title: An Overview of Enhancing Distance Learning Through Augmented and Virtual Reality Technologies

abstract: Although distance learning presents a number of interesting educational advantages as compared to in-person instruction, it is not without its downsides. We first assess the educational challenges presented by distance learning as a whole and identify 4 main challenges that distance learning currently presents as compared to in-person instruction: the lack of social interaction, reduced student engagement and focus, reduced comprehension and information retention, and the lack of flexible and customizable instructor resources. After assessing each of these challenges in-depth, we examine how AR/VR technologies might serve to address each challenge along with their current shortcomings, and finally outline the further research that is required to fully understand the potential of AR/VR technologies as they apply to distance learning.


title: Sample path generation of the stochastic volatility CGMY process and its application to path-dependent option pricing

abstract: This paper proposes the sample path generation method for the stochastic volatility version of CGMY process. We present the Monte-Carlo method for European and American option pricing with the sample path generation and calibrate model parameters to the American style S\&P 100 index options market, using the least square regression method. Moreover, we discuss path-dependent options such as Asian and Barrier options.


title: New Algorithms for Computing Field of Vision over 2D Grids

abstract: The aim of this paper is to propose new algorithms for Field of Vision (FOV) computation which improve on existing work at high resolutions. FOV refers to the set of locations that are visible from a specific position in a scene of a computer game. We summarize existing algorithms for FOV computation, describe their limitations, and present new algorithms which aim to address these limitations. We first present an algorithm which makes use of spatial data structures in a way which is new for FOV calculation. We then present a novel technique which updates a previously calculated FOV, rather than re-calculating an FOV from scratch. We compare our algorithms to existing FOV algorithms and show they provide substantial improvements to running time. Our algorithms provide the largest improvement over existing FOV algorithms at large grid sizes, thus allowing the possibility of the design of high resolution FOV-based video games.


title: FDApy: a Python package for functional data

abstract: We introduce the Python package, FDApy, as an implementation of functional data. This package provide modules for the analysis of such data. It includes classes for different dimensional data as well as irregularly sampled functional data. A simulation toolbox is also provided. It might be used to simulate different clusters of functional data. Some methodologies to handle these data are implemented, such as dimension reduction and clustering. New methods can be easily added. The package is publicly available on the Python Package Index and Github.


title: Photon pair distributions and event-mixing at low invariant mass

abstract: Neutral meson cross-section and direct photon HBT measurements in high energy pp, pA, and AA collisions rely on either subtracting or dividing out the background underlying the $M_\mathrm{inv}$ or $Q_\mathrm{inv}$ distributions. In this paper we investigate the photon pair distributions for simulated proton-proton collisions at $\sqrt{s}=13~$TeV, and study why the background in $M_\mathrm{inv}$ gets poorly described by the conventional event-mixing method implemented in experimental analyses. We show that it fails to describe the background as it gives a qualitatively different shape especially at invariant masses around the $\pi^{0}$ and $\eta$ meson mass, which is mainly caused by the correlations originating from parton fragmentation.


title: The Epoch of Giant Planet Migration Planet Search Program. I. Near-Infrared Radial Velocity Jitter of Young Sun-like Stars

abstract: We present early results from the Epoch of Giant Planet Migration program, a precise RV survey of over one hundred intermediate-age ($\sim$20$-$200 Myr) G and K dwarfs with the Habitable-Zone Planet Finder spectrograph (HPF) at McDonald Observatory's Hobby-Eberly Telescope (HET). The goals of this program are to determine the timescale and dominant physical mechanism of giant planet migration interior to the water ice line of Sun-like stars. Here, we summarize results from the first 14 months of this program, with a focus on our custom RV pipeline for HPF, a measurement of the intrinsic near-infrared RV activity of young Solar analogs, and modeling the underlying population-level distribution of stellar jitter. We demonstrate on-sky stability at the sub-2 m s$^{-1}$ level for the K2 standard HD 3765 using a least-squares matching method to extract precise RVs. Based on a subsample of 29 stars with at least three RV measurements from our program, we find a median RMS level of 34 m s$^{-1}$. This is nearly a factor of 2 lower than the median RMS level in the optical of 60 m s$^{-1}$ for a comparison sample with similar ages and spectral types as our targets. The observed near-infrared jitter measurements for this subsample are well reproduced with a log-normal parent distribution with $\mu=4.15$ and $\sigma=1.02$. Finally, by compiling RMS values from previous planet search programs, we show that near-infrared jitter for G and K dwarfs generally decays with age in a similar fashion to optical wavelengths, albeit with a shallower slope and lower overall values for ages $\lesssim$1 Gyr.


title: The GADOT Galaxy Survey: Dense Gas and Feedback in Herschel-Selected Starburst Galaxies at Redshifts 2 to 6

abstract: We report the detection of 23 OH+(1-0) absorption, emission, or P-Cygni-shaped lines and CO(9-8) emission lines in 18 Herschel-selected z=2-6 starburst galaxies with ALMA and NOEMA, taken as part of the Gas And Dust Over cosmic Time (GADOT) Galaxy Survey. We find that the CO(9-8) luminosity is higher than expected based on the far-infrared luminosity when compared to nearby star-forming galaxies. Together with the strength of the OH+ emission components, this may suggest that shock excitation of warm, dense molecular gas is more prevalent in distant massive dusty starbursts than in nearby star-forming galaxies on average, perhaps due to an impact of galactic winds on the gas. OH+ absorption is found to be ubiquitous in massive high-redshift starbursts, and is detected toward 89% of the sample. The majority of the sample shows evidence for outflows or inflows based on the velocity shifts of the OH+ absorption/emission, with a comparable occurrence rate of both at the resolution of our observations. A small subsample appears to show outflow velocities in excess of their escape velocities. Thus, starburst-driven feedback appears to be important in the evolution of massive galaxies in their most active phases. We find a correlation between the OH+ absorption optical depth and the dust temperature, which may suggest that warmer starbursts are more compact and have higher cosmic ray energy densities, leading to more efficient OH+ ion production. This is in agreement with a picture in which these high-redshift galaxies are "scaled-up" versions of the most intense nearby starbursts.


title: Observing Invisible Axions with Gravitational Waves

abstract: If the Peccei-Quinn symmetry associated to an axion has ever been restored after inflation, axion strings inevitably produce a contribution to the stochastic gravitational wave background. Combining effective field theory analysis with numerical simulations, we show that the resulting gravitational wave spectrum has logarithmic deviations from a scale invariant form with an amplitude that is significantly enhanced at low frequencies. As a result, a single ultralight axion-like particle with a decay constant larger than $10^{14}~{\rm GeV}$ and any mass between $10^{-18}~{\rm eV}$ and $10^{-28}~{\rm eV}$ leads to an observable gravitational wave spectrum and is compatible with constraints on the post-inflationary scenario from dark matter overproduction, isocurvature and dark radiation. Since the spectrum extends over a wide range of frequencies, the resulting signal could be detected by multiple experiments. We describe straightforward ways in which the Peccei-Quinn symmetry can be restored after inflation for such decay constants. We also comment on the recent possible NANOgrav signal in light of our results.


title: Spillway Preheating

abstract: In traditional models only an order one fraction of energy is transferred from the inflaton to radiation through nonperturbative resonance production in preheating immediately after inflation, due to backreaction effects. We propose a particle production mechanism that could improve the depletion of the inflaton energy density by up to four orders of magnitude. The improvement comes from the fast perturbative decays of resonantly produced daughter particles. They act as a "spillway" to drain these daughter particles, reducing their backreaction on the inflaton and keeping the resonant production effective for a longer period. Thus we dub the scenario "spillway preheating". We also show that the fraction of energy density remaining in the inflaton has a simple inverse power-law scaling in the scenario. In general, spillway preheating is a much more efficient energy dissipation mechanism, which may have other applications in model building for particle physics.


title: The Perseus ALMA Chemistry Survey (PEACHES). I. The Complex Organic Molecules in Perseus Embedded Protostars

abstract: To date, about two dozen low-mass embedded protostars exhibit rich spectra with lines of complex organic molecule (COM). These protostars seem to possess different enrichment in COMs. However, the statistics of COM abundance in low-mass protostars are limited by the scarcity of observations. This study introduces the Perseus ALMA Chemistry Survey (PEACHES), which aims at unbiasedly characterizing the chemistry of COMs toward the embedded (Class 0/I) protostars in the Perseus molecular cloud. Of 50 embedded protostars surveyed, 58% of them have emission from COMs. A 56%, 32%, and 40% of the protostars have CH$_3$OH, CH$_3$OCHO, and N-bearing COMs, respectively. The detectability of COMs depends neither on the averaged continuum brightness temperature, a proxy of the H$_2$ column density, nor on the bolometric luminosity and the bolometric temperature. For the protostars with detected COMs, CH$_3$OH has a tight correlation with CH$_3$CN, spanning more than two orders of magnitude in column densities normalized by the continuum brightness temperature, suggesting a chemical relation between CH$_3$OH and CH$_3$CN and a large chemical diversity in the PEACHES samples at the same time. A similar trend with more scatter is also found between all identified COMs, hinting at a common chemistry for the sources with COMs. The correlation between COMs is insensitive to the protostellar properties, such as the bolometric luminosity and the bolometric temperature. The abundance of larger COMs (CH$_3$OCHO and CH$_3$OCH$_3$) relative to that of smaller COMs (CH$_3$OH and CH$_3$CN) increases with the inferred gas column density, hinting at an efficient production of complex species in denser envelopes.


title: Thermodynamics of $d$-dimensional Schwarzschild black holes in the canonical ensemble

abstract: We study the thermodynamics of a $d$-dimensional Schwarzschild black hole in the canonical ensemble. This generalizes York's formalism to any number $d$ of dimensions. The canonical ensemble, characterized by a cavity of fixed radius $r$ and fixed temperature $T$ at the boundary, allows for two possible solutions in thermal equilibrium, a small and a large black hole. From the Euclidean action and the path integral approach, we obtain the free energy, the thermodynamic energy, the pressure, and the entropy, of the black hole plus cavity system. The entropy is given by the Bekenstein-Hawking area law. The heat capacity shows that the smaller black hole is in unstable equilibrium and the larger is stable. The photon sphere radius divides the stability criterion. To study perturbations, a generalized free energy function is obtained that allows to understand the possible phase transitions between classical hot flat space and the black holes. The Buchdahl radius, that appears naturally in the general relativistic study of star structure, also shows up in our context, the free energy is zero when the cavity's radius has the Buchdahl radius value. Then, if the cavity's radius is smaller than the Buchdahl radius classical hot flat space can nucleate a black hole. It is also pointed out the link between the canonical analysis performed and the direct perturbation of the path integral. Since gravitational hot flat space is a quantum system made purely of gravitons it is of interest to compare the free energies of quantum hot flat space and the stable black hole to find for which ranges of $r$ and $T$ one phase predominates over the other. Phase diagrams are displayed. The density of states at a given energy is found. Further calculations and comments are carried out, notably, a connection to thin shells in $d$ spacetime dimensions which are systems that are also apt to rigorous thermodynamics.


title: The Phantom of RAMSES user guide for galaxy simulations using Milgromian and Newtonian gravity

abstract: This document describes the general process of setting up, running, and analysing disc galaxy simulations using the freely available program Phantom of RAMSES (PoR). This implements Milgromian Dynamics (MOND) with a patch to the RAMSES grid-based $N$-body and hydrodynamical code that uses adaptive mesh refinement. We discuss the procedure of setting up isolated and interacting disc galaxy initial conditions for PoR, running the simulations, and analysing the results. This manual also concisely documents all previously developed MOND simulation codes and the results obtained with them.


title: Progenitor Mass Distribution for 22 Historic Core-Collapse Supernovae

abstract: We infer the progenitor mass distribution for 22 historic core-collapse supernovae (CCSNe) using a Bayesian hierarchical model. For this inference, we use the local star formation histories to estimate the age for each supernova (SN). These star formation histories often show multiple bursts of star formation; our model assumes that one burst is associated with the SN progenitor and the others are random bursts of star formation. The primary inference is the progenitor age distribution. Due to the limited number of historic SNe and highly uncertain star formation at young ages, we restrict our inference to the slope of the age distribution and the maximum age for CCSNe. Using single-star evolutionary models, we transform the progenitor age distribution into a progenitor mass distribution. Under these assumptions, the minimum mass for CCSNe is ${M_\textrm{min}}~=~8.60^{+0.37}_{-0.41} \ M_\odot$ and the slope of the progenitor mass distribution is $\alpha = -2.61^{+1.05}_{-1.18}$. The power-law slope for the progenitor mass distribution is consistent with the standard Salpeter initial mass function ($\alpha = -2.35$). These values are consistent with previous estimates using precursor imaging and the age-dating technique, further confirming that using stellar populations around SN and supernova remnants is a reliable way to infer the progenitor masses.


title: Hierarchical fragmentation in high redshift galaxies revealed by hydrodynamical simulations

abstract: High-redshift star-forming galaxies have very different morphologies compared to nearby ones. Indeed, they are often dominated by bright star-forming structures of masses up to $10^{8-9}$ $\mathrm{M}_\odot$ dubbed {\guillemotleft}giant clumps{\guillemotright}. However, recent observations questioned this result by showing only low-mass structures or no structure at all. We use Adaptative Mesh Refinement hydrodynamical simulations of galaxies with parsec-scale resolution to study the formation of structures inside clumpy high-redshift galaxies. We show that in very gas-rich galaxies star formation occurs in small gas clusters with masses below $10^{7-8}$ $\mathrm{M}_\odot$ that are themselves located inside giant complexes with masses up to $10^8$ and sometimes $10^9$ $\mathrm{M}_\odot$. Those massive structures are similar in mass and size to the giant clumps observed in imaging surveys, in particular with the Hubble Space Telescope. Using mock observations of simulated galaxies, we show that at very high resolution with instruments like the Atacama Large Millimeter Array or through gravitational lensing, only low-mass structures are likely to be detected, and their gathering into giant complexes might be missed. This leads to the non-detection of the giant clumps and therefore introduces a bias in the detection of these structures. We show that the simulated giant clumps can be gravitationally bound even when undetected in mocks representative for ALMA observations and HST observations of lensed galaxies. We then compare the top-down fragmentation of an initially warm disc and the bottom-up fragmentation of an initially cold disc to show that the process of formation of the clumps does not impact their physical properties.


title: The cosmology dependence of galaxy clustering and lensing from a hybrid $N$-body-perturbation theory model

abstract: We implement a model for the two-point statistics of biased tracers that combines dark matter dynamics from $N$-body simulations with an analytic Lagrangian bias expansion. Using Aemulus, a suite of $N$-body simulations built for emulation of cosmological observables, we emulate the cosmology dependence of these nonlinear spectra from redshifts $z = 0$ to $z=2$. We quantify the accuracy of our emulation procedure, which is sub-per cent at $k=1\, h {\rm Mpc}^{-1}$ for the redshifts probed by upcoming surveys and improves at higher redshifts. We demonstrate its ability to describe the statistics of complex tracer samples, including those with assembly bias and baryonic effects, reliably fitting the clustering and lensing statistics of such samples at redshift $z\simeq 0.4$ to scales of $k_{\rm max} \approx 0.6\, h\mathrm{Mpc}^{-1}$. We show that the emulator can be used for unbiased cosmological parameter inference in simulated joint clustering and galaxy--galaxy lensing analyses with data drawn from an independent $N$-body simulation. These results indicate that our emulator is a promising tool that can be readily applied to the analysis of current and upcoming datasets from galaxy surveys.


title: Anomalies in the gravitational recoil of eccentric black-hole mergers with unequal mass ratios

abstract: The radiation of linear momentum imparts a recoil (or "kick") to the center of mass of a merging black hole binary system. Recent numerical relativity calculations have shown that eccentricity can lead to an approximate 25% increase in recoil velocities for equal-mass, spinning binaries with spins lying in the orbital plane ("superkick" configurations) [U Sperhake et al. Phys. Rev. D 101 (2020) 024044 (arXiv:1910.01598)]. Here we investigate the impact of nonzero eccentricity on the kick magnitude and gravitational-wave emission of nonspinning, unequal-mass black hole binaries. We confirm that nonzero eccentricities at merger can lead to kicks which are larger by up to ~25% relative to the quasicircular case. We also find that the kick velocity $v$ has an oscillatory dependence on eccentricity, that we interpret as a consequence of changes in the angle between the infall direction at merger and the apoapsis (or periapsis) direction.


title: Explaining Excess Dipole in NVSS Data Using Superhorizon Perturbation

abstract: Many observations in recent times have shown evidence against the standard assumption of isotropy in the Big Bang model. Introducing a superhorizon scalar metric perturbation has been able to explain some of these anomalies. In this work, we probe the net velocity arising due to the perturbation. We find that this extra component does not contribute to the CMB dipole amplitude while it does contribute to the dipole in large scale structures. Thus, within this model's framework, our velocity with respect to the large scale structure is not the same as that extracted from the CMB dipole, assuming it to be of purely kinematic origin. Taking this extra velocity component into account, we study the superhorizon mode's implications for the excess dipole observed in the NRAO VLA Sky Survey (NVSS). We find that the mode can consistently explain both the CMB and NVSS observations. We also find that the model leads to small contributions to the local bulk flow and the dipole in Hubble parameter, which are consistent with observations. The model leads to several predictions which can be tested in future surveys. In particular, it implies that the observed dipole in large scale structure should be redshift dependent and should show an increase in amplitude with redshift. We also find that the Hubble parameter should show a dipole anisotropy whose amplitude must increase with redshift in the CMB frame. Similar anisotropic behaviour is expected for the observed redshift as a function of the luminosity distance.


title: Measuring time delays: I. Using a flux time series that is a linear combination of time-shifted light curves

abstract: (Abridged) Several phenomena in astrophysics generate light curves with time delays. Among these are reverberation mapping, and lensed quasars. In some systems, the measurement of the time-delay is complicated by the fact that the delayed components are unresolved and that the light curves are generated from a red-noise process. We derive the likelihood function of the observations given a model of either a combination of time-delayed light curves or a single light curve. This likelihood function is different from the auto-correlation function. We demonstrate that given a single-band light curve that is a combination of two (or more) time-shifted copies of an original light curve, generated from a red-noise probability distribution, we can test if the total-flux light curve is a composition of time-delayed copies or, alternatively, is consistent with being the original light curve. Furthermorew, in some realistic cases, it is possible to measure the time delays and flux ratios between these unresolved components even when the flux ratio is about 1/10. This method is useful for identifying lensed quasars and simultaneously measuring their time delays, and for estimating the reverberation time scales of active galactic nuclei. In a companion paper, we derive a method that uses the center-of-light position (e.g., of a lensed quasar) along with the combined flux. This allow us to identify lensed quasars and supernovae and measure their time delays, with higher fidelity compared to the flux-only method. The astrometry + flux method, however, is not suitable for quasar reverberation mapping. We also comment on the commonly used method of fitting a power-law model to a power spectrum, and present the proper likelihood function for such a fit. We test the new method on simulations and provide Python and MATLAB implementations.


title: Accurate and Efficient Simulations of Hamiltonian Mechanical Systems with Discontinuous Potentials

abstract: This article considers Hamiltonian mechanical systems with potential functions admitting jump discontinuities. The focus is on accurate and efficient numerical approximations of their solutions, which will be defined via the laws of reflection and refraction. Despite of the success of symplectic integrators for smooth mechanical systems, their construction for the discontinuous ones is nontrivial, and numerical convergence order can be impaired too. Several rather-usable numerical methods are proposed, including: a first-order symplectic integrator for general problems, a third-order symplectic integrator for problems with only one linear interface, arbitrarily high-order reversible integrators for general problems (no longer symplectic), and an adaptive time-stepping version of the previous high-order method. Interestingly, whether symplecticity leads to favorable long time performance is no longer clear due to discontinuity, as traditional Hamiltonian backward error analysis does not apply any more. Therefore, at this stage, our recommended default method is the last one. Various numerical evidence, on the order of convergence, long time performance, momentum map conservation, and consistency with the computationally-expensive penalty method, are supplied. A complex problem, namely the Sauteed Mushroom, is also proposed and numerically investigated, for which multiple bifurcations between trapped and ergodic dynamics are observed.

categories: math.NA cs.CE cs.NA math.DS physics.comp-ph


title: Twin Higgs Portal Dark Matter

abstract: Many minimal models of dark matter (DM) or canonical solutions to the hierarchy problem are either excluded or severely constrained by LHC and direct detection null results. In particular, Higgs Portal Dark Matter (HPDM) features a scalar coupling to the Higgs via a quartic interaction, and obtaining the measured relic density via thermal freeze-out gives definite direct detection predictions which are now almost entirely excluded. The Twin Higgs solves the little hierarchy problem without coloured top partners by introducing a twin sector related to the Standard Model (SM) by a discrete symmetry. We generalize HPDM to arbitrary Twin Higgs models and introduce Twin Higgs Portal Dark Matter (THPDM), which features a DM candidate with an $SU(4)$-invariant quartic coupling to the Twin Higgs scalar sector. Given the size of quadratic corrections to the DM mass, its most motivated scale is near the mass of the radial mode. In that case, DM annihilation proceeds with the full Twin Higgs portal coupling, while direct detection is suppressed by the pNGB nature of the 125 GeV Higgs. For a standard cosmological history, this results in a predicted direct detection signal for THPDM that is orders of magnitude below that of HPDM with very little dependence on the precise details of the twin sector, evading current bounds but predicting possible signals at next generation experiments. In many Twin Higgs models, twin radiation contributions to $\Delta N_\mathrm{eff}$ are suppressed by an asymmetric reheating mechanism. We study this by extending the $\nu$MTH and $X$MTH models to include THPDM and compute the viable parameter space according to the latest CMB bounds. The injected entropy dilutes the DM abundance as well, resulting in additional suppression of direct detection below the neutrino floor.


title: Supervised quantum machine learning models are kernel methods

abstract: With near-term quantum devices available and the race for fault-tolerant quantum computers in full swing, researchers became interested in the question of what happens if we replace a supervised machine learning model with a quantum circuit. While such "quantum models" are sometimes called "quantum neural networks", it has been repeatedly noted that their mathematical structure is actually much more closely related to kernel methods: they analyse data in high-dimensional Hilbert spaces to which we only have access through inner products revealed by measurements. This technical manuscript summarises and extends the idea of systematically rephrasing supervised quantum models as a kernel method. With this, a lot of near-term and fault-tolerant quantum models can be replaced by a general support vector machine whose kernel computes distances between data-encoding quantum states. Kernel-based training is then guaranteed to find better or equally good quantum models than variational circuit training. Overall, the kernel perspective of quantum machine learning tells us that the way that data is encoded into quantum states is the main ingredient that can potentially set quantum models apart from classical machine learning models.


title: Explaining the scatter in the galaxy mass-metallicity relation with gas flows

abstract: The physical origin of the scatter in the relation between galaxy stellar mass and the metallicity of the interstellar medium, i.e. the Mass-Metallicity Relation (MZR), reflects the relative importance of key processes in galaxy evolution. The \eagle cosmological hydrodynamical simulation is used to investigate the correlations between the residuals of the MZR and the residuals of the relations between stellar mass and, respectively, specific inflow, outflow and star formation rate as well as the gas fraction for central galaxies. At low redshift, all these residuals are found to be anti-correlated with the residuals of the MZR for $M_\star/\mathrm{M}_\odot \lesssim 10^{10}$. The correlations between the residuals of the MZR and the residuals of the other relations with mass are interrelated, but we find that gas fraction, specific inflow rate and specific outflow rate all have at least some independent influence on the scatter of the MZR. We find that, while for $M_\star/\mathrm{M}_\odot > 10^{10.4}$ the specific mass of the nuclear black hole is most important, for $M_\star/\mathrm{M}_\odot \lesssim 10^{10.3}$ gas fraction and specific inflow rate are the variables that correlate most strongly with the MZR scatter. The timescales involved in the residual correlations and the time that galaxies stay above the MZR are revealed to be a few Gyr. However, most galaxies that are below the MZR at $z=0$ have been below the MZR throughout their lifetimes.


title: SDSS-IV MaNGA: the "G-dwarf problem" revisited

abstract: The levels of heavy elements in stars are the product of enhancement by previous stellar generations, and the distribution of this metallicity among the population contains clues to the process by which a galaxy formed. Most famously, the "G-dwarf problem" highlighted the small number of low-metallicity G-dwarf stars in the Milky Way, which is inconsistent with the simplest picture of a galaxy formed from a "closed box" of gas. It can be resolved by treating the Galaxy as an open system that accretes gas throughout its life. This observation has classically only been made in the Milky Way, but the availability of high-quality spectral data from SDSS-IV MaNGA and the development of new analysis techniques mean that we can now make equivalent measurements for a large sample of spiral galaxies. Our analysis shows that high-mass spirals generically show a similar deficit of low-metallicity stars, implying that the Milky Way's history of gas accretion is common. By contrast, low-mass spirals show little sign of a G-dwarf problem, presenting the metallicity distribution that would be expected if such systems evolved as pretty much closed boxes. This distinction can be understood from the differing timescales for star formation in galaxies of differing masses.


title: On formal concepts of random formal contexts

abstract: In formal concept analysis, it is well-known that the number of formal concepts can be exponential in the worst case. To analyze the average case, we introduce a probabilistic model for random formal contexts and prove that the average number of formal concepts has a superpolynomial asymptotic lower bound.


title: Measuring time delays: II. Using observations of the unresolved flux and astrometry

abstract: Lensed quasars and supernovae can be used to study galaxies' gravitational potential and measure cosmological parameters. The typical image separation of objects lensed by galaxies is of the order of 0.5". Therefore, finding the ones with small separations, and measuring their time-delays using ground-based observations is challenging. We suggest a new method to identify lensed quasars and simultaneously measure their time-delays, using seeing-limited synoptic observations in which the lensed quasar images and the lensing galaxy are unresolved. We show that using the light curve of the combined flux, and the astrometric measurements of the center-of-light position of the lensed images, the lensed nature of a quasar can be identified, and its time-delay can be measured. We provide the analytic formalism to do so, taking into account the measurement errors and the fact that the power spectra of quasar light curves is red (i.e., the light curve is highly correlated). We demonstrate our method on simulated data, while its implementation to real data will be presented in future papers. Our simulations suggest that, under reasonable assumptions, the new method can detect unresolved lensed quasars and measure their time delays, even when the image separation is below 0.1", or the flux ratio between the faintest and brightest images is as low as 0.03. Python and MATLAB implementations are provided. In a companion paper, we present a method for measuring the time delay using the combined flux observations. Although the flux-only method is less powerful, it may be useful in cases in which the astrometric information is not relevant (e.g., reverberation mapping).


title: Probing core overshooting using subgiant asteroseismology: the case of KIC10273246

abstract: The size of convective cores remains uncertain, despite its substantial influence on stellar evolution, and thus on stellar ages. The seismic modeling of young subgiants can be used to obtain indirect constraints on the core structure during main sequence, thanks to the high probing potential of mixed modes. We selected the young subgiant KIC10273246, observed by Kepler, based on its mixed-mode properties. We thoroughly modeled this star, with the aim of placing constraints on the size of its main sequence convective core. We first extracted the parameters of the oscillation modes of the star using the full Kepler data set. To overcome the challenges posed by the seismic modeling of subgiants, we proposed a method which is specifically tailored for subgiants with mixed modes and consists in a nested optimization. We then applied this method to perform a detailed seismic modeling of KIC10273246. We obtained models that show good statistical agreements with the observations, both seismic and non-seismic. We showed that including core overshooting in the models significantly improves the quality of the seismic fit, optimal models being found for $\alpha_{\mathrm{ov}} = 0.15$. Higher amounts of core overshooting strongly worsen the agreement with the observations and are thus firmly ruled out. We also found that having access to two g-dominated mixed modes in young subgiants allows us to place stronger constraints on the gradient of molecular weight in the core and on the central density. This study confirms the high potential of young subgiants with mixed modes to investigate the size of main-sequence convective cores. It paves the way for a more general study including the subgiants observed with Kepler, TESS, and eventually PLATO.


title: The Highest Energy HAWC Sources are Likely Leptonic and Powered by Pulsars

abstract: The HAWC Collaboration has observed gamma rays at energies above 56 TeV from a collection of nine sources. It has been suggested that this emission could be hadronic in nature, requiring that these systems accelerate cosmic-ray protons or nuclei up to PeV-scale energies. In this paper, we instead show that the spectra of these objects favor a leptonic (inverse Compton) origin for their emission. More specifically, the gamma-ray emission from these objects can be straightforwardly accommodated within a model in which $\sim \mathcal{O}(10\%)$ of the host pulsar's spindown power is transferred into the acceleration of electrons and positrons with a power-law spectrum that extends to several hundred TeV or higher. The spectral break that is observed among these sources is naturally explained within the context of this simple model, and occurs at the energy where the timescale for energy losses matches the age of the pulsar. In contrast, this spectral feature cannot be straightforwardly accommodated in hadronic scenarios. Furthermore, hadronic models predict that these sources should produce more emission at GeV-scale energies than is observed. In light of these considerations, we conclude that HAWC's highest energy sources should be interpreted as TeV halos or pulsar wind nebulae, which produce their emission through inverse Compton scattering, and are powered by the rotational kinetic energy of their host pulsar.


title: Multimessenger constraints on the dark matter interpretation of the Fermi-LAT Galactic center excess

abstract: An excess of $\gamma$ rays in the data measured by the Fermi Large Area Telescope in the direction of the Galactic center has been reported in several publications. This excess, labeled as the Galactic center excess (GCE), is detected analyzing the data with different interstellar emission models, point source catalogs and analysis techniques. The characteristics of the GCE, recently measured with unprecedented precision, are all compatible with dark matter particles (DM) annihilating in the main halo of our Galaxy, even if other interpretations are still not excluded. We investigate the DM candidates that fit the observed GCE spectrum and spatial morphology. We assume a simple scenario with DM annihilating into a single channel but we inspect also more complicated models with two and three channels. We perform a search for a $\gamma$-ray flux from a list of 48 Milky Way dwarf spheroidal galaxies (dSphs) using state-of-the-art estimation of the DM density in these objects. Since we do not find any significant signal from the dSphs, we put upper limits on the annihilation cross section that result to be compatible with the DM candidate that fits the GCE. However, we find that the GCE DM signal is excluded by the AMS-02 $\bar{p}$ flux data for all hadronic and semi-hadronic annihilation channels unless the vertical size of the diffusion halo is smaller than 2 kpc -- which is in tension with radioactive cosmic ray fluxes and radio data. Furthermore, AMS-02 $e^+$ data rule out pure or mixed channels with a component of $e^+ e^-$. The only DM candidate that fits the GCE spectrum and is compatible with constraints obtained with the combined dSphs analysis and the AMS-02 $\bar{p}$ and $e^+$ data annihilates purely into $\mu^+\mu^-$, has a mass of 60 GeV and roughly a thermal cross section.


title: Outflow from super-Eddington flow: where it originates from and how much impact it gives?

abstract: It is widely believed that super-Eddington accretion flow can produce powerful outflow, but where it originates from and how much mass and energy are carried away to which directions? To answer to these questions, we newly perform a large-box, two-dimensional radiation hydrodynamic simulation, paying special attention lest the results should depend on adopted initial and boundary conditions. We could achieve a quasi-steady state in an unprecedentedly large range, $r=2~r_{\rm S}$-$600~r_{\rm S}$ (with $r_{\rm S}$ being the Schwarzschild radius) from the black hole. The accretion rate onto the central $10 ~M_{\odot}$ black hole is $\dot{M}_{\rm BH} \sim 180 ~L_{\rm Edd}/c^{2}$, whereas the mass outflow rate is ${\dot M}_{\rm outflow} \sim 24 ~L_{\rm Edd}/c^2$ (where $L_{\rm Edd}$ and $c$ are the Eddington luminosity and the speed of light, respectively). The ratio (${\dot M}_{\rm outflow}/{\dot M}_{\rm BH} \sim 0.14$) is much less than those reported previously. By careful inspection we find that most of outflowing gas which reach the outer boundary originates from the region at $R\lesssim140~r_{\rm S}$, while gas at $140~r_{\rm S}$-$230 ~r_{\rm S}$ forms failed outflow. Therefore, significant outflow occurs inside the trapping radius $\sim 450 ~r_{\rm S}$. The mechanical energy flux (or mass flux) reaches its maximum in the direction of $\sim 15^\circ$ ($\sim 80^\circ$) from the rotation axis. The total mechanical luminosity is $L_{\rm mec}\sim 0.16~L_{\rm Edd}$, while the isotropic X-ray luminosity varies from $L_{\rm X}^{\rm ISO}\sim 2.9~L_{\rm Edd}$, (for a face-on observer) to $\sim 2.1~L_{\rm Edd}$ (for a nearly edge-on observer). The power ratio is $L_{\rm mec}/L_{\rm X}^{\rm ISO}\sim 0.05$-$0.08$, in good agreement with the observations of Ultra-Luminous X-ray sources surrounded by optical nebulae.


title: Holographic entanglement negativity and replica symmetry breaking

abstract: Since the work of Ryu and Takayanagi, deep connections between quantum entanglement and spacetime geometry have been revealed. The negative eigenvalues of the partial transpose of a bipartite density operator is a useful diagnostic of entanglement. In this paper, we discuss the properties of the associated entanglement negativity and its R\'enyi generalizations in holographic duality. We first review the definition of the R\'enyi negativities, which contain the familiar logarithmic negativity as a special case. We then study these quantities in the random tensor network model and rigorously derive their large bond dimension asymptotics. Finally, we study entanglement negativity in holographic theories with a gravity dual, where we find that R\'enyi negativities are often dominated by bulk solutions that break the replica symmetry. From these replica symmetry breaking solutions, we derive general expressions for R\'enyi negativities and their special limits including the logarithmic negativity. In fixed-area states, these general expressions simplify dramatically and agree precisely with our results in the random tensor network model. This provides a concrete setting for further studying the implications of replica symmetry breaking in holography.


title: Enabling Dataflow Optimization for Quantum Programs

abstract: We propose an IR for quantum computing that directly exposes quantum and classical data dependencies for the purpose of optimization. The Quantum Intermediate Representation for Optimization (QIRO) consists of two dialects, one input dialect and one that is specifically tailored to enable quantum-classical co-optimization. While the first employs a perhaps more intuitive memory-semantics (quantum operations act as side-effects), the latter uses value-semantics (operations consume and produce states). Crucially, this encodes the dataflow directly in the IR, allowing for a host of optimizations that leverage dataflow analysis. We discuss how to map existing quantum programming languages to the input dialect and how to lower the resulting IR to the optimization dialect. We present a prototype implementation based on MLIR that includes several quantum-specific optimization passes. Our benchmarks show that significant improvements in resource requirements are possible even through static optimization. In contrast to circuit optimization at run time, this is achieved while incurring only a small constant overhead in compilation time, making this a compelling approach for quantum program optimization at application scale.


title: Is a Miracle-less WIMP Ruled Out?

abstract: We examine a real electroweak triplet scalar field as dark matter, abandoning the requirement that its relic abundance is determined through freeze out in a standard cosmological history (a situation which we refer to as a 'miracle-less WIMP'). We extract the bounds on such a particle from collider searches, searches for direct scattering with terrestrial targets, and searches for the indirect products of annihilation. Each type of search provides complementary information, and each is most effective in a different region of parameter space. LHC searches tend to be highly dependent on the mass of the SU(2) charged partner state, and are effective for very large or very tiny mass splitting between it and the neutral dark matter component. Direct searches are very effective at bounding the Higgs portal coupling, but ineffective once it falls below $\lambda_{\text{eff}} \lesssim 10^{-3}$. Indirect searches suffer from large astrophysical uncertainties due to the backgrounds and $J$-factors, but do provide key information for $\sim$ 100 GeV to TeV masses. Synthesizing the allowed parameter space, this example of WIMP dark matter remains viable, but only in miracle-less regimes.


title: The original Wigner's friend paradox within a realist toy model

abstract: The original Wigner's friend paradox is a gedankenexperiment involving an observer described by an external agent. The paradox highlights the tension between unitary evolution and collapse in quantum theory, and is sometimes taken as requiring a reassessment of the notion of objective reality. In this note however we present a classical toy model in which (i) The contradicting predictions at the heart of the thought experiment are reproduced (ii) Every system is in a well-defined state at all times. The toy model shows how puzzles such as Wigner's friend's experience of being in a superposition, conflicts between different agents' descriptions of the experiment, the positioning of the Heisenberg's cut and the apparent lack of objectivity of measurement outcomes can be explained within a classical model where there exists an objective state of affairs about every physical system at all times. Within the model, the debate surrounding the original Wigner's friend thought experiment and its resolution have striking similarities with arguments concerning the nature of the second law of thermodynamics. The same conclusion however does not apply to more recent extensions of the gedankenexperiment featuring multiple encapsulated observers, and shows that such extensions are indeed necessary avoid simple classical explanations.


title: Adding eccentricity to quasicircular binary-black-hole waveform models

abstract: The detection of gravitational-wave signals from coalescing eccentric binary black holes would yield unprecedented information about the formation and evolution of compact binaries in specific scenarios, such as dynamical formation in dense stellar clusters and three-body interactions. The gravitational-wave searches by the ground-based interferometers, LIGO and Virgo, rely on analytical waveform models for binaries on quasicircular orbits. Eccentric merger waveform models are less developed, and only few numerical simulations of eccentric mergers are publicly available, but several eccentric inspiral models have been developed from the post-Newtonian expansion. Here we present a novel method to convert the dominant quadrupolar mode of any circular analytical binary-black-hole model into an eccentric model. First, using numerical simulations, we examine the additional amplitude and frequency modulations of eccentric signals that are not present in their circular counterparts. Subsequently, we identify suitable analytical descriptions of those modulations and interpolate key parameters from twelve numerical simulations designated as our training dataset. This allows us to reconstruct the modulated amplitude and phase of any waveform up to mass ratio 3 and eccentricity 0.2. We find that the minimum overlap of the new model with numerical simulations is around 0.98 over all of our test dataset that are scaled to a 50M$_\odot$ black-hole binary starting at 35 Hz with aLIGO A+ design sensitivity. A Python package \pyrex easily carries out the computation of this method.


title: Fiery Cores: Bursty and Smooth Star Formation Distributions across Galaxy Centers in Cosmological Zoom-in Simulations

abstract: We present an analysis of the $R\lesssim 1.5$ kpc core regions of seven simulated Milky Way mass galaxies, from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, for a finely sampled period ($\Delta t = 2.2$ Myr) of 22 Myr at $z \approx 0$, and compare them with star formation rate (SFR) and gas surface density observations of the Milky Way's Central Molecular Zone (CMZ). Despite not being tuned to reproduce the detailed structure of the CMZ, we find that four of these galaxies are consistent with CMZ observations at some point during this 22 Myr period. The galaxies presented here are not homogeneous in their central structures, roughly dividing into two morphological classes; (a) several of the galaxies have very asymmetric gas and SFR distributions, with intense (compact) starbursts occurring over a period of roughly 10 Myr, and structures on highly eccentric orbits through the CMZ, whereas (b) others have smoother gas and SFR distributions, with only slowly varying SFRs over the period analyzed. In class (a) centers, the orbital motion of gas and star-forming complexes across small apertures ($R \lesssim 150$pc, analogously $|l|<1^\circ$ in the CMZ observations) contributes as much to tracers of star formation/dense gas appearing in those apertures, as the internal evolution of those structures does. These asymmetric/bursty galactic centers can simultaneously match CMZ gas and SFR observations, demonstrating that time-varying star formation can explain the CMZ's low star formation efficiency.


title: J\'onsson groups of various cardinalities

abstract: A group $G$ is J\'onsson if $|H| < |G|$ whenever $H$ is a proper subgroup of $G$. Using an embedding theorem of Obraztsov it is shown that there exists a J\'onsson group $G$ of infinite cardinality $\kappa$ if and only if there exists a J\'onsson algebra of cardinality $\kappa$. Thus the question as to which cardinals admit a J\'onsson group is wholly reduced to the well-studied question of which cardinals are not J\'onsson. As a consequence there exist J\'onsson groups of arbitrarily large cardinality. Another consequence is that the infinitary edge-orbit conjecture of Babai is true.


title: Understanding the uneven spread of COVID-19 in the context of the global interconnected economy

abstract: Using network analysis, this paper develops a multidimensional methodological framework for understanding the uneven (cross-country) spread of COVID-19 in the context of the global interconnected economy. The globally interconnected system of tourism mobility is modeled as a complex network, where two main stages in the temporal spread of COVID-19 are revealed and defined by the cutting-point of the 44th day from Wuhan. The first stage describes the outbreak in Asia and North America, the second one in Europe, South America, and Africa, while the outbreak in Oceania is spread along both stages. The analysis shows that highly connected nodes in the global tourism network (GTN) are infected early by the pandemic, while nodes of lower connectivity are late infected. Moreover, countries with the same network centrality as China were early infected on average by COVID-19. The paper also finds that network interconnectedness, economic openness, and transport integration are key determinants in the early global spread of the pandemic, and it reveals that the spatio-temporal patterns of the worldwide spread of COVID-19 are more a matter of network interconnectivity than of spatial proximity.


title: Average Localised Proximity: A new data descriptor with good default one-class classification performance

abstract: One-class classification is a challenging subfield of machine learning in which so-called data descriptors are used to predict membership of a class based solely on positive examples of that class, and no counter-examples. A number of data descriptors that have been shown to perform well in previous studies of one-class classification, like the Support Vector Machine (SVM), require setting one or more hyperparameters. There has been no systematic attempt to date to determine optimal default values for these hyperparameters, which limits their ease of use, especially in comparison with hyperparameter-free proposals like the Isolation Forest (IF). We address this issue by determining optimal default hyperparameter values across a collection of 246 one-class classification problems derived from 50 different real-world datasets. In addition, we propose a new data descriptor, Average Localised Proximity (ALP) to address certain issues with existing approaches based on nearest neighbour distances. Finally, we evaluate classification performance using a leave-one-dataset-out procedure, and find strong evidence that ALP outperforms IF and a number of other data descriptors, as well as weak evidence that it outperforms SVM, making ALP a good default choice.


title: Muppet: Massive Multi-task Representations with Pre-Finetuning

abstract: We propose pre-finetuning, an additional large-scale learning stage between language model pre-training and fine-tuning. Pre-finetuning is massively multi-task learning (around 50 datasets, over 4.8 million total labeled examples), and is designed to encourage learning of representations that generalize better to many different tasks. We show that pre-finetuning consistently improves performance for pretrained discriminators (e.g.~RoBERTa) and generation models (e.g.~BART) on a wide range of tasks (sentence prediction, commonsense reasoning, MRC, etc.), while also significantly improving sample efficiency during fine-tuning. We also show that large-scale multi-tasking is crucial; pre-finetuning can hurt performance when few tasks are used up until a critical point (usually above 15) after which performance improves linearly in the number of tasks.


title: The $(l,r)$-Stirling numbers: a combinatorial approach

abstract: This work deals with a new generalization of $r$-Stirling numbers using $l$-tuple of permutations and partitions called $(l,r)$-Stirling numbers of both kinds. We study various properties of these numbers using combinatorial interpretations and symmetric functions. Also, we give a limit representation of the multiple zeta function using $(l,r)$-Stirling of the first kind.


title: A Comparison of Approaches to Document-level Machine Translation

abstract: Document-level machine translation conditions on surrounding sentences to produce coherent translations. There has been much recent work in this area with the introduction of custom model architectures and decoding algorithms. This paper presents a systematic comparison of selected approaches from the literature on two benchmarks for which document-level phenomena evaluation suites exist. We find that a simple method based purely on back-translating monolingual document-level data performs as well as much more elaborate alternatives, both in terms of document-level metrics as well as human evaluation.


title: Complementary Composite Minimization, Small Gradients in General Norms, and Applications

abstract: Composite minimization is a powerful framework in large-scale convex optimization, based on decoupling of the objective function into terms with structurally different properties and allowing for more flexible algorithmic design. We introduce a new algorithmic framework for complementary composite minimization, where the objective function decouples into a (weakly) smooth and a uniformly convex term. This particular form of decoupling is pervasive in statistics and machine learning, due to its link to regularization. The main contributions of our work are summarized as follows. First, we introduce the problem of complementary composite minimization in general normed spaces; second, we provide a unified accelerated algorithmic framework to address broad classes of complementary composite minimization problems; and third, we prove that the algorithms resulting from our framework are near-optimal in most of the standard optimization settings. Additionally, we show that our algorithmic framework can be used to address the problem of making the gradients small in general normed spaces. As a concrete example, we obtain a nearly-optimal method for the standard $\ell_1$ setup (small gradients in the $\ell_{\infty}$ norm), essentially matching the bound of Nesterov (2012) that was previously known only for the Euclidean setup. Finally, we show that our composite methods are broadly applicable to a number of regression and other classes of optimization problems, where regularization plays a key role. Our methods lead to complexity bounds that are either new or match the best existing ones.


title: Translational and reorientational dynamics in deep eutectic solvents

abstract: We performed rheological measurements of the typical deep eutectic solvents (DESs) glyceline, ethaline, and reline in a very broad temperature and dynamic range, extending from the low-viscosity to the high-viscosity supercooled-liquid regime. We find that the mechanical compliance spectra can be well described by the random free-energy barrier hopping model, while the dielectric spectra on the same materials involve significant contributions arising from reorientational dynamics. The temperature-dependent viscosity and structural relaxation time, revealing non-Arrhenius behavior typical for glassy freezing, are compared to the ionic dc conductivity and relaxation times determined by broadband dielectric spectroscopy. For glyceline and ethaline we find essentially identical temperature dependences for all dynamic quantities. These findings point to a close coupling of the ionic and molecular translational and reorientational motions in these systems. However, for reline the ionic charge transport appears decoupled from the structural and reorientational dynamics, following a fractional Walden rule. Especially, at low temperatures the ionic conductivity in this DES is enhanced by about one decade compared to expectations based on the temperature dependence of the viscosity. The results for all three DESs can be understood without invoking a revolving-door mechanism previously considered as a possible charge-transport mechanism in DESs.


title: Deep Subjecthood: Higher-Order Grammatical Features in Multilingual BERT

abstract: We investigate how Multilingual BERT (mBERT) encodes grammar by examining how the high-order grammatical feature of morphosyntactic alignment (how different languages define what counts as a "subject") is manifested across the embedding spaces of different languages. To understand if and how morphosyntactic alignment affects contextual embedding spaces, we train classifiers to recover the subjecthood of mBERT embeddings in transitive sentences (which do not contain overt information about morphosyntactic alignment) and then evaluate them zero-shot on intransitive sentences (where subjecthood classification depends on alignment), within and across languages. We find that the resulting classifier distributions reflect the morphosyntactic alignment of their training languages. Our results demonstrate that mBERT representations are influenced by high-level grammatical features that are not manifested in any one input sentence, and that this is robust across languages. Further examining the characteristics that our classifiers rely on, we find that features such as passive voice, animacy and case strongly correlate with classification decisions, suggesting that mBERT does not encode subjecthood purely syntactically, but that subjecthood embedding is continuous and dependent on semantic and discourse factors, as is proposed in much of the functional linguistics literature. Together, these results provide insight into how grammatical features manifest in contextual embedding spaces, at a level of abstraction not covered by previous work.


title: Voltage-controlled antiferromagnetism in magnetic tunnel junctions

abstract: We demonstrate a voltage-controlled exchange bias effect in CoFeB/MgO/CoFeB magnetic tunnel junctions that is related to the interfacial Fe(Co)Ox formed between the CoFeB electrodes and the MgO barrier. The unique combination of interfacial antiferromagnetism, giant tunneling magnetoresistance, and sharp switching of the perpendicularly-magnetized CoFeB allows sensitive detection of the exchange bias. It is found that the exchange bias field can be isothermally controlled by magnetic fields at low temperatures. More importantly, the exchange bias can also be effectively manipulated by the electric field applied to the MgO barrier due to the voltage-controlled antiferromagnetic anisotropy in this system.


title: On the Support of a hypoelliptic diffusion on the Heisenberg group

abstract: We provide an elementary proof of the support of the law of a hypoelliptic Brownian motion on the Heisenberg group $\mathbb{H}$. We consider a control norm associated to left-invariant vector fields on $\mathbb{H}$, and describe the support in terms of the space of finite energy horizontal curves.


title: Generalized Doubly Reparameterized Gradient Estimators

abstract: Efficient low-variance gradient estimation enabled by the reparameterization trick (RT) has been essential to the success of variational autoencoders. Doubly-reparameterized gradients (DReGs) improve on the RT for multi-sample variational bounds by applying reparameterization a second time for an additional reduction in variance. Here, we develop two generalizations of the DReGs estimator and show that they can be used to train conditional and hierarchical VAEs on image modelling tasks more effectively. First, we extend the estimator to hierarchical models with several stochastic layers by showing how to treat additional score function terms due to the hierarchical variational posterior. We then generalize DReGs to score functions of arbitrary distributions instead of just those of the sampling distribution, which makes the estimator applicable to the parameters of the prior in addition to those of the posterior.


title: A Census of Blue Stragglers in Gaia DR2 Open Clusters as a Test of Population Synthesis and Mass Transfer Physics

abstract: We use photometry and proper motions from Gaia DR2 to determine the blue straggler star (BSS) populations of 16 old (1-10 Gyr), nearby ($d< 3500$ pc) open clusters. We find that the fractional number of BSS compared to RGB stars increases with age, starting near zero at 1 Gyr and flattening to $\sim 0.35$ by 4 Gyr. Fitting stellar evolutionary tracks to these BSS, we find that their mass distribution peaks at a few tenths of a solar mass above the main-sequence turnoff. BSS more than 0.5 $M_\odot$ above the turnoff make up only $\sim25$\% of the sample, and BSS more than 1.0 $M_\odot$ above the turnoff are rare. We compare this to population synthesis models of BSS formed via mass transfer using the Compact Object Synthesis and Monte Carlo Investigation Code (COSMIC). We find that standard population synthesis assumptions dramatically under-produce the number of BSS in old open clusters. We also find that these models overproduce high mass BSS relative to lower mass BSS. Expected numbers of BSS formed through dynamics do not fully account for this discrepancy. We conclude that in order to explain the observed BSS populations from Roche lobe overflow, mass-transfer from giant donors must be more stable than assumed in canonical mass-transfer prescriptions, and including non-conservative mass transfer is important in producing realistic BSS masses. Even with these modifications, it is difficult to achieve the large number of BSS observed in the oldest open clusters. We discuss some additional physics that may explain the large number of observed blue stragglers among old stellar populations.


title: Large-scale Optimization with Linear Equality Constraints using Reduced Compact Representation

abstract: For optimization problems with linear equality constraints, we prove that the (1,1) block of the inverse KKT matrix remains unchanged when projected onto the nullspace of the constraint matrix. We develop reduced compact representations of the limited-memory inverse BFGS Hessian to compute search directions efficiently when the constraint Jacobian is sparse. Orthogonal projections are implemented by a sparse QR factorization or a preconditioned LSQR iteration. In numerical experiments two proposed trust-region algorithms improve in computation times, often significantly, compared to previous implementations of related algorithms and compared to IPOPT.


title: C-for-Metal: High Performance SIMD Programming on Intel GPUs

abstract: The SIMT execution model is commonly used for general GPU development. CUDA and OpenCL developers write scalar code that is implicitly parallelized by compiler and hardware. On Intel GPUs, however, this abstraction has profound performance implications as the underlying ISA is SIMD and important hardware capabilities cannot be fully utilized. To close this performance gap we introduce C-For-Metal (CM), an explicit SIMD programming framework designed to deliver close-to-the-metal performance on Intel GPUs. The CM programming language and its vector/matrix types provide an intuitive interface to exploit the underlying hardware features, allowing fine-grained register management, SIMD size control and cross-lane data sharing. Experimental results show that CM applications from different domains outperform the best-known SIMT-based OpenCL implementations, achieving up to 2.7x speedup on the latest Intel GPU.


title: Non-tautological Hurwitz cycles

abstract: We show that various loci of stable curves of sufficiently large genus admitting degree $d$ covers of positive genus curves define non-tautological algebraic cycles on $\overline{\mathcal{M}}_{g,N}$, assuming the non-vanishing of the $d$-th Fourier coefficient of a certain modular form. Our results build on those of Graber-Pandharipande and van Zelm for degree 2 covers of elliptic curves; the main new ingredient is a method to intersect the cycles in question with boundary strata, as developed recently by Schmitt-van Zelm and the author.


title: A Paradox of Whistling Atmospherics

abstract: Analyzing paradoxes is interesting and instructive. Sometimes the analysis leads to non-trivial results. This methodological note sets out a paradox arising in the theory of propagation of electromagnetic waves in moving plasmas. The paradox is interesting in itself, and, generally speaking, it should be taken into account when analyzing geoelectromagnetic waves. The paradox is as follows: contrary to expectations, the group velocity of the waves is the same in the comoving and laboratory frames of reference. The condition for the appearance of the paradox is the quadratic dependence of the frequency on the wave number. A paradoxical property manifests itself in the theory of the propagation of radio waves (in particular, whistling atmospherics), Langmuir waves and Alfv\'en waves. From a cognitive point of view, it is interesting that the paradox can be traced in relation to de Broglie waves. An explanation of the paradox is proposed. Keywords: group velocity, moving plasma, Doppler Effect, dispersion, longitudinal waves, transverse waves.


title: Energy Non-Conservation in Quantum Mechanics

abstract: We study the conservation of energy, or lack thereof, when measurements are performed in quantum mechanics. The expectation value of the Hamiltonian of a system can clearly change when wave functions collapse in accordance with the standard textbook (Copenhagen) treatment of quantum measurement, but one might imagine that the change in energy is compensated by the measuring apparatus or environment. We show that this is not true; the change in the energy of a state after measurement can be arbitrarily large, independent of the physical measurement process. In Everettian quantum theory, while the expectation value of the Hamiltonian is conserved for the wave function of the universe (including all the branches), it is not constant within individual worlds. It should therefore be possible to experimentally measure violations of conservation of energy, and we suggest an experimental protocol for doing so.


title: Response-Time Analysis and Optimization for Probabilistic Conditional Parallel DAG Tasks

abstract: Real-time systems increasingly use multicore processors in order to satisfy thermal, power, and computational requirements. To exploit the architectural parallelism offered by the multicore processors, parallel task models, scheduling algorithms and response-time analyses with respect to real-time constraints have to be provided. In this paper, we propose a reservation-based scheduling algorithm for sporadic constrained-deadline parallel conditional DAG tasks with probabilistic execution behaviour for applications that can tolerate bounded number of deadline misses and bounded tardiness. We devise design rules and analyses to guarantee bounded tardiness for a specified bounded probability for $k$-consecutive deadline misses without enforcing late jobs to be immediately aborted.


title: Remote Learners, Home Makers: How Digital Fabrication Was Taught Online During a Pandemic

abstract: Digital fabrication courses that relied on physical makerspaces were severely disrupted by COVID-19. As universities shut down in Spring 2020, instructors developed new models for digital fabrication at a distance. Through interviews with faculty and students and examination of course materials, we recount the experiences of eight remote digital fabrication courses. We found that learning with hobbyist equipment and online social networks could emulate using industrial equipment in shared workshops. Furthermore, at-home digital fabrication offered unique learning opportunities including more iteration, machine tuning, and maintenance. These opportunities depended on new forms of labor and varied based on student living situations. Our findings have implications for remote and in-person digital fabrication instruction. They indicate how access to tools was important, but not as critical as providing opportunities for iteration; they show how remote fabrication exacerbated student inequities; and they suggest strategies for evaluating trade-offs in remote fabrication models with respect to learning objectives.


title: LDLE: Low Distortion Local Eigenmaps

abstract: We present Low Distortion Local Eigenmaps (LDLE), a manifold learning technique which constructs a set of low distortion local views of a dataset in lower dimension and registers them to obtain a global embedding. The local views are constructed using the global eigenvectors of the graph Laplacian and are registered using Procrustes analysis. The choice of these eigenvectors may vary across the regions. In contrast to existing techniques, LDLE can embed closed and non-orientable manifolds into their intrinsic dimension by tearing them apart. It also provides gluing instruction on the boundary of the torn embedding to help identify the topology of the original manifold. Our experimental results will show that LDLE largely preserved distances up to a constant scale while other techniques produced higher distortion. We also demonstrate that LDLE produces high quality embeddings even when the data is noisy or sparse.


title: A Model of Densifying Collaboration Networks

abstract: Research collaborations provide the foundation for scientific advances, but we have only recently begun to understand how they form and grow on a global scale. Here we analyze a model of the growth of research collaboration networks to explain the empirical observations that the number of collaborations scales superlinearly with institution size, though at different rates (heterogeneous densification), the number of institutions grows as a power of the number of researchers (Heaps' law) and institution sizes approximate Zipf's law. This model has three mechanisms: (i) researchers are preferentially hired by large institutions, (ii) new institutions trigger more potential institutions, and (iii) researchers collaborate with friends-of-friends. We show agreement between these assumptions and empirical data, through analysis of co-authorship networks spanning two centuries. We then develop a theoretical understanding of this model, which reveals emergent heterogeneous scaling such that the number of collaborations between institutions scale with an institution's size.


title: Singular integrals with variable kernels in dyadic settings

abstract: In this paper we explore conditions on variable symbols with respect to Haar systems, defining Calder\'on-Zygmund type operators with respect to the dyadic metrics associated to the Haar bases.We show that Petermichl's dyadic kernel can be seen as a variable kernel singular integral and we extend it to dyadic systems built on spaces of homogeneous type.


title: Supervised Momentum Contrastive Learning for Few-Shot Classification

abstract: Few-shot learning aims to transfer information from one task to enable generalization on novel tasks given a few examples. This information is present both in the domain and the class labels. In this work we investigate the complementary roles of these two sources of information by combining instance-discriminative contrastive learning and supervised learning in a single framework called Supervised Momentum Contrastive learning (SUPMOCO). Our approach avoids a problem observed in supervised learning where information in images not relevant to the task is discarded, which hampers their generalization to novel tasks. We show that (self-supervised) contrastive learning and supervised learning are mutually beneficial, leading to a new state-of-the-art on the META-DATASET - a recently introduced benchmark for few-shot learning. Our method is based on a simple modification of MOCO and scales better than prior work on combining supervised and self-supervised learning. This allows us to easily combine data from multiple domains leading to further improvements.


title: Event-Driven News Stream Clustering using Entity-Aware Contextual Embeddings

abstract: We propose a method for online news stream clustering that is a variant of the non-parametric streaming K-means algorithm. Our model uses a combination of sparse and dense document representations, aggregates document-cluster similarity along these multiple representations and makes the clustering decision using a neural classifier. The weighted document-cluster similarity model is learned using a novel adaptation of the triplet loss into a linear classification objective. We show that the use of a suitable fine-tuning objective and external knowledge in pre-trained transformer models yields significant improvements in the effectiveness of contextual embeddings for clustering. Our model achieves a new state-of-the-art on a standard stream clustering dataset of English documents.


title: Defenses Against Multi-Sticker Physical Domain Attacks on Classifiers

abstract: Recently, physical domain adversarial attacks have drawn significant attention from the machine learning community. One important attack proposed by Eykholt et al. can fool a classifier by placing black and white stickers on an object such as a road sign. While this attack may pose a significant threat to visual classifiers, there are currently no defenses designed to protect against this attack. In this paper, we propose new defenses that can protect against multi-sticker attacks. We present defensive strategies capable of operating when the defender has full, partial, and no prior information about the attack. By conducting extensive experiments, we show that our proposed defenses can outperform existing defenses against physical attacks when presented with a multi-sticker attack.


title: Nonergodic dynamics in disorder-free potentials

abstract: We review the dynamics of interacting particles in disorder-free potentials concentrating on a combination of a harmonic binding with a constant tilt. We show that a simple picture of an effective local tilt describes a variety of cases. Our examples include spinless fermions (as modeled by Heisenberg spin chain in a magnetic field), spinful fermions as well as bosons that enjoy a larger local on-site Hilbert space. We also discuss the domain-wall dynamics that reveals nonergodic features even for a relatively weak tilt as suggested by Doggen et. al. [arXiv:2012.13722]. By adding a harmonic potential on top of the static field we confirm that the surprizing regular dynamics is not entirely due to Hilbert space shuttering. It seems better explained by the inhibited transport within the domains of identically oriented spins. Once the spin-1/2 restrictions are lifted as, e.g., for bosons, the dynamics involve stronger entanglement generation. Again for domain wall melting, the effect of the harmonic potential is shown to lead mainly to an effective local tilt.


title: The Competition Between Deformation Twinning and Dislocation Slip in Deformed Face-Centered Cubic Metals

abstract: The competition between deformation twinning and dislocation slip underpins the evolution of mesoscale plasticity in face-centered cubic materials. While competition between these mechanisms is known to be related to the critical features of the generalized planar fault energy landscape, a physical theory that tracks competition over extended plasticity has yet to emerge. Here, we report a methodology to predict the mesoscale evolution of this competition in deformed crystals. Our approach implements kinetic Monte Carlo simulations to examine fault structure evolution in face-centered cubic metals using intrinsic material parameters as inputs. These results are leveraged to derive an analytical model for the evolution of the fault fraction, fault densities, and partitioning of plastic strains among deformation mechanisms. In addition, we define a competition parameter that measures the tendencies for deformation twinning and dislocation slip. In contrast to previous twinnability parameters, our derivation considers deformation history when examining mechanism competition. This contribution therefore extends the reach of deformation twinning theory beyond incipient nucleation events. These products find direct applications in work hardening and crystal plasticity models, which have previously relied on phenomenological relations to predict the mesoscale evolution of deformation twin microstructures.


title: Strong Shock in the Uniformly Expanding Universe with a Spherical Void

abstract: Propagation of strong shock wave in the expanding universe is studied using approximate analytic, and exact numerical solution of self-similar equations. Both solutions have similar properties, which change qualitatively, depending on the adiabatic powers $\gamma$. In the interval $1<\gamma<\gamma_{cr} \sim 1.16$ analytic and numeric solutions fill all the space without any voids and they are rather close to each other. At larger $\gamma>\gamma_{cr}$ a pressure becomes zero at finite radius, and a spherical void appears around the origin in both solutions. All matter is collected in thin layer behind the shock wave front. The structure of this layer qualitatively depends on $\gamma$. At the inner edge of the layer the pressure is always zero, but the density on this edge is jumping from zero to infinity at $\gamma \approx 1.4$ in both solutions.


title: A geometric approach to Lie systems: formalism of Poisson-Hopf algebra deformations

abstract: The notion of quantum algebras is merged with that of Lie systems in order to establish a new formalism called Poisson-Hopf algebra deformations of Lie systems. The procedure can be naturally applied to Lie systems endowed with a symplectic structure, the so-called Lie-Hamilton systems. This is quite a general approach, as it can be applied to any quantum deformation and any underlying manifold. One of its main features is that, Lie systems are extended to generalized systems described by involutive distributions. In this way, we obtain their new generalized (deformed) counterparts that cover, in particular, a new oscillator system with a time-dependent frequency and a position-dependent mass. Based on a recently developed procedure to construct Poisson-Hopf deformations of Lie-Hamilton systems \cite{Ballesteros5}, a novel unified approach to deformations of Lie-Hamilton systems on the real plane with a Vessiot-Guldberg Lie algebra isomorphic to $\mathfrak{sl}(2)$ is proposed. In addition, we study the deformed systems obtained from Lie-Hamilton systems associated to the oscillator algebra $\mathfrak{h}_4$, seen as a subalgebra of the 2-photon algebra $\mathfrak{h}_6$. As a particular application, we propose an epidemiological model of SISf type that uses the solvable Lie algebra $\mathfrak{b}_2$ as subalgebra of $\mathfrak{sl}(2)$, by restriction of the corresponding quantum deformed systems.


title: Inhomogeneous mean-field approach to collective excitations in disordered interacting bosons

abstract: We develop an inhomogeneous quantum mean-field approach to the behavior of collective excitations across the superfluid-Mott glass quantum phase transition in two dimensions, complementing recent quantum Monte Carlo simulations [Phys. Rev. Lett. {\bf 125}, 027002 (2020)]. In quadratic (Gaussian) approximation, the Goldstone (phase) and Higgs (amplitude) modes completely decouple. Each is described by a disordered Bogoliubov Hamiltonian which can be solved by an inhomogeneous multi-mode Bogoliubov transformation. We find that the Higgs mode is spatially localized in both phases. The corresponding scalar spectral function shows a broad peak that is noncritical in the sense that its peak frequency does not soften but remains nonzero across the quantum phase transition. In contrast, the lowest-energy Goldstone mode delocalizes in the superfluid phase, leading to a zero-frequency spectral peak. We compare these findings to the results of the quantum Monte Carlo simulations. We also relate them to general results on the localization of bosonic excitations, and we discuss the limits and generality of our approach.


title: Logical-Combinatorial Approaches in Dynamic Recognition Problems

abstract: A pattern recognition scenario, where instead of object classification into the classes by the learning set, the algorithm aims to allocate all objects to the same, the so-called normal class, is the research objective.


title: Optical polarizability of erbium-oxygen complexes in sol-gel-based silica films

abstract: For erbium-doped amorphous oxides, such as those that are used in compact lightwave devices interfaced with silicon, values of the refractive indices are commonly obtained empirically. This work, combining experimental and theoretical studies, examines silica as the exemplary host and the influence of erbium doping on the refractive index. Analysis of data is presented for the spectral refractive index in the ultraviolet to near infrared wavelength range of heavily erbium-doped silica thin films prepared by spin coating a sol-gel precursor on silicon and subsequent vacuum annealing. Effective medium Lorentz-Lorenz data analysis determines that the dopant component has a refractive index of 1.76(0.24) with wavelength dispersion constrained to within 2 percent. Considering the dopant as a localized ErO6 impurity complex, a corresponding theoretical refractive index of 1.662 is derived by calculating the optical polarizability and volume of the impurity. Data presented for room-temperature (293 K) photoluminescence in the vicinity of 1.54 micron are shown to be consistent with random variability in impurity sites. Inherent advantages of studying colloid-based materials are discussed. To the best of the authors knowledge, such a detailed study of the refractive index associated with erbium impurities in silica is being reported for the first time in the literature.


title: The Galactic Chemical Evolution of Chlorine

abstract: We measured $^{35}$Cl abundances in 52 M giants with metallicities between -0.5 $<$ [Fe/H] $<$ 0.12. Abundances and atmospheric parameters were derived using infrared spectra from CSHELL on the IRTF and from optical echelle spectra. We measured Cl abundances by fitting a H$^{35}$Cl molecular feature at 3.6985 $\mu$m with synthetic spectra. We also measured the abundances of O, Ca, Ti, and Fe using atomic absorption lines. We find that the [Cl/Fe] ratio for our stars agrees with chemical evolution models of Cl and the [Cl/Ca] ratio is broadly consistent with the solar ratio over our metallicity range. Both indicate that Cl is primarily made in core-collapse supernovae with some contributions from Type Ia SN. We suggest other potential nucleosynthesis processes, such as the $\nu$-process, are not significant producers of Cl. Finally, we also find our Cl abundances are consistent with H II and planetary nebular abundances at a given oxygen abundance, although there is scatter in the data.


title: The Neupert Effect of Flare UltraViolet and Soft X-ray Emissions

abstract: We model the Neupert effect that relates flare heating energy with the observed SXR emission. The traditional form of the Neupert effect refers to the correlation between the time-integrated HXR or microwave light curve and the SXR light curve. In this paper, instead, we use as the proxy for heating energy the ultraviolet (UV) emission at the foot-points of flare loops, and modify the model of the Neupert effect by taking into account the discrete nature of flare heating as well as cooling. In the modified empirical model, spatially resolved UV lightcurves from the transition region or upper chromosphere are each convolved with a kernel function characterizing the decay of the flare loop emission. Contributions by all loops are summed to compare with the observed total SXR emission. The model has successfully reproduced the observed SXR emission from its rise to decay. To estimate heating energies in flare loops, we also employ the UV Foot-point Calorimeter (UFC) method that infers heating rates in flare loops from these UV light curves and models evolution of flare loops with a zero-dimensional hydrodynamic code. The experiments show that a multitude of impulsive heating events do not well reproduce the observed flare SXR light curve, but a two-phase heating model leads to better agreement with observations. Comparison of the two models of the Neupert effect further allows us to calibrate the UFC method, and improve the estimate of heating rates in flare loops continuously formed by magnetic reconnection throughout the flare evolution.


title: REFORM: Fast and Adaptive Solution for Subteam Replacement

abstract: In this paper, we propose the novel problem of Subteam Replacement: given a team of people embedded in a social network to complete a certain task, and a subset of members - subteam - in this team which have become unavailable, find another set of people who can perform the subteam's role in the larger team. The ability to simultaneously replace multiple team members is highly appreciated in settings such as corporate management where team structure is highly volatile and large-scale changes are commonplace. We conjecture that a good candidate subteam should have high skill and structural similarity with the replaced subteam while sharing a similar connection with the larger team as a whole. Based on this conjecture, we propose a novel graph kernel which evaluates the goodness of candidate subteams in this holistic way freely adjustable to the need of the situation. To tackle the significant computational difficulties, we combine our kernel with a fast approximate algorithm which (a) employs effective pruning strategies, (b) exploits the similarity between candidate team structures to reduce kernel computations, and (c) features a solid theoretical bound obtained from mathematical properties of the problem. We extensively test our solution on both synthetic and real datasets to demonstrate its consistency and efficiency. Our proposed graph kernel results in more suitable replacements being proposed compared to graph kernels used in previous work, and our algorithm consistently outperforms alternative choices by finding near-optimal solutions while scaling linearly with the size of the replaced subteam.


title: The MineRL 2020 Competition on Sample Efficient Reinforcement Learning using Human Priors

abstract: Although deep reinforcement learning has led to breakthroughs in many difficult domains, these successes have required an ever-increasing number of samples, affording only a shrinking segment of the AI community access to their development. Resolution of these limitations requires new, sample-efficient methods. To facilitate research in this direction, we propose this second iteration of the MineRL Competition. The primary goal of the competition is to foster the development of algorithms which can efficiently leverage human demonstrations to drastically reduce the number of samples needed to solve complex, hierarchical, and sparse environments. To that end, participants compete under a limited environment sample-complexity budget to develop systems which solve the MineRL ObtainDiamond task in Minecraft, a sequential decision making environment requiring long-term planning, hierarchical control, and efficient exploration methods. The competition is structured into two rounds in which competitors are provided several paired versions of the dataset and environment with different game textures and shaders. At the end of each round, competitors submit containerized versions of their learning algorithms to the AIcrowd platform where they are trained from scratch on a hold-out dataset-environment pair for a total of 4-days on a pre-specified hardware platform. In this follow-up iteration to the NeurIPS 2019 MineRL Competition, we implement new features to expand the scale and reach of the competition. In response to the feedback of the previous participants, we introduce a second minor track focusing on solutions without access to environment interactions of any kind except during test-time. Further we aim to prompt domain agnostic submissions by implementing several novel competition mechanics including action-space randomization and desemantization of observations and actions.


title: Index matching computerized tomography

abstract: Computerized tomography (CT) has been used for decades by medical professionals to detect and diagnose injuries and ailments. CT scanners are based on interesting physics, but due to their bulk, cost, and safety, hands on experience with a medical CT scanner is unrealistic for undergraduate students. Therefore, operationally similar, yet small, safe, and inexpensive CT scanners are desirable teaching tools. This project details the development of a novel model CT scanning apparatus. The experimental setup presented utilizes visible light, has short data acquisition time, and operates on the same physics as its X-ray counterpart. The apparatus employs a laser and a photodiode to image a transparent material, while avoiding loss of transmitted intensity through index of refraction matching. A simple back-projection algorithm results in a 2D cross section of the scan object. We found we could collect data and reliably image samples in 15 min.


title: Property Inference From Poisoning

abstract: Property inference attacks consider an adversary who has access to the trained model and tries to extract some global statistics of the training data. In this work, we study property inference in scenarios where the adversary can maliciously control part of the training data (poisoning data) with the goal of increasing the leakage. Previous work on poisoning attacks focused on trying to decrease the accuracy of models either on the whole population or on specific sub-populations or instances. Here, for the first time, we study poisoning attacks where the goal of the adversary is to increase the information leakage of the model. Our findings suggest that poisoning attacks can boost the information leakage significantly and should be considered as a stronger threat model in sensitive applications where some of the data sources may be malicious. We describe our \emph{property inference poisoning attack} that allows the adversary to learn the prevalence in the training data of any property it chooses. We theoretically prove that our attack can always succeed as long as the learning algorithm used has good generalization properties. We then verify the effectiveness of our attack by experimentally evaluating it on two datasets: a Census dataset and the Enron email dataset. We were able to achieve above $90\%$ attack accuracy with $9-10\%$ poisoning in all of our experiments.


title: Patient Diversion Across Primary Health Centers Using Real Time Delay Predictors

abstract: In the current work, we consider diversion of childbirth patients who arrive seeking emergency admission to public primary health centers (PHCs). PHCs are the first point of contact for an Indian patient with formal medical care, and offer medical care on an outpatient basis, and limited inpatient and childbirth care. In this context, real-time prediction of the wait time of the arriving patient becomes important in order to determine whether the patient must be diverted to another PHC or not. We study this problem using a discrete event simulation that we develop of medical care operations in two PHCs in India. We approximate the labour room service at each PHC as an M/G/1 queueing system and show how the accuracy of real-time delay predictors impacts the extent of the change in operational outcomes at each PHC. We simulate patient diversion using actual delays as well as the delay estimates generated by various delay predictors based on the state of the system such as queue-length, elapsed service time, and observed delay histories. The simulation of the diversion process also incorporates travel time between the PHCs. We also propose a new delay predictor that incorporates information regarding the system state as well as the service time distribution. We compare the operational outcomes at both PHCs without diversion and with diversion using the above delay predictors. We show numerically that more accurate delay predictors lead to more equitable distribution of resources involved in provision of childbirth care across both PHCs.


title: Adaptivity without Compromise: A Momentumized, Adaptive, Dual Averaged Gradient Method for Stochastic Optimization

abstract: We introduce MADGRAD, a novel optimization method in the family of AdaGrad adaptive gradient methods. MADGRAD shows excellent performance on deep learning optimization problems from multiple fields, including classification and image-to-image tasks in vision, and recurrent and bidirectionally-masked models in natural language processing. For each of these tasks, MADGRAD matches or outperforms both SGD and ADAM in test set performance, even on problems for which adaptive methods normally perform poorly.


title: Broadband terahertz half-wave plate with multi-layers metamaterial via quantum engineering

abstract: In this paper, we employ the novel design of the metamaterial half-wave plate by using the multiple layers of the metamaterials with some specific rotation angles. The rotation angles are given by composite pulse control which is the well-known quantum control technique. A big advantage of this method can analytically calculate the rotation angles and can be easily extended to the multiple layers. The more layers of the metamaterials can present more bandwidth of half-wave plate. In this paper, we present 1, 3, 5, 7 configurations of layers of a metamaterial half-wave plate and we demonstrate that our device increasingly enhance the bandwidth of performance.


title: New Findings on GLRT Radar Detection of Nonfluctuating Targets via Phased Arrays

abstract: This paper addresses the standard generalized likelihood ratio test (GLRT) detection problem of weak signals in background noise. In so doing, we consider a nonfluctuating target embedded in complex white Gaussian noise (CWGN), in which the amplitude of the target echo and the noise power are assumed to be unknown. Important works have analyzed the performance for the referred scenario and proposed GLRT-based detectors. Such detectors are projected at an early stage (i.e., prior to the formation of a post-beamforming scalar waveform), thereby imposing high demands on hardware, processing, and data storage. From a hardware perspective, most radar systems fail to meet these strong requirements. In fact, due to hardware and computational constraints, most radars use a combination of analog and digital beamformers (sums) before any estimation or further pre-processing. The rationale behind this study is to derive a GLRT detector that meets the hardware and system requirements. In this work, we design and analyze a more practical and easy-to-implement GLRT detector, which is projected after the analog beamforming. The performance of the proposed detector is analyzed and the probabilities of detection (PD) and false alarm (PFA) are derived in closed form. Moreover, we show that in the low signal-to-noise ratio (SNR) regime, the post-beamforming GLRT detector performs better than both the classic pre-beamforming GLRT detector and the square-law detector. This finding suggests that if the signals are weak, instead of processing the signals separately, we first must to reinforce the overall signal and then assembling the system's detection statistic. At last, the SNR losses are quantified, in which the superiority of the post-beamforming GLRT detector was evidenced as the number of antennas and samples increase.


title: A Generalization of the Greene-Kleitman Duality Theorem

abstract: In this paper, we describe and prove a generalization of both the classical Greene-Kleitman duality theorem for posets and the local version proved recently by Lewis-Lyu-Pylyavskyy-Sen in studying discrete solitons, using an approach more closely linked to the approach of the classical case.


title: Blind Reconstruction of Multilayered Tissue Profiles with UWB Radar Under Bayesian Setting

abstract: In this paper, we investigate the problem of inverse electromagnetic scattering to recover multilayer human tissue profiles using ultrawideband radar systems in a Bayesian setting. We study the recovery problem in a blind setting, in which we simultaneously estimate both the dielectric/geometric properties of the one-dimensional target tissue profile and the transmitted radar waveform. To perform Bayesian parameter estimation, we propose a hybrid and adaptive Markov Chain Monte Carlo method, which combines the Slice sampling and Hamiltonian Monte Carlo approaches. The introduced sampling mechanism also incorporates the Parallel Tempering approach to escape from the local optimal regions of the complex posterior distribution. We provide empirical support through various numerical simulations for the achieved enhanced sampling efficiency compared to conventional sampling schemes. To investigate the recovery performance, we work on synthetic measurements simulating actual radar returns from multilayer tissue profiles. We derive theoretical bounds for the best achievable estimation performance in terms of normalized root mean square error and provide a comparison with the performance of our estimator.