title: Semi-Analytic Galaxies - I. Synthesis of environmental and star-forming regulation mechanisms

abstract: We present results from the semi-analytic model of galaxy formation SAG applied on the MultiDark simulation MDPL2. SAG features an updated supernova (SN) feedback scheme and a robust modelling of the environmental effects on satellite galaxies. This incorporates a gradual starvation of the hot gas halo driven by the action of ram pressure stripping (RPS), that can affect the cold gas disc, and tidal stripping (TS), which can act on all baryonic components. Galaxy orbits of orphan satellites are integrated providing adequate positions and velocities for the estimation of RPS and TS. The star formation history and stellar mass assembly of galaxies are sensitive to the redshift dependence implemented in the SN feedback model. We discuss a variant of our model that allows to reconcile the predicted star formation rate density at $z \gtrsim 3$ with the observed one, at the expense of an excess in the faint end of the stellar mass function at $z=2$. The fractions of passive galaxies as a function of stellar mass, halo mass and the halo-centric distances are consistent with observational measurements. The model also reproduces the evolution of the main sequence of star forming central and satellite galaxies. The similarity between them is a result of the gradual starvation of the hot gas halo suffered by satellites, in which RPS plays a dominant role. RPS of the cold gas does not affect the fraction of quenched satellites but it contributes to reach the right atomic hydrogen gas content for more massive satellites ($M_{\star}\gtrsim 10^{10}\,{\rm M}_{\odot}$).


title: Semi-Analytic Galaxies - II. Revealing the role of environmental and mass quenching in galaxy formation

abstract: We use the semi-analytic model of galaxy formation SAG to study the relevance of mass and environmental quenching on satellite galaxies. We find that environmental processes dominate the star formation (SF) quenching of low-mass satellites ($M_{\star} \lesssim 10^{10.5}\, {\rm M}_{\odot}$), whereas high-mass galaxies typically quench as centrals. High-mass galaxies that remain actively forming stars while being accreted are found to be mainly affected by mass quenching after their first infall. For a given stellar mass, our model predicts SF quenching to be less efficient in low-mass haloes both before and after infall, in contradiction with common interpretations of observational data. Our model supports a two-stage scenario to explain the SF quenching. Initially, the SF of satellites resembles that of centrals until the gas cooling rate is reduced to approximately half its value at infall. Then, the SF fades through secular processes that exhaust the cold gas reservoir. This reservoir is not replenished efficiently due to the action of either ram-pressure stripping (RPS) of the hot gas in low-mass satellites, or feedback from the active galactic nucleus (AGN) in high-mass satellites. The delay times for the onset of SF quenching are found to range from $\approx 3\,{\rm Gyr}$ to $\approx 1\,{\rm Gyr}$ for low-mass ($M_{\star} \approx 10^{10}\, {\rm M}_{\odot}$) and high-mass ($M_{\star} \approx 10^{11}\, {\rm M}_{\odot}$) satellites, respectively. SF fades in $\approx 1.5\,{\rm Gyr}$, largely independent of stellar mass. We find that the SF quenching of low-mass satellites supports the so-called delay-then-rapid quenching scenario. However, the SF history of $z=0$ passive satellites of any stellar mass is better described by a delay-then-fade quenching scenario.


title: Zeroth-order general Randi\`c index of $k$-generalized quasi trees

abstract: For a simple graph $G(V,E)$, the zeroth-order general Randi\' c index is defined as $^0R_{\alpha}(G)=\sum_{v\in V(G)}d(v)^{\alpha}$, where $d(v)$ is the degree of the vertex $v$ and $\alpha\ne0$ is a real number. The $k$-generalized quasi-tree is a connected graph $G$ with a subset $V_k\subset V(G)$, where $|V_k|=k$ such that $G-V_k$ is a tree, but for any subset $V_{k-1}\subset V(G)$ with cardinality $k-1$, $G-V_{k-1}$ is not a tree. In this paper, we characterize the extremal $k$-generalized quasi trees with the minimum and maximum values of the zeroth-order general Randi\' c index for $\alpha\neq 0$.


title: Almost Every Simply Typed Lambda-Term Has a Long Beta-Reduction Sequence

abstract: It is well known that the length of a beta-reduction sequence of a simply typed lambda-term of order k can be huge; it is as large as k-fold exponential in the size of the lambda-term in the worst case. We consider the following relevant question about quantitative properties, instead of the worst case: how many simply typed lambda-terms have very long reduction sequences? We provide a partial answer to this question, by showing that asymptotically almost every simply typed lambda-term of order k has a reduction sequence as long as (k-1)-fold exponential in the term size, under the assumption that the arity of functions and the number of variables that may occur in every subterm are bounded above by a constant. To prove it, we have extended the infinite monkey theorem for strings to a parametrized one for regular tree languages, which may be of independent interest. The work has been motivated by quantitative analysis of the complexity of higher-order model checking.


title: Words have bounded width in $SL(n,\mathbb{Z})$

abstract: We prove two results about width of words in $SL_n(\mathbb{Z})$. The first is that, for every $n \geq 3$, there is a constant $C(n)$ such that the width of any word in $SL_n(\mathbb{Z})$ is less than $C(n)$. The second result is that, for any word $w$, if $n$ is big enough, the width of $w$ in $SL_n(\mathbb{Z})$ is at most 87.


title: Synchronization of Two Diffusively Coupled Chaotic Parametrically Excited nonidentical Pendula

abstract: Auxiliary system approach and various nearest neighbor methods are widely used to detect generalized synchronization in non-identical coupled systems. These methods generally give contradictory results. Therefore one method alone is not sufficient to predict correct result. We show in this report that it is necessary to apply multiple methods together to come to a conclusion. These methods show a signature of generalized synchronization in diffusively coupled non-identical chaotic parametric excited pendula. But we finally find it to be the almost synchronization. It is achieved when the second Lyapunov exponent and both the system's transverse Lyapunov exponents are almost equal. The transition from asynchronous state to almost synchronization is through frequency entrainment as coupling constant is increased. Non-identity of the pendula are realized by mismatch in amplitude of parametric forcing. The frequency entrainment regime does not depend on amplitude mismatch whereas onset of almost synchronization increases with increase in mismatch. The systems


title: A complementary compact laser based neutron source

abstract: Several experiments of neutron generation using high intensity laser sources, with a power exceeding 10^19W/cm^2 via TNSA (Target Normal Sheath Acceleration) or other similar methods, have been performed in the past years in different laboratories. However, so far there is no one running neutron source based on such a technology. In the framework of the Conceptual Report Design of a new accelerator in the Eupraxia project we are studying the possibility to have a laser-based neutron source, not only by TNSA but also from self-injection schemes. We focus our attention on the applications in cultural heritage studies as well also on the complementary role that such a source can have in the framework of large facilities devoted to radiation production.


title: HoPP: Robust and Resilient Publish-Subscribe for an Information-Centric Internet of Things

abstract: This paper revisits NDN deployment in the IoT with a special focus on the interaction of sensors and actuators. Such scenarios require high responsiveness and limited control state at the constrained nodes. We argue that the NDN request-response pattern which prevents data push is vital for IoT networks. We contribute HoP-and-Pull (HoPP), a robust publish-subscribe scheme for typical IoT scenarios that targets IoT networks consisting of hundreds of resource constrained devices at intermittent connectivity. Our approach limits the FIB tables to a minimum and naturally supports mobility, temporary network partitioning, data aggregation and near real-time reactivity. We experimentally evaluate the protocol in a real-world deployment using the IoT-Lab testbed with varying numbers of constrained devices, each wirelessly interconnected via IEEE 802.15.4 LowPANs. Implementations are built on CCN-lite with RIOT and support experiments using various single- and multi-hop scenarios.


title: The role of non-spherical double counting in DFT+DMFT: total energy and structural optimization of pnictide superconductors

abstract: A simple scheme for avoiding non-spherical double counting in the combination of density func- tional theory with dynamical mean-field theory (DFT+DMFT)is developed. It is applied to total- energy calculations and structural optimization of the pnictide superconductor LaFeAsO. The results are compared to a recently proposed "exact" double-counting formulation. Both schemes bring the optimized Fe-As interatomic distance close to the experimental value. This resolves the long stand- ing controversy between DFT+DMFT and experiment for the structural optimization of LaFeAsO.


title: Privacy in Index Coding: Improved Bounds and Coding Schemes

abstract: It was recently observed in [1], that in index coding, learning the coding matrix used by the server can pose privacy concerns: curious clients can extract information about the requests and side information of other clients. One approach to mitigate such concerns is the use of $k$-limited-access schemes [1], that restrict each client to learn only part of the index coding matrix, and in particular, at most $k$ rows. These schemes transform a linear index coding matrix of rank $T$ to an alternate one, such that each client needs to learn at most $k$ of the coding matrix rows to decode its requested message. This paper analyzes $k$-limited-access schemes. First, a worst-case scenario, where the total number of clients $n$ is $2^T-1$ is studied. For this case, a novel construction of the coding matrix is provided and shown to be order-optimal in the number of transmissions. Then, the case of a general $n$ is considered and two different schemes are designed and analytically and numerically assessed in their performance. It is shown that these schemes perform better than the one designed for the case $n=2^T-1$.


title: Structure of photodissociation fronts in star-forming regions revealed by observations of high-J CO emission lines with Herschel

abstract: In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces. We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code. A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to Jup=23 in the Orion Bar (Jup=19 in NGC7023), can only originate from small structures of typical thickness of a few 1e-3 pc and at high thermal pressures (Pth~1e8 K cm-3). Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds.


title: Stability in the homology of Deligne-Mumford compactifications

abstract: Using the the theory of FS^op modules, we study the asymptotic behavior of the homology of $\overline M_{g,n}$, the Deligne--Mumford compactification of the moduli space of curves, for $n >> 0$. An FS^op module is a contravariant functor from the category of finite sets and surjections to vector spaces. Via maps that glue on marked P^1's, we give the homology of $\overline M_{g,n}$ the structure of an FS^op module and bound its degree of generation. As a consequence, we prove that the generating function $\sum_{n} \dim(H_i(\overline M_{g,n})) t^n$ is rational, and its denominator has roots in the set $\{1, 1/2, \dots, 1/p(g,i)\}$ where $p(g,i)$ is a polynomial of order $O(g^2 i^2)$. We also obtain restrictions on the decomposition of the homology of $\overline M_{g,n}$ into irreducible $S_n$ representations.


title: On Locally Decodable Index Codes

abstract: Index coding achieves bandwidth savings by jointly encoding the messages demanded by all the clients in a broadcast channel. The encoding is performed in such a way that each client can retrieve its demanded message from its side information and the broadcast codeword. In general, in order to decode its demanded message symbol, a receiver may have to observe the entire transmitted codeword. Querying or downloading the codeword symbols might involve costs to a client -- such as network utilization costs and storage requirements for the queried symbols to perform decoding. In traditional index coding solutions, this 'client aware' perspective is not considered during code design. As a result, for these codes, the number of codeword symbols queried by a client per decoded message symbol, which we refer to as 'locality', could be large. In this paper, considering locality as a cost parameter, we view index coding as a trade-off between the achievable broadcast rate (codeword length normalized by the message length) and locality, where the objective is to minimize the broadcast rate for a given value of locality and vice versa. We show that the smallest possible locality for any index coding problem is 1, and that the optimal index coding solution with locality 1 is the coding scheme based on fractional coloring of the interference graph. We propose index coding schemes with small locality by covering the side information graph using acyclic subgraphs and subgraphs with small minrank. We also show how locality can be accounted for in conventional partition multicast and cycle covering solutions to index coding. Finally, applying these new techniques, we characterize the locality-broadcast rate trade-off of the index coding problem whose side information graph is the directed 3-cycle.


title: Robust inference with knockoffs

abstract: We consider the variable selection problem, which seeks to identify important variables influencing a response $Y$ out of many candidate features $X_1, \ldots, X_p$. We wish to do so while offering finite-sample guarantees about the fraction of false positives - selected variables $X_j$ that in fact have no effect on $Y$ after the other features are known. When the number of features $p$ is large (perhaps even larger than the sample size $n$), and we have no prior knowledge regarding the type of dependence between $Y$ and $X$, the model-X knockoffs framework nonetheless allows us to select a model with a guaranteed bound on the false discovery rate, as long as the distribution of the feature vector $X=(X_1,\dots,X_p)$ is exactly known. This model selection procedure operates by constructing "knockoff copies'" of each of the $p$ features, which are then used as a control group to ensure that the model selection algorithm is not choosing too many irrelevant features. In this work, we study the practical setting where the distribution of $X$ could only be estimated, rather than known exactly, and the knockoff copies of the $X_j$'s are therefore constructed somewhat incorrectly. Our results, which are free of any modeling assumption whatsoever, show that the resulting model selection procedure incurs an inflation of the false discovery rate that is proportional to our errors in estimating the distribution of each feature $X_j$ conditional on the remaining features $\{X_k:k\neq j\}$. The model-X knockoff framework is therefore robust to errors in the underlying assumptions on the distribution of $X$, making it an effective method for many practical applications, such as genome-wide association studies, where the underlying distribution on the features $X_1,\dots,X_p$ is estimated accurately but not known exactly.


title: Cloud Quantum Computing of an Atomic Nucleus

abstract: We report a quantum simulation of the deuteron binding energy on quantum processors accessed via cloud servers. We use a Hamiltonian from pionless effective field theory at leading order. We design a low-depth version of the unitary coupled-cluster ansatz, use the variational quantum eigensolver algorithm, and compute the binding energy to within a few percent. Our work is the first step towards scalable nuclear structure computations on a quantum processor via the cloud, and it sheds light on how to map scientific computing applications onto nascent quantum devices.


title: A parallel workload has extreme variability in a production environment

abstract: Writing data in parallel is a common operation in some computing environments and a good proxy for a number of other parallel processing patterns. The duration of time taken to write data in large-scale compute environments can vary considerably. This variation comes from a number of sources, both systematic and transient. The result is a highly complex behavior that is difficult to characterize. This paper further develops the model for parallel task variability proposed in the paper "A parallel workload has extreme variability" (Henwood et. al 2016). This model is the Generalized Extreme Value (GEV) distribution. This paper further develops the systematic analysis that leads to the GEV model with the addition of a traffic congestion term. Observations of a parallel workload are presented from a High Performance Computing environment under typical production conditions, which include traffic congestion. An analysis of the workload is performed and shows the variability tends towards GEV as the order of parallelism is increased. The results are presented in the context of Amdahl's law and the predictive properties of a GEV models are discussed. A optimization for certain machine designs is also suggested.


title: Connecting the Kuramoto Model and the Chimera State

abstract: Since its discovery in 2002, the chimera state has frequently been described as a counter-intuitive, puzzling phenomenon. The Kuramoto model, in contrast, has become a celebrated paradigm useful for understanding a range of phenomena related to phase transitions, synchronization and network effects. Here we show that the chimera state can be understood as emerging naturally through a symmetry breaking bifurcation from the Kuramoto model's partially synchronized state. Our analysis sheds light on recent observations of chimera states in laser arrays, chemical oscillators, and mechanical pendula.

categories: nlin.PS math-ph math.DS math.MP nlin.AO nlin.CD


title: The Orbit of the Companion to HD 100453A: Binary-Driven Spiral Arms in a Protoplanetary Disk

abstract: HD 100453AB is a 10+/-2 Myr old binary whose protoplanetary disk was recently revealed to host a global two-armed spiral structure. Given the relatively small projected separation of the binary (1.05", or ~108 au), gravitational perturbations by the binary seemed to be a likely driving force behind the formation of the spiral arms. However, the orbit of these stars remained poorly understood, which prevented a proper treatment of the dynamical influence of the companion on the disk. We observed HD 100453AB between 2015-2017 utilizing extreme adaptive optics systems on the Very Large Telescope and Magellan Clay Telescope. We combined the astrometry from these observations with published data to constrain the parameters of the binary's orbit to a=1.06"+/-0.09", e=0.17+/-0.07, and i=32.5+/- 6.5 degrees. We utilized publicly available ALMA CO data to constrain the inclination of the disk to i~28 degrees, which is relatively co-planar with the orbit of the companion and consistent with previous estimates from scattered light images. Finally, we input these constraints into hydrodynamical and radiative transfer simulations to model the structural evolution of the disk. We find that the spiral structure and truncation of the circumprimary disk in HD 100453 are consistent with a companion-dirven origin. Furthermore, we find that the primary star's rotation, its outer disk, and the companion exhibit roughly the same direction of angular momentum, and thus the system likely formed from the same parent body of material.


title: Maximum Likelihood for Gaussian Process Classification and Generalized Linear Mixed Models under Case-Control Sampling

abstract: Modern data sets in various domains often include units that were sampled non-randomly from the population and have a latent correlation structure. Here we investigate a common form of this setting, where every unit is associated with a latent variable, all latent variables are correlated, and the probability of sampling a unit depends on its response. Such settings often arise in case-control studies, where the sampled units are correlated due to spatial proximity, family relations, or other sources of relatedness. Maximum likelihood estimation in such settings is challenging from both a computational and statistical perspective, necessitating approximations that take the sampling scheme into account. We propose a family of approximate likelihood approaches which combine composite likelihood and expectation propagation. We demonstrate the efficacy of our solutions via extensive simulations. We utilize them to investigate the genetic architecture of several complex disorders collected in case-control genetic association studies, where hundreds of thousands of genetic variants are measured for every individual, and the underlying disease liabilities of individuals are correlated due to genetic similarity. Our work is the first to provide a tractable likelihood-based solution for case-control data with complex dependency structures.


title: Central exclusive diffractive production of $p \bar{p}$ pairs in proton-proton collisions at high energies

abstract: We consider the central exclusive production of the $p\bar{p}$ in the continuum and via resonances in proton-proton collisions at high energies. We discuss the diffractive mechanism calculated within the tensor-pomeron approach including pomeron, odderon, and reggeon exchanges. The theoretical results are discussed in the context of existing WA102 and ISR experimental data and predictions for planned or current experiments at the RHIC and the LHC are presented. The distribution in ${\rm y}_{diff}$, the rapidity distance between proton and antiproton, is particularly interesting. We find a dip at ${\rm y}_{diff} = 0$ for the $p \bar{p}$ production, in contrast to the $\pi^{+}\pi^{-}$ and $K^{+}K^{-}$ production. We predict also the $p \bar{p}$ invariant mass distribution to be less steep than for the pairs of pseudoscalar mesons. We argue that these specific differences for the $p \bar{p}$ production with respect to the pseudoscalar meson pair production can be attributed to the proper treatment of the spin of produced particles. We discuss asymmetries that are due to the interference of $C=+1$ and $C=-1$ amplitudes of $p \bar{p}$ production. We have also calculated the cross section for the $pp \to pp \Lambda \overline{\Lambda}$ reaction. Here, the cross section is smaller but the characteristic feature for $d\sigma/d{\rm y}_{diff}$ is predicted to be similar to $p \bar{p}$ production. The presence of resonances in the $p \bar p$ channel may destroy the dip at ${\rm y}_{diff} = 0$. This opens the possibility to study diffractively produced resonances. We discuss the observables suited for this purpose.


title: Plasmonic and photonic crystal applications of a vector solver for 2D Kerr nonlinear waveguides

abstract: We use our vector Maxwell's nonlinear eigenmode solver to study the stationary solutions in 2D cross-section waveguides with Kerr nonlinear cores. This solver is based on the fixed power algorithm within the finite element method. First, studying nonlinear plasmonic slot waveguides, we demonstrate that, even in the low power regime, 1D studies may not provide accurate and meaningfull results compared to 2D ones. Second, we study at high powers the link between the nonlinear parameter $\gamma_{nl}$ and the change of the nonlinear propagation constant $\Delta \beta$. Third, for a specific type of photonic crystal fiber, we show that a non-trivial interplay between the band-gap edge and the nonlinearity takes place.


title: Towards dynamic interaction-based model

abstract: In this paper, we investigate how dynamic properties of reputation can influence the quality of users ranking. Reputation systems should be based on rules that can guarantee a high level of trust and help identifying unreliable units. To understand the effectiveness of dynamic properties in the evaluation of reputation, we propose our own model (DIB-RM) that is based on three factors: forgetting, cumulative, and activity period. In order to evaluate the model, we use data from StackOverflow, which also has its own reputation model. We estimate similarity of ratings between DIB-RM and the StackOverflow model so to check our hypothesis. We use two values to calculate our metric: DIB-RM reputation and $historical$ reputation. We found that $historical$ reputation gives better metric values. Our preliminary results are presented for different sets of values of the aforementioned factors in order to analyze how effectively the model can be used for modeling reputation systems.


title: Learning and Inferring a Driver's Braking Action in Car-Following Scenarios

abstract: Accurately predicting and inferring a driver's decision to brake is critical for designing warning systems and avoiding collisions. In this paper we focus on predicting a driver's intent to brake in car-following scenarios from a perception-decision-action perspective according to his/her driving history. A learning-based inference method, using onboard data from CAN-Bus, radar and cameras as explanatory variables, is introduced to infer drivers' braking decisions by combining a Gaussian mixture model (GMM) with a hidden Markov model (HMM). The GMM is used to model stochastic relationships among variables, while the HMM is applied to infer drivers' braking actions based on the GMM. Real-case driving data from 49 drivers (more than three years' driving data per driver on average) have been collected from the University of Michigan Safety Pilot Model Deployment database. We compare the GMM-HMM method to a support vector machine (SVM) method and an SVM-Bayesian filtering method. The experimental results are evaluated by employing three performance metrics: accuracy, sensitivity, specificity. The comparison results show that the GMM-HMM obtains the best performance, with an accuracy of 90%, sensitivity of 84%, and specificity of 97%. Thus, we believe that this method has great potential for real-world active safety systems.


title: Reducing Noise in Cosmological N-body Simulations with Neutrinos

abstract: We present a new method for generating initial conditions for numerical cosmological simulations in which massive neutrinos are treated as an extra set of N-body (collisionless) particles. It allows us to accurately follow the density field for both Cold Dark Matter (CDM) and neutrinos at both high and low redshifts. At high redshifts, the new method is able to reduce the shot noise in the neutrino power spectrum by a factor of more than $10^7$ compared to previous methods, where the power spectrum was dominated by shot noise at all scales. We find that our new approach also helps to reduce the noise on the total matter power spectrum on large scales, whereas on small scales the results agree with previous simulations. Our new method also allows for a systematic study of clustering of the low velocity tail of the distribution function of neutrinos. This method also allows for the study of the evolution of the overall velocity distribution as a function of the environment determined by the CDM field.


title: Quantum Spin of Elastic Wave

abstract: Unveiling intrinsic spins of propagating waves usually offers people a fundamental understanding of the geometrical and topological properties of waves from classical to quantum aspects. A great variety of research has shown that transverse waves can possess non-trivial quantum spins and topology without help of strong wave-matter interaction. However, until now we still lack essential physical insights about the spin and topological nature of longitudinal waves. Here, demonstrated by elastic waves we uncover unique quantum spins for longitudinal waves and the mixed longitudinal-transverse waves that play essential roles in topological spin-momentum locking. Based on this quantum spin perspective, several abnormal phenomena beyond pure transverse waves are attributed to the hybrid spin induced by mixed longitudinal-transverse waves. The intrinsic hybrid spin reveals the complex spin essence in elastic waves and advances our understanding about their fundamental topological properties. We also show these spin-dependent phenomena can be exploited to control the wave propagation, such as non-symmetric elastic wave excitation by spin pairs, uni-directional Rayleigh wave and spin-selected elastic wave routing. These findings are generally applicable for arbitrary waves with longitudinal and transverse components.


title: Homogeneous length functions on groups

abstract: A pseudo-length function defined on an arbitrary group $G = (G,\cdot,e, (\,)^{-1})$ is a map $\ell: G \to [0,+\infty)$ obeying $\ell(e)=0$, the symmetry property $\ell(x^{-1}) = \ell(x)$, and the triangle inequality $\ell(xy) \leqslant \ell(x) + \ell(y)$ for all $x,y \in G$. We consider pseudo-length functions which saturate the triangle inequality whenever $x=y$, or equivalently those that are homogeneous in the sense that $\ell(x^n) = n\,\ell(x)$ for all $n\in\mathbb{N}$. We show that this implies that $\ell([x,y])=0$ for all $x,y \in G$. This leads to a classification of such pseudo-length functions as pullbacks from embeddings into a Banach space. We also obtain a quantitative version of our main result which allows for defects in the triangle inequality or the homogeneity property.


title: Nonlinear water waves in shallow water in the presence of constant vorticity: A Whitham approach

abstract: Two-dimensional nonlinear gravity waves travelling in shallow water on a vertically sheared current of constant vorticity are considered. Using Euler equations, in the shallow water approximation, hyperbolic equations for the surface elevation and the horizontal velocity are derived. Using Riemann invariants of these equations, that are obtained analytically, a closed-form nonlinear evolution equation for the surface elevation is derived. A dispersive term is added to this equation using the exact linear dispersion relation. With this new single first-order partial differential equation, vorticity effects on undular bores are studied. Within the framework of weakly nonlinear waves, a KdV-type equation and a Whitham equation with constant vorticity are derived from this new model and the effect of vorticity on solitary waves and periodic waves is considered. Futhermore, within the framework of the new model and the Whitham equation a study of the effect of vorticity on the breaking time of dispersive waves and hyperbolic waves as well is carried out.


title: Multi-view Consistency as Supervisory Signal for Learning Shape and Pose Prediction

abstract: We present a framework for learning single-view shape and pose prediction without using direct supervision for either. Our approach allows leveraging multi-view observations from unknown poses as supervisory signal during training. Our proposed training setup enforces geometric consistency between the independently predicted shape and pose from two views of the same instance. We consequently learn to predict shape in an emergent canonical (view-agnostic) frame along with a corresponding pose predictor. We show empirical and qualitative results using the ShapeNet dataset and observe encouragingly competitive performance to previous techniques which rely on stronger forms of supervision. We also demonstrate the applicability of our framework in a realistic setting which is beyond the scope of existing techniques: using a training dataset comprised of online product images where the underlying shape and pose are unknown.


title: Stochastic Learning of Nonstationary Kernels for Natural Language Modeling

abstract: Natural language processing often involves computations with semantic or syntactic graphs to facilitate sophisticated reasoning based on structural relationships. While convolution kernels provide a powerful tool for comparing graph structure based on node (word) level relationships, they are difficult to customize and can be computationally expensive. We propose a generalization of convolution kernels, with a nonstationary model, for better expressibility of natural languages in supervised settings. For a scalable learning of the parameters introduced with our model, we propose a novel algorithm that leverages stochastic sampling on k-nearest neighbor graphs, along with approximations based on locality-sensitive hashing. We demonstrate the advantages of our approach on a challenging real-world (structured inference) problem of automatically extracting biological models from the text of scientific papers.


title: Phase transitions in social networks inspired by the Schelling model

abstract: We propose two models of social segregation inspired by the Schelling model. Agents in our models are nodes of evolving social networks. The total number of social connections of each node remains constant in time, though may vary from one node to the other. The first model describes a "polychromatic" society, in which colors designate different social categories of agents. The parameter $\mu$ favors/disfavors connected "monochromatic triads", i.e. connected groups of three individuals \emph{within the same social category}, while the parameter $\nu$ controls the preference of interactions between two individuals \emph{from different social categories}. The polychromatic model has several distinct regimes in $(\mu,\nu)$-parameter space. In $\nu$-dominated region, the phase diagram is characterized by the plateau in the number of the inter-color connections, where the network is bipartite, while in $\mu$-dominated region, the network looks as two weakly connected unicolor clusters. At $\mu>\mu_{crit}$ and $\nu >\nu_{crit}$ two phases are separated by a critical line, while at small values of $\mu$ and $\nu$, a gradual crossover between the two phases occurs. The second "colorless" model describes a society in which the advantage/disadvantage of forming small fully connected communities (short cycles or cliques in a graph) is controlled by a parameter $\gamma$. We analyze the topological structure of a social network in this model and demonstrate that above a critical threshold, $\gamma^+>0$, the entire network splits into a set of weakly connected clusters, while below another threshold, $\gamma^-<0$, the network acquires a bipartite graph structure. Our results propose mechanisms of formation of self-organized communities in international communication between countries, as well as in crime clans and prehistoric societies.


title: Moduli spaces of real projective structures on surfaces: Notes on a paper by V.V. Fock and A.B. Goncharov

abstract: These notes grew out of our learning and applying the methods of Fock and Goncharov concerning moduli spaces of real projective structures on surfaces with ideal triangulations. We give a self-contained treatment of Fock and Goncharov's description of the moduli space of framed marked properly convex projective structures with minimal or maximal ends, and deduce results of Marquis and Goldman as consequences. We also discuss the Poisson structure on moduli space and its relationship to Goldman's Poisson structure on the character variety.


title: $L^1$ semigroup generation for Fokker-Planck operators associated with general L\'evy driven SDEs

abstract: We prove a new generation result in $L^1$ for a large class of non-local operators with non-degenerate local terms. This class contains the operators appearing in Fokker-Planck or Kolmogorov forward equations associated with L\'evy driven SDEs, i.e. the adjoint operators of the infinitesimal generators of these SDEs. As a byproduct, we also obtain a new elliptic regularity result of independent interest. The main novelty in this paper is that we can consider very general L\'evy operators, including state-space depending coefficients with linear growth and general L\'evy measures which can be singular and have fat tails.


title: BioWorkbench: A High-Performance Framework for Managing and Analyzing Bioinformatics Experiments

abstract: Advances in sequencing techniques have led to exponential growth in biological data, demanding the development of large-scale bioinformatics experiments. Because these experiments are computation- and data-intensive, they require high-performance computing (HPC) techniques and can benefit from specialized technologies such as Scientific Workflow Management Systems (SWfMS) and databases. In this work, we present BioWorkbench, a framework for managing and analyzing bioinformatics experiments. This framework automatically collects provenance data, including both performance data from workflow execution and data from the scientific domain of the workflow application. Provenance data can be analyzed through a web application that abstracts a set of queries to the provenance database, simplifying access to provenance information. We evaluate BioWorkbench using three case studies: SwiftPhylo, a phylogenetic tree assembly workflow; SwiftGECKO, a comparative genomics workflow; and RASflow, a RASopathy analysis workflow. We analyze each workflow from both computational and scientific domain perspectives, by using queries to a provenance and annotation database. Some of these queries are available as a pre-built feature of the BioWorkbench web application. Through the provenance data, we show that the framework is scalable and achieves high-performance, reducing up to 98% of the case studies execution time. We also show how the application of machine learning techniques can enrich the analysis process.


abstract: We propose 5-brane webs for 5d $\mathcal{N}=1$ $G_2$ gauge theories. From a Higgsing of the $SO(7)$ gauge theory with a hypermultiplet in the spinor representation, we construct two types of 5-brane web configurations for the pure $G_2$ gauge theory using an O5-plane or an $\widetilde{\text{O5}}$-plane. Adding flavors to the 5-brane web for the pure $G_2$ gauge theory is also discussed. Based on the obtained 5-brane webs, we compute the partition functions for the 5d $G_2$ gauge theories using the recently suggested topological vertex formulation with an O5-plane, and we find agreement with known results.


title: Factorization Theorem Relating Euclidean and Light-Cone Parton Distributions

abstract: In a large-momentum nucleon state, the matrix element of a gauge-invariant Euclidean Wilson line operator accessible from lattice QCD can be related to the standard light-cone parton distribution function through the large-momentum effective theory (LaMET) expansion. This relation is given by a factorization theorem with a non-trivial matching coefficient. Using the operator product expansion we prove the large-momentum factorization of the quasi-parton distribution function in LaMET, and show that the more recently discussed Ioffe-time distribution approach also obeys an equivalent factorization theorem. Explicit results for the coefficients are obtained and compared at one-loop. Our proof clearly demonstrates that the matching coefficients in the $\overline{\rm MS}$ scheme depend on the large partonic momentum rather than the nucleon momentum.


title: Black Holes as Brains: Neural Networks with Area Law Entropy

abstract: Motivated by the potential similarities between the underlying mechanisms of the enhanced memory storage capacity in black holes and in brain networks, we construct an artificial quantum neural network based on gravity-like synaptic connections and a symmetry structure that allows to describe the network in terms of geometry of a d-dimensional space. We show that the network possesses a critical state in which the gapless neurons emerge that appear to inhabit a (d-1)-dimensional surface, with their number given by the surface area. In the excitations of these neurons, the network can store and retrieve an exponentially large number of patterns within an arbitrarily narrow energy gap. The corresponding micro-state entropy of the brain network exhibits an area law. The neural network can be described in terms of a quantum field, via identifying the different neurons with the different momentum modes of the field, while identifying the synaptic connections among the neurons with the interactions among the corresponding momentum modes. Such a mapping allows to attribute a well-defined sense of geometry to an intrinsically non-local system, such as the neural network, and vice versa, it allows to represent the quantum field model as a neural network.


title: Genuine Distributed Coherence

abstract: We introduce a notion of genuine distributed coherence. Such a notion is based on the possibility of concentrating on individual systems the coherence present in a distributed system, by making use of incoherent unitary transformations. We define an entropic quantifier of genuine distributed multipartite coherence for generic mixed states, and we focus on the bipartite pure-state case. In the latter case we derive necessary and sufficient conditions for the possibility of fully localizing the coherence, hence identifying the conditions for genuine distributed bipartite coherence. We analyze in detail the quantitative problem for the case of two-qubit pure states, identifying the states with the largest amount of genuine distributed coherence. Interestingly, such states do not have maximal global coherence nor maximal coherence rank.


title: Inner Bound for the Capacity Region of Noisy Channels with an Authentication Requirement

abstract: The rate regions of many variations of the standard and wire-tap channels have been thoroughly explored. Secrecy capacity characterizes the loss of rate required to ensure that the adversary gains no information about the transmissions. Authentication does not have a standard metric, despite being an important counterpart to secrecy. While some results have taken an information-theoretic approach to the problem of authentication coding, the full rate region and accompanying trade-offs have yet to be characterized. In this paper, we provide an inner bound of achievable rates with an average authentication and reliability constraint. The bound is established by combining and analyzing two existing authentication schemes for both noisy and noiseless channels. We find that our coding scheme improves upon existing schemes.


title: How nonuniform contact profiles of T cell receptors modulate thymic selection outcomes

abstract: T cell receptors (TCRs) bind foreign or self-peptides attached to major histocompatibility complex (MHC) molecules, and the strength of this interaction determines T cell activation. Optimizing the ability of T cells to recognize a diversity of foreign peptides yet be tolerant of self-peptides is crucial for the adaptive immune system to properly function. This is achieved by selection of T cells in the thymus, where immature T cells expressing unique, stochastically generated TCRs interact with a large number of self-peptide-MHC; if a TCR does not bind strongly enough to any self-peptide-MHC, or too strongly with at least one self-peptide-MHC, the T cell dies. Past theoretical work cast thymic selection as an extreme value problem, and characterized the statistical enrichment or depletion of amino acids in the post-selection TCR repertoire, showing how T cells are selected to be able to specifically recognize peptides derived from diverse pathogens, yet have limited self-reactivity. Here, we investigate how the degree of enrichment is modified by nonuniform contacts that a TCR makes with peptide-MHC. Specifically, we were motivated by recent experiments showing that amino acids at certain positions of a TCR sequence have large effects on thymic selection outcomes, and crystal structure data that reveal a nonuniform contact profile between a TCR and its peptide-MHC ligand. Using a representative TCR contact profile as an illustration, we show via simulations that the degree of enrichment now varies by position according to the contact profile, and, importantly, it depends on the implementation of nonuniform contacts during thymic selection. We explain these nontrivial results analytically. Our study has implications for understanding the selection forces that shape the functionality of the post-selection TCR repertoire.

categories: q-bio.PE cond-mat.stat-mech physics.bio-ph q-bio.QM


title: Quantum algorithm for simulating real time evolution of lattice Hamiltonians

abstract: We study the problem of simulating the time evolution of a lattice Hamiltonian, where the qubits are laid out on a lattice and the Hamiltonian only includes geometrically local interactions (i.e., a qubit may only interact with qubits in its vicinity). This class of Hamiltonians is very general and is believed to capture fundamental interactions of physics. Our algorithm simulates the time evolution of such a Hamiltonian on $n$ qubits for time $T$ up to error $\epsilon$ using $\mathcal O( nT \mathrm{polylog} (nT/\epsilon))$ gates with depth $\mathcal O(T \mathrm{polylog} (nT/\epsilon))$. Our algorithm is the first simulation algorithm that achieves gate cost quasilinear in $nT$ and polylogarithmic in $1/\epsilon$. Our algorithm also readily generalizes to time-dependent Hamiltonians and yields an algorithm with similar gate count for any piecewise slowly varying time-dependent bounded local Hamiltonian. We also prove a matching lower bound on the gate count of such a simulation, showing that any quantum algorithm that can simulate a piecewise constant bounded local Hamiltonian in one dimension to constant error requires $\tilde \Omega(nT)$ gates in the worst case. The lower bound holds even if we only require the output state to be correct on local measurements. To our best knowledge, this is the first nontrivial lower bound on the gate complexity of the simulation problem. Our algorithm is based on a decomposition of the time-evolution unitary into a product of small unitaries using Lieb-Robinson bounds. In the appendix, we prove a Lieb-Robinson bound tailored to Hamiltonians with small commutators between local terms, giving zero Lieb-Robinson velocity in the limit of commuting Hamiltonians. This improves the performance of our algorithm when the Hamiltonian is close to commuting.


title: Rotational and fine structure of open-shell molecules in nearly degenerate electronic states

abstract: An effective Hamiltonian without symmetry restriction has been developed to model the rotational and fine structure of two nearly degenerate electronic states of an open-shell molecule. In addition to the rotational Hamiltonian for an asymmetric top, this spectroscopic model includes energy separation between the two states due to difference potential and zero-point energy difference, as well as the spin-orbit (SO), Coriolis, and electron spin-molecular rotation (SR) interactions. Hamiltonian matrices are computed using orbitally and fully symmetrized case (a) and case (b) basis sets. Intensity formulae and selection rules for rotational transitions between a pair of nearly degenerate states and a nondegenerate state have also been derived using all four basis sets. It is demonstrated using real examples of free radicals that the fine structure of a single electronic state can be simulated with either a SR tensor or a combination of SO and Coriolis constants. The related molecular constants can be determined precisely only when all interacting levels are simulated simultaneously. The present study suggests that analysis of rotational and fine structure can provide quantitative insights into vibronic interactions and related effects.


title: The Unreasonable Effectiveness of Deep Features as a Perceptual Metric

abstract: While it is nearly effortless for humans to quickly assess the perceptual similarity between two images, the underlying processes are thought to be quite complex. Despite this, the most widely used perceptual metrics today, such as PSNR and SSIM, are simple, shallow functions, and fail to account for many nuances of human perception. Recently, the deep learning community has found that features of the VGG network trained on ImageNet classification has been remarkably useful as a training loss for image synthesis. But how perceptual are these so-called "perceptual losses"? What elements are critical for their success? To answer these questions, we introduce a new dataset of human perceptual similarity judgments. We systematically evaluate deep features across different architectures and tasks and compare them with classic metrics. We find that deep features outperform all previous metrics by large margins on our dataset. More surprisingly, this result is not restricted to ImageNet-trained VGG features, but holds across different deep architectures and levels of supervision (supervised, self-supervised, or even unsupervised). Our results suggest that perceptual similarity is an emergent property shared across deep visual representations.


title: A Note on a Fourier Expansion of Automorphic Functions that Nontrivially Involves the Euclidean Algorithm

abstract: We extend the Fourier expansion of automorphic functions of Piateski-Shapiro and Shalika in a way that the Euclidean algorithm appears nontrivially. Then, we discuss a potential application in the Rankin-Selberg method.


title: Axionic Landscape for Higgs Near-Criticality

abstract: The measured value of the Higgs quartic coupling $\lambda$ is peculiarly close to the critical value above which the Higgs potential becomes unstable, when extrapolated to high scales by renormalization group running. It is tempting to speculate that there is an anthropic reason behind this near-criticality. We show how an axionic field can provide a landscape of vacuum states in which $\lambda$ scans. These states are populated during inflation to create a multiverse with different quartic couplings, with a probability distribution $P$ that can be computed. If $P$ is peaked in the anthropically forbidden region of Higgs instability, then the most probable universe compatible with observers would be close to the boundary, as observed. We discuss three scenarios depending on the Higgs vacuum selection mechanism: decay by quantum tunneling; by thermal fluctuations or by inflationary fluctuations.


title: Mass Pattern of the SM Fermions: Flavor Democracy Revisited

abstract: Mass pattern of the SM fermions is one of the most important mysterious in particle physics. Flavor Democracy could shed light on this mystery. Addition of isosinglet quark and isosinglet lepton give opportunity to obtain masses of charged leptons and quarks of the 2nd and 3rd family due to small deviations of full Flavor Democracy.


title: The efficiency of community detection by most similar node pairs

abstract: Community analysis is an important way to ascertain whether or not a complex system consists of sub-structures with different properties. In this paper, we give a two level community structure analysis for the SSCI journal system by most similar co-citation pattern. Five different strategies for the selection of most similar node (journal) pairs are introduced. The efficiency is checked by the normalized mutual information technique. Statistical properties and comparisons of the community results show that both of the two level detection could give instructional information for the community structure of complex systems. Further comparisons of the five strategies indicates that, the most efficient strategy is to assign nodes with maximum similarity into the same community whether the similarity information is complete or not, while random selection generates small world local community with no inside order. These results give valuable indication for efficient community detection by most similar node pairs.


title: Formalized Conceptual Spaces with a Geometric Representation of Correlations

abstract: The highly influential framework of conceptual spaces provides a geometric way of representing knowledge. Instances are represented by points in a similarity space and concepts are represented by convex regions in this space. After pointing out a problem with the convexity requirement, we propose a formalization of conceptual spaces based on fuzzy star-shaped sets. Our formalization uses a parametric definition of concepts and extends the original framework by adding means to represent correlations between different domains in a geometric way. Moreover, we define various operations for our formalization, both for creating new concepts from old ones and for measuring relations between concepts. We present an illustrative toy-example and sketch a research project on concept formation that is based on both our formalization and its implementation.


title: A substructuring preconditioner with vertex-related interface solvers for elliptic-type equations in three dimensions

abstract: In this paper we propose a variant of the substructuring preconditioner for solving three-dimensional elliptic-type equations with strongly discontinuous coefficients. In the proposed preconditioner, we use the simplest coarse solver associated with the finite element space induced by the coarse partition, and construct vertex-related inexact interface solvers based on overlapping domain decomposition with small overlaps. This new preconditioner has an important merit: its construction and efficiency do not depend on the concrete form of the considered elliptic-type equations. % in the sense that they not only are cheap but also are easy to implement. We apply the proposed preconditioner to solve the linear elasticity problems and Maxwell's equations in three dimensions. Numerical results show that the convergence rate of PCG method with the preconditioner is nearly optimal, and also robust with respect to the (possibly large) jumps of the coefficients in the considered equations.


title: New Results on Multilevel Diversity Coding with Secure Regeneration

abstract: The problem of multilevel diversity coding with secure regeneration is revisited. Under the assumption that the eavesdropper can access the repair data for all compromised storage nodes, Shao el al. provided a precise characterization of the minimum-bandwidth-regeneration (MBR) point of the achievable normalized storage-capacity repair-bandwidth tradeoff region. In this paper, it is shown that the MBR point of the achievable normalized storage-capacity repair-bandwidth tradeoff region remains the same even if we assume that the eavesdropper can access the repair data for some compromised storage nodes (type 2 compromised nodes) but only the data contents of the remaining compromised nodes (type 1 compromised nodes), as long as the number of type 1 compromised nodes is no greater than that of type 2 compromised nodes.


title: Extremals for the Singular Moser-Trudinger Inequality via n-Harmonic Transplantation

abstract: The Moser-Trudinger embedding has been generalized in [Adimurthi A.; Sandeep K., A singular Moser-Trudinger embedding and its applications, \textit{NoDEA Nonlinear Differential Equations Appl.}, 13 (2007), no. 5-6, 585--603] to the following weighted version: if $\Omega\subset\mathbb{R}^n$ is bounded, $\omega_{n-1}$ is the $\mathcal{H}^{n-1}$ measure of the unit sphere, then for $\alpha>0$ and $\beta\in [0,n)$, $$ \sup_{u\in\mathcal{B}_1}\int_{\Omega}\frac{e^{\alpha |u|^{n/(n-1)}}}{|x|^{\beta}}\leq C \ \Leftrightarrow \ \frac{\alpha}{\alpha_n}+\frac{\beta}{n}\leq1,\qquad $$ where $\alpha_n=n\cnn$ and $\mathcal{B}_1 = \left\{ u \in W_0^{1, n}(\Omega) \ | \ \int_{\Omega} |\nabla u |^n \leq1 \right\}$. We prove that the supremum is attained on any domain $\Omega$. The paper also fills in the gaps in the proof of [Lin K.C., Extremal functions for Moser's inequality, \textit{Trans. of. Am. Math. Soc.}, 384 (1996), 2663--2671], which deals with the case $\beta=0.$


title: Gas Kinematics in FIRE Simulated Galaxies Compared to Spatially Unresolved HI Observations

abstract: The shape of a galaxy's spatially unresolved, globally integrated 21-cm emission line depends on its internal gas kinematics: galaxies with rotation-supported gas disks produce double-horned profiles with steep wings, while galaxies with dispersion-supported gas produce Gaussian-like profiles with sloped wings. Using mock observations of simulated galaxies from the FIRE project, we show that one can therefore constrain a galaxy's gas kinematics from its unresolved 21-cm line profile. In particular, we find that the kurtosis of the 21-cm line increases with decreasing $V/\sigma$, and that this trend is robust across a wide range of masses, signal-to-noise ratios, and inclinations. We then quantify the shapes of 21-cm line profiles from a morphologically unbiased sample of $\sim$2000 low-redshift, HI-detected galaxies with $M_{\rm star} = 10^{7-11} M_{\odot}$ and compare to the simulated galaxies. At $M_{\rm star} \gtrsim 10^{10} M_{\odot}$, both the observed and simulated galaxies produce double-horned profiles with low kurtosis and steep wings, consistent with rotation-supported disks. Both the observed and simulated line profiles become more Gaussian-like (higher kurtosis and less-steep wings) at lower masses, indicating increased dispersion support. However, the simulated galaxies transition from rotation to dispersion support more strongly: at $M_{\rm star} = 10^{8-10}M_{\odot}$, most of the simulations produce more Gaussian-like profiles than typical observed galaxies with similar mass, indicating that gas in the low-mass simulated galaxies is, on average, overly dispersion-supported. Most of the lower-mass simulated galaxies also have somewhat lower gas fractions than the median of the observed population. The simulations nevertheless reproduce the observed line-width baryonic Tully-Fisher relation, which is insensitive to rotation vs. dispersion support.


title: Confusing binaries: the role of stellar binaries in biasing disk properties in the Galactic Center

abstract: The population of young stars near the Supermassive black hole (SMBH) in the Galactic Center (GC) has presented an unexpected challenge to theories of star formation. Kinematics measurements of these stars have revealed a stellar disk structure (with an apparent 20% disk membership) that has provided important clues to the origin of these mysterious young stars. However many of the apparent disk properties are difficult to explain, including the low disk membership fraction and the high eccentricities, given the youth of this population. Thus far, all efforts to derive the properties of this disk have made the simplifying assumption that stars at the GC are single stars. Nevertheless, stellar binaries are prevalent in our Galaxy, and recent investigations suggested that they may also be abundant in the Galactic Center. Here we show that binaries in the disk can largely alter the apparent orbital properties of the disk. The motion of binary members around each other adds a velocity component, which can be comparable to the magnitude of the velocity around the SMBH in the GC. Thus neglecting the contribution of binaries can significantly vary the inferred stars' orbital properties. While the disk orientation is unaffected the apparent disk's 2D width is increased to about 11.2deg, similar to the observed width. For a population of stars orbiting the SMBH with zero eccentricity, unaccounted for binaries will create a wide apparent eccentricity distribution with an average of 0.23.This is consistent with the observed average eccentricity of the stars' in the disk. We suggest that this high eccentricity value, which poses a theoretical challenge, may be an artifact of binary stars. Finally our results suggest that the actual disk membership might be significantly higher than the one inferred by observations that ignore the contribution of binaries, alleviating another theoretical challenge.


title: TOPoS: IV. Chemical abundances from high-resolution observations of seven EMP stars

abstract: Extremely metal-poor stars provide us with indirect information on the first generations of massive stars. The TOPoS survey has been designed to increase the census of these stars and to provide a chemical inventory that is as detailed as possible. Seven of the most iron-poor stars have been observed with the UVES spectrograph at the ESO VLT Kueyen 8.2m telescope to refine their chemical composition. We analysed the spectra based on 1D LTE model atmospheres, but also used 3D hydrodynamical simulations of stellar atmospheres. We measured carbon in six of the seven stars: all are carbon-enhanced and belong to the low-carbon band, defined in the TOPoS II paper. We measured lithium (A(Li)=1.9) in the most iron-poor star (SDSS J1035+0641, [Fe/H] < -5.2 ). We were also able to measure Li in three stars at [Fe/H]~ -4.0, two of which lie on the Spite plateau. We confirm that SDSS J1349+1407 is extremely rich in Mg, but not in Ca. It is also very rich in Na. Several of our stars are characterised by low alpha-to-iron ratios. The lack of high-carbon band stars at low metallicity can be understood in terms of evolutionary timescales of binary systems. The detection of Li in SDSS J1035+0641 places a strong constraint on theories that aim at solving the cosmological lithium problem. The Li abundance of the two warmer stars at [Fe/H]~ -4.0 places them on the Spite plateau, while the third, cooler star, lies below. We argue that this suggests that the temperature at which Li depletion begins increases with decreasing [Fe/H]. SDSS J1349+1407 may belong to a class of Mg-rich EMP stars. We cannot assess if there is a scatter in alpha-to-iron ratios among the EMP stars or if there are several discrete populations. However, the existence of stars with low alpha-to-iron ratios is supported by our observations.


title: LINEAR: A Novel Algorithm for Reconstructing Slitless Spectroscopy from HST/WFC3

abstract: We present a grism extraction package (LINEAR) designed to reconstruct one-dimensional spectra from a collection of slitless spectroscopic images, ideally taken at a variety of orientations, dispersion directions, and/or dither positions. Our approach is to enumerate every transformation between all direct image positions (ie. a potential source) and the collection of grism images at all relevant wavelengths. This leads to solving a large, sparse system of linear equations, which we invert using the standard LSQR algorithm. We implement a number of color and geometric corrections (such as flat field, pixel-area map, source morphology, and spectral bandwidth), but assume many effects have been calibrated out (such as basic reductions, background subtraction, and astrometric refinement). We demonstrate the power of our approach with several Monte Carlo simulations and the analysis of archival data. The simulations include astrometric and photometric uncertainties, sky-background estimation, and signal-to-noise calculations. The data are G141 observations obtained with the Wide-Field Camera 3 of the Hubble Ultra-Deep Field, and show the power of our formalism by improving the spectral resolution without sacrificing the signal-to-noise (a tradeoff that is often made by current approaches). Additionally, our approach naturally accounts for source contamination, which is only handled heuristically by present softwares. We conclude with a discussion of various observations where our approach will provide much improved spectral one-dimensional spectra, such as crowded fields (star or galaxy clusters), spatially resolved spectroscopy, or surveys with strict completeness requirements. At present our software is heavily geared for Wide-Field Camera 3 IR, however we plan extend the codebase for additional instruments.


title: Predictive Leptogenesis from Minimal Lepton Flavour Violation

abstract: A predictive Leptogenesis scenario is presented based on the Minimal Lepton Flavour Violation symmetry. In the realisation with three right-handed neutrinos transforming under the same flavour symmetry of the lepton electroweak doublets, lepton masses and PMNS mixing parameters can be described according to the current data, including a large Dirac CP phase. The observed matter-antimatter asymmetry of the Universe can be achieved through Leptogenesis, with the CP asymmetry parameter $\varepsilon$ described in terms of only lepton masses, mixings and phases, plus two real parameters of the low-energy effective description. This is in contrast with the large majority of models present in the literature, where $\varepsilon$ depends on several high-energy parameters, preventing a direct connection between low-energy observables and the baryon to photon ratio today. Recovering the correct amount of baryon asymmetry in the Universe constrains the Majorana phases of the PMNS matrix within specific ranges of values: clear predictions for the neutrinoless double beta decay emerge, representing a potential smoking gun for this framework.


title: Not all stars form in clusters - measuring the kinematics of OB associations with $Gaia$

abstract: It is often stated that star clusters are the fundamental units of star formation and that most (if not all) stars form in dense stellar clusters. In this monolithic formation scenario, low density OB associations are formed from the expansion of gravitationally bound clusters following gas expulsion due to stellar feedback. $N$-body simulations of this process show that OB associations formed this way retain signs of expansion and elevated radial anisotropy over tens of Myr. However, recent theoretical and observational studies suggest that star formation is a hierarchical process, following the fractal nature of natal molecular clouds and allowing the formation of large-scale associations in-situ. We distinguish between these two scenarios by characterising the kinematics of OB associations using the Tycho-$Gaia$ Astrometric Solution catalogue. To this end, we quantify four key kinematic diagnostics: the number ratio of stars with positive radial velocities to those with negative radial velocities, the median radial velocity, the median radial velocity normalised by the tangential velocity, and the radial anisotropy parameter. Each quantity presents a useful diagnostic of whether the association was more compact in the past. We compare these diagnostics to models representing random motion and the expanding products of monolithic cluster formation. None of these diagnostics show evidence of expansion, either from a single cluster or multiple clusters, and the observed kinematics are better represented by a random velocity distribution. This result favours the hierarchical star formation model in which a minority of stars forms in bound clusters and large-scale, hierarchically-structured associations are formed in-situ.


title: Growing up in a megalopolis: Environmental effects on galaxy evolution in a supercluster at z~0.65 in UKIDSS UDS

abstract: We present a large-scale galaxy structure Cl J021734-0513 at z~0.65 discovered in the UKIDSS UDS field, made of ~20 galaxy groups and clusters, spreading over 10 Mpc. We report on a VLT/VIMOS spectroscopic follow-up program that, combined with past spectroscopy, allowed us to confirm four galaxy clusters (M200 ~ 10^14 Msol) and a dozen associated groups and star-forming galaxy overdensities. Two additional filamentary structures at z ~ 0.62 and z ~ 0.69 and foreground and background clusters at 0.6 < z < 0.7 were also confirmed along the line-of-sight. The structure subcomponents are at different formation stages. The clusters have a core dominated by passive galaxies and an established red sequence. The remaining structures are a mix of star-forming galaxy overdensities and forming groups. The presence of quiescent galaxies in the core of the latter shows that `preprocessing' has already happened before the groups fall into their more massive neighbours. Our spectroscopy allows us to derive spectral index measurements e.g. emission/absorption line equivalent widths, strength of the 4000A break, valuable to investigate the star formation history of structure members. Based on these line measurements, we select a population of `poststarburst' galaxies. These galaxies are preferentially found within the virial radius of clusters, supporting a scenario in which their recent quenching could be prompted by gas stripping by the dense intracluster medium. We derive stellar age estimates using MCMC-based spectral fitting for quiescent galaxies and find a correlation between ages and colours/stellar masses which favours a top-down formation scenario of the red sequence. A catalogue of ~650 redshifts in UDS will be released alongside the paper.


title: On the claimed X-shaped structure in the Milky Way bulge

abstract: A number of recent studies have claimed that the double red clump (RC) observed in the Milky Way bulge is a consequence of a giant X-shaped structure. In particular, Ness & Lang (2016) reported a direct detection of a faint X-shaped structure from the residual map of the Wide-Field Infrared Survey Explorer (WISE) bulge image. Here we show, however, that their result is affected substantially by whether the conventional dust extinction correction is applied or not and partly by a bulge model subtracted from the original image. We find that the residuals obtained by subtracting either ellipsoidal or boxy bulge models from the dereddened images show no obvious X-shaped structure. We further show that, even if it is real, the stellar density in the claimed X-shaped structure is way too low to be observed as a strong double RC at l=0.


title: Full-sky ray-tracing simulation of weak lensing using ELUCID simulations: exploring galaxy intrinsic alignment and cosmic shear correlations

abstract: The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail.


title: Average [O II] nebular emission associated with Mg II absorbers: Dependence on Fe II absorption

abstract: We investigate the effect of Fe II equivalent width ($W_{2600}$) and fibre size on the average luminosity of [O II]$\lambda\lambda$3727,3729 nebular emission associated with Mg II absorbers (at $0.55 \le z \le 1.3$) in the composite spectra of quasars obtained with 3 and 2 arcsec fibres in the Sloan Digital Sky Survey. We confirm the presence of strong correlations between [O II] luminosity (L$_{[\rm O~II]}$) and equivalent width ($W_{2796}$) and redshift of Mg II absorbers. However, we show L$_{[\rm O~II]}$ and average luminosity surface density suffers from fibre size effects. More importantly, for a given fibre size the average L$_{[\rm O~II]}$ strongly depends on the equivalent width of Fe II absorption lines and found to be higher for Mg II absorbers with $R \equiv$ $W_{\rm 2600}/W_{\rm 2796}$ $\ge 0.5$. In fact, we show the observed strong correlations of L$_{[\rm O~II]}$ with $W_{2796}$ and $z$ of Mg II absorbers are mainly driven by such systems. Direct [O II] detections also confirm the link between L$_{[\rm O~II]}$ and $R$. Therefore, one has to pay attention to the fibre losses and dependence of redshift evolution of Mg II absorbers on $W_{2600}$ before using them as a luminosity unbiased probe of global star formation rate density. We show that the [O II] nebular emission detected in the stacked spectrum is not dominated by few direct detections (i.e., detections $\ge 3 \sigma$ significant level). On an average the systems with $R$ $\ge 0.5$ and $W_{2796}$ $\ge 2$ \AA\ are more reddened, showing colour excess E($B-V$) $\sim$ 0.02, with respect to the systems with $R$ $< 0.5$ and most likely traces the high H I column density systems.


title: Discovery of a new fundamental plane dictating galaxy cluster evolution from gravitational lensing

abstract: In cold dark matter (CDM) cosmology, objects in the Universe have grown under the effect of gravity of dark matter. The intracluster gas in a galaxy cluster was heated when the dark-matter halo formed through gravitational collapse. The potential energy of the gas was converted to thermal energy through this process. However, this process and the thermodynamic history of the gas have not been clearly characterized in connection with with the formation and evolution of the internal structure of dark-matter halos. Here, we show that observational CLASH data of high-mass galaxy clusters lie on a plane in the three-dimensional logarithmic space of their characteristic radius $r_s$, mass $M_s$, and X-ray temperature $T_X$ with a very small orthogonal scatter. The tight correlation indicates that the gas temperature was determined at a specific cluster formation time, which is encoded in $r_s$ and $M_s$. The plane is tilted with respect to $T_X \propto M_s/r_s$, which is the plane expected in case of simplified virial equilibrium. We show that this tilt can be explained by a similarity solution, which indicates that clusters are not isolated but continuously growing through matter accretion from their outer environments. Numerical simulations reproduce the observed plane and its angle. This result holds independently of the gas physics implemented in the code, revealing the fundamental origin of this plane.


title: A theoretical study of top-mass measurements at the LHC using NLO+PS generators of increasing accuracy

abstract: In this paper we study the theoretical uncertainties in the determination of the top-quark mass using next-to-leading-order (NLO) generators interfaced to parton showers (PS) that have different levels of accuracy. Specifically we consider three generators: one that implements NLO corrections in the production dynamics, one that includes also NLO corrections in top decay in the narrow width approximation, and one that implements NLO corrections for both production and decay including finite-width and interference effects. Since our aim is to provide an assessment of the uncertainties of purely theoretical origin, we consider simplified top-mass related observables that are broadly related to those effectively used by experiments, eventually modelling experimental resolution effects with simple smearing procedures. We estimate the differences in the value of the extracted top mass that would occur due to the use of the three different NLO generators, to the variation of scales, to the choice of parton distribution functions and to the matching procedure. Furthermore, we also consider differences due to the shower and to the modelling of non-perturbative effects by interfacing our NLO generators to both Pythia8.2 and Herwig7.1, with various settings. We find very different results depending upon the adopted shower model. While with Pythia8.2 we find moderate differences between the different NLO+PS generators, with Herwig7.1 we find very large ones. Furthermore, the differences between Pythia8.2 and Herwig7.1 generators are also remarkably large.


title: Molecular reconnaissance of the $\beta$ Pictoris gas disk with the SMA: a low HCN/(CO+CO$_2$) outgassing ratio and predictions for future surveys

abstract: The exocometary origin of CO gas has been confirmed in several extrasolar Kuiper belts, with CO ice abundances consistent with Solar System comets. We here present a molecular survey of the $\beta$ Pictoris belt with the Submillimeter Array (SMA), reporting upper limits for CN, HCN, HCO$^+$, N$_2$H$^+$ and H$_2$CO, as well as for H$_2$S, CH$_3$OH, SiO and DCN from archival ALMA data. Non-detections can be attributed to rapid molecular photodissociation due to the A-star's strong UV flux. CN is the longest-lasting and most easily detectable molecule after CO in this environment. We update our NLTE excitation model to include UV fluorescence, finding it plays a key role in CO and CN excitation, and use it to turn the SMA CN/CO flux ratio constraint into an upper limit of $<2.5$% on the HCN/(CO+CO$_2$) ratio of outgassing rates. This value is consistent with, but at the low end of, the broad range observed in Solar System comets. If sublimation dominates outgassing, then this low value may be caused by decreased outgassing for the less volatile molecule HCN compared to CO. If instead UV photodesorption or collisional vaporization of unbound grains dominates outgassing, then this low ratio of rates would imply a low ice abundance ratio, which would in turn indicate a variation in cometary cyanide abundances across planetary systems. To conclude, we make predictions for future molecular surveys and show that CN and HCN should be readily detectable with ALMA around $\beta$ Pictoris for Solar-System-like exocometary compositions.


title: On the relevance of chaos for halo stars in the solar neighbourhood II

abstract: In a previous paper based on dark matter only simulations we show that, in the approximation of an analytic and static potential describing the strongly triaxial and cuspy shape of Milky Way-sized haloes, diffusion due to chaotic mixing in the neighbourhood of the Sun does not efficiently erase phase space signatures of past accretion events. In this second paper we further explore the effect of chaotic mixing using multicomponent Galactic potential models and solar neighbourhood-like volumes extracted from fully cosmological hydrodynamic simulations, thus naturally accounting for the gravitational potential associated with baryonic components, such as the bulge and disc. Despite the strong change in the global Galactic potentials with respect to those obtained in dark matter only simulations, our results confirm that a large fraction of halo particles evolving on chaotic orbits exhibit their chaotic behaviour after periods of time significantly larger than a Hubble time. In addition, significant diffusion in phase space is not observed on those particles that do exhibit chaotic behaviour within a Hubble time.


title: Optimizing deep-space optical communication under power constraints

abstract: We investigate theoretically the efficiency of deep-space optical communication in the presence of background noise. With decreasing average signal power spectral density, a scaling gap opens up between optimized simple-decoded pulse position modulation and generalized on-off keying with direct detection. The scaling of the latter follows the quantum mechanical capacity of an optical channel with additive Gaussian noise. Efficient communication is found to require a highly imbalanced distribution of instantaneous signal power. This condition can be alleviated through the use of structured receivers which exploit optical interference over multiple time bins to concentrate the signal power before the detection stage.


title: CO and dust properties in the TW Hya disk from high-resolution ALMA observations

abstract: We analyze high angular resolution ALMA observations of the TW Hya disk to place constraints on the CO and dust properties. We present new, sensitive observations of the $^{12}$CO $J = 3-2$ line at a spatial resolution of 8 AU (0\farcs14). The CO emission exhibits a bright inner core, a shoulder at $r\approx70$ AU, and a prominent break in slope at $r\approx90$ AU. Radiative transfer modeling is used to demonstrate that the emission morphology can be reasonably reproduced with a $^{12}$CO column density profile featuring a steep decrease at $r\approx15$ AU and a secondary bump peaking at $r\approx70$ AU. Similar features have been identified in observations of rarer CO isotopologues, which trace heights closer to the midplane. Substructure in the underlying gas distribution or radially varying CO depletion that affects much of the disk's vertical extent may explain the shared emission features of the main CO isotopologues. We also combine archival 1.3 mm and 870 $\mu$m continuum observations to produce a spectral index map at a spatial resolution of 2 AU. The spectral index rises sharply at the continuum emission gaps at radii of 25, 41, and 47 AU. This behavior suggests that the grains within the gaps are no larger than a few millimeters. Outside the continuum gaps, the low spectral index values of $\alpha\approx 2$ indicate either that grains up to centimeter size are present, or that the bright continuum rings are marginally optically thick at millimeter wavelengths.


title: The Need for Speed: Escape velocity and dynamical mass measurements of the Andromeda galaxy

abstract: Our nearest large cosmological neighbour, the Andromeda galaxy (M31), is a dynamical system, and an accurate measurement of its total mass is central to our understanding of its assembly history, the life-cycles of its satellite galaxies, and its role in shaping the Local Group environment. Here, we apply a novel approach to determine the dynamical mass of M31 using high velocity Planetary Nebulae (PNe), establishing a hierarchical Bayesian model united with a scheme to capture potential outliers and marginalize over tracers unknown distances. With this, we derive the escape velocity run of M31 as a function of galacto-centric distance, with both parametric and non-parametric approaches. We determine the escape velocity of M31 to be $470\pm{40}$ km s$^{-1}$ at a galacto-centric distance of 15 kpc, and also, derive the total potential of M31, estimating the virial mass and radius of the galaxy to be $0.8\pm{0.1}\times10^{12}\,M_\odot$ and $240\pm{10}$ kpc, respectively. Our M31 mass is on the low-side of the measured range, this supports the lower expected mass of the M31-Milky Way system from the timing and momentum arguments, satisfying the HI constraint on circular velocity between $10\lesssim R/\textrm{kpc}<35$, and agreeing with the stellar mass Tully-Fisher relation. To place these results in a broader context, we compare them to the key predictions of the $\Lambda{\rm CDM}$ cosmological paradigm, including the stellar-mass-halo-mass and the dark matter halo concentration-virial mass correlation, and finding it to be an outlier to this relation.


title: Jeans that fit: weighing the mass of the Milky Way analogues in the $\Lambda{\rm CDM}$ universe

abstract: The spherical Jeans equation is a widely used tool for dynamical study of gravitating systems in astronomy. Here we test its efficacy in robustly weighing the mass of Milky Way analogues, given they need not be in equilibrium or even spherical. Utilizing Milky Way stellar halos simulated in accordance with $\Lambda{\rm CDM}$ cosmology by Bullock and Johnston (2005) and analysing them under the Jeans formalism, we recover the underlying mass distribution of the parent galaxy, within distance $r/{\rm kpc}\in[10,100]$, with a bias of $\sim12\%$ and a dispersion of $\sim14\%$. Additionally, the mass profiles of triaxial dark matter halos taken from the SURFS simulation, within scaled radius $0.2<r/r_{\rm max}<3$, are measured with a bias of $\sim-2.4\%$ and a dispersion of $\sim10\%$. The obtained dispersion is not because of Poisson noise due to small particle numbers as it is twice the later. We interpret the dispersion to be due to the inherent nature of the $\Lambda{\rm CDM}$ halos, for example being aspherical and out-of-equilibrium. Hence the dispersion obtained for stellar halos sets a limit of about $12\%$ (after adjusting for random uncertainty) on the accuracy with which the mass profiles of the Milky Way-like galaxies can be reconstructed using the spherical Jeans equation. This limit is independent of the quantity and quality of the observational data. The reason for a non zero bias is not clear, hence its interpretation is not obvious at this stage.


title: LDPC Codes with Local and Global Decoding

abstract: This paper presents a theoretical study of a new type of LDPC codes motivated by practical storage applications. LDPCL codes (suffix L represents locality) are LDPC codes that can be decoded either as usual over the full code block, or locally when a smaller sub-block is accessed (to reduce latency). LDPCL codes are designed to maximize the error-correction performance vs. rate in the usual (global) mode, while at the same time providing a certain performance in the local mode. We develop a theoretical framework for the design of LDPCL codes. Our results include a design tool to construct an LDPC code with two data-protection levels: local and global. We derive theoretical results supporting this tool and we show how to achieve capacity with it. A trade-off between the gap to capacity and the number of full-block accesses is studied, and a finite-length analysis of ML decoding is performed to exemplify a trade-off between the locality capability and the full-block error-correcting capability.


title: Numerical Simulations of Multiphase Winds and Fountains from Star-Forming Galactic Disks: I. Solar Neighborhood TIGRESS Model

abstract: Gas blown away from galactic disks by supernova (SN) feedback plays a key role in galaxy evolution. We investigate outflows utilizing the solar neighborhood model of our high-resolution, local galactic disk simulation suite, TIGRESS. In our numerical implementation, star formation and SN feedback are self-consistently treated and well resolved in the multiphase, turbulent, magnetized interstellar medium. Bursts of star formation produce spatially and temporally correlated SNe that drive strong outflows, consisting of hot (T>5x10^5K) winds and warm (5050K < T < 2x10^4K) fountains. The hot gas at distance d>1kpc from the midplane has mass and energy fluxes nearly constant with d. The hot flow escapes our local Cartesian box barely affected by gravity and is expected to accelerate up to the terminal velocity of v_wind~350-500km/s. The mean mass and energy loading factors of the hot wind are 0.1 and 0.02, respectively. For warm gas, the mean outward mass flux through d=1kpc is comparable to the mean star formation rate, but only a small fraction of this gas is at velocity >50km/s. Thus, the warm outflows eventually fall back as inflows. The warm fountain flows are created by expanding hot superbubbles at d< 1kpc; at larger d neither ram pressure acceleration nor cooling transfers significant momentum or energy flux from the hot wind to the warm outflow. The velocity distribution at launching near d~1kpc better represents warm outflows than a single mass loading factor, potentially enabling development of subgrid models for warm galactic winds in arbitrary large-scale galactic potentials.


title: Tractable models of self-sustaining autocatalytic networks

abstract: Self-sustaining autocatalytic networks play a central role in living systems, from metabolism at the origin of life, simple RNA networks, and the modern cell, to ecology and cognition. A collectively autocatalytic network that can be sustained from an ambient food set is also referred to more formally as a `Reflexively Autocatalytic F-generated' (RAF) set. In this paper, we first investigate a simplified setting for studying RAFs, which are nevertheless relevant to real biochemistry and allows for a more exact mathematical analysis based on graph-theoretic concepts. This, in turn, allows for the development of efficient (polynomial-time) algorithms for questions that are computationally NP-hard in the general RAF setting. We then show how this simplified setting for RAF systems leads naturally to a more general notion of RAFs that are `generative' (they can be built up from simpler RAFs) and for which efficient algorithms carry over to this more general setting. Finally, we show how classical RAF theory can be extended to deal with ensembles of catalysts as well as the assignment of rates to reactions according to which catalysts (or combinations of catalysts) are available.


title: Model-Based Action Exploration for Learning Dynamic Motion Skills

abstract: Deep reinforcement learning has achieved great strides in solving challenging motion control tasks. Recently, there has been significant work on methods for exploiting the data gathered during training, but there has been less work on how to best generate the data to learn from. For continuous action domains, the most common method for generating exploratory actions involves sampling from a Gaussian distribution centred around the mean action output by a policy. Although these methods can be quite capable, they do not scale well with the dimensionality of the action space, and can be dangerous to apply on hardware. We consider learning a forward dynamics model to predict the result, ($x_{t+1}$), of taking a particular action, ($u$), given a specific observation of the state, ($x_{t}$). With this model we perform internal look-ahead predictions of outcomes and seek actions we believe have a reasonable chance of success. This method alters the exploratory action space, thereby increasing learning speed and enables higher quality solutions to difficult problems, such as robotic locomotion and juggling.


title: A machine learns to predict the stability of circumbinary planets

abstract: Long-period circumbinary planets appear to be as common as those orbiting single stars and have been found to frequently have orbital radii just beyond the critical distance for dynamical stability. Assessing the stability is typically done either through N-body simulations or using the stability criterion first considered by Dvorak and later developed by Holman and Wiegert: a second-order polynomial calibrated to broadly match numerical simulations. However, the polynomial is unable to capture islands of instability introduced by mean motion resonances, causing the accuracy of the criterion to approach that of a random coin-toss when close to the boundary. We show how a deep neural network (DNN) trained on N-body simulations generated with REBOUND is able to significantly improve stability predictions for circumbinary planets on initially coplanar, circular orbits. Specifically, we find that the accuracy of our DNN never drops below 86%, even when tightly surrounding the boundary of instability, and is fast enough to be practical for on-the-fly calls during likelihood evaluations typical of modern Bayesian inference. Our binary classifier DNN is made publicly available at https://github.com/CoolWorlds/orbital-stability


title: Stable Weyl points, trivial surface states and particle-hole compensation in WP2

abstract: A possible connection between extremely large magneto-resistance and the presence of Weyl points has garnered much attention in the study of topological semimetals. Exploration of these concepts in transition metal phosphide WP2 has been complicated by conflicting experimental reports. Here we combine angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT) calculations to disentangle surface and bulk contributions to the ARPES intensity, the superposition of which has plagued the determination of the electronic structure in WP2. Our results show that while the hole- and electron-like Fermi surface sheets originating from surface states have different areas, the bulk-band structure of WP2 is electron-hole-compensated in agreement with DFT. Furthermore, the detailed band structure is compatible with the presence of at least 4 temperature-independent Weyl points, confirming the topological nature of WP2 and its stability against lattice distortions.


title: Combined search for electroweak production of charginos and neutralinos in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A statistical combination of several searches for the electroweak production of charginos and neutralinos is presented. All searches use proton-proton collision data at $\sqrt{s}=$ 13 TeV, recorded with the CMS detector at the LHC in 2016 and corresponding to an integrated luminosity of 35.9 fb$^{-1}$. In addition to the combination of previous searches, a targeted analysis requiring three or more charged leptons (electrons or muons) is presented, focusing on the challenging scenario in which the difference in mass between the two least massive neutralinos is approximately equal to the mass of the Z boson. The results are interpreted in simplified models of chargino-neutralino or neutralino pair production. For chargino-neutralino production, in the case when the lightest neutralino is massless, the combination yields an observed (expected) limit at the 95% confidence level on the chargino mass of up to 650 (570) GeV, improving upon the individual analysis limits by up to 40 GeV. If the mass difference between the two least massive neutralinos is approximately equal to the mass of the Z boson in the chargino-neutralino model, the targeted search requiring three or more leptons obtains observed and expected exclusion limits of around 225 GeV on the second neutralino mass and 125 GeV on the lightest neutralino mass, improving the observed limit by about 60 GeV in both masses compared to the previous CMS result. In the neutralino pair production model, the combined observed (expected) exclusion limit on the neutralino mass extends up to 650-750 (550-750) GeV, depending on the branching fraction assumed. This extends the observed exclusion achieved in the individual analyses by up to 200 GeV. The combined result additionally excludes some intermediate gaps in the mass coverage of the individual analyses.


title: Sheaves on the alcoves and modular representations II

abstract: We relate the category of sheaves on alcoves that was constructed in "Sheaves on the alcoves and modular representations I" to the representation theory of reductive algebraic groups. In particular, we show that its indecomposable projective objects encode the simple rational characters of a reductive algebraic group in all characteristics above the Coxeter number.


abstract: We consider the set of affine alcoves associated with a root system R as a topological space and consider a certain category S of sheaves of Z-modules on this space. Here Z is the structure algebra of the root system over a field k. To any wall reflection we associate a wall crossing functor on S. In the companion article "Sheaves on the alcoves and modular representations II" we prove that S encodes the simple rational characters of the connected, simply connected algebraic group with root system R over k, in the case that k is algebraically closed with characteristic above the Coxeter number.


title: The non-triviality of the vacuum in light-front quantization: An elementary treatment

abstract: It is often stated that the vacuum is trivial when light-front (null-plane) quantization is applied to a quantum field theory, in contrast to the situation with equal-time quantization. In fact, it is has long been known that the statement is false, and that in certain cases the standard rules for light-front perturbation theory need modification. This paper gives an elementary review of these issues, including an explanation of how and when there is a failure of the elementary derivation of the rules for light-front perturbation theory.


title: Harnack inequality for a class of Kolmogorov-Fokker-Planck equations in non-divergence form

abstract: We prove invariant Harnack inequalities for certain classes of non-divergence form equations of Kolmogorov type. The operators we consider exhibit invariance properties with respect to a homogeneous Lie group structure. The coefficient matrix is assumed either to satisfy a Cordes-Landis condition on the eigenvalues, or to admit a uniform modulus of continuity.