title: Lessons from a blind study of simulated lenses: image reconstructions do not always reproduce true convergence

abstract: In the coming years, strong gravitational lens discoveries are expected to increase in frequency by two orders of magnitude. Lens-modelling techniques are being developed to prepare for the coming massive influx of new lens data, and blind tests of lens reconstruction with simulated data are needed for validation. In this paper we present a systematic blind study of a sample of 15 simulated strong gravitational lenses from the EAGLE suite of hydrodynamic simulations. We model these lenses with a free-form technique and evaluate reconstructed mass distributions using criteria based on shape, orientation, and lensed image reconstruction. Especially useful is a lensing analogue of the Roche potential in binary star systems, which we call the $\textit{lensing Roche potential}$. This we introduce in order to factor out the well-known problem of steepness or mass-sheet degeneracy. Einstein radii are on average well recovered with a relative error of ${\sim}5\%$ for quads and ${\sim}25\%$ for doubles; the position angle of ellipticity is on average also reproduced well up to $\pm10^{\circ}$, but the reconstructed mass maps tend to be too round and too shallow. It is also easy to reproduce the lensed images, but optimising on this criterion does not guarantee better reconstruction of the mass distribution.


title: An Extremely Massive Quiescent Galaxy at $z=3.493$: Evidence of Insufficiently Rapid Quenching Mechanisms in Theoretical Models

abstract: We present spectra of the most massive quiescent galaxy yet discovered at $z>3$, spectroscopically confirmed via the detection of Balmer absorption features in the $H-$ and $K-$bands of Keck/MOSFIRE. The spectra confirm a galaxy with no significant ongoing star formation, consistent with the lack of rest-frame UV flux and overall photometric spectral energy distribution. With a stellar mass of $3.1^{+0.1}_{-0.2} \times 10^{11} ~\rm{M}_\odot$ at $z = 3.493$, this galaxy is nearly three times more massive than the highest redshift spectroscopically confirmed absorption-line identified galaxy known. The star-formation history of this quiescent galaxy implies that it formed $>1000 ~\rm{M}_\odot$/yr for almost 0.5 Gyr beginning at $z\sim7.2$, strongly suggestive that it is the descendant of massive dusty star-forming galaxies at $5<z<7$ recently observed with ALMA. While galaxies with similarly extreme stellar masses are reproduced in some simulations at early times, such a lack of ongoing star formation is not seen there. This suggests the need for a more rapid quenching process than is currently prescribed, challenging our current understanding of how ultra-massive galaxies form and evolve in the early Universe.


title: Multiple populations in globular clusters: Unified efforts from stellar evolution and chemical evolution models

abstract: Recent stellar evolution models for globular clusters (GCs) in multiple population paradigm suggest that horizontal-branch (HB) morphology and mean period of type ab RR Lyrae variables are mostly determined by He & CNO abundances and relative ages for subpopulations. These parameters are also provided by chemical evolution models constructed to reproduce the Na-O anti-correlation. Therefore, a consistency check is possible between the synthetic HB and chemical evolution models. Furthermore, by combining them, a better constraint might be attained for star formation history and chemical abundances of subpopulations in GCs. We find, from such efforts made for four GCs, M4, M5, M15, and M80, that consistent results can be obtained from these two independent studies. In our unified model, He and Na abundances gradually increase over the generation, and therefore, the various extensions observed in both HB morphology and Na-O chemical pattern depend on the presence of later generation stars after the second generation. It is schematically shown that this observed diversity, however, would not be naturally explained by the models requiring dilution. Further spectroscopic observations are required, for metal-poor GCs in particular, to obtain a more detailed constraint from this approach.


title: Conformal Freeze-In of Dark Matter

abstract: We present the conformal freeze-in (COFI) scenario for dark matter production. At high energies, the dark sector is described by a gauge theory flowing towards a Banks-Zaks fixed point, coupled to the standard model via a non-renormalizable portal interaction. At the time when the dark sector is populated in the early universe, it is described by a strongly coupled conformal field theory. As the universe cools, cosmological phase transitions in the standard model sector, either electroweak or QCD, induce conformal symmetry breaking and confinement in the dark sector. One of the resulting dark bound states is stable on the cosmological time scales and plays the role of dark matter. With the Higgs portal, the COFI scenario provides a viable dark matter candidate with mass in a phenomenologically interesting 0.1-1 MeV range. With the quark portal, a dark matter candidate with mass around 1 keV is consistent with observations. Conformal bootstrap puts a non-trivial constraint on model building in this case.


title: Detection of Topological Materials with Machine Learning

abstract: Databases compiled using ab-initio and symmetry-based calculations now contain tens of thousands of topological insulators and topological semimetals. This makes the application of modern machine learning methods to topological materials possible. Using gradient boosted trees, we show how to construct a machine learning model which can predict the topology of a given existent material with an accuracy of 90%. Such predictions are orders of magnitude faster than actual ab-initio calculations. Through extensive testing of different models we determine which properties help detect topological materials. We identify the sources of our model's errors and we discuss approaches to overcome them.


title: Conformal $n$-point functions in momentum space

abstract: We present a Feynman integral representation for the general momentum-space scalar $n$-point function in any conformal field theory. This representation solves the conformal Ward identities and features an arbitrary function of $n(n-3)/2$ variables which play the role of momentum-space conformal cross-ratios. It involves $(n-1)(n-2)/2$ integrations over momenta, with the momenta running over the edges of an $(n-1)$-simplex. We provide the details in the simplest non-trivial case (4-point functions), and for this case we identify values of the operator and spacetime dimensions for which singularities arise leading to anomalies and beta functions, and discuss several illustrative examples from perturbative quantum field theory and holography.

categories: hep-th cond-mat.stat-mech cond-mat.str-el hep-ph


title: Coupling 1D stellar evolution with 3D-hydrodynamical simulations on-the-fly II: Stellar Evolution and Asteroseismic Applications

abstract: Models of stellar structure and evolution are an indispensable tool in astrophysics, yet they are known to incorrectly reproduce the outer convective layers of stars. In the first paper of this series, we presented a novel procedure to include the mean structure of 3D hydrodynamical simulations on-the-fly in stellar models, and found it to significantly improve the outer stratification and oscillation frequencies of a standard solar model. In the present work, we extend the analysis of the method; specifically how the transition point between envelope and interior affects the models. We confirm the versatility of our method by successfully repeating the entire procedure for a different grid of 3D hydro-simulations. Furthermore, the applicability of the procedure was investigated across the HR diagram and an accuracy comparable to the solar case was found. Moreover, we explored the implications on stellar evolution and find that the red-giant branch is shifted about 40 K to higher effective temperatures. Finally, we present for the first time an asteroseismic analysis based on stellar models fully utilising the stratification of 3D simulations on-the-fly. These new models significantly reduce the asteroseismic surface term for the two selected stars in the Kepler field. We extend the analysis to red giants and characterise the shape of the surface effect in this regime. Lastly, we stress that the interpolation required by our method would benefit from new 3D simulations, resulting in a finer sampling of the grid.


title: Energy of fermionic ground states with low-entanglement single-reference expansions, and tensor-based strong-coupling extensions of the coupled-cluster method

abstract: We consider a fermionic system for which there exist a single-reference configuration-interaction (CI) expansion of the ground state wave function that converges, albeit not necessarily rapidly, with respect to excitation number. We show that, if the coefficients of Slater determinants (SD) with $l\leq k$ excitations can be defined with a number of free parameters $N_{\leq k}$ bounded polynomially in $k$, the ground state energy $E$ only depends on a small fraction of all the wave function parameters, and is the solution of equations of the coupled-cluster (CC) form. This generalizes the standard CC method, for which $N_{\leq k}$ is bounded by a constant. Based on that result and low-rank tensor decompositions (LRTD), we discuss two possible extensions of the CC approach for wave functions with general polynomial bound for $N_{\leq k}$. First, one can use LRTD to represent the amplitudes of the CC cluster operator $T$ which, unlike in the CC case, is not truncated with respect to excitation number, and the energy and tensor parameters are given by a LRTD-adapted version of standard CC equations. Second, the LRTD can also be used to directly parametrize the CI coefficients, which involves different equations of the CC form. We derive those equations for up to quadruple-excitation coefficients, using a different type of excitation operator in the CC wave function ansatz, and a Hamiltonian representation in terms of excited particle and hole operators. We complete the proposed CC extensions by constructing compact tensor representations of coefficients, or $T$-amplitudes, using superpositions of tree tensor networks which take into account different possible types of entanglement between excited particles and holes. Finally, we discuss why the proposed extensions are theoretically applicable at larger coupling strengths than those treatable by the standard CC method.


title: Engineering Topological Models with a General-Purpose Symmetry-to-Hamiltonian Approach

abstract: Symmetry is at the heart of modern physics. Phases of matter are classified by symmetry breaking, topological phases are characterized by non-local symmetries, and point group symmetries are critical to our understanding of crystalline materials. Symmetries could then be used as a criterion to engineer quantum systems with targeted properties. Toward that end, we have developed a novel approach, the symmetric Hamiltonian construction (SHC), that takes as input symmetries, specified by integrals of motion or discrete symmetry transformations, and produces as output all local Hamiltonians consistent with these symmetries (see github.com/ClarkResearchGroup/qosy for our open-source code). This approach builds on the slow operator method [PRE 92, 012128]. We use our new approach to construct new Hamiltonians for topological phases of matter. Topological phases of matter are exotic quantum phases with potential applications in quantum computation. In this work, we focus on two types of topological phases of matter: superconductors with Majorana zero modes and $Z_2$ quantum spin liquids. In our first application of the SHC approach, we analytically construct a large and highly tunable class of superconducting Hamiltonians with Majorana zero modes with a given targeted spatial distribution. This result lays the foundation for potential new experimental routes to realizing Majorana fermions. In our second application, we find new $Z_2$ spin liquid Hamiltonians on the square and kagome lattices. These new Hamiltonians are not sums of commuting operators nor frustration-free and, when perturbed appropriately (in a way that preserves their $Z_2$ spin liquid behavior), exhibit level-spacing statistics that suggest non-integrability. This result demonstrates how our approach can automatically generate new spin liquid Hamiltonians with interesting properties not often seen in solvable models.


title: Power Counting the Small-$x$ Observables

abstract: We emphasize the importance of applying power counting to the small-$x$ observables, which introduces novel soft contributions usually missing and allows for a unified treatment of the Balitsky-Kovchegov (BK) evolution and various Sudakov logarithms. We use $pA \to h(p_{h\perp})X$ at forward rapidity to highlight how the power counting yields a partonic cross section with collinear and soft sectors. We show how the kinematic constraints can be obtained in the soft sector without violating the power counting. We further show how one can resum the threshold Sudakov logarithms systematically to all orders in a re-factorized framework with additional collinear-soft contributions. Direct applications to other small-$x$ processes involving heavy particles, jet (sub-)observables and EIC physics are straightforward.


title: PS1-STRM: Neural network source classification and photometric redshift catalogue for PS1 $3\pi$ DR1

abstract: The Pan-STARRS1 (PS1) $3\pi$ survey is a comprehensive optical imaging survey of three quarters of the sky in the $grizy$ broad-band photometric filters. We present the methodology used in assembling the source classification and photometric redshift (photo-z) catalogue for PS1 $3\pi$ Data Release 1, titled Pan-STARRS1 Source Types and Redshifts with Machine learning (PS1-STRM). For both main data products, we use neural network architectures, trained on a compilation of public spectroscopic measurements that has been cross-matched with PS1 sources. We quantify the parameter space coverage of our training data set, and flag extrapolation using self-organizing maps. We perform a Monte-Carlo sampling of the photometry to estimate photo-z uncertainty. The final catalogue contains $2,902,054,648$ objects. On our validation data set, for non-extrapolated sources, we achieve an overall classification accuracy of $98.1\%$ for galaxies, $97.8\%$ for stars, and $96.6\%$ for quasars. Regarding the galaxy photo-z estimation, we attain an overall bias of $\left<\Delta z_{\mathrm{norm}}\right>=0.0005$, a standard deviation of $\sigma(\Delta z_{\mathrm{norm}})=0.0322$, a median absolute deviation of $\mathrm{MAD}(\Delta z_{\mathrm{norm}})=0.0161$, and an outlier fraction of $O=1.89\%$. The catalogue will be made available as a high-level science product via the Mikulski Archive for Space Telescopes at https://doi.org/10.17909//t9-rnk7-gr88.


title: Realistic modeling of wind and supernovae shocks in star clusters: addressing ${\rm ^{22}Ne/^{20}Ne}$ and other problems in Galactic cosmic rays

abstract: Cosmic ray (CR) sources leave signatures in the isotopic abundances of CRs. Current models of Galactic CRs that consider supernovae (SNe) shocks as the main sites of particle acceleration cannot satisfactorily explain the higher ${\rm ^{22}Ne/^{20}Ne}$ ratio in CRs compared to the interstellar medium. Although stellar winds from massive stars have been invoked, their contribution relative to SNe ejecta has been taken as a free parameter. Here we present a theoretical calculation of the relative contributions of wind termination shocks (WTSs) and SNe shocks in superbubbles, based on the hydrodynamics of winds in clusters, the standard stellar mass function, and stellar evolution theory. We find that the contribution of WTSs towards the total CR production is at least $25\%$, which rises to $\gtrsim 50\%$ for young ($\lesssim 10$ Myr) clusters, and explains the observed $^{22}{\rm Ne}/^{20} {\rm Ne}$ ratio. We argue that since the progenitors of apparently isolated supernovae remnants (SNRs) are born in massive star clusters, both WTS and SNe shocks can be integrated into a combined scenario of CRs being accelerated in massive clusters. This scenario is consistent with the observed ratio of SNRs to $\gamma$-ray bright ($L_\gamma \gtrsim 10^{35}$ erg s$^{-1}$) star clusters, as predicted by star cluster mass function. Moreover, WTSs can accelerate CRs to PeV energies, and solve other longstanding problems of the standard supernova paradigm of CR acceleration.


title: An updated visual orbit of the directly-imaged exoplanet 51 Eridani b and prospects for a dynamical mass measurement with Gaia

abstract: We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$\sigma$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$).


title: A GMRT 150 MHz Search for Variables and Transients in Stripe 82

abstract: We have carried out a dedicated transient survey of 300 deg$^2$ of the SDSS Stripe 82 region using the Giant Meterwavelength Radio Telescope (GMRT) at 150 MHz. Our multi-epoch observations, together with the TGSS survey, allow us to probe variability and transient activity on four different timescales, beginning with 4 hours, and up to 4 years. Data calibration, RFI flagging, source finding and transient search were carried out in a semi-automated pipeline incorporating the SPAM recipe. This has enabled us to produce superior-quality images and carry out reliable transient search over the entire survey region in under 48 hours post-observation. Among the few thousand unique point sources found in our 5$\sigma$ single-epoch catalogs (flux density thresholds of about 24 mJy, 20 mJy, 16 mJy and 18 mJy on the respective timescales), we find $<$0.08%, 0.01%, $<$0.06% and 0.05% to be variable (beyond a significance of 4$\sigma$ and fractional variability of 30%) on timescales of 4 hours, 1 day, 1 month and 4 years respectively. This is substantially lower than that in the GHz sky, where $\sim$1% of the persistent point sources are found to be variable. Although our survey was designed to probe a superior part of the transient phase space, our transient sources did not yield any significant candidates. The transient (preferentially extragalactic) rate at 150 MHz is therefore $<$0.005 on timescales of 1 month and 4 years, and $<$0.002 on timescales of 1 day and 4 hours, beyond 7$\sigma$ detection threshold. We put these results in the perspective with the previous studies and give recommendations for future low-frequency transient surveys.


title: The Effects of Population III X-ray and Radio Backgrounds on the Cosmological 21-cm Signal

abstract: We investigate the effects of Population III (Pop III) stars and their remnants on the cosmological 21-cm global signal. By combining a semi-analytic model of Pop III star formation with a global 21-cm simulation code, we investigate how X-ray and radio emission from accreting Pop III black holes may affect both the timing and depth of the 21-cm absorption feature that follows the initial onset of star formation during the Cosmic Dawn. We compare our results to the findings of the EDGES experiment, which has reported the first detection of a cosmic 21-cm signal. In general, we find that our fiducial Pop III models, which have peak star formation rate densities of $\sim 10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ between $z \sim 10$ and $z \sim 15$, are able to match the timing of the EDGES signal quite well, in contrast to models that ignore Pop III stars. To match the unexpectedly large depth of the EDGES signal without recourse to exotic physics, we vary the parameters of emission from accreting black holes (formed as Pop III remnants) including the intrinsic strength of X-ray and radio emission as well as the local column density of neutral gas. We find that models with strong radio emission and relatively weak X-ray emission can self-consistently match the EDGES signal, though this solution requires fine-tuning. We are only able to produce signals with sharp features similar to the EDGES signal if the Pop~III IMF is peaked narrowly around $140 \, M_\odot$.


title: Detection of a low-mass stellar companion to the accelerating A2IV star HR 1645

abstract: The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions.


title: Low-rank Sachdev-Ye-Kitaev models

abstract: Motivated by recent works on atom-cavity realizations of fast scramblers, and on Cooper pairing in non-Fermi liquids, we study a family of solvable variants of the ($q=4$) Sachdev-Ye-Kitaev model in which the rank and eigenvalue distribution of the coupling matrix $J_{ij,kl}$ are tuneable. When the rank is proportional to the number of fermions, the low temperature behavior is sensitive to the eigenvalue distribution. We obtain a complete classification of the possible non-Fermi liquid quantum phases. These include two previously studied phases whose fermion scaling dimension depends continuously on the rank; we show that they are maximally chaotic, but necessitate {an extensively degenerate or negative semidefinite coupling matrix}. More generic distributions give rise to "almost Fermi liquids" with a scaling dimension $\Delta = 1/2$, but which differ from a genuine Fermi-liquid in quasi-particle decay rate, quantum Lyapunov exponent and/or specific heat.


title: Leveraging directed causal discovery to detect latent common causes

abstract: The discovery of causal relationships is a fundamental problem in science and medicine. In recent years, many elegant approaches to discovering causal relationships between two variables from observational data have been proposed. However, most of these deal only with purely directed causal relationships and cannot detect latent common causes. Here, we devise a general heuristic which takes a causal discovery algorithm that can only distinguish purely directed causal relations and modifies it to also detect latent common causes. We apply our method to two directed causal discovery algorithms, the Information Geometric Causal Inference of (Daniusis et al., 2010) and the Kernel Conditional Deviance for Causal Inference of (Mitrovic, Sejdinovic, & Teh, 2018), and extensively test on synthetic data -- detecting latent common causes in additive, multiplicative and complex noise regimes -- and on real data, where we are able to detect known common causes. In addition to detecting latent common causes, our experiments demonstrate that both the modified algorithms preserve the performance of the original in distinguishing directed causal relations.


title: Constraining black hole-galaxy scaling relations from the large-scale clustering of Active Galactic Nuclei and implied mean radiative efficiency

abstract: A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demography of supermassive black holes thus preventing accurate model comparisons and progress in this field. Here we show that the large-scale spatial distribution of local active galactic nuclei (AGN), believed to be accreting supermassive black holes, can constrain the shape and normalization of the black hole-stellar mass relation thus bypassing resolution-related observational biases. In turn, our results can set more stringent constraints on the so-called "radiative efficiency", a fundamental parameter describing the inner physics of supermassive black holes that is closely linked to their spin, geometry, and ability to release energy. The mean value of the radiative efficiency can be estimated by comparing the average total luminous output of AGN with the relic mass density locked up in quiescent supermassive black holes at galaxy centres today. For currently accepted values of the AGN obscured fractions and bolometric corrections, our newest estimates of the local supermassive black hole mass density favour mean radiative efficiencies of ~10-20\%, suggesting that the vast majority of supermassive black holes are spinning moderately to rapidly. With large-scale AGN surveys coming online, our novel methodology will enable even tighter constraints on the fundamental parameters that regulate the growth of supermassive black holes.


title: Review of Recent Heap Specification and Verification Techniques

abstract: The article provides an overview of the existing methods of dynamic memory verification; a comparative analysis is carried out; the applicability for solving problems of control, monitoring, and verification of dynamic memory is evaluated. This article is divided into eight sections. The first section introduces formal verification, followed by a section that discusses dynamic memory management problems. The third section discusses Hoare's calculus resumed by heap transformations to the stack. The fifth and sixth sections introduce the concept of dynamic memory shape analysis and the rotation of pointers. The seventh is on separation logic. The last section discusses possible areas of further research, particularly the recognition at recording level of various instances of objects; automation of proofs; "hot" code, that is, software code that updates itself when the program runs; expanding intuitiveness, for instance, on proof explanations.


title: Direct Detection of Solar Angular Momentum Loss with the Wind Spacecraft

abstract: The rate at which the solar wind extracts angular momentum from the Sun has been predicted by theoretical models for many decades, and yet we lack a conclusive measurement from in-situ observations. In this letter we present a new estimate of the time-varying angular momentum flux in the equatorial solar wind, as observed by the \textit{Wind} spacecraft from 1994-2019. We separate the angular momentum flux into contributions from the protons, alpha particles, and magnetic stresses, showing that the mechanical flux in the protons is $\sim$3 times larger than the magnetic field stresses. We observe the tendency for the angular momentum flux of fast wind streams to be oppositely signed to the slow wind streams, as noted by previous authors. From the average total flux, we estimate the global angular momentum loss rate of the Sun to be $3.3\times10^{30}$erg, which lies within the range of various MHD wind models in the literature. This angular momentum loss rate is a factor of $\sim$2 weaker than required for a Skumanich-like rotation period evolution ($\Omega_*\propto $ stellar age$^{-1/2}$), which should be considered in studies of the rotation period evolution of Sun-like stars.


title: The Proton in High Definition: Revisiting Photon-Initiated Production in High Energy Collisions

abstract: We re-examine the current state of the art for the calculation of photon-initiated processes at the LHC, as formulated in terms of a photon PDF in the proton that may be determined rather precisely from the known proton structure functions. We in particular demonstrate that a by construction more precise calculation is provided by a direct application of the structure function approach, best known from the case of Higgs Boson production via vector boson fusion. This avoids any artificial scale variation uncertainties, which can otherwise be rather significant for processes calculated within the standard approach thus far. To understand the source of these, we present a detailed comparison of the structure function approach and its relation to the photon PDF. We then provide precise predictions for the photon-initiated contribution to lepton pair production at the LHC, including the lepton pair transverse momentum distribution. Thus, by a direct application of the structure function formalism we show how the contribution from initial-state photons at the LHC may for the first time be included with high precision in a universal and straightforward way, providing a high definition picture of the photon content of the proton.


title: Classification of Critical Points in Energy Bands Based on Topology, Scaling and Symmetry

abstract: A critical point of the energy dispersion is the momentum where electron velocity vanishes. At the corresponding energy, the density of states (DOS) exhibits non-analyticity such as divergence. Critical points can be first classified as ordinary and high-order ones, and the ordinary critical points have been studied thoroughly by L\'eon van Hove. In this work, we describe and classify high-order critical points based on topology, scaling and symmetry, which are beyond L\'eon van Hove's framework. We show that high-order critical points can have power-law divergent DOS with particle-hole asymmetry, and can be realized at generic or symmetric momenta by tuning a few parameters such as twist angle, strain, pressure and/or external fields.


title: Identifying Higher Order Topology and Fractional Corner Charge Using Entanglement Spectra

abstract: We study the entanglement spectrum (ES) of two-dimensional $C_{n}$-symmetric second-order topological insulators (TIs). We show that some characteristic higher order topological observables, e.g., the filling anomaly and its associated fractional corner charge, can be determined from the ES of atomic and fragile TIs. By constructing the relationship between the configuration of Wannier orbitals and the number of protected in-gap states in the ES for different symmetric cuts in real space, we express the fractional corner charge in terms of the number of protected in-gap states of the ES. We show that our formula is robust in the presence of electron-electron interactions as long as the interactions preserve $C_{n}$ rotation symmetry and charge-conservation symmetry. Moreover, we discuss the possible signatures higher order topology in the many-body ES. Our methods allow the identification of some classes of higher order topology without requiring the usage of nested Wilson loops or nested entanglement spectra.


title: Stable limit theorems on the Poisson space

abstract: We prove limit theorems for functionals of a Poisson point process using the Malliavin calculus on the Poisson space. The target distribution is conditionally either a Gaussian vector or a Poisson random variable. The convergence is stable and our conditions are expressed in terms of the Malliavin operators. For conditionally Gaussian limits, we also obtain quantitative bounds, given for the Monge-Kantorovich transport distance in the univariate case; and for another probabilistic variational distance in higher dimension. Our work generalizes several limit theorems on the Poisson space, including the seminal works by Peccati, Sol\'e, Taqqu & Utzet for Gaussian approximations; and by Peccati for Poisson approximations; as well as the recently established fourth-moment theorem on the Poisson space of D\"obler & Peccati. We give an application to stochastic processes.


title: Rational fibered cubic fourfolds with nontrivial Brauer classes

abstract: Some classes of cubic fourfolds are birational to fibrations over $P^2$, where the fibers are rational surfaces. This relates strongly the rationality of the cubic with the rationality of these surfaces over the function field of $P^2$ and to the existence of rational sections of the fibration. It is an interesting question to understand whether the existence of a rational section is also necessary for the rationality of the cubic. Let $C_{18}$ denote the divisor inside the moduli space of smooth cubic hypersurfaces in $P^5$ whose general element is a smooth cubic containing an elliptic ruled surface associated to a sextic del Pezzo fibration. Cubics in $C_{18}$ are rational whenever the fibration has rational sections. By intersecting $C_{18}$ with the divisors $C_{14}$, $C_{26}$ and $C_{38}$, whose elements are known to be rational, we provide explicit description of these intersections and we exhibit new examples of rational cubic fourfolds whose associated good sextic del Pezzo fibration has no rational sections. We study also the divisor $C_{8}$, whose general element is a cubic containing a plane. These cubics are birational to quadric fibrations over $P^2$. Intersections between $C_8$ and the divisors $C_{26}$, $C_{38}$ give us also new classes of smooth rational cubic fourfolds whose associated quadric surface bundle doesn't have a rational section.


title: Catastrophe theory classification of Fermi surface topological transitions in two dimensions

abstract: We classify all possible singularities in the electronic dispersion of two-dimensional systems that occur when the Fermi surface changes topology, using catastrophe theory. For systems with up to seven control parameters (i.e., pressure, strain, bias voltage, etc), the theory guarantees that the singularity belongs to to one of seventeen standard types. We show that at each of these singularities the density of states diverges as a power law, with a universal exponent characteristic of the particular catastrophe, and we provide its universal ratio of amplitudes of the prefactors of energies above and below the singularity. We further show that crystal symmetry restricts which types of catastrophes can occur at the points of high symmetry in the Brillouin zone. For each of the seventeen wallpaper groups in two-dimensions, we list which catastrophes are possible at each high symmetry point.


title: Virtual Elements and curved edges

abstract: We present a new method to construct Virtual Element spaces on polygons with curved edges.


title: Quantum Log-Corrections to Swampland Conjectures

abstract: Taking the anti-de Sitter minimum of KKLT and the large volume scenario at face value, we argue for the existence of logarithmic quantum corrections to AdS swampland conjectures. If these conjectures receive such corrections, it is natural to suspect that they also arise for other swampland conjectures, in particular the dS swampland conjecture. We point out that the proposed $\log$-corrections are in accord with the implications of the recently proposed trans-Planckian censorship conjecture. We also comment on the emergence proposal in the context of both perturbative flux models and the KKLT construction.


title: Power oscillations induced by the relative Goos-Haenchen phase

abstract: By using an optical interferometer composed of a dielectric laser ellipsometer, to change the optical response of transverse electric and magnetic incident radiation, and two polarisers, to trigger the interference pattern induced by the relative Goos-Haenchen phase, we show under which conditions it is possible to optimize the laser power oscillations induced by the relative phase difference between orthogonal polarised states. The Goos-Haenchen interference can then be used to sense rotation, to test optical components, and to simulate quarter and half wave plates.


title: Fast and Automatic Periacetabular Osteotomy Fragment Pose Estimation Using Intraoperatively Implanted Fiducials and Single-View Fluoroscopy

abstract: Accurate and consistent mental interpretation of fluoroscopy to determine the position and orientation of acetabular bone fragments in 3D space is difficult. We propose a computer assisted approach that uses a single fluoroscopic view and quickly reports the pose of an acetabular fragment without any user input or initialization. Intraoperatively, but prior to any osteotomies, two constellations of metallic ball-bearings (BBs) are injected into the wing of a patient's ilium and lateral superior pubic ramus. One constellation is located on the expected acetabular fragment, and the other is located on the remaining, larger, pelvis fragment. The 3D locations of each BB are reconstructed using three fluoroscopic views and 2D/3D registrations to a preoperative CT scan of the pelvis. The relative pose of the fragment is established by estimating the movement of the two BB constellations using a single fluoroscopic view taken after osteotomy and fragment relocation. BB detection and inter-view correspondences are automatically computed throughout the processing pipeline. The proposed method was evaluated on a multitude of fluoroscopic images collected from six cadaveric surgeries performed bilaterally on three specimens. Mean fragment rotation error was 2.4 +/- 1.0 degrees, mean translation error was 2.1 +/- 0.6 mm, and mean 3D lateral center edge angle error was 1.0 +/- 0.5 degrees. The average runtime of the single-view pose estimation was 0.7 +/- 0.2 seconds. The proposed method demonstrates accuracy similar to other state of the art systems which require optical tracking systems or multiple-view 2D/3D registrations with manual input. The errors reported on fragment poses and lateral center edge angles are within the margins required for accurate intraoperative evaluation of femoral head coverage.


title: Intertwined magnetic, structural, and electronic transitions in V$_2$O$_3$

abstract: We present a coordinated study of the paramagnetic-to-antiferromagnetic, rhombohedral-to-monoclinic, and metal-to-insulator transitions in thin-film specimens of the classic Mott insulator V$_2$O$_3$ using low-energy muon spin relaxation, x-ray diffraction, and nanoscale-resolved near-field infrared spectroscopic techniques. The measurements provide a detailed characterization of the thermal evolution of the magnetic, structural, and electronic phase transitions occurring in a wide temperature range, including quantitative measurements of the high- and low-temperature phase fractions for each transition. The results reveal a stable coexistence of the high- and low-temperature phases over a broad temperature range throughout the transition. Careful comparison of temperature dependence of the different measurements, calibrated by the resistance of the sample, demonstrates that the electronic, magnetic, and structural degrees of freedom remain tightly coupled to each other during the transition process. We also find evidence for antiferromagnetic fluctuations in the vicinity of the phase transition, highlighting the important role of the magnetic degree of freedom in the metal-insulator transition.


title: Rigidity of the Torelli subgroup in $Out(F_N)$

abstract: Let $N$ be at least 4. We prove that every injective homomorphism from the Torelli subgroup into $Out(F_N)$ differs from the inclusion by a conjugation in $Out(F_N)$. This applies more generally to the following subgroups: every finite-index subgroup of $Out(F_N)$ (recovering a theorem of Farb and Handel); every subgroup that contains a finite-index subgroup of one of the groups in the Andreadakis--Johnson filtration; every subgroup that contains a power of every linearly-growing automorphism; more generally, every twist-rich subgroup (subgroups that contain sufficiently many twists in an appropriate sense). Among applications, this recovers the fact that the abstract commensurator of every group above is equal to its relative commensurator in $Out(F_N)$; it also implies that all subgroups in the Andreadakis--Johnson filtration are co-Hopfian. We also prove the same rigidity statement for subgroups of $Out(F_3)$ which contain a power of every Nielsen transformation. This shows in particular that $Out(F_3)$ and all its finite-index subgroups are co-Hopfian, extending a theorem of Farb and Handel to the $N=3$ case.


title: Kuksa: A Cloud-Native Architecture for Enabling Continuous Delivery in the Automotive Domain

abstract: Connecting vehicles to cloud platforms has enabled innovative business scenarios while raising new quality concerns, such as reliability and scalability, which must be addressed by research. Cloud-native architectures based on microservices are a recent approach to enable continuous delivery and to improve service reliability and scalability. We propose an approach for restructuring cloud platform architectures in the automotive domain into a microservices architecture. To this end, we adopted and implemented microservices patterns from literature to design the cloud-native automotive architecture and conducted a laboratory experiment to evaluate the reliability and scalability of microservices in the context of a real-world project in the automotive domain called Eclipse Kuksa. Findings indicated that the proposed architecture could handle the continuous software delivery over-the-air by sending automatic control messages to a vehicular setting. Different patterns enabled us to make changes or interrupt services without extending the impact to others. The results of this study provide evidences that microservices are a potential design solution when dealing with service failures and high payload on cloud-based services in the automotive domain.


title: Exploring the possibility of enhancing the figure-of-merit ( $>$ 2) of Na$_{0.74}$CoO$_{2}$: A combined experimental and theoretical study

abstract: Search of new thermoelectric (TE) materials with high \textit{figure-of-merit} (ZT) is always inspired the researcher in TE field. Here, we present a combined experimental and theoretical study of TE properties of Na$_{0.74}$CoO$_{2}$ compound in high-temperature region. The experimental Seebeck coefficient (S) is found to vary from 64 to 118 $\mu$V/K in the temperature range $300-620$ K. The positive values of S are indicating the dominating p-type behaviour of the compound. The observed value of thermal conductivity ($\kappa$) is $\sim$ 2.2 W/m-K at 300 K. In the temperature region $300-430$ K, the value of $\kappa$ increases up to $\sim$ 2.6 W/m-K and then decreases slowly till 620 K with the corresponding value of $\sim$ 2.4 W/m-K. We have also carried out the theoretical calculations and the best matching between experimental and calculated values of transport properties are observed in spin-polarized calculation within DFT+\textit{U} by chosen \textit{U} = 4 eV. The maximum calculated value of ZT is found to be $\sim$ 0.67 at 1200 K for p-type conduction. Our computational study suggests that the possibility of n-type behaviour of the compound which can lead to a large value of ZT at higher temperature region. Electron doping of $\sim$ 5.1$\times$10$^{20}$ cm$^{-3}$ is expected to give rise the high ZT value of $\sim$ 2.7 at 1200 K. Using these temperature-dependent ZT values, we have calculated the maximum possible values of efficiency ($\eta$) of thermoelectric generator (TEG) made by p and n-type Na$_{0.74}$CoO$_{2}$. The present study suggests that one can get the efficiency of a TE cell as high as $\sim$ 11$\%$ when the cold and hot end temperature are fixed at 300 K and 1200 K, respectively. Such high values of ZT and efficiency suggest that Na$_{0.74}$CoO$_{2}$ can be used as a potential candidate for high-temperature TE applications.


title: Attractor Mechanism and Non-Renormalization Theorem in 6D (1,0) Supergravity

abstract: We compute the macroscopic entropy of the supersymmetric rotating dyonic strings carrying linear momentum in 6D (1,0) supergravity with curvature squared corrections. Our calculation is based on Sen's entropy function formalism applied to the near-horizon geometry of the string solution taking the form of an extremal BTZ$\times S^3$. The final entropy formula states that the two independent supersymmetric completions of Riemann tensor squared contribute equally to the entropy. A further $S^3$ compactification of the 6D theory results in a matter coupled 3D supergravity model in which the quantization condition of the SU(2)$_R$ Chern-Simons level implies the horizon value of the dilaton is not modified by higher derivative interactions beyond supersymmetric curvature squared terms.


title: First M87 Event Horizon Telescope Results and the Role of ALMA

abstract: In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.


title: Teach Biped Robots to Walk via Gait Principles and Reinforcement Learning with Adversarial Critics

abstract: Controlling a biped robot to walk stably is a challenging task considering its nonlinearity and hybrid dynamics. Reinforcement learning can address these issues by directly mapping the observed states to optimal actions that maximize the cumulative reward. However, the local minima caused by unsuitable rewards and the overestimation of the cumulative reward impede the maximization of the cumulative reward. To increase the cumulative reward, this paper designs a gait reward based on walking principles, which compensates the local minima for unnatural motions. Besides, an Adversarial Twin Delayed Deep Deterministic (ATD3) policy gradient algorithm with a recurrent neural network (RNN) is proposed to further boost the cumulative reward by mitigating the overestimation of the cumulative reward. Experimental results in the Roboschool Walker2d and Webots Atlas simulators indicate that the test rewards increase by 23.50% and 9.63% after adding the gait reward. The test rewards further increase by 15.96% and 12.68% after using the ATD3_RNN, and the reason may be that the ATD3_RNN decreases the error of estimating cumulative reward from 19.86% to 3.35%. Besides, the cosine kinetic similarity between the human and the biped robot trained by the gait reward and ATD3_RNN increases by over 69.23%. Consequently, the designed gait reward and ATD3_RNN boost the cumulative reward and teach biped robots to walk better.


title: The Graphon Fourier Transform

abstract: In many network problems, graphs may change by the addition of nodes, or the same problem may need to be solved in multiple similar graphs. This generates inefficiency, as analyses and systems that are not transferable have to be redesigned. To address this, we consider graphons, which are both limit objects of convergent graph sequences and random graph models. We define graphon signals and introduce the Graphon Fourier Transform (WFT), to which the Graph Fourier Transform (GFT) is shown to converge. This result is demonstrated in two numerical experiments where, as expected, the GFT converges, hinting to the possibility of centralizing analysis and design on graphons to leverage transferability.


title: No-regret Non-convex Online Meta-Learning

abstract: The online meta-learning framework is designed for the continual lifelong learning setting. It bridges two fields: meta-learning which tries to extract prior knowledge from past tasks for fast learning of future tasks, and online-learning which deals with the sequential setting where problems are revealed one by one. In this paper, we generalize the original framework from convex to non-convex setting, and introduce the local regret as the alternative performance measure. We then apply this framework to stochastic settings, and show theoretically that it enjoys a logarithmic local regret, and is robust to any hyperparameter initialization. The empirical test on a real-world task demonstrates its superiority compared with traditional methods.


title: Hybrid data-driven physics model-based framework for enhance cyber-physical smart grid security

abstract: This paper presents a hybrid data-driven physics model-based framework for real time monitoring in smart grids. As the power grid transitions to the use of smart grid technology, it's real time monitoring becomes more vulnerable to cyber attacks like false data injections (FDI). Although smart grids cyber-physical security has an extensive scope, this paper focuses on FDI attacks, which are modeled as bad data. State of the art strategies for FDI detection in real time monitoring rely on physics model-based weighted least squares state estimation solution and statistical tests. This strategy is inherently vulnerable by the linear approximation and the companion statistical modeling error, which means it can be exploited by a coordinated FDI attack. In order to enhance the robustness of FDI detection, this paper presents a framework which explores the use of data-driven anomaly detection methods in conjunction with physics model-based bad data detection via data fusion. Multiple anomaly detection methods working at both the system level and distributed local detection level are fused. The fusion takes into consideration the confidence of the various anomaly detection methods to provide the best overall detection results. Validation considers tests on the IEEE 118 bus system.


title: Critical force in active microrheology

abstract: Soft solids like colloidal glasses exhibit a yield stress, above which the system starts to flow. The microscopic analogon in microrheology is the delocalization of a tracer particle subject to an external force exceeding a threshold value, in a glassy host. We characterize this delocalization transition based on a bifurcation analysis of the corresponding mode-coupling theory equations. A schematic model is presented first, that allows analytical progress, and the full physical model is studied numerically next. This analysis yields a continuous type A transition with a critical power law decay of the probe correlation functions with exponent $-1/2$. In order to compare with simulations with a limited duration, a finite time analysis is performed, which yields reasonable results for not-too-small wave vectors. The theoretically predicted findings are verified by Langevin dynamics simulations. For small wave vectors we find anomalous behavior for the probe position correlation function, which can be traced back to a wave vector divergence of the critical amplitude. In addition we propose and test three methods to extract the critical force from experimental data, which provide the same value of the critical force when applied to the finite-time theory or simulations.


title: Developing AlMn films for Argonne TES fabrication

abstract: The reference design for the next-generation cosmic microwave background (CMB) experiment, CMB-S4, relies on large arrays of transition edge sensor (TES) bolometers coupled to Superconducting Quantum Interference Device (SQUID)-based readout systems. Mapping the CMB to near cosmic variance limits will enable the search for signatures of inflation and constrain dark energy and neutrino physics. AlMn TESes provide simple film manufacturing and highly uniform arrays over large areas to meet the requirements of the CMB-S4 experiment. TES parameters such as critical temperature and normal resistance must be tuned to experiment specifications and can be varied based on geometry and steps in the fabrication process such as deposition layering, geometry, and baking time and temperature. Using four-terminal sensing, we measured $T_C$ and $R_N$ of AlMn 2000 ppm films and devices of varying thicknesses fabricated at Argonne National Laboratory to motivate device geometries and fabrication processes to tune $T_C$ to 150-200 mK and $R_N$ to $\sim$10 mOhms. Measurements of IV curves and time constants for the resulting devices of varying leg length were made using time-division SQUID multiplexing, and determined $T_C$, $G$, $k$, $f_{3db}$, and $R_N$. We present the results of these tests along with the geometries and fabrication steps used to tune the device parameters to the desired limits.


title: $n$-ary algebras of the first level

abstract: $n$-ary algebras of the first degeneration level are described. A detailed classification is given in the cases $n=2,3$.


title: Inference with selection, varying population size and evolving population structure: Application of ABC to a forward-backward coalescent process with interactions

abstract: Genetic data are often used to infer demographic history and changes or detect genes under selection. Inferential methods are commonly based on models making various strong assumptions: demography and population structures are supposed \textit{a priori} known, the evolution of the genetic composition of a population does not affect demography nor population structure, and there is no selection nor interaction between and within genetic strains. In this paper, we present a stochastic birth-death model with competitive interactions and asexual reproduction. We develop an inferential procedure for ecological, demographic and genetic parameters. We first show how genetic diversity and genealogies are related to birth and death rates, and to how individuals compete within and between strains. {This leads us to propose an original model of phylogenies, with trait structure and interactions, that allows multiple merging}. Second, we develop an Approximate Bayesian Computation framework to use our model for analyzing genetic data. We apply our procedure to simulated data from a toy model, and to real data by analyzing the genetic diversity of microsatellites on Y-chromosomes sampled from Central Asia human populations in order to test whether different social organizations show significantly different fertility.


title: Complex Transformer: A Framework for Modeling Complex-Valued Sequence

abstract: While deep learning has received a surge of interest in a variety of fields in recent years, major deep learning models barely use complex numbers. However, speech, signal and audio data are naturally complex-valued after Fourier Transform, and studies have shown a potentially richer representation of complex nets. In this paper, we propose a Complex Transformer, which incorporates the transformer model as a backbone for sequence modeling; we also develop attention and encoder-decoder network operating for complex input. The model achieves state-of-the-art performance on the MusicNet dataset and an In-phase Quadrature (IQ) signal dataset.


title: Pattern Formation in a Coupled Membrane-Bulk Reaction-Diffusion Model for Intracellular Polarization and Oscillations

abstract: Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a 1-D model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a 2-D model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov--Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns.


title: Symmetrisation and the Feigin-Frenkel centre

abstract: For simple Lie algebras of types B, C, and D, we provide new explicit formulas for the generators of the Feigin-Frenkel centre. These formulas make use of the symmetrisation map as well as some well-chosen symmetric invariants of $\mathfrak g$. There are some general results on the role of the symmetrisation map in the explicit description of the FF-centre. The advantage of our method is that it reduces questions about elements of the FF-centre to questions on the structure of the symmetric invariants in a type-free way. As an illustration, we deal with type G$_2$ by hand.


title: Analytical Study of the Impacts of Stochastic Load Fluctuation on the Dynamic Voltage Stability Margin Using Bifurcation Theory

abstract: This paper studies the impacts of stochastic load fluctuations, namely the fluctuation intensity and the load power variation speed, on power system dynamic voltage stability. Additionally, the trade-off relationship between the two parameters is revealed, which provides important insights regarding the potential of using energy storage to maintain voltage stability under high uncertainty. To this end, Stochastic Differential-Algebraic Equations (SDAEs) are used to model the stochastic load variation; bifurcation analysis is carried out to explain the influence of stochasticity. Numerical study and Monte Carlo simulations on the IEEE 14-bus system demonstrate that a larger fluctuation intensity or a slower load power variation speed may lead to a smaller voltage stability margin. To the best of authors' knowledge, this work uncovers the impacts of the time evolution property of the driving parameters, i.e., the load power variation speed and its trade off effect with the fluctuation intensity on the size of the dynamic voltage stability margin.


title: Facets of the Witt Vectors

abstract: We give a $K$-theoretic account of the basic properties of Witt vectors. Along the way we re-prove basic properties of the little-known Witt vector norm, give a characterization of Witt vectors in terms of algebraic $K$-theory, and a presentation of the Witt vectors that is useful for computation.


title: Radiative pulsar magnetospheres: aligned rotator

abstract: Force-free neutron star magnetospheres are nowadays well known and found through numerical simulations. Even extension to general relativity has recently been achieved. However, those solutions are by definition dissipationless, meaning that the star is unable to accelerate particles and let them radiate any photon. Interestingly, the force-free model has no free parameter however it must be superseded by a dissipative mechanism within the plasma. In this paper, we investigate the magnetosphere electrodynamics for particles moving in the radiation reaction regime, using the limit where acceleration is fully balanced by radiation, also called Aristotelian dynamics. An Ohm's law is derived, from which the dissipation rate is controlled by a one parameter family of solutions depending on the pair multiplicity~$\kappa$. The spatial extension of the dissipation zone is found self-consistently from the simulations. We show that the radiative magnetosphere of an aligned rotator tends to the force-free regime whenever the pair multiplicity becomes moderately large, $\kappa \gg 1$. However, for low multiplicity, a substantial fraction of the spindown energy goes into particle acceleration and radiation in addition to the Poynting flux, the latter remaining only dominant for large multiplicities. We show that the work done on the plasma occurs predominantly in the equatorial current sheet right outside the light-cylinder.


title: Lucene for Approximate Nearest-Neighbors Search on Arbitrary Dense Vectors

abstract: We demonstrate three approaches for adapting the open-source Lucene search library to perform approximate nearest-neighbor search on arbitrary dense vectors, using similarity search on word embeddings as a case study. At its core, Lucene is built around inverted indexes of a document collection's (sparse) term-document matrix, which is incompatible with the lower-dimensional dense vectors that are common in deep learning applications. We evaluate three techniques to overcome these challenges that can all be natively integrated into Lucene: the creation of documents populated with fake words, LSH applied to lexical realizations of dense vectors, and k-d trees coupled with dimensionality reduction. Experiments show that the "fake words" approach represents the best balance between effectiveness and efficiency. These techniques are integrated into the Anserini open-source toolkit and made available to the community.


title: Towards an Intelligent Microscope: adaptively learned illumination for optimal sample classification

abstract: Recent machine learning techniques have dramatically changed how we process digital images. However, the way in which we capture images is still largely driven by human intuition and experience. This restriction is in part due to the many available degrees of freedom that alter the image acquisition process (lens focus, exposure, filtering, etc). Here we focus on one such degree of freedom - illumination within a microscope - which can drastically alter information captured by the image sensor. We present a reinforcement learning system that adaptively explores optimal patterns to illuminate specimens for immediate classification. The agent uses a recurrent latent space to encode a large set of variably-illuminated samples and illumination patterns. We train our agent using a reward that balances classification confidence with image acquisition cost. By synthesizing knowledge over multiple snapshots, the agent can classify on the basis of all previous images with higher accuracy than from naively illuminated images, thus demonstrating a smarter way to physically capture task-specific information.


title: From Random to Regular: Variation in the Patterning of Retinal Mosaics

abstract: The various types of retinal neurons are each positioned at their respective depths within the retina where they are believed to be assembled as orderly mosaics, in which like-type neurons minimize proximity to one another. Two common statistical analyses for assessing the spatial properties of retinal mosaics include the nearest neighbor analysis, from which an index of their "regularity" is commonly calculated, and the density recovery profile derived from auto-correlation analysis, revealing the presence of an exclusion zone indicative of anti-clustering. While each of the spatial statistics derived from these analyses, the regularity index and the effective radius, can be useful in characterizing such properties of orderly retinal mosaics, they are rarely sufficient for conveying the natural variation in the self-spacing behavior of different types of retinal neurons and the extent to which that behavior generates uniform intercellular spacing across the mosaic. We consider the strengths and limitations of different spatial statistical analyses for assessing the patterning in retinal mosaics, highlighting a number of misconceptions and their frequent misuse. Rather than being diagnostic criteria for determining simply whether a population is "regular", they should be treated as descriptive statistics that convey variation in the factors that influence neuronal positioning. We subsequently apply multiple spatial statistics to the analysis of eight different mosaics in the mouse retina, demonstrating conspicuous variability in the degree of patterning present, from essentially random to notably regular. This variability in patterning has both a developmental as well as a functional significance, reflecting the rules governing the positioning of different types of neurons as the architecture of the retina is assembled (abstract truncated).


title: Better Approximate Inference for Partial Likelihood Models with a Latent Structure

abstract: Temporal Point Processes (TPP) with partial likelihoods involving a latent structure often entail an intractable marginalization, thus making inference hard. We propose a novel approach to Maximum Likelihood Estimation (MLE) involving approximate inference over the latent variables by minimizing a tight upper bound on the approximation gap. Given a discrete latent variable $Z$, the proposed approximation reduces inference complexity from $O(|Z|^c)$ to $O(|Z|)$. We use convex conjugates to determine this upper bound in a closed form and show that its addition to the optimization objective results in improved results for models assuming proportional hazards as in Survival Analysis.


title: Transparent and flexible field-effect transistors and mem-transistors with electroactive layers of solution-processed organic polyradicals

abstract: Solution processed polymers are at the core of organic electronics. Polyradicals (polymers in which each repeating unit contains an unpaired spin) are unique alternatives to their p-conjugated, semiconducting counterparts. Unique of polyradicals are tunable charge states localized at their repeating units, which enable electrically switchable charge transport regimes. Tremendous efforts were focused on polyradical memristors and batteries. Notwithstanding recent progress in doping, polyradical field-effect transistors (PR-FETs) have not been reported. Here, we show that vertical architectures, with drain-source contacts sandwiching the active layer of a strongly correlated 6-oxoverdazyl polyradical, leads to on/off ratios higher than 10+3 in p-type PR-FETs. Hole injection occurs via contact doping by tunable charge states at the polyradical-electrode interface. Transparent and flexible PR-FETs are also reported. PR-FETs are superior to existing organic FETs as they combine memristor and transistor functions in one mem-transistor device, offering a unique potential for the circuital simplification of organic electronics.


title: Conditional Law and Occupation Times of Two-Sided Sticky Brownian Motion

abstract: Sticky Brownian motion on the real line can be obtained as a weak solution of a system of stochastic differential equations. We find the conditional distribution of the process given the driving Brownian motion, both at an independent exponential time and at a fixed time t>0. As a classical problem, we find the distribution of the occupation times of a half-line, and at 0, which is the sticky point for the process.


title: Exponential Dynamical Localization for Random Word Models

abstract: We show that one-dimensional Schr{\"o}dinger operators whose potentials arise by randomly concatenating words from an underlying set exhibit exponential dynamical localization (EDL) on any compact set which trivially intersects a finite set of critical energies. We do so by first giving a new proof of spectral localization for such operators and then showing that once one has the existence of a complete orthonormal basis of eigenfunctions (with probability one), the same estimates used to prove it naturally lead to a proof of the aforementioned EDL result. The EDL statements provide new localization results for several classes of random Schr{\"o}dinger operators including random polymer models and generalized Anderson models.


title: Configuration interaction approaches for solving quantum impurity models

abstract: We develop several configuration interaction approaches for characterizing the electronic structure of an adsorbate on a metal surface (at least in model form). When one can separate adsorbate from substrate, these methods can achieve a reasonable description of adsorbate on-site electron-electron correlation in the presence of a continuum of states. While the present paper is restricted to the Anderson impurity model, there is hope that these methods can be extended to ab initio Hamiltonians, and provide insight into the structure and dynamics of molecule-metal surface interactions.


title: Solitons in the Einstein universe

abstract: We show that equations of Newtonian hydrodynamics and gravity with Einstein's cosmological constant included admit gravitostatic wave solutions propagating in the background of Einstein's static Universe. In the zero pressure limit these waves exist at an average matter density exceeding that of Einstein's Universe. They have the form of a lattice of integrable density singularities localized at the maxima of the gravitational potential. These singularities are steady-state counterparts of the so-called Zeldovich pancakes (ZP), interim wall-like structures appearing at nonlinear stages of development of gravitational instability. As the average matter density decreases, the period of the ZP lattice increases diverging at the density of Einstein's Universe. Solitary wave solutions are found at exactly the density of Einstein's Universe, and at a slightly larger density the wave may be viewed as a lattice of well-separated ZP solitons.


title: Hydrostatic mass estimates of massive galaxy clusters: a study with varying hydrodynamics flavours and non-thermal pressure support

abstract: We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but still producing clusters which broadly match observations. At the same mass resolution as BAHAMAS, for the most massive clusters ($M_{500} > 10^{15}$ M$_{\odot}$), we find changes in the SPH method produce the greatest differences in the final halo, while the subgrid models dominate at lower mass. By calculating the mass of all of the clusters using different permutations of the pressure, temperature and density profiles, created with either the true simulated data or mock spectroscopic data, we find that the spectroscopic temperature causes a bias in the hydrostatic mass estimates which increases with the mass of the cluster, regardless of the SPH flavour used. For the most massive clusters, the estimated mass of the cluster using spectroscopic density and temperature profiles is found to be as low as 50 per cent of the true mass compared to $\sim$ 90 per cent for low mass clusters. When including a correction for non-thermal pressure, the spectroscopic hydrostatic mass estimates are less biased on average and the mass dependence of the bias is reduced, although the scatter in the measurements does increase.


title: Photon generation via dynamical Casimir effect in an optomechanical cavity as a closed quantum system

abstract: We present an analytical and numerical analysis of the particle creation in an optomechanical cavity in parametric resonance. We treat both the electromagnetic field and the mirror as quantum degrees of freedom and study the dynamical evolution as a closed quantum system. We consider different initial states and investigate the spontaneous emission of photons from phonons in the mirror. We find that for initial phononic product states the evolution of the photon number can be described as a non-harmonic quantum oscillator, providing an useful tool so as to estimate the maximum and mean number of photons produced for arbitrary high energies. The efficiency of this mechanism is further analyzed for a detuned cavity as well as the possibility of stimulating the photon production by adding some initial ones to the cavity. We also find relationships for the maximum and mean entanglement between the mirror and the wall in these states. Additionally we study coherent states for the motion of the mirror to connect this model with previous results from quantum field theory with a classical mirror. Finally we study thermal states of phonons in the wall and the equilibration process that leads to a stationary distribution.


title: Cone-Adapted Shearlets and Radon Transforms

abstract: We show that the cone-adapted shearlet coefficients can be computed by means of the limited angle horizontal and vertical (affine) Radon transforms and the one-dimensional wavelet transform. This yields formulas that open new perspectives for the inversion of the Radon transform.


title: Quantitative comparison of Anderson impurity solvers applied to transport in quantum dots

abstract: We study the single impurity Anderson model (SIAM) using the equations of motion method (EOM), the non-crossing approximation (NCA), the one-crossing approximation (OCA), and Wilson's numerical renormalization group (NRG). We calculate the density of states and the linear conductance focusing on their dependence on the chemical potential and on the temperature paying special attention to the Kondo and Coulomb blockade regimes for a large range of model parameters. We report that some standard approximations based on the EOM technique display a rather unexpected poor behavior in the Coulomb blockade regime even at high temperatures. Our study offers a critical comparison between the different methods as well as a detailed compilation of the shortcomings and limitations due the approximations involved in each technique, thus allowing for a cost-benefit analysis of the different solvers that considers both numerical precision and computational performance.


title: Recovering the full dimensionality of hyperentanglement in collinear photon pairs

abstract: Exploiting hyperentanglement of photon pairs, that is, simultaneous entanglement in multiple degrees of freedom(DOFs), increases the dimensionality of Hilbert spaces for quantum information processing. However, generation of hyperentangled photon pairs collinearlly, while produces high brightness, results in a smaller Hilbert space due to the two photons being in the same spatial mode. In this letter, we point out that one can recover the full dimensionality of such hyperentanglement through a simple interference set up, similar to the time-reversed Hong-Ou-Mandel (TR-HOM) process. Different from the standard TR-HOM, we point out a critical phase condition has to be satisfied in order to recover the hyperentanglement. We theoretically analyze the realization of this approach and discuss the feasibility of generating truly hyperentangled photon pairs. Our proposed approach does not require post-selection and hence enables efficient hyper-entangled photon pairs generation for high-dimensional quantum applications.


title: Assembly of nanocube super-structures directed by surface and magnetic interactions

abstract: We model the stabilization of clusters and lattices of cuboidal particles with long-ranged magnetic dipolar and short-ranged surface interactions. Two realistic systems were considered: one with magnetization orientated in the [001] crystallographic direction, and the other with magnetization along the [111] direction. We have studied magnetic nanocubes clusters first in the limit of $T=0$~K intending to elucidate the structural genesis of low energy configurations and then analyzed finite-temperature behavior of the same systems in simulations. Our results demonstrate that dipolar coupling can stabilize nanoparticle assemblies with cubic, planar, and linear arrangements seen previously in experiments. While attractive surface energy supports the formation of super-cubes, the repulsion results in the elongated structures in the form of rods and chains. We observe the stabilization of the ferromagnetic planar arrangements of the cubes standing on their corners and in contact over edges. We illustrate that minimal energy structures depend only on the size of the assembly and balance of surface repulsion and magnetic dipolar coupling. The presented results are scalable to different particle sizes and material parameters.


title: Unsupervised Robust Disentangling of Latent Characteristics for Image Synthesis

abstract: Deep generative models come with the promise to learn an explainable representation for visual objects that allows image sampling, synthesis, and selective modification. The main challenge is to learn to properly model the independent latent characteristics of an object, especially its appearance and pose. We present a novel approach that learns disentangled representations of these characteristics and explains them individually. Training requires only pairs of images depicting the same object appearance, but no pose annotations. We propose an additional classifier that estimates the minimal amount of regularization required to enforce disentanglement. Thus both representations together can completely explain an image while being independent of each other. Previous methods based on adversarial approaches fail to enforce this independence, while methods based on variational approaches lead to uninformative representations. In experiments on diverse object categories, the approach successfully recombines pose and appearance to reconstruct and retarget novel synthesized images. We achieve significant improvements over state-of-the-art methods which utilize the same level of supervision, and reach performances comparable to those of pose-supervised approaches. However, we can handle the vast body of articulated object classes for which no pose models/annotations are available.


title: "Reply to A comment on A test of general relativity using the LARES and LAGEOS satellites and a GRACE Earth gravity model, by I. Ciufolini et al."

abstract: In 2016, we published "A test of general relativity using the LARES and LAGEOS satellites and a GRACE Earth gravity model. Measurement of Earth's dragging of inertial frames [1]", a measurement of frame-dragging, a fundamental prediction of Einstein's theory of General Relativity, using the laser-ranged satellites LARES, LAGEOS and LAGEOS 2. The formal error, or precision, of our test was about 0.2% of frame-dragging, whereas the systematic error was estimated to be about 5%. In the 2017 paper "A comment on "A test of general relativity using the LARES and LAGEOS satellites and a GRACE Earth gravity model by I. Ciufolini et al." "by L. Iorio [2] (called I2017 in the following), it was incorrectly claimed that, when comparing different Earth gravity field models, the systematic error in our test due to the Earth's even zonal harmonics of degree 6, 8, 10 could be aslarge as 15%, 6% and 36%, respectively. Furthermore, I2017 contains other, also incorrect, claims about the number of necessary significant decimal digits of the coefficients used in our test (claimed to be nine), in order to eliminate the largest uncertainties in the even zonals of degree 2 and 4, and about the non-repeatability of our test. Here we analyze and rebut those claims in I2017.


title: Objective Bayesian Analysis of a Cokriging Model for Hierarchical Multifidelity Codes

abstract: Autoregressive cokriging models have been widely used to emulate multiple computer models with different levels of fidelity. The dependence structures are modeled via Gaussian processes at each level of fidelity, where covariance structures are often parameterized up to a few parameters. The predictive distributions typically require intensive Monte Carlo approximations in previous works. This article derives new closed-form formulas to compute the means and variances of predictive distributions in autoregressive cokriging models that only depend on correlation parameters. For parameter estimation, we consider objective Bayesian analysis of such autoregressive cokriging models. We show that common choices of prior distributions, such as the constant prior and inverse correlation prior, typically lead to improper posteriors. We also develop several objective priors such as the independent reference prior and the independent Jeffreys prior that are shown to yield proper posterior distributions. This development is illustrated with a borehole function in an eight-dimensional input space and applied to an engineering application in a six-dimensional input space.


title: Joint analysis of EDGES $21$-cm line observations with standard candles and rulers in $\Lambda$CDM and non-adiabatic gCg models

abstract: A decomposed generalised Chaplygin gas (gCg) with energy flux from dark energy to dark matter, represented by a negative value for the gas parameter $\alpha$, is shown to alleviate the tension between EDGES data and the cosmological standard model. Using EDGES data and employing a Bayesian statistical analysis, the agreement with the standard model is only marginal. However, if $\alpha$ is negative enough the gCg fits remarkably well the data, even in combination with SNe Ia datasets. On the other hand, when the CMB and BAO acoustic scales are included the preferred value for $\alpha$ is near zero, implying that a small deviation from $\Lambda$CDM is predicted.


title: Localized states in global AdS

abstract: We construct both local states and scattering states with finite energy in global AdS by inserting properly regularized operators in the CFT of arbitrary conformal dimension $(\Delta)$ at an instant of time. We give the state fixed angular momentum $(\ell)$ by integrating the result over a sphere with appropriate spherical harmonics. The energy of the states and their angular resolution is computed with CFT operator methods and is independent of having an AdS interpretation. In the semiclassical limit of large conformal dimension operators, these correspond to single particles localized within subAdS scales with width $1/\sqrt{\Delta}$ in AdS units, whose subsequent evolution is controlled by bulk geodesics. Our construction allows us to place a particle in any desired geodesic. For radial geodesics, we show that the amplitude to produce the desired state can be thought of as a regularized tunneling amplitude from the boundary to the radial turning point of the radial geodesic, while for other geodesics we argue that the insertion is at the outermost radial turning point of the corresponding geodesic.


title: Modeling and simulation for varying die gap fluid extrusion

abstract: The analysis and the simulation of fluid extruded out of a die opening are of great importance in different industrial applications. Mathematical modeling and experimental studies of extrusion out of a constant die gap has been the subject of numerous publications in the literature. Extrusion Blow Molding is a polymer forming technique used to manufacture a wide range of large hollow parts such as fuel tanks used in the automotive industry. The die gap is varied during extrusion to adjust the manufactured product's final thickness and to compensate for uneven stretching during molding. In this work, mathematical modeling of varying die gap extrusion is introduced. Two models using different approaches are presented to model the time dependent effects of varying the die gap on the extrudate. In the first approach, physics based partial differential equations are used to simulate extrusion. The Navier Stokes equations are used as governing equations to simulate molten polymers. The governing equations are solved using Arbitrary Lagrangian Eulerian based Finite Element Method. A novel numerical method for defining the free surface of the extrudate is proposed. In addition, the effects of varying the die gap on the material velocity, pressure and extrudate shape are reported. In the second approach, a parameter identification scheme is used to obtain the parameters of a novel model with a predefined structure. The purpose of the proposed model is to replicate time dependent extrudate thickness of the physics based model.


title: Triple points and sign of circulation

abstract: Mach reflection generally produces a contact discontinuity whose circulation has previously been analyzed using "thermodynamic" arguments based on the Hugoniot relations across the shocks. We focus on "kinematic" techniques that avoid assumptions about the equation of state, using only jump relations for conservation of mass and momentum, but not energy. We give a new short proof for non-existence of pure (no contact) triple shocks, recovering a result of Serre. For MR with a zero-circulation but nonzero-density-jump contact we show that the incident shock must be normal. Nonexistence without contacts generalizes to two or more incident shocks if we assume all shocks are compressive. The sign of circulation across the contact has previously been controlled with entropy arguments, showing the post-Mach-stem velocity is generally smaller. We give a kinematic proof assuming compressive shocks and another condition, for example backward incident shocks, or a weak form of the Lax condition. We also show that for 2+2 interactions (two "upper" shocks with clockwise flow meeting two "lower" shocks with counterclockwise flow in a single point) the circulation sign can generally not be controlled. For $\gamma$-law pressure we show 2+2 interactions without contact must be either symmetric or antisymmetric, with symmetry favored at low Mach number and low shock strength. For full potential flow instead of the Euler equations we surprisingly find, contrary to folklore and prior results for other models, that pure triple shocks without contacts are possible, even for $\gamma$-law pressure with $1<\gamma<3$.


title: Coverage Analysis for Energy-Harvesting UAV-assisted mmWave Cellular Networks

abstract: In this paper, we jointly consider the downlink simultaneous wireless information and power transfer (SWIPT) and uplink information transmission in unmanned aerial vehicle (UAV)-assisted millimeter wave (mmWave) cellular networks, in which the user equipment (UE) locations are modeled using Poisson cluster processes (e.g., Thomas cluster processes or Mat\'ern cluster processes). Distinguishing features of mmWave communications, such as different path loss models for line-of-sight (LOS) and non-LOS (NLOS) links and directional transmissions are taken into account. In the downlink phase, the association probability, and energy coverages of different tier UAVs and ground base stations (GBSs) are investigated. Moreover, we define a successful transmission probability to jointly present the energy and signal-to-interference-plus-noise ratio (SINR) coverages and provide general expressions. In the uplink phase, we consider the scenario that each UAV receives information from its own cluster member UEs. We determine the Laplace transform of the interference components and characterize the uplink SINR coverage. In addition, we formulate the average uplink throughput, with the goal to identify the optimal time division multiplexing between the donwlink and uplink phases. Through numerical results we investigate the impact of key system parameters on the performance. We show that the network performance is improved when the cluster size becomes smaller. In addition, we analyze the optimal height of UAVs, optimal power splitting value and optimal time division multiplexing that maximizes the network performance.


title: Deep Learning at the Edge

abstract: The ever-increasing number of Internet of Things (IoT) devices has created a new computing paradigm, called edge computing, where most of the computations are performed at the edge devices, rather than on centralized servers. An edge device is an electronic device that provides connections to service providers and other edge devices; typically, such devices have limited resources. Since edge devices are resource-constrained, the task of launching algorithms, methods, and applications onto edge devices is considered to be a significant challenge. In this paper, we discuss one of the most widely used machine learning methods, namely, Deep Learning (DL) and offer a short survey on the recent approaches used to map DL onto the edge computing paradigm. We also provide relevant discussions about selected applications that would greatly benefit from DL at the edge.


title: Bottom-Up Meta-Policy Search

abstract: Despite of the recent progress in agents that learn through interaction, there are several challenges in terms of sample efficiency and generalization across unseen behaviors during training. To mitigate these problems, we propose and apply a first-order Meta-Learning algorithm called Bottom-Up Meta-Policy Search (BUMPS), which works with two-phase optimization procedure: firstly, in a meta-training phase, it distills few expert policies to create a meta-policy capable of generalizing knowledge to unseen tasks during training; secondly, it applies a fast adaptation strategy named Policy Filtering, which evaluates few policies sampled from the meta-policy distribution and selects which best solves the task. We conducted all experiments in the RoboCup 3D Soccer Simulation domain, in the context of kick motion learning. We show that, given our experimental setup, BUMPS works in scenarios where simple multi-task Reinforcement Learning does not. Finally, we performed experiments in a way to evaluate each component of the algorithm.


title: Flexible Bayesian modelling in dichotomous item response theory using mixtures of skewed item curves

abstract: Most Item Response Theory (IRT) models for dichotomous responses are based on probit or logit link functions which assume a symmetric relationship between the probability of a correct response and the latent traits of individuals submitted to a test. This assumption restricts the use of those models to the case in which all items have a symmetric behaviour. On the other hand, asymmetric models proposed in the literature impose that all the items in a test have an asymmetric behaviour. This assumption is inappropriate for great part of the tests which are, in general, composed by both symmetric and asymmetric items. Furthermore, a straightforward extension of the existing models in the literature would require a prior selection of the items' symmetry/asymmetry status. This paper proposes a Bayesian IRT model that accounts for symmetric and asymmetric items in a flexible though parsimonious way. That is achieved by assigning a finite mixture prior to the skewness parameter, with one of the mixture components being a point-mass at zero. This allows for analyses under both model selection and model averaging approaches. Asymmetric item curves are designed through the centred skew normal distribution, which has a particularly appealing parametrisation in terms of parameter interpretation and computational efficiency. An efficient MCMC algorithm is proposed to perform Bayesian inference and its performance is investigated in some simulated examples. Finally, the proposed methodology is applied to a data set from a large scale educational exam in Brazil.


title: The Bitlet Model: Defining a Litmus Test for the Bitwise Processing-in-Memory Paradigm

abstract: This paper describes an analytical modeling tool called Bitlet that can be used, in a parameterized fashion, to understand the affinity of workloads to processing-in-memory (PIM) as opposed to traditional computing. The tool uncovers interesting trade-offs between operation complexity (cycles required to perform an operation through PIM) and other key parameters, such as system memory bandwidth, data transfer size, the extent of data alignment, and effective memory capacity involved in PIM computations. Despite its simplicity, the model has already proven useful. In the future, we intend to extend and refine Bitlet to further increase its utility.


title: GCI detection from raw speech using a fully-convolutional network

abstract: Glottal Closure Instants (GCI) detection consists in automatically detecting temporal locations of most significant excitation of the vocal tract from the speech signal. It is used in many speech analysis and processing applications, and various algorithms have been proposed for this purpose. Recently, new approaches using convolutional neural networks have emerged, with encouraging results. Following this trend, we propose a simple approach that performs a mapping from the speech waveform to a target signal from which the GCIs are obtained by peak-picking. However, the ground truth GCIs used for training and evaluation are usually extracted from EGG signals, which are not perfectly reliable and often not available. To overcome this problem, we propose to train our network on high-quality synthetic speech with perfect ground truth. The performances of the proposed algorithm are compared with three other state-of-the-art approaches using publicly available datasets, and the impact of using controlled synthetic or real speech signals in the training stage is investigated. The experimental results demonstrate that the proposed method obtains similar or better results than other state-of-the-art algorithms and that using large synthetic datasets with many speakers offers a better generalization ability than using a smaller database of real speech and EGG signals.


title: Fourier Analysis, Computing, and Image Formation for Spotlight Synthetic Aperture Radar

abstract: This article is written to serve as an introduction and survey of imaging with synthetic aperture radar (SAR). The reader will benefit from having some familiarity with harmonic analysis, electromagnetic radiation, and inverse problems. After an overview of the SAR problem and some main concepts, the SAR imaging problem is contextualized in terms of classical harmonic analysis. Within this context, we consider partial Fourier sums of off-centered Fourier data and correspondingly the convolutional kernels resulting from conventional SAR image formation techniques. Following this, we revisit imaging of random complex signals from frequency data as in SAR, providing simpler derivations of some previous results and extending these ideas to the continuous setting. These concepts are tied in with the derived convolutional kernels, and it is deduced how good an image approximation is when it is obtained from only a small band of high frequency Fourier coefficients. Finally, regularization methods are presented to improve the quality of SAR images. Corresponding MATLAB software is made available for reproducibility of most figures and to facilitate further exploration of the methods presented here.