title: Weakly-Supervised Deep Learning for Domain Invariant Sentiment Classification

abstract: The task of learning a sentiment classification model that adapts well to any target domain, different from the source domain, is a challenging problem. Majority of the existing approaches focus on learning a common representation by leveraging both source and target data during training. In this paper, we introduce a two-stage training procedure that leverages weakly supervised datasets for developing simple lift-and-shift-based predictive models without being exposed to the target domain during the training phase. Experimental results show that transfer with weak supervision from a source domain to various target domains provides performance very close to that obtained via supervised training on the target domain itself.


title: Turbulence -- Obstacle Interactions in the Lagrangian Framework: Applications for Stochastic Modeling in Canopy Flows

abstract: Lagrangian stochastic models are widely used to predict and analyze turbulent dispersion in complex environments, such as in various terrestrial and marine canopy flows. However, due to a lack of empirical data, it is still not understood how particular features of highly inhomogeneous canopy flows affect the Lagrangian statistics. In this work, we study Lagrangian short time statistics by analyzing empirical Lagrangian trajectories in sub-volumes of space that are small in comparison with the canopy height. For the analysis we used 3D Lagrangian trajectories measured in a dense canopy flow model in a wind-tunnel, using an extended version of real-time 3D particle tracking velocimetry (3D-PTV). One of our key results is that the random turbulent fluctuations due to the intense dissipation were more dominant than the flow's inhomogeneity in affecting the short-time Lagrangian statistics. This amounts to a so-called quasi-homogeneous regime of Lagrangian statistics at small scales. Using the Lagrangian dataset we calculate the Lagrangian autocorrelation function and the second-order Lagrangian structure-function, and extract associated parameters, namely a Lagrangian velocity decorrelation timescale, $T_i$, and the Kolmogorov constant, $C_0$. We demonstrate that in the quasi-homogeneous regime, both these functions are well represented using a second-order Lagrangian stochastic model that was designed for homogeneous flows. Furthermore, we show that the spatial variations of the Lagrangian separation of scales, $T_i/\tau_\eta$, and the Kolmogorov constant, $C_0$, cannot be explained by the variation of the Reynolds number, $Re_\lambda$, in space, and that $T_i/\tau_\eta$ was small as compared with homogeneous turbulence predictions at similar $Re_\lambda$. We thus hypothesize that this occurred due to the so-called "wake production", and show empirical results supporting our hypothesis.


title: Distribution Density, Tails, and Outliers in Machine Learning: Metrics and Applications

abstract: We develop techniques to quantify the degree to which a given (training or testing) example is an outlier in the underlying distribution. We evaluate five methods to score examples in a dataset by how well-represented the examples are, for different plausible definitions of "well-represented", and apply these to four common datasets: MNIST, Fashion-MNIST, CIFAR-10, and ImageNet. Despite being independent approaches, we find all five are highly correlated, suggesting that the notion of being well-represented can be quantified. Among other uses, we find these methods can be combined to identify (a) prototypical examples (that match human expectations); (b) memorized training examples; and, (c) uncommon submodes of the dataset. Further, we show how we can utilize our metrics to determine an improved ordering for curriculum learning, and impact adversarial robustness. We release all metric values on training and test sets we studied.


title: Minimal Radius Enclosing Polyellipsoids

abstract: In this paper we analyze the extension of the classical smallest enclosing disk problem to the case of the location of a polyellipsoid to fully cover a set of demand points in $\mathbb{R}^d$. We prove that the problem is polynomially solvable in fixed dimension and analyze mathematical programming formulations for it. We also consider some geometric approaches for the problem in case the foci of the polyellipsoids are known. Extensions of the classical algorithm by Elzinga-Hearn are also derived for this new problem. Moreover, we address several extensions of the problem, as the case where the foci of the enclosing polyellipsoid are not given and have to be determined among a potential set of points or the induced covering problems when instead of polyellipsoids, one uses ordered median polyellipsoids. For these problems we also present Mixed Integer (Non) Linear Programming strategies that lead to efficient ways to solve it.


title: zfit: scalable pythonic fitting

abstract: Statistical modeling is a key element in many scientific fields and especially in High-Energy Physics (HEP) analysis. The standard framework to perform this task in HEP is the C++ ROOT/RooFit toolkit; with Python bindings that are only loosely integrated into the scientific Python ecosystem. In this paper, zfit, a new alternative to RooFit written in pure Python, is presented. Most of all, zfit provides a well defined high-level API and workflow for advanced model building and fitting, together with an implementation on top of TensorFlow, allowing a transparent usage of CPUs and GPUs. It is designed to be extendable in a very simple fashion, allowing the usage of cutting-edge developments from the scientific Python ecosystem in a transparent way. The main features of zfit are introduced, and its extension to data analysis, especially in the context of HEP experiments, is discussed.


title: Higgs vacuum (in)stability during inflation: the dangerous relevance of de Sitter departure and Planck-suppressed operators

abstract: The measured Standard Model parameters lie in a range such that the Higgs potential, once extrapolated up to high scales, develops a minimum of negative energy density. This has important cosmological implications. In particular, during inflation, quantum fluctuations could have pushed the Higgs field beyond its potential barrier, triggering the formation of anti-de Sitter regions, with fatal consequences for our universe. By requiring that this did not happen, one can in principle connect (and constrain) Standard Model parameters with the energy scale of inflation. In this context, we highlight the sensitivity of the fate of our vacuum to seemingly irrelevant physics. In particular, the departure of inflation from an exact de Sitter phase, as well as Planck-suppressed derivative operators, can, already and surprisingly, play a decisive role in (de)stabilizing the Higgs during inflation. Furthermore, in the stochastic dynamics, we quantify the impact of the amplitude of the noise differing from the one of a massless field, as well as of going beyond the slow-roll approximation by using a phase-space approach. On a general ground, our analysis shows that relating the period of inflation to precision particle physics requires a knowledge of these "irrelevant" effects.


title: Integral characterization for Poincar\'e half-maps in planar linear systems

abstract: The intrinsic nature of a problem usually suggests a first suitable method to deal with it. Unfortunately, the apparent ease of application of these initial approaches may make their possible flaws seem to be inherent to the problem and often no alternative ways to solve it are searched for. For instance, since linear systems of differential equations are easy to integrate, Poincar\'e half-maps for piecewise linear systems are always studied by using the direct integration of the system in each zone of linearity. However, this approach is accompanied by two important defects: due to the different spectra of the involved matrices, many cases and strategies must be considered and, since the flight time appears as a new variable, nonlinear complicated equations arise. This manuscript is devoted to present a novel theory to characterize Poincar\'e half-maps in planar linear systems that avoids the computation of their solutions and the problems it causes. This new perspective rests on the use of line integrals of a specific conservative vector field which is orthogonal to the flow of the linear system. Besides the obvious mathematical interest, this approach is attractive because it allows to simplify the study of piecewise-linear systems and deal with open problems in this field.


title: Longitudinal short-distance constraints for the hadronic light-by-light contribution to $(g-2)_\mu$ with large-$N_c$ Regge models

abstract: While the low-energy part of the hadronic light-by-light (HLbL) tensor can be constrained from data using dispersion relations, for a full evaluation of its contribution to the anomalous magnetic moment of the muon $(g-2)_\mu$ also mixed- and high-energy regions need to be estimated. Both can be addressed within the operator product expansion (OPE), either for configurations where all photon virtualities become large or one of them remains finite. Imposing such short-distance constraints (SDCs) on the HLbL tensor is thus a major aspect of a model-independent approach towards HLbL scattering. Here, we focus on longitudinal SDCs, which concern the amplitudes containing the pseudoscalar-pole contributions from $\pi^0$, $\eta$, $\eta'$. Since these conditions cannot be fulfilled by a finite number of pseudoscalar poles, we consider a tower of excited pseudoscalars, constraining their masses and transition form factors from Regge theory, the OPE, and phenomenology. Implementing a matching of the resulting expressions for the HLbL tensor onto the perturbative QCD quark loop, we are able to further constrain our calculation and significantly reduce its model dependence. We find that especially for the $\pi^0$ the corresponding increase of the HLbL contribution is much smaller than previous prescriptions in the literature would imply. Overall, we estimate that longitudinal SDCs increase the HLbL contribution by $\Delta a_\mu^\text{LSDC}=13(6)\times 10^{-11}$. This number does not include the contribution from the charm quark, for which we find $a_\mu^{c\text{-quark}}=3(1)\times 10^{-11}$.


title: Thresholds versus fractional expectation-thresholds

abstract: Proving a conjecture of Talagrand, a fractional version of the 'expectation-threshold' conjecture of Kalai and the second author, we show for any increasing family $F$ on a finite set $X$ that $p_c (F) =O( q_f (F) \log \ell(F))$, where $p_c(F)$ and $q_f(F)$ are the threshold and 'fractional expectation-threshold' of $F$, and $\ell(F)$ is the largest size of a minimal member of $F$. This easily implies several heretofore difficult results and conjectures in probabilistic combinatorics, including thresholds for perfect hypergraph matchings (Johansson--Kahn--Vu), bounded-degree spanning trees (Montgomery), and bounded-degree spanning graphs (new). We also resolve (and vastly extend) the 'axial' version of the random multi-dimensional assignment problem (earlier considered by Martin--M\'{e}zard--Rivoire and Frieze--Sorkin). Our approach builds on a recent breakthrough of Alweiss, Lovett, Wu and Zhang on the Erd\H{o}s--Rado 'Sunflower Conjecture'.


title: Experimental validation of the $1/\tau$ -scaling entropy generation in finite-time thermodynamics with dry air

abstract: The second law of thermodynamics can be described as the non-decreasing of the entropy in the irreversible thermodynamic process. Such phenomenon can be quantitatively evaluated with the irreversible entropy generation (IEG), which was recently found to follow a $1/\tau$ scaling for the system under a long contact time $\tau$ with the thermal bath. This scaling, predicted in many finite-time thermodynamic models, is of great potential in the optimization of heat engines, yet remains lack of direct experimental validation. In this letter, we design an experimental apparatus to test such scaling by compressing dry air in a temperature-controlled water bath. More importantly, we quantitatively verify the optimized control protocol to reduce the IEG. Such optimization shall bring new insight to the practical design of heat engine cycles.


title: An Ultimate Approach of Mitigating Attacks in RPL Based Low Power Lossy Networks

abstract: The Routing Protocol for Low-Power and Lossy Networks (RPL) is the existing routing protocol for Internet of Things (IoT). RPL is a proactive,lightweight, Distance Vector protocol which offers security against various forms of routing attacks. Still, there are various attacks(as rank, version attacks and many more ) which is possible in this network due to problem of unauthenticated or unencrypted control frames, centralized root controller, compromised or unauthenticated devices and many more ways. There are various solutions present in the literature but every solution has its pros and cons. There is no appropriate system framework till now which completely solves these all issues. So, we present an ultimate approach to mitigate these RPL attacks more efficiently and effectively. We use IDS based system for internal attacks and a mini-firewall for removing the external attacks. In IDS based approach, we use intrusion detection system at multiple locations for analyzing the behaviour of nodes. The final decision whether the node is attacker or not depends on mainly three things as: trust between the neighbouring nodes, local decision by multiple sink nodes and global decision by root node. We also use some blockchain features in this framework for better internal security. We use some threshold values and rules in mini-firewall for removing external attacks. In this paper, we provide the proposed approach and theoretical analysis of this approach which provide better protection from these attacks than any other method.


title: LISA and the Existence of a Fast-Merging Double Neutron Star Formation Channel

abstract: Using a Milky Way double neutron star (DNS) merger rate of 210 Myr$^{-1}$, as derived by the Laser Interferometer Gravitational-Wave Observatory (LIGO), we demonstrate that the Laser Interferometer Space Antenna (LISA) will detect on average 240 (330) DNSs within the Milky Way for a 4-year (8-year) mission with a signal-to-noise ratio greater than 7. Even adopting a more pessimistic rate of 42 Myr$^{-1}$, as derived by the population of Galactic DNSs, we find a significant detection of 46 (65) Milky Way DNSs. These DNSs can be leveraged to constrain formation scenarios. In particular, traditional NS-discovery methods using radio telescopes are unable to detect DNSs with $P_{\rm orb}$ $\lesssim$1 hour (merger times $\lesssim$10 Myr). If a fast-merging channel exists that forms DNSs at these short orbital periods, LISA affords, perhaps, the only opportunity to observationally characterize these systems; we show that toy models for possible formation scenarios leave unique imprints on DNS orbital eccentricities, which may be measured by LISA for values as small as $\sim$10$^{-2}$.


title: An Empirical Study of Generation Order for Machine Translation

abstract: In this work, we present an empirical study of generation order for machine translation. Building on recent advances in insertion-based modeling, we first introduce a soft order-reward framework that enables us to train models to follow arbitrary oracle generation policies. We then make use of this framework to explore a large variety of generation orders, including uninformed orders, location-based orders, frequency-based orders, content-based orders, and model-based orders. Curiously, we find that for the WMT'14 English $\to$ German translation task, order does not have a substantial impact on output quality, with unintuitive orderings such as alphabetical and shortest-first matching the performance of a standard Transformer. This demonstrates that traditional left-to-right generation is not strictly necessary to achieve high performance. On the other hand, results on the WMT'18 English $\to$ Chinese task tend to vary more widely, suggesting that translation for less well-aligned language pairs may be more sensitive to generation order.


title: Compact almost automorphic solutions for semilinear parabolic evolution equations

abstract: In this paper, using the subvariant functional method due to Favard \cite{Favard}, we prove the existence of aunique compact almost automorphic solution for a class of semilinear evolution equations in Banach spaces. More specifically, we improve the assumptions in \cite{CieuEzz}, we show that the almost automorphy of the coefficients in a weaker sense (Stepanov almost automorphy of order $1\leq p <\infty$) is enough to obtain solutions that are almost automorphic in a strong sense (Bochner almost automorphy). We distinguish two cases, $ p=1 $ and $ p>1$. Moreover, we propose to study a class of reaction-diffusion problems.


title: Learning to Manipulate Deformable Objects without Demonstrations

abstract: In this paper we tackle the problem of deformable object manipulation through model-free visual reinforcement learning (RL). In order to circumvent the sample inefficiency of RL, we propose two key ideas that accelerate learning. First, we propose an iterative pick-place action space that encodes the conditional relationship between picking and placing on deformable objects. The explicit structural encoding enables faster learning under complex object dynamics. Second, instead of jointly learning both the pick and the place locations, we only explicitly learn the placing policy conditioned on random pick points. Then, by selecting the pick point that has Maximal Value under Placing (MVP), we obtain our picking policy. This provides us with an informed picking policy during testing, while using only random pick points during training. Experimentally, this learning framework obtains an order of magnitude faster learning compared to independent action-spaces on our suite of deformable object manipulation tasks with visual RGB observations. Finally, using domain randomization, we transfer our policies to a real PR2 robot for challenging cloth and rope coverage tasks, and demonstrate significant improvements over standard RL techniques on average coverage.


title: Identifying Exo-Earth Candidates in Direct Imaging Data through Bayesian Classification

abstract: Future space telescopes may be able to directly image $\sim$10 - 100 planets with sizes and orbits consistent with habitable surface conditions ("exo-Earth candidates" or EECs), but observers will face difficulty in distinguishing these from the potentially hundreds of non-habitable "false positives" which will also be detected. To maximize the efficiency of follow-up observations, a prioritization scheme must be developed to determine which planets are most likely to be EECs. In this paper, we present a Bayesian method for estimating the likelihood that any directly imaged extrasolar planet is a true exo-Earth candidate by interpreting the planet's apparent magnitude and separation in light of existing exoplanet statistics. As a specific application of this general framework, we use published estimates of the discovery yield of future space-based direct imaging mission concepts to conduct "mock surveys" in which we compute the likelihood that each detected planet is an EEC. We find that it will be difficult to determine which planets are EECs with $>50\%$ confidence using single-band photometry immediately upon their detection. The best way to reduce this ambiguity would be to constrain the planet's orbit by revisiting the system multiple times or through a radial velocity precursor survey. Astrometric or radial velocity constraints on the planet's mass would offer a lesser benefit. Finally, we show that a Bayesian approach to prioritizing targets would improve the follow-up efficiency of a direct imaging survey versus a blind approach using the same data. For example, the prioritized approach could reduce the amount of integration time required for the spectral detection (or rejection) of water absorption in most EECs by a factor of two.


title: Electric-magnetic duality in the quantum double models of topological orders with gapped boundaries

abstract: We generalize the Electric-magnetic (EM) duality in the quantum double (QD) models to the case of topological orders with gapped boundaries. We also map the QD models with boundaries to the Levin-Wen (LW) models with boundaries. To this end, we Fourier transform and rewrite the extended QD model with a finite gauge group $G$ on a trivalent lattice with a boundary. Gapped boundary conditions of the model before the transformation are known to be characterized by the subgroups $K \subseteq G$. We find that after the transformation, the boundary conditions are then characterized by the Frobenius algebras $A_{G,K}$ in $\mathrm{Rep}_G$. An $A_{G,K}$ is the dual space of the quotient of the group algebra of $G$ over that of $K$, and $\mathrm{Rep}_G$ is the category of the representations of $G$. The EM duality on the boundary is revealed by mapping the $K$'s to $A_{G,K}$'s. We also show that our transformed extended QD model can be mapped to an extended LW model on the same lattice via enlarging the Hilbert space of the extended LW model. Moreover, our transformed extended QD model elucidates the phenomenon of anyon splitting in anyon condensation.


title: Dual Separated Variables and Scalar Products

abstract: Separation of variables (SoV) is an extremely efficient and elegant technique for analysing physical systems but its application to integrable spin chains was limited until recently to the simplest su(2) cases. In this paper we continue developing the SoV program for higher-rank spin chains and demonstrate how to derive the measure for the su(3) case. Our results are a natural consequence of factorisability of the wave function and functional orthogonality relations following from the interplay between Baxter equations for Q-functions and their dual.

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


title: Multilevel evolutionary developmental optimization (MEDO): A theoretical framework for understanding preferences and selection dynamics

abstract: What is motivation and how does it work? Where do goals come from and how do they vary within and between species and individuals? Why do we prefer some things over others? MEDO is a theoretical framework for understanding these questions in abstract terms, as well as for generating and evaluating specific hypotheses that seek to explain goal-oriented behavior. MEDO views preferences as selective pressures influencing the likelihood of particular outcomes. With respect to biological organisms, these patterns must compete and cooperate in shaping system evolution. To the extent that shaping processes are themselves altered by experience, this enables feedback relationships where histories of reward and punishment can impact future motivation. In this way, various biases can undergo either amplification or attenuation, resulting in preferences and behavioral orientations of varying degrees of inter-temporal and inter-situational stability. MEDO specifically models all shaping dynamics in terms of natural selection operating on multiple levels--genetic, neural, and cultural--and even considers aspects of development to themselves be evolutionary processes. Thus, MEDO reflects a kind of generalized Darwinism, in that it assumes that natural selection provides a common principle for understanding the emergence of complexity within all dynamical systems in which replication, variation, and selection occur. However, MEDO combines this evolutionary perspective with economic decision theory, which describes both the preferences underlying individual choices, as well as the preferences underlying choices made by engineers in designing optimized systems. In this way, MEDO uses economic decision theory to describe goal-oriented behaviors as well as the interacting evolutionary optimization processes from which they emerge. (Please note: this manuscript was written and finalized in 2012.)


title: Highly-scalable, physics-informed GANs for learning solutions of stochastic PDEs

abstract: Uncertainty quantification for forward and inverse problems is a central challenge across physical and biomedical disciplines. We address this challenge for the problem of modeling subsurface flow at the Hanford Site by combining stochastic computational models with observational data using physics-informed GAN models. The geographic extent, spatial heterogeneity, and multiple correlation length scales of the Hanford Site require training a computationally intensive GAN model to thousands of dimensions. We develop a hierarchical scheme for exploiting domain parallelism, map discriminators and generators to multiple GPUs, and employ efficient communication schemes to ensure training stability and convergence. We developed a highly optimized implementation of this scheme that scales to 27,500 NVIDIA Volta GPUs and 4584 nodes on the Summit supercomputer with a 93.1% scaling efficiency, achieving peak and sustained half-precision rates of 1228 PF/s and 1207 PF/s.


title: G2SAT: Learning to Generate SAT Formulas

abstract: The Boolean Satisfiability (SAT) problem is the canonical NP-complete problem and is fundamental to computer science, with a wide array of applications in planning, verification, and theorem proving. Developing and evaluating practical SAT solvers relies on extensive empirical testing on a set of real-world benchmark formulas. However, the availability of such real-world SAT formulas is limited. While these benchmark formulas can be augmented with synthetically generated ones, existing approaches for doing so are heavily hand-crafted and fail to simultaneously capture a wide range of characteristics exhibited by real-world SAT instances. In this work, we present G2SAT, the first deep generative framework that learns to generate SAT formulas from a given set of input formulas. Our key insight is that SAT formulas can be transformed into latent bipartite graph representations which we model using a specialized deep generative neural network. We show that G2SAT can generate SAT formulas that closely resemble given real-world SAT instances, as measured by both graph metrics and SAT solver behavior. Further, we show that our synthetic SAT formulas could be used to improve SAT solver performance on real-world benchmarks, which opens up new opportunities for the continued development of SAT solvers and a deeper understanding of their performance.


title: Proactive Optimization with Machine Learning: Femto-caching with Future Content Popularity

abstract: Optimizing resource allocation with predicted information has shown promising gain in boosting network performance and improving user experience. Earlier research efforts focus on optimizing proactive policies under the assumption of knowing the future information. Recently, various techniques have been proposed to predict the required information, and the prediction results were then treated as the true value in the optimization, i.e., "first-predict-then-optimize". In this paper, we introduce a proactive optimization framework for anticipatory resource allocation, where the future information is implicitly predicted under the same objective with the policy optimization in a single step. An optimization problem is formulated to integrate the implicit prediction and the policy optimization, based on the conditional distribution of the future information given the historical observations. To solve such a problem, we transform it equivalently to a problem depending on the joint distribution of future and historical information. Then, we resort to unsupervised learning with neural networks to learn the proactive policy as a function of the past observations via stochastic optimization. We take proactive caching and bandwidth allocation at base stations as a concrete example, where the objective function is the conditional expectation of successful offloading probability taken over the future popularity given the historically observed popularity. We use simulation to validate the proposed framework and compare it with the "first-predict-then-optimize" strategy and a heuristic "end-to-end" optimization strategy with supervised learning.


title: Entanglement in coupled kicked tops with chaotic dynamics

abstract: The entanglement of eigenstates in two coupled, classically chaotic kicked tops is studied in dependence of their interaction strength. The transition from the non-interacting and unentangled system towards full random matrix behavior is governed by a universal scaling parameter. Using suitable random matrix transition ensembles we express this transition parameter as a function of the subsystem sizes and the coupling strength for both unitary and orthogonal symmetry classes. The universality is confirmed for the level spacing statistics of the coupled kicked tops and a perturbative description is in good agreement with numerical results. The statistics of Schmidt eigenvalues and entanglement entropies of eigenstates is found to follow a universal scaling as well. Remarkably this is not only the case for large subsystems of equal size but also if one of them is much smaller. For the entanglement entropies a perturbative description is obtained, which can be extended to large couplings and provides very good agreement with numerical results. Furthermore, the transition of the statistics of the entanglement spectrum towards the random matrix limit is demonstrated for different ratios of the subsystem sizes.


title: Energy-momentum tensor and helicity for gauge fields coupled to a pseudo-scalar inflaton

abstract: We study the energy-momentum tensor and helicity of gauge fields coupled through $g \phi F \tilde{F}/4$ to a pseudo-scalar field $\phi$ driving inflation. Under the assumption of a constant time derivative of the background inflaton, we compute analitically divergent and finite terms of the energy density and helicity of gauge fields for any value of the coupling $g$. We introduce a suitable adiabatic expansion for mode functions of physical states of the gauge fields which correctly reproduces ultraviolet divergences in average quantities and identify corresponding counterterms. Our calculations shed light on the accuracy and the range of validity of approximated analytic estimates of the energy density and helicity terms previously existed in the literature in the strongly coupled regime only, i.e. for $g \dot \phi/(2H) \gg 1$. We discuss the implications of our analytic calculations for the backreaction of quantum fluctuations onto the inflaton evolution.


title: The first steps of Interstellar Phosphorus Chemistry

abstract: Phosphorus-bearing species are an essential key to form life on Earth, however they have barely been detected in the interstellar medium. Since only PN and PO have been identified so far towards star-forming regions, the chemical formation pathways of P-bearing molecules are not easy to constrain and are thus highly debatable. An important factor still missing in the chemical models is the initial elemental abundance of phosphorus, i.e. the depletion level of P. In order to overcome this problem, we study P-bearing species in diffuse/translucent clouds. In these objects phosphorus is mainly in the gas phase and therefore the elemental initial abundance needed in our chemical simulations corresponds to the cosmic one and is thus well constrained. An advanced chemical model with an updated P-chemistry network has been used. Single-pointing observations were performed with the IRAM 30m telescope towards the line of sight to the blazar B0355+508 aiming for the (2-1) transitions of PN, PO, HCP and CP. This line of sight incorporates five diffuse/translucent clouds. The (2-1) transitions of the PN, PO, HCP and CP were not detected. We report detections of the (1-0) lines of $\mathrm{^{13}CO}$, HNC and CN along with a first detection of $\mathrm{C^{34}S}$. We have reproduced the observations of HNC, CN, CS and CO in every cloud by applying typical conditions for diffuse/translucent clouds. According to our best-fit model, the most abundant P-bearing species are HCP and CP ($\sim10^{-10}$), followed by PN, PO and $\mathrm{PH_3}$ ($\sim10^{-11}$). We show that the production of P-bearing species is favoured towards translucent rather than diffuse clouds, where the environment provides a stronger shielding from the interstellar radiation. Based on our improved model, the (1-0) transitions of HCP, CP, PN and PO are expected to be detectable with estimated intensities up to $\sim200$ mK.


title: Gaps between primes

abstract: We discuss recent advances on weak forms of the Prime $k$-tuple Conjecture, and its role in proving new estimates for the existence of small gaps between primes and the existence of large gaps between primes.


title: The KLEVER Survey: Spatially resolved metallicity maps and gradients in a sample of 1.2 < z < 2.5 lensed galaxies

abstract: We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2<z<2.5 galaxies by means of a full coverage of the YJ, H and K near-infrared bands. Detailed metallicity maps and gradients are derived for a sub-sample of 28 galaxies from reconstructed source plane emission line maps, exploiting the variety of different emission line diagnostics provided by the broad wavelength coverage of the survey. About 85% of these galaxies are characterised by metallicity gradients shallower than 0.05 dex/kpc and 89% are consistent with flat slope within 3$\sigma$ (67% within 1$\sigma$), suggesting a mild evolution with cosmic time. In the context of cosmological simulations and chemical evolution models, the presence of efficient feedback mechanisms and/or extended star formation profiles on top of the classical "inside-out" scenario of mass assembly is generally required to reproduce the observed flatness of the metallicity gradients beyond z$\sim$1 . Three galaxies with significantly (> 3$\sigma$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation.


title: Hyperboloidal framework for the Kerr spacetime

abstract: Motivated by the need of a robust geometrical framework for the calculation of long, and highly accurate waveforms for extreme-mass-ratio inspirals, this work presents an extensive study of the hyperboloidal formalism for the Kerr spacetime and the Teukolsky equation. In a first step, we introduce a generic coordinate system foliating the Kerr spacetime into hypersurfaces of constant time extending between the black-hole horizon and future null infinity, while keeping track of the underlying degrees of freedom. Then, we express the Teukolsky equation in terms of these generic coordinates with focus on applications in both the time and frequency domains. Specifically, we derive a wave-like equation in $2+1$ dimensions, whose unique solution follows directly from the prescription of initial data (no external boundary conditions). Moreover, we extend the hyperboloidal formulation into the frequency domain. A comparison with the standard form of the Teukolsky equations allows us to express the regularisation factors in terms of the hyperboloidal degrees of freedom. In the second part, we discuss several hyperboloidal gauges for the Kerr solution. Of particular importance, this paper introduces the minimal gauge. The resulting expressions for the Kerr metric and underlying equations are simple enough for eventual (semi)-analytical studies. Despite the simplicity, the gauge has a very rich structure as it naturally leads to two possible limits to extremality, namely the standard extremal Kerr spacetime and its near-horizon geometry. When applied to the Teukolsky equation in the frequency domain, we show that the minimal gauge actually provides the spacetime counterpart of the well-known Leaver's formalism. Finally, we recast the hyperboloidal gauges for the Kerr spacetime available in the literature within the framework introduced here.


title: Adversarial Machine Learning Phases of Matter

abstract: We study the robustness of machine learning approaches to adversarial perturbations, with a focus on supervised learning scenarios. We find that typical phase classifiers based on deep neural networks are extremely vulnerable to adversarial perturbations: adding a tiny amount of carefully crafted noises into the original legitimate examples will cause the classifiers to make incorrect predictions at a notably high confidence level. Through the lens of activation maps, we find that some important underlying physical principles and symmetries remain to be adequately captured for classifiers with even near-perfect performance. This explains why adversarial perturbations exist for fooling these classifiers. In addition, we find that, after adversarial training, the classifiers will become more consistent with physical laws and consequently more robust to certain kinds of adversarial perturbations. Our results provide valuable guidance for both theoretical and experimental future studies on applying machine learning techniques to condensed matter physics.


title: Multi-Layer Restricted Boltzmann Machine Representation of 1D Quantum Many-Body Wave Functions

abstract: We consider representing two classes of 1D quantum wave functions of spin systems, including the AKLT and CFT correlator wave functions, in terms of multi-layer restricted Boltzmann machines. In our prescription, the AKLT wave function can be exactly represented by a 2-layer restricted Boltzmann machine with five hidden spins per visible spin. The construction can be generalized to prove that any MPS wave function on $N$ unit cells with finite bond dimension can be approximated by a 2-layer restricted Boltzmann machine with $\mathcal{O}(N)$ hidden spins within an error which scales linearly with $N$. The Haldane-Shastry wave function or a chiral boson CFT correlator wave function, as any Jastrow type of wave functions, can be exactly written as a 1-layer Boltzmann machine with $\mathcal{O}(N^2)$ hidden spins and $N$ visible spins. Applying the cumulant expansion, we further find that the chiral boson CFT correlator wave function (with small vertex operator conformal dimension $\alpha$, i.e., $\alpha<0.1$) can be approximated, within 99.9\% accuracy up to 22 visible spins, by a 1-layer RBM with $\mathcal{O}(N)$ hidden spins. The cumulant expansion also leads us to a physically inspiring result in which the hidden spins of the restricted Boltzmann machine can be interpreted as the conformal tower of the chiral boson CFT on the cylinder.


title: Trialities of minimally supersymmetric 2d gauge theories

abstract: We study dynamics of two-dimensional N=(0,1) supersymmetric gauge theories. In particular, we propose that there is an infrared triality between certain triples of theories with orthogonal and symplectic gauge groups. The proposal is supported by matching of anomalies and elliptic genera. This triality can be viewed as a (0,1) counterpart of the (0,2) triality proposed earlier by two of the authors and A. Gadde. We also describe the relation between global anomalies in gauge theoretic and sigma-model descriptions, filling in a gap in the present literature.


title: Combining Sterile Neutrino Fits to Short Baseline Data with IceCube Data

abstract: Recent global fits to short-baseline neutrino oscillation data have been performed finding preference for a sterile neutrino solution (3+1) over null. In the most recent iteration, it was pointed out that an unstable sterile neutrino (3+1+decay) may be a better description of the data. This is due to the fact that this model significantly reduces the tension between appearance and disappearance datasets. In this work, we add a one-year IceCube dataset to the global fit obtaining new results for the standard 3+1 and 3+1+decay sterile neutrino scenarios. We find that the 3+1+decay model provides a better fit than the 3+1, even in the presence of IceCube, with reduced appearance to disappearance tension. The 3+1+decay model is a 5.4$\sigma$ improvement over the null hypothesis and a 2.8$\sigma$ improvement over the standard 3+1 model.


title: On the Gas Content, Star Formation Efficiency, and Environmental Quenching of Massive Galaxies in Proto-Clusters at z~2.0-2.5

abstract: We present ALMA Band 6 (nu=233GHz, lambda=1.3mm) continuum observations towards 68 'normal' star-forming galaxies within two Coma-like progenitor structures at z=2.10 and 2.47, from which ISM masses are derived, providing the largest census of molecular gas mass in overdense environments at these redshifts. Our sample comprises galaxies with a stellar mass range of 1x10^9M_sun - 4x10^11M_sun with a mean M_*~6x10^10M_sun. Combining these measurements with multiwavelength observations and SED modeling, we characterize the gas mass fraction and the star formation efficiency, and infer the impact of the environment on galaxies' evolution. Most of our detected galaxies (~70%) have star formation efficiencies and gas fractions similar to those found for coeval field galaxies and in agreement with the field scaling relations. However, we do find that the proto-clusters contain an increased fraction of massive, gas-poor galaxies, with low gas fractions (f_gas<6-10%) and red rest-frame ultraviolet/optical colors typical of post-starburst and passive galaxies. The relatively high abundance of passive galaxies suggests an accelerated evolution of massive galaxies in proto-cluster environments. The large fraction of quenched galaxies in these overdense structures also implies that environmental quenching takes place during the early phases of cluster assembly, even before virialization. From our data, we derive a quenching efficiency of E_q~0.45 and an upper limit on the quenching timescale of T_q<1Gyr.


title: MUSE Analysis of Gas around Galaxies (MAGG) -- I: Survey design and the environment of a near pristine gas cloud at z~3.5

abstract: We present the design, methods, and first results of the MUSE Analysis of Gas around Galaxies (MAGG) survey, a large programme on the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope (VLT) which targets 28 z > 3.2 quasars to investigate the connection between optically-thick gas and galaxies at z~3-4. MAGG maps the environment of 52 strong absorption line systems at z > 3, providing the first statistical sample of galaxies associated with gas-rich structures in the early Universe. In this paper, we study the galaxy population around a very metal poor gas cloud at z~3.5 towards the quasar J124957.23-015928.8. We detect three Lyman alpha emitters within <200km/s of the cloud redshift, at projected separations <185 kpc (physical). The presence of star-forming galaxies near a very metal-poor cloud indicates that metal enrichment is still spatially inhomogeneous at this redshift. Based on its very low metallicity and the presence of nearby galaxies, we propose that the most likely scenario for this LLS is that it lies within a filament which may be accreting onto a nearby galaxy. Taken together with the small number of other LLSs studied with MUSE, the observations to date show a range of different environments near strong absorption systems. The full MAGG survey will significantly expand this sample and enable a statistical analysis of the link between gas and galaxies to pin down the origin of these diverse environments at z~3-4.


title: Mediator assisted cooling in quantum annealing

abstract: We show a significant reduction of errors for an architecture of quantum annealers (QA) where bosonic modes mediate the interaction between qubits. These systems have a large redundancy in the subspace of solutions, supported by arbitrarily large bosonic occupations. We explain how this redundancy leads to a mitigation of errors when the bosonic modes operate in the ultrastrong coupling regime. Numerical simulations also predict a large increase of qubit coherence for a specific annealing problem with mediated interactions. We provide evidences that noise reduction occurs in more general types of quantum computers with similar architectures.


title: Strong Supercooling as a Consequence of Renormalization Group Consistency

abstract: Classically scale-invariant models are attractive not only because they may offer a solution to the long-standing gauge hierarchy problem, but also due to their role in facilitating strongly supercooled cosmic phase transitions. In this paper, we investigate the interplay between these two aspects. We do so in the context of the electroweak phase transition (EWPT) in the minimal scale-invariant theory. We find that the amount of supercooling generally decreases for increasing scalar couplings. However, the stabilization of the electroweak scale against the Planck scale requires the absence of Landau poles in the respective energy range. Scalar couplings at the TeV scale can therefore not become larger than $\mathcal{O}(10^{-1})$. As a consequence, all fully consistent parameter points predict the EWPT not to complete before the QCD transition, which then eventually triggers the generation of the electroweak scale. We also discuss the potential of the model to give rise to an observable gravitational wave signature, as well as the possibility to accommodate a dark matter candidate.


title: BART: Denoising Sequence-to-Sequence Pre-training for Natural Language Generation, Translation, and Comprehension

abstract: We present BART, a denoising autoencoder for pretraining sequence-to-sequence models. BART is trained by (1) corrupting text with an arbitrary noising function, and (2) learning a model to reconstruct the original text. It uses a standard Tranformer-based neural machine translation architecture which, despite its simplicity, can be seen as generalizing BERT (due to the bidirectional encoder), GPT (with the left-to-right decoder), and many other more recent pretraining schemes. We evaluate a number of noising approaches, finding the best performance by both randomly shuffling the order of the original sentences and using a novel in-filling scheme, where spans of text are replaced with a single mask token. BART is particularly effective when fine tuned for text generation but also works well for comprehension tasks. It matches the performance of RoBERTa with comparable training resources on GLUE and SQuAD, achieves new state-of-the-art results on a range of abstractive dialogue, question answering, and summarization tasks, with gains of up to 6 ROUGE. BART also provides a 1.1 BLEU increase over a back-translation system for machine translation, with only target language pretraining. We also report ablation experiments that replicate other pretraining schemes within the BART framework, to better measure which factors most influence end-task performance.


title: Newly discovered dwarf galaxies in the MATLAS low density fields

abstract: We present the photometric properties of 2210 newly identified dwarf galaxy candidates in the MATLAS fields. The Mass Assembly of early Type gaLAxies with their fine Structures (MATLAS) deep imaging survey mapped $\sim$142 deg$^2$ of the sky around nearby isolated early type galaxies using MegaCam on the Canada-France-Hawaii Telescope, reaching surface brightnesses of $\sim$ 28.5 - 29 in the g-band. The dwarf candidates were identified through a direct visual inspection of the images and by visually cleaning a sample selected using a partially automated approach, and were morphologically classified at the time of identification. Approximately 75% of our candidates are dEs, indicating that a large number of early type dwarfs also populate low density environments, and 23.2% are nucleated. Distances were determined for 13.5% of our sample using pre-existing $z_{spec}$ measurements and HI detections. We confirm the dwarf nature for 99% of this sub-sample based on a magnitude cut $M_g$ = -18. Additionally, most of these ($\sim$90%) have relative velocities suggesting that they form a satellite population around nearby massive galaxies rather than an isolated field sample. Assuming that the candidates over the whole survey are satellites of the nearby galaxies, we demonstrate that the MATLAS dwarfs follow the same scaling relations as dwarfs in the Local Group as well as the Virgo and Fornax clusters. We also find that the nucleated fraction increases with $M_g$, and find evidence of a morphology-density relation for dwarfs around isolated massive galaxies.


title: Towards Scalable Continuous-Time Trajectory Optimization for Multi-Robot Navigation

abstract: Scalable multi-robot transition is essential for ubiquitous adoption of robots. As a step towards it, a computationally efficient decentralized algorithm for continuous-time trajectory optimization in multi-robot scenarios based upon model predictive control is introduced. The robots communicate only their current states and goals rather than sharing their whole trajectory; using this data each robot predicts a continuous-time trajectory for every other robot exploiting optimal control based motion primitives that are corrected for spatial inter-robot interactions using least squares. A non linear program (NLP) is formulated for collision avoidance with the predicted trajectories of other robots. The NLP is condensed by using time as a parametrization resulting in an unconstrained optimization problem and can be solved in a fast and efficient manner. Additionally, the algorithm resizes the robot to accommodate it's trajectory tracking error. The algorithm was tested in simulations on Gazebo with aerial robots. Early results indicate that the proposed algorithm is efficient for upto forty homogeneous robots and twenty one heterogeneous robots occupying 20\% of the available space.


title: On the possibilities of classical nova identifications among historical Far Eastern guest star observations

abstract: More than 100 guest star observations have been obtained by Chinese, Korean, Japanese and Vietnamese astronomers between ~600 BCE and ~1690 CE. Comparing the coordinates from the information given in old texts for eight supernova recoveries with modern supernova remnant positions, we estimate a typical positional accuracy of the order of 0.3 to 7 deg for these supernovae. These values turn out to be also a start for the expected deviation angle between a classical nova observed as a guest star and its modern counterpart among known cataclysmic variables (CVs). However, there are considerable disagreements among modern authors in the interpretation of ancient Far Eastern texts, emphasizing the need to consult again the original historic sources, in order to improve the positioning reliability. We also discuss the typical amplitudes of well observed classical novae and find that modern counterparts of nova guest stars should be V = 18 mag and thus easily observable. In this context we also consider the "hibernation scenario" and conclude that it is impossible to decide from currently-available observations whether hibernation is common. In addition to the limiting magnitude around 2 mag for ancient guest star detections mentioned in the literature, we consider the possibility that also fainter guest stars (4-5 mag) could have been detected by ancient observers and give arguments in favor of this possibility. For these limits we compare the expected nova detection rate of ancient naked-eye observers with that during modern times, and conclude that they coincide in order of magnitude, which implies that, indeed, a considerable number of classical nova remnants should be hidden among the Far Eastern guest star reports. Finally, we present a statistical analysis of the probability of casual misidentifications based on frequency and galactic distribution of CVs in the AAVSO-VSX catalogue.


title: On a Conjecture of Nagy on Extremal Densities

abstract: We disprove a conjecture of Nagy on the maximum number of copies N(G,H) of a fixed graph G in a large graph H with prescribed edge density. Nagy conjectured that for all G, the quantity N(G,H) is asymptotically maximised by either a quasi-star or a quasi-clique. We show this is false for infinitely many graphs, the smallest of which has 6 vertices and 6 edges. We also propose some new conjectures for the behaviour of N(G,H), and present some evidence for them.


title: Ordered Memory

abstract: Stack-augmented recurrent neural networks (RNNs) have been of interest to the deep learning community for some time. However, the difficulty of training memory models remains a problem obstructing the widespread use of such models. In this paper, we propose the Ordered Memory architecture. Inspired by Ordered Neurons (Shen et al., 2018), we introduce a new attention-based mechanism and use its cumulative probability to control the writing and erasing operation of the memory. We also introduce a new Gated Recursive Cell to compose lower-level representations into higher-level representation. We demonstrate that our model achieves strong performance on the logical inference task (Bowman et al., 2015)and the ListOps (Nangia and Bowman, 2018) task. We can also interpret the model to retrieve the induced tree structure, and find that these induced structures align with the ground truth. Finally, we evaluate our model on the Stanford SentimentTreebank tasks (Socher et al., 2013), and find that it performs comparatively with the state-of-the-art methods in the literature.


title: Angular momentum and morphological sequence of massive galaxies through Dark Sage

abstract: We study the present-day connection between galaxy morphology and angular momentum using the {\sc Dark Sage} semi-analytic model of galaxy formation. For galaxies between $ 10^{11}-10^{12} \mathrm{M}_{\odot}$ in stellar mass, the model successfully predicts the observed trend whereby galaxies with more prominent disks exhibit higher {\em stellar} disk specific angular momentum ($j_{\rm stellar, disk}$) at fixed stellar mass. However, when we include the gas in the disk, bulge-dominated galaxies have the highest {\em total} disk specific angular momentum ($j_{\rm total, disk}$). We attribute this to a large contribution from an extended disk of cold gas in typical bulge-dominated galaxies. We find the relationship between $j_{\rm dark matter}$ and morphology to be quite complex. Surprisingly, in this stellar mass range, not only do bulge-dominated galaxies tend to live in halos with higher $j_{\rm dark matter}$ than disk-dominated galaxies, but intermediate galaxies (those with roughly equal fractions of bulge and disk mass) have the lowest $j_{\rm dark matter}$ of all. Yet, when controlling for halo mass, rather than stellar mass, the relationship between $j_{\rm dark matter}$ and morphology vanishes. Based on these results, halo mass rather than angular momentum is the main driver of the predicted morphology sequence at high masses. In fact, in our stellar mass range, disk-dominated galaxies live in dark matter halos that are roughly 1/10th the mass of their bulge-dominated counterparts.


title: On the Count Probability of Many Correlated Symmetric Events

abstract: We consider $N$ events that are defined on a common probability space. Those events shell have a common probability function that is symmetric with respect to interchanging the events. We ask for the probability distribution of the number of events that occur. If the probability of a single event is proportional to $1/N$ the resulting count probability is Poisson distributed in the limit of $N\rightarrow \infty$ for independent events. In this paper we calculate the characteristic function of the limiting count probability distribution for events that are correlated up to an arbitrary but finite order.


title: A hierarchical mean field model of interacting spins

abstract: We consider a system of hierarchical interacting spins under dynamics of spin-flip type with a ferromagnetic mean field interaction, scaling with the hierarchical distance, coupled with a system of linearly interacting hierarchical diffusions of Ornstein-Uhlenbeck type. In particular, the diffusive variables enter in the spin-flip rates, effectively acting as dynamical magnetic fields. In absence of the diffusions, the spin-flip dynamics can be thought of as a modification of the Curie--Weiss model. We study the mean field and the two-level hierarchical model, in the latter case restricting to a subcritical regime, corresponding to high temperatures, obtaining macroscopic limits at different spatio-temporal scales and studying the phase transitions in the system. We also formulate a generalization of our results to the $k$-th level hierarchical case, for any $k$ finite, in the subcritical regime. We finally address the supercritical regime, in the zero-temperature limit, for the two-level hierarchical case, proceeding heuristically with the support of numerics.


title: A z=0 Multi-wavelength Galaxy Synthesis I: A WISE and GALEX Atlas of Local Galaxies

abstract: We present an atlas of ultraviolet and infrared images of ~15,750 local (d < 50 Mpc) galaxies, as observed by NASA's WISE and GALEX missions. These maps have matched resolution (FWHM 7.5'' and 15''), matched astrometry, and a common procedure for background removal. We demonstrate that they agree well with resolved intensity measurements and integrated photometry from previous surveys. This atlas represents the first part of a program (the z=0 Multi-wavelength Galaxy Synthesis) to create a large, uniform database of resolved measurements of gas and dust in nearby galaxies. The images and associated catalogs are publicly available at the NASA/IPAC Infrared Science Archive. This atlas allows us estimate local and integrated star formation rates (SFRs) and stellar masses (M$_\star$) across the local galaxy population in a uniform way. In the appendix, we use the population synthesis fits of Salim et al. (2016, 2018) to calibrate integrated M$_\star$ and SFR estimators based on GALEX and WISE. Because they leverage an SDSS-base training set of >100,000 galaxies, these calibrations have high precision and allow us to rigorously compare local galaxies to Sloan Digital Sky Survey results. We provide these SFR and M$_\star$ estimates for all galaxies in our sample and show that our results yield a "main sequence" of star forming galaxies comparable to previous work. We also show the distribution of intensities from resolved galaxies in NUV-to-WISE1 vs. WISE1-to-WISE3 space, which captures much of the key physics accessed by these bands.


title: Stochastic entropy production in diffusive systems

abstract: Computing the stochastic entropy production associated with the evolution of a stochastic dynamical system is a well-established problem. In a small number of cases such as the Ornstein-Uhlenbeck process, of which we give a complete exposition, the distribution of entropy production can be obtained analytically, but in general it is much harder. A recent development in solving the Fokker-Planck equation, in which the solution is written as a product of positive functions, enables the distribution to be obtained approximately, with the assistance of simple numerical techniques. Using examples in one and higher dimension, we demonstrate how such a framework is very convenient for the computation of stochastic entropy production in diffusion processes.


title: Locally recoverable codes on surfaces

abstract: A linear error correcting code is a subspace of a finite-dimensional space over a finite field with a fixed coordinate system. Such a code is said to be locally recoverable with locality $r$ if, for every coordinate, its value at a codeword can be deduced from the value of (certain) $r$ other coordinates of the codeword. These codes have found many recent applications, e.g., to distributed cloud storage. We will discuss the problem of constructing good locally recoverable codes and present some constructions using algebraic surfaces that improve previous constructions and sometimes provide codes that are optimal in a precise sense. The main conceptual contribution of this paper is to consider surfaces fibered over a curve in such a way that each recovery set is constructed from points in a single fiber. This allows us to use the geometry of the fiber to guarantee the local recoverability and use the global geometry of the surface to get a hold on the standard parameters of our codes. We look in detail at situations where the fibers are rational or elliptic curves and provide many examples applying our methods.


title: Contrast data mining for the MSSM from strings

abstract: We apply techniques from data mining to the heterotic orbifold landscape in order to identify new MSSM-like string models. To do so, so-called contrast patterns are uncovered that help to distinguish between areas in the landscape that contain MSSM-like models and the rest of the landscape. First, we develop these patterns in the well-known $\mathbb{Z}_6$-II orbifold geometry and then we generalize them to all other $\mathbb{Z}_N$ orbifold geometries. Our contrast patterns have a clear physical interpretation and are easy to check for a given string model. Hence, they can be used to scale down the potentially interesting area in the landscape, which significantly enhances the search for MSSM-like models. Thus, by deploying the knowledge gain from contrast mining into a new search algorithm we create many novel MSSM-like models, especially in corners of the landscape that were hardly accessible by the conventional search algorithm, for example, MSSM-like $\mathbb{Z}_6$-II models with $\Delta(54)$ flavor symmetry.


title: Vector Leptoquarks Beyond Tree Level

abstract: Models with massive vector leptoquarks, resulting from an $SU(4)$ gauge symmetry spontaneously broken at the TeV scale, are of great phenomenological interest given the current "anomalies" in semileptonic $B$ decays. We analyze the relations between low- and high-energy observables in such class of models to next-to-leading order accuracy in the $SU(4)$ gauge coupling $g_4$. For large values of $g_4$, motivated by recent $B$-physics data, one-loop corrections are sizeable. The main effect is an enhanced contribution at low-energy, at fixed on-shell couplings. This result has important implications for current and future high-energy searches of vector leptoquark models.


title: Enskog kinetic theory for multicomponent granular suspensions

abstract: The Navier--Stokes transport coefficients of multicomponent granular suspensions at moderate densities are obtained in the context of the (inelastic) Enskog kinetic theory. The suspension is modeled as an ensemble of solid particles where the influence of the interstitial gas on grains is via a viscous drag force plus a stochastic Langevin-like term defined in terms of a background temperature. In the absence of spatial gradients, it is shown first that the system reaches a homogeneous steady state where the energy lost by inelastic collisions and viscous friction is compensated for by the energy injected by the stochastic force. Once the homogeneous steady state is characterized, a \emph{normal} solution to the set of Enskog equations is obtained by means of the Chapman--Enskog expansion around the \emph{local} version of the homogeneous state. To first-order in spatial gradients, the Chapman--Enskog solution allows us to identify the Navier--Stokes transport coefficients associated with the mass, momentum, and heat fluxes. In addition, the first-order contributions to the partial temperatures and the cooling rate are also calculated. Explicit forms for the diffusion coefficients, the shear and bulk viscosities, and the first-order contributions to the partial temperatures and the cooling rate are obtained in steady-state conditions by retaining the leading terms in a Sonine polynomial expansion. The results show that the dependence of the transport coefficients on inelasticity is clearly different from that found in its granular counterpart (no gas phase). The present work extends previous theoretical results for \emph{dilute} multicomponent granular suspensions [Khalil and Garz\'o, Phys. Rev. E \textbf{88}, 052201 (2013)] to higher densities.


title: Improved spectral convergence rates for graph Laplacians on epsilon-graphs and k-NN graphs

abstract: In this paper we improve the spectral convergence rates for graph-based approximations of Laplace-Beltrami operators constructed from random data. We utilize regularity of the continuum eigenfunctions and strong pointwise consistency results to prove that spectral convergence rates are the same as the pointwise consistency rates for graph Laplacians. In particular, for an optimal choice of the graph connectivity $\varepsilon$, our results show that the eigenvalues and eigenvectors of the graph Laplacian converge to those of the Laplace-Beltrami operator at a rate of $O(n^{-1/(m+4)})$, up to log factors, where $m$ is the manifold dimension and $n$ is the number of vertices in the graph. Our approach is general and allows us to analyze a large variety of graph constructions that include $\varepsilon$-graphs and $k$-NN graphs.


title: Proper $r$-harmonic functions on the Thurston geometries

abstract: For any positive natural number $r\in\N^+$ we construct new explicit proper $r$-harmonic functions on the celebrated $3$-dimensional Thurston geometries $\Sol$, $\Nil$, $\SL2$, $\H^2\times\rn$ and $\s^2\times\rn$.


title: Adiabatic Quantum Estimation: A Numerical Study of the Heisenberg XX Model with Antisymmetric Exchange

abstract: In this paper, we address the adiabatic technique for quantum estimation of the azimuthal orientation of a magnetic field. Exactly solving a model consisting of a two-qubit system, where one of which is driven by a static magnetic field while the other is coupled with the magnetic field rotating adiabatically, we obtain the analytical expression of the quantum Fisher information (QFI). We investigate how the two-qubit system can be used to probe the azimuthal direction of the field and analyze the roles of the intensities of the magnetic fields, Dzyaloshinskii-Moriya interaction, spin-spin coupling coefficient, and the polar orientation of the rotating field on the precision of the estimation. In particular, it is illustrated that the QFI trapping or saturation may occur if the qubit is subjected to a strong rotating field. Moreover, we discuss how the azimuthal direction of the rotating field can be estimated using only the qubit not affected by that field and investigate the conditions under which this strategy is more efficient than use of the qubit locally interacting with the adiabatically rotating field. Interestingly, in the one-qubit scenario, it was found that when the rotating field is weak, the best estimation is achieved by subjecting the probe to a static magnetic field.


title: Practical Repetition-Aware Grammar Compression

abstract: The goal of grammar compression is to construct a small sized context free grammar which uniquely generates the input text data. Among grammar compression methods, RePair is known for its good practical compression performance. MR-RePair was recently proposed as an improvement to RePair for constructing small-sized context free grammar for repetitive text data. However, a compact encoding scheme has not been discussed for MR-RePair. We propose a practical encoding method for MR-RePair and show its effectiveness through comparative experiments. Moreover, we extend MR-RePair to run-length context free grammar and design a novel variant for it called RL-MR-RePair. We experimentally demonstrate that a compression scheme consisting of RL-MR-RePair and the proposed encoding method show good performance on real repetitive datasets.


title: Planetary Astronomy-Understanding the Origin of the Solar System

abstract: There is a vibrant and effective planetary science community in Canada. We do research in the areas of meteoritics, asteroid and trans-Neptunian object orbits and compositions, and space weather, and are involved in space probe missions to study planetary surfaces and interiors. For Canadian planetary scientists to deliver the highest scientific impact possible, we have several recommendations. Our top recommendation is to join LSST and gain access to the full data releases by hosting a data centre, which could be done by adding to the CADC, which is already highly involved in hosting planetary data and supporting computational modelling for orbital studies. We also support MSE, which can provide spectroscopy and thus compositional information for thousands of small bodies. We support a Canadian-led microsatellite, POEP, which will provide small body sizes by measuring occultations. We support the idea of piggybacking space weather instruments on other astronomical space probes to provide data for the space weather community. Many Canadian planetary scientists are involved in space probe missions, but through haphazard and temporary arrangements like co-appointments at US institutions, so we would like the community to support Canadian researchers to participate in these large, international missions.


title: Classical nucleation theory predicts the shape of the nucleus in homogeneous solidification

abstract: Macroscopic models of nucleation provide powerful tools for understanding activated phase transition processes. These models do not provide atomistic insights and can thus sometime lack material-specific descriptions. Here we provide a comprehensive framework for constructing a continuum picture from an atomistic simulation of homogeneous nucleation. We use this framework to determine the shape of the equilibrium solid nucleus that forms inside bulk liquid for a Lennard-Jones potential. From this shape, we then extract the anisotropy of the solid-liquid interfacial free energy, by performing a reverse Wulff construction in the space of spherical harmonic expansions. We find that the shape of the nucleus is nearly spherical and that its anisotropy can be perfectly described using classical models.


title: BSM Physics from Enlarged Gauge Symmetry: the 331 Model, a case of study

abstract: We discuss the most relevant features of a BSM model with extra gauge symmetry, the so called 331 model. The gauge group is $SU(3)_c\times SU(3)_L \times U(1)_X$ impling the presence of extra gauge bosons, both charged and neutral, as well as extra/exotic fermions and an enlarged scalar sector. We present the relevant phenomenology of doubly-charged gauge bosons, which are a distinctive feature of a version of the 331 model, and discuss the role of BSM phenomenology as a tool for testing GUT (inspired) theories.


title: The Quantum Alternating Operator Ansatz on Maximum k-Vertex Cover

abstract: The Quantum Alternating Operator Ansatz is a generalization of the Quantum Approximate Optimization Algorithm (QAOA) designed for finding approximate solutions to combinatorial optimization problems with hard constraints. In this paper, we study Maximum $k$-Vertex Cover under this ansatz due to its modest complexity, while still being more complex than the well studied problems of Max-Cut and Max E3-LIN2. Our approach includes (i) a performance comparison between easy-to-prepare classical states and Dicke states as starting states, (ii) a performance comparison between two $XY$-Hamiltonian mixing operators: the ring mixer and the complete graph mixer, (iii) an analysis of the distribution of solutions via Monte Carlo sampling, and (iv) the exploration of efficient angle selection strategies. Our results are: (i) Dicke states improve performance compared to easy-to-prepare classical states, (ii) an upper bound on the simulation of the complete graph mixer, (iii) the complete graph mixer improves performance relative to the ring mixer, (iv) numerical results indicating the standard deviation of the distribution of solutions decreases exponentially in $p$ (the number of rounds in the algorithm), requiring an exponential number of random samples to find a better solution in the next round, and (v) a correlation of angle parameters which exhibit high quality solutions that behave similarly to a discretized version of the Quantum Adiabatic Algorithm.


title: A Novel Procedure for the Assessment of the Seismic Performance of Frame Structures by Means of Limit Analysis

abstract: Limit analysis is a computationally efficient tool to assess the resistance and the failure mode of structures but does not provide any information on the displacement capacity, which is one of the concepts which most affects the seismic safety. Therefore, since many researchers did not consider limit analysis as a possible tool for the seismic assessment of structures, its widespread employment has been prevented. In this paper this common belief is questioned and the authors show that limit analysis can be useful in the evaluation of the seismic performance of frame structures. In particular, to overcome the limitation on the possibility to evaluate the displacements of a structure based on a limit analysis approach, an approximated capacity curve is reconstructed. The latter is based on a limit analysis strategy, which takes into account the second order effects, and evaluates the displacement capacity considering a post-peak softening branch and a threshold on the allowed plastic rotations. Then, based on this simplified capacity curve, an equivalent single degree of freedom system is defined in order to assess the seismic performance of frame structures. The proposed simplified strategy is implemented in a dedicated software and the obtained results are validated with well-established approaches based on nonlinear static analyses, showing the reliability and the computational efficiency of this methodology


title: Four-Field Mixed Finite Element Methods for Incompressible Nonlinear Elasticity

abstract: We introduce conformal mixed finite element methods for $2$D and $3$D incompressible nonlinear elasticity in terms of displacement, displacement gradient, the first Piola-Kirchhoff stress tensor, and pressure, where finite elements for the $\mathrm{curl}$ and the $\mathrm{div}$ operators are used to discretize strain and stress, respectively. These choices of elements follow from the strain compatibility and the momentum balance law. Some inf-sup conditions are derived to study the stability of methods. By considering $96$ choices of simplicial finite elements of degree less than or equal to $2$ in $2$D and $3$D, we conclude that $28$ choices in $2$D and $6$ choices in $3$D satisfy these inf-sup conditions. The performance of stable finite element choices are numerically studied. Although the proposed methods are computationally more expensive than the standard two-field methods for incompressible elasticity, they are potentially useful for accurate approximations of strain and stress as they are independently computed in the solution process.


title: Parametric Interpolation Framework for Scalar Conservation Laws

abstract: In this paper we present a novel framework for obtaining high-order numerical methods for scalar conservation laws in one-space dimension for both the homogeneous and non-homogeneous case. The numerical schemes for these two settings are somewhat different in the presence of shocks, however at their core they both rely heavily on the solution curve being represented parametrically. By utilizing high-order parametric interpolation techniques we succeed to obtain fifth order accuracy ( in space ) everywhere in the computation domain, including the shock location itself. In the presence of source terms a slight modification is required, yet the spatial order is maintained but with an additional temporal error appearing. We provide a detailed discussion of a sample scheme for non-homogeneous problems which obtains fifth order in space and fourth order in time even in the presence of shocks.


title: Explicit Diagonalization of Pair Interaction Models

abstract: We provide a general method for constructing bosonic Bogoliubov transformations that diagonalize a general class of quadratic Hamiltonians. These Hamiltonians describe the pair interaction models. Bogoliubov transformations are constructed algebraically, and the resulting Hamiltonians become the second quantizations of explicit one-particle Hamiltonians. Moreover, an explicit formula for the ground state energies is given. Our method systematically diagonalizes various models of quantum field theory, including a model of a harmonic oscillator coupled to a Bose field and the Pauli-Fierz models in the dipole approximation.


title: Does Speech enhancement of publicly available data help build robust Speech Recognition Systems?

abstract: Automatic speech recognition (ASR) systems play a key role in many commercial products including voice assistants. Typically, they require large amounts of clean speech data for training which gives an undue advantage to large organizations which have tons of private data. In this paper, we have first curated a fairly big dataset using publicly available data sources. Thereafter, we tried to investigate if we can use publicly available noisy data to train robust ASR systems. We have used speech enhancement to clean the noisy data first and then used it together with its cleaned version to train ASR systems. We have found that using speech enhancement gives 9.5\% better word error rate than training on just noisy data and 9\% better than training on just clean data. It's performance is also comparable to the ideal case scenario when trained on noisy and its clean version.


title: Stability estimates in a partial data inverse boundary value problem for biharmonic operators at high frequencies

abstract: We study the inverse boundary value problems of determining a potential in the Helmholtz type equation for the perturbed biharmonic operator from the knowledge of the partial Cauchy data set. Our geometric setting is that of a domain whose inaccessible portion of the boundary is contained in a hyperplane, and we are given the Cauchy data set on the complement. The uniqueness and logarithmic stability for this problem were established in [37] and [7], respectively. We establish stability estimates in the high frequency regime, with an explicit dependence on the frequency parameter, under mild regularity assumptions on the potentials, sharpening those of [7].


title: Centrosymmetric Stochastic Matrices

abstract: We consider the convex set $\Gamma_{m,n}$ of $m\times n$ stochastic matrices and the convex set $\Gamma_{m,n}^\pi\subset \Gamma_{m,n}$ of $m\times n$ centrosymmetric stochastic matrices (stochastic matrices that are symmetric under rotation by 180 degrees). For $\Gamma_{m,n}$, we demonstrate a Birkhoff theorem for its extreme points and create a basis from certain $(0,1)$-matrices. For $\Gamma_{m,n}^\pi$, we characterize its extreme points and create bases, whose construction depends on the parity of $m$, using our basis construction for stochastic matrices. For each of $\Gamma_{m,n}$ and $\Gamma_{m,n}^\pi$, we further characterize their extreme points in terms of their associated bipartite graphs, we discuss a graph parameter called the fill and compute it for the various basis elements, and we examine the number of vertices of the faces of these sets. We provide examples illustrating the results throughout.


title: Inequalities for indices of coincidence and entropies

abstract: We consider a probability distribution depending on a real parameter $x$. As functions of $x$, the R\'enyi entropy and the Tsallis entropy can be expressed in terms of the associated index of coincidence $S(x)$. We establish recurrence relations and inequalities for $S(x),$ which can be used in order to get information concerning the two entropies.


title: The moduli space of multi-scale differentials

abstract: We construct a compactification of the moduli spaces of abelian differentials on Riemann surfaces with prescribed zeroes and poles. This compactification, called the moduli space of multi-scale differentials, is a complex orbifold with normal crossing boundary. Locally, our compactification can be described as the normalization of an explicit blowup of the incidence variety compactification, which was defined in [BCGGM18] as the closure of the stratum of abelian differentials in the closure of the Hodge bundle. We also define families of projectivized multi-scale differentials, which gives a proper Deligne-Mumford stack, and our compactification is the orbifold corresponding to it. Moreover, we perform a real oriented blowup of the unprojectivized moduli space of multi-scale differentials such that the $\mathrm{GL}_2(\mathbb R)$-action in the interior of the moduli space extends continuously to the boundary.


title: Efficient Identification in Linear Structural Causal Models with Instrumental Cutsets

abstract: One of the most common mistakes made when performing data analysis is attributing causal meaning to regression coefficients. Formally, a causal effect can only be computed if it is identifiable from a combination of observational data and structural knowledge about the domain under investigation (Pearl, 2000, Ch. 5). Building on the literature of instrumental variables (IVs), a plethora of methods has been developed to identify causal effects in linear systems. Almost invariably, however, the most powerful such methods rely on exponential-time procedures. In this paper, we investigate graphical conditions to allow efficient identification in arbitrary linear structural causal models (SCMs). In particular, we develop a method to efficiently find unconditioned instrumental subsets, which are generalizations of IVs that can be used to tame the complexity of many canonical algorithms found in the literature. Further, we prove that determining whether an effect can be identified with TSID (Weihs et al., 2017), a method more powerful than unconditioned instrumental sets and other efficient identification algorithms, is NP-Complete. Finally, building on the idea of flow constraints, we introduce a new and efficient criterion called Instrumental Cutsets (IC), which is able to solve for parameters missed by all other existing polynomial-time algorithms.


title: X-ray diffractive imaging of controlled gas-phase molecules: Toward imaging of dynamics in the molecular frame

abstract: We report experimental results on the diffractive imaging of three-dimensionally aligned 2,5-diiodothiophene molecules. The molecules were aligned by chirped near-infrared laser pulses, and their structure was probed at a photon energy of 9.5 keV ($\lambda\approx130 \text{pm}$) provided by the Linac Coherent Light Source. Diffracted photons were recorded on the CSPAD detector and a two-dimensional diffraction pattern of the equilibrium structure of 2,5-diiodothiophene was recorded. The retrieved distance between the two iodine atoms agrees with the quantum-chemically calculated molecular structure to within 5 %. The experimental approach allows for the imaging of intrinsic molecular dynamics in the molecular frame, albeit this requires more experimental data which should be readily available at upcoming high-repetition-rate facilities.


title: Relative C*-simplicity and characterizations for normal subgroups

abstract: The notion of a plump subgroup was recently introduced by Amrutam. This is a relativized version of Powers' averaging property, and it is known that Powers' averaging property is equivalent to C*-simplicity. With this in mind, we introduce a relativized notion of C*-simplicity, and show that for normal subgroups it is equivalent to plumpness, along with several other characterizations.


title: Asymptotically unbiased estimation of physical observables with neural samplers

abstract: We propose a general framework for the estimation of observables with generative neural samplers focusing on modern deep generative neural networks that provide an exact sampling probability. In this framework, we present asymptotically unbiased estimators for generic observables, including those that explicitly depend on the partition function such as free energy or entropy, and derive corresponding variance estimators. We demonstrate their practical applicability by numerical experiments for the 2d Ising model which highlight the superiority over existing methods. Our approach greatly enhances the applicability of generative neural samplers to real-world physical systems.


title: Quantifying the Semantic Core of Gender Systems

abstract: Many of the world's languages employ grammatical gender on the lexeme. For example, in Spanish, the word for 'house' (casa) is feminine, whereas the word for 'paper' (papel) is masculine. To a speaker of a genderless language, this assignment seems to exist with neither rhyme nor reason. But is the assignment of inanimate nouns to grammatical genders truly arbitrary? We present the first large-scale investigation of the arbitrariness of noun-gender assignments. To that end, we use canonical correlation analysis to correlate the grammatical gender of inanimate nouns with an externally grounded definition of their lexical semantics. We find that 18 languages exhibit a significant correlation between grammatical gender and lexical semantics.


title: Quasi normal modes for Heavy Vector Mesons in a Finite Density Plasma

abstract: There are strong indications that ultra-relativistic heavy ion collisions, produced in accelarators, lead to the formation of a new state of matter: the quark gluon plasma (QGP). This deconfined QCD matter is expected to exist just for very short times after the collision. All the information one can get about the plasma is obtained from the particles that reach the detectors. Among them, heavy vector mesons are particularly important. The abundance of $ c {\bar c}$ and $b {\bar b}$ states produced in a heavy ion collision is a source of information about the plasma. In contrast to the light mesons, that completely dissociate when the plasma is formed, heavy mesons presumably undergo a partial thermal dissociation. The dissociation degree depends on the temperature and also on the presence of magnetic fields and on the density (chemical potential). So, in order to get information about the plasma out of the quarkonium abundance data, one needs to resort to models that provide the dependence of the dissociation degree on these factors. Holographic phenomenological models provide a nice description for charmonium and bottomonium quasi-states in a plasma. In particular, quasi-normal modes associated with quarkonia states have been studied recently for a plasma with magnetic fields. Here we extend this analysis of quasinormal models to the case when charmonium and bottomonium are inside a plasma with finite density. The complex frequencies obtained are then compared with a Breit Wigner approximation for the peaks of the corresponding thermal spectral functions, in order to investigate the quantitative agreement of the different descriptions of quarkonium quasi-states.


title: Bending of light in the Universe filled with quintessential dark energy

abstract: As a local effect of dynamical dark energy, bending of light in the presence of a spherically symmetric and static black hole surrounded by quintessence has been studied. Having in mind recent observational data, we have treated the problem as a deviation from Kottler space-time. This deviation is measured by a perturbation parameter $\varepsilon$ included in the equation of state parameter of quintessence as $\omega_q=-1+\frac{1}{3}\varepsilon$. Here, the deflection angle is calculated and then the result is compared with \cite{Arakida:2011th} in the limit $\varepsilon\rightarrow 0$ where the quintessence behaves like the cosmological constant. It is shown that unlike the cosmological constant, the effect of quintessence on the photon energy equation can not be absorbed into the definition of impact parameter. Moreover in this paper, we generalize the Kiselev black hole to the case that there is a modified Chaplygin gas as the dark energy component of the universe and show that the resulted metric can be reduced to the Kiselev metric by adjusting some arbitrary parameters.


title: Equivalence groupoid and group classification of a class of variable-coefficient Burgers equations

abstract: We study admissible transformations and Lie symmetries for a class of variable-coefficient Burgers equations. We combine the advanced methods of splitting into normalized subclasses and of mappings between classes that are generated by families of point transformations parameterized by arbitrary elements of the original classes. A nontrivial differential constraint on the arbitrary elements of the class of variable-coefficient Burgers equations leads to its partition into two subclasses, which are related to normalized classes via families of point transformations parameterized by subclasses' arbitrary elements. One of the mapped classes is proved to be normalized in the extended generalized sense, and its effective extended generalized equivalence group is found. Using the mappings between classes and the algebraic method of group classification, we carry out the group classification of the initial class with respect to its equivalence groupoid.


title: Assessment of climate change effects on mountain ecosystems through a cross-site analysis in the Alps and Apennines

abstract: Mountain ecosystems are sensitive indicators of climate change. Long-term studies may be extremely useful in assessing the responses of high-elevation ecosystems to climate change and other anthropogenic drivers. Mountain research sites within the LTER (Long-Term Ecosystem Research) network are representative of various types of ecosystems and span a wide bioclimatic and elevational range. Here, we present a synthesis and a review of the main results from long-term ecological studies in mountain ecosystems at 20 LTER sites in Italy, Switzerland and Austria. We analyzed a set of key climate parameters, such as temperature and snow cover duration, in relation to vascular species composition, plant traits, abundance patterns, pedoclimate, nutrient dynamics in soils and water, phenology and composition of freshwater biota. The overall results highlight the rapid response of mountain ecosystems to climate change. As temperatures increased, vegetation cover in alpine and subalpine summits increased as well. Years with limited snow cover duration caused an increase in soil temperature and microbial biomass during the growing season. Effects on freshwater ecosystems were observed, in terms of increases in solutes, decreases in nitrates and changes in plankton phenology and benthos communities. This work highlights the importance of comparing and integrating long-term ecological data collected in different ecosystems, for a more comprehensive overview of the ecological effects of climate change. Nevertheless, there is a need for i) adopting co-located monitoring site networks to improve our ability to obtain sound results from cross-site analysis, ii) carrying out further studies, with fine spatial and temporal resolutions to improve understanding of responses to extreme events, and iii) increasing comparability and standardizing protocols across networks to clarify local from global patterns.


title: The Tradeoff Between Coverage and Computation in Wireless Networks

abstract: We consider a distributed edge computing scenario consisting of several wireless nodes that are located over an area of interest. Specifically, some of the "master" nodes are tasked to sense the environment (e.g., by acquiring images or videos via cameras) and process the corresponding sensory data, while the other nodes are assigned as "workers" to help the computationally-intensive processing tasks of the masters. A new tradeoff that has not been previously explored in the existing literature arises in such a formulation: On one hand, one wishes to allocate as many master nodes as possible to cover a large area for accurate monitoring. On the other hand, one also wishes to allocate as many worker nodes as possible to maximize the computation rate of the sensed data. It is in the context of this tradeoff that this work is presented. By utilizing the basic physical layer principles of wireless communication systems, we formulate and analyze the tradeoff between the coverage and computation performance of spatial networks. We also present an algorithm to find the optimal tradeoff and demonstrate its performance through numerical simulations.


title: Weight of Evidence as a Basis for Human-Oriented Explanations

abstract: Interpretability is an elusive but highly sought-after characteristic of modern machine learning methods. Recent work has focused on interpretability via $\textit{explanations}$, which justify individual model predictions. In this work, we take a step towards reconciling machine explanations with those that humans produce and prefer by taking inspiration from the study of explanation in philosophy, cognitive science, and the social sciences. We identify key aspects in which these human explanations differ from current machine explanations, distill them into a list of desiderata, and formalize them into a framework via the notion of $\textit{weight of evidence}$ from information theory. Finally, we instantiate this framework in two simple applications and show it produces intuitive and comprehensible explanations.


title: An NIP structure which does not interpret an infinite group but whose Shelah expansion interprets an infinite field

abstract: We describe one.


title: Correlations Between Conduction Electrons in Dense Plasmas

abstract: Most treatments of electron-electron correlations in dense plasmas either ignore them entirely (random phase approximation) or neglect the role of ions (jellium approximation). In this work, we go beyond both these approximations to derive a new formula for the electron-electron static structure factor which properly accounts for the contributions of both ionic structure and quantum-mechanical dynamic response in the electrons. The result can be viewed as a natural extension of the quantum Ornstein-Zernike theory of ionic and electronic correlations, and it is suitable for dense plasmas in which the ions are classical and the conduction electrons are quantum-mechanical. The corresponding electron-electron pair distribution functions are compared with the results of path integral Monte Carlo simulations, showing good agreement whenever no strong electron resonance states are present. We construct approximate potentials of mean force which describe the effective screened interaction between electrons. Significant deviations from Debye-H\"uckel screening are present at temperatures and densities relevant to high energy density experiments involving warm and hot dense plasmas. The presence of correlations between conduction electrons is likely to influence the electron-electron contribution to the electron and thermal conductivity. It is expected that excitation processes involving the conduction electrons (e.g., free-free absorption) will also be affected.


title: An outer shade of Pal: Abundance analysis of the outer halo globular cluster Palomar 13

abstract: At a Galactocentric distance of 27 kpc, Pal 13 is an old globular cluster (GC) in the outer halo. We present a chemical abundance analysis of this remote system from high-resolution spectra obtained with Keck/HIRES. Owing to the low signal-to-noise ratio of the data, our analysis is based on a coaddition of the spectra of 18 member stars. We are able to determine integrated abundance ratios for 16 species of 14 elements, of $\alpha$-elements (Mg,Si,Ca,Ti), Fe-peak (Sc,Mn,Cr,Ni,Cu,and Zn), and n-capture elements (Y,Ba). While the mean Na abundance is found to be slightly enhanced and halo-like, our method does not allow us to probe an abundance spread that would be expected in this light element if multiple populations are present in Pal 13. We find a metal-poor mean metallicity of $-1.91\pm0.05$ (statistical) $\pm$ 0.22 (systematic), confirming that Pal 13 is a typical metal-poor representative of the outer halo. While there are some differences between individual $\alpha$-elements, such as halo-like Mg and Si versus the mildly lower Ca and Ti abundances, the mean [$\alpha$/Fe] of 0.34$\pm$0.06 is consistent with the marginally lower $\alpha$ component of the halo field and GC stars. We discuss our results in the context of other objects in the outer halo and consider which of these objects were likely accreted. We also discuss the properties of their progenitors. While chemically, Pal 13 is similar to Gaia-Enceladus and some of its GCs, this is not supported by its kinematic properties. Moreover, its chemodynamical similarity with NGC 5466, a progeny of the Sequoia accretion event, might indicate a common origin in this progenitor. However, the ambiguities in the full abundance space of this comparison emphasize the difficulties in unequivocally labeling a single GC as an accreted object, let alone assigning it to a single progenitor. (Abridged)


title: A 100-kiloparsec wind feeding the circumgalactic medium of a massive compact galaxy

abstract: Ninety per cent of baryons are located outside galaxies, either in the circumgalactic or intergalactic medium. Theory points to galactic winds as the primary source of the enriched and massive circumgalactic medium. Winds from compact starbursts have been observed to flow to distances somewhat greater than ten kiloparsecs, but the circumgalactic medium typically extends beyond a hundred kiloparsecs. Here we report optical integral field observations of the massive but compact galaxy SDSS J211824.06+001729.4. The oxygen [O II] lines at wavelengths of 3726 and 3729 angstroms reveal an ionized outflow spanning 80 by 100 square kiloparsecs, depositing metal-enriched gas at 10,000 kelvin through an hourglass-shaped nebula that resembles an evacuated and limb-brightened bipolar bubble. We also observe neutral gas phases at temperatures of less than 10,000 kelvin reaching distances of 20 kiloparsecs and velocities of around 1,500 kilometres per second. This multi-phase outflow is probably driven by bursts of star formation, consistent with theory.