title: Antiviral Drug-Membrane Permeability: the Viral Envelope and Cellular Organelles

abstract: To shorten the time required to find effective new drugs, like antivirals, a key parameter to consider is membrane permeability, as a compound intended for an intracellular target with poor permeability will have low efficacy. Here, we present a computational model that considers both drug characteristics and membrane properties for the rapid assessment of drugs permeability through the coronavirus envelope and various cellular membranes. We analyze 79 drugs that are considered as potential candidates for the treatment of SARS-CoV-2 and determine their time of permeation in different organelle membranes grouped by viral baits and mammalian processes. The computational results are correlated with experimental data, present in the literature, on bioavailability of the drugs, showing a negative correlation between fast permeation and most promising drugs. This model represents an important tool capable of evaluating how permeability affects the ability of compounds to reach both intended and unintended intracellular targets in an accurate and rapid way. The method is general and flexible and can be employed for a variety of molecules, from small drugs to nanoparticles, as well to a variety of biological membranes.


title: Mirostat: A Neural Text Decoding Algorithm that Directly Controls Perplexity

abstract: Neural text decoding is important for generating high-quality texts using language models. To generate high-quality text, popular decoding algorithms like top-k, top-p (nucleus), and temperature-based sampling truncate or distort the unreliable low probability tail of the language model. Though these methods generate high-quality text after parameter tuning, they are ad hoc. Not much is known about the control they provide over the statistics of the output, which is important since recent reports show text quality is highest for a specific range of likelihoods. Here, first we provide a theoretical analysis of perplexity in top-k, top-p, and temperature sampling, finding that cross-entropy behaves approximately linearly as a function of p in top-p sampling whereas it is a nonlinear function of k in top-k sampling, under Zipfian statistics. We use this analysis to design a feedback-based adaptive top-k text decoding algorithm called mirostat that generates text (of any length) with a predetermined value of perplexity, and thereby high-quality text without any tuning. Experiments show that for low values of k and p in top-k and top-p sampling, perplexity drops significantly with generated text length, which is also correlated with excessive repetitions in the text (the boredom trap). On the other hand, for large values of k and p, we find that perplexity increases with generated text length, which is correlated with incoherence in the text (confusion trap). Mirostat avoids both traps: experiments show that cross-entropy has a near-linear relation with repetition in generated text. This relation is almost independent of the sampling method but slightly dependent on the model used. Hence, for a given language model, control over perplexity also gives control over repetitions. Experiments with human raters for fluency, coherence, and quality further verify our findings.


title: Defining Pointwise Lower Scalar Curvature Bounds for $C^0$ Metrics with Regularization by Ricci Flow

abstract: We survey some recent work using Ricci flow to create a class of local definitions of weak lower scalar curvature bounds that is well defined for $C^0$ metrics. We discuss several properties of these definitions and explain some applications of this approach to questions regarding uniform convergence of metrics with scalar curvature bounded below. Finally, we consider the relationship between this approach and some other generalized notions of lower scalar curvature bounds.


title: Theory of Crystallization versus Vitrification

abstract: The competition between crystallization and vitrification in glass-forming materials manifests as a non-monotonic behavior in the time-temperature transformation (TTT) diagrams, which quantify the time scales for crystallization as a function of temperature. We develop a coarse-grained lattice model, the Arrow-Potts model, to explore the physics behind this competition. Using Monte Carlo simulations, the model showcases non-monotonic TTT diagrams resulting in polycrystalline structures, with two distinct regimes limited by either crystal nucleation or growth. At high temperatures, crystallization is limited by nucleation and results in the growth of compact crystal grains. At low temperatures, crystal growth is influenced by glassy dynamics, and proceeds through dynamically heterogeneous and hierarchical relaxation pathways producing fractal and ramified crystals. To explain these phenomena, we combine the Kolmogorov-Johnson-Mehl-Avrami theory with the field theory of nucleation, a random walk theory for crystal growth, and the dynamical facilitation theory for glassy dynamics. The unified theory yields an analytical formula relating crystallization timescale to the nucleation and growth rates through universal exponents governing glassy dynamics of the model. We show that the formula with the universal exponents yields excellent agreement with the Monte Carlo simulation data and thus, it also accounts for the non-monotonic TTT diagrams produced by the model. Both the model and theory can be used to understand structural ordering in various glassy systems including bulk metallic glass alloys, organic molecules, and colloidal suspensions.


title: A survey for variable young stars with small telescopes: III -- Warm spots on the active star V1598Cyg

abstract: Magnetic spots on low-mass stars can be traced and characterised using multi-band photometric light curves. Here we analyse an extensive data set for one active star, V1598Cyg, a known variable K dwarf which is either pre-main sequence and/or in a close binary system. Our light curve contains 2854 photometric data points, mostly in $V$, $R_c$, $I_c$, but also in $U$, $B$ and $H\alpha$, with a total baseline of about 4yr, obtained with small telescopes as part of the HOYS project. We find that V1598Cyg is a very fast rotator with a period of 0.8246 days and varying amplitudes in all filters, best explained as a signature of strong magnetic activity and spots. We fit the photometric amplitudes in $V$, $R_c$, $I_c$ and use them to estimate spot properties, using a grid-based method that is also propagating uncertainties. We verify the method on a partial data set with high cadence and all five broad-band filters. The method yields spot temperatures and fractional spot coverage with typical uncertainties of 100K and 3-4%, respectively. V1598Cyg consistently exhibits spots that are a few hundred degrees warmer than the photosphere, most likely indicating that the light curve is dominated by chromospheric plage. The spot activity varies over our observing baseline, with a typical time scale of 0.5-1yr, which we interpret as the typical spot lifetime. Combining our light curve with archival data, we find a six year cycle in the average brightness, that is probably a sign of a magnetic activity cycle.


title: maxsmooth: Rapid Maximally Smooth Function Fitting With Applications in Global 21-cm Cosmology

abstract: Maximally Smooth Functions (MSFs) are a form of constrained functions in which there are no inflection points or zero crossings in high order derivatives. Consequently, they have applications to signal recovery in experiments where signals of interest are expected to be non-smooth features masked by larger smooth signals or foregrounds. They can also act as a powerful tool for diagnosing the presence of systematics. The constrained nature of MSFs makes fitting these functions a non-trivial task. We introduce maxsmooth, an open source package that uses quadratic programming to rapidly fit MSFs. We demonstrate the efficiency and reliability of maxsmooth by comparison to commonly used fitting routines and show that we can reduce the fitting time by approximately two orders of magnitude. We introduce and implement with maxsmooth Partially Smooth Functions, which are useful for describing elements of non-smooth structure in foregrounds. This work has been motivated by the problem of foreground modelling in 21-cm cosmology. We discuss applications of maxsmooth to 21-cm cosmology and highlight this with examples using data from the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) and the Large-aperture Experiment to Detect the Dark Ages (LEDA) experiments. We demonstrate the presence of a sinusoidal systematic in the EDGES data with a log-evidence difference of $86.19\pm0.12$ when compared to a pure foreground fit. MSFs are applied to data from LEDA for the first time in this paper and we identify the presence of sinusoidal systematics. maxsmooth is pip installable and available for download at: https://github.com/htjb/maxsmooth


title: Equivalence theorems for compound design problems with application in mixed models

abstract: In the present paper we consider design criteria which depend on several designs simultaneously. We formulate equivalence theorems based on moment matrices (if criteria depend on designs via moment matrices) or with respect to the designs themselves (for finite design regions). We apply the obtained optimality conditions to the multiple-group random coefficient regression models and illustrate the results by simple examples.


title: Families of Gr\"obner Degenerations, Grassmannians and Universal Cluster Algebras

abstract: Let $V$ be the weighted projective variety defined by a weighted homogeneous ideal $J$ and $C$ a maximal cone in the Gr\"obner fan of $J$ with $m$ rays. We construct a flat family over $\mathbb A^m$ that assembles the Gr\"obner degenerations of $V$ associated with all faces of $C$. This is a multi-parameter generalization of the classical one-parameter Gr\"obner degeneration associated to a weight. We explain how our family can be constructed from Kaveh-Manon's recent work on the classification of toric flat families over toric varieties: it is the pull-back of a toric family defined by a Rees algebra with base $X_C$ (the toric variety associated to $C$) along the universal torsor $\mathbb A^m \to X_C$. We apply this construction to the Grassmannians ${\rm Gr}(2,\mathbb C^n)$ with their Pl\"ucker embeddings and the Grassmannian ${\rm Gr}\big(3,\mathbb C^6\big)$ with its cluster embedding. In each case, there exists a unique maximal Gr\"obner cone whose associated initial ideal is the Stanley-Reisner ideal of the cluster complex. We show that the corresponding cluster algebra with universal coefficients arises as the algebra defining the flat family associated to this cone. Further, for ${\rm Gr}(2,\mathbb C^n)$ we show how Escobar-Harada's mutation of Newton-Okounkov bodies can be recovered as tropicalized cluster mutation.


title: Transient Jupiter Co-orbitals from Solar System Sources

abstract: We demonstrate dynamical pathways from main-belt asteroid and Centaur orbits to those in co-orbital motion with Jupiter, including the retrograde (inclination $i>90^o$) state. We estimate that at any given time, there should be $\sim1$ kilometer-scale or larger escaped asteroid in a transient direct (prograde) orbit with semimajor axis near that of Jupiter's ($a\simeq a_J$), with proportionally more smaller objects as determined by their size distribution. Most of these objects would be in the horseshoe dynamical state, which are hard to detect due to their moderate eccentricities (spending most of their time beyond 5 AU) and longitudes relative to Jupiter being spread nearly all over the sky. We also show that $\approx$1% of the transient asteroid co-orbital population is on retrograde orbits with Jupiter. This population, like the recently identified asteroid (514107) 2015 BZ$_{509}$, can spend millions of years with $a\simeq a_J$ including tens or hundreds of thousands of years formally in the retrograde 1:-1 co-orbital resonance. Escaping near-Earth asteroids (NEAs) are thus likely the precursors to the handful of known high-inclination objects with $a\simeq a_J$. We compare the production of jovian co-orbitals from escaping NEAs with those from incoming Centaurs. We find that temporary direct co-orbitals are likely dominated by Centaur capture, but we only find production of (temporary) retrograde jovian co-orbitals (including very long-lived ones) from the NEA source. We postulate that the primordial elimination of the inner Solar System's planetesimal population could provide a supply route for a metastable outer Solar System reservoir for the high-inclination Centaurs.


title: On Loss Functions and Recurrency Training for GAN-based Speech Enhancement Systems

abstract: Recent work has shown that it is feasible to use generative adversarial networks (GANs) for speech enhancement, however, these approaches have not been compared to state-of-the-art (SOTA) non GAN-based approaches. Additionally, many loss functions have been proposed for GAN-based approaches, but they have not been adequately compared. In this study, we propose novel convolutional recurrent GAN (CRGAN) architectures for speech enhancement. Multiple loss functions are adopted to enable direct comparisons to other GAN-based systems. The benefits of including recurrent layers are also explored. Our results show that the proposed CRGAN model outperforms the SOTA GAN-based models using the same loss functions and it outperforms other non-GAN based systems, indicating the benefits of using a GAN for speech enhancement. Overall, the CRGAN model that combines an objective metric loss function with the mean squared error (MSE) provides the best performance over comparison approaches across many evaluation metrics.


title: Objective frequentist uncertainty quantification for atmospheric CO$_2$ retrievals

abstract: The steadily increasing amount of atmospheric carbon dioxide (CO$_2$) is affecting the global climate system and threatening the long-term sustainability of Earth's ecosystem. In order to better understand the sources and sinks of CO$_2$, NASA operates the Orbiting Carbon Observatory-2 & 3 satellites to monitor CO$_2$ from space. These satellites make passive radiance measurements of the sunlight reflected off the Earth's surface in different spectral bands, which are then inverted in an ill-posed inverse problem to obtain estimates of the atmospheric CO$_2$ concentration. In this work, we propose a new CO$_2$ retrieval method that uses known physical constraints on the state variables and direct inversion of the target functional of interest to construct well-calibrated frequentist confidence intervals based on convex programming. We compare the method with the current operational retrieval procedure, which uses prior knowledge in the form of probability distributions to regularize the problem. We demonstrate that the proposed intervals consistently achieve the desired frequentist coverage, while the operational uncertainties are poorly calibrated in a frequentist sense both at individual locations and over a spatial region in a realistic simulation experiment. We also study the influence of specific nuisance state variables on the length of the proposed intervals and identify certain key variables that can greatly reduce the final uncertainty given additional deterministic or probabilistic constraints, and develop a principled framework to incorporate such information into our method.


title: The spontaneous exchange bias effect on La$_{1.5}$Ba$_{0.5}$CoMnO$_{6}$

abstract: La$_{1.5}$Ba$_{0.5}$CoMnO$_{6}$ is a re-entrant cluster glass material exhibiting a robust negative exchange bias (EB) effect even after being cooled from an unmagnetized state down to low temperature in zero magnetic field. Here we thoroughly investigate this phenomena by performing magnetization as a function of applied field [$M(H)$] measurements at several different temperatures and maximum applied magnetic fields ($H_m$). The spontaneous EB (SEB) effect is observed below 20 K, and shows a maximum value for $H_m$ = 75 kOe. The effect is greatly enhanced when the $M(H)$ curves are measured after the system is cooled in the presence of a magnetic field. The asymmetry of the $M(H)$ curves here investigated can be well described by a recently proposed model based on the unconventional relaxation of the SG-like moments during the hysteresis cycle.


title: Powerful ionized gas outflows in the interacting radio galaxy 4C +29.30

abstract: We investigate the ionized gas excitation and kinematics in the inner $4.3 \times 6.2$ kpc$^{2}$ of the merger radio galaxy 4C +29.30. Using optical integral field spectroscopy with the Gemini North Telescope, we present flux distributions, line-ratio maps, peak velocities and velocity dispersion maps as well as channel maps with a spatial resolution of $\approx 955$ pc. We observe high blueshifts of up to $\sim -650$ km s$^{-1}$, in a region $\sim 1''$ south of the nucleus (the southern knot, SK), which also presents high velocity dispersions ($\sim 250$ km s$^{-1}$), which we attribute to an outflow. A possible redshifted counterpart is observed north from the nucleus (the northern knot, NK). We propose that these regions correspond to a bipolar outflow possibly due to the interaction of the radio jet with the ambient gas. We estimate a total ionized gas mass outflow rate of $\dot{M}_{out} = 25.4 \substack{+11.5 \\ -7.5}$ M$_\odot$ yr$^{-1}$ with a kinetic power of $\dot{E} = 8.1 \substack{+10.7 \\ -4.0} \times 10^{42}$ erg s$^{-1}$, which represents $5.8 \substack{+7.6 \\ -2.9} \%$ of the AGN bolometric luminosity. These values are higher than usually observed in nearby active galaxies with the same bolometric luminosities and could imply a significant impact of the outflows in the evolution of the host galaxy. The excitation is higher in the NK (that correlates with extended X-ray emission, indicating the presence of hotter gas) than in the SK, supporting a scenario in which an obscuring dust lane is blocking part of the AGN radiation to reach the southern region of the galaxy.


title: Higgs inflation in Einstein-Cartan gravity

abstract: We study inflation driven by the Higgs field in the Einstein-Cartan formulation of gravity. In this theory, the presence of the Holst and Nieh-Yan terms with the Higgs field non-minimally coupled to them leads to three additional coupling constants. For a broad range of parameters, we find that inflation is both possible and consistent with observations. In most cases, the spectral index is given by $n_s=1-2/N_\star$ (with $N_\star$ the number of e-foldings) whereas the tensor-to-scalar ratio $r$ can vary between about $10^{-10}$ and $1$. Thus, there are scenarios of Higgs inflation in the Einstein-Cartan framework for which the detection of gravitational waves from inflation is possible in the near future. In certain limits, the known models of Higgs inflation in the metric and Palatini formulations of gravity are reproduced. Finally, we discuss the robustness of inflationary dynamics against quantum corrections due to the scalar and fermion fields.


title: Neural Network-based Reconstruction in Compressed Sensing MRI Without Fully-sampled Training Data

abstract: Compressed Sensing MRI (CS-MRI) has shown promise in reconstructing under-sampled MR images, offering the potential to reduce scan times. Classical techniques minimize a regularized least-squares cost function using an expensive iterative optimization procedure. Recently, deep learning models have been developed that model the iterative nature of classical techniques by unrolling iterations in a neural network. While exhibiting superior performance, these methods require large quantities of ground-truth images and have shown to be non-robust to unseen data. In this paper, we explore a novel strategy to train an unrolled reconstruction network in an unsupervised fashion by adopting a loss function widely-used in classical optimization schemes. We demonstrate that this strategy achieves lower loss and is computationally cheap compared to classical optimization solvers while also exhibiting superior robustness compared to supervised models. Code is available at https://github.com/alanqrwang/HQSNet.


title: Moments of the doubly truncated selection elliptical distributions with emphasis on the unified multivariate skew-$t$ distribution

abstract: In this paper, we compute doubly truncated moments for the selection elliptical (SE) class of distributions, which includes some multivariate asymmetric versions of well-known elliptical distributions, such as, the normal, Student's t, slash, among others. We address the moments for doubly truncated members of this family, establishing neat formulation for high order moments as well as for its first two moments. We establish sufficient and necessary conditions for the existence of these truncated moments. Further, we propose optimized methods able to deal with extreme setting of the parameters, partitions with almost zero volume or no truncation which are validated with a brief numerical study. Finally, we present some results useful in interval censoring models. All results has been particularized to the unified skew-t (SUT) distribution, a complex multivariate asymmetric heavy-tailed distribution which includes the extended skew-t (EST), extended skew-normal (ESN), skew-t (ST) and skew-normal (SN) distributions as particular and limiting cases.


title: Incommensurate magnetic states induced by ordering competition in $\mathrm{Ba_{1-x}Na_xFe_2As_2}$

abstract: Quantum criticality nearby a certain magnetic phase transition beneath the superconducting dome of $\mathrm{Ba_{1-x}Na_xFe_2As_2}$ is attentively studied by virtue of a phenomenological theory in conjunction with renormalization group approach. We report that ordering competition between magnetic and superconducting fluctuations is capable of coaxing incommensurate (IC) magnetic states to experience distinct fates depending upon their spin configurations. The $C_2$-symmetry IC magnetic stripe with perpendicular magnetic helix dominates over other $C_2$-symmetry magnetic competitors and hints at a potential candidate for the unknown $C_2$-symmetry magnetic state. Amongst $C_4$-symmetry IC magnetic phases, IC charge spin density wave is substantiated to be superior, shedding light on the significant intertwining of charge and spin degrees of freedom. Meanwhile, ferocious fluctuations render a sharp fall of superfluid density alongside dip of critical temperature as well as intriguing behavior of London penetration depth.


title: Homogeneous Floquet time crystal from weak ergodicity breaking

abstract: Recent works on observation of discrete time-crystalline signatures throw up major puzzles on the necessity of localization for stabilizing such out-of-equilibrium phases. Motivated by these studies, we delve into a clean interacting Floquet system, whose quasi-spectrum conforms to the ergodic Wigner-Dyson distribution, yet with an unexpectedly robust, long-lived time-crystalline dynamics in the absence of disorder or fine-tuning. We relate such behavior to a measure zero set of nonthermal Floquet eigenstates with long-range spatial correlations, which coexist with otherwise thermal states at near-infinite temperature and develop a high overlap with a family of translationally invariant, symmetry-broken initial conditions. This resembles the notion of "dynamical scars" that remain robustly localized throughout a thermalizing Floquet spectrum with fractured structure. We dub such a long-lived discrete time crystal formed in partially nonergodic systems, "scarred discrete time crystal" which is distinct by nature from those stabilized by either many-body localization or prethermalization mechanism.


title: Adversarial Robustness for Machine Learning Cyber Defenses Using Log Data

abstract: There has been considerable and growing interest in applying machine learning for cyber defenses. One promising approach has been to apply natural language processing techniques to analyze logs data for suspicious behavior. A natural question arises to how robust these systems are to adversarial attacks. Defense against sophisticated attack is of particular concern for cyber defenses. In this paper, we develop a testing framework to evaluate adversarial robustness of machine learning cyber defenses, particularly those focused on log data. Our framework uses techniques from deep reinforcement learning and adversarial natural language processing. We validate our framework using a publicly available dataset and demonstrate that our adversarial attack does succeed against the target systems, revealing a potential vulnerability. We apply our framework to analyze the influence of different levels of dropout regularization and find that higher dropout levels increases robustness. Moreover 90% dropout probability exhibited the highest level of robustness by a significant margin, which suggests unusually high dropout may be necessary to properly protect against adversarial attacks.


title: Cosmological Model Parameter Dependence of the Matter Power Spectrum Covariance from the DEUS-PUR $Cosmo$ Simulations

abstract: Future galaxy surveys will provide accurate measurements of the matter power spectrum across an unprecedented range of scales and redshifts. The analysis of these data will require one to accurately model the imprint of non-linearities of the matter density field. In particular, these induce a non-Gaussian contribution to the data covariance that needs to be properly taken into account to realise unbiased cosmological parameter inference analyses. Here, we study the cosmological dependence of the matter power spectrum covariance using a dedicated suite of N-body simulations, the Dark Energy Universe Simulation - Parallel Universe Runs (DEUS-PUR) {\it Cosmo}. These consist of 512 realizations for 10 different cosmologies where we vary the matter density $\Omega_m$, the amplitude of density fluctuations $\sigma_8$, the reduced Hubble parameter $h$ and a constant dark energy equation of state $w$ by approximately $10\%$. We use these data to evaluate the first and second derivatives of the power spectrum covariance with respect to a fiducial $\Lambda$CDM cosmology. We find that the variations can be as large as $150\%$ depending on the scale, redshift and model parameter considered. By performing a Fisher matrix analysis we explore the impact of different choices in modelling the cosmological dependence of the covariance. Our results suggest that fixing the covariance to a fiducial cosmology can significantly affect the recovered parameter errors and that modelling the cosmological dependence of the variance while keeping the correlation coefficient fixed can alleviate the impact of this effect.


title: Comparison of recoil polarization in the $^{12}{\rm C}(\vec{e},{e}'\vec{p})$ process for protons extracted from $s$ and $p$ shell

abstract: We present first measurements of the double ratio of the polarization transfer components $(P^{\prime}_{\!x} \!/ P^{\prime}_{\!z} )_p/ (P^{\prime}_{\!x} \!/ P^{\prime}_{\!z} )_s$ for knock-out protons from $s$ and $p$ shells in $^{12}{\rm C}$ measured by the $^{12}{\rm C}(\vec{e},{e}'\vec{p}\,)$ reaction in quasi-elastic kinematics. The data are compared to theoretical predictions in relativistic distorted-wave impulse approximation. Our results show that differences between $s$- and $p$-shell protons, observed when compared at the same initial momentum (missing momentum) largely disappear when the comparison is done at the same proton virtuality. We observe no density-dependent medium modifications for protons from $s$ and $p$ shells with the same virtuality in spite of the large differences in the respective nuclear densities.


title: Investigation of Phase Distortion on Perceived Speech Quality for Hearing-impaired Listeners

abstract: Phase serves as a critical component of speech that influences the quality and intelligibility. Current speech enhancement algorithms are beginning to address phase distortions, but the algorithms focus on normal-hearing (NH) listeners. It is not clear whether phase enhancement is beneficial for hearing-impaired (HI) listeners. We investigated the influence of phase distortion on speech quality through a listening study, in which NH and HI listeners provided speech-quality ratings using the MUSHRA procedure. In one set of conditions, the speech was mixed with babble noise at 4 different signal-to-noise ratios (SNRs) from -5 to 10 dB. In another set of conditions, the SNR was fixed at 10 dB and the noisy speech was presented in a simulated reverberant room with T60s ranging from 100 to 1000 ms. The speech level was kept at 65 dB SPL for NH listeners and amplification was applied for HI listeners to ensure audibility. Ideal ratio masking (IRM) was used to simulate speech enhancement. Two objective metrics (i.e., PESQ and HASQI) were utilized to compare subjective and objective ratings. Results indicate that phase distortion has a negative impact on perceived quality for both groups and PESQ is more closely correlated with human ratings.


title: Presentation and Analysis of a Multimodal Dataset for Grounded Language Learning

abstract: Grounded language acquisition -- learning how language-based interactions refer to the world around them -- is amajor area of research in robotics, NLP, and HCI. In practice the data used for learning consists almost entirely of textual descriptions, which tend to be cleaner, clearer, and more grammatical than actual human interactions. In this work, we present the Grounded Language Dataset (GoLD), a multimodal dataset of common household objects described by people using either spoken or written language. We analyze the differences and present an experiment showing how the different modalities affect language learning from human in-put. This will enable researchers studying the intersection of robotics, NLP, and HCI to better investigate how the multiple modalities of image, text, and speech interact, as well as show differences in the vernacular of these modalities impact results.


title: The EFT Likelihood for Large-Scale Structure in Redshift Space

abstract: We study the EFT likelihood for biased tracers in redshift space, for which the bias expansion of the galaxy velocity field $\mathbf{v}_g$ plays a fundamental role. The equivalence principle forbids stochastic contributions to $\mathbf{v}_g$ to survive at small $k$. Therefore, at leading order in derivatives the form of the likelihood ${\cal P}[\tilde{\delta}_g|\delta,\!\mathbf{v}]$ to observe a redshift-space galaxy overdensity $\tilde{\delta}_g(\tilde{\mathbf{x}})$ given a rest-frame matter and velocity fields $\delta(\mathbf{x})$, $\mathbf{v}(\mathbf{x})$ is fixed by the rest-frame noise. If this noise is Gaussian with constant power spectrum, ${\cal P}[\tilde{\delta}_g|\delta,\!\mathbf{v}]$ is also a Gaussian in the difference between $\tilde{\delta}_g(\tilde{\mathbf{x}})$ and its bias expansion: redshift-space distortions only make the covariance depend on $\delta(\mathbf{x})$ and $\mathbf{v}(\mathbf{x})$. We then show how to match this result to perturbation theory, and that one can consistently neglect the field-dependent covariance if the bias expansion is stopped at second order in perturbations. We discuss qualitatively how this affects numerical implementations of the EFT-based forward modeling, and how the picture changes when the survey window function is taken into account.


title: Analytic marginalization of $N(z)$ uncertainties in tomographic galaxy surveys

abstract: We present a new method to marginalize over uncertainties in redshift distributions, $N(z)$, within tomographic cosmological analyses applicable to current and upcoming photometric galaxy surveys. We allow for arbitrary deviations from the best-guess $N(z)$ governed by a general covariance matrix describing the uncertainty in our knowledge of redshift distributions. In principle, this is marginalization over hundreds or thousands of new parameters describing potential deviations as a function of redshift and tomographic bin. However, by linearly expanding the theory predictions around a fiducial model, this marginalization can be performed analytically, resulting in a modified data covariance matrix that effectively downweights the modes of the data vector that are more sensitive to redshift distribution variations. We showcase this method by applying it to the galaxy clustering measurements from the Hyper Suprime-Cam first data release. We illustrate how to marginalize over sample-variance of the calibration sample and a large general systematic uncertainty in photometric estimation methods, and explore the impact of priors imposing smoothness in the redshift distributions.


title: Reliable Broadcast in Practical Networks: Algorithm and Evaluation

abstract: Reliable broadcast is an important primitive to ensure that a source node can reliably disseminate a message to all the non-faulty nodes in an asynchronous and failure-prone networked system. Byzantine Reliable Broadcast protocols were first proposed by Bracha in 1987, and have been widely used in fault-tolerant systems and protocols. Several recent protocols have improved the round and bit complexity of these algorithms. Motivated by the constraints in practical networks, we revisit the problem. In particular, we use cryptographic hash functions and erasure coding to reduce communication and computation complexity and simplify the protocol design. We also identify the fundamental trade-offs of Byzantine Reliable Broadcast protocols with respect to resilience (number of nodes), local computation, round complexity, and bit complexity. Finally, we also design and implement a general testing framework for similar communication protocols. We evaluate our protocols using our framework. The results demonstrate that our protocols have superior performance in practical networks.


title: Analytical construction of the projectile motion trajectory in midair

abstract: A classic problem of the motion of a projectile thrown at an angle to the horizon is studied. Air resistance force is taken into account. The quadratic law for the resistance force is used. An analytic approach applies for the investigation. Equations of the projectile motion are solved analytically for an arbitrarily large period of time. The constructed analytical solutions are universal. As a limit case of motion, the vertical asymptote formula is obtained. The value of the vertical asymptote is calculated directly from the initial conditions of motion. There is no need to study the problem numerically. The found analytical solutions are highly accurate over a wide range of parameters. The motion of a baseball, a tennis ball and shuttlecock of badminton are presented as examples.


title: A First-Principle Model for Polarization Swings during Reconnection-Powered Flares

abstract: We show that magnetic reconnection in a magnetically-dominated fast-cooling plasma can naturally produce bright flares accompanied by rotations in the synchrotron polarization vector. With particle-in-cell simulations of reconnection, we find that flares are powered by efficient particle acceleration at the interface of merging magnetic flux ropes, or "plasmoids". The accelerated particles stream through the post-merger plasmoid towards the observer, thus progressively illuminating regions with varying plane-of-sky field direction, and so leading to a rotation in the observed polarization vector. Our results provide evidence for magnetic reconnection as the physical cause of high-energy flares from the relativistic jets of blazars (which recent observations have shown to be frequently associated with polarization rotations), and provide a first-principle physical mechanism for such flares.


title: Ultraviolet-Based Science in the Solar System: Advances and Next Steps

abstract: We review the importance of recent UV observations of solar system targets and discuss the need for further measurements, instrumentation and laboratory work in the coming decade. In the past decade, numerous important advances have been made in solar system science using ultraviolet (UV) spectroscopic techniques. Formerly used nearly exclusively for studies of giant planet atmospheres, planetary exospheres and cometary emissions, UV imaging spectroscopy has recently been more widely applied. The geyser-like plume at Saturn's moon Enceladus was discovered in part as a result of UV stellar occultation observations, and this technique was used to characterize the plume and jets during the entire Cassini mission. Evidence for a similar style of activity has been found at Jupiter's moon Europa using Hubble Space Telescope (HST) UV emission and absorption imaging. At other moons and small bodies throughout the solar system, UV spectroscopy has been utilized to search for activity, probe surface composition, and delineate space weathering effects; UV photometric studies have been used to uncover regolith structure. Insights from UV imaging spectroscopy of solar system surfaces have been gained largely in the last 1-2 decades, including studies of surface composition, space weathering effects (e.g. radiolytic products) and volatiles on asteroids (e.g. [2][39][48][76][84]), the Moon (e.g. [30][46][49]), comet nuclei (e.g. [85]) and icy satellites (e.g. [38][41-44][45][47][65]). The UV is sensitive to some species, minor contaminants and grain sizes often not detected in other spectral regimes. In the coming decade, HST observations will likely come to an end. New infrastructure to bolster future UV studies is critically needed. These needs include both developmental work to help improve future UV observations and laboratory work to help interpret spacecraft data. UV instrumentation will be a critical tool on missions to a variety of targets in the coming decade, especially for the rapidly expanding application of UV reflectance investigations of atmosphereless bodies.


title: Hybrid Deep Learning Gaussian Process for Diabetic Retinopathy Diagnosis and Uncertainty Quantification

abstract: Diabetic Retinopathy (DR) is one of the microvascular complications of Diabetes Mellitus, which remains as one of the leading causes of blindness worldwide. Computational models based on Convolutional Neural Networks represent the state of the art for the automatic detection of DR using eye fundus images. Most of the current work address this problem as a binary classification task. However, including the grade estimation and quantification of predictions uncertainty can potentially increase the robustness of the model. In this paper, a hybrid Deep Learning-Gaussian process method for DR diagnosis and uncertainty quantification is presented. This method combines the representational power of deep learning, with the ability to generalize from small datasets of Gaussian process models. The results show that uncertainty quantification in the predictions improves the interpretability of the method as a diagnostic support tool. The source code to replicate the experiments is publicly available at https://github.com/stoledoc/DLGP-DR-Diagnosis.


title: Return-Oriented Programming in RISC-V

abstract: RISC-V is an open-source hardware ISA based on the RISC design principles, and has been the subject of some novel ROP mitigation technique proposals due to its open-source nature. However, very little work has actually evaluated whether such an attack is feasible assuming a typical RISC-V implementation. We show that RISC-V ROP can be used to perform Turing complete calculation and arbitrary function calls by leveraging gadgets found in a version of the GNU libc library. Using techniques such as self-modifying ROP chains and algorithmic ROP chain generation, we demonstrate the power of RISC-V ROP by creating a compiler that converts code of arbitrary complexity written in a popular Turing-complete language into RISC-V ROP chains.


title: Supervised Neural Networks for RFI Flagging

abstract: Neural network (NN) based methods are applied to the detection of radio frequency interference (RFI) in post-correlation,post-calibration time/frequency data. While calibration doesaffect RFI for the sake of this work a reduced dataset inpost-calibration is used. Two machine learning approachesfor flagging real measurement data are demonstrated usingthe existing RFI flagging technique AOFlagger as a groundtruth. It is shown that a single layer fully connects networkcan be trained using each time/frequency sample individuallywith the magnitude and phase of each polarization and Stokesvisibilities as features. This method was able to predict aBoolean flag map for each baseline to a high degree of accuracy achieving a Recall of 0.69 and Precision of 0.83 and anF1-Score of 0.75.


title: Aggregate Analytic Window Query over Spatial Data

abstract: Analytic window query is a commonly used query in the relational databases. It answers the aggregations of data over a sliding window. For example, to get the average prices of a stock for each day. However, it is not supported in the spatial databases. Because the spatial data are not in a one-dimension space, there is no straightforward way to extend the original analytic window query to spatial databases. But these queries are useful and meaningful. For example, to find the average number of visits for all the POIs in the circle with a fixed radius for each POI as the centre. In this paper, we define the aggregate analytic window query over spatial data and propose algorithms for grid index and tree-index. We also analyze the complexity of the algorithms to prove they are efficient and practical.


title: Monte Carlo Simulations of BFSS and IKKT Matrix Models

abstract: In this thesis, we studied the bosonic BFSS and IKKT matrix models using Monte Carlo simulations. First, we explored some toy models to check the validity of the numerical simulations. Then we simulated the BFSS matrix model using Hamiltonian Monte Carlo (HMC) algorithm. In the BFSS matrix model, we used the Polyakov loop as an order parameter to investigate the large-N behaviour of this model at different temperatures. Our simulations confirmed that the model exhibits a confinement-deconfinement phase transition as the temperature of the system is varied. Besides the Polyakov loop, other observables such as internal energy and extent of space were also computed. In the bosonic IKKT model, we studied the spontaneous symmetry breaking (SSB) of SO(10) symmetry using the moment of inertia tensor and found that there is no SSB of SO(10) symmetry in this model. Besides the eigenvalues of the moment of inertia tensor, other observable such as extent of spacetime was also computed. We also studied the simulation theory of the phase-quenched IKKT model.


title: MeV-scale performance of water-based and pure liquid scintillator detectors

abstract: This paper presents studies of the performance of water-based liquid scintillator in both 1-kt and 50-kt detectors. Performance is evaluated in comparison to both pure water Cherenkov detectors and a nominal model for pure scintillator detectors. Performance metrics include energy, vertex, and angular resolution, along with a metric for ability to separate the Cherenkov from the scintillation signal, as being representative of various particle identification capabilities that depend on the Cherenkov / scintillation ratio. We also modify the time profile of scintillation light to study the same performance metrics as a function of rise and decay time. We go on to interpret these results in terms of their impact on certain physics goals, such as solar neutrinos and the search for Majorana neutrinos. This work supports and validates previous results, and the assumptions made therein, by using a more complete detector model and full reconstruction. We confirm that a high-coverage, 50-kt detector would be capable of better than 10 (1)\% precision on the CNO neutrino flux with a WbLS (pure LS) target in 5 years of data taking. A 1-kt LS detector, with a conservative 50\% fiducial volume of 500~t, can achieve a better than 5\% detection. Using the liquid scintillator model, we find a sensitivity into the normal hierarchy region for Majorana neutrinos, with half life sensitivity of $T^{0\nu\beta\beta}_{1/2} > 1.4 \times 10^{28}$ years at 90\% CL for 10 years of data taking with a Te-loaded target.


title: Linear Independence Of Some Irrational Numbers

abstract: This note presents an analytic technique for proving the linear independence of certain small subsets of real numbers over the rational numbers. The applications of this test produce simple linear independence proofs for the subsets of triples $\{1, e, \pi\}$, $\{1, e, \pi^{-1}\}$, and $\{1, \pi^r, \pi^s\}$, where $1\leq r<s $ are fixed integers.


title: Dark Matter Spectra from the Electroweak to the Planck Scale

abstract: We compute the decay spectrum for dark matter (DM) with masses above the scale of electroweak symmetry breaking, all the way to the Planck scale. For an arbitrary hard process involving a decay to the unbroken standard model, we determine the prompt distribution of stable states including photons, neutrinos, positrons, and antiprotons. These spectra are a crucial ingredient in the search for DM via indirect detection at the highest energies as being probed in current and upcoming experiments including IceCube, HAWC, CTA, and LHAASO. Our approach improves considerably on existing methods. For example, we include all relevant electroweak interactions. The importance of these effects grow with DM mass, and by an EeV our spectra can differ by orders of magnitude from existing results.


title: Pairing in magic-angle twisted bilayer graphene: role of phonon and plasmon umklapp

abstract: Identifying the microscopic mechanism for superconductivity in magic-angle twisted bilayer graphene (MATBG) is an outstanding open problem. While MATBG exhibits a rich phase-diagram, driven partly by the strong interactions relative to the electronic bandwidth, its single-particle properties are unique and likely play an important role in some of the phenomenological complexity. Some of the salient features include an electronic bandwidth smaller than the characteristic phonon bandwidth and a non-trivial structure of the underlying Bloch wavefunctions. We perform a theoretical study of the cooperative effects due to phonons and plasmons on pairing in order to disentangle the distinct role played by these modes on superconductivity. We consider a variant of MATBG with an enlarged number of fermion flavors, $N \gg 1$, where the study of pairing instabilities reduces to the conventional (weak-coupling) Eliashberg framework. In particular, we show that certain umklapp processes involving mini-optical phonon modes, which arise physically as a result of the folding of the original acoustic branch of graphene due to the moir\'e superlattice structure, contribute significantly towards enhancing pairing. We also investigate the role played by the dynamics of the screened Coulomb interaction on pairing, which leads to an enhancement in a narrow window of fillings, and study the effect of external screening due to a metallic gate on superconductivity. At strong coupling the dynamical pairing interaction leaves a spectral mark in the single particle tunneling density of states. We thus predict such features will appear at specific frequencies of the umklapp phonons corresponding to the sound velocity of graphene times an integer multiple of the Brillouin zone size.


title: Discrete and higher-form symmetries in SCFTs from wrapped M5-branes

abstract: We analyze topological mass terms of BF type arising in supersymmetric M-theory compactifications to $AdS_5$. These describe spontaneously broken higher-form gauge symmetries in the bulk. Different choices of boundary conditions for the BF terms yield dual field theories with distinct global discrete symmetries. We discuss in detail these symmetries and their 't Hooft anomalies for 4d $\mathcal N = 1$ SCFTs arising from M5-branes wrapped on a Riemann surface without punctures, including theories from M5-branes at a $\mathbb Z_2$ orbifold singularity. The anomaly polynomial is computed via inflow and contains background fields for discrete global 0-, 1-, and 2-form symmetries and continuous 0-form symmetries, as well as axionic background fields. The latter are properly interpreted in the context of anomalies in the space of coupling constants.


title: Considerations for Atmospheric Retrieval of High-Precision Brown Dwarf Spectra

abstract: Isolated brown dwarfs provide remarkable laboratories for understanding atmospheric physics in the low-irradiation regime, and can be observed more precisely than exoplanets. As such, they provide a glimpse into the future of high-SNR observations of exoplanets. In this work, we investigate several new considerations that are important for atmospheric retrievals of high-quality thermal emission spectra of sub-stellar objects. We pursue this using an adaptation of the HyDRA atmospheric retrieval code. We propose a parametric pressure-temperature (P-T) profile for brown dwarfs consisting of multiple atmospheric layers, parameterised by the temperature change across each layer. This model allows the steep temperature gradient of brown dwarf atmospheres to be accurately retrieved while avoiding commonly-encountered numerical artefacts. The P-T model is especially flexible in the photosphere, which can reach a few tens of bar for T-dwarfs. We demonstrate an approach to include model uncertainties in the retrieval, focusing on uncertainties introduced by finite spectral and vertical resolution in the atmospheric model used for retrieval (~8\% in the present case). We validate our retrieval framework by applying it to a simulated data set and then apply it to the HST/WFC3 spectrum of the T-dwarf 2MASS J2339+1352. We retrieve sub-solar abundances of H2O and CH4 in the object at ~0.1 dex precision. Additionally, we constrain the temperature structure to within ~100 K in the photosphere. Our results demonstrate the promise of high-SNR spectra to provide high-precision abundance estimates of sub-stellar objects.


title: Delayed Rebounds in the Two-Ball Bounce Problem

abstract: In the classroom demonstration of a tennis ball dropped on top of a basketball, the surprisingly high bounce of the tennis ball is typically explained using momentum conservation for elastic collisions, with the basketball-floor collision treated as independent from the collision between the two balls. This textbook explanation is extended to inelastic collisions by including a coefficient of restitution. This independent contact model (ICM), as reviewed in this paper, is accurate for a wide variety of cases, even when the collisions are not truly independent. However, it is easy to explore situations that are not explained by the ICM, such as swapping the tennis ball for a ping-pong ball. In this paper, we study the conditions that lead to a "delayed rebound effect," a small first bounce followed by a higher second bounce, using techniques accessible to an undergraduate student. The dynamical model is based on the familiar solution of the damped harmonic oscillator. We focus on making the equations of motion dimensionless for numerical simulation, and reducing the number of parameters and initial conditions to emphasize universal behavior. The delayed rebound effect is found for a range of parameters, most commonly in cases where the first bounce is lower than the ICM prediction.


title: Boosted Neutrinos and Relativistic Dark Particles as Messengers from Reheating

abstract: Usually information from early eras such as reheating is hard to come by. In this paper we argue that, given the right circumstances, right-handed sterile neutrinos decaying to left-handed active ones at relatively late times can carry information from reheating by propagating freely over the thermal history. For not too small mixing angles, suitable right-handed neutrino masses are around ${\cal O}$(MeV-GeV). We identify the typical spectra and argue that they provide information on the ratio of the inflaton mass to the reheating temperature. This primordial neutrino signal can be strong enough that it can be detected in IceCube. More speculatively, for a reheating temperature and inflaton mass satisfying $T_R={\cal O}(1-100)\, {\rm MeV}$, and $m_\phi={\cal O}(10^{16-19})\,$GeV they may even explain the observed PeV events. Also more general relativistic dark particles can play the role of such messengers, potentially not only allowing for the PeV events but also alleviating the $H_0$-tension .


title: Using Diffuse Scattering to Observe X-Ray-Driven Nonthermal Melting

abstract: We present results from the SPring-8 Angstrom Compact free electron LAser (SACLA) XFEL facility, using a high intensity ($\sim\!10^{20}\,$W/cm$^2$) X-ray pump X-ray probe scheme to observe changes in the ionic structure of silicon induced by X-ray heating of the electrons. By avoiding Laue spots in the scattering signal from a single crystalline sample, we observe a rapid rise in diffuse scattering, which we attribute to a loss of lattice order and a transition to a liquid state within 100 fs of irradiation, a timescale which agrees well with first principles simulations, but is faster than that predicted by purely inertial behavior. This method is capable of observing liquid scattering without masking or filtering of signal from the ambient solid, allowing the liquid structure to be measured throughout and beyond the phase change.


title: Contrasting lattice geometry dependent versus independent quantities: ramifications for Berry curvature, energy gaps, and dynamics

abstract: In the tight-binding description of electronic, photonic, or cold atomic dynamics in a periodic lattice potential, particle motion is described in terms of hopping amplitudes and potentials on an abstract network of discrete sites corresponding to physical orbitals in the lattice. The physical attributes of the orbitals, including their locations in three-dimensional space, are independent pieces of information. In this paper we identify a notion of geometry-independence: any physical quantity or observable that depends only on the tight-binding parameters (and not on the explicit information about the orbital geometry) is said to be "geometry-independent." The band structure itself, and for example the Chern numbers of the bands in a two-dimensional system, are geometry-independent, while the Bloch-band Berry curvature is geometry-dependent. Careful identification of geometry-dependent versus independent quantities can be used as a novel principle for constraining a variety of results. By extending the notion of geometry-independence to certain classes of interacting systems, where the many-body energy gap is evidently geometry-independent, we shed new light on a hypothesized relation between many-body energy gaps of fractional Chern insulators and the uniformity of Bloch band Berry curvature in the Brillouin zone. We furthermore explore the geometry-dependence of semiclassical wave packet dynamics, and use this principle to show how two different types of Hall response measurements may give markedly different results due to the fact that one is geometry-dependent, while the other is geometry-independent. Similar considerations apply for anomalous thermal Hall response, in both electronic and spin systems.


title: Gravitational positivity bounds

abstract: We study the validity of positivity bounds in the presence of a massless graviton, assuming the Regge behavior of the amplitude. Under this assumption, the problematic $t$-channel pole is canceled with the UV integral of the imaginary part of the amplitude in the dispersion relation, which gives rise to finite corrections to the positivity bounds. We find that low-energy effective field theories (EFT) with "wrong" sign are generically allowed. The allowed amount of the positivity violation is determined by the Regge behavior. This violation is suppressed by $M_{\rm pl}^{-2}\alpha'$ where $\alpha'$ is the scale of Reggeization. This implies that the positivity bounds can be applied only when the cutoff scale of EFT is much lower than the scale of Reggeization. We then obtain the positivity bounds on scalar-tensor EFT at one-loop level. Implications of our results on the degenerate higher-order scalar-tensor (DHOST) theory are also discussed.


title: $\chi$aro$\nu$: a tool for neutrino flux generation from WIMPs

abstract: Indirect searches for signatures of corpuscular dark matter have been performed using all cosmic messengers: gamma rays, cosmic rays, and neutrinos. The search for dark matter from neutrinos is of particular importance since they are the only courier that can reach detectors from dark matter processes in dense environments, such as the core of the Sun or Earth, or from the edge of the observable Universe. In this work, we introduce $\chi$aro$\nu$, a software package that, in the spirit of its mythological Greek namesake $\chi \acute\alpha \rho \omega \nu $, bridges the dark sector and Standard Model by predicting neutrino fluxes from different celestial dark matter agglomerations. The flux at the point of production is either computed internally by $\chi$aro$\nu$ or is taken from user supplied tables. $\chi$aro$\nu$ then propagates this flux through vacuum or dense media and returns the expected neutrino yield at an observer's location. In developing $\chi$aro$\nu$, we have revisited and updated the production of neutrinos in dense media, updated the propagation of high-energy neutrinos, and studied the sources of uncertainty in neutrino transport. This package is coupled to a new calculation that includes electroweak corrections resulting in the most up-to-date and complete repository of neutrino fluxes from dark matter decay and annihilation over the energy range of 1 GeV to 10 PeV coming from the Earth, the Sun, and the Galactic halo.


title: Evidence of attraction between charge-carriers in a doped Mott insulator

abstract: Recent progress in optically trapped ultracold atomic gases is now making it possible to access microscopic observables in doped Mott insulators, which are the parent states of high-temperature superconductors. This makes it possible to address longstanding questions about the temperature scales at which attraction between charge carriers are present, and their mechanism. Controllable theoretical results for this problem are not available at low temperature due to the sign problem. In this work, we employ worm-algorithm Monte Carlo to obtain completely unbiased results for two charge carriers in a Mott insulator. Our method gives access to lower temperatures than what is currently possible in experiments, and provides evidence for attraction between dopants at a temperature scale that is now feasible in ultracold atomic systems. We also report on spin-correlations in the presence of charge carriers, which are directly comparable to experiments.


title: Delayed Photons from Binary Evolution Help Reionize the Universe

abstract: High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (fesc) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is 70%. Moreover, simulations so far yield values of fesc falling only on the lower end of the roughly 10-20% range, the amount presumed necessary to reionize the Universe. Analyzing a high-resolution (4 pc) cosmological radiation hydrodynamic simulation we study how fesc changes when we include two different products of binary stellar evolution - stars stripped of their hydrogen envelopes and massive blue stragglers. Both produce significant amounts of ionizing photons 10-200 Myr after each starburst. We find the relative importance of these photons are amplified with respect to escaped ionizing photons, because peaks in star formation rates (SFRs) and fesc are often out of phase by this 10-200 Myr. Additionally, low mass, bursty galaxies emit Lyman continuum radiation primarily from binary products when SFRs are low. Observations of these galaxies by the James Webb Space Telescope could provide crucial information on the evolution of binary stars as a function of redshift. Overall, including stripped stars and massive blue stragglers increases our photon-weighted mean escape fraction by around 13% and 10%, respectively, resulting in a mean fesc of 17%. Our results emphasize that using updated stellar population synthesis models with binary stellar evolution provides a more sound physical basis for stellar reionization.


title: Impact of fermionic operators on the Higgs width measurement

abstract: The off-shell Higgs production in pp -> ZZ at LHC provides at present the most direct measurement of the Higgs width in the absence of beyond the standard model contributions. Here, we analyze the impact of anomalous Z couplings to fermions on the Higgs width determination. We show that, despite these couplings being strongly constrained by the available electroweak precision data, they can substantially affect the Higgs width determination at the LHC Runs 2 and 3. Conversely, in larger integrated luminosities runs, such as those foreseen at the high luminosity LHC and high energy LHC setups, the effect of such anomalous interactions in the Higgs width measurement is minimal.


title: Solving grain size inconsistency between ALMA polarization and VLA continuum in the Ophiuchus IRS 48 protoplanetary disk

abstract: The protoplanetary disk around Ophiuchus IRS 48 shows an azimuthally asymmetric dust distribution in (sub-)millimeter observations, which is interpreted as a vortex, where millimeter/centimeter-sized particles are trapped at the location of the continuum peak. In this paper, we present 860 $\mu$m ALMA observations of polarized dust emission of this disk. The polarized emission was detected toward a part of the disk. The polarization vectors are parallel to the disk minor axis, and the polarization fraction was derived to be $1-2$\%. These characteristics are consistent with models of self-scattering of submillimeter-wave emission, which indicate a maximum grain size of $\sim100$ $\mu$m. However, this is inconsistent with the previous interpretation of millimeter/centimeter dust particles being trapped by a vortex. To explain both, ALMA polarization and previous ALMA and VLA observations, we suggest that the thermal emission at 860 $\mu$m wavelength is optically thick ($\tau_{\rm abs}\sim7.3$) at the dust trap with the maximum observable grain size of $\sim100$ $\mu$m rather than an optically thin case with $\sim$ cm dust grains. We note that we cannot rule out that larger dust grains are accumulated near the midplane if the 860 $\mu$m thermal emission is optically thick.


title: Fully exclusive heavy quark-antiquark pair production from a colourless initial state at NNLO in QCD

abstract: We present a local subtraction scheme for computing next-to-next-to-leading order QCD corrections to the production of a massive quark-antiquark pair from a colourless initial state. The subtraction terms are built following the CoLoRFulNNLO method and refined in such a way that their integration gives rise to compact, fully analytic expressions. All ingredients necessary for a numerical implementation of our subtraction scheme are provided in detail. As an example, we calculate the fully differential decay rate of the Standard Model Higgs boson to massive bottom quarks at next-to-next-to-leading order accuracy in perturbative QCD.


title: Anomalies and accidental symmetries: charging the scalar leptoquark under $L_\mu$-$L_\tau$

abstract: While the $S_3$ scalar leptoquark presents a possible tree-level explanation of the $b \to s \ell \ell$ flavour anomalies, it suffers from two conceptual problems which are often disregarded by model-builders. Firstly, the quantum numbers of the $S_3$ allow for a renormalisable diquark operator that would trigger rapid proton decay unless its coupling were tuned away. Secondly, one expects the leptoquark to have generic couplings to leptons, which require tuning to avoid stringent experimental bounds on lepton flavour violation. By gauging a $U(1)$ current that acts as $L_\mu - L_\tau$ on the Standard Model (SM) fermions, and under which the leptoquark has charge $-1$, one can remedy both these problems. The additional $U(1)$, which is spontaneously broken at some high scale, is associated with a massive $Z^\prime$ gauge boson and a scalar SM singlet $\Phi$, which play no direct role in mediating the anomalous $B$ meson decays. By computing one- and two-loop mass corrections, we show that this pair of particles can be hidden away at much higher mass scales without destabilising either the Higgs or the leptoquark masses. The only low-energy relic of gauging $L_\mu - L_\tau$ is thus the accidental global symmetry structure of the lagrangian. On the other hand, we find quite generally that an $S_3$ leptoquark that mediates the $b \to s \ell \ell$ anomalies cannot be much heavier than a few TeV without itself inducing large Higgs mass corrections.


title: Revisiting the AdS Boson Stars: the Mass-Charge Relations

abstract: Motivated by the recent progress in solving the large charge sector of conformal field theories, we revisit the mass-charge relation of boson stars asymptotic to global AdS. We construct and classify a large number of electrically charged boson star solutions in a toy model and two supergravity models arising from the $SU(3)$ and $U(1)^4$ truncation of $D=4$ $SO(8)$ gauged maximal supergravity. We find a simple ansatz for the chemical potential that can fit the numerical data in striking accuracy for the full range of charge. Combining with the first law of thermodynamics, we can then evaluate the mass as a function of the charge and obtain the free energy in the fixed charge ensemble. We show that in the toy model, the ground state can be either the extremal RN black hole or the boson stars depending on the parameter region. For the $SU(3)$ truncation, there always exists a boson star that has smaller free energy than the extremal RN black hole, in contrast to the $U(1)^4$ model where the extremal RN black hole is always the ground state. In all models, for boson star solutions with arbitrarily large charge, we show that the large charge expansion of the mass reproduces the same structure exhibited in the CFT side.


title: The LPM effect in sequential bremsstrahlung: nearly complete results for QCD

abstract: The splitting processes of bremsstrahlung and pair production in a medium are coherent over large distances in the very high energy limit, which leads to a suppression known as the Landau-Pomeranchuk-Migdal (LPM) effect. We continue study of the case when the coherence lengths of two consecutive splitting processes overlap (which is important for understanding corrections to standard treatments of the LPM effect in QCD), avoiding soft-emission approximations. Previous work has computed overlap effects for double splitting $g \to gg \to ggg$. To make use of those results, one also needs calculations of related virtual loop corrections to single splitting $g \to gg$ in order to cancel severe (power-law) infrared (IR) divergences. This paper provides calculations of nearly all such processes involving gluons and discusses how to organize the results to demonstrate the cancellation. In the soft emission limit, our results reproduce the known double-log behavior of earlier authors who worked in leading-log approximation. We also present a first (albeit numerical and not yet analytic) investigation of sub-leading, single IR logarithms. Ultraviolet divergences appearing in our calculations correctly renormalize the coupling $\alpha_{\rm s}$ in the usual LPM result for leading-order $g \to gg$.


title: Probing Quantum Speed Limits with Ultracold Gases

abstract: Quantum Speed Limits (QSLs) rule the minimum time for a quantum state to evolve into a distinguishable state in an arbitrary physical process. These fundamental results constrain a notion of distance travelled by the quantum state, known as the Bures angle, in terms of the speed of evolution set by nonadiabatic energy fluctuations. We theoretically propose how to measure QSLs in an ultracold quantum gas confined in a time-dependent harmonic trap. In this highly-dimensional system of continuous variables, quantum tomography is prohibited. Yet, QSLs can be probed whenever the dynamics is self-similar by measuring as a function of time the cloud size of the ultracold gas. This makes possible to determine the Bures angle and energy fluctuations, as we discuss for various ultracold atomic systems.


title: Mapping the inhomogeneous Universe with Standard Sirens: Degeneracy between inhomogeneity and modified gravity theories

abstract: The detection of gravitational waves (GWs) and an accompanying electromagnetic (E/M) counterpart have been suggested as a future probe for cosmology and theories of gravity. In this paper, we present calculations of the luminosity distance of sources taking into account inhomogeneities in the matter distribution that are predicted in numerical simulations of structure formation. In addition, we show that inhomogeneities resulting from clustering of matter can mimic certain classes of modified gravity theories, or other effects that dampen GW amplitudes, and deviations larger than $\delta \nu \sim \mathcal{O}(0.1)\ (99\%\ \rm{C.L.})$ to the extra friction term $\nu$, from zero, would be necessary to distinguish them. For these, we assume mock GWs sources, with known redshift, based on binary population synthesis models, between redshifts $z=0$ and $z=5$. We show that future GW detectors, like Einstein Telescope or Cosmic Explorer, will be needed for strong constraints on the inhomogeneity parameters and breaking the degeneracy between modified gravity effects and matter anisotropies by measuring $\nu$ at $5 \%$ and $1 \%$ level with $100$ and $350$ events respectively.


title: Precision Higgs Couplings in Neutral Naturalness Models: an Effective Field Theory Approach

abstract: The Higgs sector in neutral naturalness models provides a portal to the hidden sectors, and thus measurements of Higgs couplings at current and future colliders play a central role in constraining the parameter space of the model. We investigate a class of neutral naturalness models, in which the Higgs boson is a pseudo-Goldstone boson from the universal SO(N)/SO(N-1) coset structure. Integrating out the radial mode from the spontaneous global symmetry breaking, we obtain various dimension-six operators in the Standard Model effective field theory, and calculate the low energy Higgs effective potential with radiative corrections included. We perform a chi-square fit to the Higgs coupling precision measurements at current and future colliders and show that the new physics scale could be explored up to 2.7 (2.8) TeV without (with) the Higgs invisible decay channels at future Higgs factories.


title: Mass and rate of hierarchical black hole mergers in young, globular and nuclear star clusters

abstract: Hierarchical mergers are one of the distinctive signatures of binary black hole (BBH) formation through dynamical evolution. Here, we present a fast semi-analytic approach to simulate hierarchical mergers in nuclear star clusters (NSCs), globular clusters (GCs) and young star clusters (YSCs). Hierarchical mergers are more common in NSCs than they are in both GCs and YSCs, because of the different escape velocity. The mass distribution of hierarchical BBHs strongly depends on the properties of first-generation BBHs, such as their progenitor's metallicity. In our fiducial model, we form black holes (BHs) with masses up to $\sim{}10^3$ M$_\odot$ in NSCs and up to $\sim{}10^2$ M$_\odot$ in both GCs and YSCs. When escape velocities in excess of 100 km~s$^{-1}$ are considered, BHs with mass $>10^3$ M$_\odot$ are allowed to form in NSCs. Hierarchical mergers lead to the formation of BHs in the pair instability mass gap and intermediate-mass BHs, but only in metal-poor environments. The local BBH merger rate in our models ranges from $\sim{}10$ to $\sim{} 60$ Gpc$^{-3}$ yr$^{-1}$; hierarchical BBHs in NSCs account for $\sim{}10^{-2}- 0.2$ Gpc$^{-3}$ yr$^{-1}$, with a strong upper limit of $\sim{}10$ Gpc$^{-3}$ yr$^{-1}$. When comparing our models with the second gravitational-wave transient catalog, we find that multiple formation channels are favored to reproduce the observed BBH population.


title: Cosmic inflation as a renormalization-group flow: the running of power spectra in quantum gravity

abstract: We study the running of power spectra in inflationary cosmology as a renormalization-group flow from the de Sitter fixed point. The beta function is provided by the equations of the background metric. The spectra of the scalar and tensor fluctuations obey RG evolution equations with vanishing anomalous dimensions in the superhorizon limit. By organizing the perturbative expansion in terms of leading and subleading logs, we calculate the spectral indices, their runnings, the runnings of the runnings, etc., to the next-to-leading log order in quantum gravity with fakeons (i.e., the theory $R+R^{2}+C^{2}$ with the fakeon prescription/projection for $C^{2}$). We show that these quantities are related to the spectra in a universal way. We also compute the first correction to the relation $r=-8n_{T}$ and provide a number of quantum gravity predictions that can be hopefully tested in the forthcoming future.


title: Constraining Delay Time Distribution of Binary Neutron Star Mergers from Host Galaxy Properties

abstract: Gravitational wave (GW) observatories are discovering binary neutron star mergers (BNSMs), and in at least one event we were able to track it down in multiple wavelengths of light, which allowed us to identify the host galaxy. Using a catalogue of local galaxies with inferred star formation histories and adopting a BNSM delay time distribution (DTD) model, we investigate the dependence of BNSM rate on an array of galaxy properties. Compared to the intrinsic property distribution of galaxies, that of BNSM host galaxies is skewed toward galaxies with redder colour, lower specific star formation rate, higher luminosity, and higher stellar mass, reflecting the tendency of higher BNSM rates in more massive galaxies. We introduce a formalism to efficiently make forecast on using host galaxy properties to constrain DTD models. We find comparable constraints from the dependence of BNSM occurrence distribution on galaxy colour, specific star formation rate, and stellar mass, all better than those from dependence on $r$-band luminosity. The tightest constraints come from using individual star formation histories of host galaxies, which reduces the uncertainties on DTD parameters by a factor of three or more. Substantially different DTD models can be differentiated with about 10 BNSM detections. To constrain DTD parameters at 10\% precision level requires about one hundred detections, achievable with GW observations on a decade time scale.


title: A Renyi Quantum Null Energy Condition: Proof for Free Field Theories

abstract: The Quantum Null Energy Condition (QNEC) is a lower bound on the stress-energy tensor in quantum field theory that has been proved quite generally. It can equivalently be phrased as a positivity condition on the second null shape derivative of the relative entropy $S_{\text{rel}}(\rho||\sigma)$ of an arbitrary state $\rho$ with respect to the vacuum $\sigma$. The relative entropy has a natural one-parameter family generalization, the Sandwiched Renyi divergence $S_n(\rho||\sigma)$, which also measures the distinguishability of two states for arbitrary $n\in[1/2,\infty)$. A Renyi QNEC, a positivity condition on the second null shape derivative of $S_n(\rho||\sigma)$, was conjectured in previous work. In this work, we study the Renyi QNEC for free and superrenormalizable field theories in spacetime dimension $d>2$ using the technique of null quantization. In the above setting, we prove the Renyi QNEC in the case $n>1$ for arbitrary states. We also provide counterexamples to the Renyi QNEC for $n<1$.


title: Dark Energy Survey Year 1 Results: Constraining Baryonic Physics in the Universe

abstract: Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3$\times$2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modeling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude ($Q_1$) generally reflects the strength of baryon feedback. With the upper limit of $Q_1$ prior being bound by the Illustris feedback scenarios, we reach $\sim 20\%$ improvement in the constraint of $S_8=\sigma_8(\Omega_{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3$\times$2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative $Q_1$ prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than $2 \sigma$.


title: Double Copy for Celestial Amplitudes

abstract: Celestial amplitudes which use conformal primary wavefunctions rather than plane waves as external states offer a novel opportunity to study properties of amplitudes with manifest conformal covariance and give insight into a potential holographic celestial CFT at the null boundary of asymptotically flat space. Since translation invariance is obscured in the conformal basis, features of amplitudes that heavily rely on it appear to be lost. Among these are the remarkable relations between gauge theory and gravity amplitudes known as the double copy. Nevertheless, properties of amplitudes reflecting fundamental aspects of the perturbative regime of quantum field theory are expected to survive a change of basis. Here we show that there exists a well-defined procedure for a celestial double copy. This requires a generalization of the usual squaring of numerators which entails first promoting them to generalized differential operators acting on external wavefunctions, and then squaring them. We demonstrate this procedure for three and four point celestial amplitudes, and give an argument for its validity to all multiplicities.


title: Quantum Geometric Contributions to the BKT Transition: Beyond Mean Field Theory

abstract: We study quantum geometric contributions to the Berezinskii-Kosterlitz-Thouless (BKT) transition temperature, $T_{\mathrm{BKT}}$, in the presence of fluctuations beyond BCS theory. Because quantum geometric effects become progressively more important with stronger pairing attraction, a full understanding of 2D multi-orbital superconductivity requires the incorporation of preformed pairs. We find it is through the effective mass of these pairs that quantum geometry enters the theory and this suggests that the quantum geometric effects are present in the non-superconducting pseudogap phase as well. Increasing these geometric contributions tends to raise $T_{\mathrm{BKT}}$ which then competes with fluctuation effects that generally depress it. We argue that a way to physically quantify the magnitude of these geometric terms is in terms of the ratio of the pairing onset temperature $T^*$ to $T_{\mathrm{BKT}}$. Our paper calls attention to an experimental study demonstrating how both temperatures and, thus, their ratio may be currently accessible. They can be extracted from the same voltage-current measurements which are generally used to establish BKT physics. We use these observations to provide rough preliminary estimates of the magnitude of the geometric contributions in, for example, magic angle twisted bilayer graphene.

categories: cond-mat.supr-con cond-mat.mes-hall cond-mat.quant-gas cond-mat.str-el


title: Relative Acceleration Noise Mitigation for Nanocrystal Matter-wave Interferometry: Application to Entangling Masses via Quantum Gravity

abstract: Matter wave interferometers with large momentum transfers, irrespective of specific implementations, will face a universal dephasing due to relative accelerations between the interferometric mass and the associated apparatus. Here we propose a solution that works even without actively tracking the relative accelerations: putting both the interfering mass and its associated apparatus in a freely falling capsule, so that the strongest inertial noise components vanish due to the equivalence principle. In this setting, we investigate two of the most important remaining noise sources: (a) the non-inertial jitter of the experimental setup and (b) the gravity-gradient noise. We show that the former can be reduced below desired values by appropriate pressures and temperatures, while the latter can be fully mitigated in a controlled environment. We finally apply the analysis to a recent proposal for testing the quantum nature of gravity [S. Bose et. al. Phys. Rev. Lett 119, 240401 (2017)] through the entanglement of two masses undergoing interferometry. We show that the relevant entanglement witnessing is feasible with achievable levels of relative acceleration noise.


title: Dynamic Defense Against Byzantine Poisoning Attacks in Federated Learning

abstract: Federated learning, as a distributed learning that conducts the training on the local devices without accessing to the training data, is vulnerable to Byzatine poisoning adversarial attacks. We argue that the federated learning model has to avoid those kind of adversarial attacks through filtering out the adversarial clients by means of the federated aggregation operator. We propose a dynamic federated aggregation operator that dynamically discards those adversarial clients and allows to prevent the corruption of the global learning model. We assess it as a defense against adversarial attacks deploying a deep learning classification model in a federated learning setting on the Fed-EMNIST Digits, Fashion MNIST and CIFAR-10 image datasets. The results show that the dynamic selection of the clients to aggregate enhances the performance of the global learning model and discards the adversarial and poor (with low quality models) clients.


title: Regression-based imputation of explanatory discrete missing data

abstract: Imputation of missing values is a strategy for handling non-responses in surveys or data loss in measurement processes, which may be more effective than ignoring them. When the variable represents a count, the literature dealing with this issue is scarce. Likewise, if problems of over- or under-dispersion are observed, generalisations of the Poisson distribution are recommended for carrying out imputation. In order to assess the performance of various regression models in the imputation of a discrete variable compared to classical counting models, this work presents a comprehensive simulation study considering a variety of scenarios and real data. To do so we compared the results of estimations using only complete data, and using imputations based on the Poisson, negative binomial, Hermite, and COMPoisson distributions, and the ZIP and ZINB models for excesses of zeros. The results of this work reveal that the COMPoisson distribution provides in general better results in any dispersion scenario, especially when the amount of missing information is large. When the variable presenting missing values is a count, the most widely used method is to assume that a classical Poisson model is the best alternative to impute the missing counts; however, in real-life research this assumption is not always correct, and it is common to find count variables exhibiting overdispersion or underdispersion, for which the Poisson model is no longer the best to use in imputation. In several of the scenarios considered the performance of the methods analysed differs, something which indicates that it is important to analyse dispersion and the possible presence of excess zeros before deciding on the imputation method to use. The COMPoisson model performs well as it is flexible regarding the handling of counts with characteristics of over- and under-dispersion, as well as with equidispersion.


title: A Flexible and Modular Body-Machine Interface for Individuals Living with Severe Disabilities

abstract: This paper presents a control interface to translate the residual body motions of individuals living with severe disabilities, into control commands for body-machine interaction. A custom, wireless, wearable multi-sensor network is used to collect motion data from multiple points on the body in real-time. The solution proposed successfully leverage electromyography gesture recognition techniques for the recognition of inertial measurement units-based commands (IMU), without the need for cumbersome and noisy surface electrodes. Motion pattern recognition is performed using a computationally inexpensive classifier (Linear Discriminant Analysis) so that the solution can be deployed onto lightweight embedded platforms. Five participants (three able-bodied and two living with upper-body disabilities) presenting different motion limitations (e.g. spasms, reduced motion range) were recruited. They were asked to perform up to 9 different motion classes, including head, shoulder, finger, and foot motions, with respect to their residual functional capacities. The measured prediction performances show an average accuracy of 99.96% for able-bodied individuals and 91.66% for participants with upper-body disabilities. The recorded dataset has also been made available online to the research community. Proof of concept for the real-time use of the system is given through an assembly task replicating activities of daily living using the JACO arm from Kinova Robotics.


title: Computation of free boundary minimal surfaces via extremal Steklov eigenvalue problems

abstract: Recently Fraser and Schoen showed that the solution of a certain extremal Steklov eigenvalue problem on a compact surface with boundary can be used to generate a free boundary minimal surface, i.e., a surface contained in the ball that has (i) zero mean curvature and (ii) meets the boundary of the ball orthogonally (doi:10.1007/s00222-015-0604-x). In this paper, we develop numerical methods that use this connection to realize free boundary minimal surfaces. Namely, on a compact surface, $\Sigma$, with genus $\gamma$ and $b$ boundary components, we maximize $\sigma_j(\Sigma,g) \ L(\partial \Sigma, g)$ over a class of smooth metrics, $g$, where $\sigma_j(\Sigma,g)$ is the $j$-th nonzero Steklov eigenvalue and $L(\partial \Sigma, g)$ is the length of $\partial \Sigma$. Our numerical method involves (i) using conformal uniformization of multiply connected domains to avoid explicit parameterization for the class of metrics, (ii) accurately solving a boundary-weighted Steklov eigenvalue problem in multi-connected domains, and (iii) developing gradient-based optimization methods for this non-smooth eigenvalue optimization problem. For genus $\gamma =0$ and $b=2,\dots, 9, 12, 15, 20$ boundary components, we numerically solve the extremal Steklov problem for the first eigenvalue. The corresponding eigenfunctions generate a free boundary minimal surface, which we display in striking images. For higher eigenvalues, numerical evidence suggests that the maximizers are degenerate, but we compute local maximizers for the second and third eigenvalues with $b=2$ boundary components and for the third and fifth eigenvalues with $b=3$ boundary components.


title: Epidemic plateau in critical SIR dynamics with non-trivial initial conditions

abstract: Containment measures implemented by some countries to suppress the spread of COVID-19 have resulted in a slowdown of the epidemic characterized by time series of daily infections plateauing over extended periods of time. We prove that such a dynamical pattern is compatible with critical Susceptible-Infected-Removed (SIR) dynamics. In traditional analyses of the critical SIR model, the critical dynamical regime is started from a single infected node. The application of containment measures to an ongoing epidemic, however, has the effect to make the system enter in its critical regime with a number of infected individuals potentially large. We describe how such non-trivial starting conditions affect the critical behavior of the SIR model. We perform a theoretical and large-scale numerical investigation of the model. We show that the expected outbreak size is an increasing function of the initial number of infected individuals, while the expected duration of the outbreak is a non-monotonic function of the initial number of infected individuals. Also, we precisely characterize the magnitude of the fluctuations associated with the size and duration of the outbreak in critical SIR dynamics with non-trivial initial conditions. Far from heard immunity, fluctuations are much larger than average values, thus indicating that predictions of plateauing time series may be particularly challenging.


title: Time-dependent variational principle of mixed matrix product states in the thermodynamic limit

abstract: We describe a time evolution algorithm for quantum spin chains whose Hamiltonians are composed of an infinite uniform left and right bulk part, and an arbitrary finite region in between. The left and right bulk parts are allowed to be different from each other. The algorithm is based on the time-dependent variational principle (TDVP) of matrix product states. It is inversion-free and very simple to adapt from an existing TDVP code for finite systems. The importance of working in the projective Hilbert space is highlighted. We study the quantum Ising model as a benchmark and an illustrative example. The spread of information after a local quench is studied in both the ballistic and the diffusive case. We also offer a derivation of TDVP directly from symplectic geometry.


title: Generative Classifiers as a Basis for Trustworthy Image Classification

abstract: With the maturing of deep learning systems, trustworthiness is becoming increasingly important for model assessment. We understand trustworthiness as the combination of explainability and robustness. Generative classifiers (GCs) are a promising class of models that are said to naturally accomplish these qualities. However, this has mostly been demonstrated on simple datasets such as MNIST and CIFAR in the past. In this work, we firstly develop an architecture and training scheme that allows GCs to operate on a more relevant level of complexity for practical computer vision, namely the ImageNet challenge. Secondly, we demonstrate the immense potential of GCs for trustworthy image classification. Explainability and some aspects of robustness are vastly improved compared to feed-forward models, even when the GCs are just applied naively. While not all trustworthiness problems are solved completely, we observe that GCs are a highly promising basis for further algorithms and modifications. We release our trained model for download in the hope that it serves as a starting point for other generative classification tasks, in much the same way as pretrained ResNet architectures do for discriminative classification.


title: Strong coupling diagnostics for multi-mode open systems

abstract: We present a new method to diagnose strong coupling in multi-mode open systems. Our method presents a non-trivial extension of exceptional point (EP) analysis employed for such systems; specifically, we show how eigenvectors can not only reproduce all the features predicted by EPs but are also able to identify the physical modes that hybridize in different regions of the strong coupling regime. As a demonstration, we apply this method to study hybridization physics in a three-mode optomechanical system and determine the parameter regime for efficient sideband cooling of the mechanical oscillator in the presence of reservoir correlations.


title: Metamodel Based Forward and Inverse Design for Passive Vibration Suppression

abstract: Aperiodic metamaterials represent a class of structural systems that are composed of different building blocks (cells), instead of a self-repeating chain of the same unit cells. Optimizing aperiodic cellular structural systems thus presents high-dimensional problems that are challenging to solve using purely high-fidelity structural optimization approaches. Specialized analytical modeling along with metamodel based optimization can provide a more tractable alternative solution approach. To this end, this paper presents a design automation framework applied to a 1D metamaterial system, namely a drill string, where vibration suppression is of utmost importance. The drill string comprises a set of nonuniform rings attached to the outer surface of a longitudinal rod. As such, the resultant system can now be perceived as an aperiodic 1D metamaterial with each ring/gap representing a cell. Despite being a 1D system, the simultaneous consideration of multiple DoF (i.e., torsional, axial, and lateral motions) poses significant computational challenges. Therefore, a transfer matrix method (TMM) is employed to analytically determine the frequency response of the drill string. A suite of neural networks (ANN) is trained on TMM samples (which present minute-scale computing costs per evaluation), to model the frequency response. ANN-based optimization is then performed to minimize mass subject to constraints on the gap between consecutive resonance peaks in one case, and minimizing this gap in the second case, leading to crucial improvements over baselines. Further novel contribution occurs through the development of an inverse modeling approach that can instantaneously produce the 1D metamaterial design with minimum mass for a given desired non-resonant frequency range. This is accomplished by using invertible neural networks, and results show promising alignment with forward solutions.


title: Extreme-K categorical samples problem

abstract: With histograms as its foundation, we develop Categorical Exploratory Data Analysis (CEDA) under the extreme-$K$ sample problem, and illustrate its universal applicability through four 1D categorical datasets. Given a sizable $K$, CEDA's ultimate goal amounts to discover by data's information content via carrying out two data-driven computational tasks: 1) establish a tree geometry upon $K$ populations as a platform for discovering a wide spectrum of patterns among populations; 2) evaluate each geometric pattern's reliability. In CEDA developments, each population gives rise to a row vector of categories proportions. Upon the data matrix's row-axis, we discuss the pros and cons of Euclidean distance against its weighted version for building a binary clustering tree geometry. The criterion of choice rests on degrees of uniformness in column-blocks framed by this binary clustering tree. Each tree-leaf (population) is then encoded with a binary code sequence, so is tree-based pattern. For evaluating reliability, we adopt row-wise multinomial randomness to generate an ensemble of matrix mimicries, so an ensemble of mimicked binary trees. Reliability of any observed pattern is its recurrence rate within the tree ensemble. A high reliability value means a deterministic pattern. Our four applications of CEDA illuminate four significant aspects of extreme-$K$ sample problems.


title: A new framework for the computation of Hessians

abstract: We investigate the computation of Hessian matrices via Automatic Differentiation, using a graph model and an algebraic model. The graph model reveals the inherent symmetries involved in calculating the Hessian. The algebraic model, based on Griewank and Walther's state transformations, synthesizes the calculation of the Hessian as a formula. These dual points of view, graphical and algebraic, lead to a new framework for Hessian computation. This is illustrated by developing edge_pushing, a new truly reverse Hessian computation algorithm that fully exploits the Hessian's symmetry. Computational experiments compare the performance of edge_pushing on sixteen functions from the CUTE collection against two algorithms available as drivers of the software ADOL-C, and the results are very promising.


title: Uniqueness in Cauchy problems for diffusive real-valued strict local martingales

abstract: For a real-valued one dimensional diffusive strict local martingale,, we provide a set of smooth functions in which the Cauchy problem has a unique classical solution under a local H\"older condition. Under the weaker Engelbert-Schmidt conditions, we provide a set in which the Cauchy problem has a unique weak solution. We exemplify our results using quadratic normal volatility models and the two dimensional Bessel process.


title: Hybrid reflections from multiple x-ray scattering in epitaxial bismuth telluride topological insulator films

abstract: Epitaxial films of bismuth telluride topological insulators have received increasing attention due to potential applications in spintronic and quantum computation. One of the most important properties of epitaxial films is the presence of interface defects due to lateral lattice mismatch since electrically active defects can drastically compromise device performance. By describing hybrid reflections in hexagonal bismuth telluride films on cubic substrates, in-plane lattice mismatches were characterized with accuracy at least 20 times better than using other X-ray diffraction methods, providing clear evidence of 0.007\% lateral lattice mismatch, consistent with stress relaxation associated with van der Waals gaps in the film structure.