title: The algebraic $\alpha$-Ford tree under evolution

abstract: Null models of binary phylogenetic trees are useful for testing hypotheses on real world phylogenies. In this paper we consider phylogenies as binary trees without edge lengths together with a sampling measure and encode them as algebraic measure trees. This allows to describe the degree of similarity between actual and simulated phylogenies by focusing on the sample shape of subtrees and their subtree masses. We describe the annealed law of the statistics of subtree masses of null models, namely the branching tree, the coalescent tree, and the comb tree in more detail. Finally, we use methods from martingale problems to characterize evolving phylogenetic trees in the diffusion limit.


title: Higher Kazhdan projections, $\ell_2$-Betti numbers and Baum-Connes conjectures

abstract: We introduce higher-dimensional analogs of Kazhdan projections in matrix algebras over group $C^*$-algebras and Roe algebras. These projections are constructed in the framework of cohomology with coefficients in unitary representations and in certain cases give rise to non-trivial $K$-theory classes. We apply the higher Kazhdan projections to establish a relation between $\ell_2$-Betti numbers of a group and surjectivity of different Baum-Connes type assembly maps.


title: Non-Markovian trajectories involving future in the semi-classical path integral expression

abstract: Semiclassical path integral expression for a quantum system coupled to a harmonic bath is derived based on the stationary phase condition. It is discovered that the system path is non-Markovian. Most strikingly, the system path not only couples to its past (as in the Langevin equation), but also to its future, i.e. the equation of motion for the system is an integro-differential equation that involves all times. Numerical tests are performed to confirm that the future-involved term is indeed necessary. Because of the future-non-Markovian nature of the equation, the numerical solution cannot be obtained by iterative methods. Instead, root search algorithms must be employed.


title: A Survey of Constrained Gaussian Process Regression: Approaches and Implementation Challenges

abstract: Gaussian process regression is a popular Bayesian framework for surrogate modeling of expensive data sources. As part of a broader effort in scientific machine learning, many recent works have incorporated physical constraints or other a priori information within Gaussian process regression to supplement limited data and regularize the behavior of the model. We provide an overview and survey of several classes of Gaussian process constraints, including positivity or bound constraints, monotonicity and convexity constraints, differential equation constraints provided by linear PDEs, and boundary condition constraints. We compare the strategies behind each approach as well as the differences in implementation, concluding with a discussion of the computational challenges introduced by constraints.


title: Managing Consensus-Based Cooperative Task Allocation for IIoT Networks

abstract: Current IoT services include industry-oriented services, which often require objects to run more than one task. However, the exponential growth of objects in IoT poses the challenge of distributing and managing task allocation among objects. One of the main goals of task allocation is to improve the quality of information and maximize the tasks to be performed. Although there are approaches that optimize and manage the dynamics of nodes, not all consider the quality of information and the distributed allocation over the cluster service. This paper proposes the mechanism CONTASKI for task allocation in IIoT networks to distribute tasks among objects. It relies on collaborative consensus to allocate tasks and similarity capabilities to know which objects can play in accomplishing those tasks. CONTASKI was evaluated on NS-3 and achieved 100% of allocated tasks incases with 75 and 100 nodes, and, on average, more than 80% clusters performed tasks in a low response time.


title: Interval parking functions

abstract: Interval parking functions (IPFs) are a generalization of ordinary parking functions in which each car is willing to park only in a fixed interval of spaces. Each interval parking function can be expressed as a pair $(a,b)$, where $a$ is a parking function and $b$ is a dual parking function. We say that a pair of permutations $(x,y)$ is \emph{reachable} if there is an IPF $(a,b)$ such that $x,y$ are the outcomes of $a,b$, respectively, as parking functions. Reachability is reflexive and antisymmetric, but not in general transitive. We prove that its transitive closure, the \emph{pseudoreachability order}, is precisely the bubble-sort order on the symmetric group $\Sym_n$, which can be expressed in terms of the normal form of a permutation in the sense of du~Cloux; in particular, it is isomorphic to the product of chains of lengths $2,\dots,n$. It is thus seen to be a special case of Armstrong's sorting order, which lies between the Bruhat and (left) weak orders.


title: Domain Adaptation with Morphologic Segmentation

abstract: We present a novel domain adaptation framework that uses morphologic segmentation to translate images from arbitrary input domains (real and synthetic) into a uniform output domain. Our framework is based on an established image-to-image translation pipeline that allows us to first transform the input image into a generalized representation that encodes morphology and semantics - the edge-plus-segmentation map (EPS) - which is then transformed into an output domain. Images transformed into the output domain are photo-realistic and free of artifacts that are commonly present across different real (e.g. lens flare, motion blur, etc.) and synthetic (e.g. unrealistic textures, simplified geometry, etc.) data sets. Our goal is to establish a preprocessing step that unifies data from multiple sources into a common representation that facilitates training downstream tasks in computer vision. This way, neural networks for existing tasks can be trained on a larger variety of training data, while they are also less affected by overfitting to specific data sets. We showcase the effectiveness of our approach by qualitatively and quantitatively evaluating our method on four data sets of simulated and real data of urban scenes. Additional results can be found on the project website available at http://jonathank.de/research/eps/ .


title: Evidence for a Rotational Component in the Circumgalactic Medium of Nearby Galaxies

abstract: We present results of a study comparing the relative velocity of $\rm Ly\alpha$ absorbers to the rotation velocity of nearby galaxy disks in the local universe ($z \leq 0.03$). We have obtained rotation curves via long-slit spectroscopy of eight galaxies with the Southern African Large Telescope, and combine this dataset with an additional 16 galaxies with data from the literature. Each galaxy appears within $3R_{\rm vir}$ of a QSO sightline with archival Cosmic Origin Spectrograph (COS) spectra. We study the velocity orientation of absorbers with respect to nearby galaxy's rotation, and compare with results from both the Steidel et al. (2002) monolithic halo model and a new cylindrical Navarro-Frenk-White galaxy halo model to interpret these data in the context of probing 3D galaxy halos via 1D QSO absorption-line spectroscopy. Relative to these models we find that up to $59\pm5\%$ of $\rm Ly\alpha$ absorbers have velocities consistent with co-rotation. We find the $\rm Ly\alpha$ co-rotation fraction to decrease with galaxy luminosity ($L^*$) and impact parameter in a model-independent fashion. We report that both anti-rotating absorbers and those found near luminous galaxies ($L \gtrsim 0.5 L^*$) mostly have low Doppler $b$-parameters ($b \lesssim 50$ km $\rm s^{-1}$). Absorbers consistent with co-rotation show a wide range of Doppler $b$-parameters. Finally, we find a strong anticorrelation between co-rotation fraction and galaxy inclination, which is at odds with recent metal-line kinematic studies and suggests the kinematic and geometric distribution of the circumgalactic medium is complex and multiphase.


title: The Sample Complexity of Teaching-by-Reinforcement on Q-Learning

abstract: We study the sample complexity of teaching, termed as "teaching dimension" (TDim) in the literature, for the teaching-by-reinforcement paradigm, where the teacher guides the student through rewards. This is distinct from the teaching-by-demonstration paradigm motivated by robotics applications, where the teacher teaches by providing demonstrations of state/action trajectories. The teaching-by-reinforcement paradigm applies to a wider range of real-world settings where a demonstration is inconvenient, but has not been studied systematically. In this paper, we focus on a specific family of reinforcement learning algorithms, Q-learning, and characterize the TDim under different teachers with varying control power over the environment, and present matching optimal teaching algorithms. Our TDim results provide the minimum number of samples needed for reinforcement learning, and we discuss their connections to standard PAC-style RL sample complexity and teaching-by-demonstration sample complexity results. Our teaching algorithms have the potential to speed up RL agent learning in applications where a helpful teacher is available.


title: Black hole collapse and bounce in effective loop quantum gravity

abstract: We derive effective equations with loop quantum gravity corrections for the Lema\^itre-Tolman-Bondi family of space-times, and use these to study quantum gravity effects in the Oppenheimer-Snyder collapse model. For this model, after the formation of a black hole with an apparent horizon, quantum gravity effects become important in the space-time region where the energy density and space-time curvature scalars become comparable to the Planck scale. These quantum gravity effects first stop the collapse of the dust matter field when its energy density reaches the Planck scale, and then cause the dust field to begin slowly expanding. Due to this continued expansion, the matter field will eventually extend beyond the apparent horizon, at which point the horizon disappears and there is no longer a black hole. There are no singularities anywhere in this space-time. In addition, in the limit that edge effects are neglected, we show that the dynamics for the interior of the star of uniform energy density follow the loop quantum cosmology effective Friedman equation for the spatially flat Friedman-Lema\^itre-Robertson-Walker space-time. Finally, we estimate the lifetime of the black hole, as measured by a distant observer, to be $\sim (GM)^2/\ell_{\rm Pl}$.


title: Digital Contact Tracing Using IP Colocation

abstract: The spread of an infectious disease through a population can be modeled using a network or a graph. In digital advertising, internet device graphs are graph data sets that organize identifiers produced by mobile phones, PCs, TVs, and tablets as they access media on the internet. Characterized by immense scale, they have become ubiquitous as they enable targeted advertising, content customization and tracking. This paper posits that internet device graphs, in particular those based on IP colocation, can provide significant utility in predicting and modeling the spread of infectious disease. Starting the week of March 16th, 2020, in the United States, many individuals began to `shelter-in-place' as schools and workplaces across the nation closed because of the COVID-19 pandemic. This paper quantifies the effect of the shelter-in-place orders on a large scale internet device graph with more than a billion nodes by studying the graph before and after orders went into effect. The effects are clearly visible. The structure of the graph suggests behavior least conducive to transmission of infection occurred in the US between April 12th and 19th, 2020. This paper also discusses the utility of device graphs for i) contact tracing, ii) prediction of `hot spots', iii) simulation of infectious disease spread, and iv) delivery of advertisement-based warnings to potentially exposed individuals. The paper also posits an overarching question: can systems and datasets amassed by entities in the digital ad ecosystem aid in the fight against COVID-19?


title: Submodular Maximization in Clean Linear Time

abstract: In this paper, we provide the first deterministic algorithm that achieves the tight $1-1/e$ approximation guarantee for submodular maximization under a cardinality (size) constraint while making a number of queries that scales only linearly with the size of the ground set $n$. To complement our result, we also show strong information-theoretic lower bounds. More specifically, we show that when the maximum cardinality allowed for a solution is constant, no algorithm making a sub-linear number of function evaluations can guarantee any constant approximation ratio. Furthermore, when the constraint allows the selection of a constant fraction of the ground set, we show that any algorithm making fewer than $\Omega(n/\log(n))$ function evaluations cannot perform better than an algorithm that simply outputs a uniformly random subset of the ground set of the right size. We then provide a variant of our deterministic algorithm for the more general knapsack constraint, which is the first linear-time algorithm that achieves $1/2$-approximation guarantee for this constraint. Finally, we extend our results to the general case of maximizing a monotone submodular function subject to the intersection of a $p$-set system and multiple knapsack constraints. We extensively evaluate the performance of our algorithms on multiple real-life machine learning applications, including movie recommendation, location summarization, twitter text summarization and video summarization.


title: Aerial Platforms with Reconfigurable Smart Surfaces for 5G and Beyond

abstract: Aerial platforms are expected to deliver enhanced and seamless connectivity in the fifth generation (5G) wireless networks and beyond (B5G). This is generally achievable by supporting advanced onboard communication features embedded in heavy and energy-intensive equipment. Alternatively, reconfigurable smart surfaces (RSS), which smartly exploit/recycle signal reflections in the environment, are increasingly being recognized as a new wireless communication paradigm to improve communication links. In fact, their reduced cost, low power use, light weight, and flexible deployment make them an attractive candidate for integration with 5G/B5G technologies. In this article, we discuss comprehensive approaches to the integration of RSS in aerial platforms. First, we present a review of RSS technology, its operations and types of communication. Next, we describe how RSS can be used in aerial platforms, and we propose a control architecture workflow. Then, several potential use cases are presented and discussed. Finally, associated research challenges are identified.


title: Improving Piecewise Linear Snow Density Models through Hierarchical Spatial and Orthogonal Functional Smoothing

abstract: Snow density estimates as a function of depth are used for understanding climate processes, evaluating water accumulation trends in polar regions, and estimating glacier mass balances. The common and interpretable physically-derived differential equation models for snow density are piecewise linear as a function of depth (on a transformed scale); thus, they can fail to capture important data features. Moreover, the differential equation parameters show strong spatial autocorrelation. To address these issues, we allow the parameters of the physical model, including random change points over depth, to vary spatially. We also develop a framework for functionally smoothing the physically-motivated model. To preserve inference on the interpretable physical model, we project the smoothing function into the physical model's spatially varying null space. The proposed spatially and functionally smoothed snow density model better fits the data while preserving inference on physical parameters. Using this model, we find significant spatial variation in the parameters that govern snow densification.


title: Characterising continuum variability in the radio-loud narrow-line Seyfert 1 galaxy IRAS 17020+4544

abstract: We present results of temporal and spectral analyses on four XMM-Newton EPIC pn observations of IRAS 17020+4544, a narrow-line Seyfert 1 galaxy with evidence of a radio jet. Analysis of the light curves reveals that this radio-loud source does not behave like the bulk population of its radio-quiet counterparts. A trend of spectral hardening with increased flux is found. Variability is found to increase with energy, though it decreases as the spectrum hardens. The first 40 ks of the most recent observation behave uniquely among the epochs, exhibiting a softer spectral state than at any other time. Possible non-stationarity at low energies is found, with no such effect present at higher energies, suggesting at least two distinct spectral components. A reverberation signature is confirmed, with the lag-frequency, lag-energy, and covariance spectra changing significantly during the soft-state epoch. The temporal analysis suggests a variable power-law in the presence of a reflection component, thus motivating such a fit for the 0.3-10 keV EPIC pn spectra from all epochs. We find an acceptable spectral fit using the timing-motivated parameters and report the detection of a broad Fe K emission line, requiring an additional model component beyond the reflection spectrum. We discuss links between this source and other narrow-line Seyfert 1 sources that show evidence of jet activity, finding similarities among this currently very limited sample of interesting objects.


title: Drag enhancement in a dusty Kolmogorov flow

abstract: Particles suspended in a fluid exert feedback forces that can significantly impact the flow, altering the turbulent drag and velocity fluctuations. We study flow modulation induced by particles heavier than the carrier fluid in the framework of an Eulerian two-way coupled model, where particles are represented by a continuum density transported by a compressible velocity field, exchanging momentum with the fluid phase. We implement the model in direct numerical simulations of the turbulent Kolmogorov flow, a simplified setting allowing for studying the momentum balance and the turbulent drag in the absence of boundaries. We show that the amplitude of the mean flow and the turbulence intensity are reduced by increasing particle mass loading with the consequent enhancement of the friction coefficient. Surprisingly, turbulence suppression is stronger for particles of smaller inertia. We understand such a result by mapping the equations for dusty flow, in the limit of vanishing inertia, to a Newtonian flow with an effective forcing reduced by the increase in fluid density due to the presence of particles. We also discuss the negative feedback produced by turbophoresis which mitigates the effects of particles, especially with larger inertia, on the turbulent flow.


title: What's in the Image? Explorable Decoding of Compressed Images

abstract: The ever-growing amounts of visual contents captured on a daily basis necessitate the use of lossy compression methods in order to save storage space and transmission bandwidth. While extensive research efforts are devoted to improving compression techniques, every method inevitably discards information. Especially at low bit rates, this information often corresponds to semantically meaningful visual cues, so that decompression involves significant ambiguity. In spite of this fact, existing decompression algorithms typically produce only a single output, and do not allow the viewer to explore the set of images that map to the given compressed code. In this work we propose the first image decompression method to facilitate user-exploration of the diverse set of natural images that could have given rise to the compressed input code, thus granting users the ability to determine what could and what could not have been there in the original scene. Specifically, we develop a novel deep-network based decoder architecture for the ubiquitous JPEG standard, which allows traversing the set of decompressed images that are consistent with the compressed JPEG file. To allow for simple user interaction, we develop a graphical user interface comprising several intuitive exploration tools, including an automatic tool for examining specific solutions of interest. We exemplify our framework on graphical, medical and forensic use cases, demonstrating its wide range of potential applications.


title: High order methods for acoustic scattering: Coupling Farfield Expansions ABC with Deferred-Correction methods

abstract: Arbitrary high order numerical methods for time-harmonic acoustic scattering problems originally defined on unbounded domains are constructed. This is done by coupling recently developed high order local absorbing boundary conditions (ABCs) with finite difference methods for the Helmholtz equation. These ABCs are based on exact representations of the outgoing waves by means of farfield expansions. The finite difference methods, which are constructed from a deferred-correction (DC) technique, approximate the Helmholtz equation and the ABCs, with the appropriate number of terms, to any desired order. As a result, high order numerical methods with an overall order of convergence equal to the order of the DC schemes are obtained. A detailed construction of these DC finite difference schemes is presented. Additionally, a rigorous proof of the consistency of the DC schemes with the Helmholtz equation and the ABCs in polar coordinates is also given. The results of several numerical experiments corroborate the high order convergence of the novel method.


title: Gravity Probe Spin: Prospects for measuring general-relativistic precession of intrinsic spin using a ferromagnetic gyroscope

abstract: An experimental test at the intersection of quantum physics and general relativity is proposed: measurement of relativistic frame dragging and geodetic precession using intrinsic spin of electrons. The behavior of intrinsic spin in spacetime dragged and warped by a massive rotating body is an experimentally open question, hence the results of such a measurement could have important theoretical consequences. Such a measurement is possible by using mm-scale ferromagnetic gyroscopes in orbit around the Earth. Under conditions where the rotational angular momentum of a ferromagnet is sufficiently small, a ferromagnet's angular momentum is dominated by atomic electron spins and is predicted to exhibit macroscopic gyroscopic behavior. If such a ferromagnetic gyroscope is sufficiently isolated from the environment, rapid averaging of quantum uncertainty via the spin-lattice interaction enables readout of the ferromagnetic gyroscope dynamics with sufficient sensitivity to measure both the Lense-Thirring (frame dragging) and de Sitter (geodetic precession) effects due to the Earth.


title: Two-photon interference: the Hong-Ou-Mandel effect

abstract: Nearly 30 years ago, two-photon interference was observed, marking the beginning of a new quantum era. Indeed, two-photon interference has no classical analogue, giving it a distinct advantage for a range of applications. The peculiarities of quantum physics may now be used to our advantage to outperform classical computations, securely communicate information, simulate highly complex physical systems and increase the sensitivity of precise measurements. This separation from classical to quantum physics has motivated physicists to study two-particle interference for both fermionic and bosonic quantum objects. So far, two-particle interference has been observed with massive particles, among others, such as electrons and atoms, in addition to plasmons, demonstrating the extent of this effect to larger and more complex quantum systems. A wide array of novel applications to this quantum effect is to be expected in the future. This review will thus cover the progress and applications of two-photon (two-particle) interference over the last three decades.


title: Cross-Cultural Similarity Features for Cross-Lingual Transfer Learning of Pragmatically Motivated Tasks

abstract: Much work in cross-lingual transfer learning explored how to select better transfer languages for multilingual tasks, primarily focusing on typological and genealogical similarities between languages. We hypothesize that these measures of linguistic proximity are not enough when working with pragmatically-motivated tasks, such as sentiment analysis. As an alternative, we introduce three linguistic features that capture cross-cultural similarities that manifest in linguistic patterns and quantify distinct aspects of language pragmatics: language context-level, figurative language, and the lexification of emotion concepts. Our analyses show that the proposed pragmatic features do capture cross-cultural similarities and align well with existing work in sociolinguistics and linguistic anthropology. We further corroborate the effectiveness of pragmatically-driven transfer in the downstream task of choosing transfer languages for cross-lingual sentiment analysis.


title: Fast Free-text Authentication via Instance-based Keystroke Dynamics

abstract: Keystroke dynamics study the way in which users input text via their keyboards. Having the ability to differentiate users, typing behaviors can unobtrusively form a component of a behavioral biometric recognition system to improve existing account security. Keystroke dynamics systems on free-text data have previously required 500 or more characters to achieve reasonable performance. In this paper, we propose a novel instance-based graph comparison algorithm called the instance-based tail area density (ITAD) metric to reduce the number of keystrokes required to authenticate users. Additionally, commonly used features in the keystroke dynamics literature, such as monographs and digraphs, are all found to be useful in informing who is typing. The usefulness of these features for authentication is determined using a random forest classifier and validated across two publicly available datasets. Scores from the individual features are fused to form a single matching score. With the fused matching score and our ITAD metric, we achieve equal error rates (EERs) for 100 and 200 testing digraphs of 9.7% and 7.8% for the Clarkson II dataset, improving upon state-of-the-art of 35.3% and 15.3%.


title: The Z\gamma{} transverse-momentum spectrum at NNLO+N3LL

abstract: We consider the transverse-momentum ($p_T$) distribution of $Z\gamma$ pairs produced in hadronic collisions. Logarithmically enhanced contributions at small $p_T$ are resummed to all orders in QCD perturbation theory and combined with the fixed-order prediction. We achieve the most advanced prediction for the $Z\gamma$ $p_T$ spectrum by matching next-to-next-to-next-to-leading logarithmic (N$^3$LL) resummation to the integrated cross section at next-to-next-to-leading order (NNLO). By considering $\ell^+\ell^-\gamma$ production at the fully differential level, including spin correlations, interferences and off-shell effects, arbitrary cuts can be applied to the leptons and the photon. We present results at the LHC in presence of fiducial cuts and find agreement with the $13$\,TeV ATLAS data at the few-percent level.


title: Elastodynamic single-sided homogeneous Green's function representation: Theory and numerical examples

abstract: The homogeneous Green's function is the difference between an impulse response and its time-reversal. According to existing representation theorems, the homogeneous Green's function associated with source-receiver pairs inside a medium can be computed from measurements at a boundary enclosing the medium. However, in many applications such as seismic imaging, time-lapse monitoring, medical imaging, non-destructive testing, etc., media are only accessible from one side. A recent development of wave theory has provided a representation of the homogeneous Green's function in an elastic medium in terms of wavefield recordings at a single (open) boundary. Despite its single-sidedness, the elastodynamic homogeneous Green's function representation accounts for all orders of scattering inside the medium. We present the theory of the elastodynamic single-sided homogeneous Green's function representation and illustrate it with numerical examples for 2D laterally-invariant media. For propagating waves, the resulting homogeneous Green's functions match the exact ones within numerical precision, demonstrating the accuracy of the theory. In addition, we analyse the accuracy of the single-sided representation of the homogeneous Green's function for evanescent wave tunnelling.


title: Quasiregular curves: H\"older continuity and higher integrability

abstract: We show that a $K$-quasiregular $\omega$-curve from a Euclidean domain to a Euclidean space with respect to a covector $\omega$ is locally $(1/K)(\lVert \omega\rVert/|\omega|_{\ell_1})$-H\"older continuous. We also show that quasiregular curves enjoy higher integrability.


title: Photoacoustic fluctuation imaging: theory and application to blood flow imaging

abstract: Photoacoustic fluctuation imaging, which exploits randomness in photoacoustic generation, provides enhanced images in terms of resolution and visibility, as compared to conventional photoacoustic images. While a few experimental demonstrations of photoacoustic fluctuation imaging have been reported, it has to date not been described theoretically. In the first part of this work, we propose a theory relevant to fluctuations induced either by random illumination patterns or by random distributions of absorbing particles. The theoretical predictions are validated by Monte Carlo finite-difference time-domain simulations of photoacoustic generation in random particle media. We provide a physical insight into why visibility artefacts are absent from second-order fluctuation images. In the second part, we demonstrate experimentally that harnessing randomness induced by the flow of red blood cells produce photoacoustic fluctuation images free of visibility artefacts. As a first proof of concept, we obtain two-dimensional images of blood vessel phantoms. Photoacoustic fluctuation imaging is finally applied in vivo to obtain 3D images of the vascularization in a chicken embryo.


title: Play with One's Feelings: A Study on Emotion Awareness for Player Experience

abstract: Affective interaction between players of video games can elicit rich and varying patterns of emotions. In multiplayer activities that take place in a common space (such as sports and board games), players are generally aware of the emotions of their teammates or opponents as they can directly observe their behavioral patterns, facial expressions, head pose, body stance and so on. Players of online video games, however, are not generally aware of the other players' emotions given the limited channels of direct interaction among them (e.g. via emojis or chat boxes). It also turns out that the impact of real-time emotionawareness on play is still unexplored in the space of online digital games. Motivated by this lack of empirical knowledge on the role of the affect of others to one's gameplay performance in this paper we investigate the degrees to which the expression of manifested emotions of an opponent can affect the emotions of the player and consequently his gameplay behavior. In this initial study, we test our hypothesis on a two-player adversarial car racing game. We perform a comprehensive user study to evaluate the emotions, behaviors, and attitudes of players in emotion aware versus emotion agnostic game versions. Our findings suggest that expressing the emotional state of the opponent through an emoji in real-time affects the emotional state and behavior of players that can consequently affect their playing experience.


title: Contamination of TEM Holders Quantified and Mitigated with Open-Hardware, High-Vacuum Bakeout System

abstract: Hydrocarbon contamination plagues high-resolution and analytical electron microscopy by depositing carbonaceous layers onto surfaces during electron irradiation, which can render carefully prepared specimens useless. Increased specimen thickness degrades resolution with beam broadening alongside loss of contrast. The large inelastic cross-section of carbon hampers accurate atomic species detection. Oxygen and water molecules pose problems of lattice damage by chemically etching the specimen during imaging. These constraints on high-resolution and spectroscopic imaging demand clean, high-vacuum microscopes with dry pumps. Here, we present an open-hardware design of a high-vacuum manifold for transmission electron microscopy (TEM) holders to mitigate hydrocarbon and residual species exposure. We quantitatively show that TEM holders are inherently dirty and introduce a range of unwanted chemical species. Overnight storage in our manifold reduces contaminants by 1-2 orders of magnitude and promotes 2-4 times faster vacuum recovery. A built-in bakeout system further reduces contaminants partial pressure to below $10^{-10}$ Torr (~4 orders of magnitude down from ambient storage) and alleviates monolayer adsorption during a typical TEM experiment. We determine that bakeout of TEM holder with specimen held therein is the optimal cleaning method. Our high-vacuum manifold design is published with open-source blueprints, parts list, and cost.


title: Falconer-type estimates for dot products

abstract: We present a family of sharpness examples for Falconer-type single dot product results. In particular, for $d\geq 2,$ for any $s<\frac{d+1}{2},$ we construct a Borel probability measure $\mu$ satisfying the energy estimate $I_s(\mu)<\infty,$ yet the estimate \begin{equation} (\mu \times \mu)\{(x,y):1\leq x\cdot y \leq 1+\epsilon\} \leq C\epsilon \end{equation} does not hold with constants independent of $\epsilon$. It is known (\cite{EIT11}) that such an estimate always holds with $C$ independent of $\epsilon$ if $I_{\frac{d+1}{2}}(\mu)<\infty$. Thus our estimate proves the sharpness of the dimensional threshold in this result and generalizes similar results (\cite{Mat95}, \cite{IS16}) established in the case when the dot product $x \cdot y$ is replaced by the Euclidean distance function $|x-y|$, or, more generally, ${||x-y||}_K$, the distance that comes from the norm induced by a symmetric convex body $K$ with a smooth boundary and non-vanishing curvature. Our constructions are partially based on ideas that come from discrete incidence theory.


title: On the existence of global smooth solutions to the parabolic-elliptic Keller-Segel system with irregular initial data

abstract: We consider the parabolic-elliptic Keller-Segel system \[ \left\{ \begin{aligned} u_t &= \Delta u - \chi \nabla \cdot (u \nabla v), \\ 0 &= \Delta v - v + u \end{aligned} \right. \tag{$\star$} \] in a smooth bounded domain $\Omega \subseteq \mathbb{R}^n$, $n\in\mathbb{N}$, with Neumann boundary conditions. We look at both chemotactic attraction ($\chi > 0$) and repulsion ($\chi < 0$) scenarios in two and three dimensions. The key feature of interest for the purposes of this paper is under which conditions said system still admits global classical solutions due to the smoothing properties of the Laplacian even if the initial data is very irregular. Regarding this, we show for initial data $\mu \in \mathcal{M}_+(\overline{\Omega})$ that, if either \begin{align} \bullet\;& n = 2,\, \chi < 0\, \text{ or } \\ \bullet\;& n = 2,\, \chi > 0\, \text{ and the initial mass is small or } \\ \bullet\;& n = 3,\, \chi < 0\, \text{ and } \mu = f \in L^p(\Omega), p > 1 \end{align} holds, it is still possible to construct global classical solutions to ($\star$), which are continuous in $t = 0$ in the vague topology on $M_+(\overline{\Omega})$.


title: The Topological Symmetric Orbifold

abstract: We analyse topological orbifold conformal field theories on the symmetric product of a complex surface M. By exploiting the mathematics literature we show that a canonical quotient of the operator ring has structure constants given by Hurwitz numbers. This proves a conjecture in the physics literature on extremal correlators. Moreover, it allows to leverage results on the combinatorics of the symmetric group to compute more structure constants explicitly. We recall that the full orbifold chiral ring is given by a symmetric orbifold Frobenius algebra. This construction enables the computation of topological genus zero and genus one correlators, and to prove the vanishing of higher genus contributions. The efficient description of all topological correlators sets the stage for a proof of a topological AdS/CFT correspondence. Indeed, we propose a concrete mathematical incarnation of the proof, relating Gromow-Witten theory in the bulk to the quantum cohomology of the Hilbert scheme on the boundary.


title: Understanding and Mitigating Exploding Inverses in Invertible Neural Networks

abstract: Invertible neural networks (INNs) have been used to design generative models, implement memory-saving gradient computation, and solve inverse problems. In this work, we show that commonly-used INN architectures suffer from exploding inverses and are thus prone to becoming numerically non-invertible. Across a wide range of INN use-cases, we reveal failures including the non-applicability of the change-of-variables formula on in- and out-of-distribution (OOD) data, incorrect gradients for memory-saving backprop, and the inability to sample from normalizing flow models. We further derive bi-Lipschitz properties of atomic building blocks of common architectures. These insights into the stability of INNs then provide ways forward to remedy these failures. For tasks where local invertibility is sufficient, like memory-saving backprop, we propose a flexible and efficient regularizer. For problems where global invertibility is necessary, such as applying normalizing flows on OOD data, we show the importance of designing stable INN building blocks.


title: LiDARsim: Realistic LiDAR Simulation by Leveraging the Real World

abstract: We tackle the problem of producing realistic simulations of LiDAR point clouds, the sensor of preference for most self-driving vehicles. We argue that, by leveraging real data, we can simulate the complex world more realistically compared to employing virtual worlds built from CAD/procedural models. Towards this goal, we first build a large catalog of 3D static maps and 3D dynamic objects by driving around several cities with our self-driving fleet. We can then generate scenarios by selecting a scene from our catalog and "virtually" placing the self-driving vehicle (SDV) and a set of dynamic objects from the catalog in plausible locations in the scene. To produce realistic simulations, we develop a novel simulator that captures both the power of physics-based and learning-based simulation. We first utilize ray casting over the 3D scene and then use a deep neural network to produce deviations from the physics-based simulation, producing realistic LiDAR point clouds. We showcase LiDARsim's usefulness for perception algorithms-testing on long-tail events and end-to-end closed-loop evaluation on safety-critical scenarios.


title: Foundations for Bayesian inference with engineered likelihood functions for robust amplitude estimation

abstract: We present mathematical and conceptual foundations for the task of robust amplitude estimation using engineered likelihood functions (ELFs), a framework introduced in Wang et al. [PRX Quantum 2, 010346 (2021)] that uses Bayesian inference to enhance the rate of information gain in quantum sampling. These ELFs, which are obtained by choosing tunable parameters in a parametrized quantum circuit to minimize the expected posterior variance of an estimated parameter, play an important role in estimating the expectation values of quantum observables. We give a thorough characterization and analysis of likelihood functions arising from certain classes of quantum circuits and combine this with the tools of Bayesian inference to give a procedure for picking optimal ELF tunable parameters. Finally, we present numerical results to demonstrate the performance of ELFs.


title: Minimizing estimation runtime on noisy quantum computers

abstract: The number of measurements demanded by hybrid quantum-classical algorithms such as the variational quantum eigensolver (VQE) is prohibitively high for many problems of practical value. For such problems, realizing quantum advantage will require methods which dramatically reduce this cost. Previous quantum algorithms that reduce the measurement cost (e.g. quantum amplitude and phase estimation) require error rates that are too low for near-term implementation. Here we propose methods that take advantage of the available quantum coherence to maximally enhance the power of sampling on noisy quantum devices, reducing measurement number and runtime compared to the standard sampling method of the variational quantum eigensolver (VQE). Our scheme derives inspiration from quantum metrology, phase estimation, and the more recent "alpha-VQE" proposal, arriving at a general formulation that is robust to error and does not require ancilla qubits. The central object of this method is what we call the "engineered likelihood function" (ELF), used for carrying out Bayesian inference. We show how the ELF formalism enhances the rate of information gain in sampling as the physical hardware transitions from the regime of noisy intermediate-scale quantum computers into that of quantum error corrected ones. This technique speeds up a central component of many quantum algorithms, with applications including chemistry, materials, finance, and beyond. Similar to VQE, we expect small-scale implementations to be realizable on today's quantum devices.


title: Helical Edge States and Quantum Phase Transitions in Tetralayer Graphene

abstract: Helical conductors with spin-momentum locking are promising platforms for Majorana fermions. Here we report observation of two topologically distinct phases supporting helical edge states in charge neutral Bernal-stacked tetralayer graphene in Hall bar and Corbino geometries. As the magnetic field B and out-of-plane displacement field D are varied, we observe a phase diagram consisting of an insulating phase and two metallic phases, with 0, 1 and 2 helical edge states, respectively. These phases are accounted for by a theoretical model that relates their conductance to spin-polarization plateaus. Transitions between them arise from a competition among inter-layer hopping, electrostatic and exchange interaction energies. Our work highlights the complex competing symmetries and the rich quantum phases in few-layer graphene.


title: A One-Pass Private Sketch for Most Machine Learning Tasks

abstract: Differential privacy (DP) is a compelling privacy definition that explains the privacy-utility tradeoff via formal, provable guarantees. Inspired by recent progress toward general-purpose data release algorithms, we propose a private sketch, or small summary of the dataset, that supports a multitude of machine learning tasks including regression, classification, density estimation, near-neighbor search, and more. Our sketch consists of randomized contingency tables that are indexed with locality-sensitive hashing and constructed with an efficient one-pass algorithm. We prove competitive error bounds for DP kernel density estimation. Existing methods for DP kernel density estimation scale poorly, often exponentially slower with an increase in dimensions. In contrast, our sketch can quickly run on large, high-dimensional datasets in a single pass. Exhaustive experiments show that our generic sketch delivers a similar privacy-utility tradeoff when compared to existing DP methods at a fraction of the computation cost. We expect that our sketch will enable differential privacy in distributed, large-scale machine learning settings.


title: On Chow-weight homology of motivic complexes and its relation to motivic homology

abstract: We study in detail the so-called Chow-weight homology of Voevodsky motivic complexes and relate it to motivic homology. We generalize earlier results and prove that the vanishing of higher motivic homology groups of a motif $M$ implies similar vanishing for its Chow-weight homology along with effectivity properties of the higher terms of its weight complex $t(M)$ and of higher Deligne weight quotients of its cohomology. Applying this statement to motives with compact support we obtain a similar relation between the vanishing of Chow groups and the cohomology with compact support of varieties. Moreover, we prove that if higher motivic homology groups of a geometric motif or a variety over a universal domain are torsion (in a certain "range") then the exponents of these groups are uniformly bounded. To prove our main results we study Voevodsky slices of motives. Since the slice functors do not respect the compactness of motives, the results of the previous Chow-weight homology paper are not sufficient for our purposes; this is our main reason to extend them to ($w_{Chow}$-bounded below) motivic complexes.


title: A Cochain Level Proof of Adem Relations in the Mod 2 Steenrod Algebra

abstract: In 1947, N.E. Steenrod defined the Steenrod Squares, which are mod 2 cohomology operations, using explicit cochain formulae for cup-i products of cocycles. He later recast the construction in more general homological terms, using group homology and acyclic model methods, rather than explicit cochain formulae, to define mod p operations for all primes p. Steenrod's student J. Adem applied the homological point of view to prove fundamental relations, known as the Adem relations, in the algebra of cohomology operations generated by the Steenrod operations. In this paper we give a proof of the mod 2 Adem relations at the cochain level. Specifically, given a mod 2 cocycle, we produce explicit cochain formulae whose coboundaries are the Adem relations among compositions of Steenrod Squares applied to the cocycle, using Steenrod's original cochain definition of the Square operations.


title: A Note on the Global Convergence of Multilayer Neural Networks in the Mean Field Regime

abstract: In a recent work, we introduced a rigorous framework to describe the mean field limit of the gradient-based learning dynamics of multilayer neural networks, based on the idea of a neuronal embedding. There we also proved a global convergence guarantee for three-layer (as well as two-layer) networks using this framework. In this companion note, we point out that the insights in our previous work can be readily extended to prove a global convergence guarantee for multilayer networks of any depths. Unlike our previous three-layer global convergence guarantee that assumes i.i.d. initializations, our present result applies to a type of correlated initialization. This initialization allows to, at any finite training time, propagate a certain universal approximation property through the depth of the neural network. To achieve this effect, we introduce a bidirectional diversity condition.


title: Spectral deconvolution of unitarily invariant matrix models

abstract: The present paper implements a complex analytic method to recover the spectrum of a matrix perturbed by either the addition or the multiplication of a random matrix noise, under the assumption that the distribution of the noise is unitarily invariant. This method, introduced by Arizmendi, Tarrago and Vargas in arXiv:1711.08871, is done in two steps : the first step consists in a fixed point method to compute the Stieltjes transform of the desired distribution in a certain domain, and the second step is a classical deconvolution by a Cauchy distribution, whose parameter depends on the intensity of the noise. We also provide explicit bounds for the mean squared error of the first step.


title: ATP synthase: evolution, energetics, and membrane interactions

abstract: The synthesis of ATP, life's 'universal energy currency', is the most prevalent chemical reaction in biological systems, and is responsible for fueling nearly all cellular processes, from nerve impulse propagation to DNA synthesis. ATP synthases, the family of enzymes that carry out this endless task, are nearly as ubiquitous as the energy-laden molecule they are responsible for making. The F-type ATP synthase (F-ATPase) is found in every domain of life, and is believed to predate the divergence of these lineages over 1.5 billion years ago. These enzymes have therefore facilitated the survival of organisms in a wide range of habitats, ranging from the deep-sea thermal vents to the human intestine. In this review, we present an overview of the current knowledge of the structure and function of F-type ATPases, highlighting several adaptations that have been characterized across taxa. We emphasize the importance of studying these features within the context of the enzyme's particular lipid environment: Just as the interactions between an organism and its physical environment shape its evolutionary trajectory, ATPases are impacted by the membranes within which they reside. We argue that a comprehensive understanding of the structure, function, and evolution of membrane proteins -- including ATP synthase -- requires such an integrative approach.


title: Directional Pruning of Deep Neural Networks

abstract: In the light of the fact that the stochastic gradient descent (SGD) often finds a flat minimum valley in the training loss, we propose a novel directional pruning method which searches for a sparse minimizer in or close to that flat region. The proposed pruning method does not require retraining or the expert knowledge on the sparsity level. To overcome the computational formidability of estimating the flat directions, we propose to use a carefully tuned $\ell_1$ proximal gradient algorithm which can provably achieve the directional pruning with a small learning rate after sufficient training. The empirical results demonstrate the promising results of our solution in highly sparse regime (92% sparsity) among many existing pruning methods on the ResNet50 with the ImageNet, while using only a slightly higher wall time and memory footprint than the SGD. Using the VGG16 and the wide ResNet 28x10 on the CIFAR-10 and CIFAR-100, we demonstrate that our solution reaches the same minima valley as the SGD, and the minima found by our solution and the SGD do not deviate in directions that impact the training loss. The code that reproduces the results of this paper is available at https://github.com/donlan2710/gRDA-Optimizer/tree/master/directional_pruning.


title: AWAC: Accelerating Online Reinforcement Learning with Offline Datasets

abstract: Reinforcement learning (RL) provides an appealing formalism for learning control policies from experience. However, the classic active formulation of RL necessitates a lengthy active exploration process for each behavior, making it difficult to apply in real-world settings such as robotic control. If we can instead allow RL algorithms to effectively use previously collected data to aid the online learning process, such applications could be made substantially more practical: the prior data would provide a starting point that mitigates challenges due to exploration and sample complexity, while the online training enables the agent to perfect the desired skill. Such prior data could either constitute expert demonstrations or sub-optimal prior data that illustrates potentially useful transitions. While a number of prior methods have either used optimal demonstrations to bootstrap RL, or have used sub-optimal data to train purely offline, it remains exceptionally difficult to train a policy with offline data and actually continue to improve it further with online RL. In this paper we analyze why this problem is so challenging, and propose an algorithm that combines sample efficient dynamic programming with maximum likelihood policy updates, providing a simple and effective framework that is able to leverage large amounts of offline data and then quickly perform online fine-tuning of RL policies. We show that our method, advantage weighted actor critic (AWAC), enables rapid learning of skills with a combination of prior demonstration data and online experience. We demonstrate these benefits on simulated and real-world robotics domains, including dexterous manipulation with a real multi-fingered hand, drawer opening with a robotic arm, and rotating a valve. Our results show that incorporating prior data can reduce the time required to learn a range of robotic skills to practical time-scales.


title: Convection and Spin-Up During Common Envelope Evolution: The Formation of Short-Period Double White Dwarfs

abstract: The formation channels and predicted populations of double-white dwarfs (DWDs) are important because a subset will evolve to be gravitational-wave sources and/or progenitors of Type Ia supernovae. Given the observed population of short-period DWDs, we calculate the outcomes of common envelope evolution when convective effects are included. For each observed white dwarf in a DWD system, we identify all progenitor stars with an equivalent proto-WD core mass from a comprehensive suite of stellar evolution models. With the second observed white dwarf as the companion, we calculate the conditions under which convection can accommodate the energy released as the orbit decays, including (if necessary), how much the envelope must spin-up during the common envelope phase. The predicted post-CE final separations closely track the observed DWD orbital parameter space, further strengthening the view that convection is a key ingredient in common envelope evolution.


title: Gradient Free Minimax Optimization: Variance Reduction and Faster Convergence

abstract: Many important machine learning applications amount to solving minimax optimization problems, and in many cases there is no access to the gradient information, but only the function values. In this paper, we focus on such a gradient-free setting, and consider the nonconvex-strongly-concave minimax stochastic optimization problem. In the literature, various zeroth-order (i.e., gradient-free) minimax methods have been proposed, but none of them achieve the potentially feasible computational complexity of $\mathcal{O}(\epsilon^{-3})$ suggested by the stochastic nonconvex minimization theorem. In this paper, we adopt the variance reduction technique to design a novel zeroth-order variance reduced gradient descent ascent (ZO-VRGDA) algorithm. We show that the ZO-VRGDA algorithm achieves the best known query complexity of $\mathcal{O}(\kappa(d_1 + d_2)\epsilon^{-3})$, which outperforms all previous complexity bound by orders of magnitude, where $d_1$ and $d_2$ denote the dimensions of the optimization variables and $\kappa$ denotes the condition number. In particular, with a new analysis technique that we develop, our result does not rely on a diminishing or accuracy-dependent stepsize usually required in the existing methods. To our best knowledge, this is the first study of zeroth-order minimax optimization with variance reduction. Experimental results on the black-box distributional robust optimization problem demonstrates the advantageous performance of our new algorithm.


title: Solving polynomials with ordinary differential equations

abstract: In this work we consider a given root of a family of n-degree polynomials as a one-variable function that depends only on the independent term. Then we prove that this function satisfies several ordinary differential equations (ODE). More concretely, it satisfies several simple separated variables ODE, a first order generalized Abel ODE of degree n-1 and an (n-1)-th order linear ODE. Although some of our results are not new, our approach is simple and self-contained. For n=2, 3 and 4 we recover, from these ODE, the classical formulas for solving these polynomials.


title: Building One-Shot Semi-supervised (BOSS) Learning up to Fully Supervised Performance

abstract: Reaching the performance of fully supervised learning with unlabeled data and only labeling one sample per class might be ideal for deep learning applications. We demonstrate for the first time the potential for building one-shot semi-supervised (BOSS) learning on Cifar-10 and SVHN up to attain test accuracies that are comparable to fully supervised learning. Our method combines class prototype refining, class balancing, and self-training. A good prototype choice is essential and we propose a technique for obtaining iconic examples. In addition, we demonstrate that class balancing methods substantially improve accuracy results in semi-supervised learning to levels that allow self-training to reach the level of fully supervised learning performance. Rigorous empirical evaluations provide evidence that labeling large datasets is not necessary for training deep neural networks. We made our code available at https://github.com/lnsmith54/BOSS to facilitate replication and for use with future real-world applications.


title: Enhancing nonlinear damping by parametric-direct internal resonance

abstract: Mechanical sources of nonlinear damping play a central role in modern physics, from solid-state physics to thermodynamics. The microscopic theory of mechanical dissipation [M. I . Dykman, M. A. Krivoglaz, Physica Status Solidi (b) 68, 111 (1975)] suggests that nonlinear damping of a resonant mode can be strongly enhanced when it is coupled to a vibration mode that is close to twice its resonance frequency. To date, no experimental evidence of this enhancement has been realized. In this letter, we experimentally show that nanoresonators driven into parametric-direct internal resonance provide supporting evidence for the microscopic theory of nonlinear dissipation. By regulating the drive level, we tune the parametric resonance of a graphene nanodrum over a range of 40-70 MHz to reach successive two-to-one internal resonances, leading to a nearly two-fold increase of the nonlinear damping. Our study opens up an exciting route towards utilizing modal interactions and parametric resonance to realize resonators with engineered nonlinear dissipation over wide frequency range.


title: Byzantine-Robust Learning on Heterogeneous Datasets via Bucketing

abstract: In Byzantine robust distributed or federated learning, a central server wants to train a machine learning model over data distributed across multiple workers. However, a fraction of these workers may deviate from the prescribed algorithm and send arbitrary messages. While this problem has received significant attention recently, most current defenses assume that the workers have identical data. For realistic cases when the data across workers are heterogeneous (non-iid), we design new attacks which circumvent current defenses, leading to significant loss of performance. We then propose a simple bucketing scheme that adapts existing robust algorithms to heterogeneous datasets at a negligible computational cost. We also theoretically and experimentally validate our approach, showing that combining bucketing with existing robust algorithms is effective against challenging attacks. Our work is the first to establish guaranteed convergence for the non-iid Byzantine robust problem under realistic assumptions.


title: A generalization of the Levinson theorem about the asymptotic value of the scattering phase shift

abstract: In quantum scattering theory, there exists a relationship between the difference in the scattering phase shifts at threshold and infinity and the number of bound states, which is established by the Levinson theorem. The presence of Castillejo, Dalitz and Dyson poles in the scattering amplitude, as well as Jaffe and Low primitives, corresponding to zeros of $D$ function on the unitary cut, modify the Levinson theorem. The asymptotic value of the scattering phase shift is shown to be determined by the number of bound states, the number of Castillejo, Dalitz and Dyson poles, and the number of primitives. Some consequences of the generalized theorem with respect to properties of nucleon-nucleon interactions are discussed.


title: Semantic Curiosity for Active Visual Learning

abstract: In this paper, we study the task of embodied interactive learning for object detection. Given a set of environments (and some labeling budget), our goal is to learn an object detector by having an agent select what data to obtain labels for. How should an exploration policy decide which trajectory should be labeled? One possibility is to use a trained object detector's failure cases as an external reward. However, this will require labeling millions of frames required for training RL policies, which is infeasible. Instead, we explore a self-supervised approach for training our exploration policy by introducing a notion of semantic curiosity. Our semantic curiosity policy is based on a simple observation -- the detection outputs should be consistent. Therefore, our semantic curiosity rewards trajectories with inconsistent labeling behavior and encourages the exploration policy to explore such areas. The exploration policy trained via semantic curiosity generalizes to novel scenes and helps train an object detector that outperforms baselines trained with other possible alternatives such as random exploration, prediction-error curiosity, and coverage-maximizing exploration.


title: Galaxy bias and primordial non-Gaussianity: insights from galaxy formation simulations with IllustrisTNG

abstract: We study the impact that large-scale perturbations of (i) the matter density and (ii) the primordial gravitational potential with local primordial non-Gaussianity (PNG) have on galaxy formation using the IllustrisTNG model. We focus on the linear galaxy bias $b_1$ and the coefficient $b_\phi$ of the scale-dependent bias induced by PNG, which describe the response of galaxy number counts to these two types of perturbations, respectively. We perform our study using separate universe simulations, in which the effect of the perturbations is mimicked by changes to the cosmological parameters: modified cosmic matter density for $b_1$ and modified amplitude $\mathcal{A}_s$ of the primordial scalar power spectrum for $b_\phi$. We find that the widely used universality relation $b_\phi = 2\delta_c(b_1 - 1)$ is a poor description of the bias of haloes and galaxies selected by stellar mass $M_*$, which is instead described better by $b_\phi(M_*) = 2\delta_c(b_1(M_*) - p)$ with $p \in [0.4, 0.7]$. This is explained by the different impact that matter overdensities and local PNG have on the median stellar-to-halo-mass relation. A simple model of this impact allows us to describe the stellar mass dependence of $b_1$ and $b_\phi$ fairly well. Our results also show a nontrivial relation between $b_1$ and $b_\phi$ for galaxies selected by color and black hole mass accretion rate. Our results provide refined priors on $b_\phi$ for local PNG constraints and forecasts using galaxy clustering. Given that the widely used universality relation underpredicts $b_\phi(M_*)$, existing analyses may underestimate the true constraining power on local PNG.


title: Constraining Magnetic Fields in the Circumgalactic Medium

abstract: We study the properties of magnetic fields in the circumgalactic medium (CGM) of $z<1$ galaxies by correlating Faraday rotation measures of $\sim 1,000$ high-redshift radio sources with the foreground galaxy number density estimated from the DESI Legacy Imaging Surveys. This method enables us to extract signals of rotation measures contributed by intervening gas around multiple galaxies. Our results show that there is no detectable correlation between the distribution of rotation measures and the number of foreground galaxies, contrary to several previous results. Utilizing the non-detection signals, we estimate $3\sigma$ upper limits to the rotation measures from the CGM of $\sim 20 \rm \ rad/m^{2}$ within~50 kpc and $\sim 10 \rm \ rad / m^{2}$ at separations of $100$ kpc. By adopting a column density distribution of ionized gas obtained from absorption line measurements, we further estimate the strengths of coherent magnetic fields parallel to the line of sight of $<\rm 2 \ \mu G$ in the CGM. We show that the estimated upper limits of rotation measures and magnetic field strengths are sufficient to constrain outputs of recent galaxy magneto-hydrodynamic simulations. Finally, we discuss possible causes for the inconsistency between our results and previous works.


title: Two regularized energy-preserving finite difference methods for the logarithmic Klein-Gordon equation

abstract: We present and analyze two regularized finite difference methods which preserve energy of the logarithmic Klein-Gordon equation (LogKGE). In order to avoid singularity caused by the logarithmic nonlinearity of the LogKGE, we propose a regularized logarithmic Klein-Gordon equation (RLogKGE) with a small regulation parameter $0<\varepsilon\ll1$ to approximate the LogKGE with the convergence order $O(\varepsilon)$. By adopting the energy method, the inverse inequality, and the cut-off technique of the nonlinearity to bound the numerical solution, the error bound $O(h^{2}+\frac{\tau^{2}}{\varepsilon^{2}})$ of the two schemes with the mesh size $h$, the time step $\tau$ and the parameter $\varepsilon$. Numerical results are reported to support our conclusions.


title: Generalizations of prime submodules over non-commutative rings

abstract: Throughout this paper, $R$ is an associative ring (not necessarily commutative) with identity and $M$ is a right $R$-module with unitary. In this paper, we introduce a new concept of $\phi$-prime submodule over an associative ring with identity. Thus we define the concept as following: Assume that $S(M)$ is the set of all submodules of $M$ and $\phi:S(M)\rightarrow S(M)\cup\{\emptyset\}$ is a function. For every $Y\in S(M)$ and ideal $I$ of $R,$ a proper submodule $X$ of $M$ is called $\phi$-prime, if $YI\subseteq X$ and $YI\nsubseteq\phi(X),$ then $Y\subseteq X$ or $I\subseteq(X:_{R}M)$. Then we examine the properties of $\phi$-prime submodules and characterize it when $M$ is a multiplication module.


title: The CARMENES search for exoplanets around M dwarfs. Variability of the He I line at 10830 \AA

abstract: The He I infrared (IR) triplet at 10830 \AA is known as an activity indicator in solar-type stars and has become a primary diagnostic in exoplanetary transmission spectroscopy. He I lines are a tracer of the stellar extreme-ultraviolet irradiation from the transition region and corona. We study the variability of the He I IR triplet lines in a spectral time series of 319 M~dwarf stars that was obtained with the CARMENES high-resolution optical and near-infrared spectrograph at Calar Alto. We detect He I IR line variability in 18% of our sample stars, all of which show H$\alpha$ in emission. Therefore, we find detectable He I IR variability in 78% of the sub-sample of stars with H$\alpha$ emission. Detectable variability is strongly concentrated in the latest spectral sub-types, where the He I IR lines during quiescence are typically weak. The fraction of stars with detectable He I IR variation remains lower than 10% for stars earlier than M3.0 V, while it exceeds 30% for the later spectral sub-types. Flares are accompanied by particularly pronounced line variations, including strongly broadened lines with red and blue asymmetries. However, we also find evidence for enhanced He I IR absorption, which is potentially associated with increased high-energy irradiation levels at flare onset. Generally, He I IR and H$\alpha$ line variations tend to be correlated, with H$\alpha$ being the most sensitive indicator in terms of pseudo-equivalent width variation. This makes the He I IR triplet a favourable target for planetary transmission spectroscopy.


title: The shape and simplicity biases of adversarially robust ImageNet-trained CNNs

abstract: Increasingly more similarities between human vision and convolutional neural networks (CNNs) have been revealed in the past few years. Yet, vanilla CNNs often fall short in generalizing to adversarial or out-of-distribution (OOD) examples which humans demonstrate superior performance. Adversarial training is a leading learning algorithm for improving the robustness of CNNs on adversarial and OOD data; however, little is known about the properties, specifically the shape bias and internal features learned inside adversarially-robust CNNs. In this paper, we perform a thorough, systematic study to understand the shape bias and some internal mechanisms that enable the generalizability of AlexNet, GoogLeNet, and ResNet-50 models trained via adversarial training. We find that while standard ImageNet classifiers have a strong texture bias, their R counterparts rely heavily on shapes. Remarkably, adversarial training induces three simplicity biases into hidden neurons in the process of "robustifying" CNNs. That is, each convolutional neuron in R networks often changes to detecting (1) pixel-wise smoother patterns, i.e., a mechanism that blocks high-frequency noise from passing through the network; (2) more lower-level features i.e. textures and colors (instead of objects);and (3) fewer types of inputs. Our findings reveal the interesting mechanisms that made networks more adversarially robust and also explain some recent findings e.g., why R networks benefit from a much larger capacity (Xie et al. 2020) and can act as a strong image prior in image synthesis (Santurkar et al. 2019).


title: Tree-level Interference in VBF production of Vh

abstract: Vector boson scattering is a well known probe of electroweak symmetry breaking. Here we study a related process of two electroweak vector bosons scattering into a vector boson and a Higgs boson ($VV \rightarrow Vh, V=W,Z$). This process exhibits tree level interference and grows with energy if the Higgs couplings to electroweak bosons deviate from their Standard Model values. Therefore, this process is particularly sensitive to the relative sign of the ratio of the coupling between the Higgs and the $W$ and $Z$, $\lambda_{WZ}$. In this work we show that a high energy lepton collider is well suited to study this process through vector boson fusion, estimate the potential sensitivity to this ratio, and show that a relatively modest amount of data can exclude $\lambda_{WZ} \simeq -1$.


title: Fallback Rates from Partial Tidal Disruption Events

abstract: A tidal disruption event (TDE) occurs when a star plunges through a supermassive black hole's tidal radius, at which point the star's self-gravity is overwhelmed by the tidal gravity of the black hole. In a partial TDE, where the star does not reach the full disruption radius, only a fraction of the star's mass is tidally stripped while the rest remains intact in the form of a surviving core. Analytical arguments have recently suggested that the temporal scaling of the fallback rate of debris to the black hole asymptotes to $t^{-9/4}$ for partial disruptions, effectively independently of the mass of the intact core. We present hydrodynamical simulations that verify the existence of this predicted, $t^{-9/4}$ scaling. We also define a break timescale -- the time at which the fallback rate transitions from a $t^{-5/3}$ scaling to the characteristic $t^{-9/4}$ scaling -- and measure this break timescale as a function of the impact parameter and the surviving core mass. These results deepen our understanding of the properties and breadth of possible fallback curves expected from TDEs and will therefore facilitate more accurate interpretation of data from wide-field surveys.


title: Visualization of two-dimensional transition dipole moment texture in momentum space using high-harmonic generation spectroscopy

abstract: Highly nonlinear optical phenomena can provide access to properties of electronic systems which are otherwise difficult to access through conventional linear optical spectroscopies. In particular, high harmonic generation (HHG) in crystalline solids is strikingly different from that in atomic gases, and it enables us to access electronic properties such as the band structure, Berry curvature, and valence electron density. Here, we show that polarization-resolved HHG measurements can be used to probe the transition dipole moment (TDM) texture in momentum space in two dimensional semiconductors. TDM is directly related to the internal structure of the electronic system and governs the optical properties. We study HHG in black phosphorus, which offers a simple two-band system, with bandgap resonant excitation. We observed a unique crystal-orientation dependence of the HHG yields and polarizations and succeeded in reconstructing the TDM texture related to the inter-atomic bonding structure. Our results demonstrate the potential of high harmonic spectroscopy for probing electronic wavefunctions in crystalline solids.


title: Lyman continuum leakage in faint star-forming galaxies at redshift z=3-3.5 probed by gamma-ray bursts

abstract: We present the observations of Lyman continuum (LyC) emission in the afterglow spectra of GRB 191004B at $z=3.5055$, together with those of the other two previously known LyC-emitting long gamma-ray bursts (LGRB) (GRB 050908 at $z=3.3467$, and GRB 060607A at $z=3.0749$), to determine their LyC escape fraction and compare their properties. From the afterglow spectrum of GRB 191004B we determine a neutral hydrogen column density at the LGRB redshift of $\log(N_{\rm HI}/cm^{-2})= 17.2 \pm 0.15$, and negligible extinction ($A_{\rm V}=0.03 \pm 0.02$ mag). The only metal absorption lines detected are CIV and SiIV. In contrast to GRB 050908 and GRB 060607A, the host galaxy of GRB 191004B displays significant Ly$\alpha$ emission. From its Ly$\alpha$ emission and the non-detection of Balmer emission lines we constrain its star-formation rate (SFR) to $1 \leq$ SFR $\leq 4.7$ M$_{\odot}\ yr^{-1}$. We fit the Ly$\alpha$ emission with a shell model and find parameters values consistent with the observed ones. The absolute LyC escape fractions we find for GRB 191004B, GRB 050908 and GRB 060607A are of $0.35^{+0.10}_{-0.11}$, $0.08^{+0.05}_{-0.04}$ and $0.20^{+0.05}_{-0.05}$, respectively. We compare the LyC escape fraction of LGRBs to the values of other LyC emitters found from the literature, showing that LGRB afterglows can be powerful tools to study LyC escape for faint high-redshift star-forming galaxies. Indeed we could push LyC leakage studies to much higher absolute magnitudes. The host galaxies of the three LGRB presented here have all $M_{\rm 1600} > -19.5$ mag, with the GRB 060607A host at $M_{\rm 1600} > -16$ mag. LGRB hosts may therefore be particularly suitable for exploring the ionizing escape fraction in galaxies that are too faint or distant for conventional techniques. Furthermore the time investment is very small compared to galaxy studies. [Abridged]


title: Black Holes, Moduli, and Long-Range Forces

abstract: It is well known that an identical pair of extremal Reissner-Nordstrom black holes placed a large distance apart will exert no force on each other. In this paper, I establish that the same result holds in a very large class of two-derivative effective theories containing an arbitrary number of gauge fields and moduli, where the appropriate analog of an extremal Reissner-Nordstrom black hole is a charged, spherically symmetric black hole with vanishing surface gravity or vanishing horizon area. Analogous results hold for black branes.


title: CosTuuM: polarized thermal dust emission by magnetically oriented spheroidal grains

abstract: We present the new open source C++-based Python library CosTuuM that can be used to generate infrared absorption and emission coefficients for arbitrary mixtures of spheroidal dust grains that are (partially) aligned with a magnetic field. We outline the algorithms underlying the software, demonstrate the accuracy of our results using benchmarks from literature, and use our tool to investigate some commonly used approximative recipes. We find that the linear polarization fraction for a partially aligned dust grain mixture can be accurately represented by an appropriate linear combination of perfectly aligned grains and grains that are randomly oriented, but that the commonly used picket fence alignment breaks down for short wavelengths. We also find that for a fixed dust grain size, the absorption coefficients and linear polarization fraction for a realistic mixture of grains with various shapes cannot both be accurately represented by a single representative grain with a fixed shape, but that instead an average over an appropriate shape distribution should be used. Insufficient knowledge of an appropriate shape distribution is the main obstacle in obtaining accurate optical properties. CosTuuM is available as a standalone Python library and can be used to generate optical properties to be used in radiative transfer applications.


title: Finding eV-scale Light Relics with Cosmological Observables

abstract: Cosmological data provide a powerful tool in the search for physics beyond the Standard Model (SM). An interesting target are light relics, new degrees of freedom which decoupled from the SM while relativistic. Nearly massless relics contribute to the radiation energy budget, and are commonly parametrized as variations in the effective number $N_{\rm eff}$ of neutrino species. Additionally, relics with masses greater than $10^{-4}$ eV become non-relativistic before today, and thus behave as matter instead of radiation. This leaves an imprint in the clustering of the large-scale structure of the universe, as light relics have important streaming motions, mirroring the case of massive neutrinos. Here we forecast how well current and upcoming cosmological surveys can probe light massive relics (LiMRs). We consider minimal extensions to the SM by both fermionic and bosonic relic degrees of freedom. By combining current and upcoming cosmic-microwave-background and large-scale-structure surveys, we forecast the significance at which each LiMR, with different masses and temperatures, can be detected. We find that a very large coverage of parameter space will be attainable by upcoming experiments, opening the possibility of exploring uncharted territory for new physics beyond the SM.


title: A Black-Hole Feedback Valve in Massive Galaxies

abstract: Star formation in the universe's most massive galaxies proceeds furiously early in time but then nearly ceases. Plenty of hot gas remains available but does not cool and condense into star-forming clouds. Active galactic nuclei (AGN) release enough energy to inhibit cooling of the hot gas, but energetic arguments alone do not explain why quenching of star formation is most effective in high-mass galaxies. In fact, optical observations show that quenching is more closely related to a galaxy's central stellar velocity dispersion ($\sigma_v$) than to any other characteristic. Here, we show that high $\sigma_v$ is critical to quenching because a deep central potential well maximizes the efficacy of AGN feedback. In order to remain quenched, a galaxy must continually sweep out the gas ejected from its aging stars. Supernova heating can accomplish this task as long as the AGN sufficiently reduces the gas pressure of the surrounding circumgalactic medium (CGM). We find that CGM pressure acts as the control knob on a valve that regulates AGN feedback and suggest that feedback power self-adjusts so that it suffices to lift the CGM out of the galaxy's potential well. Supernova heating then drives a galactic outflow that remains homogeneous if $\sigma_v \gtrsim 240 \, {\rm km \, s^{-1}}$. AGN feedback can effectively quench galaxies with a comparable velocity dispersion, but feedback in galaxies with a much lower velocity dispersion tends to result in convective circulation and accumulation of multiphase gas within the galaxy.


title: Stochastic modelling of star-formation histories II: star-formation variability from molecular clouds and gas inflow

abstract: A key uncertainty in galaxy evolution is the physics regulating star formation, ranging from small-scale processes related to the life-cycle of molecular clouds within galaxies to large-scale processes such as gas accretion onto galaxies. We study the imprint of such processes on the time-variability of star formation with an analytical approach tracking the gas mass of galaxies ("regulator model"). Specifically, we quantify the strength of the fluctuation in the star-formation rate (SFR) on different timescales, i.e. the power spectral density (PSD) of the star-formation history, and connect it to gas inflow and the life-cycle of molecular clouds. We show that in the general case the PSD of the SFR has three breaks, corresponding to the correlation time of the inflow rate, the equilibrium timescale of the gas reservoir of the galaxy, and the average lifetime of individual molecular clouds. On long and intermediate timescales (relative to the dynamical timescale of the galaxy), the PSD is typically set by the variability of the inflow rate and the interplay between outflows and gas depletion. On short timescales, the PSD shows an additional component related to the life-cycle of molecular clouds, which can be described by a damped random walk with a power-law slope of $\beta\approx2$ at high frequencies with a break near the average cloud lifetime. We discuss star-formation "burstiness" in a wide range of galaxy regimes, study the evolution of galaxies about the main sequence ridgeline, and explore the applicability of our method for understanding the star-formation process on cloud-scale from galaxy-integrated measurements.


title: Hydrodynamic attractor of a hybrid viscous fluid in Bjorken flow

abstract: The nonequilibrium evolution in a boost-invariant Bjorken flow of a hybrid viscous fluid model containing two interacting components with different viscosities, such that they represent strongly and weakly self-coupled sectors, is shown to be characterized by a hydrodynamic attractor which has an early-time behavior that is reminiscent of the so-called bottom-up thermalization scenario in heavy-ion collisions. The hydrodynamization times for the two sectors can differ strongly, with details depending on the curve realized on the two-dimensional attractor surface, which might account for different scenarios for small and large systems in nuclear collisions. The total system behaves like a single viscous fluid with a dynamically determined effective shear viscosity.


title: The Impact of Different Parameterizations on the Interpretation of CP Violation in Neutrino Oscillations

abstract: CP violation in the lepton mass matrix will be probed with good precision in upcoming experiments. The amount of CP violation present in oscillations can be quantified in numerous ways and is typically parameterized by the complex phase $\delta_{\rm PDG}$ in the standard PDG definition of the lepton mixing matrix. There are additional parameterizations of the lepton mixing matrix as well. Through various examples, we explore how, given the current data, different parameterizations can lead to different conclusions when working with parameterization dependent variables, such as $\delta$. We demonstrate how the smallness of $|U_{e3}|$ governs the scale of these results. We then demonstrate how $\delta$ can be misleading and argue that the Jarlskog is the cleanest means of presenting the amount of CP violation in the lepton sector. We also confirm that, among the different parameterizations considered, the standard PDG parameterization has a number of convenient features.


title: Event Horizon Wrinklification

abstract: The possible existence of stable black holes with entropies larger than the corresponding Schwarzschild black hole has been discussed extensively. The recently proposed "rough" black holes provide a concrete example of this. The fear is that, in accordance with the Second Law of thermodynamics, the familiar smooth-skinned black holes might spontaneously "wrinklify" into such an object. We show that this fear is to some extent justified, in the sense that AdS black holes with more entropy than the AdS-Schwarzschild black hole of the same mass do exist.


title: Searching for Dark Matter Signals in Timing Spectra at Neutrino Experiments

abstract: The sensitivity to dark matter signals at neutrino experiments is fundamentally challenged by the neutrino rates, as they leave similar signatures in their detectors. As a way to improve the signal sensitivity, we investigate a dark matter search strategy which utilizes the timing and energy spectra to discriminate dark matter from neutrino signals at low-energy, pulsed-beam neutrino experiments. This strategy was proposed in our companion paper arXiv:1906.10745, which we apply to potential searches at COHERENT, JSNS$^2$, and CCM. These experiments are not only sources of neutrinos but also high intensity sources of photons. The dark matter candidate of interest comes from the relatively prompt decay of a dark sector gauge boson which may replace a Standard-Model photon, so the delayed neutrino events can be suppressed by keeping prompt events only. Furthermore, prompt neutrino events can be rejected by a cut in recoil energy spectra, as their incoming energy is relatively small and bounded from above while dark matter may deposit a sizable energy beyond it. We apply the search strategy of imposing a combination of energy and timing cuts to the existing CsI data of the COHERENT experiment as a concrete example, and report a mild excess beyond known backgrounds. We then investigate the expected sensitivity reaches to dark matter signals in our benchmark experiments.


title: Simplified fast detector simulation in MadAnalysis 5

abstract: We introduce a new simplified fast detector simulator in the MadAnalysis 5 platform. The Python-like interpreter of the programme has been augmented by new commands allowing for a detector parametrisation through smearing and efficiency functions. On run time, an associated C++ code is automatically generated and executed to produce reconstructed-level events. In addition, we have extended the MadAnalysis 5 recasting infrastructure to support our detector emulator, and we provide predefined LHC detector configurations. We have compared predictions obtained with our approach to those resulting from the usage of the Delphes 3 software, both for Standard Model processes and a few new physics signals. Results generally agree to a level of about 10% or better, the largest differences in the predictions stemming from the different strategies that are followed to model specific detector effects. Equipped with these new functionalities, MadAnalysis 5 now offers a new user-friendly way to include detector effects when analysing collider events, the simulation of the detector and the analysis being both handled either through a set of intuitive Python commands or directly within the C++ core of the platform.


title: SAGAN-II : Molecular gas content of giant radio galaxies

abstract: Radio galaxies with jets of relativistic particles are usually hosted by massive elliptical galaxies with active nuclei powered by accretion of interstellar matter onto a supermassive black hole. In some rare cases (<5%), their jets drive the overall structure to sizes larger than 700 kpc, and they are called giant radio galaxies (GRGs). A very small fraction of the population of such radio galaxies contains molecular and atomic gas in the form of rings or discs that can fuel star formation. The origin of this gas is not well known; it has sometimes been associated with a minor merger with a gas-rich disc galaxy (e.g. Centaurus A) or cooling of material from a hot X-ray atmosphere (e.g. cooling flows). The giant radio jets might be the extreme evolution of these objects, and they can teach us about the radio galaxy evolution. We selected 12 targets from a catalogue of 820 GRGs that are likely to be in a gas-accretion and star formation phase. The targets were selected from the mid-infrared to contain heated dust. We report here the results of IRAM-30m observations, the molecular gas content, and the star formation efficiency, and we discuss the origin of the gas and disc morphology. Three out of our 12 targets are detected, and for the others, we report significant upper limits. We combine our three detections and upper limits with four additional detected GRGs from the literature to discuss the results. Most of the GRG targets belong to the main sequence, and a large fraction are in the passive domain. Their star formation efficiency is comparable to normal galaxies, except for two galaxies that are deficient in molecular gas with a short (~200Myr) depletion time, and a quiescent gas-rich giant spiral galaxy. In general, the depletion time is much longer than the lifetime of the giant radio jet.


title: Reheating-Induced Axion Dark Matter After Low Scale Inflation

abstract: A kinetic mixing between the axion and the inflaton allows for a production of axion dark matter even if the inflationary Hubble scale is smaller than the zero-temperature axion mass. We analyze the axion dynamics in this recently discovered "inflaxion" framework, and present a new cosmological scenario where the axion drifts away from its vacuum during the reheating epoch, giving rise to the observed dark matter abundance. We discuss the implications for both the QCD axion and axion-like particles.


title: The Unreasonable Effectiveness of Higher-Derivative Supergravity in AdS$_4$ Holography

abstract: We study four-derivative corrections to four-dimensional $\mathcal{N}=2$ minimal gauged supergravity and show that they are controlled by two real constants. The solutions of the equations of motion in the two-derivative theory are not modified by the higher-derivative corrections. We use this to derive a general formula for the regularized on-shell action for any asymptotically locally AdS$_4$ solution of the theory and show how the higher-derivative corrections affect black hole thermodynamic quantities in a universal way. We employ our results in the context of holography to derive explicit expressions for the subleading corrections in the large $N$ expansion of supersymmetric partition functions on various compact manifolds for a large class of three-dimensional SCFTs.


title: Two-component mass models of the lensing galaxy in the quadruply imaged supernova iPTF16geu

abstract: The first resolved, multiply imaged supernova Type Ia, iPTF16geu, was observed 4 years ago, five decades after such systems were first envisioned. Because of the unique properties of the source, these systems hold a lot of promise for the study of galaxy structure and cosmological parameters. However, this very first example presented modelers with a few puzzles. It was expected that to explain image fluxes a contribution from microlensing by stars would be required, but to accommodate the magnitude of microlensing, the density slope of the elliptical power law lens model had to be quite shallow, $\rho_{2D} \propto r^{-0.7}$. Furthermore, the center of mass had to be displaced from that of observed light by ~0.1 kpc, and the position angle of light distribution was misaligned with that of mass by ~40 degrees. In this paper we present mass models that resolve the first two problems, and suggest a resolution of the third. Motivated by observations of local ellipticals, and some recent analysis of galaxy-scale lenses, our mass models consist of two offset (baryonic) mass components. The resulting mass distributions have a single centroid, but are lopsided, and have isodensity contours that are not purely elliptical and not self-similar with radius. For many of our models the microlensing requirements are modest, and the ring formed by the extended supernova host galaxy resembles the observed one.


title: Temperature enhancement of thermal Hall conductance quantization

abstract: The quest for non-Abelian quasiparticles has inspired decades of experimental and theoretical efforts, where the scarcity of direct probes poses a key challenge. Among their clearest signatures is a thermal Hall conductance with quantized half-integer value in natural units $ \pi^2 k_B^2 T /3 h$ ($T$ is temperature, $h$ the Planck constant, $k_B$ the Boltzmann constant). Such a value was indeed recently observed in a quantum-Hall system and a magnetic insulator. We show that a non-topological "thermal metal" phase that forms due to quenched disorder may disguise as a non-Abelian phase by well approximating the trademark quantized thermal Hall response. Remarkably, the quantization here improves with temperature, in contrast to fully gapped systems. We provide numerical evidence for this effect and discuss its possible implications for the aforementioned experiments.


title: Nonrelativistic model of tetraquarks and predictions for their masses from fits to charmed and bottom meson data

abstract: We investigate a nonrelativistic model of tetraquarks, which are assumed to be compact and to consist of diquark-antidiquark pairs. We fit, for the first time, basically all currently known values for the measured masses of 45 mesons, including both charmed and bottom mesons, to the model and predict masses of tetraquarks as well as diquarks. In particular, we find masses of four axial-vector diquarks, i.e., $qc$, $cc$, $qb$, and $bb$, where $q = u,d$, and 24 ground-state tetraquarks, including both heavy-light tetraquarks ($qc\overline{qc}$ and $qb\overline{qb}$) and heavy tetraquarks ($cc\overline{cc}$ and $bb\overline{bb}$). In general, our results for the masses of $qb\overline{qb}$, $cc\overline{cc}$, and $bb\overline{bb}$ are largely comparable with other reported results, whereas our results for the masses of $qc\overline{qc}$ are slightly larger than what has been found earlier. Finally, we identify some of the obtained predictions for masses of tetraquarks with masses of experimental tetraquark candidates, and especially, we find that $\psi(4660)$, $Z_b(10610)$, and $Z_b(10650)$ could be described by the model.


title: Retrieving scattering clouds and disequilibrium chemistry in the atmosphere of HR 8799e

abstract: Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.


title: Accurately Weighing Neutrinos with Cosmological Surveys

abstract: A promising avenue to measure the total, and potentially individual, mass of neutrinos consists of leveraging cosmological datasets, such as the cosmic microwave background and surveys of the large-scale structure of the universe. In order to obtain unbiased estimates of the neutrino mass, however, many effects ought to be included. Here we forecast, via a Markov Chain Monte Carlo likelihood analysis, whether measurements by two galaxy surveys: DESI and {\it Euclid}, when added to the CMB-S4 experiment, are sensitive to two effects that can alter neutrino-mass measurements. The first is the slight difference in the suppression of matter fluctuations that each neutrino-mass hierarchy generates, at fixed total mass. The second is the growth-induced scale-dependent bias (GISDB) of haloes produced by massive neutrinos. We find that near-future surveys can distinguish hierarchies with the same total mass only at the $1\sigma$ level; thus, while these are poised to deliver a measurement of the sum of neutrino masses, they cannot significantly discern the mass of each individual neutrino in the foreseeable future. We further find that neglecting the GISDB induces up to a $1\sigma$ overestimation of the total neutrino mass, and we show how to absorb this effect via a redshift-dependent parametrization of the scale-independent bias.


title: Density Deconvolution with Normalizing Flows

abstract: Density deconvolution is the task of estimating a probability density function given only noise-corrupted samples. We can fit a Gaussian mixture model to the underlying density by maximum likelihood if the noise is normally distributed, but would like to exploit the superior density estimation performance of normalizing flows and allow for arbitrary noise distributions. Since both adjustments lead to an intractable likelihood, we resort to amortized variational inference. We demonstrate some problems involved in this approach, however, experiments on real data demonstrate that flows can already out-perform Gaussian mixtures for density deconvolution.