title: Modularity, Code Specialization, and Zero-Cost Abstractions for Program Verification

abstract: For all the successes in verifying low-level, efficient, security-critical code, little has been said or studied about the structure, architecture and engineering of such large-scale proof developments. We present the design, implementation and evaluation of a set of language-based techniques that allow the programmer to modularly write and verify code at a high level of abstraction, while retaining control over the compilation process and producing high-quality, zero-overhead, low-level code suitable for integration into mainstream software. We implement our techniques within the F* proof assistant, and specifically its shallowly-embedded Low* toolchain that compiles to C. Through our evaluation, we establish that our techniques were critical in scaling the popular HACL* library past 100,000 lines of verified source code, and brought about significant gains in proof engineer productivity. The exposition of our methodology converges on one final, novel case study: the streaming API, a finicky API that has historically caused many bugs in high-profile software. Using our approach, we manage to capture the streaming semantics in a generic way, and apply it ``for free'' to over a dozen use-cases. Six of those have made it into the reference implementation of the Python programming language, replacing the previous CVE-ridden code.


title: Generating images from caption and vice versa via CLIP-Guided Generative Latent Space Search

abstract: In this research work we present CLIP-GLaSS, a novel zero-shot framework to generate an image (or a caption) corresponding to a given caption (or image). CLIP-GLaSS is based on the CLIP neural network, which, given an image and a descriptive caption, provides similar embeddings. Differently, CLIP-GLaSS takes a caption (or an image) as an input, and generates the image (or the caption) whose CLIP embedding is the most similar to the input one. This optimal image (or caption) is produced via a generative network, after an exploration by a genetic algorithm. Promising results are shown, based on the experimentation of the image Generators BigGAN and StyleGAN2, and of the text Generator GPT2


title: Online Learning with Simple Predictors and a Combinatorial Characterization of Minimax in 0/1 Games

abstract: Which classes can be learned properly in the online model? -- that is, by an algorithm that at each round uses a predictor from the concept class. While there are simple and natural cases where improper learning is necessary, it is natural to ask how complex must the improper predictors be in such cases. Can one always achieve nearly optimal mistake/regret bounds using "simple" predictors? In this work, we give a complete characterization of when this is possible, thus settling an open problem which has been studied since the pioneering works of Angluin (1987) and Littlestone (1988). More precisely, given any concept class C and any hypothesis class H, we provide nearly tight bounds (up to a log factor) on the optimal mistake bounds for online learning C using predictors from H. Our bound yields an exponential improvement over the previously best known bound by Chase and Freitag (2020). As applications, we give constructive proofs showing that (i) in the realizable setting, a near-optimal mistake bound (up to a constant factor) can be attained by a sparse majority-vote of proper predictors, and (ii) in the agnostic setting, a near-optimal regret bound (up to a log factor) can be attained by a randomized proper algorithm. A technical ingredient of our proof which may be of independent interest is a generalization of the celebrated Minimax Theorem (von Neumann, 1928) for binary zero-sum games. A simple game which fails to satisfy Minimax is "Guess the Larger Number", where each player picks a number and the larger number wins. The payoff matrix is infinite triangular. We show this is the only obstruction: if a game does not contain triangular submatrices of unbounded sizes then the Minimax Theorem holds. This generalizes von Neumann's Minimax Theorem by removing requirements of finiteness (or compactness), and captures precisely the games of interest in online learning.


title: A Novel Use of Discrete Wavelet Transform Features in the Prediction of Epileptic Seizures from EEG Data

abstract: This paper demonstrates the predictive superiority of discrete wavelet transform (DWT) over previously used methods of feature extraction in the diagnosis of epileptic seizures from EEG data. Classification accuracy, specificity, and sensitivity are used as evaluation metrics. We specifically show the immense potential of 2 combinations (DWT-db4 combined with SVM and DWT-db2 combined with RF) as compared to others when it comes to diagnosing epileptic seizures either in the balanced or the imbalanced dataset. The results also highlight that MFCC performs less than all the DWT used in this study and that, The mean-differences are statistically significant respectively in the imbalanced and balanced dataset. Finally, either in the balanced or the imbalanced dataset, the feature extraction techniques, the models, and the interaction between them have a statistically significant effect on the classification accuracy.


title: The Privatization of AI Research(-ers): Causes and Potential Consequences -- From university-industry interaction to public research brain-drain?

abstract: The private sector is playing an increasingly important role in basic Artificial Intelligence (AI) R&D. This phenomenon, which is reflected in the perception of a brain drain of researchers from academia to industry, is raising concerns about a privatisation of AI research which could constrain its societal benefits. We contribute to the evidence base by quantifying transition flows between industry and academia and studying its drivers and potential consequences. We find a growing net flow of researchers from academia to industry, particularly from elite institutions into technology companies such as Google, Microsoft and Facebook. Our survival regression analysis reveals that researchers working in the field of deep learning as well as those with higher average impact are more likely to transition into industry. A difference-in-differences analysis of the effect of switching into industry on a researcher's influence proxied by citations indicates that an initial increase in impact declines as researchers spend more time in industry. This points at a privatisation of AI knowledge compared to a counterfactual where those high-impact researchers had remained in academia. Our findings highlight the importance of strengthening the public AI research sphere in order to ensure that the future of this powerful technology is not dominated by private interests.


title: Heterogeneous Graph based Deep Learning for Biomedical Network Link Prediction

abstract: Multi-scale biomedical knowledge networks are expanding with emerging experimental technologies that generates multi-scale biomedical big data. Link prediction is increasingly used especially in bipartite biomedical networks to identify hidden biological interactions and relationshipts between key entities such as compounds, targets, gene and diseases. We propose a Graph Neural Networks (GNN) method, namely Graph Pair based Link Prediction model (GPLP), for predicting biomedical network links simply based on their topological interaction information. In GPLP, 1-hop subgraphs extracted from known network interaction matrix is learnt to predict missing links. To evaluate our method, three heterogeneous biomedical networks were used, i.e. Drug-Target Interaction network (DTI), Compound-Protein Interaction network (CPI) from NIH Tox21, and Compound-Virus Inhibition network (CVI). Our proposed GPLP method significantly outperforms over the state-of-the-art baselines. In addition, different network incompleteness is analysed with our devised protocol, and we also design an effective approach to improve the model robustness towards incomplete networks. Our method demonstrates the potential applications in other biomedical networks.


title: Combined theoretical study of the $D^+ \to \pi^+ \eta \eta$ and $D^+ \to \pi^+ \pi^0 \eta $ reactions

abstract: We study the $D^+ \to \pi^+ \eta \eta $ and $D^+ \to \pi^+ \pi^0 \eta$ reactions, which are single Cabibbo suppressed and can proceed both through internal and external emission. The primary mechanisms at quark level are considered, followed by hadronization to produce three mesons in the $D^+$ decay, and after that the final state interaction of these mesons leads to the production of the $a_0(980)$ resonance, seen in the $\pi^+ \eta$, $\pi^0 \eta$ mass distributions. The theory has three unknown parameters to determine the shape of the distributions and the ratio between the $D^+ \to \pi^+ \eta \eta$ and $D^+ \to \pi^+ \pi^0 \eta$ rates. This ratio restricts much the sets of parameters but there is still much freedom leading to different shapes in the mass distributions. We call for a measurement of these mass distributions that will settle the reaction mechanism, while at the same time provide relevant information on the way that the $a_0(980)$ resonance is produced in the reactions.


title: Insights into ventilation hysteresis shift due to flow unsteadiness in ventilated supercavitation

abstract: Understanding ventilation strategy of a supercavity is important for designing high-speed underwater vehicles wherein an artificial gas pocket is created behind a flow separation device for drag reduction. Our study investigates the effect of flow unsteadiness on the ventilation requirements to form (CQf) and collapse (CQc) a supercavity. Imposing flow unsteadiness on the incoming flow has shown an increment in higher CQf at low free stream velocity and lower CQf at high free stream velocity. High-speed imaging reveals distinctly different behaviors in the recirculation region for low and high freestream velocity under unsteady flows. At low free stream velocities, the recirculation region formed downstream of a cavitator shifted vertically with flow unsteadiness, resulting in lower bubble collision and coalescence probability, which is critical for the supercavity formation process. The recirculation region negligibly changed with flow unsteadiness at high free stream velocity and less ventilation is required to form a supercavity compared to that of the steady incoming flow. Such a difference is attributed to the increased transverse Reynolds stress that aids bubble collision in a confined space of the recirculation region. CQc is found to heavily rely on the vertical component of the flow unsteadiness and the free stream velocity. Interfacial instability located upper rear of the supercavity develops noticeably with flow unsteadiness and additional bubbles formed by the distorted interface shed from the supercavity, resulting in an increased CQc. Further analysis on the quantification of such additional bubble leakage rate indicates that the development and amplitude of the interfacial instability accounts for the variation of CQc under a wide range of flow unsteadiness. Our study provides some insights on the design of a ventilation strategy for supercavitating vehicles in practice.


title: The conforming virtual element method for polyharmonic and elastodynamics problems: a review

abstract: In this paper, we review recent results on the conforming virtual element approximation of polyharmonic and elastodynamics problems. The structure and the content of this review is motivated by three paradigmatic examples of applications: classical and anisotropic Cahn-Hilliard equation and phase field models for brittle fracture, that are briefly discussed in the first part of the paper. We present and discuss the mathematical details of the conforming virtual element approximation of linear polyharmonic problems, the classical Cahn-Hilliard equation and linear elastodynamics problems.


title: A late times approximation for the SYK spectral form factor

abstract: We find a late times approximation for the SYK spectral form factor from a large $N$ steepest descent version of the path integral over two replica collective fields. Main ingredients are a suitable uv regularization of the two replica kinetic operator, the property of its Fourier transform and some spectral analysis of the four point function two replica ladder kernel.

categories: hep-th cond-mat.dis-nn cond-mat.str-el quant-ph


title: Internal-wave billiards in trapezoids and similar tables

abstract: We call internal-wave billiard the dynamical system of a point particle that moves freely inside a planar domain (the table) and is reflected by its boundary according to this rule: reflections are standard Fresnel reflections but with the pretense that the boundary at any collision point is either horizontal or vertical (relative to a predetermined direction representing gravity). These systems are point particle approximations for the motion of internal gravity waves in closed containers, hence the name. For a class of tables similar to rectangular trapezoids, but with the slanted leg replaced by a general curve with downward concavity, we prove that the dynamics has only three asymptotic regimes: (1) there exist a global attractor and a global repellor, which are periodic and might coincide; (2) there exists a beam of periodic trajectories, whose boundary (if any) comprises an attractor and a repellor for all the other trajectories; (3) all trajectories are dense (that is, the system is minimal). Furthermore, in the prominent case where the table is an actual trapezoid, we study the sets in parameter space relative to the three regimes. We prove in particular that the set for (1) has positive measure (giving a rigorous proof of the existence of Arnol'd tongues for internal-wave billiards), whereas the sets for (2) and (3) are non-empty but have measure zero.


title: Low-energy decomposition results over finite fields

abstract: We prove various low-energy decomposition results, showing that we can decompose a finite set $A\subset \mathbb{F}_p$ satisfying $|A|<p^{5/8}$, into $A = S\sqcup T$ so that, for a non-degenerate quadratic $f\in \mathbb{F}_p[x,y]$, we have \[ |\{(s_1,s_2,s_3,s_4)\in S^4 : s_1 + s_2 = s_3 + s_4\}| \ll |A|^{3 - \frac15 + \varepsilon} \] and \[ |\{(t_1,t_2,t_3,t_4)\in T^4 : f(t_1, t_2) = f(t_3, t_4)\}|\ll |A|^{3 - \frac15 + \varepsilon}\,. \] Variations include extending this result to large $A$ and a low-energy decomposition involving additive energy of images of rational functions. This gives a quantitative improvement to a result of Roche-Newton, Shparlinski and Winterhof as well as a generalisation of a result of Rudnev, Shkredov and Stevens. We consider applications to conditional expanders, exponential sum estimates and the finite field Littlewood problem. In particular, we improve results of Mirzaei, Swaenepoel and Winterhof and Garcia.


title: Calibration of the Surface Array of the Pierre Auger Observatory

abstract: The Pierre Auger Observatory is designed to study cosmic rays of the highest energies ($>10^{19}$ eV). The ground array of the Observatory will consist of 1600 water Cherenkov detectors deployed over 3000 km^2. The remoteness and large number of detectors require a robust, automatic self-calibration procedure. It relies on the measurement of the average charge collected by a photomultiplier tube from the Cherenkov light produced by a vertical and central through-going muon determined to 5 - 10% at the detector via a novel rate-based technique and to 3% precision through analysis of histograms of the charge distribution. The parameters needed for the calibration are measured every minute, allowing for an accurate determination of the signals recorded from extensive air showers produced by primary cosmic rays. The method also enables stable and uniform triggering conditions to be achieved.


title: Construction of Nahm data and BPS monopoles with continuous symmetries

abstract: We study solutions to Nahm's equations with continuous symmetries and, under certain (mild) hypotheses, we classify the corresponding Ans\"atze. Using our classification, we construct novel Nahm data, and prescribe methods for generating further solutions. Finally, we use these results to construct new BPS monopoles with spherical symmetry.


title: Johnson-Nyquist Noise Effects in Neutron Electric-Dipole-Moment Experiments

abstract: Magnetic Johnson-Nyquist noise (JNN) originating from metal electrodes, used to create a static electric field in neutron electric-dipole-moment (nEDM) experiments, may limit the sensitivity of measurements. We present here the first dedicated study on JNN applied to a large-scale long-measurement-time experiment with the implementation of a co-magnetometry. In this study, we derive surface- and volume-averaged root-mean-square normal noise amplitudes at a certain frequency bandwidth for a cylindrical geometry. In addition, we model the source of noise as a finite number of current dipoles and demonstrate a method to simulate temporal and three-dimensional spatial dependencies of JNN. The calculations are applied to estimate the impact of JNN on measurements with the new apparatus, n2EDM, at the Paul Scherrer Institute. We demonstrate that the performances of the optically pumped $^{133}$Cs magnetometers and $^{199}$Hg co-magnetometers, which will be used in the apparatus, are not limited by JNN. Further, we find that in measurements deploying a co-magnetometer system, the impact of JNN is negligible for nEDM searches down to a sensitivity of $4\,\times\,10^{-28}\,e\cdot{\rm cm}$ in a single measurement; therefore, the use of economically and mechanically favored solid aluminum electrodes is possible.


title: Capacity and quantum geometry of parametrized quantum circuits

abstract: To harness the potential of noisy intermediate-scale quantum devices, it is paramount to find the best type of circuits to run hybrid quantum-classical algorithms. Key candidates are parametrized quantum circuits that can be effectively implemented on current devices. Here, we evaluate the capacity and trainability of these circuits using the geometric structure of the parameter space via the effective quantum dimension, which reveals the expressive power of circuits in general as well as of particular initialization strategies. We assess the expressive power of various popular circuit types and find striking differences depending on the type of entangling gates used. Particular circuits are characterized by scaling laws in their expressiveness. We identify a transition in the quantum geometry of the parameter space, which leads to a decay of the quantum natural gradient for deep circuits. For shallow circuits, the quantum natural gradient can be orders of magnitude larger in value compared to the regular gradient; however, both of them can suffer from vanishing gradients. By tuning a fixed set of circuit parameters to randomized ones, we find a region where the circuit is expressive, but does not suffer from barren plateaus, hinting at a good way to initialize circuits. We show an algorithm that prunes redundant parameters of a circuit without affecting its effective dimension. Our results enhance the understanding of parametrized quantum circuits and can be immediately applied to improve variational quantum algorithms.


title: Maxwell Fields in Boosted Kerr Black Holes

abstract: The spacetime of a boosted Bondi-Sachs rotating black hole is considered as a proper background to examine electromagnetic configurations connected to analytic solutions of Maxwell equations. In our analysis, we first use the Bondi-Sachs transformations in order to bring the boosted rotating black hole metric into the Kerr-Schild form, from which zero angular momentum observers (ZAMOs) are constructed via the ADM formalism. In Kerr-Schild coordinates we obtain the Killing fields as sources of Maxwell electrodynamics, and we fix a ZAMO in order to evaluate the components of the electric and magnetic fields, from which we obtain nonsingular patterns of an eventual momentum-energy emission of a boosted Kerr-Schild black hole. Distinct patterns are examined and discussed in the case of variations of the boost parameter $\gamma$. We extend our analysis by considering the nonsingular electromagnetic emission in the framework of a boosted Bondi-Sachs rotating black hole, as it moves at relativistic speeds. We also discuss possible mechanisms that may resemble magnetospheres of rotating boosted black holes and give rise to hydromagnetic flows from accretion discs and to the production of jets.


title: The Milky Way Project: Probing Star Formation with First Results on Yellowballs from DR2

abstract: Yellowballs (YBs) were first discovered during the Milky Way Project citizen-science initiative (MWP; Simpson et al. 2012). MWP users noticed compact, yellow regions in Spitzer Space Telescope mid-infrared (MIR) images of the Milky Way plane and asked professional astronomers to explain these "yellow balls." Follow-up work by Kerton et al. (2015) determined that YBs likely trace compact photo-dissociation regions associated with massive and intermediate-mass star formation. YBs were included as target objects in a version of the Milky Way Project launched in 2016 (Jayasinghe et al. 2016), which produced a listing of over 6000 YB locations. We have measured distances, cross-match associations, physical properties, and MIR colors of ~500 YBs within a pilot region covering the l= 30 - 40 degrees, b= +/- 1 degree region of the Galactic plane. We find 20-30% of YBs in our pilot region contain high-mass star formation capable of becoming expanding H II regions that produce MIR bubbles. A majority of YBs represent intermediate-mass star-forming regions whose placement in evolutionary diagrams suggest they are still actively accreting, and may be precursors to optically-revealed Herbig Ae/Be nebulae. Many of these intermediate-mass YBs were missed by surveys of massive star-formation tracers and thus this catalog provides information for many new sites of star formation. Future work will expand this pilot region analysis to the entire YB catalog.


title: Private Linear Transformation: The Individual Privacy Case

abstract: This paper considers the single-server Private Linear Transformation (PLT) problem when individual privacy is required. In this problem, there is a user that wishes to obtain $L$ linear combinations of a $D$-subset of messages belonging to a dataset of $K$ messages stored on a single server. The goal is to minimize the download cost while keeping the identity of every message required for the computation individually private. The individual privacy requirement implies that, from the perspective of the server, every message is equally likely to belong to the $D$-subset of messages that constitute the support set of the required linear combinations. We focus on the setting in which the matrix of coefficients pertaining to the required linear combinations is the generator matrix of a Maximum Distance Separable code. We establish lower and upper bounds on the capacity of PLT with individual privacy, where the capacity is defined as the supremum of all achievable download rates. We show that our bounds are tight under certain divisibility conditions. In addition, we present lower bounds on the capacity of the settings in which the user has a prior side information about a subset of messages.


title: Interpolated family of non group-like simple integral fusion rings of Lie type

abstract: This paper is motivated by the quest of a non-group irreducible finite index depth 2 maximal subfactor. We compute the generic fusion rules of the Grothendieck ring of Rep(PSL(2,q)), q prime-power, by applying a Verlinde-like formula on the generic character table. We then prove that this family of fusion rings R_q interpolates to all integers q>=2, providing (when q is not prime-power) the first example of infinite family of non group-like simple integral fusion rings. Furthermore, they pass all the known criteria of (unitary) categorification. This provides infinitely many serious candidates for solving the famous open problem of whether there exists an integral fusion category which is not weakly group-theoretical. We prove that a complex categorification (if any) of an interpolated fusion ring R_q (with q non prime-power) cannot be braided, and so its Drinfeld center must be simple. In general, this paper proves that a non-pointed simple fusion category is non-braided if and only if its Drinfeld center is simple; and also that every simple integral fusion category is weakly group-theoretical if and only if every simple integral modular fusion category is pointed.


title: W*-rigidity paradigms for embeddings of II$_1$ factors

abstract: We undertake a systematic study of W*-rigidity paradigms for the embeddability relation $\hookrightarrow$ between separable II$_1$ factors and its stable version $\hookrightarrow_s$, obtaining large families of non stably isomorphic II$_1$ factors that are mutually embeddable and families of II$_1$ factors that are mutually non stably embeddable. We provide an augmentation functor $G \mapsto H_G$ from the category of groups into icc groups, so that $L(H_{G_1}) \hookrightarrow_s L(H_{G_2})$ iff $G_1 \hookrightarrow G_2$. We construct complete intervals of II$_1$ factors, including a strict chain of II$_1$ factors $(M_k)_{k\in \mathbb{Z}}$ with the property that if $N$ is any II$_1$ factor with $M_i \hookrightarrow_s N$ and $N \hookrightarrow_s M_{j}$, then $N \cong M_k^t$ for some $i \leq k \leq j$ and $t > 0$.


title: Private Linear Transformation: The Joint Privacy Case

abstract: We introduce the problem of Private Linear Transformation (PLT). This problem includes a single (or multiple) remote server(s) storing (identical copies of) $K$ messages and a user who wants to compute $L$ linear combinations of a $D$-subset of these messages by downloading the minimum amount of information from the server(s) while protecting the privacy of the entire set of $D$ messages. This problem generalizes the Private Information Retrieval and Private Linear Computation problems. In this work, we focus on the single-server case. For the setting in which the coefficient matrix of the required $L$ linear combinations generates a Maximum Distance Separable (MDS) code, we characterize the capacity -- defined as the supremum of all achievable download rates, for all parameters $K, D, L$. In addition, we present lower and/or upper bounds on the capacity for the settings with non-MDS coefficient matrices and the settings with a prior side information.


title: Ubiquitous Dynamical Time Asymmetry in Measurements on Materials and Biological Systems

abstract: Many measurements on soft condensed matter (e.g., biological and materials) systems track low-dimensional observables projected from the full system phase space as a function of time. Examples are dynamic structure factors, spectroscopic and rheological response functions, and time series of distances derived from optical tweezers, single-molecule spectroscopy and molecular dynamics simulations. In many such systems the projection renders the reduced dynamics non-Markovian and the observable is not prepared in, or initially sampled from and averaged over, a stationary distribution. We prove that such systems always exhibit non-equilibrium, time asymmetric dynamics. That is, they evolve in time with a broken time-translation invariance in a manner closely resembling aging dynamics. We identify the entropy associated with the breaking of time-translation symmetry that is a measure of the instantaneous thermodynamic displacement of latent, hidden degrees of freedom from their stationary state. Dynamical time asymmetry is a general phenomenon, independent of the underlying energy surface, and is frequently even visible in measurements on systems that have fully reached equilibrium. This finding has fundamental implications for the interpretation of many experiments on, and simulations of, biological and materials systems.

categories: cond-mat.stat-mech cond-mat.soft math-ph math.MP math.PR physics.bio-ph


title: Pulling the Higgs and Top needles from the jet stack with Feature Extended Supervised Tagging

abstract: Jet tagging has become an essential tool for new physics searches at the high-energy frontier. For jets that contain energetic charged leptons we introduce Feature Extended Supervised Tagging (FEST) which, in addition to jet substructure, considers the features of the charged lepton within the jet. With this method we build dedicated taggers to discriminate among boosted $H \to \ell \nu q \bar q$, $t \to \ell \nu b$, and QCD jets (with $\ell$ an electron or muon). The taggers have an impressive performance, allowing for overall light jet rejection factors of $10^4-10^5$, for top quark / Higgs boson efficiencies of $0.5$. The taggers are also excellent in the discrimination of Higgs bosons from top quarks and vice versa, for example rejecting top quarks by factors of $100-300$ for Higgs boson efficiencies of $0.5$. We demonstrate the potential of these taggers to improve the sensitivity to new physics by using as example a search for a new $Z'$ boson decaying into $Z H$, in the fully-hadronic final state.


title: Noncommutative-geometry wormholes with isotropic pressure

abstract: The strategy adopted in the original Morris-Thorne wormhole was to retain complete control over the geometry at the expense of certain engineering considerations. The purpose of this paper is to obtain several complete wormhole solutions by assuming a noncommutative-geometry background with a concomitant isotropic-pressure condition. This condition allows us to consider a cosmological setting with a perfect-fluid equation of state. An extended form of the equation generalizes the first solution and subsequently leads to the generalized Chaplygin-gas model and hence to a third solution. The solutions obtained extend several previous results. This paper also reiterates the need for a noncommutative-geometry background to account for the enormous radial tension that is characteristic of Morris-Thorne wormholes.


title: Sharp estimates for spanning trees

abstract: We prove the following sharp estimate for the number of spanning trees of a graph in terms of its vertex-degrees: a simple graph $G$ on $n$ vertices has at most $(1/n^{2}) \prod_{v \in V(G)} (d(v)+1)$ spanning trees. This result is tight (for complete graphs), and improves earlier estimates of Alon from 1990 and Kostochka from 1995 by a factor of about $1/n$ (for dense graphs). We additionally show that an analogous bound holds for the weighted spanning tree enumerator of a (nonnegatively) weighted graph as well.


title: Keep the Gradients Flowing: Using Gradient Flow to Study Sparse Network Optimization

abstract: Training sparse networks to converge to the same performance as dense neural architectures has proven to be elusive. Recent work suggests that initialization is the key. However, while this direction of research has had some success, focusing on initialization alone appears to be inadequate. In this paper, we take a broader view of training sparse networks and consider the role of regularization, optimization, and architecture choices on sparse models. We propose a simple experimental framework, Same Capacity Sparse vs Dense Comparison (SC-SDC), that allows for a fair comparison of sparse and dense networks. Furthermore, we propose a new measure of gradient flow, Effective Gradient Flow (EGF), that better correlates to performance in sparse networks. Using top-line metrics, SC-SDC and EGF, we show that default choices of optimizers, activation functions and regularizers used for dense networks can disadvantage sparse networks. Based upon these findings, we show that gradient flow in sparse networks can be improved by reconsidering aspects of the architecture design and the training regime. Our work suggests that initialization is only one piece of the puzzle and taking a wider view of tailoring optimization to sparse networks yields promising results.


title: Reed-Muller Subcodes: Machine Learning-Aided Design of Efficient Soft Recursive Decoding

abstract: Reed-Muller (RM) codes are conjectured to achieve the capacity of any binary-input memoryless symmetric (BMS) channel, and are observed to have a comparable performance to that of random codes in terms of scaling laws. On the negative side, RM codes lack efficient decoders with performance close to that of a maximum likelihood decoder for general parameters. Also, they only admit certain discrete sets of rates. In this paper, we focus on subcodes of RM codes with flexible rates that can take any code dimension from 1 to n, where n is the blocklength. We first extend the recursive projection-aggregation (RPA) algorithm proposed recently by Ye and Abbe for decoding RM codes. To lower the complexity of our decoding algorithm, referred to as subRPA in this paper, we investigate different ways for pruning the projections. We then derive the soft-decision based version of our algorithm, called soft-subRPA, that is shown to improve upon the performance of subRPA. Furthermore, it enables training a machine learning (ML) model to search for \textit{good} sets of projections in the sense of minimizing the decoding error rate. Training our ML model enables achieving very close to the performance of full-projection decoding with a significantly reduced number of projections. For instance, our simulation results on a (64,14) RM subcode show almost identical performance for full-projection decoding and pruned-projection decoding with 15 projections picked via training our ML model. This is equivalent to lowering the complexity by a factor of more than 4 without sacrificing the decoding performance.


title: The GEM Benchmark: Natural Language Generation, its Evaluation and Metrics

abstract: We introduce GEM, a living benchmark for natural language Generation (NLG), its Evaluation, and Metrics. Measuring progress in NLG relies on a constantly evolving ecosystem of automated metrics, datasets, and human evaluation standards. Due to this moving target, new models often still evaluate on divergent anglo-centric corpora with well-established, but flawed, metrics. This disconnect makes it challenging to identify the limitations of current models and opportunities for progress. Addressing this limitation, GEM provides an environment in which models can easily be applied to a wide set of tasks and in which evaluation strategies can be tested. Regular updates to the benchmark will help NLG research become more multilingual and evolve the challenge alongside models. This paper serves as the description of the data for which we are organizing a shared task at our ACL 2021 Workshop and to which we invite the entire NLG community to participate.


title: Universal systole bounds for arithmetic locally symmetric spaces

abstract: The systole of a closed Riemannian manifold is the minimal length of a non-contractible closed loop. We give a uniform lower bound for the systole for large classes of simple arithmetic locally symmetric orbifolds. We establish new bounds for the translation length of a semisimple element x in SL_n(R) in terms of its associated Mahler measure. We use these geometric methods to prove the existence of extensions of number fields in which fixed sets of primes have certain prescribed splitting behavior.


title: The simplest erasing substitution

abstract: In this work, we begin the study of a new class of dynamical systems determined by interval maps generated by the symbolic action of erasing substitution rules. We do this by discussing in some detail the geometric, analytical, dynamical and arithmetic properties of a particular example, which has the virtue of being arguably the simplest and that at the same time produces interesting properties and new challenging problems.


title: A proposal for realising quantum scars in the tilted 1D Fermi-Hubbard model

abstract: Motivated by recent observations of ergodicity breaking due to Hilbert space fragmentation in 1D Fermi-Hubbard chains with a tilted potential [Scherg et al., arXiv:2010.12965], we show that the same system also hosts quantum many-body scars in a regime $U\approx \Delta \gg J$ at electronic filling factor $\nu=1$. We numerically demonstrate that the scarring phenomenology in this model is similar to other known realisations such as Rydberg atom chains, including persistent dynamical revivals and ergodicity-breaking many-body eigenstates. At the same time, we show that the mechanism of scarring in the Fermi-Hubbard model is different from other examples in the literature: the scars originate from a subgraph, representing a free spin-1 paramagnet, which is weakly connected to the rest of the Hamiltonian's adjacency graph. Our work demonstrates that correlated fermions in tilted optical lattices provide a platform for understanding the interplay of many-body scarring and other forms of ergodicity breaking, such as localisation and Hilbert space fragmentation.


title: Anomalous gravitomagnetic moment and non-universality of the axial vortical effect at finite temperature

abstract: The coupling between the spin of a massive Dirac fermion and the angular momentum of the medium, i.e. the gravitomagnetic moment, is shown here to be renormalized by QED interactions at finite temperature. This means that the anomalous gravitomagnetic moment (AGM) does not vanish, and implies that thermal effects can break the Einstein equivalence principle in quantum field theory, as argued previously. We also show that the AGM causes radiative corrections to the axial current of massive fermions induced by vorticity in quantum relativistic fluids, similarly to the previous findings for massless fermions. The radiative QCD effects on the AGM should significantly affect the production of polarized hadrons in heavy-ion collisions.


title: Trace of $\Lambda(1405)$ resonance in low energy $K^{-}+\, ^{3}\mathrm{He}\rightarrow (\pi^{0}\Sigma^{0})+d$ reaction

abstract: In present work, we investigated $K^{-}+\, ^{3}\mathrm{He}$ reaction at low energies. The coupled-channel Faddeev AGS equations were solved for $\bar{K}Nd-\pi\Sigma{d}$ three-body system in momentum representation to extract the scattering amplitudes. To trace the signature of the $\Lambda$(1405) resonance in the $\pi\Sigma$ invariant mass, the deuteron energy spectrum for $K^{-}+\,^{3}\mathrm{He}\rightarrow\pi\Sigma{d}$ reaction was obtained. Different types of $\bar{K}N-\pi\Sigma$ potentials based on phenomenological and chiral SU(3) approaches were used. As a remarkable result of this investigation, it was found that the deuteron energy spectrum, reflecting the $\Lambda$(1405) mass distribution and width, depends quite sensitively on the $\bar{K}N-\pi\Sigma$ model of interaction. Hence accurate measurements of the $\pi\Sigma$ mass distribution have the potential to discriminate between possible mechanisms at work in the formation of the $\Lambda$(1405).


title: Learning domain-agnostic visual representation for computational pathology using medically-irrelevant style transfer augmentation

abstract: Suboptimal generalization of machine learning models on unseen data is a key challenge which hampers the clinical applicability of such models to medical imaging. Although various methods such as domain adaptation and domain generalization have evolved to combat this challenge, learning robust and generalizable representations is core to medical image understanding, and continues to be a problem. Here, we propose STRAP (Style TRansfer Augmentation for histoPathology), a form of data augmentation based on random style transfer from non-medical style source such as artistic paintings, for learning domain-agnostic visual representations in computational pathology. Style transfer replaces the low-level texture content of an image with the uninformative style of randomly selected style source image, while preserving the original high-level semantic content. This improves robustness to domain shift and can be used as a simple yet powerful tool for learning domain-agnostic representations. We demonstrate that STRAP leads to state-of-the-art performance, particularly in the presence of domain shifts, on two particular classification tasks in computational pathology.


title: Anomalous light-induced broadening of the spin-noise resonance in cesium vapor

abstract: We uncover a highly nontrivial dependence of the spin-noise (SN) resonance broadening induced by the intense probe beam. The measurements were performed by probing the cell with cesium vapor at the wavelengths of the transition ${6}^2S_{1/2} \leftrightarrow {6}^2P_{3/2}$ ($\mathrm{D}_2$ line) with the unresolved hyperfine structure of the excited state. The light-induced broadening of the SN resonance was found to differ strongly at different slopes of the $\mathrm{D}_2$ line and, generally, varied nonmonotonically with light power. We discuss the effect in terms of the phenomenological Bloch equations for the spin fluctuations and demonstrate that the SN broadening behavior strongly depends on the relation between the pumping and excited-level decay rates, the spin precession, and decoherence rates. To reconcile the puzzling experimental results, we propose that the degree of optical perturbation of the spin-system is controlled by the route of the excited-state relaxation of the atom or, in other words, that the act of optical excitation of the atom does not necessarily break down completely its ground-state coherence and continuity of the spin precession. Spectral asymmetry of the effect, in this case, is provided by the position of the "closed" transition $F = 4 \leftrightarrow F' = 5$ at the short-wavelength side of the line. This hypothesis, however, remains to be proven by microscopic calculations.


title: A double copy for asymptotic symmetries in the self-dual sector

abstract: We give a double copy construction for the symmetries of the self-dual sectors of Yang-Mills (YM) and gravity, in the light-cone formulation. We find an infinite set of double copy constructible symmetries. We focus on two families which correspond to the residual diffeomorphisms on the gravitational side. For the first one, we find novel non-perturbative double copy rules in the bulk. The second family has a more striking structure, as a non-perturbative gravitational symmetry is obtained from a perturbatively defined symmetry on the YM side. At null infinity, we find the YM origin of the subset of extended Bondi-Metzner-Sachs (BMS) symmetries that preserve the self-duality condition. In particular, holomorphic large gauge YM symmetries are double copied to holomorphic supertranslations. We also identify the single copy of superrotations with certain non-gauge YM transformations that to our knowledge have not been previously presented in the literature.


title: Assessing Time-Varying Causal Effect Moderation in the Presence of Cluster-Level Treatment Effect Heterogeneity

abstract: The micro-randomized trial (MRT) is a sequential randomized experimental design to empirically evaluate the effectiveness of mobile health (mHealth) intervention components that may be delivered at hundreds or thousands of decision points. MRTs have motivated a new class of causal estimands, termed "causal excursion effects", for which semiparametric inference can be conducted via a weighted, centered least squares criterion (Boruvka et al., 2018). Existing methods assume between-subject independence and non-interference. Deviations from these assumptions often occur. In this paper, causal excursion effects are revisited under potential cluster-level treatment effect heterogeneity and interference, where the treatment effect of interest may depend on cluster-level moderators. Utility of the proposed methods is shown by analyzing data from a multi-institution cohort of first year medical residents in the United States.


title: Exploiting dynamic quantum circuits in a quantum algorithm with superconducting qubits

abstract: The execution of quantum circuits on real systems has largely been limited to those which are simply time-ordered sequences of unitary operations followed by a projective measurement. As hardware platforms for quantum computing continue to mature in size and capability, it is imperative to enable quantum circuits beyond their conventional construction. Here we break into the realm of dynamic quantum circuits on a superconducting-based quantum system. Dynamic quantum circuits involve not only the evolution of the quantum state throughout the computation, but also periodic measurements of a subset of qubits mid-circuit and concurrent processing of the resulting classical information within timescales shorter than the execution times of the circuits. Using noisy quantum hardware, we explore one of the most fundamental quantum algorithms, quantum phase estimation, in its adaptive version, which exploits dynamic circuits, and compare the results to a non-adaptive implementation of the same algorithm. We demonstrate that the version of real-time quantum computing with dynamic circuits can offer a substantial and tangible advantage when noise and latency are sufficiently low in the system, opening the door to a new realm of available algorithms on real quantum systems.


title: Forces within hadrons on the light front

abstract: In this work, we find the light front densities for momentum and forces, including pressure and shear forces, within hadrons. This is achieved by deriving relativistically correct expressions relating these densities to the gravitational form factors $A(t)$ and $D(t)$ associated with the energy momentum tensor. The derivation begins from the fundamental definition of density in a quantum field theory, namely the expectation value of a local operator within a spatially-localized state. We find that it is necessary to use the light front formalism to define a density that corresponds to internal hadron structure. When using the instant form formalism, it is impossible to remove the spatial extent of the hadron wave function from any density, and -- even within instant form dynamics -- one does not obtain a Breit frame Fourier transform for a properly defined density. Within the front formalism, we derive new expressions for various mechanical properties of hadrons, including the mechanical radius, as well as for stability conditions. The multipole ansatz for the form factors is used as an example to illustrate all of these findings.


title: Popular differences for matrix patterns

abstract: The following combinatorial conjecture arises naturally from recent ergodic-theoretic work of Ackelsberg, Bergelson, and Best. Let $M_1$, $M_2$ be $k\times k$ integer matrices, $G$ be a finite abelian group of order $N$, and $A\subseteq G^k$ with $|A|\ge\alpha N^k$. If $M_1$, $M_2$, $M_1-M_2$, and $M_1+M_2$ are automorphisms of $G^k$, is it true that there exists a popular difference $d \in G^k\setminus\{0\}$ such that \[\#\{x \in G^k: x, x+M_1d, x+M_2d, x+(M_1+M_2)d \in A\} \ge (\alpha^4-o(1))N^k.\] We show that this conjecture is false in general, but holds for $G = \mathbb{F}_p^n$ with $p$ an odd prime given the additional spectral condition that no pair of eigenvalues of $M_1M_2^{-1}$ (over $\overline{\mathbb{F}}_p$) are negatives of each other. In particular, the "rotated squares" pattern does not satisfy this eigenvalue condition, and we give a construction of a set of positive density in $(\mathbb{F}_5^n)^2$ for which that pattern has no nonzero popular difference. This is in surprising contrast to three-point patterns, which we handle over all compact abelian groups and which do not require an additional spectral condition.


title: Agent Incentives: A Causal Perspective

abstract: We present a framework for analysing agent incentives using causal influence diagrams. We establish that a well-known criterion for value of information is complete. We propose a new graphical criterion for value of control, establishing its soundness and completeness. We also introduce two new concepts for incentive analysis: response incentives indicate which changes in the environment affect an optimal decision, while instrumental control incentives establish whether an agent can influence its utility via a variable X. For both new concepts, we provide sound and complete graphical criteria. We show by example how these results can help with evaluating the safety and fairness of an AI system.


title: Electrically Switchable Broadband Photonic Bound States in the Continuum

abstract: The question of how to continuously manipulate a photonic system between a radiating state and a bound state is an important challenge in photonics, as its solution promises broad technological relevance for optical sensors, modulators, switches and displays. Existing approaches utilise the inherent wave-nature of electromagnetic fields to their advantage, and are commonly identified as bound states in the continuum (BICs), as they resemble singular bound states embedded in a band of radiating states. Although quasi-BICs have been demonstrated for numerous symmetric periodic photonic crystals, their existence so far has been limited to narrow spectral ranges, and their application in switches or modulators needs large changes of the refractive index. Here, we show that the incidence of two guided symmetric substrate waves of opposite phase onto a symmetric periodic film waveguide enables the excitation of self-stabilising BICs which are electrically switchable in a broad optical range, with small changes of the refractive index. An experimental verification of the concept shows a switching contrast C=700 at a wavelength of 532 nm and C=1000 at a wavelength of 632.8 nm.


title: Report of the Workshop on Program Synthesis for Scientific Computing

abstract: Program synthesis is an active research field in academia, national labs, and industry. Yet, work directly applicable to scientific computing, while having some impressive successes, has been limited. This report reviews the relevant areas of program synthesis work for scientific computing, discusses successes to date, and outlines opportunities for future work. This report is the result of the Workshop on Program Synthesis for Scientific Computing was held virtually on August 4-5 2020 (https://prog-synth-science.github.io/2020/).


title: The regularized free fall I -- Index computations

abstract: Main results are, firstly, a generalization of the Conley-Zehnder index from ODEs to the delay equation at hand and, secondly, the equality of the Morse index and the clockwise normalized Conley-Zehnder index.


title: Spin misalignment of black hole binaries from young star clusters: implications for the origin of gravitational waves events

abstract: Recent studies indicate that the progenitors of merging black hole (BH) binaries from young star clusters can undergo a common envelope phase just like isolated binaries. If the stars emerge from the common envelope as naked cores, tidal interactions can efficiently synchronize their spins before they collapse into BHs. Contrary to the isolated case, these binary BHs can also undergo dynamical interactions with other BHs in the cluster before merging. The interactions can tilt the binary orbital plane, leading to spin-orbit misalignment. We estimate the spin properties of merging binary BHs undergoing this scenario by combining up-to-date binary population synthesis and accurate few-body simulations. We show that post-common envelope binary BHs are likely to undergo only a single encounter, due to the high binary recoil velocity and short coalescence times. Adopting conservative limits on the binary-single encounter rates, we obtain a local BH merger rate density of ~6.6 yr^-1 Gpc^-3. Assuming low (<0.2) natal BH spins, this scenario reproduces the trends in the distributions of effective spin Xeff and precession parameters Xp inferred from GWTC-2, including the peaks at (Xeff, Xp) ~ (0.1, 0.2) and the tail at negative Xeff values.


title: From Culture to Clothing: Discovering the World Events Behind A Century of Fashion Images

abstract: Fashion is intertwined with external cultural factors, but identifying these links remains a manual process limited to only the most salient phenomena. We propose a data-driven approach to identify specific cultural factors affecting the clothes people wear. Using large-scale datasets of news articles and vintage photos spanning a century, we present a multi-modal statistical model to detect influence relationships between happenings in the world and people's choice of clothing. Furthermore, on two image datasets we apply our model to improve the concrete vision tasks of visual style forecasting and photo timestamping. Our work is a first step towards a computational, scalable, and easily refreshable approach to link culture to clothing.


title: Exact Langevin Dynamics with Stochastic Gradients

abstract: Stochastic gradient Markov Chain Monte Carlo algorithms are popular samplers for approximate inference, but they are generally biased. We show that many recent versions of these methods (e.g. Chen et al. (2014)) cannot be corrected using Metropolis-Hastings rejection sampling, because their acceptance probability is always zero. We can fix this by employing a sampler with realizable backwards trajectories, such as Gradient-Guided Monte Carlo (Horowitz, 1991), which generalizes stochastic gradient Langevin dynamics (Welling and Teh, 2011) and Hamiltonian Monte Carlo. We show that this sampler can be used with stochastic gradients, yielding nonzero acceptance probabilities, which can be computed even across multiple steps.


title: Generacion de voces artificiales infantiles en castellano con acento costarricense

abstract: This article evaluates a first experience of generating artificial children's voices with a Costa Rican accent, using the technique of statistical parametric speech synthesis based on Hidden Markov Models. The process of recording the voice samples used for learning the models, the fundamentals of the technique used and the subjective evaluation of the results through the perception of a group of people is described. The results show that the intelligibility of the results, evaluated in isolated words, is lower than the voices recorded by the group of participating children. Similarly, the detection of the age and gender of the speaking person is significantly affected in artificial voices, relative to recordings of natural voices. These results show the need to obtain larger amounts of data, in addition to becoming a numerical reference for future developments resulting from new data or from processes to improve results in the same technique.


title: 1-form Symmetries of 4d N=2 Class S Theories

abstract: We determine the 1-form symmetry group for any 4d N = 2 class S theory constructed by compactifying a 6d N=(2,0) SCFT on a Riemann surface with arbitrary regular untwisted and twisted punctures. The 6d theory has a group of mutually non-local dimension-2 surface operators, modulo screening. Compactifying these surface operators leads to a group of mutually non-local line operators in 4d, modulo screening and flavor charges. Complete specification of a 4d theory arising from such a compactification requires a choice of a maximal subgroup of mutually local line operators, and the 1-form symmetry group of the chosen 4d theory is identified as the Pontryagin dual of this maximal subgroup. We also comment on how to generalize our results to compactifications involving irregular punctures. Finally, to complement the analysis from 6d, we derive the 1-form symmetry from a Type IIB realization of class S theories.


title: Magnetically Modified Spherical Accretion in GRMHD: Reconnection-Driven Convection and Jet Propagation

abstract: We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulence in a magnetically arrested (or a nearly arrested) state. Electromagnetic jets with outflow efficiencies ~10-200% occupy the polar regions, reaching several hundred gravitational radii before they dissipate due to the kink instability. The jet directions fluctuate in time and can be tilted by as much as ~30 degrees with respect to black hole spin, but this tilt does not depend strongly on the tilt of the initial magnetic field. A jet forms even when there is no initial net vertical magnetic flux since turbulent, horizon-scale fluctuations can generate a net vertical field locally. Peak jet power is obtained for an initial magnetic field tilted by 40-80 degrees with respect to the black hole spin because this maximizes the amount of magnetic flux that can reach the black hole. These simulations may be a reasonable model for low luminosity black hole accretion flows such as Sgr A* or M87.


title: Observation of a prethermal discrete time crystal

abstract: The conventional framework for defining and understanding phases of matter requires thermodynamic equilibrium. Extensions to non-equilibrium systems have led to surprising insights into the nature of many-body thermalization and the discovery of novel phases of matter, often catalyzed by driving the system periodically. The inherent heating from such Floquet drives can be tempered by including strong disorder in the system, but this can also mask the generality of non-equilibrium phases. In this work, we utilize a trapped-ion quantum simulator to observe signatures of a non-equilibrium driven phase without disorder: the prethermal discrete time crystal (PDTC). Here, many-body heating is suppressed not by disorder-induced many-body localization, but instead via high-frequency driving, leading to an expansive time window where non-equilibrium phases can emerge. We observe a number of key features that distinguish the PDTC from its many-body-localized disordered counterpart, such as the drive-frequency control of its lifetime and the dependence of time-crystalline order on the energy density of the initial state. Floquet prethermalization is thus presented as a general strategy for creating, stabilizing and studying intrinsically out-of-equilibrium phases of matter.


title: Information Content of Higher-Order Galaxy Correlation Functions

abstract: The shapes of galaxy N-point correlation functions can be used as standard rulers to constrain the distance-redshift relationship and thence the expansion rate of the Universe. The cosmological density fields traced by late-time galaxy formation are initially nearly Gaussian, and hence all the cosmological information can be extracted from their 2-Point Correlation Function (2PCF) or its Fourier-space analog the power spectrum. Subsequent nonlinear evolution under gravity, as well as halo and then galaxy formation, generate higher-order correlation functions. Since the mapping of the initial to the final density field is, on large scales, invertible, it is often claimed that the information content of the initial field's power spectrum is equal to that of all the higher-order functions of the final, nonlinear field. This claim implies that reconstruction of the initial density field from the nonlinear field renders analysis of higher-order correlation functions of the latter superfluous. We here show that this claim is false when the N-point functions are used as standard rulers. Constraints available from joint analysis of the galaxy power spectrum and bispectrum (Fourier-space analog of the 3-Point Correlation Function) can, in some cases, exceed those offered by the initial power spectrum even when the reconstruction is perfect. We provide a mathematical justification for this claim and also demonstrate it using a large suite of N-body simulations. In particular, we show that for the z = 0 real-space matter field in the limit of vanishing shot noise, taking modes up to k_max = 0.2 h/Mpc, using the bispectrum alone offers a factor of two reduction in the variance on the cosmic distance scale relative to that available from the power spectrum.


title: Mapping stellar surfaces II: An interpretable Gaussian process model for light curves

abstract: The use of Gaussian processes (GPs) as models for astronomical time series datasets has recently become almost ubiquitous, given their ease of use and flexibility. GPs excel in particular at marginalization over the stellar signal in cases where the variability due to starspots rotating in and out of view is treated as a nuisance, such as in exoplanet transit modeling. However, these effective models are less useful in cases where the starspot signal is of primary interest since it is not obvious how the parameters of the GP model are related to the physical properties of interest, such as the size, contrast, and latitudinal distribution of the spots. Instead, it is common practice to explicitly model the effect of individual starspots on the light curve and attempt to infer their properties via optimization or posterior inference. Unfortunately, this process is degenerate, ill-posed, and often computationally intractable when applied to stars with more than a few spots and/or to ensembles of many light curves. In this paper, we derive a closed-form expression for the mean and covariance of a Gaussian process model that describes the light curve of a rotating, evolving stellar surface conditioned on a given distribution of starspot sizes, contrasts, and latitudes. We demonstrate that this model is correctly calibrated, allowing one to robustly infer physical parameters of interest from one or more stellar light curves, including the typical radii and the mean and variance of the latitude distribution of starspots. Our GP has far-ranging implications for understanding the variability and magnetic activity of stars from both light curves and radial velocity (RV) measurements, as well as for robustly modeling correlated noise in both transiting and RV exoplanet searches. Our implementation is efficient, user-friendly, and open source, available as the Python package starry-process.


title: The Cosmological Trajectories Method: Modelling cosmic structure formation in the non-linear regime

abstract: We introduce a novel approach, the Cosmological Trajectories Method (CTM), to model nonlinear structure formation in the Universe by expanding gravitationally-induced particle trajectories around the Zel'dovich approximation. A new Beyond Zel'dovich approximation is presented, which expands the CTM to leading second-order in the gravitational interaction and allows for post-Born gravitational scattering. In the Beyond Zel'dovich approximation we derive the exact expression for the matter clustering power spectrum. This is calculated to leading order and is available in the CTM MODULE. We compare the Beyond Zel'dovich approximation power spectrum and correlation function to other methods including 1-loop Standard Perturbation Theory (SPT), 1-loop Lagrangian Perturbation Theory (LPT) and Convolution Lagrangian Perturbation Theory (CLPT). We find that the Beyond Zel'dovich approximation power spectrum performs well, matching simulations to within $\pm{10}\%$, on mildly non-linear scales, and at redshifts above $z=1$ it outperforms the Zel'dovich approximation. We also find that the Beyond Zel'dovich approximation models the BAO peak in the correlation function at $z=0$ more accurately, to within $\pm{5}\%$ of simulations, than the Zel'dovich approximation, SPT 1-loop and CLPT.


title: Divergence of the Gr\"uneisen ratio at symmetry-enhanced first-order quantum phase transitions

abstract: Studies of the Gr\"uneisen ratio, i.e., the ratio between thermal expansion and specific heat, have become a powerful tool in the context of quantum criticality, since it was shown theoretically that the Gr\"uneisen ratio displays characteristic power-law divergencies upon approaching the transition point of a continuous quantum phase transition. Here we show that the Gr\"uneisen ratio also diverges at a symmetry-enhanced first-order quantum phase transition, albeit with mean-field exponents, as the enhanced symmetry implies the vanishing of a mode gap which is finite away from the transition. We provide explicit results for simple pseudo-spin models, both with and without Goldstone modes in the stable phases, and discuss implications.


title: Critical anomalous metals near superconductivity in models with random interactions

abstract: Anomalous metals are observed in numerous experiments on disordered two-dimensional systems proximate to superconductivity. A characteristic feature of an anomalous metal is that its low temperature conductivity has a weakly temperature dependent value, significantly higher than that of a disordered Fermi liquid. We propose a dynamical mean-field model of an anomalous metal: interacting electrons similar in structure to that of the well-studied universal Hamiltonian of mesoscopic metallic grains, but with independent random interactions between pairs of sites, involving Cooper pair hopping and spin exchange. We find evidence for critical anomalous phases or points between a superconducting phase and a disordered Fermi liquid phase in this model. Our results are obtained by a renormalization group analysis in a weak coupling limit, and a complementary solution at large $M$ when the spin symmetry is generalized to USp($M$). The large $M$ limit describes the anomalous metal by fractionalization of the electron into spinons, holons, and doublons, with these partons forming critical non-Fermi liquids in the Sachdev-Ye-Kitaev class. We compute the low temperature conductivity in the large $M$ limit, and find temperature-independent values moderately enhanced from that in the disordered metal.


title: Living with Neighbors. III. The Origin of the Spin$-$Orbit Alignment of Galaxy Pairs: A Neighbor versus the Large-scale Structure

abstract: Recent observations revealed a coherence between the spin vector of a galaxy and the orbital motion of its neighbors. We refer to the phenomenon as "the spin$-$orbit alignment (SOA)" and explore its physical origin via the IllustrisTNG simulation. This is the first study to utilize a cosmological hydrodynamic simulation to investigate the SOA of galaxy pairs. In particular, we identify paired galaxies at $z = 0$ having the nearest neighbor with mass ratios from 1/10 to 10 and calculate the spin$-$orbit angle for each pair. Our results are as follows. (a) There exists a clear preference for prograde orientations (i.e., SOA) for galaxy pairs, qualitatively consistent with observations. (b) The SOA is significant for both baryonic and dark matter spins, being the strongest for gas and the weakest for dark matter. (c) The SOA is stronger for less massive targets and for targets having closer neighbors. (d) The SOA strengthens for galaxies in low-density regions, and the signal is dominated by central$-$satellite pairs in low-mass halos. (e) There is an explicit dependence of the SOA on the duration of interaction with its current neighbor. Taken together, we propose that the SOA witnessed at $z = 0$ has been developed mainly by interactions with a neighbor for an extended period of time, rather than tidal torque from the ambient large-scale structure.


title: Early Science from POSSUM: Shocks, turbulence, and a massive new reservoir of ionised gas in the Fornax cluster

abstract: We present the first Faraday rotation measure (RM) grid study of an individual low-mass cluster -- the Fornax cluster -- which is presently undergoing a series of mergers. Exploiting commissioning data for the POlarisation Sky Survey of the Universe's Magnetism (POSSUM) covering a $\sim34$ square degree sky area using the Australian Square Kilometre Array Pathfinder (ASKAP), we achieve an RM grid density of $\sim25$ RMs per square degree from a 280 MHz band centred at 887 MHz, which is similar to expectations for forthcoming GHz-frequency all-sky surveys. We thereby probe the extended magnetoionic structure of the cluster in unprecedented detail. We find that the scatter in the Faraday RM of confirmed background sources is increased by $16.8\pm2.4$ rad m$^{-2}$ within 1 degree (360 kpc) projected distance to the cluster centre, which is 2--4 times more extended than the presently-detectable X-ray-emitting intracluster medium (ICM). The Faraday-active plasma is more massive than the X-ray-emitting ICM, with an average density that broadly matches expectations for the Warm-Hot Intergalactic Medium. The morphology of the Faraday depth enhancement exhibits the classic morphology of an astrophysical bow shock on the southwest side of the main Fornax cluster, and an extended, swept-back wake on the northeastern side. Our favoured explanation is an ongoing merger between the main cluster and a sub-cluster to the southwest. The shock's Mach angle and stand-off distance lead to a self-consistent transonic merger speed with Mach 1.06. The region hosting the Faraday depth enhancement shows a decrement in both total and polarised intensity. We fail to identify a satisfactory explanation for this; further observations are warranted. Generally, our study illustrates the scientific returns that can be expected from all-sky grids of discrete sources generated by forthcoming all-sky radio surveys.


title: The imitation game: Proca stars that can mimic the Schwarzschild shadow

abstract: Can a dynamically robust bosonic star (BS) produce an (effective) shadow that mimics that of a black hole (BH)? The BH shadow is linked to the existence of light rings (LRs). For free bosonic fields, yielding mini-BSs, it is known that these stars can become ultra-compact - i.e., possess LRs - but only for perturbatively unstable solutions. We show this remains the case even when different self-interactions are considered. However, an effective shadow can arise in a different way: if BSs reproduce the existence of an innermost stable circular orbit (ISCO) for timelike geodesics (located at $r_{\rm ISCO}=6M$ for a Schwarzschild BH of mass M), the accretion flow morphology around BHs is mimicked and an effective shadow arises in an astrophysical environment. Even though spherical BSs may accommodate stable timelike circular orbits all the way down to their centre, we show the angular velocity along such orbits may have a maximum away from the origin, at $R_{\Omega}$; this scale was recently observed to mimic the BH's ISCO in some scenarios of accretion flow. Then: (i) for free scalar fields or with quartic self-interactions, $R_{\Omega}\neq 0$ only for perturbatively unstable BSs; (ii) for higher scalar self-interactions, e.g. axionic, $R_{\Omega}\neq 0$ is possible for perturbatively stable BSs, but no solution with $R_{\Omega}=6M$ was found in the parameter space explored; (iii) but for free vector fields, yielding Proca stars (PSs), perturbatively stable solutions with $R_{\Omega}\neq 0$ exist, and indeed $R_{\Omega}=6M$ for a particular solution. Thus, dynamically robust spherical PSs can mimic the shadow of a (near-)equilibrium Schwarzschild BH with the same M, in an astrophysical environment, despite the absence of a LR, at least under some observation conditions, as we confirm by comparing the lensing of such PSs and Schwarzschild BHs.


title: Substructure at High Speed I: Inferring the Escape Velocity in the Presence of Kinematic Substructure

abstract: The local escape velocity provides valuable inputs to the mass profile of the Galaxy, and requires understanding the tail of the stellar speed distribution. Following Leonard $\&$ Tremaine (1990), various works have since modeled the tail of the stellar speed distribution as $\propto (v_{\rm{esc}} -v)^k$, where $v_{\rm{esc}}$ is the escape velocity, and $k$ is the slope of the distribution. In such studies, however, these two parameters were found to be largely degenerate and often a narrow prior is imposed on $k$ in order to constrain $v_{\rm{esc}}$. Furthermore, the validity of the power law form is likely to break down in the presence of multiple kinematic substructures. In this paper, we introduce a strategy that for the first time takes into account the presence of kinematic substructure. We model the tail of the velocity distribution as a sum of multiple power laws without imposing strong priors. Using mock data, we show the robustness of this method in the presence of kinematic structure that is similar to the recently-discovered Gaia Sausage. In a companion paper, we present the new measurement of the escape velocity and subsequently the mass of the Milky Way using Gaia DR2 data.


title: Lifetimes of local excitations in disordered dipolar quantum systems

abstract: When a strongly disordered system of interacting quantum dipoles is locally excited, the excitation relaxes on some (potentially very long) timescale. We analyze this relaxation process, both for electron glasses with strong Coulomb interactions - in which particle-hole dipoles are emergent excitations - and for systems (e.g., quantum magnets or ultracold dipolar molecules) made up of microscopic dipoles. We consider both energy relaxation rates ($T_1$ times) and dephasing rates ($T_2$ times), and their dependence on frequency, temperature, and polarization. Systems in both two and three dimensions are considered, along with the dimensional crossover in quasi-two dimensional geometries. A rich set of scaling laws is found.

categories: cond-mat.dis-nn cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el quant-ph


title: APOGEE discovery of a chemically atypical star disrupted from NGC 6723 and captured by the Milky Way bulge

abstract: The central (`bulge') region of the Milky Way is teeming with a significant fraction of mildly metal-deficient stars with atmospheres that are strongly enriched in cyanogen ($^{12}$C$^{14}$N). Some of these objects, which are also known as nitrogen-enhanced stars, are hypothesised to be relics of the ancient assembly history of the Milky Way. Although the chemical similarity of nitrogen-enhanced stars to the unique chemical patterns observed in globular clusters has been observed, a direct connection between field stars and globular clusters has not yet been proven. In this work, we report on high-resolution, near-infrared spectroscopic observations of the bulge globular cluster NGC 6723, and the serendipitous discovery of a star, 2M18594405$-$3651518, located outside the cluster (near the tidal radius) but moving on a similar orbit, providing the first clear piece of evidence of a star that was very likely once a cluster member and has recently been ejected. Its nitrogen abundance ratio ([N/Fe]$\gtrsim + 0.94$) is well above the typical Galactic field-star levels, and it exhibits noticeable enrichment in the heavy $s$-process elements (Ce, Nd, and Yb), along with moderate carbon enrichment; all characteristics are known examples in globular clusters. This result suggests that some of the nitrogen-enhanced stars in the bulge likely originated from the tidal disruption of globular clusters.


title: Monte Carlo simulations of multiple populations in globular clusters: constraints on the cooling flow vs. accretion scenario using million bodies simulations

abstract: I simulate the evolution of a stellar system hosting two stellar populations whose initial set up is defined according to the two main scenarios proposed for the origin of multiple populations in Galactic globular clusters: (i) formation of a second generation from a cooling flow of pristine+polluted gas and (ii) accretion of polluted gas onto the proto-stellar disks of a fraction of low-mass stars. For this purpose, Monte Carlo simulations containing from $10^{5}$ up to $3\cdot 10^{6}$ particles have been run including the effect of stellar evolution, binary interactions, external tidal field and a detailed modelling of the proto-stellar disk structure. The early accretion of gas onto proto-stellar disks is unable to produce discrete populations and to alter the chemical composition of a significant ($>10\%$) fraction of stars unless a disk lifetime larger ($t_{disk}\sim20~Myr$) than that predicted by models is assumed. Moreover, in this scenario the mixing timescale of the two populations is too short to reproduce the observed segregation of the chemically enriched population. On the other hand, simulations run within the cooling flow scenario can evolve after a Hubble time into stellar systems with a first-to-second population mass ratio similar to that observed in globular clusters, provided that an initial filling-factor $r_{h}/r_{J}>0.15$ is adopted. However, in the weak tidal field regime a radial segregation of the second population stronger than what observed in Milky Way globular clusters at large Galactocentric distances is predicted. This discrepancy disappears in simulations following eccentric orbits in a realistic axisymmetric potential.


title: Gravitational wave cosmology with extreme mass-ratio inspirals

abstract: The Laser Interferometer Space Antenna (LISA) will open the mHz frequency window of the gravitational wave (GW) landscape. Among all the new GW sources expected to emit in this frequency band, extreme mass-ratio inspirals (EMRIs) constitute a unique laboratory for astrophysics and fundamental physics. Here we show that EMRIs can also be used to extract relevant cosmological information, complementary to both electromagnetic (EM) and other GW observations. By using the loudest EMRIs (SNR$>$100) detected by LISA as dark standard sirens, statistically matching their sky localisation region with mock galaxy catalogs, we find that constraints on $H_0$ can reach $\sim$1.1% ($\sim$3.6%) accuracy, at the 90% credible level, in our best (worst) case scenario. By considering a dynamical dark energy (DE) cosmological model, with $\Lambda$CDM parameters fixed by other observations, we further show that in our best (worst) case scenario $\sim$5.9% ($\sim$12.3%) relative uncertainties at the 90% credible level can be obtained on $w_0$, the DE equation of state parameter. Besides being relevant in their own right, EMRI measurements will be affected by different systematics compared to both EM and ground-based GW observations. Cross validation with complementary cosmological measurements will therefore be of paramount importance, especially if convincing evidence of physics beyond $\Lambda$CDM emerges from future observations.


title: The Galactic Faraday rotation sky 2020

abstract: This work gives an update to existing reconstructions of the Galactic Faraday rotation sky by processing almost all Faraday rotation data sets available at the end of the year 2020. Observations of extra-Galactic sources in recent years have, among other regions, further illuminated the previously under-constrained southern celestial sky, as well as parts of the inner disc of the Milky Way. This has culminated in an all-sky data set of 55,190 data points, which is a significant expansion on the 41,330 used in previous works, hence making an updated separation of the Galactic component a promising venture. The increased source density allows us to present our results in a resolution of about $1.3\cdot 10^{-2}\, \mathrm{deg}^2$ ($46.8\,\mathrm{arcmin}^2$), which is a twofold increase compared to previous works. As for previous Faraday rotation sky reconstructions, this work is based on information field theory, a Bayesian inference scheme for field-like quantities which handles noisy and incomplete data. In contrast to previous reconstructions, we find a significantly thinner and pronounced Galactic disc with small-scale structures exceeding values of several thousand $\mathrm{rad}\,\mathrm{m}^{-2}$. The improvements can mainly be attributed to the new catalog of Faraday data, but are also supported by advances in correlation structure modeling within numerical information field theory. We furthermore give a detailed discussion on statistical properties of the Faraday rotation sky and investigate correlations to other data sets.


title: Probing the Halo Gas Distribution in the Inner Galaxy with Fermi Bubble Observations

abstract: The hot halo gas distribution in the inner Milky Way (MW) contains key fossil records of the past energetic feedback processes in the Galactic center. Here we adopt a variety of spherical and disk-like MW halo gas models as initial conditions in a series of simulations to investigate the formation of the Fermi bubbles in the jet-shock scenario. The simulation results are compared directly with relevant X-ray and gamma-ray observations of the Fermi bubbles to constrain the halo gas distribution in the inner Galaxy before the Fermi bubble event. Our best-fit gas density distribution can be described by a power law in radius $n_{\rm e}(r)=0.01(r/1 \text{~kpc})^{-1.5}$ cm$^{-3}$. Our study can not determine if there is an inner density core, which if exists, should be very small with size $r_{c} \lesssim 0.5$ kpc. When extrapolating to large radii $r\sim 50-90$ kpc, our derived density distribution lies appreciably below the recently estimated gas densities from ram-pressure stripping calculations, suggesting that the halo gas density profile either flattens out or has one or more discontinuities within $10 \lesssim r \lesssim 50$ kpc. Some of these discontinuities may be related to the eROSITA bubbles, and our derived gas density profile may correspond to the hot gas distribution in the inner eROSITA bubbles about $5$ Myr ago.


title: Skills-based on technological knowledge in the digital economy activity

abstract: This research seeks to measure the impact of people with technological knowledge on regional digital economic activity and the implications of prosperous cities' contagion effect on neighbouring ones. The focus of this study is quantitative, cross-sectional, and its design is correlational-causal. This study covers seven micro-regions of Minas Gerais in Brazil, organized in 89 municipalities, with 69% urban population and 31% rural. The data used consisted of 4,361 observations obtained in the Brazilian government's public repositories, organized into panel data, and analysed using partial least squares, micro-regional spatial regression, and identification patterns with machine learning. The confirmatory analysis of the regression test establishes a significant impact between the CE's technological knowledge and the digital economic activity AED through a predictive value of R2 = .749, \b{eta} = .867, p = .000 (value t = 18,298). With high notoriety among the variables, public and private university institutions (IUPP), professors with doctorates and masters (DCNT), and information technology occupations (CBO). A geographic concentration of companies that demand technology-based skills had effects by slowing down the development of small municipalities, suggesting the development of new government technology initiatives that support new business models based on technological knowledge.


title: An abstract theory of physical measurements

abstract: The question of what should be meant by a measurement is tackled from a mathematical perspective whose physical interpretation is that a measurement is a fundamental process via which a finite amount of classical information is produced. This translates into an algebraic and topological definition of measurement space that caters for the distinction between quantum and classical measurements and allows a notion of observer to be derived.


title: On the relation between Lya absorbers and local galaxy filaments

abstract: The intergalactic medium (IGM) is believed to contain the majority of baryons in the universe and to trace the same dark matter structure as galaxies, forming filaments and sheets. Lya absorbers, which sample the neutral component of the IGM, have been extensively studied at low and high redshift, but the exact relation between Lya absorption, galaxies and the large-scale structure is observationally not well-constrained. In this study, we aim at characterising the relation between Lya absorbers and nearby overdense cosmological structures (galaxy filaments) at recession velocities Delta v \leq 6700 km/s by using archival observational data from various instruments. We analyse 587 intervening Lya absorbers in the spectra of 302 extragalactic background sources obtained with the COS installed on the HST. We combine the absorption-line information with galaxy data of five local galaxy filaments originally mapped by Courtois et al. (2013). Along the 91 sightlines that pass close to a filament, we identify 215 (227) Lya absorption systems (components). Among these, 74 Lya systems are aligned in position and velocity with the galaxy filaments, indicating that these absorbers and the galaxies trace the same large-scale structure. The filament-aligned Lya absorbers have a ~90 percent higher rate of incidence (dN/dz=189 for log N HI \geq 13.2) and a mildly shallower slope (-beta = -1.47) of the column density distribution function than the general Lya population at z=0, reflecting the filaments' matter overdensity. The strongest Lya absorbers are preferentially found near galaxies or close to the axis of a filament, although there is substantial scatter in this relation. Our sample of absorbers clusters more strongly around filament axes than a randomly distributed sample would do (as confirmed by a KS test), but the clustering signal is less pronounced than for the galaxies in the filaments.


title: Implications for First-Order Cosmological Phase Transitions from the Third LIGO-Virgo Observing Run

abstract: We place constraints on the normalized energy density in gravitational waves from first-order strong phase transitions using data from Advanced LIGO and Virgo's first, second and third observing runs. First, adopting a broken power law model, we place $95 \%$ confidence level upper limits simultaneously on the gravitational-wave energy density at 25 Hz from unresolved compact binary mergers, $\Omega_{\rm CBC} < 6.1 \times 10^{-9}$, and strong first-order phase transitions, $\Omega_{\rm BPL} < 4.4 \times 10^{-9}$. The inclusion of the former is necessary since we expect this astrophysical signal to be the foreground of any detected spectrum. We then consider two more complex phenomenological models, limiting at 25 Hz the gravitational-wave background due to bubble collisions to $\Omega_{\rm pt} < 5.0\times 10^{-9}$ and the background due to sound waves to $\Omega_{\rm pt} < 5.8\times10^{-9}$ at $95 \%$ confidence level for phase transitions occurring at temperatures above $10^8$GeV.


title: The Demographics of Wide-Separation Planets

abstract: I begin this review by first defining what is meant by exoplanet demographics, and then motivating why we would like as broad a picture of exoplanet demographics as possible. I then outline the methodology and pitfalls to measuring exoplanet demographics in practice. I next review the methods of detecting exoplanets, focusing on the ability of these methods to detect wide separation planets. For the purposes of this review, I define wide separation as separations beyond the 'snow line' of the protoplanetary disk, which is at roughly $\sim3$ au for a sunlike star. I note that this definition is somewhat arbitrary, and the practical boundary depends on the host star mass, planet mass and radius, and detection method. I review the approximate scaling relations for the signal-to-noise ratio for the detectability of exoplanets as a function of the relevant physical parameters, including the host star properties. I provide abroad overview of what has already been learned from the transit, radial velocity, direct imaging, and microlensing methods. I outline the challenges to synthesizing the demographics using different methods and discuss some preliminary first steps in this direction. Finally, I describe future prospects for providing a nearly complete statistical census of exoplanets.


title: A survey of some recent applications of optimal transport methods to econometrics

abstract: This paper surveys recent applications of methods from the theory of optimal transport to econometric problems.


title: Dynamical Evolution of Multiple-Population Globular Clusters

abstract: We have carried out a set of Monte Carlo simulations to study a number of fundamental aspects of the dynamical evolution of multiple stellar populations in globular clusters with different initial masses, fractions of second generation (2G) stars, and structural properties. Our simulations explore and elucidate: 1) the role of early and long-term dynamical processes and stellar escape in the evolution of the fraction of 2G stars and the link between the evolution of the fraction of 2G stars and various dynamical parameters; 2) the link between the fraction of 2G stars inside the cluster and in the population of escaping stars during a cluster's dynamical evolution; 3) the dynamics of the spatial mixing of the first-generation (1G) and 2G stars and the details of the structural properties of the two populations as they evolve toward mixing; 4) the implications of the initial differences between the spatial distribution of 1G and 2G stars for the evolution of the anisotropy in the velocity distribution and the expected radial profile of the 1G and 2G anisotropy for clusters at different stages of their dynamical history; 5) the variation of the degree of energy equipartition of the 1G and the 2G populations as a function of the distance from the cluster's centre and the cluster's evolutionary phase.


title: Archimedes' principle for ideal gas

abstract: We prove Archimedes' principle for a macroscopic ball in ideal gas consisting of point particles with non-zero mass. The main result is an asymptotic theorem, as the number of point particles goes to infinity and their total mass remains constant. We also show that, asymptotically, the gas has an exponential density as a function of height. We find the asymptotic inverse temperature of the gas. We derive an accurate estimate of the volume of the phase space using the local central limit theorem.


title: Initial condition assessment for reaction-diffusion glioma growth models: A translational MRI/histology (in)validation study

abstract: Diffuse gliomas are highly infiltrative tumors whose early diagnosis and follow-up usually rely on magnetic resonance imaging (MRI). However, the limited sensitivity of this technique makes it impossible to directly assess the extent of the glioma cell invasion, leading to sub-optimal treatment planing. Reaction-diffusion growth models have been proposed for decades to extrapolate glioma cell infiltration beyond margins visible on MRI and predict its spatial-temporal evolution. These models nevertheless require an initial condition, that is the tumor cell density values at every location of the brain at diagnosis time. Several works have proposed to relate the tumor cell density function to abnormality outlines visible on MRI but the underlying assumptions have never been verified so far. In this work we propose to verify these assumptions by stereotactic histological analysis of a non-operated brain with glioblastoma using a tailored 3D-printed slicer. Cell density maps are computed from histological slides using a deep learning approach. The density maps are then registered to a postmortem MR image and related to an MR-derived geodesic distance map to the tumor core. The relation between the edema outlines visible on T2 FLAIR MRI and the distance to the core is also investigated. Our results suggest that (i) the previously suggested exponential decrease of the tumor cell density with the distance to the tumor core is not unreasonable but (ii) the edema outlines may in general not correspond to a cell density iso-contour and (iii) the commonly adopted tumor cell density value at these outlines is likely overestimated. These findings highlight the limitations of using conventional MRI to derive glioma cell density maps and point out the need of validating other methods to initialize reaction-diffusion growth models and make them usable in clinical practice.

categories: physics.med-ph cs.LG eess.IV physics.bio-ph


title: A tidally tilted sectoral dipole pulsation mode in the eclipsing binary TIC 63328020

abstract: We report the discovery of the third tidally tilted pulsator, TIC 63328020. Observations with the TESS satellite reveal binary eclipses with an orbital period of 1.1057 d, and $\delta$ Scuti-type pulsations with a mode frequency of 21.09533 d$^{-1}$. This pulsation exhibits a septuplet of orbital sidelobes as well as a harmonic quintuplet. Using the oblique pulsator model, the primary oscillation is identified as a sectoral dipole mode with $l = 1, |m| = 1$. We find the pulsating star to have $M_1 \simeq 2.5\, {\rm M}_\odot$, $R_1 \simeq 3 \, {\rm R}_\odot$, and $T_{\rm eff,1} \simeq 8000$ K, while the secondary has $M_2 \simeq 1.1 \, {\rm M}_\odot$, $R_2 \simeq 2 \, {\rm R}_\odot$, and $T_{\rm eff,2} \simeq 5600$ K. Both stars appear to be close to filling their respective Roche lobes. The properties of this binary as well as the tidally tilted pulsations differ from the previous two tidally tilted pulsators, HD74423 and CO Cam, in important ways. We also study the prior history of this system with binary evolution models and conclude that extensive mass transfer has occurred from the current secondary to the primary.


title: Direct numerical simulation of quasi-two-dimensional MHD turbulent shear flows

abstract: Direct numerical simulations (DNS) are performed to study the turbulent shear flow of an electrically conducting fluid in a cylindrical container. The flow is driven by the interaction between the radial electric currents ($I$) injected through a large number of small electrodes at the bottom wall and an axial magnetic field. All the numerical parameters, including the geometry of the container, the total injcected currents and the magnetic field, are in line with the experiment performed in J. Fluid Mech. 456, 137-159. First, witth laminar Hartmann layers, three dimensional simulations recover experimentally measured quantities (global angular momentum, velocity profiles). The variation laws for the wall shear stresses, the energy spectra and visualizations of flow structures near the side wall highlight separation and turbulence within the side wall layers. Furthermore, a parametric analysis of the flow reveals that Ekman recirculations have significant influence on the vortex size, the free shear layer, and the global dissipation. Second, we recover the scaling laws of the cutoff scale that separate the large quasi-two-dimensional scales from the small three-dimensional ones (J. Fluid Mech. 118, 507-518), and thus establish their validity in sheared MHD turbulence. Furthermore, we find that three-componentality are and the three-dimensionality appear concurrently and that both the two-dimensional cutoff frequency and the mean energy associated to the axial component of velocity scale with $N_t$, respectively as $0.063N_t^{0.37}$ and $0.126 N_t^{-0.92}$.


title: A Historical and Statistical Studyof the Software Vulnerability Landscape

abstract: Understanding the landscape of software vulnerabilities is key for developing effective security solutions. Fortunately, the evaluation of vulnerability databases that use a framework for communicating vulnerability attributes and their severity scores, such as the Common Vulnerability Scoring System (CVSS), can help shed light on the nature of publicly published vulnerabilities. In this paper, we characterize the software vulnerability landscape by performing a historical and statistical analysis of CVSS vulnerability metrics over the period of 2005 to 2019 through using data from the National Vulnerability Database. We conduct three studies analyzing the following: the distribution of CVSS scores (both empirical and theoretical), the distribution of CVSS metric values and how vulnerability characteristics change over time, and the relative rankings of the most frequent metric value over time. Our resulting analysis shows that the vulnerability threat landscape has been dominated by only a few vulnerability types and has changed little during the time period of the study. The overwhelming majority of vulnerabilities are exploitable over the network. The complexity to successfully exploit these vulnerabilities is dominantly low; very little authentication to the target victim is necessary for a successful attack. And most of the flaws require very limited interaction with users. However on the positive side, the damage of these vulnerabilities is mostly confined within the security scope of the impacted components. A discussion of lessons that could be learned from this analysis is presented.


title: Apollo: Transferable Architecture Exploration

abstract: The looming end of Moore's Law and ascending use of deep learning drives the design of custom accelerators that are optimized for specific neural architectures. Architecture exploration for such accelerators forms a challenging constrained optimization problem over a complex, high-dimensional, and structured input space with a costly to evaluate objective function. Existing approaches for accelerator design are sample-inefficient and do not transfer knowledge between related optimizations tasks with different design constraints, such as area and/or latency budget, or neural architecture configurations. In this work, we propose a transferable architecture exploration framework, dubbed Apollo, that leverages recent advances in black-box function optimization for sample-efficient accelerator design. We use this framework to optimize accelerator configurations of a diverse set of neural architectures with alternative design constraints. We show that our framework finds high reward design configurations (up to 24.6% speedup) more sample-efficiently than a baseline black-box optimization approach. We further show that by transferring knowledge between target architectures with different design constraints, Apollo is able to find optimal configurations faster and often with better objective value (up to 25% improvements). This encouraging outcome portrays a promising path forward to facilitate generating higher quality accelerators.