title: Phenomenology of $\tau^- \rightarrow \pi^- \nu_{\tau} \gamma$ using LCSR

abstract: We present the study of radiative tau decay ($\tau^- \rightarrow \pi^- \nu_\tau \gamma$), computing the structure dependent contribution using Light Cone Sum Rules. This decay includes the same form factors as the radiative pion decay with the crucial difference that the momentum transfer squared, $t$, between the pion-photon system is positive, which makes these form factors timelike and also as $t$ can now take values upto $m_\tau^2$ , it can produce real hadronic resonances. The analytical form for these form factors has been calculated using light cone sum rules method and the invariant mass-spectrum in the $\pi-\gamma$ system and the decay width are presented. The structure dependent parameter, $\gamma$, the ratio of the axial to vector form factor is found to be in good agreement with the experimental determination.


title: Cardiovascular deconditioning during long-term spaceflight through multiscale modeling

abstract: Human spaceflight has been fascinating man for centuries, representing the intangible need to explore the unknown, challenge new frontiers, advance technology and push scientific boundaries further. A key area of importance is cardiovascular deconditioning, that is, the collection of hemodynamic changes - from blood volume shift and reduction to altered cardiac function - induced by sustained presence in microgravity. A thorough grasp of the 0G adjustment point per se is important from a physiological viewpoint and fundamental for astronauts' safety and physical capability on long spaceflights. However, hemodynamic details of cardiovascular deconditioning are incomplete, inconsistent and poorly measured to date; thus a computational approach can be quite valuable. We present a validated 1D-0D multiscale model to study the cardiovascular response to long-term 0G spaceflight in comparison to the 1G supine reference condition. Cardiac work, oxygen consumption and contractility indexes, as well as central mean and pulse pressures were reduced, augmenting the cardiac deconditioning scenario. Exercise tolerance of a spaceflight traveler was found to be comparable to an untrained person with a sedentary lifestyle. At the capillary-venous level significant waveform alterations were observed which can modify the regular perfusion and average nutrient supply at the cellular level. The present study suggests special attention should be paid to future long spaceflights which demand prompt physical capacity at the time of restoration of partial gravity (e.g., Moon/Mars landing). Since spaceflight deconditioning has features similar to accelerated aging understanding deconditioning mechanisms in microgravity are also relevant to the understanding of aging physiology on Earth.


title: Transcription-driven genome organization: a model for chromosome structure and the regulation of gene expression tested through simulations

abstract: Current models for the folding of the human genome see a hierarchy stretching down from chromosome territories, through A/B compartments and TADs (topologically-associating domains), to contact domains stabilized by cohesin and CTCF. However, molecular mechanisms underlying this folding, and the way folding affects transcriptional activity, remain obscure. Here we review physical principles driving proteins bound to long polymers into clusters surrounded by loops, and present a parsimonious yet comprehensive model for the way the organization determines function. We argue that clusters of active RNA polymerases and their transcription factors are major architectural features; then, contact domains, TADs, and compartments just reflect one or more loops and clusters. We suggest tethering a gene close to a cluster containing appropriate factors -- a transcription factory -- increases the firing frequency, and offer solutions to many current puzzles concerning the actions of enhancers, super-enhancers, boundaries, and eQTLs (expression quantitative trait loci). As a result, the activity of any gene is directly influenced by the activity of other transcription units around it in 3D space, and this is supported by Brownian-dynamics simulations of transcription factors binding to cognate sites on long polymers.


title: Low-temperature thermal conductivity of Co$_{1-x}$M$_x$Si (M=Fe, Ni) alloys

abstract: We study the low-temperature electrical and thermal conductivity of CoSi and Co$_{1-x}$M$_x$Si alloys (M = Fe, Ni; $x \leq$ 0.06). Measurements show that the low-temperature electrical conductivity of Co$_{1-x}$Fe$_{x}$Si alloys decreases at $x > $ 0.01 by an order of magnitude compared with that of pure CoSi. It was expected that both the lattice and electronic contributions to thermal conductivity would decrease in the alloys. However, our experimental results revealed that at temperatures below 20K the thermal conductivity of Fe- and Ni-containing alloys is several times larger than that of pure CoSi. We discuss possible mechanisms of the thermal conductivity enhancement. The most probable one is related to the dominant scattering of phonons by charge carriers. We propose a simple theoretical model that takes into account the complex semimetallic electronic structure of CoSi with nonequivalent valleys, and show that it explains well the increase of the lattice thermal conductivity with increasing disorder and the linear temperature dependence of the thermal conductivity in the Co$_{1-x}$Fe$_x$Si alloys below 20K.


title: Large deviation expansions for the coefficients of random walks on the general linear group

abstract: Let $(g_n)_{n\geq 1}$ be a sequence of independent and identically distributed elements of the general linear group $GL(d, \mathbb R)$. Consider the random walk $G_n: = g_n \ldots g_1$. Under suitable conditions, we establish Bahadur-Rao-Petrov type large deviation expansion for the coefficients $\langle f, G_n v \rangle$, where $f \in (\mathbb R^d)^*$ and $v \in \mathbb R^d$. In particular, our result implies the large deviation principle with an explicit rate function, thus improving significantly the large deviation bounds established earlier. Moreover, we establish Bahadur-Rao-Petrov type large deviation expansion for the coefficients $\langle f, G_n v \rangle$ under the changed measure. Toward this end we prove the H\"{o}lder regularity of the stationary measure corresponding to the Markov chain $G_n v /|G_n v|$ under the changed measure, which is of independent interest. In addition, we also prove local limit theorems with large deviations for the coefficients of $G_n$.


title: EigenGame: PCA as a Nash Equilibrium

abstract: We present a novel view on principal component analysis (PCA) as a competitive game in which each approximate eigenvector is controlled by a player whose goal is to maximize their own utility function. We analyze the properties of this PCA game and the behavior of its gradient based updates. The resulting algorithm -- which combines elements from Oja's rule with a generalized Gram-Schmidt orthogonalization -- is naturally decentralized and hence parallelizable through message passing. We demonstrate the scalability of the algorithm with experiments on large image datasets and neural network activations. We discuss how this new view of PCA as a differentiable game can lead to further algorithmic developments and insights.


title: Nuclear parton distribution functions with uncertainties in a general mass variable flavor number scheme

abstract: In this article we obtain a new set of nuclear parton distribution functions (nuclear PDFs) at next-to-leading order and next-to-next-to-leading order accuracy in perturbative QCD. The common nuclear deep-inelastic scattering (DIS) data analyzed in our study are complemented by the available charged-current neutrino DIS data with nuclear targets and data from Drell-Yan cross-section measurements for several nuclear targets. In addition, the most recent DIS data from the Jefferson Lab CLAS and Hall C experiments are also added to our data sample. For these specific datasets, we consider the impact of target mass corrections and higher twist effects which are expected to be important in the region of large $x$ and intermediate-to-low $Q^2$. Our analysis is based on a publicly available open-source tool {\tt APFEL}, which has been modified to be applicable for our analysis of nuclear PDFs. Heavy quark contributions to nuclear DIS are considered within the framework of the {\tt FONLL} general-mass variable-flavor-number scheme. The most recent {\tt CT18} PDFs are used as baseline proton PDFs. The uncertainties of nuclear PDFs are determined using the standard Hessian approach. The results of our global QCD analysis are compared with existing nuclear PDF sets and with the fitted cross-sections, for which our set of nuclear PDFs provides a very good description.


title: Injective Hyperbolicity for quotients of balls and polydisks

abstract: In this article we study the injective Kobayashi metric on complex surfaces.


title: Berry-Esseen bounds and moderate deviations for the norm, entries and spectral radius of products of positive random matrices

abstract: Let $(g_{n})_{n\geq 1}$ be a sequence of independent and identically distributed positive random $d\times d$ matrices and consider the matrix product $G_n: = g_n \ldots g_1$. Under suitable conditions, we establish the Berry-Esseen bounds on the rate of convergence in the central limit theorem and moderate deviation expansions of Cram\'er type, for the matrix norm $\| G_n \|$ of $G_n$, for its $(i,j)$-th entry $G_n^{i,j}$, and the and for its spectral radius $\rho(G_n)$.


title: Arc Flow Formulations Based on Dynamic Programming: Theoretical Foundations and Applications

abstract: Network flow formulations are among the most successful tools to solve optimization problems. Such formulations correspond to determining an optimal flow in a network. One particular class of network flow formulations is the arc flow, where variables represent flows on individual arcs of the network. For $\mathcal{NP}$-hard problems, polynomial-sized arc flow models typically provide weak linear relaxations and may have too much symmetry to be efficient in practice. Instead, arc flow models with a pseudo-polynomial size usually provide strong relaxations and are efficient in practice. The interest in pseudo-polynomial arc flow formulations has grown considerably in the last twenty years, in which they have been used to solve many open instances of hard problems. A remarkable advantage of pseudo-polynomial arc flow models is the possibility to solve practical-sized instances directly by a Mixed Integer Linear Programming solver, avoiding the implementation of complex methods based on column generation. In this survey, we present theoretical foundations of pseudo-polynomial arc flow formulations, by showing a relation between their network and Dynamic Programming (DP). This relation allows a better understanding of the strength of these formulations, through a link with models obtained by Dantzig-Wolfe decomposition. The relation with DP also allows a new perspective to relate state-space relaxation methods for DP with arc flow models. We also present a dual point of view to contrast the linear relaxation of arc flow models with that of models based on paths and cycles. To conclude, we review the main solution methods and applications of arc flow models based on DP in several domains such as cutting, packing, scheduling, and routing.


title: An In-Plane $J_2$-Invariance Condition and Control Algorithm for Highly Elliptical Satellite Formations

abstract: Phased formations of satellites provide an important means of keeping multiple satellites in close proximity without becoming dangerously close to one another. In order to minimize the amount of propellant necessary to keep a set of satellites in a phased formation, this paper presents a new condition for in-plane linearized secular J2-invariance in highly eccentric orbits. A maintenance strategy outlined here combats the growth of out-of-plane motion, and the delta-v cost of this strategy is analyzed. For context, this paper also discusses relevant formations that benefit from this condition and maintenance strategy. A highly elliptical perching formation called the boomerang perch is also discussed.


title: Dynamic Facial Asset and Rig Generation from a Single Scan

abstract: The creation of high-fidelity computer-generated (CG) characters used in film and gaming requires intensive manual labor and a comprehensive set of facial assets to be captured with complex hardware, resulting in high cost and long production cycles. In order to simplify and accelerate this digitization process, we propose a framework for the automatic generation of high-quality dynamic facial assets, including rigs which can be readily deployed for artists to polish. Our framework takes a single scan as input to generate a set of personalized blendshapes, dynamic and physically-based textures, as well as secondary facial components (e.g., teeth and eyeballs). Built upon a facial database consisting of pore-level details, with over $4,000$ scans of varying expressions and identities, we adopt a self-supervised neural network to learn personalized blendshapes from a set of template expressions. We also model the joint distribution between identities and expressions, enabling the inference of the full set of personalized blendshapes with dynamic appearances from a single neutral input scan. Our generated personalized face rig assets are seamlessly compatible with cutting-edge industry pipelines for facial animation and rendering. We demonstrate that our framework is robust and effective by inferring on a wide range of novel subjects, and illustrate compelling rendering results while animating faces with generated customized physically-based dynamic textures.


title: SiMPLISTIC: A Novel Pairwise Potential for Implicit Solvent Lipid Simulations with Single-site Models

abstract: Implicit solvent, coarse-grained models with pairwise interactions can access the largest length and time scales in molecular dynamics simulations, owing to the absence of interactions with a huge number of solvent particles, the smaller number of interaction sites in the model molecules, and the lack of fast sub-molecular degrees of freedom. In this paper, we describe a maximally coarse-grained model for lipids in implicit water. The model is called SiMPLISTIC, which abbreviates for Single-site Model with Pairwise interaction for Lipids in Implicit Solvent with Tuneable Intrinsic Curvature. SiMPLISTIC lipids rapidly self-assemble into realistic non-lamellar and lamellar phases such as inverted micelles and bilayers, the spontaneous curvature of the phase being determined by a single free parameter of the model. Model membrane simulations with the lamellar lipids show satisfactory fluid and gel phases with no interdigitation or tilt. The model lipids follow rigid body dynamics suggested by empirical studies, and generate bilayer elastic properties consistent with experiments and other simulations. SiMPLISTIC can also simulate mixtures of lipids that differ in their packing parameter or length, the latter leading to the phenomenon of hydrophobic mismatch driven domain formation. The model has a large scope due to its speed, conceptual and computational simplicity, and versatility. Applications may range from large-scale simulations for academic and industrial research on various lipid-based systems, such as lyotropic liquid crystals, biological and biomimetic membranes, vectors for drug and gene delivery, to fast, lightweight, interactive simulations for gaining insights into self-assembly, lipid polymorphism, biomembrane organization etc.


title: ISAAQ -- Mastering Textbook Questions with Pre-trained Transformers and Bottom-Up and Top-Down Attention

abstract: Textbook Question Answering is a complex task in the intersection of Machine Comprehension and Visual Question Answering that requires reasoning with multimodal information from text and diagrams. For the first time, this paper taps on the potential of transformer language models and bottom-up and top-down attention to tackle the language and visual understanding challenges this task entails. Rather than training a language-visual transformer from scratch we rely on pre-trained transformers, fine-tuning and ensembling. We add bottom-up and top-down attention to identify regions of interest corresponding to diagram constituents and their relationships, improving the selection of relevant visual information for each question and answer options. Our system ISAAQ reports unprecedented success in all TQA question types, with accuracies of 81.36%, 71.11% and 55.12% on true/false, text-only and diagram multiple choice questions. ISAAQ also demonstrates its broad applicability, obtaining state-of-the-art results in other demanding datasets.


title: Uniform potential density for rational points on algebraic groups and elliptic K3 surfaces

abstract: A collection of varieties satisfies uniform potential density if each of them contains a dense subset of rational points after extending its ground field by a bounded degree. In this paper, we prove that uniform potential density holds for connected algebraic groups of a fixed dimension over fields of characteristic zero as well as elliptic K3 surfaces over number fields.


title: On the singular Weinstein conjecture and the existence of escape orbits for $b$-Beltrami fields

abstract: Motivated by Poincar\'e's orbits going to infinity in the (restricted) three-body (see [26] and [6]), we investigate the generic existence of heteroclinic-like orbits in a neighbourhood of the critical set of a $b$-contact form. This is done by using the singular counterpart [3] of Etnyre--Ghrist's contact/Beltrami correspondence [9], and genericity results concerning eigenfunctions of the Laplacian established by Uhlenbeck [29]. Specifically, we analyze the $b$-Beltrami vector fields on $b$-manifolds of dimension $3$ and prove that for a generic asymptotically exact $b$-metric they exhibit escape orbits. We also show that a generic asymptotically symmetric $b$-Beltrami vector field on an asymptotically flat $b$-manifold has a generalized singular periodic orbit and at least $4$ escape orbits. Generalized singular periodic orbits are trajectories of the vector field whose $\alpha$- and $\omega$-limit sets intersect the critical surface. These results are a first step towards proving the singular Weinstein conjecture.


title: Absence of true localization in many-body localized phases

abstract: We have recently shown that the logarithmic growth of the entanglement entropy following a quantum quench in a many-body localized (MBL) phase is accompanied by a slow growth of the number entropy, $S_N\sim\ln\ln t$. Here we provide an in-depth numerical study of $S_N(t)$ for the disordered Heisenberg chain and show that this behavior is not transient and persists even for very strong disorder. Calculating the truncated R\'enyi number entropy $S_N^{(\alpha)}(t)=(1-\alpha)^{-1}\ln\sum_n p^\alpha(n)$ for $\alpha\ll 1$ and $p(n)>p_c$ -- which is sensitive to large number fluctuations occurring with low probability -- we demonstrate that the particle number distribution $p(n)$ in one half of the system has a continuously growing tail. This indicates a slow but steady increase of the number of particles crossing between the partitions in the interacting case, and is in sharp contrast to Anderson localization, for which we show that $S_N^{(\alpha\to 0)}(t)$ saturates for any cutoff $p_c>0$. We show, furthermore, that the growth of $S_N$ is $\mathit not$ the consequence of rare states or rare regions but rather represents typical behavior. These findings provide strong evidence that the interacting system is never fully localized even for very strong but finite disorder.


title: Where to stand when playing darts?

abstract: This paper analyzes the question of where one should stand when playing darts. If one stands at distance $d>0$ and aims at $a\in \mathbb{R}^n$, then the dart (modelled by a random vector $X$ in $\mathbb{R}^n$) hits a random point given by $a+dX$. Next, given a payoff function $f$, one considers $$ \sup_a Ef(a+dX) $$ and asks if this is decreasing in $d$; i.e., whether it is better to stand closer rather than farther from the target. Perhaps surprisingly, this is not always the case and understanding when this does or does not occur is the purpose of this paper. We show that if $X$ has a so-called selfdecomposable distribution, then it is always better to stand closer for any payoff function. This class includes all stable distributions as well as many more. On the other hand, if the payoff function is $\cos(x)$, then it is always better to stand closer if and only if the characteristic function $|\phi_X(t)|$ is decreasing on $[0,\infty)$. We will then show that if there are at least two point masses, then it is not always better to stand closer using $\cos(x)$. If there is a single point mass, one can find a different payoff function to obtain this phenomenon. Another large class of darts $X$ for which there are bounded continuous payoff functions for which it is not always better to stand closer are distributions with compact support. This will be obtained by using the fact that the Fourier transform of such distributions has a zero in the complex plane. This argument will work whenever there is a complex zero of the Fourier transform. Finally, we analyze if the property of it being better to stand closer is closed under convolution and/or limits.


title: Deep learning for time series classification

abstract: Time series analysis is a field of data science which is interested in analyzing sequences of numerical values ordered in time. Time series are particularly interesting because they allow us to visualize and understand the evolution of a process over time. Their analysis can reveal trends, relationships and similarities across the data. There exists numerous fields containing data in the form of time series: health care (electrocardiogram, blood sugar, etc.), activity recognition, remote sensing, finance (stock market price), industry (sensors), etc. Time series classification consists of constructing algorithms dedicated to automatically label time series data. The sequential aspect of time series data requires the development of algorithms that are able to harness this temporal property, thus making the existing off-the-shelf machine learning models for traditional tabular data suboptimal for solving the underlying task. In this context, deep learning has emerged in recent years as one of the most effective methods for tackling the supervised classification task, particularly in the field of computer vision. The main objective of this thesis was to study and develop deep neural networks specifically constructed for the classification of time series data. We thus carried out the first large scale experimental study allowing us to compare the existing deep methods and to position them compared other non-deep learning based state-of-the-art methods. Subsequently, we made numerous contributions in this area, notably in the context of transfer learning, data augmentation, ensembling and adversarial attacks. Finally, we have also proposed a novel architecture, based on the famous Inception network (Google), which ranks among the most efficient to date.


title: Charge symmetry broken complex coacervation

abstract: Liquid-liquid phase separation has emerged as one of the important paradigms in the chemical physics as well as biophysics of charged macromolecular systems. We elucidate an equilibrium phase separation mechanism based on charge regulation, i.e., protonation-deprotonation equilibria controlled by pH, in an idealized macroion system which can serve as a proxy for simple coacervation. First, a low-density density-functional calculation reveals the dominance of two-particle configurations coupled by ion adsorption on neighboring macroions. Then a binary cell model, solved on the Debye-H\"uckel as well as the full nonlinear Poisson-Boltzmann level, unveils the charge-symmetry breaking as inducing the phase separation between low- and high-density phases as a function of pH. These results can be identified as a charge symmetry broken complex coacervation between chemically identical macroions.

categories: cond-mat.soft cond-mat.stat-mech physics.bio-ph physics.chem-ph


title: The Horizontal Branch morphology of the globular cluster NGC 1261 using AstroSat

abstract: We present the results obtained from the UV photometry of the globular cluster NGC 1261 using Far-UV (FUV) and Near-UV (NUV) images acquired with the Ultraviolet Imaging Telescope (UVIT) onboard the Astrosat satellite. We utilized the UVIT data combined with HST, GAIA, and ground-based optical photometric data to construct the different UV colour-magnitude diagrams (CMDs). We detected blue HB (BHB), and two extreme HB (EHB) stars in FUV, whereas full HB, i.e., red HB (RHB), BHB as well as EHB is detected in NUV CMDs. The 2 EHB stars, identified in both NUV and FUV, are confirmed members of the cluster. The HB stars form a tight sequence in UV-optical CMDs, which is almost aligned with Padova isochrones. This study sheds light on the significance of UV imaging to probe the HB morphology in GCs.


title: Spin fluctuations in quantized transport of magnetic topological insulators

abstract: In magnetic topological insulators, quantized electronic transport is interwined with spontaneous magnetic ordering, as magnetization controls band gaps, hence band topology, through the exchange interaction. We show that considering the exchange gaps at the mean-field level is inadequate to predict phase transitions between electronic states of distinct topology. Thermal spin fluctuations disturbing the magnetization can act as frozen disorders that strongly scatter electrons, reducing the onset temperature of quantized transport appreciably even in the absence of structural impurities. This effect, which has hitherto been overlooked, provides an alternative explanation of recent experiments on intrinsic magnetic topological insulators.


title: Understanding tables with intermediate pre-training

abstract: Table entailment, the binary classification task of finding if a sentence is supported or refuted by the content of a table, requires parsing language and table structure as well as numerical and discrete reasoning. While there is extensive work on textual entailment, table entailment is less well studied. We adapt TAPAS (Herzig et al., 2020), a table-based BERT model, to recognize entailment. Motivated by the benefits of data augmentation, we create a balanced dataset of millions of automatically created training examples which are learned in an intermediate step prior to fine-tuning. This new data is not only useful for table entailment, but also for SQA (Iyyer et al., 2017), a sequential table QA task. To be able to use long examples as input of BERT models, we evaluate table pruning techniques as a pre-processing step to drastically improve the training and prediction efficiency at a moderate drop in accuracy. The different methods set the new state-of-the-art on the TabFact (Chen et al., 2020) and SQA datasets.


title: Universal limiting transition temperature for the high $T_\mathrm{c}$ superconductors

abstract: Since their discovery three decades ago it has emerged that the physics of high-$T_\mathrm{c}$ cuprate superconductors is characterised by multiple temperature scales, and a phenomenology that deviates significantly from the conventional paradigm of superconductivity. Below the "pseudogap" temperature, $T^*$, a range of experiments indicate a reduction in the density of states. Lower in temperature, $T_\mathrm{c}$ marks the onset of a bulk zero-resistance state. Intermediate between these, numerous other temperature boundaries and cross-overs have been identified, often postulated to be associated with fluctuations of some kind. However, for the most part there is little consensus either over their definitions or the physical mechanisms at play. One of these temperature scales is the Nernst onset temperature, $T_\mathrm{onset}$, below which thermoelectric phenomena characteristic of superconductivity are observed. Depending on the material, the difference between $T_\mathrm{onset}$ and $T_\mathrm{c}$ ranges from almost nothing to over 100K. In this paper, we identify $T_\mathrm{onset}$ from published experimental data according to a consistent definition; this reveals a remarkable consistency of behaviour across the whole high-$T_\mathrm{c}$ family, despite the appreciable variations in other parameters. Our analysis suggests that in the cuprates there is an inherent limit to $T_\mathrm{c} \sim 135$K. We compare this behaviour with other strongly-correlated superconductors and propose a unified picture of superconducting fluctuations close to the Mott state.


title: DASGIL: Domain Adaptation for Semantic and Geometric-aware Image-based Localization

abstract: Long-Term visual localization under changing environments is a challenging problem in autonomous driving and mobile robotics due to season, illumination variance, etc. Image retrieval for localization is an efficient and effective solution to the problem. In this paper, we propose a novel multi-task architecture to fuse the geometric and semantic information into the multi-scale latent embedding representation for visual place recognition. To use the high-quality ground truths without any human effort, the effective multi-scale feature discriminator is proposed for adversarial training to achieve the domain adaptation from synthetic virtual KITTI dataset to real-world KITTI dataset. The proposed approach is validated on the Extended CMU-Seasons dataset and Oxford RobotCar dataset through a series of crucial comparison experiments, where our performance outperforms state-of-the-art baselines for retrieval-based localization and large-scale place recognition under the challenging environment.


title: An uncountable Mackey-Zimmer theorem

abstract: The Mackey-Zimmer theorem classifies ergodic group extensions $X$ of a measure-preserving system $Y$ by a compact group $K$, by showing that such extensions are isomorphic to a group skew-product $X \equiv Y \rtimes_\rho H$ for some closed subgroup $H$ of $K$. An analogous theorem is also available for ergodic homogeneous extensions $X$ of $Y$, namely that they are isomorphic to a homogeneous skew-product $Y \rtimes_\rho H/M$. These theorems have many uses in ergodic theory, for instance playing a key role in the Host-Kra structural theory of characteristic factors of measure-preserving systems. The existing proofs of the Mackey-Zimmer theorem require various "countability", "separability", or "metrizability" hypotheses on the group $\Gamma$ that acts on the system, the base space $Y$, and the group $K$ used to perform the extension. In this paper we generalize the Mackey-Zimmer theorem to "uncountable" settings in which these hypotheses are omitted, at the cost of making the notion of a measure-preserving system and a group extension more abstract. However, this abstraction is partially counteracted by the use of a "canonical model" for abstract measure-preserving systems developed in a companion paper. In subsequent work we will apply this theorem to also obtain uncountable versions of the Host-Kra structural theory.


title: D3C: Reducing the Price of Anarchy in Multi-Agent Learning

abstract: In multiagent systems, the complex interaction of fixed incentives can lead agents to outcomes that are poor (inefficient) not only for the group, but also for each individual. Price of anarchy is a technical, game-theoretic definition that quantifies the inefficiency arising in these scenarios -- it compares the welfare that can be achieved through perfect coordination against that achieved by self-interested agents at a Nash equilibrium. We derive a differentiable, upper bound on a price of anarchy that agents can cheaply estimate during learning. Equipped with this estimator, agents can adjust their incentives in a way that improves the efficiency incurred at a Nash equilibrium. Agents do so by learning to mix their reward (equiv. negative loss) with that of other agents by following the gradient of our derived upper bound. We refer to this approach as D3C. In the case where agent incentives are differentiable, D3C resembles the celebrated Win-Stay, Lose-Shift strategy from behavioral game theory, thereby establishing a connection between the global goal of maximum welfare and an established agent-centric learning rule. In the non-differentiable setting, as is common in multiagent reinforcement learning, we show the upper bound can be reduced via evolutionary strategies, until a compromise is reached in a distributed fashion. We demonstrate that D3C improves outcomes for each agent and the group as a whole on several social dilemmas including a traffic network exhibiting Braess's paradox, a prisoner's dilemma, and several multiagent domains.


title: Stratified Global MHD Models of Accretion Disks in Semi-Detached Binaries

abstract: We present results of the first global magnetohydrodynamic (MHD) simulations of accretion disks fed by Roche lobe overflow, including vertical stratification, in order to investigate the roles of spiral shocks, magnetorotational instability (MRI), and the accretion stream on disk structure and evolution. Our models include a simple treatment of gas thermodynamics, with orbital Mach numbers at the inner edge of the disk $M_{\rm in}$ of 5 and 10. We find mass accretion rates to vary considerably on all time scales, with only the Mach 5 model reaching a clear quasi-stationary state. For Mach 10, the model undergoes an outside-in magnetically-driven accretion event occurring on a time scale of $\sim10$ orbital periods of the binary. Both models exhibit spiral shocks inclined with respect to the binary plane, with their position and inclination changing rapidly. However, the time-averaged location of these shocks in the equatorial plane is well-fit by simple linear models. MRI turbulence in the disk generates toroidal magnetic field patterns (butterfly diagrams) that are in some cases irregular, perhaps due to interaction with spiral structure. While many of our results are in good agreement with local studies, we find some features (most notably those related to spiral shocks) can only be captured in global models such as studied here. Thus, while global studies remain computationally expensive -- even as idealized models -- they are essential (along with more sophisticated treatment of radiation transport and disk thermodynamics) for furthering our understanding of accretion in binary systems.


title: Interpreting Graph Neural Networks for NLP With Differentiable Edge Masking

abstract: Graph neural networks (GNNs) have become a popular approach to integrating structural inductive biases into NLP models. However, there has been little work on interpreting them, and specifically on understanding which parts of the graphs (e.g. syntactic trees or co-reference structures) contribute to a prediction. In this work, we introduce a post-hoc method for interpreting the predictions of GNNs which identifies unnecessary edges. Given a trained GNN model, we learn a simple classifier that, for every edge in every layer, predicts if that edge can be dropped. We demonstrate that such a classifier can be trained in a fully differentiable fashion, employing stochastic gates and encouraging sparsity through the expected $L_0$ norm. We use our technique as an attribution method to analyze GNN models for two tasks -- question answering and semantic role labeling -- providing insights into the information flow in these models. We show that we can drop a large proportion of edges without deteriorating the performance of the model, while we can analyse the remaining edges for interpreting model predictions.


title: Understanding Self-supervised Learning with Dual Deep Networks

abstract: We propose a novel theoretical framework to understand contrastive self-supervised learning (SSL) methods that employ dual pairs of deep ReLU networks (e.g., SimCLR). First, we prove that in each SGD update of SimCLR with various loss functions, including simple contrastive loss, soft Triplet loss and InfoNCE loss, the weights at each layer are updated by a \emph{covariance operator} that specifically amplifies initial random selectivities that vary across data samples but survive averages over data augmentations. To further study what role the covariance operator plays and which features are learned in such a process, we model data generation and augmentation processes through a \emph{hierarchical latent tree model} (HLTM) and prove that the hidden neurons of deep ReLU networks can learn the latent variables in HLTM, despite the fact that the network receives \emph{no direct supervision} from these unobserved latent variables. This leads to a provable emergence of hierarchical features through the amplification of initially random selectivities through contrastive SSL. Extensive numerical studies justify our theoretical findings. Code is released in https://github.com/facebookresearch/luckmatters/tree/master/ssl.


title: Independent factorization of the last zero arcsine law for Bessel processes with drift

abstract: We show that the last zero before time $t$ of a recurrent Bessel process with drift starting at $0$ has the same distribution as the product of an independent right censored exponential random variable and a beta random variable. This extends a recent result of Schulte-Geers and Stadje (2017) from Brownian motion with drift to recurrent Bessel processes with drift. Our proof is intuitive and direct while avoiding heavy computations. For this we develop a novel additive decomposition for the square of a Bessel process with drift that may be of independent interest.


title: Ball packings for links

abstract: The ball number of a link $L$, denoted by $ball(L)$, is the minimum number of solid balls (not necessarily of the same size) needed to realize a necklace representing $L$. In this paper, we show that $ball(L)\leq 5 cr(L)$ where $cr(L)$ denotes the crossing number of $L$. To this end, we use Lorentz geometry applied to ball packings. The well-known Koebe-Andreev-Thurston circle packing Theorem is also an important brick for the proof. Our approach yields to an algorithm to construct explicitly the desired necklace representation of $L$ in the 3-dimensional space.


title: Emergent Social Learning via Multi-agent Reinforcement Learning

abstract: Social learning is a key component of human and animal intelligence. By taking cues from the behavior of experts in their environment, social learners can acquire sophisticated behavior and rapidly adapt to new circumstances. This paper investigates whether independent reinforcement learning (RL) agents in a multi-agent environment can learn to use social learning to improve their performance. We find that in most circumstances, vanilla model-free RL agents do not use social learning. We analyze the reasons for this deficiency, and show that by imposing constraints on the training environment and introducing a model-based auxiliary loss we are able to obtain generalized social learning policies which enable agents to: i) discover complex skills that are not learned from single-agent training, and ii) adapt online to novel environments by taking cues from experts present in the new environment. In contrast, agents trained with model-free RL or imitation learning generalize poorly and do not succeed in the transfer tasks. By mixing multi-agent and solo training, we can obtain agents that use social learning to gain skills that they can deploy when alone, even out-performing agents trained alone from the start.


title: Tunneling transport of unitary fermions across the superfluid transition

abstract: We investigate the transport of a Fermi gas with unitarity-limited interactions across the superfluid phase transition, probing its response to a direct current (dc) drive through a tunnel junction. As the superfluid critical temperature is crossed from below, we observe the evolution from a highly nonlinear to an Ohmic conduction characteristics, associated with the critical breakdown of the Josephson dc current induced by pair condensate depletion. Moreover, we reveal a large and dominant anomalous contribution to resistive currents, which reaches its maximum at the lowest attained temperature, fostered by the tunnel coupling between the condensate and phononic Bogoliubov-Anderson excitations. Increasing the temperature, while the zeroing of supercurrents marks the transition to the normal phase, the conductance drops considerably but remains much larger than that of a normal, uncorrelated Fermi gas tunneling through the same junction. We attribute such enhanced transport to incoherent tunneling of sound modes, which remain weakly damped in the collisional hydrodynamic fluid of unpaired fermions at unitarity.


title: Utilizing Transfer Learning and a Customized Loss Function for Optic Disc Segmentation from Retinal Images

abstract: Accurate segmentation of the optic disc from a retinal image is vital to extracting retinal features that may be highly correlated with retinal conditions such as glaucoma. In this paper, we propose a deep-learning based approach capable of segmenting the optic disc given a high-precision retinal fundus image. Our approach utilizes a UNET-based model with a VGG16 encoder trained on the ImageNet dataset. This study can be distinguished from other studies in the customization made for the VGG16 model, the diversity of the datasets adopted, the duration of disc segmentation, the loss function utilized, and the number of parameters required to train our model. Our approach was tested on seven publicly available datasets augmented by a dataset from a private clinic that was annotated by two Doctors of Optometry through a web portal built for this purpose. We achieved an accuracy of 99.78\% and a Dice coefficient of 94.73\% for a disc segmentation from a retinal image in 0.03 seconds. The results obtained from comprehensive experiments demonstrate the robustness of our approach to disc segmentation of retinal images obtained from different sources.


title: The Future of the Correlated Electron Problem

abstract: A central problem in modern condensed matter physics is the understanding of materials with strong electron correlations. Despite extensive work, the essential physics of many of these systems is not understood and there is very little ability to make predictions in this class of materials. In this manuscript we share our personal views on the major open problems in the field of correlated electron systems. We discuss some possible routes to make progress in this rich and fascinating field. This manuscript is the result of the vigorous discussions and deliberations that took place at Johns Hopkins University during a three-day workshop January 27, 28, and 29, 2020 that brought together six senior scientists and 46 more junior scientists. Our hope, is that the topics we have presented will provide inspiration for others working in this field and motivation for the idea that significant progress can be made on very hard problems if we focus our collective energies.


title: Wavenumber-explicit convergence of the $hp$-FEM for the full-space heterogeneous Helmholtz equation with smooth coefficients

abstract: A convergence theory for the $hp$-FEM applied to a variety of constant-coefficient Helmholtz problems was pioneered in the papers [Melenk-Sauter, 2010], [Melenk-Sauter, 2011], [Esterhazy-Melenk, 2012], [Melenk-Parsania-Sauter, 2013]. This theory shows that, if the solution operator is bounded polynomially in the wavenumber $k$, then the Galerkin method is quasioptimal provided that $hk/p \leq C_1$ and $p\geq C_2 \log k$, where $C_1$ is sufficiently small, $C_2$ is sufficiently large, and both are independent of $k,h,$ and $p$. The significance of this result is that if $hk/p= C_1$ and $p=C_2\log k$, then quasioptimality is achieved with the total number of degrees of freedom proportional to $k^d$; i.e., the $hp$-FEM does not suffer from the pollution effect. This paper proves the analogous quasioptimality result for the heterogeneous (i.e. variable-coefficient) Helmholtz equation, posed in $\mathbb{R}^d$, $d=2,3$, with the Sommerfeld radiation condition at infinity, and $C^\infty$ coefficients. We also prove a bound on the relative error of the Galerkin solution in the particular case of the plane-wave scattering problem. These are the first ever results on the wavenumber-explicit convergence of the $hp$-FEM for the Helmholtz equation with variable coefficients.


title: Maximum power heat engines and refrigerators in the fast-driving regime

abstract: We study the optimization of the performance of arbitrary periodically driven thermal machines. Within the assumption of fast modulation of the driving parameters, we derive the optimal cycle that universally maximizes the extracted power of heat engines, the cooling power of refrigerators, and in general any linear combination of the heat currents. We denote this optimal solution as ``generalized Otto cycle'' since it shares the basic structure with the standard Otto cycle, but it is characterized by a greater number of fast strokes. We bound this number in terms of the dimension of the Hilbert space of the system used as working fluid. The generality of these results allows for a widespread range of applications, such as reducing the computational complexity for numerical approaches, or obtaining the explicit form of the optimal protocols when the system-baths interactions are characterized by a single thermalization scale. In this case, we compare the thermodynamic performance of a collection of optimally driven non-interacting and interacting qubits. Remarkably, for refrigerators the non-interacting qubits perform almost as well as the interacting ones, while in the heat engine case there is a many-body advantage both in the maximum power, and in the efficiency at maximum power. Additionally, we illustrate our general results studying the paradigmatic model of a qutrit-based heat engine. Our results strictly hold in the semiclassical case in which no coherence is generated by the driving, and finally we discuss the non-commuting case.


title: Nearly Minimax Optimal Reinforcement Learning for Discounted MDPs

abstract: We study the reinforcement learning problem for discounted Markov Decision Processes (MDPs) under the tabular setting. We propose a model-based algorithm named UCBVI-$\gamma$, which is based on the \emph{optimism in the face of uncertainty principle} and the Bernstein-type bonus. We show that UCBVI-$\gamma$ achieves an $\tilde{O}\big({\sqrt{SAT}}/{(1-\gamma)^{1.5}}\big)$ regret, where $S$ is the number of states, $A$ is the number of actions, $\gamma$ is the discount factor and $T$ is the number of steps. In addition, we construct a class of hard MDPs and show that for any algorithm, the expected regret is at least $\tilde{\Omega}\big({\sqrt{SAT}}/{(1-\gamma)^{1.5}}\big)$. Our upper bound matches the minimax lower bound up to logarithmic factors, which suggests that UCBVI-$\gamma$ is nearly minimax optimal for discounted MDPs.


title: Twisted bilayer graphene.VI. An Exact Diagonalization Study of Twisted Bilayer Graphene at Non-Zero Integer Fillings

abstract: Using exact diagonalization, we study the projected Hamiltonian with Coulomb interaction in the 8 flat bands of first magic angle twisted bilayer graphene. Employing the U(4) (U(4)$\times$U(4)) symmetries in the nonchiral (chiral) flat band limit, we reduced the Hilbert space to an extent which allows for study around $\nu=\pm 3,\pm2,\pm1$ fillings. In the first chiral limit $w_0/w_1=0$ where $w_0$ ($w_1$) is the $AA$ ($AB$) stacking hopping, we find that the ground-states at these fillings are extremely well-described by Slater determinants in a so-called Chern basis, and the exactly solvable charge $\pm1$ excitations found in [arXiv:2009.14200] are the lowest charge excitations up to system sizes $8\times8$ (for restricted Hilbert space) in the chiral-flat limit. We also find that the Flat Metric Condition (FMC) used in [arXiv:2009.11301,2009.11872,2009.12376,2009.13530,2009.14200] for obtaining a series of exact ground-states and excitations holds in a large parameter space. For $\nu=-3$, the ground state is the spin and valley polarized Chern insulator with $\nu_C=\pm1$ at $w_0/w_1\lesssim0.9$ (0.3) with (without) FMC. At $\nu=-2$, we can only numerically access the valley polarized sector, and we find a spin ferromagnetic phase when $w_0/w_1\gtrsim0.5t$ where $t\in[0,1]$ is the factor of rescaling of the actual TBG bandwidth, and a spin singlet phase otherwise, confirming the perturbative calculation [arXiv:2009.13530]. The analytic FMC ground state is, however, predicted in the intervalley coherent sector which we cannot access [arXiv:2009.13530]. For $\nu=-3$ with/without FMC, when $w_0/w_1$ is large, the finite-size gap $\Delta$ to the neutral excitations vanishes, leading to phase transitions. Further analysis of the ground state momentum sectors at $\nu=-3$ suggests a competition among (nematic) metal, momentum $M_M$ ($\pi$) stripe and $K_M$-CDW orders at large $w_0/w_1$.


title: Recoverable Systems

abstract: Motivated by the established notion of storage codes, we consider sets of infinite sequences over a finite alphabet such that every $k$-tuple of consecutive entries is uniquely recoverable from its $l$-neighborhood in the sequence. We address the problem of finding the maximum growth rate of the set, which we term capacity, as well as constructions of explicit families that approach the optimal rate. The techniques that we employ rely on the connection of this problem with constrained systems. In the second part of the paper we consider a modification of the problem wherein the entries in the sequence are viewed as random variables over a finite alphabet that follow some joint distribution, and the recovery condition requires that the Shannon entropy of the $k$-tuple conditioned on its $l$-neighborhood be bounded above by some $\epsilon>0.$ We study properties of measures on infinite sequences that maximize the metric entropy under the recoverability condition. Drawing on tools from ergodic theory, we prove some properties of entropy-maximizing measures. We also suggest a procedure of constructing an $\epsilon$-recoverable measure from a corresponding deterministic system.


title: Quantifying social organization and political polarization in online platforms

abstract: Optimism about the Internet's potential to bring the world together has been tempered by concerns about its role in inflaming the 'culture wars'. Via mass selection into like-minded groups, online society may be becoming more fragmented and polarized, particularly with respect to partisan differences. However, our ability to measure the social makeup of online communities, and in turn understand the social organization of online platforms, is limited by the pseudonymous, unstructured, and large-scale nature of digital discussion. We develop a neural embedding methodology to quantify the positioning of online communities along social dimensions by leveraging large-scale patterns of aggregate behaviour. Applying our methodology to 5.1B Reddit comments made in 10K communities over 14 years, we measure how the macroscale community structure is organized with respect to age, gender, and U.S. political partisanship. Examining political content, we find Reddit underwent a significant polarization event around the 2016 U.S. presidential election, and remained highly polarized for years afterward. Contrary to conventional wisdom, however, individual-level polarization is rare; the system-level shift in 2016 was disproportionately driven by the arrival of new and newly political users. Political polarization on Reddit is unrelated to previous activity on the platform, and is instead temporally aligned with external events. We also observe a stark ideological asymmetry, with the sharp increase in 2016 being entirely attributable to changes in right-wing activity. Our methodology is broadly applicable to the study of online interaction, and our findings have implications for the design of online platforms, understanding the social contexts of online behaviour, and quantifying the dynamics and mechanisms of online polarization.


title: On Secondary Atomization and blockage of surrogate cough droplets in single and multi-layer face masks

abstract: By now it is well-understood that the usage of facemasks provides protection from transmission of viral loads through exhalation and inhalation of respiratory droplets. Therefore, during the current Covid-19 pandemic the usage of face masks is strongly recommended by health officials. Although three-layer masks are generally advised for usage, many commonly available or homemade masks contain only single and double layers. In this study, we show through detailed physics based analyses and high speed imaging that high momentum cough droplets on impingement on single- and double-layer masks can lead to significant partial penetration and more importantly atomization into numerous much smaller daughter droplets, thereby increasing the total population of the aerosol, which can remain suspended for a longer time. The possibility of secondary atomization of high momentum cough droplets due to impingement, hydrodynamic focusing and extrusion through the microscale pores in the fibrous network of the mask has not been explored before. However, this unique mode of aerosol generation poses a finite risk of infection as shown in this work. We also demonstrate that in single layer masks close to 70 % of a given droplet volume is atomized and only 30 % is trapped within the fibers. The entrapped volume is close to 90 % for double layer masks which still allows some atomization into smaller droplets. We however found that a triple-layer surgical mask permits negligible penetration and hence should be effective in preventing disease transmission.


title: The Lyman Continuum Escape Survey: Connecting Time-Dependent [OIII] and [OII] Line Emission with Lyman Continuum Escape Fraction in Simulations of Galaxy Formation

abstract: Escaping Lyman continuum photons from galaxies likely reionized the intergalactic medium at redshifts $z\gtrsim6$. However, the Lyman continuum is not directly observable at these redshifts and secondary indicators of Lyman continuum escape must be used to estimate the budget of ionizing photons. Observationally, at redshifts $z\sim2-3$ where the Lyman continuum is observationally accessible, surveys have established that many objects that show appreciable Lyman continuum escape fractions $f_{esc}$ also show enhanced [OIII]/[OII] (O$_{32}$) emission line ratios. Here, we use radiative transfer analyses of cosmological zoom-in simulations of galaxy formation to study the physical connection between $f_{esc}$ and O$_{32}$. Like the observations, we find that the largest $f_{esc}$ values occur at elevated O$_{32}\sim3-10$ and that the combination of high $f_{esc}$ and low O$_{32}$ is extremely rare. While high $f_{esc}$ and O$_{32}$ often are observable concurrently, the timescales of the physical origin for the processes are very different. Large O$_{32}$ values fluctuate on short ($\sim$1 Myr) timescales during the Wolf-Rayet-powered phase after the formation of star clusters, while channels of low absorption are established over tens of megayears by collections of supernovae. We find that while there is no direct causal relation between $f_{esc}$ and O$_{32}$, high $f_{esc}$ most often occurs after continuous input from star formation-related feedback events that have corresponding excursions to large O$_{32}$ emission. These calculations are in agreement with interpretations of observations that large $f_{esc}$ tends to occur when O$_{32}$ is large, but large O$_{32}$ does not necessarily imply efficient Lyman continuum escape.


title: Static response and Love numbers of Schwarzschild black holes

abstract: We derive the quadratic action for the physical degrees of freedom of massless spin-0, spin-1, and spin-2 perturbations on a Schwarzschild--(A)dS background in arbitrary dimensions. We then use these results to compute the static response of asymptotically flat Schwarzschild black holes to external fields. Our analysis reproduces known facts about black hole Love numbers, in particular that they vanish for all types of perturbation in four spacetime dimensions, but also leads to new results. For instance, we find that neutral Schwarzschild black holes polarize in the presence of an electromagnetic background in any number of spacetime dimensions except four. Moreover, we calculate for the first time black hole Love numbers for vector-type gravitational perturbations in higher dimensions and find that they generically do not vanish. Along the way, we shed some light on an apparent discrepancy between previous results in the literature, and clarify some aspects of the matching between perturbative calculations of static response on a Schwarzschild background and the point-particle effective theory


title: Space Telescope and Optical Reverberation Mapping Project. XII. Broad-Line Region Modeling of NGC 5548

abstract: We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$\beta$, C IV, and Ly$\alpha$ broad emission lines. We find an extended disk-like H$\beta$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$\alpha$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$\alpha$ emission arising at smaller radii than the H$\beta$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$\beta$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole.


title: Black hole perturbations of massive and partially massless spin-2 fields in (anti) de Sitter spacetime

abstract: We provide a systematic and comprehensive derivation of the linearized dynamics of massive and partially massless spin-2 particles in a Schwarzschild (anti) de Sitter black hole background, in four and higher spacetime dimensions. In particular, we show how to obtain the quadratic actions for the propagating modes and recast the resulting equations of motion in a Schr\"odinger-like form. In the case of partially massless fields in Schwarzschild de Sitter spacetime, we study the isospectrality between modes of different parity. In particular, we prove isospectrality analytically for modes with multipole number $L=1$ in four spacetime dimensions, providing the explicit form of the underlying symmetry. We show that isospectrality between partially massless modes of different parity is broken in higher-dimensional Schwarzschild de Sitter spacetimes.


title: Matter Power Spectrum Emulator for f(R) Modified Gravity Cosmologies

abstract: Testing a subset of viable cosmological models beyond General Relativity (GR), with implications for cosmic acceleration and the Dark Energy associated with it, is within the reach of Rubin Observatory Legacy Survey of Space and Time (LSST) and a part of its endeavor. Deviations from GR-w(z)CDM models can manifest in the growth rate of structure and lensing, as well as in screening effects on non-linear scales. We explore the constraining power of small-scale deviations predicted by the f(R) Hu-Sawicki Modified Gravity (MG) candidate, by emulating this model with COLA (COmoving Lagrangian Acceleration) simulations. We present the experimental design, data generation, and interpolation schemes in cosmological parameters and across redshifts for the emulation of the boost in the power spectra due to Modified Gravity effects. Three preliminary applications of the emulator highlight the sensitivity to cosmological parameters, Fisher forecasting and Markov Chain Monte Carlo inference for a fiducial cosmology. This emulator will play an important role for future cosmological analysis handling the formidable amount of data expected from Rubin Observatory LSST.


title: The Electro-Weak Phase Transition at Colliders: Confronting Theoretical Uncertainties and Complementary Channels

abstract: We explore and contrast the capabilities of future colliders to probe the nature of the electro-weak phase transition. We focus on the real singlet scalar field extension of the Standard Model, representing the most minimal, yet most elusive, framework that can enable a strong first-order electro-weak phase transition. By taking into account the theoretical uncertainties and employing the powerful complementarity between gauge and Higgs boson pair channels in the searches for new scalar particles, we find that a 100 TeV proton collider has the potential to confirm or falsify a strong first-order transition. Our results hint towards this occurring relatively early in its lifetime. Furthermore, by extrapolating down to 27 TeV, we find that a lower-energy collider may also probe a large fraction of the parameter space, if not all. Such early discoveries would allow for precise measurements of the new phenomena to be obtained at future colliders and would pave the way to definitively verify whether this is indeed the physical remnant of a scalar field that catalyses a strong first-order transition.


title: Magnetic Topological Quantum Chemistry

abstract: For over 100 years, the group-theoretic characterization of crystalline solids has provided the foundational language for diverse problems in physics and chemistry. However, the group theory of crystals with commensurate magnetic order has remained incomplete for the past 70 years, due to the complicated symmetries of magnetic crystals. In this work, we complete the 100-year-old problem of crystalline group theory by deriving the small corepresentations, momentum stars, compatibility relations, and magnetic elementary band corepresentations of the 1,421 magnetic space groups (MSGs), which we have made freely accessible through tools on the Bilbao Crystallographic Server. We extend Topological Quantum Chemistry to the MSGs to form a complete, real-space theory of band topology in magnetic and nonmagnetic crystalline solids - Magnetic Topological Quantum Chemistry (MTQC). Using MTQC, we derive the complete set of symmetry-based indicators of electronic band topology, for which we identify symmetry-respecting bulk and anomalous surface and hinge states.

categories: cond-mat.mes-hall cond-mat.str-el math-ph math.MP


title: Diabolical touching point in the magnetic energy levels of topological nodal-line metals

abstract: For three-dimensional metals, Landau levels disperse as a function of the magnetic field and the momentum wavenumber parallel to the field. In this two-dimensional parameter space, it is shown that two conically-dispersing Landau levels can touch at a diabolical point -- a Landau-Dirac point. The conditions giving rise to Landau-Dirac points are shown to be magnetic breakdown (field-induced quantum tunneling) and certain crystallographic spacetime symmetry. Both conditions are realizable in topological nodal-line metals, as we exemplify with CaP$_3$. A Landau-Dirac point reveals itself in anomalous batman-like peaks in the magnetoresistance, as well as in the onset of optical absorption linearly evolving to zero frequency as a function of the field magnitude/orientation.


title: #Election2020: The First Public Twitter Dataset on the 2020 US Presidential Election

abstract: The integrity of democratic political discourse is at the core to guarantee free and fair elections. With social media often dictating the tones and trends of politics-related discussion, it is of paramount important to be able to study online chatter, especially in the run up to important voting events, like in the case of the upcoming November 3, 2020 U.S. Presidential Election. Limited access to social media data is often the first barrier to impede, hinder, or slow down progress, and ultimately our understanding of online political discourse. To mitigate this issue and try to empower the Computational Social Science research community, we decided to publicly release a massive-scale, longitudinal dataset of U.S. politics- and election-related tweets. This multilingual dataset that we have been collecting for over one year encompasses hundreds of millions of tweets and tracks all salient U.S. politics trends, actors, and events between 2019 and 2020. It predates and spans the whole period of Republican and Democratic primaries, with real-time tracking of all presidential contenders of both sides of the isle. After that, it focuses on presidential and vice-presidential candidates. Our dataset release is curated, documented and will be constantly updated on a weekly-basis, until the November 3, 2020 election and beyond. We hope that the academic community, computational journalists, and research practitioners alike will all take advantage of our dataset to study relevant scientific and social issues, including problems like misinformation, information manipulation, interference, and distortion of online political discourse that have been prevalent in the context of recent election events in the United States and worldwide. Our dataset is available at: https://github.com/echen102/us-pres-elections-2020


title: Direct Detection of Atomic Dark Matter in White Dwarfs

abstract: Dark matter could have a dissipative asymmetric subcomponent in the form of atomic dark matter (aDM). This arises in many scenarios of dark complexity, and is a prediction of neutral naturalness, such as the Mirror Twin Higgs model. We show for the first time how White Dwarf cooling provides strong bounds on aDM. In the presence of a small kinetic mixing between the dark and SM photon, stars are expected to accumulate atomic dark matter in their cores, which then radiates away energy in the form of dark photons. In the case of white dwarfs, this energy loss can have a detectable impact on their cooling rate. We use measurements of the white dwarf luminosity function to tightly constrain the kinetic mixing parameter between the dark and visible photons, across many orders of magnitude in DM mass, down to values of $\epsilon \sim 10^{-12}$. Using this method we can constrain scenarios in which aDM constitutes fractions as small as $10^{-3}$ of the total dark matter density. Our methods are highly complementary to other methods of probing aDM, especially in scenarios where the aDM is arranged in a dark disk, which can make direct detection extremely difficult but actually slightly enhances our cooling constraints.


title: HOLISMOKES -- IV. Efficient mass modeling of strong lenses through deep learning

abstract: Modelling the mass distributions of strong gravitational lenses is often necessary to use them as astrophysical and cosmological probes. With the high number of lens systems ($>10^5$) expected from upcoming surveys, it is timely to explore efficient modeling approaches beyond traditional MCMC techniques that are time consuming. We train a CNN on images of galaxy-scale lenses to predict the parameters of the SIE mass model ($x,y,e_x,e_y$, and $\theta_E$). To train the network, we simulate images based on real observations from the HSC Survey for the lens galaxies and from the HUDF as lensed galaxies. We tested different network architectures, the effect of different data sets, and using different input distributions of $\theta_E$. We find that the CNN performs well and obtain with the network trained with a uniform distribution of $\theta_E$ $>0.5"$ the following median values with $1\sigma$ scatter: $\Delta x=(0.00^{+0.30}_{-0.30})"$, $\Delta y=(0.00^{+0.30}_{-0.29})" $, $\Delta \theta_E=(0.07^{+0.29}_{-0.12})"$, $\Delta e_x = -0.01^{+0.08}_{-0.09}$ and $\Delta e_y = 0.00^{+0.08}_{-0.09}$. The bias in $\theta_E$ is driven by systems with small $\theta_E$. Therefore, when we further predict the multiple lensed image positions and time delays based on the network output, we apply the network to the sample limited to $\theta_E>0.8"$. In this case, the offset between the predicted and input lensed image positions is $(0.00_{-0.29}^{+0.29})"$ and $(0.00_{-0.31}^{+0.32})"$ for $x$ and $y$, respectively. For the fractional difference between the predicted and true time delay, we obtain $0.04_{-0.05}^{+0.27}$. Our CNN is able to predict the SIE parameters in fractions of a second on a single CPU and with the output we can predict the image positions and time delays in an automated way, such that we are able to process efficiently the huge amount of expected lens detections in the near future.


title: Evolution of the density PDF in star forming clouds: the role of gravity

abstract: We derive an analytical theory of the PDF of density fluctuations in supersonic turbulence in the presence of gravity in star-forming clouds. The theory is based on a rigorous derivation of a combination of the Navier-Stokes continuity equations for the fluid motions and the Poisson equation for the gravity. It extends upon previous approaches first by including gravity, second by considering the PDF as a dynamical system, not a stationary one. We derive the transport equations of the density PDF, characterize its evolution and determine the density threshold above which gravity strongly affects and eventually dominates the dynamics of turbulence. We demonstrate the occurence of {\it two} power law tails in the PDF, with two characteristic exponents, corresponding to two different stages in the balance between turbulence and gravity. Another important result of this study is to provide a procedure to relate the observed {\it column density} PDFs to the corresponding {\it volume density} PDFs. This allows to infer, from the observation of column-densities, various physical parameters characterizing molecular clouds, notably the virial parameter. Furthermore, the theory offers the possibility to date the clouds in units of ${t}_{\rm coll}$, the time since a statistically significant fraction of the cloud started to collapse. The theoretical results and diagnostics reproduce very well numerical simulations and observations of star-forming clouds. The theory provides a sound theoretical foundation and quantitative diagnostics to analyze observations or numerical simulations of star-forming regions and to characterize the evolution of the density PDF in various regions of molecular clouds.


title: EvryFlare III: Temperature Evolution and Habitability Impacts of Dozens of Superflares Observed Simultaneously by Evryscope and TESS

abstract: Superflares may provide the dominant source of biologically relevant UV radiation to rocky habitable zone M-dwarf planets (M-Earths), altering planetary atmospheres and conditions for surface life. The combined line and continuum flare emission has usually been approximated by a 9000 K blackbody. If superflares are hotter, then the UV emission may be 10X higher than predicted from the optical. However, it is unknown for how long M-dwarf superflares reach temperatures above 9000 K. Only a handful of M-dwarf superflares have been recorded with multi-wavelength high-cadence observations. We double the total number of events in the literature using simultaneous Evryscope and TESS observations to provide the first systematic exploration of the temperature evolution of M-dwarf superflares. We also increase the number of superflaring M-dwarfs with published time-resolved blackbody evolution by ~10X. We measure temperatures at 2 min cadence for 42 superflares from 27 K5-M5 dwarfs. We find superflare peak temperatures (defined as the mean of temperatures corresponding to flare FWHM) increase with flare energy and impulse. We find the amount of time flares emit at temperatures above 14,000 K depends on energy. We discover 43% of the flares emit above 14,000 K, 23% emit above 20,000 K and 5% emit above 30,000 K. The largest and hottest flare briefly reached 42,000 K. Some do not reach 14,000 K. During superflares, we estimate M-Earths orbiting <200 Myr stars typically receive a top-of-atmosphere UV-C flux of ~120 W m^-2 and up to 10^3 W m^-2, 100-1000X the time-averaged XUV flux from Proxima Cen.


title: Light charged fermions in quantum gravity

abstract: We study the impact of quantum gravity on a system of chiral fermions that are charged under an Abelian gauge group. Under the impact of quantum gravity, a finite value of the gauge coupling could be generated and in turn drive four-fermion interactions to criticality. We find indications that the gravity-gauge-fermion interplay protects the lightness of fermions for a large enough number of fermions. On the other hand, for a smaller number of fermions, chiral symmetry may be broken, which would be in tension with the observation of light fermions.


title: The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. III. Substructures in Protostellar Disks

abstract: The prevalence of substructures in $\sim1-10$ Myr old protoplanetary disks, which are often linked to planet formation, has raised the question of how early such features form, and as a corollary, how early planet formation begins. Here we present observations of seven protostellar disks (aged $\sim0.1-1$ Myr) from the VLA/ALMA Nascent Disk and Multiplicity Survey of Orion Protostars (VANDAM: Orion) that show clear substructures, thereby demonstrating that these features can form early in the lifetimes of disks. We use simple analytic models as well as detailed radiative transfer modeling to characterize their structure. In particular we show that at least four of the sources have relatively massive envelopes, indicating that they are particularly young, likely the youngest disks with substructures known to-date. Several of these disks also have emission from an inner disk that is offset from the center of the ring structure. Given the size of the cleared out regions of the disk, it is unclear, however, whether these features are related to planet formation, or rather if they are signposts of close-separation binary formation at early times.


title: Gauging the Kitaev chain

abstract: We gauge the fermion parity symmetry of the Kitaev chain. While the bulk of the model becomes an Ising chain of gauge-invariant spins in a tilted field, near the boundaries the global fermion parity symmetry survives gauging, leading to local gauge-invariant Majorana operators. In the absence of vortices, the Higgs phase exhibits fermionic symmetry-protected topological (SPT) order distinct from the Kitaev chain. Moreover, the deconfined phase can be stable even in the presence of vortices. We also undertake a comprehensive study of a gently gauged model which interpolates between the ordinary and gauged Kitaev chains. This showcases rich quantum criticality and illuminates the topological nature of the Higgs phase. Even in the absence of superconducting terms, gauging leads to an SPT phase which is intrinsically gapless due to an emergent anomaly.

categories: cond-mat.str-el cond-mat.quant-gas hep-lat hep-th quant-ph


title: Displaced Vertex signatures of a pseudo-Goldstone sterile neutrino

abstract: Low-scale models of neutrino mass generation often feature sterile neutrinos with masses in the GeV-TeV range, which can be produced at colliders through their mixing with the Standard Model neutrinos. We consider an alternative scenario in which the sterile neutrino is produced in the decay of a heavier particle, such that its production cross section does not depend on the active-sterile neutrino mixing angles. The mixing angles can be accessed through the decays of the sterile neutrino, provided that they lead to observable displaced vertices. We present an explicit realization of this scenario in which the sterile neutrino is the supersymmetric partner of a pseudo-Nambu-Goldstone boson, and is produced in the decays of higgsino-like neutralinos and charginos. The model predicts the active-sterile neutrino mixing angles in terms of a small number of parameters. We show that a sterile neutrino with a mass between a few 10 GeV and 200 GeV can lead to observable displaced vertices at the LHC, and outline a strategy for reconstructing experimentally its mixing angles.


title: Biases and Cosmic Variance in Molecular Gas Abundance Measurements

abstract: Recent deep millimeter-wave surveys attempt to measure the carbon monoxide (CO) luminosity function and mean molecular gas density through blind detections of CO emission lines. While the cosmic star formation rate density is now constrained in fields hundreds of arcmin$^2$ or more, molecular gas studies have been limited to $\le50$ arcmin$^2$. These small fields result in significant biases that have not been accounted for in published results. To quantify these biases, we assign CO luminosities to halos in cosmological simulations to produce mock observations for a range of field sizes. We find that fields $\lesssim 10$ arcmin$^2$ alter the recovered shape of the luminosity function, causing underestimates of the number of bright objects. Our models suggest that current surveys are sensitive enough to detect sources responsible for approximately half of the cosmic molecular gas density at high redshift. However, uncertainties in the gas density measurement are large, and cosmic variance may double the uncertainty claimed in these surveys. As a result, the field size needed to detect redshift evolution in the molecular gas at high confidence may be more than an order of magnitude larger than what current surveys have achieved. Shot power intensity mapping measurements are particularly sensitive to Poisson variance and require yet larger areas to constrain the gas density or its evolution. We provide a simple prescription for approximating uncertainty in total CO emission as a function of survey area and redshift, for both direct detection and intensity mapping surveys.


title: Experimental Evidence that Zn Impurities Pin Pair-Density-Wave Order in La$_{2-x}$Ba$_x$CuO$_4$

abstract: Both Zn-doping and $c$-axis magnetic fields have been observed to increase the spin stripe order in La$_{2-x}$Ba$_x$CuO$_4$ with $x$ close to 1/8. For $x=0.095$, the applied magnetic field also causes superconducting layers to decouple, presumably by favoring pair-density-wave order that consequently frustrates interlayer Josephson coupling. Here we show that introducing 1% Zn also leads to an initial onset of two-dimensional (2D) superconductivity, followed by 3D superconductivity at lower temperatures, even in zero field. We infer that the Zn pins pair-density-wave order locally, establishing the generality of such behavior.


title: Charge density waves and their transitions in anisotropic quantum Hall systems

abstract: In recent experiments, external anisotropy has been a useful tool to tune different phases and study their competitions. In this paper, we look at the quantum Hall charge density wave states in the $N=2$ Landau level. Without anisotropy, there are two first-order phase transitions between the Wigner crystal, the $2$-electron bubble phase, and the stripe phase. By adding mass anisotropy, our analytical and numerical studies show that the $2$-electron bubble phase disappears and the stripe phase significantly enlarges its domain in the phase diagram. Meanwhile, a regime of stripe crystals that may be observed experimentally is unveiled after the bubble phase gets out. Upon increase of the anisotropy, the energy of the phases at the transitions becomes progressively smooth as a function of the filling. We conclude that all first-order phase transitions are replaced by continuous phase transitions, providing a possible realisation of continuous quantum crystalline phase transitions.


title: Closure Traces: Novel Calibration-Insensitive Quantities for Radio Astronomy

abstract: Closure phases and closure amplitudes have proven critical to modern radio interferometry due to their insensitivity to the uncertain station gains. We present the first set of closure quantities constructed from parallel-hand and cross-hand visibilities that are insensitive to both station gains and to polarimetric leakage. These complex "closure traces" are a natural extension of closure amplitudes and closure phases, are independent of all station-based linear corruptions of the polarized visibilities, and are complete in the sense that they contain all remaining information present in the visibility data. Products of closure traces on so-called "conjugate" quadrangles are sensitive only to structure in the source polarization fraction -- independent of variations in the Stokes $I$ structure -- and thereby provide unambiguous probes of polarization in astronomical sources.


title: Unified treatment of synchronization patterns in generalized networks with higher-order, multilayer, and temporal interactions

abstract: When describing complex interconnected systems, one often has to go beyond the standard network description to account for generalized interactions. Here, we establish a unified framework to simplify the stability analysis of cluster synchronization patterns for a wide range of generalized networks, including hypergraphs, multilayer networks, and temporal networks. The framework is based on finding a simultaneous block diagonalization of the matrices encoding the synchronization pattern and the network topology. As an application, we use simultaneous block diagonalization to unveil an intriguing type of chimera states that appear only in the presence of higher-order interactions. The unified framework established here can be extended to other dynamical processes and can facilitate the discovery of emergent phenomena in complex systems with generalized interactions.


title: Self-calibration of weak lensing systematic effects using combined two- and three-point statistics

abstract: We investigate the prospects for using the weak lensing bispectrum alongside the power spectrum to control systematic uncertainties in a Euclid-like survey. Three systematic effects are considered: the intrinsic alignment of galaxies, uncertainties in the means of tomographic redshift distributions, and multiplicative bias in the measurement of the shear signal. We find that the bispectrum is very effective in mitigating these systematic errors. Varying all three systematics simultaneously, a joint power spectrum and bispectrum analysis reduces the area of credible regions for the cosmological parameters $\Omega_\mathrm{m}$ and $\sigma_8$ by a factor of 90 and for the two parameters of a time-varying dark energy equation of state by a factor of almost 20, compared with the baseline approach of using the power spectrum alone and of imposing priors consistent with the accuracy requirements specified for Euclid. We also demonstrate that including the bispectrum self-calibrates all three systematic effects to the stringent levels required by the forthcoming generation of weak lensing surveys, thereby reducing the need for external calibration data.


title: Compensating PDE actuator and sensor dynamics using Sylvester equation

abstract: We consider the problem of stabilizing PDE-ODE cascade systems in which the input is applied to the PDE system whose output drives the ODE system. We also consider the dual problem of constructing an observer for ODE-PDE cascade systems in which the output of the ODE system drives the PDE system, whose output is measured. The PDE in these problems is stable and the ODE is unstable. While the ODE system models the plant in both the problems, the PDE system models the actuator in the stabilization problem and the sensor in the dual problem. In the literature, these problems have been solved for specific PDE models using the backstepping approach. In contrast, in the present work we consider these problems in an abstract framework by letting the PDE system be any regular linear system. Using a state transformation obtained by solving a Sylvester equation with unbounded operators, we first diagonalize the state operator corresponding to the cascade systems. We then solve the stabilization problem and the dual estimation problem, provided they are solvable, by solving certain finite-dimensional counterparts. We also derive necessary and sufficient conditions for verifying the solvability of these problems. We show that the controller which solves the stabilization problem is robust to certain unbounded perturbations. We illustrate our theory by designing a stabilizing controller for a PDE-ODE cascade in which the PDE is a 1D diffusion equation and an observer for a ODE-PDE cascade in which the PDE is a 1D wave equation.


title: Competition of Spinon Fermi Surface and Heavy Fermi Liquids states from the Periodic Anderson to the Hubbard model

abstract: We study a model of correlated electrons coupled by tunnelling to a layer of itinerant metallic electrons, which allows to interpolate from a frustrated limit favorable to spin liquid states to a Kondo-lattice limit favorable to interlayer coherent heavy metallic states. We study the competition of the spinon fermi surface state and the interlayer coherent heavy Kondo metal that appears with increasing tunnelling. Employing a slave rotor mean-field approach, we obtain a phase diagram and describe two regimes where the spin liquid state is destroyed by weak interlayer tunnelling, (i) the Kondo limit in which the correlated electrons can be viewed as localized spin moments and (ii) near the Mott metal-insulator-transition where the spinon Fermi surface transitions continuously into a Fermi liquid. We study the shape of LDOS spectra of the putative spin liquid layer in the heavy Fermi liquid phase and describe the temperature dependence of its width arising from quasiparticle interactions and disorder effects throughout this phase diagram, in an effort to understand recent STM experiments of the candidate spin liquid 1T-TaSe$_2$ residing on metallic 1H-TaSe$_2$. Comparison of the shape and temperature dependence of the theoretical and experimental LDOS suggest that this system is either close to the localized Kondo limit, or in an intermediate coupling regime where the Kondo coupling and the Heisenberg exchange interaction are comparable.


title: On the rate of crustal failures in young magnetars

abstract: The activity of magnetars is powered by their intense and dynamic magnetic fields and has been proposed as the trigger to extragalactic Fast Radio Bursts. Here we estimate the frequency of crustal failures in young magnetars, by computing the magnetic stresses in detailed magneto-thermal simulations including Hall drift and Ohmic dissipation. The initial internal topology at birth is poorly known but is likely to be much more complex than a dipole. Thus, we explore a wide range of initial configurations, finding that the expected rate of crustal failures varies by orders of magnitude depending on the initial magnetic configuration. Our results show that this rate scales with the crustal magnetic energy, rather than with the often used surface value of the dipolar component related to the spin-down torque. The estimated frequency of crustal failures for a given dipolar component can vary by orders of magnitude for different initial conditions, depending on how much magnetic energy is distributed in the crustal non-dipolar components, likely dominant in newborn magnetars. The quantitative reliability of the expected event rate could be improved by a better treatment of the magnetic evolution in the core and the elastic/plastic crustal response, here not included. Regardless of that, our results are useful inputs in modelling the outburst rate of young Galactic magnetars, and their relation with the Fast Radio Bursts in our and other galaxies.


title: Composition determination of cosmic rays from the muon content of the showers

abstract: The origin and nature of ultra high energy cosmic rays remains being a mystery. However, great progress has been made in recent years due to the observations performed by the Pierre Auger Observatory and Telescope Array. In particular, it is believed that the composition information of the cosmic rays as a function of the energy can play a fundamental role for the understanding of their origin. The best indicators for primary mass composition are the muon content of extensive air shower and the atmospheric depth of the shower maximum. In this work we consider a maximum likelihood method to perform mass composition analyses based on the number of muons measured by underground muon detectors. The analyses are based on numerical simulations of the showers. The effects introduced by the detectors and the methods used to reconstruct the experimental data are also taken into account through a dedicated simulation that uses as input the information of the simulated showers. In order to illustrate the use of the method, we consider AMIGA (Auger Muons and Infill for the Ground Array), the low energy extension of the Pierre Auger Observatory that directly measures the muonic content of extensive air showers. We also study in detail the impact of the use of different high energy hadronic interaction models in the composition analyses performed. It is found that differences of a few percent between the predicted number of muons have a significant impact on composition determination.


title: The CAMELS project: Cosmology and Astrophysics with MachinE Learning Simulations

abstract: We present the Cosmology and Astrophysics with MachinE Learning Simulations --CAMELS-- project. CAMELS is a suite of 4,233 cosmological simulations of $(25~h^{-1}{\rm Mpc})^3$ volume each: 2,184 state-of-the-art (magneto-)hydrodynamic simulations run with the AREPO and GIZMO codes, employing the same baryonic subgrid physics as the IllustrisTNG and SIMBA simulations, and 2,049 N-body simulations. The goal of the CAMELS project is to provide theory predictions for different observables as a function of cosmology and astrophysics, and it is the largest suite of cosmological (magneto-)hydrodynamic simulations designed to train machine learning algorithms. CAMELS contains thousands of different cosmological and astrophysical models by way of varying $\Omega_m$, $\sigma_8$, and four parameters controlling stellar and AGN feedback, following the evolution of more than 100 billion particles and fluid elements over a combined volume of $(400~h^{-1}{\rm Mpc})^3$. We describe the simulations in detail and characterize the large range of conditions represented in terms of the matter power spectrum, cosmic star formation rate density, galaxy stellar mass function, halo baryon fractions, and several galaxy scaling relations. We show that the IllustrisTNG and SIMBA suites produce roughly similar distributions of galaxy properties over the full parameter space but significantly different halo baryon fractions and baryonic effects on the matter power spectrum. This emphasizes the need for marginalizing over baryonic effects to extract the maximum amount of information from cosmological surveys. We illustrate the unique potential of CAMELS using several machine learning applications, including non-linear interpolation, parameter estimation, symbolic regression, data generation with Generative Adversarial Networks (GANs), dimensionality reduction, and anomaly detection.


title: Energy-participation quantization of Josephson circuits

abstract: Superconducting microwave circuits incorporating nonlinear devices, such as Josephson junctions, are one of the leading platforms for emerging quantum technologies. Increasing circuit complexity further requires efficient methods for the calculation and optimization of the spectrum, nonlinear interactions, and dissipation in multi-mode distributed quantum circuits. Here, we present a method based on the energy-participation ratio (EPR) of a dissipative or nonlinear element in an electromagnetic mode. The EPR, a number between zero and one, quantifies how much of the energy of a mode is stored in each element. It obeys universal constraints--valid regardless of the circuit topology and nature of the nonlinear elements. The EPR of the elements are calculated from a unique, efficient electromagnetic eigenmode simulation of the linearized circuit, including lossy elements. Their set is the key input to the determination of the quantum Hamiltonian of the system. The method provides an intuitive and simple-to-use tool to quantize multi-junction circuits. It is especially well-suited for finding the Hamiltonian and dissipative parameters of weakly anharmonic systems, such as transmon qubits coupled to resonators, or Josephson transmission lines. We experimentally tested this method on a variety of Josephson circuits, and demonstrated agreement within several percents for nonlinear couplings and modal Hamiltonian parameters, spanning five-orders of magnitude in energy, across a dozen samples.


title: Nonsmooth Exact Penalization Second-order Methods for incompressible Bingham flows

abstract: We consider the exact penalization of the incompressibility condition $div(u)=0$ for the velocity field of a Bingham fluid in terms of the $L^1$-norm. This penalization procedure results in a nonsmooth optimization problem for which we propose an algorithm using generalized second-order information. Our method solves the resulting nonsmooth problem by considering the steepest descent direction and extra generalized second-order information associated to the nonsmooth term. This method has the advantage that the divergence-free property is enforced by the descent direction proposed by the method without the need of build-in divergence-free approximation schemes. The inexact penalization approach, given by the $L^2$-norm, is also considered in our discussion and comparison.


title: Pix2Prof: fast extraction of sequential information from galaxy imagery via a deep natural language 'captioning' model

abstract: We present 'Pix2Prof', a deep learning model that can eliminate any manual steps taken when extracting galaxy profiles. We argue that a galaxy profile of any sort is conceptually similar to a natural language image caption. This idea allows us to leverage image captioning methods from the field of natural language processing, and so we design Pix2Prof as a float sequence 'captioning' model suitable for galaxy profile inference. We demonstrate the technique by approximating a galaxy surface brightness (SB) profile fitting method that contains several manual steps. Pix2Prof processes $\sim$1 image per second on an Intel Xeon E5 2650 v3 CPU, improving on the speed of the manual interactive method by more than two orders of magnitude. Crucially, Pix2Prof requires no manual interaction, and since galaxy profile estimation is an embarrassingly parallel problem, we can further increase the throughput by running many Pix2Prof instances simultaneously. In perspective, Pix2Prof would take under an hour to infer profiles for $10^5$ galaxies on a single NVIDIA DGX-2 system. A single human expert would take approximately two years to complete the same task. Automated methodology such as this will accelerate the analysis of the next generation of large area sky surveys expected to yield hundreds of millions of targets. In such instances, all manual approaches -- even those involving a large number of experts -- will be impractical.


title: 'Interaction-Free' Channel Discrimination

abstract: In this work, we investigate the question, which objects one can discriminate perfectly by 'interaction-free' measurements. To this end, we interpret the Elitzur-Vaidman bomb-tester experiment as a quantum channel discrimination problem and generalize the notion of 'interaction-free' measurement to arbitrary quantum channels. Our main result is a necessary and sufficient criterion for when it is possible or impossible to discriminate quantum channels in an 'interaction-free' manner (i.e., such that the discrimination error probability and the 'interaction' probability can be made arbitrarily small). For the case where our condition holds, we devise an explicit protocol with the property that both probabilities approach zero with an increasing number of channel uses, $N$. More specifically, the 'interaction' probability in our protocol decays as $\frac{1}{N}$ and we show that this rate is the optimal achievable one. Furthermore, our protocol only needs at most one ancillary qubit and might thus be implementable in near-term experiments. For the case where our condition does not hold, we prove an inequality that quantifies the trade-off between the error probability and the 'interaction' probability.