title: SINk: A Middleware for Synchronization of Heterogeneous Software Interfaces

abstract: Software is everywhere. The increasing requirement of supporting a wide variety of domains has rapidly increased the complexity of software systems, making them hard to maintain and the training process harder for end-users, which in turn ultimately demanded the development of user-friendly application software with simple interfaces that makes them easy, especially for non-professional personnel, to employ, and interact with. However, due to the lack of source code access for third-party software and the lack of non-graphical interfaces such as web-services or RMI (Remote Method Invocation) access to application functionality, synchronization between heterogeneous closed-box software interfaces and a user-friendly version of those interfaces has become a major challenge. In this paper, we design SINk, a middleware that enables synchronization of multiple heterogeneous software applications, using only graphical interface, without the need for source code access or access to the entire platform's control. SINk helps with synchronization of closed-box industry software, where in fact the only possible way of communication is through software interfaces. It leverages efficient client sever architecture, socket based protocol, adaptation to resolution changes, and parameter mapping mechanisms to transfer control events to ensure the real-time requirements of synchronization among multiple interfaces are met. Our proof-of-concept evaluation shows there is in fact potential usage of our middleware in a wide variety of domains.


title: Pseudoscalar-exchange contribution to $(g-2)_{\mu}$ from rational approximants

abstract: We present our recent results on the hadronic light-by-light pseudoscalar-exchange contribution to the anomalous magnetic moment of the muon using rational approximants. Our work provides a generalization of Pad\'e approximants to the bivariate case to describe the most general doubly-virtual transition form factor, which is required in the calculation. This method provides a powerful tool to systematically implement the known experimental data on transition form factors in the space-like region and low energies in a model-independent way. Given the lack of experimental data on the doubly-virtual transition form factor, we make use of the pseduoscalar decays into lepton pairs. We find an interesting and puzzling situation which calls for new experimental measurements to clarify the present state.


abstract: In this paper we introduce the notion of $\Delta$-weak character amenable Banach algebras and investigate $\Delta$-weak character amenability of certain Banach algebras such as projective tensor product $A\widehat{\otimes}B$, Lau product $A\times_TB$ for every Banach algebra $ A$ and $B$, where $T$ be a homomorphism from $B$ into $A$, abstract Segal algebras and module extension Banach algebras.


title: Toroidal dipole excitations in metamolecules formed by interacting plasmonic nanorods

abstract: We show how the elusive toroidal dipole moment appears as a radiative excitation eigenmode in a metamolecule resonator that is formed by pairs of plasmonic nanorods. We analyze one such nanorod configuration - a toroidal metamolecule. We find that the radiative interactions in the toroidal metamolecule can be qualitatively represented by a theoretical model based on an electric point dipole arrangement. Both a finite-size rod model and the point dipole approximation demonstrate how the toroidal dipole moment is subradiant and difficult to excite by incident light. By means of breaking the geometric symmetry of the metamolecule, the toroidal mode can be excited by linearly polarized light and appears as a Fano resonance dip in the forward scattered light. We provide simple optimization protocols for maximizing the toroidal dipole mode excitation. This opens up possibilities for simplified control and driving of metamaterial arrays consisting of toroidal dipole unit-cell resonators.


title: Stochastically Transitive Models for Pairwise Comparisons: Statistical and Computational Issues

abstract: There are various parametric models for analyzing pairwise comparison data, including the Bradley-Terry-Luce (BTL) and Thurstone models, but their reliance on strong parametric assumptions is limiting. In this work, we study a flexible model for pairwise comparisons, under which the probabilities of outcomes are required only to satisfy a natural form of stochastic transitivity. This class includes parametric models including the BTL and Thurstone models as special cases, but is considerably more general. We provide various examples of models in this broader stochastically transitive class for which classical parametric models provide poor fits. Despite this greater flexibility, we show that the matrix of probabilities can be estimated at the same rate as in standard parametric models. On the other hand, unlike in the BTL and Thurstone models, computing the minimax-optimal estimator in the stochastically transitive model is non-trivial, and we explore various computationally tractable alternatives. We show that a simple singular value thresholding algorithm is statistically consistent but does not achieve the minimax rate. We then propose and study algorithms that achieve the minimax rate over interesting sub-classes of the full stochastically transitive class. We complement our theoretical results with thorough numerical simulations.


title: Theoretical Uncertainties in Proton Lifetime Estimates

abstract: We recapitulate the primary sources of theoretical uncertainties in proton lifetime estimates in renormalizable, four-dimensional and non-supersymmetric grand unifications that represent the most conservative framework in which this question may be addressed at the perturbative level. We point out that many of these uncertainties are so severe that there are only very few scenarios in which an NLO approach, as crucial as it is for a real testability of any specific model, is actually sensible. Among these, the most promising seems to be the minimal renormalizable SO(10) GUT whose high-energy gauge symmetry is spontaneously broken by the adjoint and the five-index antisymmetric irreducible representations.


title: The mathematics of causal sets

abstract: The causal set approach to quantum gravity is based on the hypothesis that the underlying structure of spacetime is that of a random partial order. We survey some of the interesting mathematics that has arisen in connection with the causal set hypothesis, and describe how the mathematical theory can be translated to the application area. We highlight a number of open problems of interest to those working in causal set theory.


title: PERCH: Perception via Search for Multi-Object Recognition and Localization

abstract: In many robotic domains such as flexible automated manufacturing or personal assistance, a fundamental perception task is that of identifying and localizing objects whose 3D models are known. Canonical approaches to this problem include discriminative methods that find correspondences between feature descriptors computed over the model and observed data. While these methods have been employed successfully, they can be unreliable when the feature descriptors fail to capture variations in observed data; a classic cause being occlusion. As a step towards deliberative reasoning, we present PERCH: PErception via SeaRCH, an algorithm that seeks to find the best explanation of the observed sensor data by hypothesizing possible scenes in a generative fashion. Our contributions are: i) formulating the multi-object recognition and localization task as an optimization problem over the space of hypothesized scenes, ii) exploiting structure in the optimization to cast it as a combinatorial search problem on what we call the Monotone Scene Generation Tree, and iii) leveraging parallelization and recent advances in multi-heuristic search in making combinatorial search tractable. We prove that our system can guaranteedly produce the best explanation of the scene under the chosen cost function, and validate our claims on real world RGB-D test data. Our experimental results show that we can identify and localize objects under heavy occlusion--cases where state-of-the-art methods struggle.


title: Explicit Densities of Multidimensional L\'evy Walks

abstract: We provide explicit formulas for asymptotic densities of the 2- and 3-dimensional ballistic L\'evy walks. It turns out that in the 3D case the densities are given by elementary functions. The densities of the 2D L\'evy walks are expressed in terms of hypergeometric functions and the right-side Riemann-Liouville fractional derivative which allows to efficiently evaluate them numerically. The theoretical results agree with Monte-Carlo simulations. The obtained functions solve certain differential equations with the fractional material derivative.


title: Coherent Quantization using Coloured Surfaces

abstract: In this note, we revisit the quantization of Lie bialgebras described by the second author, placing it in the more general framework of the quantization of moduli spaces developed in our previous work. In particular, we show that embeddings of quilted surfaces (which are compatible with the choice of skeleton) induce morphisms between the corresponding quantized moduli spaces of flat connections. As an application, we describe quantizations of both the variety of Lagrangian subalgebras and the de-Cocini Procesi wonderful compactification, which are compatible with the action of the (quantized) Poisson Lie group.


title: A fast algorithm for simulating multiphase flows through periodic geometries of arbitrary shape

abstract: This paper presents a new boundary integral equation (BIE) method for simulating particulate and multiphase flows through periodic channels of arbitrary smooth shape in two dimensions. The authors consider a particular system---multiple vesicles suspended in a periodic channel of arbitrary shape---to describe the numerical method and test its performance. Rather than relying on the periodic Green's function as classical BIE methods do, the method combines the free-space Green's function with a small auxiliary basis, and imposes periodicity as an extra linear condition. As a result, we can exploit existing free-space solver libraries, quadratures, and fast algorithms, and handle a large number of vesicles in a geometrically complex channel. Spectral accuracy in space is achieved using the periodic trapezoid rule and product quadratures, while a first-order semi-implicit scheme evolves particles by treating the vesicle-channel interactions explicitly. New constraint-correction formulas are introduced that preserve reduced areas of vesicles, independent of the number of time steps taken. By using two types of fast algorithms, (i) the fast multipole method (FMM) for the computation of the vesicle-vesicle and the vesicle-channel hydrodynamic interaction, and (ii) a fast direct solver for the BIE on the fixed channel geometry, the computational cost is reduced to $O(N)$ per time step where $N$ is the spatial discretization size. Moreover, the direct solver inverts the wall BIE operator at $t = 0$, stores its compressed representation and applies it at every time step to evolve the vesicle positions, leading to dramatic cost savings compared to classical approaches. Numerical experiments illustrate that a simulation with $N=128, 000$ can be evolved in less than a minute per time step on a laptop.


title: Hierarchical gravitational fragmentation. I. Collapsing cores within collapsing clouds

abstract: We investigate the Hierarchical Gravitational Fragmentation scenario through numerical simulations of the prestellar stages of the collapse of a marginally gravitationally unstable isothermal sphere immersed in a strongly gravitationally unstable, uniform background medium. The core developes a Bonnor-Ebert (BE)-like density profile, while at the time of singularity (the protostar) formation the envelope approaches a singular-isothermal-sphere (SIS)-like $r^-2$ density profile. However, these structures are never hydrostatic. In this case, the central flat region is characterized by an infall speed, while the envelope is characterized by a uniform speed. This implies that the hydrostatic SIS initial condition leading to Shu's classical inside-out solution is not expected to occur, and therefore neither should the inside-out solution. Instead, the solution collapses from the outside-in, naturally explaining the observation of extended infall velocities. The core, defined by the radius at which it merges with the background, has a time-variable mass, and evolves along the locus of the ensemble of observed prestellar cores in a plot of $M/M_{BE}$ vs. $M$, where $M$ is the core's mass and $M_{BE}$ is the critical Bonnor-Ebert mass, spanning the range from the "stable" to the "unstable" regimes, even though it is collapsing at all times. We conclude that the presence of an unstable background allows a core to evolve dynamically from the time when it first appears, even when it resembles a pressure-confined, stable BE-sphere. The core can be thought of as a ram-pressure confined BE-sphere, with an increasing mass due to the accretion from the unstable background.


title: Integrated optical addressing of an ion qubit

abstract: The long coherence times and strong Coulomb interactions afforded by trapped ion qubits have enabled realizations of the necessary primitives for quantum information processing (QIP), and indeed the highest-fidelity quantum operations in any qubit to date. But while light delivery to each individual ion in a system is essential for general quantum manipulations and readout, experiments so far have employed optical systems cumbersome to scale to even a few tens of qubits. Here we demonstrate lithographically defined nanophotonic waveguide devices for light routing and ion addressing fully integrated within a surface-electrode ion trap chip. Ion qubits are addressed at multiple locations via focusing grating couplers emitting through openings in the trap electrodes to ions trapped 50 $\mu$m above the chip; using this light we perform quantum coherent operations on the optical qubit transition in individual $^{88}$Sr$^+$ ions. The grating focuses the beam to a diffraction-limited spot near the ion position with a 2 $\mu$m 1/$e^2$-radius along the trap axis, and we measure crosstalk errors between $10^{-2}$ and $4\times10^{-4}$ at distances 7.5-15 $\mu$m from the beam center. Owing to the scalability of the planar fabrication employed, together with the tight focusing and stable alignment afforded by optics integration within the trap chip, this approach presents a path to creating the optical systems required for large-scale trapped-ion QIP.


title: Minimal physical requirements for crystal growth self-poisoning

abstract: Self-poisoning is a kinetic trap that can impair or prevent crystal growth in a wide variety of physical settings. Here we use dynamic mean-field theory and computer simulation to argue that poisoning is ubiquitous because its emergence requires only the notion that a molecule can bind in two (or more) ways to a crystal; that those ways are not energetically equivalent; and that the associated binding events occur with sufficiently unequal probability. If these conditions are met then the steady-state growth rate is in general a non-monotonic function of the thermodynamic driving force for crystal growth, which is the characteristic of poisoning. Our results also indicate that relatively small changes of system parameters could be used to induce recovery from poisoning.


title: Mean-field limit of generalized Hawkes processes

abstract: We generalize multivariate Hawkes processes mainly by including a dependence with respect to the age of the process, i.e. the delay since the last point. Within this class, we investigate the limit behaviour, when n goes to infinity, of a system of n mean-field interacting age-dependent Hawkes processes. We prove that such a system can be approximated by independent and identically distributed age dependent point processes interacting with their own mean intensity. This result generalizes the study performed by Delattre, Fournier and Hoffmann (2015). In continuity with the work of Chevallier et al. (2015), the second goal of this paper is to give a proper link between these generalized Hawkes processes as microscopic models of individual neurons and the age-structured system of partial differential equations introduced by Pakdaman, Perthame and Salort (2010) as macroscopic model of neurons.


title: A classification of torsors over Laurent polynomial rings

abstract: Let R\_n be the ring of Laurent polynomials in n variables over a field k of characteristic zero and let K\_n be its fraction field.Given a linear algebraic k-group $G$, we show that a K\_n-torsor under G which is unramified with respect to X=Spec(R\_n)extends to a unique toral R\_n-torsor under G. This result, in turn, allows us to classify all G-torsors over R\_n.


title: Sparse Polynomial Systems with many Positive Solutions from Bipartite Simplicial Complexes

abstract: Consider a regular triangulation of the convex-hull $P$ of a set $\mathcal A$ of $n$ points in $\mathbb R^d$, and a real matrix $C$ of size $d \times n$. A version of Viro's method allows to construct from these data an unmixed polynomial system with support $\mathcal A$ and coefficient matrix $C$ whose number of positive solutions is bounded from below by the number of $d$-simplices which are positively decorated by $C$. We show that all the $d$-simplices of a triangulation can be positively decorated if and only if the triangulation is balanced, which in turn is equivalent to the fact that its dual graph is bipartite. This allows us to identify, among classical families, monomial supports which admit maximally positive systems, i.e. systems all toric complex solutions of which are real and positive. These families give some evidence in favor of a conjecture due to Bihan. We also use this technique in order to construct fewnomial systems with many positive solutions. This is done by considering a simplicial complex with bipartite dual graph included in a regular triangulation of the cyclic polytope.


title: On the Bound States of Schr\"odinger Operators with $\delta$-interactions on Conical Surfaces

abstract: In dimension greater than or equal to three, we investigate the spectrum of a Schr{\"o}dinger operator with a $\delta$-interaction supported on a cone whose cross section is the sphere of co-dimension two. After decomposing into fibers, we prove that there is discrete spectrum only in dimension three and that it is generated by the axisymmetric fiber. We get that these eigenvalues are non-decreasing functions of the aperture of the cone and we exhibit the precise logarithmic accumulation of the discrete spectrum below the threshold of the essential spectrum.


title: Quantum dissipative dynamics of a bistable system in the sub-Ohmic to super-Ohmic regime

abstract: We investigate the quantum dynamics of a multilevel bistable system coupled to a bosonic heat bath beyond the perturbative regime. We consider different spectral densities of the bath, in the transition from sub-Ohmic to super-Ohmic dissipation, and different cutoff frequencies. The study is carried out by using the real-time path integral approach of the Feynman-Vernon influence functional. We find that, in the crossover dynamical regime characterized by damped \emph{intrawell} oscillations and incoherent tunneling, the short time behavior and the time scales of the relaxation starting from a nonequilibrium initial condition depend nontrivially on the spectral properties of the heat bath.


title: A Topologically Driven Glass in Ring Polymers

abstract: The static and dynamic properties of ring polymers in concentrated solutions remains one of the last deep unsolved questions in Polymer Physics. At the same time, the nature of the glass transition in polymeric systems is also not well understood. In this work we study a novel glass transition in systems made of circular polymers by exploiting the topological constraints that are conjectured to populate concentrated solutions of rings. We show that such rings strongly inter-penetrate through one another, generating an extensive network of topological interactions that dramatically affects their dynamics. We show that a kinetically arrested state can be induced by randomly pinning a small fraction of the rings. This occurs well above the classical glass transition temperature at which microscopic mobility is lost. Our work demonstrates both the existence of long-lived inter-ring penetrations and also realises a novel, topologically-induced, glass transition.


title: Two-loop Integrand Decomposition into Master Integrals and Surface Terms

abstract: Loop amplitudes are conveniently expressed in terms of master integrals whose coefficients carry the process dependent information. Similarly before integration, the loop integrands may be expressed as a linear combination of propagator products with universal numerator-tensors. Such a decomposition is an important input for the numerical unitarity approach, which constructs integrand coefficients from on-shell tree amplitudes. We present a new method to organise multi-loop integrands into a direct sum of terms that integrate to zero (surface terms) and remaining master integrands. This decomposition facilitates a general, numerical unitarity approach for multi-loop amplitudes circumventing analytic integral reduction. Vanishing integrals are well known as integration-by-parts identities. Our construction can be viewed as an explicit construction of a complete set of integration-by-parts identities excluding doubled propagators. Interestingly, a class of 'horizontal' identities is singled out which hold as well for altered propagator powers.


title: Counting Restricted Dyck Paths Through Random Walks

abstract: We show connection between Dyck paths with peaks of bounded height and random walks. The correspondence between a certain class of random walks and such Dyck paths allows us to develop a probabilistic perspective on Chebyshev polynomials.


title: Interacting binaries W Serpentids and Double Periodic Variables

abstract: W Serpentids and Double Periodic Variables (DPVs) are candidates for close interacting binaries in a non-conservative evolutionary stage; while W Serpentids are defined by high-excitation ultraviolet emission lines present during most orbital phases, and by usually showing variable orbital periods, DPVs are characterized by a long photometric cycle lasting roughly 33 times the (practically constant) orbital period. We report the discovery of 7 new Galactic DPVs, increasing the number of known DPVs in our Galaxy by 50%. We find that DPVs are tangential-impact systems, i.e. their primaries have radii barely larger than the critical Lubow-Shu radius. These systems are expected to show transient discs, but we find that they host stable discs with radii smaller than the tidal radius. Among tangential-impact systems including DPVs and semi-detached Algols, only DPVs have primaries with masses between 7 and 10 $M_{\odot}$. We find that DPVs are in a Case-B mass transfer stage with donor masses between 1 and 2 M$_{\odot}$ and with primaries resembling Be stars. W Serpentids are impact and non-impact systems, their discs extend until the last non-intersecting orbit and show a larger range of stellar mass and mass ratio than DPVs. Infrared photometry reveals significant color excesses in many DPVs and W Serpentids, usually larger for the latter ones, suggesting variable amounts of circumstellar matter.


title: Ultra high energy cosmic rays: the highest energy frontier

abstract: Ultra-high energy cosmic rays (UHECRs) are the highest energy messengers of the present universe, with energies up to $10^{20}$ eV. Studies of astrophysical particles (nuclei, electrons, neutrinos and photons) at their highest observed energies have implications for fundamental physics as well as astrophysics. The primary particles interact in the atmosphere and generate extensive air showers. Analysis of those showers enables one not only to estimate the energy, direction and most probable mass of the primary cosmic particles, but also to obtain information about the properties of their hadronic interactions at an energy more than one order of magnitude above that accessible with the current highest energy human-made accelerator. In this contribution we will review the state-of-the-art in UHECRs detection. We will present the leading experiments Pierre Auger Observatory and Telescope Array and discuss the cosmic ray energy spectrum, searches for directional anisotropy, studies of mass composition, the determination of the number of shower muons (which is sensitive to the shower hadronic interactions) and the proton-air cross section.


title: Model-Independent Dark Energy Equation of State from Unanchored Baryon Acoustic Oscillations

abstract: Ratios of line of sight baryon acoustic oscillation (BAO) peaks at two redshifts only depend upon the average dark energy equation of states between those redshifts, as the dependence on anchors such as the BAO scale or the Hubble constant is canceled in a ratio. As a result, BAO ratios provide a probe of dark energy which is independent of both the cosmic distance ladder and the early evolution of universe. In this note, we use ratios to demonstrate that the known tension between the Lyman alpha forest BAO measurement and other probes arises entirely from recent (0.57<z<2.34) cosmological expansion. Using ratios of the line of sight Lyman alpha forest and BOSS CMASS BAO scales, we show that there is already more than 3 sigma tension with the standard LambdaCDM cosmological model which implies that either (i) The BOSS Lyman alpha forest measurement of the Hubble parameter was too low as a result of a statistical fluctuation or systematic error or else (ii) the dark energy equation of state falls steeply at high redshift.


title: The site frequency spectrum for general coalescents

abstract: General genealogical processes such as $\Lambda$- and $\Xi$-coalescents, which respectively model multiple and simultaneous mergers, have important applications in studying marine species, strong positive selection, recurrent selective sweeps, strong bottlenecks, large sample sizes, and so on. Recently, there has been significant progress in developing useful inference tools for such general models. In particular, inference methods based on the site frequency spectrum (SFS) have received noticeable attention. Here, we derive a new formula for the expected SFS for general $\Lambda$- and $\Xi$-coalescents, which leads to an efficient algorithm. For time-homogeneous coalescents, the runtime of our algorithm for computing the expected SFS is $O(n^2)$, where $n$ is the sample size. This is a factor of $n^2$ faster than the state-of-the-art method. Furthermore, in contrast to existing methods, our method generalizes to time-inhomogeneous $\Lambda$- and $\Xi$-coalescents with measures that factorize as $\Lambda(dx)/\zeta(t)$ and $\Xi(dx)/\zeta(t)$, respectively, where $\zeta$ denotes a strictly positive function of time. The runtime of our algorithm in this setting is $O(n^3)$. We also obtain general theoretical results for the identifiability of the $\Lambda$ measure when $\zeta$ is a constant function, as well as for the identifiability of the function $\zeta$ under a fixed $\Xi$ measure.


title: Decomposability of Finitely Generated Torsion-free Nilpotent Groups

abstract: We describe an algorithm for deciding whether or not a given finitely generated torsion-free nilpotent group is decomposable as the direct product of nontrivial subgroups.


title: Global Antineutrino Modeling for a Web Application

abstract: Antineutrinos stream freely from rapidly decaying fission products within the cores of nuclear reactors and from long-lived natural radioactivity within the rocky layers of the Earth. These global antineutrinos produce detectable signals in large ultra-clear volumes of water- or hydrocarbon-based target liquids, which are viewed by inward-facing photomultiplier tubes. Detected antineutrinos provide information about their shrouded sources and about the fundamental properties of neutrinos themselves. This paper presents the input data, formulae, and plots resulting from the calculations, which, in addition to the time-dependent reaction rates and energy spectra, model the directions of the antineutrinos from IAEA-registered nuclear power reactors and of the neutrinos from $^8$B decay in the Sun. The model includes estimates of the steady state reaction rates and energy spectra of the antineutrinos from the crust and mantle of the Earth. Results are available for any location near the surface of the Earth and comprise both quasi-elastic scattering on free protons and elastic scattering on atomic electrons. This paper compares model results for two underground locations, the Boulby Mine in the United Kingdom and the Morton Salt Mine in the United States. Operational nuclear power reactors are within about $20$ kilometers of these mines, making them candidate sites for antineutrino detectors capable of identifying, monitoring, and locating remote nuclear activity. The model, which is implemented in a web application at https://geoneutrinos.org/reactors/, provides references for the input data and the formulae, as well as an interactive calculator of the significance of the rate of any of the neutrino sources relative to other sources taken as background.


title: Optimal Slater-determinant approximation of fermionic wave functions

abstract: We study the optimal Slater-determinant approximation of an $N$-fermion wave function analytically. That is, we seek the Slater-determinant (constructed out of $N$ orthonormal single-particle orbitals) wave function having largest overlap with a given $N$-fermion wave function. Some simple lemmas have been established and their usefulness is demonstrated on some structured states, such as the Greenberger-Horne-Zeilinger state. In the simplest nontrivial case of three fermions in six orbitals, which the celebrated Borland-Dennis discovery is about, the optimal Slater approximation wave function is proven to be built out of the natural orbitals in an interesting way. We also show that the Hadamard inequality is useful for finding the optimal Slater approximation of some special target wave functions.


title: Nonlinear Evolution of a Baroclinic Wave and Imbalanced Dissipation

abstract: We consider the nonlinear evolution of an unstable baroclinic wave in a regime of rotating stratified flow that is of relevance to interior circulation in the oceans and in the atmosphere---a regime characterized by small large-scale Rossby and Froude numbers, a small vertical to horizontal aspect ratio, and no bounding horizontal surfaces. Using high-resolution simulations of the non-hydrostatic Boussinesq equations and companion integrations of the balanced quasi-geostrophic equations, we present evidence for a local route to dissipation of balanced energy directly through interior turbulent cascades. Analysis of simulations presented in this study suggest that a developing baroclinic instability can lead to secondary instabilities that can cascade a small fraction of the energy forward to unbalanced scales. Mesoscale shear and strain resulting from the hydrostatic geostrophic baroclinic instability drive frontogenesis. The fronts in turn support ageostrophic secondary circulation and instabilities. These two processes acting together lead to a quick rise in dissipation rate which then reaches a peak and begins to fall slowly when frontogenesis slows down; eventually balanced and imbalanced modes decouple. Dissipation of balanced energy by imbalanced processes scales exponentially with Rossby number of the base flow. We expect that this scaling will hold more generally than for the specific setup we consider given the fundamental nature of the dynamics involved. A break is seen in the total energy spectrum at small scales: While a steep $k^{-3}$ geostrophic scaling (where $k$ is the three-dimensional wavenumber) is seen at intermediate scales, the smaller scales display a shallower $k^{-5/3}$ scaling, reminiscent of the atmospheric spectra of Nastrom & Gage. At higher Ro, the vertical shear spectrum has a minimum, like in some relevant obervations


title: From the weak Bruhat order to crystal posets

abstract: We investigate the ways in which fundamental properties of the weak Bruhat order on a Weyl group can be lifted (or not) to a corresponding highest weight crystal graph, viewed as a partially ordered set; the latter projects to the weak order via the key map. First, a crystal theoretic analogue of the statement that any two reduced expressions for the same Coxeter group element are related by Coxeter moves is proven for all lower intervals in a simply or doubly laced crystal. On the other hand, it is shown that no finite set of moves exists, even in type A, for arbitrary crystal graph intervals. In fact, it is shown that there are relations of arbitrarily high degree amongst crystal operators that are not implied by lower degree relations. Second, for crystals associated to Kac-Moody algebras it is shown for lower intervals that the Mobius function is always 0, 1, or -1, and in finite type this is also proven for upper intervals, with a precise formula given in each case. Moreover, the order complex for each of these intervals is proven to be homotopy equivalent to a ball or to a sphere of some dimension, despite often not being shellable. For general intervals, examples are constructed with arbitrarily large Mobius function, again even in type A. Any interval having Mobius function other than 0, 1, or -1 is shown to contain within it a relation amongst crystal operators that is not implied by the relations giving rise to the local structure of the crystal, making precise a tight relationship between the Mobius function and these somewhat unexpected relations appearing in crystals. New properties of the key map are also derived. The key is shown to be determined entirely by the edge-colored poset-theoretic structure of the crystal, and a recursive algorithm is given for calculating it.


title: Observing Near-Earth Objects with the James Webb Space Telescope

abstract: The James Webb Space Telescope (JWST) has the potential to enhance our understanding of near-Earth objects (NEOs). We present results of investigations into the observability of NEOs given the nominal observing requirements of JWST on elongation (85-135 degrees) and non-sidereal rates ($<$30mas/s). We find that approximately 75% of NEOs can be observed in a given year. However, observers will need to wait for appropriate observing windows. We find that JWST can easily execute photometric observations of meter-sized NEOs which will enhance our understanding of the small NEO population.


title: Explicit upper bounds for the spectral distance of two trace class operators

abstract: Given two trace class operators A and B on a separable Hilbert space we provide an upper bound for the Hausdorff distance of their spectra involving only the distance of A and B in operator norm and the singular values of A and B. By specifying particular asymptotics of the singular values our bound reproduces or improves existing bounds for the spectral distance. The proof is based on lower and upper bounds for determinants of trace class operators of independent interest.


title: Melting the core of giant planets: impact on tidal dissipation

abstract: Giant planets are believed to host central dense rocky/icy cores that are key actors in the core-accretion scenario for their formation. In the same time, some of their components are unstable in the temperature and pressure regimes of central regions of giant planets and only ab-initio EOS computations can address the question of the state of matter. In this framework, several works demonstrated that erosion and redistribution of core materials in the envelope must be taken into account. These complex mechanisms thus may deeply modify giant planet interiors for which signatures of strong tidal dissipation have been obtained for Jupiter and Saturn. The best candidates to explain this dissipation are the viscoelastic dissipation in the central dense core and turbulent friction acting on tidal inertial waves in their fluid convective envelope. In this work, we study the consequences of the possible melting of central regions for the efficiency of each of these mechanisms.


title: The impact of model detail on power grid resilience measures

abstract: Extreme events represent a challenge to natural as well as man-made systems. For critical infrastructure like power grids, we need to understand their resilience against large disturbances. Recently, new measures of the resilience of dynamical systems have been developed in the complex system literature. Basin stability and survivability respectively assess the asymptotic and transient behavior of a system when subjected to arbitrary, localized but large perturbations. To employ these methods to assess the resilience of power grids, we need to choose a model of the power grid. So far the most popular model that has been studied is the classical swing equation model for the frequency response of generators and motors. In this paper we study a more sophisticated model of synchronous machines that also takes voltage dynamics into account, and compare it to the previously studied model. This model has been found to give an accurate picture of the long term evolution of synchronous machines in the engineering literature for post fault studies. We find evidence that some stable fix points of the swing equation become unstable when we add voltage dynamics. If this occurs the asymptotic behavior of the system can be dramatically altered, and basin stability estimates obtained with the swing equation can be dramatically wrong. We also find that the survivability does not change significantly when taking the voltage dynamics into account. Further, the limit cycle type asymptotic behaviour is strongly correlated with transient voltages that violate typical operational voltage bounds. Thus, transient voltage bounds are dominated by transient frequency bounds and play no large role for realistic parameters.


title: Complex molecules in W51 North region

abstract: We present Submillimeter Array (SMA) molecular line observations in two 2 GHz-wide bands centered at 217.5 and 227.5 GHz, toward the massive star forming region W51 North. We identified 84 molecular line transitions from 17 species and their isotopologues. The molecular gas distribution of these lines mainly peaks in the continuum position of W51 North, and has a small tail extending to the west, probably associated with W51 d2. In addition to the commonly detected nitrogen and oxygen-bearing species, we detected a large amount of transitions of the Acetone (CH$_3$COCH$_3$) and Methyl Formate (CH$_3$OCHO), which may suggest that these molecules are present in an early evolutionary stage of the massive stars. We also found that W51 North is an ethanol-rich source. There is no obvious difference in the molecular gas distributions between the oxygen-bearing and nitrogen-bearing molecules. Under the assumption of Local Thermodynamic Equilibrium (LTE), with the XCLASS tool, the molecular column densities, and rotation temperatures are estimated. We found that the oxygen-bearing molecules have considerable higher column densities and fractional abundances than the nitrogen-bearing molecules. The rotation temperatures range from 100 to 200 K, suggesting that the molecular emission could be originated from a warm environment. Finally, based on the gas distributions, fractional abundances and the rotation temperatures, we conclude that CH$_3$OH, C$_2$H$_5$OH, CH$_3$COCH$_3$ and CH$_3$CH$_2$CN might be synthesized on the grain surface, while gas phase chemistry is responsible for the production of CH$_3$OCH$_3$, CH$_3$OCHO and CH$_2$CHCN.


title: Four and one more: The formation history and total mass of globular clusters in the Fornax dSph

abstract: We have determined the detailed star formation history and total mass of the globular clusters in the Fornax dwarf spheroidal using archival HST WFPC2 data. Colour magnitude diagrams are constructed in the F555W and F814W bands and corrected for the effect of Fornax field star contamination, after which we use the routine Talos to derive the quantitative star formation history as a function of age and metallicity. The star formation history of the Fornax globular clusters shows that Fornax 1, 2, 3 and 5 are all dominated by ancient (>10 Gyr) populations. Cluster Fornax 1,2 and 3 display metallicities as low as [Fe/H]=-2.5 while Fornax 5 is slightly more metal-rich at [Fe/H]=-1.8, consistent with resolved and unresolved metallicity tracers. Conversely, Fornax 4 is dominated by a more metal-rich~([Fe/H]=-1.2) and younger population at 10 Gyr, inconsistent with the other clusters. A lack of stellar populations overlapping with the main body of Fornax argues against the nucleus cluster scenario for Fornax 4. The combined stellar mass in globular clusters as derived from the SFH is (9.57$\pm$0.93)$\times$10$^{5}$ M$_{\odot}$ which corresponds to 2.5$\pm$0.2 percent of the total stellar mass in Fornax. The mass of the four most metal-poor clusters can be further compared to the metal-poor Fornax field to yield a mass fraction of 19.6$\pm$3.1 percent. Therefore, the SFH results provide separate supporting evidence for the unusually high mass fraction of the GCs compared to the Fornax field population.


title: Reply to "Comment on `Axion induced oscillating electric dipole moments' " [1]

abstract: We respond to a paper of Flambaum, et.al. [Phys. Rev. D95, no. 5, 058701 (2017)], claiming there is no effective induced oscillating electric dipole moment, e.g., for the electron, arising from interaction with an oscillating cosmic axion background via the anomaly. The relevant Feynman amplitude, Fig.(1), as computed by Flambaum et.al., becomes a total divergence, and vanishes. Contrary to this result, we obtained a nonvanishing amplitude, that yields physical electric dipole radiation for an electron (or any magnetic dipole moment) immersed in a cosmic axion field. We argue that the Flambaum et.al. counter-claim is incorrect, and is based upon a misunderstanding of a physics choice vs. gauge choice, and an assumption that electric dipoles be defined only by coupling to static (constant in time) electric fields.


title: Star formation and molecular hydrogen in dwarf galaxies: a non-equilibrium view

abstract: We study the connection of star formation to atomic (HI) and molecular hydrogen (H$_2$) in isolated, low metallicity dwarf galaxies with high-resolution ($m_{\rm gas}$ = 4 M$_\odot$, $N_{\rm ngb}$ = 100) SPH simulations. The model includes self-gravity, non-equilibrium cooling, shielding from an interstellar radiation field, the chemistry of H$_2$ formation, H$_2$-independent star formation, supernova feedback and metal enrichment. We find that the H$_2$ mass fraction is sensitive to the adopted dust-to-gas ratio and the strength of the interstellar radiation field, while the star formation rate is not. Star formation is regulated by stellar feedback, keeping the gas out of thermal equilibrium for densities $n <$ 1 cm$^{-3}$. Because of the long chemical timescales, the H$_2$ mass remains out of chemical equilibrium throughout the simulation. Star formation is well-correlated with cold ( T $\leqslant$ 100 K ) gas, but this dense and cold gas - the reservoir for star formation - is dominated by HI, not H$_2$. In addition, a significant fraction of H$_2$ resides in a diffuse, warm phase, which is not star-forming. The ISM is dominated by warm gas (100 K $<$ T $\leqslant 3\times 10^4$ K) both in mass and in volume. The scale height of the gaseous disc increases with radius while the cold gas is always confined to a thin layer in the mid-plane. The cold gas fraction is regulated by feedback at small radii and by the assumed radiation field at large radii. The decreasing cold gas fractions result in a rapid increase in depletion time (up to 100 Gyrs) for total gas surface densities $\Sigma_{\rm HI+H_2} \lesssim$ 10 M$_\odot$pc$^{-2}$, in agreement with observations of dwarf galaxies in the Kennicutt-Schmidt plane.


title: Size evolution of normal and compact galaxies in the EAGLE simulation

abstract: We present the evolution of galaxy sizes, from redshift 2 to 0, for actively star forming and passive galaxies in the cosmological hydrodynamical 1003 cMpc3 simulation of the EAGLE project. We find that the sizes increase with stellar mass , but that the relation weakens with increasing redshift. Separating galaxies by their star formation activity, we find that passive galaxies are typically smaller than active galaxies at fixed stellar mass. These trends are consistent with those found in observations and the level of agreement between the predicted and observed size - mass relation is of order 0.1 dex for z < 1 and 0.2-0.3 dex from redshift 1 to 2. We use the simulation to compare the evolution of individual galaxies to that of the population as a whole. While the evolution of the size-stellar mass relation for active galaxies provides a good proxy for the evolution of individual galaxies, the evolution of individual passive galaxies is not well represented by the observed size - mass relation due to the evolving number density of passive galaxies. Observations of z \approx 2 galaxies have revealed an abundance of massive red compact galaxies, that depletes below z \approx 1. We find that a similar population forms naturally in the simulation. Comparing these galaxies to their z = 0 descendants, we find that all compact galaxies grow in size due to the high-redshift stars migrating outwards. Approximately 60% of the compact galaxies increase in size further due to renewed star formation and/or mergers.


title: Sterile neutrino Dark Matter production from scalar decay in a thermal bath

abstract: We calculate the production rate of singlet fermions from the decay of neutral or charged scalar fields in a hot plasma. We find that there are considerable thermal corrections when the temperature of the plasma exceeds the mass of the decaying scalar. We give analytic expressions for the temperature-corrected production rates in the regime where the decay products are relativistic. We also study the regime of non-relativistic decay products numerically. Our results can be used to determine the abundance and momentum distribution of Dark Matter particles produced in scalar decays. The inclusion of thermal corrections helps to improve predictions for the free streaming of the Dark Matter particles, which is crucial to test the compatibility of a given model with cosmic structure formation. With some modifications, our results may be generalised to the production of other Dark Matter candidates in scalar decays.


title: Lyman-Alpha Emitting Galaxies as a Probe of Reionization: Large-Scale Bubble Morphology and Small-Scale Absorbers

abstract: The visibility of LyA emitting galaxies during the Epoch of Reionization is controlled by both diffuse HI patches in large-scale bubble morphology and small-scale absorbers. To investigate the impact on LyA photons, we apply a novel combination of analytic and numerical calculations to three scenarios: (i) the `bubble' model, where only diffuse HI outside ionized bubbles is present; (ii) the `web' model, where HI exists only in overdense self-shielded gas; and (iii) the more realistic 'web-bubble' model, which contains both. Our analysis confirms that there is a degeneracy between the ionization structure of the intergalactic medium (IGM) and the HI fraction inferred from LyA surveys, as the three models suppress LyA flux equally with very different HI fractions. We argue that a joint analysis of the LyA luminosity function and the rest-frame equivalent width distribution/LyA fraction can break this degeneracy and provide constraints on the reionization history and its topology. We further show that constraints can improve if we consider the full shape of the M_UV-dependent redshift evolution of the LyA fraction of Lyman break galaxies. Contrary to conventional wisdom, we find that (i) a drop of LyA fraction larger for UV-faint than for UV-bright galaxies can be reproduced with web and web-bubble models and therefore does not provide exclusive evidence of patchy reionization, and (ii) the IGM-transmission PDF is unimodal for bubble models and bimodal in web models. We further highlight the importance of galaxy-absorber cross-correlation. Comparing our models to observations, the neutral fraction at z~7 is likely to be of order of tens of per cent when interpreted with bubble or web-bubble models. Alternatively, we obtain a conservative lower limit ~1% in the web models, if we allow for a drop in the photoionization rate by a factor of ~100 from the post-reionized universe. [abridged]


title: Space Telescope and Optical Reverberation Mapping Project. III. Optical Continuum Emission and Broad-Band Time Delays in NGC 5548

abstract: We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission.


title: Quasi-bound states of massive scalar fields in the Kerr black-hole spacetime: Beyond the hydrogenic approximation

abstract: Rotating black holes can support quasi-stationary (unstable) bound-state resonances of massive scalar fields in their exterior regions. These spatially regular scalar configurations are characterized by instability timescales which are much longer than the timescale $M$ set by the geometric size (mass) of the central black hole. It is well-known that, in the small-mass limit $\alpha\equiv M\mu\ll1$ (here $\mu$ is the mass of the scalar field), these quasi-stationary scalar resonances are characterized by the familiar hydrogenic oscillation spectrum: $\omega_{\text{R}}/\mu=1-\alpha^2/2{\bar n}^2_0$, where the integer $\bar n_0(l,n;\alpha\to0)=l+n+1$ is the principal quantum number of the bound-state resonance (here the integers $l=1,2,3,...$ and $n=0,1,2,...$ are the spheroidal harmonic index and the resonance parameter of the field mode, respectively). As it depends only on the principal resonance parameter $\bar n_0$, this small-mass ($\alpha\ll1$) hydrogenic spectrum is obviously degenerate. In this paper we go beyond the small-mass approximation and analyze the quasi-stationary bound-state resonances of massive scalar fields in rapidly-spinning Kerr black-hole spacetimes in the regime $\alpha=O(1)$. In particular, we derive the non-hydrogenic (and, in general, non-degenerate) resonance oscillation spectrum ${{\omega_{\text{R}}}/{\mu}}=\sqrt{1-(\alpha/{\bar n})^2}$, where $\bar n(l,n;\alpha)=\sqrt{(l+1/2)^2-2m\alpha+2\alpha^2}+1/2+n$ is the generalized principal quantum number of the quasi-stationary resonances. This analytically derived formula for the characteristic oscillation frequencies of the composed black-hole-massive-scalar-field system is shown to agree with direct numerical computations of the quasi-stationary bound-state resonances.


title: How stellar feedback simultaneously regulates star formation and drives outflows

abstract: We present an analytic model for how momentum deposition from stellar feedback simultaneously regulates star formation and drives outflows in a turbulent interstellar medium (ISM). Because the ISM is turbulent, a given patch of ISM exhibits sub-patches with a range of surface densities. The high-density patches are 'pushed' by feedback, thereby driving turbulence and self-regulating local star formation. Sufficiently low-density patches, however, are accelerated to above the escape velocity before the region can self-adjust and are thus vented as outflows. In the turbulent-pressure-supported regime, when the gas fraction is $\gtrsim 0.3$, the ratio of the turbulent velocity dispersion to the circular velocity is sufficiently high that at any given time, of order half of the ISM has surface density less than the critical value and thus can be blown out on a dynamical time. The resulting outflows have a mass-loading factor ($\eta \equiv M_{\rm out}/M_{\star}$) that is inversely proportional to the gas fraction times the circular velocity. At low gas fractions, the star formation rate needed for local self-regulation, and corresponding turbulent Mach number, decline rapidly; the ISM is 'smoother', and it is actually more difficult to drive winds with large mass-loading factors. Crucially, our model predicts that stellar-feedback-driven outflows should be suppressed at $z \lesssim 1$ in $M_{\star} \gtrsim 10^{10} M_{\odot}$ galaxies. This mechanism allows massive galaxies to exhibit violent outflows at high redshifts and then 'shut down' those outflows at late times, thereby enabling the formation of a smooth, extended thin stellar disk. We provide simple fitting functions for $\eta$ that should be useful for sub-resolution and semi-analytic models. [abridged]


title: The Concentration Dependence of the Galaxy-Halo Connection: Modeling Assembly Bias with Abundance Matching

abstract: Empirical methods for connecting galaxies to their dark matter halos have become essential for interpreting measurements of the spatial statistics of galaxies. In this work, we present a novel approach for parameterizing the degree of concentration dependence in the abundance matching method. This new parameterization provides a smooth interpolation between two commonly used matching proxies: the peak halo mass and the peak halo maximal circular velocity. This parameterization controls the amount of dependence of galaxy luminosity on halo concentration at a fixed halo mass. Effectively this interpolation scheme enables abundance matching models to have adjustable assembly bias in the resulting galaxy catalogs. With the new 400 Mpc/h DarkSky Simulation, whose larger volume provides lower sample variance, we further show that low-redshift two-point clustering and satellite fraction measurements from SDSS can already provide a joint constraint on this concentration dependence and the scatter within the abundance matching framework.


title: Model-independent determination of the axial mass parameter in quasielastic antineutrino-nucleon scattering

abstract: Understanding the charged current quasielestic (CCQE) neutrino-nucleus interaction is important for precision studies of neutrino oscillations. The theoretical description of the interaction depends on the combination of a nuclear model with the knowledge of form factors. While the former has received considerable attention, the latter, in particular the axial form factor, is implemented using the historical dipole model. Instead, we use a model-independent approach, presented in a previous study, to analyze the muon antineutrino CCQE mineral oil data published by the MiniBooNE collaboration. We combine the cross section for scattering of antineutrinos off protons in carbon and hydrogen, using the same axial form factor for both. The extracted value of the axial mass parameter $m_A = 0.84^{+0.12}_{-0.04} \pm {0.11} \, {\rm GeV}$ is in very good agreement with the model-independent value extracted from MiniBooNE's neutrino data. Going beyond a one-parameter description of the axial form factor, we extract values of the axial form factor in the range of $Q^2=0.1...1.0$ GeV$^2$, finding a very good agreement with the analogous extraction from the neutrino data. We discuss the implications of these results.


title: Comparing Experiments to the Fault-Tolerance Threshold

abstract: Achieving error rates that meet or exceed the fault-tolerance threshold is a central goal for quantum computing experiments, and measuring these error rates using randomized benchmarking is now routine. However, direct comparison between measured error rates and thresholds is complicated by the fact that benchmarking estimates average error rates while thresholds reflect worst-case behavior when a gate is used as part of a large computation. These two measures of error can differ by orders of magnitude in the regime of interest. Here we facilitate comparison between the experimentally accessible average error rates and the worst-case quantities that arise in current threshold theorems by deriving relations between the two for a variety of physical noise sources. Our results indicate that it is coherent errors that lead to an enormous mismatch between average and worst case, and we quantify how well these errors must be controlled to ensure fair comparison between average error probabilities and fault-tolerance thresholds.


title: On the universality of interstellar filaments: theory meets simulations and observations

abstract: Filaments are ubiquitous in the universe. Recent observations have revealed that stars and star clusters form preferentially along dense filaments. Understanding the formation and properties of filaments is therefore a crucial step in understanding star formation. Here we perform three-dimensional high-resolution magnetohydrodynamical simulations that follow the evolution of molecular clouds and the formation of filaments and stars. We apply a filament detection algorithm and compare simulations with different combinations of physical ingredients: gravity, turbulence, magnetic fields and jet/outflow feedback. We find that gravity-only simulations produce significantly narrower filament profiles than observed, while simulations that include turbulence produce realistic filament properties. For these turbulence simulations, we find a remarkably universal filament width of 0.10 +/- 0.02 pc, which is independent of the star formation history of the clouds. We derive a theoretical model that provides a physical explanation for this characteristic filament width, based on the sonic scale (lambda_sonic) of molecular cloud turbulence. Our derivation provides lambda_sonic as a function of the cloud diameter L, the velocity dispersion sigma_v, the gas sound speed c_s, and the ratio of thermal to magnetic pressure, plasma beta. For typical cloud conditions in the Milky Way spiral arms, we find lambda_sonic = 0.04-0.16 pc, in excellent agreement with the filament width of 0.05-0.15 pc from observations. Consistent with the theoretical model assumptions, we find that the velocity dispersion inside the filaments is subsonic and supersonic outside. We further explain the observed p=2 scaling of the filament density profile, rho ~ r^(-p) with the collision of two planar shocks forming a filament at their intersection.

categories: astro-ph.SR astro-ph.GA astro-ph.IM physics.data-an physics.flu-dyn


title: Characterization of decohering quantum systems: Machine learning approach

abstract: Adaptive data collection and analysis, where data are being fed back to update the measurement settings, can greatly increase speed, precision, and reliability of the characterization of quantum systems. However, decoherence tends to make adaptive characterization difficult. As an example, we consider two coupled discrete quantum systems. When one of the systems can be controlled and measured, the standard method to characterize another, with an unknown frequency $\omega_{\rm r}$, is swap spectroscopy. Here, adapting measurements can provide estimates whose error decreases exponentially in the number of measurement shots rather than as a power law in conventional swap spectroscopy. However, when the decoherence time is so short that an excitation oscillating between the two systems can only undergo less than a few tens of vacuum Rabi oscillations, this approach can be marred by a severe limit on accuracy unless carefully designed. We adopt machine learning techniques to search for efficient policies for the characterization of decohering quantum systems. We find, for instance, that when the system undergoes more than 2 Rabi oscillations during its relaxation time $T_1$, $O(10^3)$ measurement shots are sufficient to reduce the squared error of the Bayesian initial prior of the unknown frequency $\omega_{\rm r}$ by a factor $O(10^4)$ or larger. We also develop policies optimized for extreme initial parameter uncertainty and for the presence of imperfections in the readout.


title: Evidence for the Thermal Sunyaev-Zel'dovich Effect Associated with Quasar Feedback

abstract: Using a radio-quiet subsample of the Sloan Digital Sky Survey spectroscopic quasar catalogue, spanning redshifts 0.5-3.5, we derive the mean millimetre and far-infrared quasar spectral energy distributions (SEDs) via a stacking analysis of Atacama Cosmology Telescope and Herschel-Spectral and Photometric Imaging REceiver data. We constrain the form of the far-infrared emission and find 3$\sigma$-4$\sigma$ evidence for the thermal Sunyaev-Zel'dovich (SZ) effect, characteristic of a hot ionized gas component with thermal energy $(6.2 \pm 1.7)\times 10^{60}$ erg. This amount of thermal energy is greater than expected assuming only hot gas in virial equilibrium with the dark matter haloes of $(1-5)\times 10^{12}h^{-1}$M$_\odot$ that these systems are expected to occupy, though the highest quasar mass estimates found in the literature could explain a large fraction of this energy. Our measurements are consistent with quasars depositing up to $(14.5 \pm 3.3)~\tau_8^{-1}$ per cent of their radiative energy into their circumgalactic environment if their typical period of quasar activity is $\tau_8\times~10^8$ yr. For high quasar host masses, $\sim10^{13}h^{-1}$M$_\odot$, this percentage will be reduced. Furthermore, the uncertainty on this percentage is only statistical and additional systematic uncertainties enter at the 40 per cent level. The SEDs are dust dominated in all bands and we consider various models for dust emission. While sufficiently complex dust models can obviate the SZ effect, the SZ interpretation remains favoured at the 3$\sigma$-4$\sigma$ level for most models.


title: Non-Global Correlations in Collider Physics

abstract: Despite their importance for precision QCD calculations, correlations between in- and out-of-jet regions of phase space have never directly been observed. These so-called non-global effects are present generically whenever a collider physics measurement is not explicitly dependent on radiation throughout the entire phase space. In this paper, we introduce a novel procedure based on mutual information, which allows us to isolate these non-global correlations between measurements made in different regions of phase space. We study this procedure both analytically and in Monte Carlo simulations in the context of observables measured on hadronic final states produced in $e^+e^-$ collisions, though it is more widely applicable. The procedure exploits the sensitivity of soft radiation at large angles to non-global correlations, and we calculate these correlations through next-to-leading logarithmic accuracy. The bulk of these non-global correlations are found to be described in Monte Carlo simulation. They increase by the inclusion of non-perturbative effects, which we show can be incorporated in our calculation through the use of a model shape function. This procedure illuminates the source of non-global correlations and has connections more broadly to fundamental quantities in quantum field theory.


title: Testing the Molecular-Hydrogen Kennicutt-Schmidt Law in the Low-Density Environments of Extended Ultraviolet Disk Galaxies

abstract: Studying star formation beyond the optical radius of galaxies allows us to test empirical relations in extreme conditions with low average gas density and low molecular fraction. Previous studies discovered galaxies with extended ultraviolet (XUV) disks, which often contain star forming regions with lower Halpha-to-far-UV (FUV) flux ratios compared to inner disk star forming regions. However, most previous studies lack measurements of molecular gas, which is presumably the component of the interstellar medium out of which stars form. We analyzed published CO measurements and upper limits for fifteen star forming regions in the XUV or outer disk of three nearby spiral galaxies and a new CO upper limit from the IRAM 30 m telescope in one star forming region at r = 3.4 r_25 in the XUV disk of NGC 4625. We found that the star forming regions are in general consistent with the same molecular-hydrogen Kennicutt-Schmidt law that applies within the optical radius, independent of whether we used Halpha or FUV as the star formation rate (SFR) tracer. However, a number of the CO detections are significantly offset towards higher SFR surface density for their molecular hydrogen surface density. Deeper CO data may enable us to use the presence or absence of molecular gas as an evolutionary probe to break the degeneracy between age and stochastic sampling of the initial mass function as the explanation for the low Halpha-to-FUV flux ratios in XUV disks.


title: CLASH-VLT: Environment-driven evolution of galaxies in the z=0.209 cluster Abell 209

abstract: The analysis of galaxy properties and the relations among them and the environment, can be used to investigate the physical processes driving galaxy evolution. We study the cluster A209 by using the CLASH-VLT spectroscopic data combined with Subaru photometry, yielding to 1916 cluster members down to a stellar mass of 10^{8.6} Msun. We determine: i) the stellar mass function of star-forming and passive galaxies; ii) the intra-cluster light and its properties; iii) the orbits of low- and high-mass passive galaxies; and iv) the mass-size relation of ETGs. The stellar mass function of the star-forming galaxies does not depend on the environment, while the slope found for passive galaxies becomes flatter in the densest region. The color distribution of the intra-cluster light is consistent with the color of passive members. The analysis of the dynamical orbits shows that low-mass passive galaxies have tangential orbits, avoiding small pericenters around the BCG. The mass-size relation of low-mass passive ETGs is flatter than that of high mass galaxies, and its slope is consistent with that of field star-forming galaxies. Low-mass galaxies are also more compact within the scale radius of 0.65 Mpc. The ratio between stellar and number density profiles shows a mass segregation in the center. The comparative analysis of the stellar and total density profiles indicates that this effect is due to dynamical friction. Our results are consistent with a scenario in which the "environmental quenching" of low-mass galaxies is due to mechanisms such as harassment out to R200, starvation and ram-pressure stripping at smaller radii, as supported by the analysis of the mass function, of the dynamical orbits and of the mass-size relation of passive early-types in different regions. Our analyses support the idea that the intra-cluster light is formed through the tidal disruption of subgiant galaxies.


title: Dust Evolution Can Produce Scattered Light Gaps in Protoplanetary Disks

abstract: Recent imaging of protoplanetary disks with high resolution and contrast have revealed a striking variety of substructure. Of particular interest are cases where near-infrared scattered light images show evidence for low-intensity annular "gaps." The origins of such structures are still uncertain, but the interaction of the gas disk with planets is a common interpretation. We study the impact that the evolution of the solid material can have on the observable properties of disks in a simple scenario without any gravitational or hydrodynamical disturbances to the gas disk structure. Even with a smooth and continuous gas density profile, we find that the scattered light emission produced by small dust grains can exhibit ring-like depressions similar to those presented in recent observations. The physical mechanisms responsible for these features rely on the inefficient fragmentation of dust particles. The occurrence and position of the proposed "gap" features depend most strongly on the dust-to-gas ratio, the fragmentation threshold velocity, the strength of the turbulence, and the age of the disk, and should be generic (at some radius) for typically adopted disk parameters. The same physical processes can affect the thermal emission at optically thin wavelengths ($\sim$1 mm), although the behavior can be more complex; unlike for disk-planet interactions, a "gap" should not be present at these longer wavelengths.


title: Thermal History of the Universe After Inflation

abstract: When did the universe thermalize? In this talk I review the status of this issue and its importance in establishing the expected properties of dark matter, the growth of large-scale structure, and the viability of inflation models when confronted with CMB observations. I also present a novel approach to tackling the theoretical challenges surrounding inflationary (p)reheating, which seeks to extend past work on the Effective Field Theory of Inflation to the time of reheating.


title: Competition between intrinsic and extrinsic effects in the quenching of the superconducting state in FeSeTe thin films

abstract: We report the first experimental observation of the quenching of the superconducting state in current-voltage characteristics of an iron-based superconductor, namely, in FeSeTe thin films. Based on available theoretical models, our analysis suggests the presence of an intrinsic flux-flow electronic instability along with non-negligible extrinsic thermal effects. The coexistence and competition of these two mechanisms classify the observed instability as halfway between those of low-temperature and of high-temperature superconductors, where thermal effects are respectively largely negligible or predominant.


title: $S$-duality of $u(1)$ gauge theory with $\theta =\pi$ on non-orientable manifolds: Applications to topological insulators and superconductors

abstract: Electric-magnetic duality ($S$-duality) is a well-known property of pure $u(1)$ gauge theory in 3+1 dimensions. In this paper, we investigate the compatibility of this duality with time-reversal symmetry. We consider two theories obtained by coupling a Dirac fermion with an "inverted" sign of the mass $m$ to a $u(1)$ gauge field. Time-reversal in the two theories is implemented respectively via the $T$ and $CT$ symmetries of the Dirac fermion. It was recently conjectured (C. Wang and T. Senthil (arXiv:1505.03520), and M. Metlitski and A.Vishwanath (arXiv:1505.05142)) that in the $|m| \to \infty$ limit these two theories are $S$-dual to each other. We provide support for this conjecture by studying partition functions of the two theories on non-orientable manifolds as a way to probe the realization of time-reversal. Upon integrating out the Dirac fermion, topological terms in the actions of the two theories are generated. While on an orientable manifold topological terms in both theories reduce to a $\theta$-term with $\theta = \pi$, on a non-orientable manifold they are distinct. We explicitly compute partition functions of the two theories on the manifold $\mathbb{RP}^4$ and show that they are equal; this result combined with certain physical arguments is sufficient to establish the duality. The two theories can be viewed as a gauged topological insulator in class AII and a gauged topological superconductor in class AIII, and the bulk duality allows us to derive previously conjectured non-trivial symmetric gapped surface states of these phases.


title: The Impact of Environmental Fluctuations on Evolutionary Fitness Functions

abstract: The concept of fitness as a measure for a species's success in natural selection is central to the theory of evolution. We here investigate how reproduction rates which are not constant but vary in response to environmental fluctuations, influence a species' prosperity and thereby its fitness. Interestingly, we find that not only larger growth rates but also reduced sensitivities to environmental changes substantially increase the fitness. Thereby, depending on the noise level of the environment, it might be an evolutionary successful strategy to minimize this sensitivity rather than to optimize the reproduction speed. Also for neutral evolution, where species with exactly the same properties compete, variability in the growth rates plays a crucial role. The time for one species to fixate is strongly reduced in the presence of environmental noise. Hence, environmental fluctuations constitute a possible explanation for effective population sizes inferred from genetic data that often are much smaller than the census population size.


title: SPT 2-Channel Kondo Model in the Structure of Normal Metal/Quantum Dot/DIII-class Topological Superconductor

abstract: We investigate the Kondo effect in a structure which is constructed by embedding one quantum dot between a normal metal and a \emph{DIII}-class topological superconductor supporting Majorana doublets at its ends. It is observed that Kondo correlation occurs between the localized state in the dot and two continuum states simultaneously, i.e., the continuum state in the metal and the continuum Andreev reflection state between the metal and topological superconductor. As a result, the Kondo model Hamiltonian is topologically protected by the $SU(2)_s\rtimes Z_2^T$ symmetry. Besides, the Kondo temperature has an opportunity to present three forms in turn, following the enhancement of the coupling between the dot and Majorana doublet. This phenomenon exactly reflects the special role of Majorana doublet in tuning the Kondo effect.


title: The Catalog of Earth-Like Exoplanet Survey TArgets (CELESTA): A Database of Habitable Zones around Nearby Stars

abstract: Locating planets in circumstellar Habitable Zones is a priority for many exoplanet surveys. Space-based and ground-based surveys alike require robust toolsets to aid in target selection and mission planning. We present the Catalog of Earth-Like Exoplanet Survey Targets (CELESTA), a database of Habitable Zones around 36,000 nearby stars. We calculated stellar parameters, including effective temperatures, masses, and radii, and we quantified the orbital distances and periods corresponding to the circumstellar Habitable Zones. We gauged the accuracy of our predictions by contrasting CELESTA's computed parameters to observational data. We ascertain a potential return on investment by computing the number of Habitable Zones probed for a given survey duration. A versatile framework for extending the functionality of CELESTA into the future enables ongoing comparisons to new observations, and recalculations when updates to Habitable Zone models, stellar temperatures, or parallax data become available. We expect to upgrade and expand CELESTA using data from the Gaia mission as the data becomes available.


title: The dust content of damped Lyman-alpha systems in the Sloan Digital Sky Survey

abstract: The dust-content of damped Lyman-alpha systems (DLAs) is an important observable for understanding their origin and the neutral gas reservoirs of galaxies. While the average colour-excess of DLAs, E(B-V), is known to be <15 milli-magnitudes (mmag), both detections and non-detections with ~2 mmag precision have been reported. Here we find 3.2-sigma statistical evidence for DLA dust-reddening of 774 Sloan Digital Sky Survey (SDSS) quasars by comparing their fitted spectral slopes to those of ~7000 control quasars. The corresponding E(B-V) is 3.0 +/- 1.0 mmag, assuming a Small Magellanic Cloud (SMC) dust extinction law, and it correlates strongly (3.5-sigma) with the metal content, characterised by the SiII1526 absorption-line equivalent width, providing additional confidence that the detection is due to dust in the DLAs. Evolution of E(B-V) over the redshift range 2.1 < z < 4.0 is limited to <2.5 mmag per unit redshift (1-sigma), consistent with the known, mild DLA metallicity evolution. There is also no apparent relationship with neutral hydrogen column density, N(HI), though the data are consistent with a mean E(B-V)/N(HI) = (3.5 +/- 1.0) x 10^{-24} mag cm^2, approximately the ratio expected from the SMC scaled to the lower metallicities typical of DLAs. We implement the SDSS selection algorithm in a portable code to assess the potential for systematic, redshift-dependent biases stemming from its magnitude and colour-selection criteria. The effect on the mean E(B-V) is negligible (<5 per cent) over the entire redshift range of interest. Given the broad potential usefulness of this implementation, we make it publicly available.


title: On the Use of the Number Count of Blue Horizontal--Branch Stars to Infer the Dominant Building Blocks of the Milky Way Halo

abstract: The formation of the Milky Way stellar halo is thought to be the result of merging and accretion of building blocks such as dwarf galaxies and massive globular clusters. Recently, Deason et al. (2015) suggested that the Milky Way outer halo formed mostly from big building blocks, such as dwarf spheroidal galaxies, based on the similar number ratio of blue straggler (BS) stars to blue horizontal-branch (BHB) stars. Here we demonstrate, however, that this result is seriously biased by not taking into detailed consideration on the formation mechanism of BHB stars from helium enhanced second-generation population. In particular, the high BS-to-BHB ratio observed in the outer halo fields is most likely due to a small number of BHB stars provided by GCs rather than to a large number of BS stars. This is supported by our dynamical evolution model of GCs which shows preferential removal of first generation stars in GCs. Moreover, there are a sufficient number of outer halo GCs which show very high BS-to-BHB ratio. Therefore, the BS-to-BHB number ratio is not a good indicator to use in arguing that more massive dwarf galaxies are the main building blocks of the Milky Way outer halo. Several lines of evidence still suggest that GCs can contribute a significant fraction of the outer halo stars.


title: Radiative plateau inflation

abstract: We describe how monomial chaotic inflation becomes compatible with the latest CMB data thanks to radiative corrections producing a plateau. The interactions of the inflation with other fields, required for reheating, can flatten the potential and moderate the production of primordial gravitational waves, keeping these below the current upper bound. We show that the appearance of a plateau requires that the inflaton couples to fermions and to another scalar or a gauge group. We give concrete examples of minimal particle physics models leading to plateaus for quadratic and quartic chaotic inflation. We also provide a three-parameter model-independent description of radiatively corrected inflation that is amenable to CMB analyses.


title: Noiseless manipulation of helical edge state transport by a quantum magnet

abstract: The current through a helical edge state of a quantum-spin-Hall insulator may be fully transmitted through a magnetically gapped region due to a combination of spin-transfer torque and spin pumping [Meng {\em et al.}, Phys. Rev. B {\bf 90}, 205403 (2014)]. Using a scattering approach, we here argue that in such a system the current is effectively carried by electrons with energies below the magnet-induced gap and well below the Fermi energy. This has striking consequences, such as the absence of shot noise, an exponential suppression of thermal noise, and an obstruction of thermal transport. For two helical edges covered by the same quantum magnet, the device can act as robust noiseless current splitter.


title: Formation of globular clusters in atomic-cooling halos via rapid gas condensation and fragmentation during the epoch of reionization

abstract: We investigate the formation of metal-poor globular clusters (GCs) at the center of two dark matter halos with $M_{\rm halo} \sim4\times10^7 M_\odot$ at $z>10$ using cosmological radiation-hydrodynamics simulations. We find that very compact ($\lesssim$ 1 pc) and massive ($\sim6\times10^5 M_\odot$) clusters form rapidly when pristine gas collapses isothermally with the aid of efficient Ly$\alpha$ emission during the transition from molecular-cooling halos to atomic-cooling halos. Because the local free-fall time of dense star-forming gas is very short ($\ll 1\,{\rm Myr}$), a large fraction of the collapsed gas is turned into stars before stellar feedback processes blow out the gas and shut down star formation. Although the early stage of star formation is limited to a small region of the central star-forming disk, we find that the disk quickly fragments due to metal enrichment from supernovae. Sub-clusters formed in the fragmented clouds eventually merge with the main cluster at the center. The simulated clusters closely resemble the local GCs in mass and size but show a metallicity spread that is much wider than found in the local GCs. We discuss a role of pre-enrichment by Pop III and II stars as a potential solution to the latter issue. Although not without shortcomings, it is encouraging that a naive blind (not tuned) cosmological simulation presents a possible channel for the formation of at least some massive GCs.


title: Nonsingular transformations and dimension spaces

abstract: For any adic transformation $T$ defined on the path space $X$ of an ordered Bratteli diagram, endowed with a Markov measure $\mu$, we construct an explicit dimension space (which corresponds to a matrix values random walk on $\mathbb{Z}$) whose Poisson boundary can be identified as a $\mathbb{Z}$-space with the dynamical system $(X,\mu,T)$. We give a couple of examples to show how dimension spaces can be used in the study of nonsingular transformations.


title: Update on the GRB universal scaling E$_{\rm{X,iso}}$-E$_{\rm{\gamma,iso}}$-E$_{\rm{pk}}$ with ten years of $Swift$ data

abstract: From a comprehensive statistical analysis of $Swift$ X-ray light-curves of gamma-ray bursts (GRBs) collected from December 2004 to the end of 2010, we found a three-parameter correlation between the isotropic energy emitted in the rest frame 1-10$^4$ keV energy band during the prompt emission (E$_{\rm{\gamma,iso}}$), the rest frame peak of the prompt emission energy spectrum (E$_{\rm{pk}}$), and the X-ray energy emitted in the rest frame 0.3-30 keV observed energy band (E$_{\rm{X,iso}}$), computed excluding the contribution of the flares. In this paper, we update this correlation with the data collected until June 2014, expanding the sample size with $\sim$35% more objects, where the number of short GRBs doubled. With this larger sample we confirm the existence of a universal correlation that connects the prompt and afterglow properties of long and short GRBs. We show that this correlation does not depend on the X-ray light-curve morphology and that further analysis is necessary to firmly exclude possible biases derived by redshift measurements. In addition we discuss about the behavior of the peculiar objects as ultra-long GRBs and we propose the existence of an intermediate group between long and short GRBs. Interestingly, two GRBs with uncertain classification fall into this category. Finally, we discuss the physics underlying this correlation, in the contest of the efficiency of conversion of the prompt $\gamma$-ray emission energy into the kinetic energy of the afterglow, the photosferic model, and the cannonball model.


title: Explicit formulas for infinitely many Shimura curves in genus 4

abstract: In this paper we construct infinitely many Shimura curves contained in the locus of Jacobians of genus four curves. All Jacobians in these families are ${\mathbb Z}/3$ covers of varying elliptic curves that appear in a geometric construction of Pirola, and include an example of a Shimura-Teichmuller curve that parameterizes Jacobians that are suitable ${\mathbb Z}/6$ covers of ${\mathbb P}^1$. We compute explicitly the period matrices of the Shimura curves we construct using the original construction of Shimura for moduli spaces of abelian varieties with automorphisms.


title: Atmospheric Cherenkov Gamma-ray Telescopes

abstract: The stereoscopic imaging atmospheric Cherenkov technique, developed in the 1980s and 1990s, is now used by a number of existing and planned gamma-ray observatories around the world. It provides the most sensitive view of the very high energy gamma-ray sky (above 30 GeV), coupled with relatively good angular and spectral resolution over a wide field-of-view. This Chapter summarizes the details of the technique, including descriptions of the telescope optical systems and cameras, as well as the most common approaches to data analysis and gamma-ray reconstruction.


title: Uniform semiclassical approximations for one-dimensional fermionic systems

abstract: A thorough account is given of the derivation of uniform semiclassical approximations to the particle and kinetic energy densities of N noninteracting bounded fermions in one dimension. The employed methodology allows the inclusion of non-perturbative effects via an infinite resummation of the Poisson summation formula.


title: Structure estimation for mixed graphical models in high-dimensional data

abstract: Undirected graphical models are a key component in the analysis of complex observational data in a large variety of disciplines. In many of these applications one is interested in estimating the undirected graphical model underlying a distribution over variables with different domains. Despite the pervasive need for such an estimation method, to date there is no such method that models all variables on their proper domain. We close this methodological gap by combining a new class of mixed graphical models with a structure estimation approach based on generalized covariance matrices. We report the performance of our methods using simulations, illustrate the method with a dataset on Autism Spectrum Disorder (ASD) and provide an implementation as an R-package.


title: Does Unitary Measurement Theory Lead to an Everettian Interpretation?

abstract: Quantum-mechanical interpretation-related implications of the theory of unitary premeasurement [1] on complete measurement (objectification or collapse included) are investigated in the present article with a view to give an affirmative answer to the question in the title. It is argued that both Bohr's and von Neumann's ideas lead to those of Everett. Hence, the latter can be, in some sense, considered to be a continuation and elaboration of both former approaches. The importance of the idea of relativeness in Everett's theory is emphasized. To free the relative-state theory from its roots both of classicalness and of subjective observation in the argument of this study, the general or unfolded version of Everett's theory is sketched.


title: Borel Complexity and Potential Canonical Scott Sentences

abstract: We define and investigate HC-forcing invariant formulas of set theory, whose interpretations in the hereditarily countable sets are well behaved under forcing extensions. This leads naturally to a notion of cardinality ||Phi|| for sentences Phi of $L_{\omega_1,\omega}$, which counts the number of sentences of $L_{\infty,\omega}$ that, in some forcing extension, become a canonical Scott sentence of a model of Phi. We show this cardinal bounds the complexity of (Mod(Phi), iso), the class of models of Phi with universe omega, by proving that (Mod(Phi),iso) is not Borel reducible to (Mod(Psi),iso) whenever ||Psi|| < ||Phi||. Using these tools, we analyze the complexity of the class of countable models of four complete, first-order theories T for which (Mod(T),iso) is properly analytic, yet admit very different behavior. We prove that both `Binary splitting, refining equivalence relations' and Koerwien's example of an eni-depth 2, omega-stable theory have (Mod(T),iso) non-Borel, yet neither is Borel complete. We give a slight modification of Koerwien's example that also is omega-stable, eni-depth 2, but is Borel complete. Additionally, we prove that I_{\infty,\omega}(Phi)<\beth_{\omega_1} whenever (Mod(Phi),iso) is Borel.


title: Dense Gas Fraction and Star Formation Efficiency Variations in the Antennae Galaxies

abstract: We use the CARMA millimeter interferometer to map the Antennae Galaxies (NGC4038/39), tracing the bulk of the molecular gas via the 12CO(1-0) line and denser molecular gas via the high density transitions HCN(1-0), HCO+(1-0), CS(2-1), and HNC(1-0). We detect bright emission from all tracers in both the two nuclei and three locales in the overlap region between the two nuclei. These three overlap region peaks correspond to previously identified "supergiant molecular clouds". We combine the CARMA data with Herschel infrared (IR) data to compare observational indicators of the star formation efficiency (SFR/H2~IR/CO), dense gas fraction (HCN/CO), and dense gas star formation efficiency (IR/HCN). Regions within the Antennae show ratios consistent with those seen for entire galaxies, but these ratios vary by up to a factor of 6 within the galaxy. The five detected regions vary strongly in both their integrated intensities and these ratios. The northern nucleus is the brightest region in mm-wave line emission, while the overlap region is the brightest part of the system in the IR. We combine the CARMA and Herschel data with ALMA CO data to report line ratio patterns for each bright point. CO shows a declining spectral line energy distribution, consistent with previous studies. HCO+(1-0) emission is stronger than HCN(1-0) emission, perhaps indicating either more gas at moderate densities or higher optical depth than is commonly seen in more advanced mergers.


title: Server Placement with Shared Backups for Disaster-Resilient Clouds

abstract: A key strategy to build disaster-resilient clouds is to employ backups of virtual machines in a geo-distributed infrastructure. Today, the continuous and acknowledged replication of virtual machines in different servers is a service provided by different hypervisors. This strategy guarantees that the virtual machines will have no loss of disk and memory content if a disaster occurs, at a cost of strict bandwidth and latency requirements. Considering this kind of service, in this work, we propose an optimization problem to place servers in a wide area network. The goal is to guarantee that backup machines do not fail at the same time as their primary counterparts. In addition, by using virtualization, we also aim to reduce the amount of backup servers required. The optimal results, achieved in real topologies, reduce the number of backup servers by at least 40%. Moreover, this work highlights several characteristics of the backup service according to the employed network, such as the fulfillment of latency requirements.


title: Protein Structure Prediction by Protein Alignments

abstract: Proteins are the basic building blocks of life. They usually perform functions by folding to a particular structure. Understanding the folding process could help the researchers to understand the functions of proteins and could also help to develop supplemental proteins for people with deficiencies and gain more insight into diseases associated with troublesome folding proteins. Experimental methods are both expensive and time consuming. In this thesis I introduce a new machine learning based method to predict the protein structure. The new method improves the performance from two directions: creating accurate protein alignments and predicting accurate protein contacts. First, I present an alignment framework MRFalign which goes beyond state-of-the-art methods and uses Markov Random Fields to model a protein family and align two proteins by aligning two MRFs together. Compared to other methods, that can only model local-range residue correlation, MRFs can model long-range residue interactions and thus, encodes global information in a protein. Secondly, I present a Group Graphical Lasso method for contact prediction that integrates joint multi-family Evolutionary Coupling analysis and supervised learning to improve accuracy on proteins without many sequence homologs. Different from single-family EC analysis that uses residue co-evolution information in only the target protein family, our joint EC analysis uses residue co-evolution in both the target family and its related families, which may have divergent sequences but similar folds. Our method can also integrate supervised learning methods to further improve accuracy. We evaluate the performance of both methods including each of its components on large public benchmarks. Experiments show that our methods can achieve better accuracy than existing state-of-the-art methods under all the measurements on most of the protein classes.


title: A characterization of modified mock theta functions

abstract: We give a characterization of modified (in the sense of Zwegers) mock theta functions, parallel to that of ordinary theta functions. Namely, modified mock theta functions are characterized by their analyticity properties, elliptic transformation properties, and by being annihilated by certain second order differential operators.


title: NYTRO: When Subsampling Meets Early Stopping

abstract: Early stopping is a well known approach to reduce the time complexity for performing training and model selection of large scale learning machines. On the other hand, memory/space (rather than time) complexity is the main constraint in many applications, and randomized subsampling techniques have been proposed to tackle this issue. In this paper we ask whether early stopping and subsampling ideas can be combined in a fruitful way. We consider the question in a least squares regression setting and propose a form of randomized iterative regularization based on early stopping and subsampling. In this context, we analyze the statistical and computational properties of the proposed method. Theoretical results are complemented and validated by a thorough experimental analysis.


title: Surface waves in granular phononic crystals

abstract: The existence of surface elastic waves at a mechanically free surface of granular phononic crystals is studied. The granular phononic crystals are made of spherical particles distributed periodically on a simple cubic lattice. It is assumed that the particles are interacting by means of normal, shear and bending contact rigidities. First, Rayleigh-type surface acoustic waves, where the displacement of the particles takes place in the sagittal plane while the particles possess one rotational and two translational degrees of freedom, are analyzed. Second, shear-horizontal-type waves, where the displacement of the particles is normal to the sagittal plane while the particles possess one translational and two rotational degrees of freedom are studied. The existence of zero-group velocity surface acoustic waves of Rayleigh-type is theoretically predicted and interpreted. A comparison with surface waves predicted by the Cosserat theory is performed, and its limitations are established.


title: New Phases and Dissociation-Recombination of Hydrogen Deuteride to 3.4 Megabar

abstract: We present infrared absorption studies of solid hydrogen deuteride to pressures as high as 3.4 megabar in a diamond anvil cell and temperatures in the range 5 to 295 K. Above 198 GPa the sample transforms to a mixture of HD ,H2 and D2, interpreted as a process of dissociation and recombination. Three new phases-lines are observed, two of which differ remarkably from those of the high-pressure homonuclear species, but none are metallic. The time-dependent spectral changes are analyzed to determine the molecular concentrations as a function of time; the nucleon exchange achieves steady state concentrations in ~20 hours at ~200 GPa.


title: Moduli Interpretations for Noncongruence Modular Curves

abstract: We define the notion of a $G$-structure for elliptic curves, where $G$ is a finite 2-generated group. When $G$ is abelian, a $G$-structure is the same as a classical congruence level structure. There is a natural action of $\text{SL}_2(\mathbb{Z})$ on these level structures. If $\Gamma$ is a stabilizer of this action, then the quotient of the upper half plane by $\Gamma$ parametrizes isomorphism classes of elliptic curves equipped with $G$-structures. When $G$ is "sufficiently" nonabelian, the stabilizers $\Gamma$ are noncongruence. As a result we realize noncongruence modular curves as moduli spaces of elliptic curves equipped with nonabelian $G$-structures. As applications we describe links to the Inverse Galois Problem, and show how our moduli interpretations explains the bad primes for the Unbounded Denominators Conjecture, and allows us to translate the conjecture into the language of geometry and Galois theory.