title: Rapid and Robust Parameter Inference for Binary Mergers

abstract: The detection rate for compact binary mergers has grown as the sensitivity of the global network of ground based gravitational wave detectors has improved, now reaching the stage where robust automation of the analyses is essential. Automated low-latency algorithms have been developed that send out alerts when candidate signals are detected. The alerts include sky maps to facilitate electromagnetic follow up observations, along with probabilities that the system might contain a neutron star, and hence be more likely to generate an electromagnetic counterpart. Data quality issues, such as loud noise transients (glitches), can adversely affect the low-latency algorithms, causing false alarms and throwing off parameter estimation. Here a new analysis method is presented that is robust against glitches, and capable of producing fully Bayesian parameter inference, including sky maps and mass estimates, in a matter of minutes. Key elements of the method are wavelet-based de-noising, penalized maximization of the likelihood during the initial search, rapid sky localization using pre-computed inner products, and heterodyned likelihoods for full Bayesian inference.


title: Ultrafast Frustration-Breaking and Magnetophononic Driving of Singlet Excitations in a Quantum Magnet

abstract: Ideal magnetic frustration forms the basis for the emergence of exotic quantum spin states that are entirely nonmagnetic. Such singlet spin states are the defining feature of the Shastry-Sutherland model, and of its faithful materials realization in the quantum antiferromagnet SrCu$_2$(BO$_3$)$_2$. To address these states on ultrafast timescales, despite their lack of any microscopic order parameter, we introduce a nonlinear magnetophononic mechanism to alter the quantum spin dynamics by driving multiple optical phonon modes coherently and simultaneously. We apply intense terahertz pulses to create a nonequilibrium modulation of the magnetic interactions that breaks the ideal frustration of SrCu$_2$(BO$_3$)$_2$, such that previously forbidden physics can be driven in a coherent manner. Specifically, this driving populates a purely magnetic excitation, the singlet branch of the two-triplon bound state, by resonance with the difference frequency of two pumped phonons. Our results demonstrate how light-driven phonons can be used for the ultrafast and selective manipulation of interactions in condensed matter, even at frequencies far from those of the pump spectrum, offering valuable additional capabilities for the dynamical control of quantum many-body phenomena.


title: Control of Stochastic Quantum Dynamics by Differentiable Programming

abstract: Control of the stochastic dynamics of a quantum system is indispensable in fields such as quantum information processing and metrology. However, there is no general ready-made approach to the design of efficient control strategies. Here, we propose a framework for the automated design of control schemes based on differentiable programming ($\partial \mathrm{P}$). We apply this approach to the state preparation and stabilization of a qubit subjected to homodyne detection. To this end, we formulate the control task as an optimization problem where the loss function quantifies the distance from the target state, and we employ neural networks (NNs) as controllers. The system's time evolution is governed by a stochastic differential equation (SDE). To implement efficient training, we backpropagate the gradient information from the loss function through the SDE solver using adjoint sensitivity methods. As a first example, we feed the quantum state to the controller and focus on different methods of obtaining gradients. As a second example, we directly feed the homodyne detection signal to the controller. The instantaneous value of the homodyne current contains only very limited information on the actual state of the system, masked by unavoidable photon-number fluctuations. Despite the resulting poor signal-to-noise ratio, we can train our controller to prepare and stabilize the qubit to a target state with a mean fidelity of around 85%. We also compare the solutions found by the NN to a hand-crafted control strategy.


title: Positivity in Multi-Field EFTs

abstract: We discuss the general method for obtaining full positivity bounds on multi-field effective field theories (EFTs). While the leading order forward positivity bounds are commonly derived from the elastic scattering of two (superposed) external states, we show that for a generic EFT containing 3 or more low-energy modes, this approach only gives incomplete bounds. We then identify the allowed parameter space as the dual to a spectrahedron, constructed from crossing symmetries of the amplitude, and show that finding the optimal bounds for a given number of modes is equivalent to a geometric problem: finding the extremal rays of a spectrahedron. We show how this is done analytically for simple cases, and numerically formulated as semidefinite programming (SDP) problems for more complicated cases. We demonstrate this approach with a number of well-motivated examples in particle physics and cosmology, including EFTs of scalars, vectors, fermions and gravitons. In all these cases, we find that the SDP approach leads to results that either improve the previous ones or are completely new. We also find that the SDP approach is numerically much more efficient.


title: Lyman-$\alpha$ emission from a WISE-selected optically faint powerful radio galaxy M151304.72-252439.7 at $z$ = 3.132

abstract: We report the detection of a large ($\sim90$ kpc) and luminous $\mathrm{Ly\alpha}$ nebula [$L\mathrm{_{Ly\alpha}}$ = $(6.80\pm0.08)\times 10^{44}$] $\rm{\,erg\,s^{-1}}$ around an optically faint (r$>23$ mag) radio galaxy M1513-2524 at $z\mathrm{_{em}}$=3.132. The double-lobed radio emission has an extent of 184 kpc, but the radio core, i.e., emission associated with the active galactic nucleus (AGN) itself, is barely detected. This object was found as part of our survey to identify high-$z$ quasars based on Wide-field Infrared Survey Explorer (WISE) colors. The optical spectrum has revealed $\mathrm{Ly\alpha}$, NV, CIV and HeII emission lines with a very weak continuum. Based on long-slit spectroscopy and narrow band imaging centered on the $\mathrm{Ly\alpha}$ emission, we identify two spatial components: a "compact component" with high velocity dispersion ($\sim 1500$$\rm{\,km\,s^{-1}}$) seen in all three lines, and an "extended component", having low velocity dispersion (i.e., 700-1000$\rm{\,km\,s^{-1}}$). The emission line ratios are consistent with the compact component being in photoionization equilibrium with an AGN. We also detect spatially extended associated $\mathrm{Ly\alpha}$ absorption, which is blue-shifted within 250-400$\rm{\,km\,s^{-1}}$ of the $\mathrm{Ly\alpha}$ peak. The probability of $\mathrm{Ly\alpha}$ absorption detection in such large radio sources is found to be low ($\sim$10%) in the literature. M1513-2524 belongs to the top few percent of the population in terms of $\mathrm{Ly\alpha}$ and radio luminosities. Deep integral field spectroscopy is essential for probing this interesting source and its surroundings in more detail.


title: The growth of H II regions around massive stars: the role of metallicity and dust

abstract: Gas metallicity $Z$ and the related dust-to-gas ratio $f_\textrm{d}$ can influence the growth of HII regions via metal line cooling and UV absorption. We model these effects in star-forming regions containing massive stars. We compute stellar feedback from photoionization and radiation pressure (RP) using Monte Carlo radiative transfer coupled with hydrodynamics, including stellar and diffuse radiation fields. We follow a $10^5$ M$_\odot$ turbulent cloud with $Z/$Z$_\odot$ = 2, 1, 0.5, 0.1 and $f_\textrm{d} = 0.01 Z/$Z$_\odot$ with a cluster-sink particle method for star formation. The models evolve for at least 1.5 Myr under feedback. Lower $Z$ results in higher temperatures and therefore larger HII regions. For $Z \ge$Z$_\odot$, radiation pressure $P_\textrm{rad}$ can dominate locally over the gas pressure $P_\textrm{gas}$ in the inner half-parsec around sink particles. Globally, the ratio of $P_\textrm{rad}/P_\textrm{gas}$ is around 1 (2 Z$_\odot$), 0.3 (Z$_\odot$), 0.1 (0.5 Z$_\odot$), and 0.03 (0.1 Z$_\odot$). In the solar model, excluding RP results in an ionized volume several times smaller than the fiducial model with both mechanisms. Excluding RP and UV attenuation by dust results in a larger ionized volume than the fiducial case. That is, UV absorption hinders growth more than RP helps it. The radial expansion velocity of ionized gas reaches $+$15 km/s outwards, while neutral gas has inward velocities for most of the runtime, except for 0.1 Z$_\odot$ which exceeds $+$4 km/s. $Z$ and $f_\textrm{d}$ do not significantly alter the star formation efficiency, rate, or cluster half-mass radius, with the exception of 0.1 Z$_\odot$ due to the earlier expulsion of neutral gas.


title: Searching for energy-resolved quasi-periodic oscillations in AGN

abstract: X-ray quasi-periodic oscillations (QPOs) in AGN allow us to probe and understand the nature of accretion in highly curved space-time, yet the most robust form of detection (i.e. repeat detections over multiple observations) has been limited to a single source to-date, with only tentative claims of single observation detections in several others. The association of those established AGN QPOs with a specific spectral component has motivated us to search the XMM-Newton archive and analyse the energy-resolved lightcurves of 38 bright AGN. We apply a conservative false alarm testing routine folding in the uncertainty and covariance of the underlying broad-band noise. We also explore the impact of red-noise leak and the assumption of various different forms (power-law, broken power-law and lorentzians) for the underlying broad-band noise. In this initial study, we report QPO candidates in 6 AGN (7 including one tentative detection in MRK~766) from our sample of 38, which tend to be found at characteristic energies and, in four cases, at the same frequency across at least two observations, indicating they are highly unlikely to be spurious in nature.


title: The rapid transition from star-formation to AGN dominated rest-frame UV light at z ~ 4

abstract: With the advent of deep optical-to-near-infrared extragalactic imaging on the degree scale, samples of high-redshift sources are being selected that contain both bright star-forming (SF) galaxies and faint active galactic nuclei (AGN). In this study we investigate the transition between SF and AGN-dominated systems at $z \simeq 4$ in the rest-frame UV. We find a rapid transition to AGN-dominated sources bright-ward of $M_{\rm UV} \simeq -23.2$. The effect is observed in the rest-frame UV morphology and size-luminosity relation, where extended clumpy systems become point-source dominated, and also in the available spectra for the sample. These results allow us to derive the rest-frame UV luminosity function for the SF and AGN-dominated sub-samples. We find the SF-dominated LF is best fit with a double-power law, with a lensed Schechter function being unable to explain the existence of extremely luminous SF galaxies at $M_{\rm UV} \simeq -23.5$. If we identify AGN-dominated sources according to a point-source morphology criterion we recover the relatively flat faint-end slope of the AGN LF determined in previous studies. If we instead separate the LF according to the current spectroscopic AGN fraction, we find a steeper faint-end slope of $\alpha = -1.83 \pm 0.11$. Using a simple model to predict the rest-frame AGN LF from the $z = 4 $ galaxy LF we find that the increasing impact of host galaxy light on the measured morphology of faint AGN can explain our observations.


title: Dating individual quasars with the HeII proximity effect

abstract: Constraints on the time-scales of quasar activity are key to understanding the formation and growth of supermassive black holes (SMBHs), quasar triggering mechanisms, and possible feedback effects on their host galaxies. However, observational estimates of this so-called quasar lifetime are highly uncertain (t_Q~10^4-10^9 yr), because most methods are indirect and involve many model-dependent assumptions. Direct evidence of earlier activity is gained from the higher ionization state of the intergalactic medium (IGM) in the quasar environs, observable as enhanced Ly$\alpha$ transmission in the so-called proximity zone. Due to the ~30 Myr equilibration time-scale of HeII in the z~3 IGM, the size of the HeII proximity zone depends on the time the quasar had been active before our observation t_on<t_Q, enabling up to $\pm$0.2 dex precise measurements of individual quasar on-times that are comparable to the e-folding time-scale t_S~44 Myr of SMBH growth. Here we present the first statistical sample of 13 quasars whose accurate and precise systemic redshifts allow for measurements of sufficiently precise HeII quasar proximity zone sizes between ~2 and ~15 proper Mpc from science-grade Hubble Space Telescope (HST) spectra. Comparing these sizes to predictions from cosmological hydrodynamical simulations post-processed with one-dimensional radiative transfer, we infer a broad range of quasar on-times from t_on<1 Myr to t_on>30 Myr that does not depend on quasar luminosity, black hole mass, or Eddington ratio. These results point to episodic quasar activity over a long duty cycle, but do not rule out substantial SMBH growth during phases of radiative inefficiency or obscuration.


title: Stromlo Stellar Tracks: non-solar scaled abundances for massive stars

abstract: We present the Stromlo Stellar Tracks, a set of stellar evolutionary tracks, computed by modifying the Modules for Experiments in Stellar Astrophysics (MESA) 1D stellar evolution package, to fit the Galactic Concordance abundances for hot ($T > 8000$ K) massive ($\geq 10M_\odot$) Main-Sequence (MS) stars. Until now, all stellar evolution tracks are computed at solar, scaled-solar, or alpha-element enhanced abundances, and none of these models correctly represent the Galactic Concordance abundances at different metallicities. This paper is the first implementation of Galactic Concordance abundances to the stellar evolution models. The Stromlo tracks cover massive stars ($10\leq M/M_\odot \leq 300$) with varying rotations ($v/v_{\rm crit} = 0.0, 0.2, 0.4$) and a finely sampled grid of metallicities ($-2.0 \leq {\rm [Z/H]} \leq +0.5$; $\Delta {\rm [Z/H]} = 0.1$) evolved from the pre-main sequence to the end of $^{12}$Carbon burning. We find that the implementation of Galactic Concordance abundances is critical for the evolution of main-sequence, massive hot stars in order to estimate accurate stellar outputs (L, T$_{\rm eff}$, $g$), which, in turn, have a significant impact on determining the ionizing photon luminosity budgets. We additionally support prior findings of the importance that rotation plays on the evolution of massive stars and their ionizing budget. The evolutionary tracks for our Galactic Concordance abundance scaling provide a more empirically motivated approach than simple uniform abundance scaling with metallicity for the analysis of HII regions and have considerable implications in determining nebular emission lines and metallicity. Therefore, it is important to refine the existing stellar evolutionary models for comprehensive high-redshift extragalactic studies. The Stromlo tracks are publicly available to the astronomical community online.


title: Mass accretion rates and multi-scale halo environment in cold and warm dark matter cosmologies

abstract: We study the evolving environment dependence of mass accretion by dark haloes in simulations of cold and warm dark matter (CDM and WDM) cosmologies. The latter allows us to probe the nature of halo growth at scales below the WDM half-mode mass, which form an extreme regime of nonlinear collisionless dynamics and offer an excellent test-bed for ideas relating to hierarchical growth. As environmental proxies, we use the local halo-centric matter density $\delta$ and tidal anisotropy $\alpha$, as well as large-scale halo bias $b_1$. Our analysis, while reproducing known trends for environment-dependent accretion in CDM, as well as the comparison between accretion in CDM and WDM, reveals several interesting new features. As expected from excursion set models, WDM haloes have higher specific accretion rates, dominated by the accretion of diffuse mass, as compared to CDM haloes. For low-mass WDM haloes, we find that the environment-dependence of both diffuse mass accretion as well as accretion by mergers is almost fully explained by $\alpha$. For the other cases, $\delta$ plays at least a comparable role. We detect, for the first time, a significant and evolving assembly bias due to diffuse mass accretion for low-mass CDM and WDM haloes (after excluding splashback objects), with a $z=0$ strength higher than with almost all known secondary variables and largely explained by $\alpha$. Our results place constraints on semi-analytical merger tree algorithms, which in turn could affect the predictions of galaxy evolution models based on them.


title: Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) XII. Extended [C II] Structure (Merger or Outflow) in a z = 6.72 Red Quasar

abstract: We present ALMA [C II] 158 $\mu$m line and far-infrared (FIR) continuum emission observations toward HSC J120505.09$-$000027.9 (J1205$-$0000) at $z = 6.72$ with the beam size of $\sim 0''.8 \times 0''.5$ (or 4.1 kpc $\times$ 2.6 kpc), the most distant red quasar known to date. Red quasars are modestly reddened by dust, and are thought to be in rapid transition from an obscured starburst to an unobscured normal quasar, driven by powerful active galactic nucleus (AGN) feedback which blows out a cocoon of interstellar medium (ISM). The FIR continuum of J1205$-$0000 is bright, with an estimated luminosity of $L_{\rm FIR} \sim 3 \times 10^{12}~L_\odot$. The [C II] line emission is extended on scales of $r \sim 5$ kpc, greater than the FIR continuum. The line profiles at the extended regions are complex and broad (FWHM $\sim 630-780$ km s$^{-1}$). Although it is not practical to identify the nature of this extended structure, possible explanations include (i) companion/merging galaxies and (ii) massive AGN-driven outflows. For the case of (i), the companions are modestly star-forming ($\sim 10~M_\odot$ yr$^{-1}$), but are not detected by our Subaru optical observations ($y_{\rm AB,5\sigma} = 24.4$ mag). For the case of (ii), our lower-limit to the cold neutral outflow rate is $\sim 100~M_\odot$ yr$^{-1}$. The outflow kinetic energy and momentum are both much smaller than what predicted in energy-conserving wind models, suggesting that the AGN feedback in this quasar is not capable of completely suppressing its star formation.


title: Modeling compact binary signals and instrumental glitches in gravitational wave data

abstract: Transient non-gaussian noise in gravitational wave detectors, commonly referred to as glitches, pose challenges for inference of the astrophysical properties of detected signals when the two are coincident in time. Current analyses aim towards modeling and subtracting the glitches from the data using a flexible, morphology-independent model in terms of sine-gaussian wavelets before the signal source properties are inferred using templates for the compact binary signal. We present a new analysis of gravitational wave data that contain both a signal and glitches by simultaneously modeling the compact binary signal in terms of templates and the instrumental glitches using sine-gaussian wavelets. The model for the glitches is generic and can thus be applied to a wide range of glitch morphologies without any special tuning. The simultaneous modeling of the astrophysical signal with templates allows us to efficiently separate the signal from the glitches, as we demonstrate using simulated signals injected around real O2 glitches in the two LIGO detectors. We show that our new proposed analysis can separate overlapping glitches and signals, estimate the compact binary parameters, and provide ready-to-use glitch-subtracted data for downstream inference analyses.


title: AT2017gfo: Bayesian inference and model selection of multi-component kilonovae and constraints on the neutron star equation of state

abstract: The joint detection of the gravitational wave GW170817, of the short $\gamma$-ray burst GRB170817A and of the kilonova AT2017gfo, generated by the the binary neutron star merger observed on August 17, 2017, is a milestone in multimessenger astronomy and provides new constraints on the neutron star equation of state. We perform Bayesian inference and model selection on AT2017gfo using semi-analytical, multi-components models that also account for non-spherical ejecta. Observational data favor anisotropic geometries to spherically symmetric profiles, with a log-Bayes' factor of ${\sim}10^{4}$, and favor multi-component models against single-component ones. The best fitting model is an anisotropic three-component composed of dynamical ejecta plus neutrino and viscous winds. Using the dynamical ejecta parameters inferred from the best-fitting model and numerical-relativity relations connecting the ejecta properties to the binary properties, we constrain the binary mass ratio to $q<1.54$ and the reduced tidal parameter to $120<\tilde\Lambda<1110$. Finally, we combine the predictions from AT2017gfo with those from GW170817, constraining the radius of a neutron star of $1.4~{\rm M}_\odot$ to $12.2\pm0.5~{\rm km}$ ($1\sigma$ level). This prediction could be further strengthened by improving kilonova models with numerical-relativity information.


title: Following up the Kepler field: Masses of Targets for transit timing and atmospheric characterization

abstract: We identify a set of planetary systems observed by Kepler that merit transit timing variation (TTV) analysis given the orbital periods of transiting planets, the uncertainties for their transit times and the number of transits observed during the Kepler mission. We confirm the planetary nature of 4 KOIs within multicandidate systems. We forward model each of the planetary systems identified to determine which systems are likely to yield mass constraints that may be significantly improved upon with follow-up transit observations. We find projected TTVs diverge by more than 90 minutes after 6000 days in 27 systems, including 22 planets with orbital periods exceeding 25 days. Such targets would benefit the most from additional transit timing data. TTV follow-up could push exoplanet characterization to lower masses, at greater orbital periods and at cooler equilibrium temperatures than is currently possible from the Kepler dataset alone. Combining TTVs and recently revised stellar parameters, we characterize an ensemble of homogeneously selected planets and identify planets in the Kepler field with large enough estimated transmission annuli for atmospheric characterization with JWST.


title: Boson stars and solitons confined in a Minkowski box

abstract: We consider the static charged black hole bomb system, originally designed for a (uncharged) rotating superradiant system by Press and Teukolsky. A charged scalar field confined in a Minkowski cavity with a Maxwell gauge field has a quantized spectrum of normal modes that can fit inside the box. Back-reacting non-linearly these normal modes, we find the hairy solitons, a.k.a boson stars (depending on the chosen $U(1)$ gauge), of the theory. The scalar condensate is totally confined inside the box and, outside it, we have the Reissner-Nordstr\"om solution. The Israel junction conditions at the box surface layer determine the stress tensor that the box must have to confine the scalar hair. Some of these horizonless hairy solutions exist for any value of the scalar field charge and not only above the natural critical charges of the theory (namely, the critical charges for the onset of the near-horizon and superradiant instabilities of the Reissner-Nordstr\"om black hole). However, the ground state solutions have a non-trivial intricate phase diagram with a main and a secondary family of solitons (some with a Chandrasekhar mass limit but others without) and there are a third and a fourth critical scalar field charges where the soliton spectra changes radically. Most of these intricate properties are not captured by a higher order perturbative analysis of the problem where we simply back-react a normal mode of the system.


title: Stochastic Optimization for Vaccine and Testing Kit Allocation for the COVID-19 Pandemic

abstract: The pandemic caused by the SARS-CoV-2 virus has exposed many flaws in the decision-making strategies used to distribute resources to combat global health crises. In this paper, we leverage reinforcement learning and optimization to improve upon the allocation strategies for various resources. In particular, we consider a problem where a central controller must decide where to send testing kits to learn about the uncertain states of the world (active learning); then, use the new information to construct beliefs about the states and decide where to allocate resources. We propose a general model coupled with a tunable lookahead policy for making vaccine allocation decisions without perfect knowledge about the state of the world. The lookahead policy is compared to a population-based myopic policy which is more likely to be similar to the present strategies in practice. Each vaccine allocation policy works in conjunction with a testing kit allocation policy to perform active learning. Our simulation results demonstrate that an optimization-based lookahead decision making strategy will outperform the presented myopic policy.


title: The Evolution of the Lyman-Alpha Luminosity Function During Reionization

abstract: The time frame in which hydrogen reionization occurred is highly uncertain, but can be constrained by observations of Lyman-alpha (Ly$\alpha$) emission from distant sources. Neutral hydrogen in the intergalactic medium (IGM) attenuates Ly$\alpha$~photons emitted by galaxies. As reionization progressed the IGM opacity decreased, increasing Ly$\alpha$~visibility. The galaxy Ly$\alpha$~luminosity function (LF) is thus a useful tool to constrain the timeline of reionization. In this work, we model the Ly$\alpha$~LF as a function of redshift, $z=5-10$, and average IGM neutral hydrogen fraction, $\overline{x}_\textsc{hi}$. We combine the Ly$\alpha$~luminosity probability distribution obtained from inhomogeneous reionization simulations with a model for the UV LF to model the Ly$\alpha$~LF. As the neutral fraction increases, the average number density of Ly$\alpha$~emitting galaxies decreases, and are less luminous, though for $\overline{x}_\textsc{hi} \lesssim 0.4$ there is only a small decrease of the Ly$\alpha$~LF. We use our model to infer the IGM neutral fraction at $z=6.6, 7.0, 7.3$ from observed Ly$\alpha$~LFs. We conclude that there is a significant increase in the neutral fraction with increasing redshift: $\overline{x}_\textsc{hi}(z=6.6)=0.08^{+ 0.08}_{- 0.05}, \, \overline{x}_\textsc{hi}(z=7.0)=0.28 \pm 0.05$ and $\overline{x}_\textsc{hi}(z=7.3)=0.83^{+ 0.06}_{- 0.07}$. We predict trends in the Ly$\alpha$~luminosity density and Schechter parameters as a function of redshift and the neutral fraction. We find that the Ly$\alpha$~luminosity density decreases as the universe becomes more neutral. Furthermore, as the neutral fraction increases, the faint-end slope of the Ly$\alpha$~LF steepens, and the characteristic Ly$\alpha$~luminosity shifts to lower values, concluding that the evolving shape of the Ly$\alpha$~LF -- not just its integral -- is an important tool to study reionization.


title: Conformal upper bounds for the volume spectrum

abstract: In this paper, we prove upper bounds for the volume spectrum of a Riemannian manifold that depend only on the volume, dimension and a conformal invariant.


title: Semantic Video Segmentation for Intracytoplasmic Sperm Injection Procedures

abstract: We present the first deep learning model for the analysis of intracytoplasmic sperm injection (ICSI) procedures. Using a dataset of ICSI procedure videos, we train a deep neural network to segment key objects in the videos achieving a mean IoU of 0.962, and to localize the needle tip achieving a mean pixel error of 3.793 pixels at 14 FPS on a single GPU. We further analyze the variation between the dataset's human annotators and find the model's performance to be comparable to human experts.


title: An implicit constitutive relation in which the stress and the linearized strain appear linearly, for describing the small displacement gradient response of elastic solids

abstract: In this short note we develop a constitutive relation that is linear in both the Cauchy stress and the linearized strain, by linearizing implicit constitutive relations between the stress and the deformation gradient that have been put into place to describe the response of elastic bodies (see \cite{rajagopal2003implicit}), by assuming that the displacement gradient is small. These implicit equations include the classical linearized elastic constitutive approximation as well as constitutive relations that imply limiting strain, as special subclasses.


title: A New Method for Simulating Photoprocesses in Astrochemical Models

abstract: We propose a new model for treating solid-phase photoprocesses in interstellar ice analogues. In this approach, photoionization and photoexcitation are included in more detail, and the production of electronically-excited (suprathermal) species is explicitly considered. In addition, we have included non-thermal, non-diffusive chemistry to account for the low-temperature characteristic of cold cores. As an initial test of our method, we have simulated two previous experimental studies involving the UV irradiation of pure solid O$_2$. In contrast to previous solid-state astrochemical model calculations which have used gas-phase photoabsorption cross-sections, we have employed solid-state cross-sections in our calculations. This method allows the model to be tested using well-constrained experiments rather than poorly constrained gas-phase abundances in ISM regions. Our results indicate that inclusion of non-thermal reactions and suprathermal species allows for reproduction of low-temperature solid-phase photoprocessing that simulate interstellar ices within cold ($\sim$ 10 K) dense cores such as TMC-1.

categories: astro-ph.GA cond-mat.mtrl-sci physics.chem-ph physics.comp-ph physics.space-ph


title: The hot spots conjecture can be false: Some numerical examples

abstract: The hot spots conjecture is only known to be true for special geometries. It can be shown numerically that the hot spots conjecture can fail to be true for easy to construct bounded domains with one hole. The underlying eigenvalue problem for the Laplace equation with Neumann boundary condition is solved with boundary integral equations yielding a non-linear eigenvalue problem. Its discretization via the boundary element collocation method in combination with the algorithm by Beyn yields highly accurate results both for the first non-zero eigenvalue and its corresponding eigenfunction which is due to superconvergence. Additionally, it can be shown numerically that the ratio between the maximal/minimal value inside the domain and its maximal/minimal value on the boundary can be larger than $1+10^{-3}$. Finally, numerical examples for easy to construct domains with up to five holes are provided which fail the hot spots conjecture as well.


title: Surgery on Aut(F2)

abstract: We study a geometric construction of certain finite index subgroups of Aut(F2).


title: Machine Learning for User Partitioning and Phase Shifters Design in RIS-Aided NOMA Networks

abstract: A novel reconfigurable intelligent surface (RIS) aided non-orthogonal multiple access (NOMA) downlink transmission framework is proposed. We formulate a long-term stochastic optimization problem that involves a joint optimization of NOMA user partitioning and RIS phase shifting, aiming at maximizing the sum data rate of the mobile users (MUs) in NOMA downlink networks. To solve the challenging joint optimization problem, we invoke a modified object migration automation (MOMA) algorithm to partition the users into equal-size clusters. To optimize the RIS phase-shifting matrix, we propose a deep deterministic policy gradient (DDPG) algorithm to collaboratively control multiple reflecting elements (REs) of the RIS. Different from conventional training-then-testing processing, we consider a long-term self-adjusting learning model where the intelligent agent is capable of learning the optimal action for every given state through exploration and exploitation. Extensive numerical results demonstrate that: 1) The proposed RIS-aided NOMA downlink framework achieves an enhanced sum data rate compared with the conventional orthogonal multiple access (OMA) framework. 2) The proposed DDPG algorithm is capable of learning a dynamic resource allocation policy in a long-term manner. 3) The performance of the proposed RIS-aided NOMA framework can be improved by increasing the granularity of the RIS phase shifts. The numerical results also show that reducing the granularity of the RIS phase shifts and increasing the number of REs are two efficient methods to improve the sum data rate of the MUs.


title: Improving Portuguese Semantic Role Labeling with Transformers and Transfer Learning

abstract: The Natural Language Processing task of determining "Who did what to whom" is called Semantic Role Labeling. For English, recent methods based on Transformer models have allowed for major improvements in this task over the previous state of the art. However, for low resource languages, like Portuguese, currently available semantic role labeling models are hindered by scarce training data. In this paper, we explore a model architecture with only a pre-trained Transformer-based model, a linear layer, softmax and Viterbi decoding. We substantially improve the state-of-the-art performance in Portuguese by over 15 F1. Additionally, we improve semantic role labeling results in Portuguese corpora by exploiting cross-lingual transfer learning using multilingual pre-trained models, and transfer learning from dependency parsing in Portuguese, evaluating the various proposed approaches empirically.


title: Reconstructing Patchy Reionization with Deep Learning

abstract: The precision anticipated from next-generation cosmic microwave background (CMB) surveys will create opportunities for characteristically new insights into cosmology. Secondary anisotropies of the CMB will have an increased importance in forthcoming surveys, due both to the cosmological information they encode and the role they play in obscuring our view of the primary fluctuations. Quadratic estimators have become the standard tools for reconstructing the fields that distort the primary CMB and produce secondary anisotropies. While successful for lensing reconstruction with current data, quadratic estimators will be sub-optimal for the reconstruction of lensing and other effects at the expected sensitivity of the upcoming CMB surveys. In this paper we describe a convolutional neural network, ResUNet-CMB, that is capable of the simultaneous reconstruction of two sources of secondary CMB anisotropies, gravitational lensing and patchy reionization. We show that the ResUNet-CMB network significantly outperforms the quadratic estimator at low noise levels and is not subject to the lensing-induced bias on the patchy reionization reconstruction that would be present with a straightforward application of the quadratic estimator.


title: Learn by Guessing: Multi-Step Pseudo-Label Refinement for Person Re-Identification

abstract: Unsupervised Domain Adaptation (UDA) methods for person Re-Identification (Re-ID) rely on target domain samples to model the marginal distribution of the data. To deal with the lack of target domain labels, UDA methods leverage information from labeled source samples and unlabeled target samples. A promising approach relies on the use of unsupervised learning as part of the pipeline, such as clustering methods. The quality of the clusters clearly plays a major role in methods performance, but this point has been overlooked. In this work, we propose a multi-step pseudo-label refinement method to select the best possible clusters and keep improving them so that these clusters become closer to the class divisions without knowledge of the class labels. Our refinement method includes a cluster selection strategy and a camera-based normalization method which reduces the within-domain variations caused by the use of multiple cameras in person Re-ID. This allows our method to reach state-of-the-art UDA results on DukeMTMC-Market1501 (source-target). We surpass state-of-the-art for UDA Re-ID by 3.4% on Market1501-DukeMTMC datasets, which is a more challenging adaptation setup because the target domain (DukeMTMC) has eight distinct cameras. Furthermore, the camera-based normalization method causes a significant reduction in the number of iterations required for training convergence.


title: First Integrals and symmetries of nonholonomic systems

abstract: In nonholonomic mechanics, the presence of constraints in the velocities breaks the well-under\-stood link between symmetries and first integrals of holonomic systems, expressed in Noether's Theorem. However there is a known special class of first integrals of nonholonomic systems generated by vector fields tangent to the group orbits, called {\it horizontal gauge momenta}, that suggest that some version of this link should still hold. In this paper we prove that, under certain conditions on the symmetry Lie group, the (nonholonomic) momentum map is conserved along the nonholonomic dynamics, thus extending Noether Theorem to the nonholonomic framework. Our analysis leads to a constructive method, with fundamental consequences to the integrability of some nonholonomic systems as well as their hamiltonization. We apply our results to three paradigmatic examples: the snakeboard, a solid of revolution rolling without sliding on a plane and a heavy homogeneous ball that rolls without sliding inside a convex surface of revolution. In particular, for the snakeboard we show the existence of a new horizontal gauge momentum that reveals new aspects of its integrability.


title: Generative Adversarial Networks with Physical Evaluators for Spray Simulation of Pintle Injector

abstract: Due to the adjustable geometry, pintle injectors are specially suitable for the liquid rocket engines which require a widely throttleable range. While applying the conventional computational fluid dynamics approaches to simulate the complex spray phenomena in the whole range still remains to be a great challenge. In this paper, a novel deep learning approach used to simulate instantaneous spray fields under continuous operating conditions is explored. Based on one specific type of neural networks and the idea of physics constraint, a Generative Adversarial Networks with Physics Evaluators (GAN-PE) framework is proposed. The geometry design and mass flux information are embedded as the inputs. After the adversarial training between the generator and discriminator, the generated field solutions are fed into the two physics evaluators. In this framework, mass conversation evaluator is designed to improve the training robustness and convergence. And the spray angle evaluator, which is composed of a down-sampling CNN and theoretical model, guides the networks generating the spray solutions more according with the injection conditions. The characterization of the simulated spray, including the spray morphology, droplet distribution and spray angle, is well predicted. The work suggests a great potential of the prior physics knowledge employment in the simulation of instantaneous flow fields.


title: Proper splittings and reduced solutions of matrix equations

abstract: In this article we apply proper splittings of matrices to develop an iterative process to approximate solutions of matrix equations of the form TX = W. Moreover, by using the partial order induced by positive semidefinite matrices, we obtain equivalent conditions to the convergence of this process. We also include some speed comparison results of the convergence of this method. In addition, for all matrix T we propose a proper splitting based on the polar decomposition of T.


title: Observing the host galaxies of high-redshift quasars with JWST: predictions from the BlueTides simulation

abstract: The bright emission from high-redshift quasars completely conceals their host galaxies in the rest-frame ultraviolet/optical, with detection of the hosts in these wavelengths eluding even the Hubble Space Telescope (HST) using detailed point spread function (PSF) modelling techniques. In this study we produce mock images of a sample of z=7 quasars extracted from the BlueTides simulation, and apply Markov Chain Monte Carlo-based PSF modelling to determine the detectability of their host galaxies with the James Webb Space Telescope (JWST). While no statistically significant detections are made with HST, we predict that at the same wavelengths and exposure times JWST NIRCam imaging will detect ~50% of quasar host galaxies. We investigate various observational strategies, and find that NIRCam wide-band imaging in the long-wavelength filters results in the highest fraction of successful quasar host detections, detecting >80% of the hosts of bright quasars in exposure times of 5 ks. Exposure times of ~5 ks are required to detect the majority of host galaxies in the NIRCam wide-band filters, however even 10 ks exposures with MIRI result in <30% successful host detections. We find no significant trends between galaxy properties and their detectability. The PSF modelling can accurately recover the host magnitudes, radii, and spatial distribution of the larger-scale emission, when accounting for the central core being contaminated by residual quasar flux. Care should be made when interpreting the host properties measured using PSF modelling.


title: Forward Modeling Helioseismic Signatures of One- and Two-cell Meridional Circulation

abstract: Using a 3D global solver of the linearized Euler equations, we model acoustic oscillations over background velocity flow fields of single-cell meridional circulation with deep and shallow return flows as well as double-cell meridional circulation with strong and weak reversals. The velocities are generated using a mean-field hydrodynamic and dynamo model -- moving through the regimes with minimal parameter changes; counter-rotation near the base of the tachocline is induced by sign inversion of the non-diffusive action of turbulent Reynolds stresses ($\Lambda$-effect) due to the radial inhomogeneity of the Coriolis number. By mimicking the stochastic excitation of resonant modes in the convective interior, we simulate realization noise present in solar observations. Using deep-focusing to analyze differences in travel-time signatures between the four regimes, as well as comparing to solar observations, we show that current helioseismology techniques may offer important insights about the location and strength of the return flow, however, that it may not currently be possible to definitively distinguish between profiles of single-cell or double-cell meridional circulation.


title: Importance of magnetic fields in highly eccentric discs with applications to tidal disruption events

abstract: Whether tidal disruption events (TDEs) circularise or accrete directly as a highly eccentric disc is the subject of current research and appears to depend sensitively on the disc thermodynamics. In a previous paper we applied the theory of eccentric discs to TDE discs using an $\alpha-$prescription for the disc stress, which leads to solutions that exhibit extreme, potentially unphysical, behaviour. In this paper we further explore the dynamical vertical structure of highly eccentric discs using alternative stress models that are better motivated by the behaviour of magnetic fields in eccentric discs. We find that the presence of a coherent magnetic field has a stabilising effect on the dynamics and can significantly alter the behaviour of highly eccentric radiation dominated discs. We conclude that magnetic fields are important for the evolution of TDE discs.


title: Dynamics and Rheology of Ring-Linear Blend Semidilute Solutions in Extensional Flow: Single Molecule Experiments

abstract: Ring polymers exhibit unique flow properties due to their closed chain topology. Despite recent progress, we have not yet achieved a full understanding of the nonequilibrium flow behavior of rings in nondilute solutions where intermolecular interactions greatly influence chain dynamics. In this work, we directly observe the dynamics of DNA rings in semidilute ring-linear polymer blends using single molecule techniques. We systematically investigate ring polymer relaxation dynamics from high extension and transient and steady-state stretching dynamics in planar extensional flow for a series of ring-linear blends with varying ring fraction. Our results show multiple molecular sub-populations for ring relaxation in ring-linear blends, as well as large conformational fluctuations for rings in steady extensional flow, even long after the initial transient stretching process has subsided. We further quantify the magnitude and characteristic timescales of ring conformational fluctuations as a function of blend composition. Interestingly, we find that the magnitude of ring conformational fluctuations follows a non-monotonic response with increasing ring fraction, first increasing at low ring fraction and then substantially decreasing at large ring fraction in ring-linear blends. A unique set of ring polymer conformations are observed during the transient stretching process, which highlights the prevalence of molecular individualism and supports the notion of complex intermolecular interactions in ring-linear polymer blends. Together with results from molecular simulations, our results suggest that ring conformational fluctuations arise due to ring-linear threading and intermolecular hydrodynamic interactions (HI). Taken together, our results provide a new molecular understanding of ring polymer dynamics in ring-linear blends in nonequilibrium flow.


title: The damped harmonic oscillator at the classical limit of the Snyder-de Sitter space

abstract: Valtancoli in his paper entitled [P. Valtancoli, Canonical transformations, and minimal length J. Math. Phys. 56, 122107 (2015)] has shown how the deformation of the canonical transformations can be made compatible with the deformed Poisson brackets. Based on this work and through an appropriate canonical transformation, we solve the problem of one dimensional (1D) damped harmonic oscillator at the classical limit of the Snyder-de Sitter (SdS) space. We show that the equations of the motion can be described by trigonometric functions with frequency and period depending on the deformed and the damped parameters. We eventually discuss the influences of these parameters on the motion of the system.


title: Detecting a network of hijacked journals by its archive

abstract: This study describes a method to detect hijacked journals based on the analysis of the archives of clone journals. This approach is most effective in discovering a network of hijacked journals that have the same organizer(s). Analysis of the archives of clone journals allowed to detect 62 URLs of hijacked journals. It also provided the possibility to predict two clone websites before they became operational. This study shows that most detected hijacked journals represent a network of clone journals organized by one or several fraudulent individuals. The information and content of nine legitimate journals were compromised in international and national scientometric databases.


title: Alkali metals in white dwarf atmospheres as tracers of ancient planetary crusts

abstract: White dwarfs that accrete the debris of tidally disrupted asteroids provide the opportunity to measure the bulk composition of the building blocks, or fragments, of exoplanets. This technique has established a diversity in compositions comparable to what is observed in the solar system, suggesting that the formation of rocky planets is a generic process. Whereas the relative abundances of lithophile and siderophile elements within the planetary debris can be used to investigate whether exoplanets undergo differentiation, the composition studies carried out so far lack unambiguous tracers of planetary crusts. Here we report the detection of lithium in the atmospheres of four cool (<5,000 K) and old (cooling ages 5-10 Gyr) metal-polluted white dwarfs, where one also displays photospheric potassium. The relative abundances of these two elements with respect to sodium and calcium strongly suggest that all four white dwarfs have accreted fragments of planetary crusts. We detect an infrared excess in one of the systems, indicating that accretion from a circumstellar debris disk is on-going. The main-sequence progenitor mass of this star was $4.8\pm0.2 M_\odot$, demonstrating that rocky, differentiated planets may form around short-lived B-type stars.


title: Fast Flavor Oscillations of Astrophysical Neutrinos with $1,\,2,\,\ldots,\,\infty$ Crossings

abstract: In the early Universe, as well as in supernovae and merging neutron stars, neutrinos have such high densities that they affect each other and exhibit collective flavor oscillations. A crucial ingredient for fast collective flavor oscillations is that the electron lepton number (ELN) distribution changes its sign as a function of direction, i.e., has a zero crossing. We present a study in two dimensions and show how fast flavor oscillations depend on the ELN and its crossings. We show that a large number of crossings can inhibit flavor oscillations. This may be a natural self-limiting mechanism that stabilizes the flavor content of the dense neutrino gas in a vast majority of scenarios, especially the early Universe, where the angular distributions for all flavors are very similar and crossings occur mainly due to fluctuations.


title: Ability of unbounded pairs of observers to achieve quantum advantage in random access codes with a single pair of qubits

abstract: Complications in preparing and preserving quantum correlations stimulate recycling of a single quantum resource in information processing and communication tasks multiple times. Here, we consider a scenario involving multiple independent pairs of observers acting with unbiased inputs on a single pair of spatially separated qubits sequentially. In this scenario, we address whether more than one pair of observers can demonstrate quantum advantage in some specific $2 \rightarrow 1$ and $3 \rightarrow 1$ random access codes. Interestingly, we not only address these in the affirmative, but also illustrate that unbounded pairs can exhibit quantum advantage. Furthermore, these results remain valid even when all observers perform suitable projective measurements and an appropriate separable state is initially shared.


title: Reddit Entity Linking Dataset

abstract: We introduce and make publicly available an entity linking dataset from Reddit that contains 17,316 linked entities, each annotated by three human annotators and then grouped into Gold, Silver, and Bronze to indicate inter-annotator agreement. We analyze the different errors and disagreements made by annotators and suggest three types of corrections to the raw data. Finally, we tested existing entity linking models that are trained and tuned on text from non-social media datasets. We find that, although these existing entity linking models perform very well on their original datasets, they perform poorly on this social media dataset. We also show that the majority of these errors can be attributed to poor performance on the mention detection subtask. These results indicate the need for better entity linking models that can be applied to the enormous amount of social media text.


title: A Survey on Embedding Dynamic Graphs

abstract: Embedding static graphs in low-dimensional vector spaces plays a key role in network analytics and inference, supporting applications like node classification, link prediction, and graph visualization. However, many real-world networks present dynamic behavior, including topological evolution, feature evolution, and diffusion. Therefore, several methods for embedding dynamic graphs have been proposed to learn network representations over time, facing novel challenges, such as time-domain modeling, temporal features to be captured, and the temporal granularity to be embedded. In this survey, we overview dynamic graph embedding, discussing its fundamentals and the recent advances developed so far. We introduce the formal definition of dynamic graph embedding, focusing on the problem setting and introducing a novel taxonomy for dynamic graph embedding input and output. We further explore different dynamic behaviors that may be encompassed by embeddings, classifying by topological evolution, feature evolution, and processes on networks. Afterward, we describe existing techniques and propose a taxonomy for dynamic graph embedding techniques based on algorithmic approaches, from matrix and tensor factorization to deep learning, random walks, and temporal point processes. We also elucidate main applications, including dynamic link prediction, anomaly detection, and diffusion prediction, and we further state some promising research directions in the area.


title: Absorption of Microwaves by Random-Anisotropy Magnets

abstract: Microscopic model of the interaction of spins with a microwave field in a random-anisotropy magnet has been developed. Numerical results show that microwave absorption occurs in a broad range of frequencies due to the distribution of Imry-Ma domains on sizes and effective anisotropy. That distribution is also responsible for the weak dependence of the absorption on the damping. At a fixed frequency of the ac-field the spatial pattern of the amplitude of spin oscillations exhibits maxima corresponding to the oscillations of resonant magnetic domains. The dependence of the peak absorption frequency on the ratio of the magnitude of per-spin random anisotropy to the strength of the ferromagnetic exchange agrees with the scaling derived from the Imry-Ma argument. The effect noticeably increases in low-dimensional systems, suggesting that materials comprised of microscopic amorphous leaves and wires can be promising candidates for enhanced microwave absorption.


title: Parallel Scaling of the Regionally-Implicit Discontinuous Galerkin Method with Quasi-Quadrature-Free Matrix Assembly

abstract: In this work we investigate the parallel scalability of the numerical method developed in Guthrey and Rossmanith [The regionally implicit discontinuous Galerkin method: Improving the stability of DG-FEM, SIAM J. Numer. Anal. (2019)]. We develop an implementation of the regionally-implicit discontinuous Galerkin (RIDG) method in DoGPack, which is an open source C++ software package for discontinuous Galerkin methods. Specifically, we develop and test a hybrid OpenMP and MPI parallelized implementation of DoGPack with the goal of exploring the efficiency and scalability of RIDG in comparison to the popular strong stability-preserving Runge-Kutta discontinuous Galerkin (SSP-RKDG) method. We demonstrate that RIDG methods are able to hide communication latency associated with distributed memory parallelism, due to the fact that almost all of the work involved in the method is highly localized to each element, producing a localized prediction for each region. We demonstrate the enhanced efficiency and scalability of the of the RIDG method and compare it to SSP-RKDG methods and show extensibility to very high order schemes. The two-dimensional scaling study is performed on machines at the Institute for Cyber-Enabled Research at Michigan State University, using up to 1440 total cores on Intel(R) Xeon(R) Gold 6148 CPU @ 2.40GHz CPUs. The three dimensional scaling study is performed on Livermore Computing clusters at at Lawrence Livermore National Laboratory, using up to 28672 total cores on Intel Xeon CLX-8276L CPUs with Omni-Path interconnects.


title: Open Fishchain in N=4 Supersymmetric Yang-Mills Theory

abstract: We consider a cusped Wilson line with J insertions of scalar fields in N=4 SYM and prove that in a certain limit the Feynman graphs are integrable to all loop orders. We identify the integrable system as a quantum fishchain with open boundary conditions. The existence of the boundary degrees of freedom results in the boundary reflection operator acting non-trivially on the physical space. We derive the Baxter equation for Q-functions and provide the quantisation condition for the spectrum. This allows us to find the non-perturbative spectrum numerically.

categories: hep-th cond-mat.stat-mech math-ph math.MP nlin.SI


title: Perelman's entropy on ancient Ricci flows

abstract: In [ZY2], the second author proved Perelman's assertion, namely, for an ancient Ricci flow with bounded and nonnegative curvature operator, bounded entropy is equivalent to noncollapsing on all scales. In this paper, we continue this discussion. It turns out that the curvature operator nonnegativity is not a necessary condition, and we need only to assume a consequence of Hamilton's trace Harnack. Furthermore, we show that this condition holds for steady Ricci solitons with nonnegative Ricci curvature.


title: On the Fourier transform of a quantitative trait: Implications for compressive sensing

abstract: This paper explores the genotype-phenotype relationship. It outlines conditions under which the dependence of a quantitative trait on the genome might be predictable, based on measurement of a limited subset of genotypes. It uses the theory of real-valued Boolean functions in a systematic way to translate trait data into the Fourier domain. Important trait features, such as the roughness of the trait landscape or the modularity of a trait have a simple Fourier interpretation. Ruggedness at a gene location corresponds to high sensitivity to mutation, while a modular organization of gene activity reduces such sensitivity. Traits where rugged loci are rare will naturally compress gene data in the Fourier domain, leading to a sparse representation of trait data, concentrated in identifiable, low-level coefficients. This Fourier representation of a trait organizes epistasis in a form which is isometric to the trait data. As Fourier matrices are known to be maximally incoherent with the standard basis, this permits employing compressive sensing techniques to work from data sets that are relatively small -- sometimes even of polynomial size -- compared to the exponentially large sets of possible genomes. This theory provides a theoretical underpinning for systematic use of Boolean function machinery to dissect the dependency of a trait on the genome and environment.


title: An integrated abundance model for estimating county-level prevalence of opioid misuse in Ohio

abstract: Opioid misuse is a national epidemic and a significant drug related threat to the United States. While the scale of the problem is undeniable, estimates of the local prevalence of opioid misuse are lacking, despite their importance to policy-making and resource allocation. This is due, in part, to the challenge of directly measuring opioid misuse at a local level. In this paper, we develop a Bayesian hierarchical spatio-temporal abundance model that integrates indirect county-level data on opioid-related outcomes with state-level survey estimates on prevalence of opioid misuse to estimate the latent county-level prevalence and counts of people who misuse opioids. A simulation study shows that our integrated model accurately recovers the latent counts and prevalence. We apply our model to county-level surveillance data on opioid overdose deaths and treatment admissions from the state of Ohio. Our proposed framework can be applied to other applications of small area estimation for hard to reach populations, which is a common occurrence with many health conditions such as those related to illicit behaviors.


title: Recent advances in superlattice frequency multipliers

abstract: Semiconductor superlattice multipliers have emerged as a nonlinear medium capable to generate radiation in a wide frequency range. This property facilitates the potential of sources suitable for sensing and spectroscopy applications. In this study, we further investigate the consequences on harmonic generation in a superlattice multiplier after excitation by an input signal oscillating at different frequencies. Here we provide a rigorous description of our theoretical model including a semiclassical Boltzmann approach to nonlinear miniband transport and non-equilibrium Greens functions calculations treating scattering processes under forward and reverse bias. To fully exploit the features of this radiation source, we focus on the effects of elastic scattering and systemat-ic imperfections in the superlattice structure which lead to asymmetric current flow.


title: Extremal fluctuations driving the relaxation in glassy energy landscapes

abstract: Cooperative events requiring anomalously large fluctuations are a defining characteristic for the onset of glassy relaxation across many materials. The importance of such intermittent events has been noted in systems as diverse as superconductors, metallic glasses, gels, colloids, and granular piles. Here, we show that prohibiting the attainment of new record-high energy fluctuations -- by explicitly imposing a ``lid'' on the fluctuation spectrum -- impedes further relaxation in the glassy phase. This lid allows us to directly measure the impact of record events on the evolving system in extensive simulations of aging in such vastly distinct glass formers as spin glasses and tapped granular piles. Interpreting our results in terms of a dynamics of records succeeds in explaining the ubiquity of both, the logarithmic decay of the energy and the memory effects encoded in the scaling of two-time correlation functions of aging systems.


title: Searching for evidence of algorithmic randomness and incomputability in the output of quantum random number generators

abstract: Ideal quantum random number generators (QRNGs) can produce algorithmically random and thus incomputable sequences, in contrast to pseudo-random number generators. However, the verification of the presence of algorithmic randomness and incomputability is a nontrivial task. We present the results of a search for algorithmic randomness and incomputability in the output from two different QRNGs, performed by applying tests based on the Solovay-Strassen test of primality and the Chaitin-Schwartz theorem. The first QRNG uses measurements of quantum vacuum fluctuations. The second QRNG is based on polarization measurements on entangled single photons; for this generator, we use looped (and thus highly compressible) strings that also allow us to assess the ability of the tests to detect repeated bit patterns. Compared to a previous search for algorithmic randomness, our study increases statistical power by almost 3 orders of magnitude.


title: Explainable Neural Network-based Modulation Classification via Concept Bottleneck Models

abstract: While RFML is expected to be a key enabler of future wireless standards, a significant challenge to the widespread adoption of RFML techniques is the lack of explainability in deep learning models. This work investigates the use of CB models as a means to provide inherent decision explanations in the context of DL-based AMC. Results show that the proposed approach not only meets the performance of single-network DL-based AMC algorithms, but provides the desired model explainability and shows potential for classifying modulation schemes not seen during training (i.e. zero-shot learning).


title: Multiperiodic orbits from interacting soft spots in cyclically-sheared amorphous solids

abstract: When an amorphous solid is deformed cyclically, it may reach a steady state in which the paths of constituent particles trace out closed loops that repeat in each driving cycle. A remarkable variant has been noticed in simulations where the period of particle motions is a multiple of the period of driving, but the reasons for this behavior have remained unclear. Motivated by mesoscopic features of displacement fields in experiments on jammed solids, we propose and analyze a simple model of interacting soft spots -- locations where particles rearrange under stress and that resemble two-level systems with hysteresis. We show that multiperiodic behavior can arise among just three or more soft spots that interact with each other, but in all cases it requires frustrated interactions, illuminating this otherwise elusive type of interaction. We suggest directions for seeking this signature of frustration in experiments and for achieving it in designed systems.


title: Search for axion-like dark matter using solid-state nuclear magnetic resonance

abstract: We report the results of an experimental search for ultralight axion-like dark matter in the mass range 162 neV to 166 neV. The detection scheme of our Cosmic Axion Spin Precession Experiment (CASPEr) is based on a precision measurement of $^{207}$Pb solid-state nuclear magnetic resonance in a polarized ferroelectric crystal. Axion-like dark matter can exert an oscillating torque on $^{207}$Pb nuclear spins via the electric-dipole moment coupling $g_d$, or via the gradient coupling $g_{\text{aNN}}$. We calibrated the detector and characterized the excitation spectrum and relaxation parameters of the nuclear spin ensemble with pulsed magnetic resonance measurements in a 4.4 T magnetic field. We swept the magnetic field near this value and searched for axion-like dark matter with Compton frequency within a 1 MHz band centered at 39.65 MHz. Our measurements place the upper bounds $|g_d|<9.5\times10^{-4}\,\text{GeV}^{-2}$ and $|g_{\text{aNN}}|<2.8\times10^{-1}\,\text{GeV}^{-1}$ (95% confidence level) in this frequency range. The constraint on $g_d$ corresponds to an upper bound of $1.0\times 10^{-21}\,\text{e}\cdot\text{cm}$ on the amplitude of oscillations of the neutron electric dipole moment, and $4.3\times 10^{-6}$ on the amplitude of oscillations of CP-violating $\theta$ parameter of quantum chromodynamics. Our results demonstrate the feasibility of using solid-state nuclear magnetic resonance to search for axion-like dark matter in the nano-electronvolt mass range.


title: (In)equality distance patterns and embeddability into Hilbert spaces

abstract: We show that compact Riemannian manifolds, regarded as metric spaces with their global geodesic distance, cannot contain a number of rigid structures such as (a) arbitrarily large regular simplices or (b) arbitrarily long sequences of points equidistant from pairs of points preceding them in the sequence. All of this provides evidence that Riemannian metric spaces admit what we term loose embeddings into finite-dimensional Euclidean spaces: continuous maps that preserve both equality as well as inequality. We also prove a local-to-global principle for Riemannian-metric-space loose embeddability: if every finite subspace thereof is loosely embeddable into a common $\mathbb{R}^N$, then the metric space as a whole is loosely embeddable into $\mathbb{R}^N$ in a weakened sense.


title: A Low-Frequency Torsion Pendulum with Interferometric Readout

abstract: We describe a torsion pendulum with a large mass-quadrupole moment and a resonant frequency of 2.8 mHz, whose angle is measured using a modified Michelson interferometer. The system achieved noise levels of $\sim200\ \text{prad}/\sqrt{\text{Hz}}$ between 0.2-30 Hz and $\sim10\ \text{prad}/\sqrt{\text{Hz}}$ above 100 Hz. Such a system can be applied to a broad range of fields from the study of rotational seismic motion and elastogravity signals to gravitational wave observation and tests of gravity.


title: Viewing angle effects in quasar application to cosmology

abstract: The symmetry axes of active galactic nuclei (AGN) are randomly distributed in space but highly inclined sources are heavily obscured and are not seen as quasars with broad emission lines. The obscuring torus geometry determines the average viewing angle, and if the torus geometry changes with the redshift, this average viewing angle will also change. Thus the ratio between the isotropic luminosity and observed luminosity may change systematically with redshift. Therefore, if we use quasars to measure the luminosity distance by evaluating the isotropic absolute luminosity and measuring the observed flux, we can have a redshift-dependent bias which can propagate to cosmological parameters. We propose a toy model for testing the effect of viewing angle uncertainty on measurement of the luminosity distance. The model is based on analytical description of the obscuring torus applied to one-parameter observational data. It illustrates the possible change of the torus covering factor between the two chosen redshift ranges. We have estimated the possible error on specific cosmological parameters (H0,Omega_m) for the flat Lambda-CDM cosmology if a method is calibrated at low redshift and applied to the higher redshift. The errors on cosmological parameters due to potential dependence of viewing angle on redshift are found to be potentially significant, and the effect will have to be accommodated in the future in all quasar-based cosmological methods. A careful systematic study of AGN mean viewing angle across redshift is necessary, with the use of appropriate samples and models which uniquely determine the inclination of each source.


title: Regression Discontinuity Design with Many Thresholds

abstract: Numerous empirical studies employ regression discontinuity designs with multiple cutoffs and heterogeneous treatments. A common practice is to normalize all the cutoffs to zero and estimate one effect. This procedure identifies the average treatment effect (ATE) on the observed distribution of individuals local to existing cutoffs. However, researchers often want to make inferences on more meaningful ATEs, computed over general counterfactual distributions of individuals, rather than simply the observed distribution of individuals local to existing cutoffs. This paper proposes a consistent and asymptotically normal estimator for such ATEs when heterogeneity follows a non-parametric function of cutoff characteristics in the sharp case. The proposed estimator converges at the minimax optimal rate of root-n for a specific choice of tuning parameters. Identification in the fuzzy case, with multiple cutoffs, is impossible unless heterogeneity follows a finite-dimensional function of cutoff characteristics. Under parametric heterogeneity, this paper proposes an ATE estimator for the fuzzy case that optimally combines observations to maximize its precision.


title: Escape and absorption probabilities for obliquely reflected Brownian motion in a quadrant

abstract: We consider an obliquely reflected Brownian motion $Z$ with positive drift in a quadrant stopped at time $T$, where $T:=\inf \{ t>0 : Z(t)=(0,0) \}$ is the first hitting time of the origin. Such a process can be defined even in the non-standard case where the reflection matrix is not completely-$\mathcal{S}$. We show that in this case the process has two possible behaviors: either it tends to infinity or it hits the corner (origin) in a finite time. Given an arbitrary starting point $(u,v)$ in the quadrant, we consider the escape (resp. absorption) probabilities $\mathbb{P}_{(u,v)}[T=\infty]$ (resp. $\mathbb{P}_{(u,v)}[T<\infty]$). We establish the partial differential equations and the oblique Neumann boundary conditions which characterize the escape probability and provide a functional equation satisfied by the Laplace transform of the escape probability. We then give asymptotics for the absorption probability in the simpler case where the starting point in the quadrant is $(u,0)$. We exhibit a remarkable geometric condition on the parameters which characterizes the case where the absorption probability has a product form and is exponential. We call this new criterion the dual skew symmetry condition due to its natural connection with the skew symmetry condition for the stationary distribution. We then obtain an explicit integral expression for the Laplace transform of the escape probability. We conclude by presenting exact asymptotics for the escape probability at the origin.


title: A Block Bidiagonalization Method for Fixed-Accuracy Low-Rank Matrix Approximation

abstract: We present randUBV, a randomized algorithm for matrix sketching based on the block Lanzcos bidiagonalization process. Given a matrix $\bf{A}$, it produces a low-rank approximation of the form ${\bf UBV}^T$, where $\bf{U}$ and $\bf{V}$ have orthonormal columns in exact arithmetic and $\bf{B}$ is block bidiagonal. In finite precision, the columns of both ${\bf U}$ and ${\bf V}$ will be close to orthonormal. Our algorithm is closely related to the randQB algorithms of Yu, Gu, and Li (2018) in that the entries of $\bf{B}$ are incrementally generated and the Frobenius norm approximation error may be efficiently estimated. Our algorithm is therefore suitable for the fixed-accuracy problem, and so is designed to terminate as soon as a user input error tolerance is reached. Numerical experiments suggest that the block Lanczos method is generally competitive with or superior to algorithms that use power iteration, even when $\bf{A}$ has significant clusters of singular values.


title: Rank functions on triangulated categories

abstract: We introduce the notion of a rank function on a triangulated category $\mathcal{C}$ which generalizes the Sylvester rank function in the case when $\mathcal{C}=\operatorname{Perf}(A)$ is the perfect derived category of a ring $A$. We show that rank functions are closely related to functors into simple triangulated categories and classify Verdier quotients into simple triangulated categories in terms of particular rank functions called localizing. If $\mathcal{C}=\operatorname{Perf}(A)$ as above, localizing rank functions also classify finite homological epimorphisms from $A$ into differential graded skew-fields or, more generally, differential graded Artining rings. To establish these results, we develop the theory of derived localization of differential graded algebras at thick subcategories of their perfect derived categories. This is a far-reaching generalization of Cohn's matrix localization of rings and has independent interest.


title: Broadband optical power limiter metalens

abstract: In recent years, the need for high-power laser is of great interest for different applications ranging from direct-laser processing, light detection, medicine, and lighting. However, high-power lasers with high intensities give rise to fundamental problems for optical detectors and imaging systems with low threshold damage, which still need reliable solutions. Here, we report and numerically demonstrate a hybrid system that synergically combines a broadband optical power limiter with a transmittance difference between on-state (70^o C) and off-state (25^o C) about 62.5%, and a diffraction-limited broadband metalens from 1534 nm to 1664 nm. Such a metalens power limiter could be used in any system requiring an intermediate focal plane in the optical path to the detector from damage by exposure to high-intensity lasers.


title: Force probe simulations using an adaptive resolution scheme

abstract: Molecular simulations of the forced unfolding and refolding of biomolecules or molecular complexes allow to gain important kinetic, structural and thermodynamic information about the folding process and the underlying energy landscape. In force probe molecular dynamics (FPMD) simulations, one pulls one end of the molecule with a constant velocity in order to induce the relevant conformational transitions. Since the extended configuration of the system has to fit into the simulation box together with the solvent such simulations are very time consuming. Here, we apply a hybrid scheme in which the solute is treated with atomistic resolution and the solvent molecules far away from the solute are described in a coarse-grained manner. We use the adaptive resolution scheme (AdResS) that has very successfully been applied to various examples of equilibrium simulations. We perform FPMD simulations using AdResS on a well studied system, a dimer formed from mechanically interlocked calixarene capsules. The results of the multiscale simulations are compared to all-atom simulations of the identical system and we observe that the size of the region in which atomistic resolution is required depends on the pulling velocity, i.e. the particular non-equilibrium situation. For large pulling velocities a larger all atom region is required. Our results show that multiscale simulations can be applied also in the strong non-equilibrium situations that the system experiences in FPMD simulations.


title: Robust Maximum Entropy Behavior Cloning

abstract: Imitation learning (IL) algorithms use expert demonstrations to learn a specific task. Most of the existing approaches assume that all expert demonstrations are reliable and trustworthy, but what if there exist some adversarial demonstrations among the given data-set? This may result in poor decision-making performance. We propose a novel general frame-work to directly generate a policy from demonstrations that autonomously detect the adversarial demonstrations and exclude them from the data set. At the same time, it's sample, time-efficient, and does not require a simulator. To model such adversarial demonstration we propose a min-max problem that leverages the entropy of the model to assign weights for each demonstration. This allows us to learn the behavior using only the correct demonstrations or a mixture of correct demonstrations.


title: Probing hadronic interactions with measurements at ultra-high energies with the Pierre Auger Observatory

abstract: The characteristics of an extensive air shower derive from both the mass of the primary ultra-high-energy cosmic ray that seeds its development and the properties of the hadronic interactions that feed it. With its hybrid detector design, the Pierre Auger Observatory measures both the longitudinal development of showers in the atmosphere and the lateral distribution of particles arriving at the ground, from which a number of parameters are calculated and compared with predictions from current hadronic interaction models tuned to LHC data. At present, a tension exists concerning the production of muons, in that the measured abundance exceeds all predictions. This discrepancy, measured up to center-of-mass energies of $\sim$ 140 TeV, is irresolvable through mass composition arguments, constrained by measurements of the depth of the electromagnetic-shower maximum. Here, we discuss a compilation of hadronically-sensitive shower observables and their comparisons with model predictions and conclude with a brief discussion of what measurements with the new detectors of the AugerPrime upgrade will bring to the table.


title: Metropolis-Hastings with Averaged Acceptance Ratios

abstract: Markov chain Monte Carlo (MCMC) methods to sample from a probability distribution $\pi$ defined on a space $(\Theta,\mathcal{T})$ consist of the simulation of realisations of Markov chains $\{\theta_{n},n\geq1\}$ of invariant distribution $\pi$ and such that the distribution of $\theta_{i}$ converges to $\pi$ as $i\rightarrow\infty$. In practice one is typically interested in the computation of expectations of functions, say $f$, with respect to $\pi$ and it is also required that averages $M^{-1}\sum_{n=1}^{M}f(\theta_{n})$ converge to the expectation of interest. The iterative nature of MCMC makes it difficult to develop generic methods to take advantage of parallel computing environments when interested in reducing time to convergence. While numerous approaches have been proposed to reduce the variance of ergodic averages, including averaging over independent realisations of $\{\theta_{n},n\geq1\}$ simulated on several computers, techniques to reduce the "burn-in" of MCMC are scarce. In this paper we explore a simple and generic approach to improve convergence to equilibrium of existing algorithms which rely on the Metropolis-Hastings (MH) update, the main building block of MCMC. The main idea is to use averages of the acceptance ratio w.r.t. multiple realisations of random variables involved, while preserving $\pi$ as invariant distribution. The methodology requires limited change to existing code, is naturally suited to parallel computing and is shown on our examples to provide substantial performance improvements both in terms of convergence to equilibrium and variance of ergodic averages. In some scenarios gains are observed even on a serial machine.


title: Partial Identification in Nonseparable Binary Response Models with Endogenous Regressors

abstract: This paper considers (partial) identification of a variety of counterfactual parameters in binary response models with possibly endogenous regressors. Our framework allows for nonseparable index functions with multi-dimensional latent variables, and does not require parametric distributional assumptions. We leverage results on hyperplane arrangements and cell enumeration from the literature on computational geometry in order to provide a tractable means of computing the identified set. We demonstrate how various functional form, independence, and monotonicity assumptions can be imposed as constraints in our optimization procedure to tighten the identified set. Finally, we apply our method to study the effects of health insurance on the decision to seek medical treatment.


title: Quantitative Corner Case Feature Analysis of Hybrid Automata with ForFET$^{SMT}$

abstract: The analysis and verification of hybrid automata (HA) models against rich formal properties can be a challenging task. Existing methods and tools can mainly reason whether a given property is satisfied or violated. However, such qualitative answers might not provide sufficient information about the model behaviors. This paper presents the ForFET$^{SMT}$ tool which can be used to reason quantitatively about such properties. It employs feature automata and can evaluate quantitative property corners of HA. ForFET$^{SMT}$ uses two third-party formal verification tools as its backbone: the SpaceEx reachability tool and the SMT solver dReach/dReal. Herein, we describe the design and implementation of ForFET$^{SMT}$ and present its functionalities and modules. To improve the usability of the tool for non-expert users, we also provide a list of quantitative property templates.


title: Three lens space summands from the Poincar\'e homology sphere

abstract: A regular fiber of the Seifert fibering of the Poincar\'e homology sphere admits a Dehn surgery to $L(2,1)\#L(3,2)\#L(5,4)$. We prove that this is the only knot in the Poincar\'e homology sphere with a surgery to a connected sum of more than two lens spaces.


title: CFAR-Based Interference Mitigation for FMCW Automotive Radar Systems

abstract: In this paper, constant false alarm rate (CFAR) detector-based approaches are proposed for interference mitigation of Frequency modulated continuous wave (FMCW) radars. The proposed methods exploit the fact that after dechirping and low-pass filtering operations the targets' beat signals of FMCW radars are composed of exponential sinusoidal components while interferences exhibit short chirp waves within a sweep. The spectra of interferences in the time-frequency ($t$-$f$) domain are detected by employing a 1-D CFAR detector along each frequency bin and then the detected map is dilated as a mask for interference suppression. They are applicable to the scenarios in the presence of multiple interferences. Compared to the existing methods, the proposed methods reduce the power loss of useful signals and are very computationally efficient. Their interference mitigation performances are demonstrated through both numerical simulations and experimental results.


title: Spontaneous fractional Chern insulators in transition metal dichalcogenides Moire superlattices

abstract: Moir{\'e} superlattice realized in two-dimensional heterostructures offers an exciting platform to access strongly-correlated electronic states. In this work, we study transition metal dichalcogenides (TMD) Moir{\'e} superlattices with time-reversal-symmetry and nontrivial spin{/valley}-Chern numbers. Utilizing realistic material parameters and the method of exact diagonalization, we find that at a certain twisting angle and fractional filling, gapped fractional topological states, i.e., fractional Chern insulators, are naturally {stabilized} by simply introducing the Coulomb repulsion. In contrast to fractional quantum Hall systems, where the time-reversal symmetry has to be broken explicitly, these fractional states break the time-reversal symmetry spontaneously. {We show that the Chern number contrasting in the opposite valleys imposes a strong constraint on the nature of fractional Chern insulator and the associated low energy excitations.} We also propose to realize the non-abelian Moore-Read state in TMD Moir{\'e} superlattice sandwiched between nonlinear dielectric media.


title: Prior Knowledge Input to Improve LSTM Auto-encoder-based Characterization of Vehicular Sensing Data

abstract: Precision in event characterization in connected vehicles has become increasingly important with the responsive connectivity that is available to the modern vehicles. Event characterization via vehicular sensors are utilized in safety and autonomous driving applications in vehicles. While characterization systems have been shown to be capable of predicting the risky driving patterns, precision of such systems still remains an open issue. The major issues against the driving event characterization systems need to be addressed in connected vehicle settings, which are the heavy imbalance and the event infrequency of the driving data and the existence of the time-series detection systems that are optimized for vehicular settings. To overcome the problems, we introduce the application of the prior-knowledge input method to the characterization systems. Furthermore, we propose a recurrent-based denoising auto-encoder network to populate the existing data for a more robust training process. The results of the conducted experiments show that the introduction of knowledge-based modelling enables the existing systems to reach significantly higher accuracy and F1-score levels. Ultimately, the combination of the two methods enables the proposed model to attain 14.7\% accuracy boost over the baseline by achieving an accuracy of 0.96.


title: SpotPatch: Parameter-Efficient Transfer Learning for Mobile Object Detection

abstract: Deep learning based object detectors are commonly deployed on mobile devices to solve a variety of tasks. For maximum accuracy, each detector is usually trained to solve one single specific task, and comes with a completely independent set of parameters. While this guarantees high performance, it is also highly inefficient, as each model has to be separately downloaded and stored. In this paper we address the question: can task-specific detectors be trained and represented as a shared set of weights, plus a very small set of additional weights for each task? The main contributions of this paper are the following: 1) we perform the first systematic study of parameter-efficient transfer learning techniques for object detection problems; 2) we propose a technique to learn a model patch with a size that is dependent on the difficulty of the task to be learned, and validate our approach on 10 different object detection tasks. Our approach achieves similar accuracy as previously proposed approaches, while being significantly more compact.


title: Hedging with Bitcoin Futures: The Effect of Liquidation Loss Aversion and Aggressive Trading

abstract: We consider the hedging problem where a futures position can be automatically liquidated by the exchange without notice. We derive a semi-closed form for an optimal hedging strategy with dual objectives - to minimise both the variance of the hedged portfolio and the probability of liquidations due to insufficient collateral. The optimal solution depends on the statistical characteristics of the spot and futures extreme returns and parameters that characterise the hedger by loss aversion, choice of leverage and collateral management. An empirical analysis of bitcoin shows that the optimal strategy combines superior hedge effectiveness with a reduction in the probability of liquidation. We compare the performance of seven major direct and inverse hedging instruments traded on five different exchanges, based on minute-level data. We also link this performance to novel speculative trading metrics, which differ markedly between venues.


title: Voltage-controlled long-range propagation of indirect excitons in van der Waals heterostructure

abstract: Indirect excitons (IXs), also known as interlayer excitons, can form the medium for excitonic devices whose operation is based on controlled propagation of excitons. A proof of principle for excitonic devices was demonstrated in GaAs heterostructures where the operation of excitonic devices is limited to low temperatures. IXs in van der Waals transition-metal dichalcogenide (TMD) heterostructures are characterized by high binding energies making IXs robust at room temperature and offering an opportunity to create excitonic devices operating at high temperatures suitable for applications. However, a characteristic feature of TMD heterostructures is the presence of moir\'e superlattice potentials, which are predicted to cause modulations of IX energy reaching tens of meV. These in-plane energy landscapes can lead to IX localization, making IX propagation fundamentally different in TMD and GaAs heterostructures and making uncertain if long-range IX propagation can be realized in TMD heterostructures. In this work, we realize long-range IX propagation with the $1/e$ IX luminescence decay distances reaching 13 microns in a MoSe$_2$/WSe$_2$ heterostructure. We trace the IX luminescence along the IX propagation path. We also realize control of the long-range IX propagation: the IX luminescence signal in the drain of an excitonic transistor is controlled within 40 times by gate voltage. These data show that the long-range IX propagation is possible in TMD heterostructures with the predicted moir\'e superlattice potentials.


title: Finding the Sequence of Largest Small n-Polygons by Numerical Optimization

abstract: LSP(n), the largest small polygon with n vertices, is the polygon of unit diameter that has maximal area A(n). It is known that for all odd values $n \geq 3$, LSP(n) is the regular n-polygon; however, this statement is not valid for even values of n. Finding the polygon LSP(n) and A(n) for even values $n \geq 6$ has been a long-standing challenge. In this work, we develop high-precision numerical solution estimates of A(n) for even values $n \geq 4$, using the Mathematica model development environment and the IPOPT local nonlinear optimization solver engine. First, we present a revised (tightened) LSP model that greatly assists the efficient solution of the model-class considered. This is followed by numerical results for an illustrative sequence of even values of n, up to $n \leq 1000$. Our results are in close agreement with, or surpass, the best results reported in all earlier studies. Most of these earlier works addressed special cases up to $n \leq 20$, while others obtained numerical optimization results for a range of values from $6 \leq n \leq 100$. For completeness, we also calculate numerically optimized results for a selection of odd values of n, up to $n \leq 999$: these results can be compared to the corresponding theoretical (exact) values. The results obtained are used to provide regression model-based estimates of the optimal area sequence {A(n)}, for all even and odd values n of interest, thereby essentially solving the entire LSP model-class numerically, with demonstrably high precision.


title: A Research Ecosystem for Secure Computing

abstract: Computing devices are vital to all areas of modern life and permeate every aspect of our society. The ubiquity of computing and our reliance on it has been accelerated and amplified by the COVID-19 pandemic. From education to work environments to healthcare to defense to entertainment - it is hard to imagine a segment of modern life that is not touched by computing. The security of computers, systems, and applications has been an active area of research in computer science for decades. However, with the confluence of both the scale of interconnected systems and increased adoption of artificial intelligence, there are many research challenges the community must face so that our society can continue to benefit and risks are minimized, not multiplied. Those challenges range from security and trust of the information ecosystem to adversarial artificial intelligence and machine learning. Along with basic research challenges, more often than not, securing a system happens after the design or even deployment, meaning the security community is routinely playing catch-up and attempting to patch vulnerabilities that could be exploited any minute. While security measures such as encryption and authentication have been widely adopted, questions of security tend to be secondary to application capability. There needs to be a sea-change in the way we approach this critically important aspect of the problem: new incentives and education are at the core of this change. Now is the time to refocus research community efforts on developing interconnected technologies with security "baked in by design" and creating an ecosystem that ensures adoption of promising research developments. To realize this vision, two additional elements of the ecosystem are necessary - proper incentive structures for adoption and an educated citizenry that is well versed in vulnerabilities and risks.


title: An identity concerning the Riemann-zeta function

abstract: For a certain function $J(s)$ we prove that the identity $$\frac{\zeta(2s)}{\zeta(s)}-\left(s-\frac{1}{2}\right)J(s)=\frac{\zeta(2s+1)}{\zeta(s+1/2)}, $$ holds in the half-plane Re$(s)>1/2$ and both sides of the equality are analytic in this half-plane.


title: Federated Learning-Based Risk-Aware Decision toMitigate Fake Task Impacts on CrowdsensingPlatforms

abstract: Mobile crowdsensing (MCS) leverages distributed and non-dedicated sensing concepts by utilizing sensors imbedded in a large number of mobile smart devices. However, the openness and distributed nature of MCS leads to various vulnerabilities and consequent challenges to address. A malicious user submitting fake sensing tasks to an MCS platform may be attempting to consume resources from any number of participants' devices; as well as attempting to clog the MCS server. In this paper, a novel approach that is based on horizontal federated learning is proposed to identify fake tasks that contain a number of independent detection devices and an aggregation entity. Detection devices are deployed to operate in parallel with each device equipped with a machine learning (ML) module, and an associated training dataset. Furthermore, the aggregation module collects the prediction results from individual devices and determines the final decision with the objective of minimizing the prediction loss. Loss measurement considers the lost task values with respect to misclassification, where the final decision utilizes a risk-aware approach where the risk is formulated as a function of the utility loss. Experimental results demonstrate that using federated learning-driven illegitimate task detection with a risk aware aggregation function improves the detection performance of the traditional centralized framework. Furthermore, the higher performance of detection and lower loss of utility can be achieved by the proposed framework. This scheme can even achieve 100%detection accuracy using small training datasets distributed across devices, while achieving slightly over an 8% increase in detection improvement over traditional approaches.


title: Heterologous autoimmunity and prokaryotic immune defense

abstract: Some prokaryotes possess CRISPR-Cas systems that provide adaptive immunity to viruses guided by DNA segments called spacers acquired from invading phage. However, the patchy incidence and limited memory breadth of CRISPR-Cas systems suggest that their fitness benefits are offset by costs. Here, we propose that cross-reactive CRISPR targeting can lead to heterologous autoimmunity, whereby foreign spacers guide self-targeting in a spacer-length dependent fashion. Balancing antiviral defense against autoimmunity predicts a scaling relation between spacer length and CRISPR repertoire size. We find evidence for this scaling through comparative analysis of sequenced prokaryotic genomes, and show that this association also holds at the level of CRISPR types. In contrast, the scaling is absent in strains with nonfunctional CRISPR loci. Finally, we demonstrate that stochastic spacer loss can explain variations around the scaling relation, even between strains of the same species. Our results suggest that heterologous autoimmunity is a selective factor shaping the evolution of CRISPR-Cas systems.


title: PSF Estimation in Crowded Astronomical Imagery as a Convolutional Dictionary Learning Problem

abstract: We present a new algorithm for estimating the Point Spread Function (PSF) in wide-field astronomical images with extreme source crowding. Robust and accurate PSF estimation in crowded astronomical images dramatically improves the fidelity of astrometric and photometric measurements extracted from wide-field sky monitoring imagery. Our radically new approach utilizes convolutional sparse representations to model the continuous functions involved in the image formation. This approach avoids the need to detect and precisely localize individual point sources that is shared by existing methods. In experiments involving simulated astronomical imagery, it significantly outperforms the recent alternative method with which it is compared.