title: FIESTA 4: optimized Feynman integral calculations with GPU support

abstract: This paper presents a new major release of the program FIESTA (Feynman Integral Evaluation by a Sector decomposiTion Approach). The new release is mainly aimed at optimal performance at large scales when one is increasing the number of sampling points in order to reduce the uncertainty estimates. The release now supports graphical processor units (GPU) for the numerical integration, methods to optimize cluster-usage, as well as other speed, memory, and stability improvements.


title: Inferring the Distances of Fast Radio Bursts Through Associated 21-cm Absorption

abstract: The distances of Fast Radio Burst (FRB) sources are currently unknown. We show that the 21-cm absorption line of hydrogen can be used to infer the redshifts of FRB sources, and determine whether they are Galactic or extragalactic. We calculate a probability of $\sim 10\%$ for the host galaxy of an FRB to exhibit a 21-cm absorption feature of equivalent width $\gtrsim 10 ~\mathrm{km}~\mathrm{s}^{-1}$. Arecibo, along with several future radio observatories, should be capable of detecting such associated 21-cm absorption signals for strong bursts of $\gtrsim {\rm several}~\mathrm{Jy}$ peak flux densities.


title: Moral hazard under ambiguity

abstract: In this paper, we extend the Holmstro\"om and Milgrom problem [47] by adding uncertainty about the volatility of the output for both the Agent and the Principal. We study more precisely the impact of the "Nature" playing against the Agent and the Principal by choosing the worst possible volatility of the output. We solve the first--best and the second--best problems associated with this framework and we show that optimal contracts are in a class of contracts similar to [14, 15], linear with respect to the output and its quadratic variation. We compare our results with the classical problem in [47].


title: The Gini Coefficient as a Morphological Measurement of Strongly Lensed Galaxies in the Image Plane

abstract: Characterization of the morphology of strongly lensed galaxies is challenging because images of such galaxies are typically highly distorted. Lens modeling and source plane reconstruction is one approach that can provide reasonably undistorted images from which morphological measurements can be made, although at the expense of a highly spatially variable telescope PSF when mapped back to the source plane. Unfortunately, modeling the lensing mass is a time and resource intensive process, and in many cases there are too few constraints to precisely model the lensing mass. If, however, useful morphological measurements could be made in the image plane rather than the source plane, it would bypass this issue and obviate the need for a source reconstruction process for some applications. We examine the use of the Gini coefficient as one such measurement. Because it depends on the cumulative distribution of the light of a galaxy, but not the relative spatial positions, the fact that surface brightness is conserved by lensing means that the Gini coefficient may be well-preserved by strong gravitational lensing. Through simulations, we test the extent to which the Gini coefficient is conserved, including by effects due to PSF convolution and pixelization, to determine whether it is invariant enough under lensing to be used as a measurement of galaxy morphology that can be made in the image plane.


title: Bayesian parameter estimation for effective field theories

abstract: We present procedures based on Bayesian statistics for estimating, from data, the parameters of effective field theories (EFTs). The extraction of low-energy constants (LECs) is guided by theoretical expectations in a quantifiable way through the specification of Bayesian priors. A prior for natural-sized LECs reduces the possibility of overfitting, and leads to a consistent accounting of different sources of uncertainty. A set of diagnostic tools are developed that analyze the fit and ensure that the priors do not bias the EFT parameter estimation. The procedures are illustrated using representative model problems, including the extraction of LECs for the nucleon mass expansion in SU(2) chiral perturbation theory from synthetic lattice data.


title: Modular Invariants of a Vector and a Covector: a proof of a conjecture of Bonnaf\'e and Kemper

abstract: Consider a finite dimensional vector space $V$ over a finite field $\mathbb{F}_q$. We give a minimal generating set for the ring of invariants $\mathbb{F}_q[V \oplus V^*]^{\text{GL}(V)}$, and show that this ring is a Gorenstein ring but is not a complete intersection. These results confirm a conjecture of Bonnaf\'e and Kemper.


title: Probing proton fluctuations with asymmetric rapidity correlations

abstract: Intrinsic fluctuations of the proton saturation momentum generate asymmetric rapidity distributions on an event-by-event basis. We argue that the asymmetric component, $\left<a_1^2\right>$, of the orthogonal polynomial decomposition of the two-particle rapidity correlation function is a sensitive probe to this distribution of fluctuations. We present a simple model connecting the experimentally measured $\left<a_1^2\right>$ to the variance, $\sigma$, of the distribution of the logarithm of the proton saturation scale. We find that $\sigma\approx 0. 5-1$ describes the asymmetric component of the rapidity correlations recently measured by the ATLAS collaboration.


title: Chiral anomaly without relativity

abstract: Perspectives article on J. Xiong et al., Science 350, 413 (2015).


title: Weak index pairs and the Conley index for discrete multivalued dynamical systems

abstract: Motivated by the problem of reconstructing dynamics from samples we revisit the Conley index theory for discrete multivalued dynamical systems. We introduce a new, less restrictive definition of the isolating neighbourhood. It turns out that then the main tool for the construction of the index, i.e. the index pair, is no longer useful. In order to overcome this obstacle we use the concept of weak index pairs.


title: Inclusion Matrices and the MDS Conjecture

abstract: Let F_q be a finite field of order q with characteristic p. An arc is an ordered family of at least k vectors in (F_q)^k in which every subfamily of size k is a basis of (F_q)^k. The MDS conjecture, which was posed by Segre in 1955, states that if k <= q, then an arc in (F_q)^k has size at most q+1, unless q is even and k=3 or k=q-1, in which case it has size at most q+2. We propose a conjecture which would imply that the MDS conjecture is true for almost all values of k when q is odd. We prove our conjecture in two cases and thus give simpler proofs of the MDS conjecture when k <= p, and if q is not prime, for k <= 2p-2. To accomplish this, given an arc G of (F_q)^k and a nonnegative integer n, we construct a matrix M_G^{\uparrow n}, which is related to an inclusion matrix, a well-studied object in combinatorics. Our main results relate algebraic properties of the matrix M_G^{\uparrow n} to properties of the arc G and may provide new tools in the computational classification of large arcs.


title: $B$-Rigidity of flag $2$-spheres without $4$-belt

abstract: Associated to every finite simplicial complex $K$, there is a moment-angle complex $\mathcal {Z}_{K}$. In this paper, we use some algebraic invariants to solve the $B$-rigidity problem for some special simplicial compelexes.


title: Three-point correlator of twist-2 light-ray operators in N=4 SYM in BFKL approximation

abstract: We present calculation of the correlation function of three twist-2 operators in the BFKL limit. The calculation is performed in N = 4 SYM but the result is valid in other gauge theories such as QCD. The obtained leading order structure constant is exact for any number of colors.


title: Photoinduced suppression of the ferroelectric instability in PbTe

abstract: The interactions between electrons and phonons drive a large array of technologically relevant material properties including ferroelectricity, thermoelectricity, and phase-change behaviour. In the case of many group IV-VI, V, and related materials, these interactions are strong and the materials exist near electronic and structural phase transitions. Their close proximity to phase instability produces a fragile balance among the various properties. The prototypical example is PbTe whose incipient ferroelectric behaviour has been associated with large phonon anharmonicity and thermoelectricity. Experimental measurements on PbTe reveal anomalous lattice dynamics, especially in the soft transverse optical phonon branch. This has been interpreted in terms of both giant anharmonicity and local symmetry breaking due to off-centering of the Pb ions. The observed anomalies have prompted renewed theoretical and computational interest, which has in turn revived focus on the extent that electron-phonon interactions drive lattice instabilities in PbTe and related materials. Here, we use Fourier-transform inelastic x-ray scattering (FT-IXS) to show that photo-injection of free carriers stabilizes the paraelectric state. With support from constrained density functional theory (CDFT) calculations, we find that photoexcitation weakens the long-range forces along the cubic direction tied to resonant bonding and incipient ferroelectricity. This demonstrates the importance of electronic states near the band edges in determining the equilibrium structure.


title: Biexciton-mediated superradiant photon blockade

abstract: The photon blockade is a hallmark of quantum light transport through a single two-level system that can accomodate only one photon. Here, we theoretically show that two-photon transmission can be suppressed even for a seemingly classical system with large number of quantum dots in a cavity when the biexciton nonlinearity is taken into account. We reveal the nonmonotonous dependence of the second-order correlation function of the transmitted photons on the biexciton binding energy. The blockade is realized by proper tuning the biexciton resonance that controls the collective superradiant modes.


title: Computing optimal interfacial structure of ordered phases

abstract: We propose a general framework of computing interfacial structure. If an ordered phase is involved, the interfacial structure can be obtained by simply minimizing the free energy with compatible boundary conditions. The framework is applied to Landau- Brazovskii model and works efficiently.


title: A Continuous Max-Flow Approach to Cyclic Field Reconstruction

abstract: Reconstruction of an image from noisy data using Markov Random Field theory has been explored by both the graph-cuts and continuous max-flow community in the form of the Potts and Ishikawa models. However, neither model takes into account the particular cyclic topology of specific intensity types such as the hue in natural colour images, or the phase in complex valued MRI. This paper presents \textit{cyclic continuous max-flow} image reconstruction which models the intensity being reconstructed as having a fundamentally cyclic topology. This model complements the Ishikawa model in that it is designed with image reconstruction in mind, having the topology of the intensity space inherent in the model while being readily extendable to an arbitrary intensity resolution.


title: A Feynman Integral and its Recurrences and Associators

abstract: We determine closed and compact expressions for the epsilon-expansion of certain Gaussian hypergeometric functions expanded around half-integer values by explicitly solving for their recurrence relations. This epsilon-expansion is identified with the normalized solution of the underlying Fuchs system of four regular singular points. We compute its regularized zeta series (giving rise to two independent associators) whose ratio gives the epsilon-expansion at a specific value. Furthermore, we use the well known one-loop massive bubble integral as an example to demonstrate how to obtain all-order epsilon-expansions for Feynman integrals and how to construct representations for Feynman integrals in terms of generalized hypergeometric functions. We use the method of differential equations in combination with the recently established general solution for recurrence relations with non-commutative coefficients.


title: General normal forms for any additive logic

abstract: In this article, we define general normal forms for any logic that has propositional part and whose non-propositional connectives distribute over the finite disjunctions. We do not require the non-propositional connectives to be closed on the set of formulas, so our normal forms cover logics with partial connectives too. We also show that most of the known normal forms in the literature are in fact particular cases of our general forms. These general normal forms are natural improvement of the distributive normal forms of J. Hintikka and their modal analogues, e.g. [Anderson] and [Fine].


title: Higher Hybrid Bottomonia in an Extended Potential Model

abstract: Using our extension of the quark potential model to hybrid mesons that fits well to the available lattice results, we now calculate the masses, radii, wave functions at origin, leptonic and two photon decay widths, E1 and M1 radiative transitions for a significant number of bottomonium mesons. These mesons include both conventional and hybrid ones with radial and angular excitations. Our numerical solutions of the Schrodinger equation are related to QCD through the Born-Oppenheimer approach. Relativistic corrections in masses and decay widths are also calculated by applying the leading order perturbation theory. The calculated results are compared with available experimental data and the theoretical results by other groups. We also identify the states of $\Upsilon(10860)$, $\Upsilon(11020)$, and $Y_b(10890)$ mesons by comparing their experimental masses and decay widths with our results.


title: Expected Regularized Total Variation of Brownian Motion

abstract: We introduce a notion of regularized total variation on an interval for continuous functions with unbounded variation. The definition of regularized total variation is obtained from that of total variation by subtracting a penalty for the size of the partition used to estimate the variation. We present an explicit construction of a partition achieving the regularized total variation, and use this construction to estimate the expected regularized total variation of Brownian motion on an interval.


title: Delayed-feedback chimera states: Forced multiclusters and stochastic resonance

abstract: A nonlinear oscillator model with negative time-delayed feedback is studied numeri- cally under external deterministic and stochastic forcing. It is found that in the unforced system complex partial synchronization patterns like chimera states as well as salt-and-pepper like soli- tary states arise on the route from regular dynamics to spatio-temporal chaos. The control of the dynamics by external periodic forcing is demonstrated by numerical simulations. It is shown that one-cluster and multi-cluster chimeras can be achieved by adjusting the external forcing frequency to appropriate resonance conditions. If a stochastic component is superimposed to the determin- istic external forcing, chimera states can be induced in a way similar to stochastic resonance, they appear, therefore, in regimes where they do not exist without noise.


title: Fisher Matrix Optimization of Cosmic Microwave Background Interferometers

abstract: We describe a method for forecasting errors in interferometric measurements of polarization of the cosmic microwave background (CMB) radiation, based on the use of the Fisher matrix calculated from the visibility covariance and relation matrices. In addition to noise and sample variance, the method can account for many kinds of systematic error by calculating an augmented Fisher matrix, including parameters that characterize the instrument along with the cosmological parameters to be estimated. The method is illustrated with examples of gain errors and errors in polarizer orientation. The augmented Fisher matrix approach is applicable to a much wider range of problems beyond CMB interferometry.


title: The quantum (non-Abelian) Potts model and its exact solution

abstract: We generalize the classical one dimensional Potts model to the case where the symmetry group is a non-Abelian finite group. It turns out that this new model has a quantum nature in that its spectrum of energy eigenstates consists of entangled states. We determine the complete energy spectrum, i.e. the ground states and all the excited states with their degeneracy structure. We calculate the partition function by two different algebraic and combinatorial methods. We also determine the entanglement properties of its ground states.


title: Homogeneous M2 duals

abstract: Motivated by the search for new gravity duals to M2 branes with $N>4$ supersymmetry --- equivalently, M-theory backgrounds with Killing superalgebra $\mathfrak{osp}(N|4)$ for $N>4$ --- we classify homogeneous M-theory backgrounds with symmetry Lie algebra $\mathfrak{so}(n) \oplus \mathfrak{so}(3,2)$ for $n=5,6,7$. We find that there are no new backgrounds with $n=6,7$ but we do find a number of new (to us) backgrounds with $n=5$. All backgrounds are metrically products of the form $\operatorname{AdS}_4 \times P^7$, with $P$ riemannian and homogeneous under the action of $\operatorname{SO}(5)$, or $S^4 \times Q^7$ with $Q$ lorentzian and homogeneous under the action of $\operatorname{SO}(3,2)$. At least one of the new backgrounds is supersymmetric (albeit with only $N=2$) and we show that it can be constructed from a supersymmetric Freund--Rubin background via a Wick rotation. Two of the new backgrounds have only been approximated numerically. (The second version of this paper includes an appendix by Alexander~S.~Haupt, closing a gap in our original analysis.)


title: Are Game Platforms suitable for Parkinson Disease patients?

abstract: Parkinson's Disease (PD) is a progressive neurodegenerative movement disorder that affects more that 6 million people worldwide. Motor dysfunction gradually increases as the disease progress. It is usually mild in the early stages of the disease but it relentlessly progresses to a severe or very severe disability that is characterized by increasing degrees of bradykinesia, hypokinesia, muscle rigidity, loss of postural reflexes and balance control as well as freezing of gait. In addition to a line of treatment based on dopaminergic PD-specific drugs, attending neurologists strongly recommend regular exercise combined with physiotherapy. However, the routine of traditional rehabilitation often create boredom and loss of interest. Opportunities to liven up a daily exercise schedule may well take the form of character-based virtual reality games which engage the player to physically train in a non-linear and looser fashion, providing an experience that varies from one game loop the next. Such "exergames", a word that results from the amalgamation of the words "exercise" and "game" challenge patients into performing movements of varying complexity in a playful and immersive virtual environment. In fact, today's game consoles using controllers like Nintendo's Wii, Sony PlayStation Eye and the Microsoft Kinect sensor present new opportunities to infuse motivation and variety to an otherwise mundane physiotherapy routine. But are these controllers and the games built for them appropriate for PD patients? In this paper we present some of these approaches and discuss their suitability for these patients mainly on the basis of demands made on balance, agility and gesture precision.


title: Analysis of Intel's Haswell Microarchitecture Using The ECM Model and Microbenchmarks

abstract: This paper presents an in-depth analysis of Intel's Haswell microarchitecture for streaming loop kernels. Among the new features examined is the dual-ring Uncore design, Cluster-on-Die mode, Uncore Frequency Scaling, core improvements as new and improved execution units, as well as improvements throughout the memory hierarchy. The Execution-Cache-Memory diagnostic performance model is used together with a generic set of microbenchmarks to quantify the efficiency of the microarchitecture. The set of microbenchmarks is chosen such that it can serve as a blueprint for other streaming loop kernels.


title: Combining features of the Unreal and Unity Game Engines to hone development skills

abstract: Two of the most popular game engines today, Unreal Engine v4.x and Unity Game Engine v5.x have recently adopted competitive and very appealing pricing structures for individual game developers and small teams. One may lean towards one or the other game engine based on various criteria: existing familiarization / vested interest, steepness of learning curve, quality, richness, variability and pricing of add-on assets, initial cost of basic ownership versus subscribing to updates and acquiring needed modules, etc. In this paper we present complementary features of the Unreal and Unity game engines that can be combined to enhance understanding game development and to hone relevant skills. Such combinations of different game engine features can increase the added value for our undergraduate informatics engineering students taking a "Game Technologies" class.


title: A Size-Free CLT for Poisson Multinomials and its Applications

abstract: An $(n,k)$-Poisson Multinomial Distribution (PMD) is the distribution of the sum of $n$ independent random vectors supported on the set ${\cal B}_k=\{e_1,\ldots,e_k\}$ of standard basis vectors in $\mathbb{R}^k$. We show that any $(n,k)$-PMD is ${\rm poly}\left({k\over \sigma}\right)$-close in total variation distance to the (appropriately discretized) multi-dimensional Gaussian with the same first two moments, removing the dependence on $n$ from the Central Limit Theorem of Valiant and Valiant. Interestingly, our CLT is obtained by bootstrapping the Valiant-Valiant CLT itself through the structural characterization of PMDs shown in recent work by Daskalakis, Kamath, and Tzamos. In turn, our stronger CLT can be leveraged to obtain an efficient PTAS for approximate Nash equilibria in anonymous games, significantly improving the state of the art, and matching qualitatively the running time dependence on $n$ and $1/\varepsilon$ of the best known algorithm for two-strategy anonymous games. Our new CLT also enables the construction of covers for the set of $(n,k)$-PMDs, which are proper and whose size is shown to be essentially optimal. Our cover construction combines our CLT with the Shapley-Folkman theorem and recent sparsification results for Laplacian matrices by Batson, Spielman, and Srivastava. Our cover size lower bound is based on an algebraic geometric construction. Finally, leveraging the structural properties of the Fourier spectrum of PMDs we show that these distributions can be learned from $O_k(1/\varepsilon^2)$ samples in ${\rm poly}_k(1/\varepsilon)$-time, removing the quasi-polynomial dependence of the running time on $1/\varepsilon$ from the algorithm of Daskalakis, Kamath, and Tzamos.

categories: cs.DS cs.GT cs.LG math.PR math.ST stat.TH


title: Raman resonance in iron-based superconductors: The magnetic scenario

abstract: We perform theoretical analysis of polarization-sensitive Raman spectroscopy on NaFe$_{1-x}$Co$_x$As, EuFe$2$As$_2$, SrFe$_2$As$_2$, and Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$, focusing on two features seen in the $B_{1g}$ symmetry channel (in one Fe unit cell notation): the strong temperature dependence of the static, uniform Raman response in the normal state and the existence of a collective mode in the superconducting state. We show that both features can be explained by the coupling of fermions to pairs of magnetic fluctuations via the Aslamazov-Larkin process. We first analyze magnetically-mediated Raman intensity at the leading two-loop order and then include interactions between pairs of magnetic fluctuations. We show that the full Raman intensity in the $B_{1g}$ channel can be viewed as the result of the coupling of light to Ising-nematic susceptibility via Aslamazov-Larkin process. We argue that the singular temperature dependence in the normal state is the combination of the temperature dependencies of the Aslamazov-Larkin vertex and of Ising-nematic susceptibility. We discuss two scenarios for the resonance below $T_c$. In one, the resonance is due to the development of a pole in the fully renormalized Ising-nematic susceptibility. The other is the orbital excitonic scenario, in which spin fluctuations generate an attractive interaction between low-energy fermions.


title: Unifying distillation and privileged information

abstract: Distillation (Hinton et al., 2015) and privileged information (Vapnik & Izmailov, 2015) are two techniques that enable machines to learn from other machines. This paper unifies these two techniques into generalized distillation, a framework to learn from multiple machines and data representations. We provide theoretical and causal insight about the inner workings of generalized distillation, extend it to unsupervised, semisupervised and multitask learning scenarios, and illustrate its efficacy on a variety of numerical simulations on both synthetic and real-world data.


title: Hessian and graviton propagator of the proper vertex

abstract: The proper spin-foam vertex amplitude is obtained from the EPRL vertex by projecting out all but a single gravitational sector, in order to achieve correct semi-classical behavior. In this paper we calculate the gravitational two-point function predicted by the proper spin-foam vertex to lowest order in the vertex expansion. We find the same answer as in the EPRL case in the `continuum spectrum' limit, so that the theory is consistent with the predictions of linearized gravity in the regime of small curvature. The method for calculating the two-point function is similar to that used in prior works: we cast it in terms of an action integral and to use stationary phase methods. Thus, the calculation of the Hessian matrix plays a key role. Once the Hessian is calculated, it is used not only to calculate the two-point function, but also to calculate the coefficient appearing in the semi-classical limit of the proper vertex amplitude itself. This coefficient is the effective discrete "measure factor" encoded in the spin-foam model. Through a non-trivial cancellation of different factors, we find that this coefficient is the same as the coefficient in front of the term in the asymptotics of the EPRL vertex corresponding to the selected gravitational sector.


title: Adjoint Methods for Guiding Adaptive Mesh Refinement in Wave Propagation Problems

abstract: One difficulty in developing numerical methods for hyperbolic systems of conservation laws is the fact that solutions often contain regions where much higher resolution is required than elsewhere in the domain, particularly since the solution may contain discontinuities or other localized features. The Clawpack software deals with this issue by using block-structured adaptive mesh refinement to selectively refine around propagating waves. For problems where only a target area of the total solution is of interest, a method that allows identifying and refining the grid only in regions that influence this target area would significantly reduce the computational cost of finding a solution. In this work, we show that solving the time-dependent adjoint equation and using a suitable inner product with the forward solution allows more precise refinement of the relevant waves. We present acoustics examples in one and two dimensions and a tsunami propagation example. To perform these simulations, the use of the adjoint method has been integrated into the adaptive mesh refinement strategy of the open source Clawpack and GeoClaw software. We also present results that show that the accuracy of the solution is maintained and the computational time required is significantly reduced through the integration of the adjoint method into AMR.


title: Effective field theory of dissipative fluids

abstract: We develop an effective field theory for dissipative fluids which governs the dynamics of long-lived gapless modes associated with conserved quantities. The resulting theory gives a path integral formulation of fluctuating hydrodynamics which systematically incorporates nonlinear interactions of noises. The dynamical variables are mappings between a "fluid spacetime" and the physical spacetime and an essential aspect of our formulation is to identify the appropriate symmetries in the fluid spacetime. The theory applies to nonlinear disturbances around a general density matrix. For a thermal density matrix, we require an additional $Z_2$ symmetry, to which we refer as the local KMS condition. This leads to the standard constraints of hydrodynamics, as well as a nonlinear generalization of the Onsager relations. It also leads to an emergent supersymmetry in the classical statistical regime, and a higher derivative deformation of supersymmetry in the full quantum regime.

categories: hep-th cond-mat.stat-mech gr-qc hep-ph math-ph math.MP


title: The Evolution of Stellar Rotation and the hydrogen atmospheres of habitable-zone Terrestrial Planets

abstract: Terrestrial planets formed within gaseous protoplanetary disks can accumulate significant hydrogen envelopes. The evolution of such an atmosphere due to XUV driven evaporation depends on the activity evolution of the host star, which itself depends sensitively on its rotational evolution, and therefore on its initial rotation rate. In this letter, we derive an easily applicable method for calculating planetary atmosphere evaporation that combines models for a hydrostatic lower atmosphere and a hydrodynamic upper atmosphere. We show that the initial rotation rate of the central star is of critical importance for the evolution of planetary atmospheres and can determine if a planet keeps or loses its primordial hydrogen envelope. Our results highlight the need for a detailed treatment of stellar activity evolution when studying the evolution of planetary atmospheres.


title: Solitons and vortices in nonlinear potential wells

abstract: We consider self-trapping of topological modes governed by the one- and two-dimensional (1D and 2D) nonlinear-Schrodinger/Gross-Pitaevskii equation with effective single- and double-well (DW) nonlinear potentials induced by spatial modulation of the local strength of the self-defocusing nonlinearity. This setting, which may be implemented in optics and Bose-Einstein condensates, aims to extend previous studies, which dealt with single-well nonlinear potentials. In the 1D setting, we find several types of symmetric, asymmetric and antisymmetric states, focusing on scenarios of the spontaneous symmetry breaking. The single-well model is extended by including rocking motion of the well, which gives rise to Rabi oscillations between the fundamental and dipole modes. Analysis of the 2D single-well setting gives rise to stable modes in the form of ordinary dipoles, vortex-antivortex dipoles (VADs), and vortex triangles (VTs), which may be considered as produced by spontaneous breaking of the axial symmetry. The consideration of the DW configuration in 2D reveals diverse types of modes built of components trapped in the two wells, which may be fundamental states and vortices with topological charges m = 1 and 2, as well as VADs (with m = 0) and VTs (with m = 2).


title: Nonperturbative Regulator for Chiral Gauge Theories?

abstract: We propose a nonperturbative gauge invariant regulator for d-dimensional chiral gauge theories on the lattice. The method involves simulating domain wall fermions in d + 1 dimensions with quantum gauge fields that reside on one d-dimensional surface and are extended into the bulk via gradient flow. The result is a theory of gauged fermions plus mirror fermions, where the mirror fermions couple to the gauge fields via a form factor that becomes exponentially soft with the separation between domain walls. The resultant theory has a local d-dimensional interpretation only if the chiral fermion representation is anomaly free. A physical realization of this construction would imply the existence of mirror fermions in the standard model that are invisible except for interactions induced by vacuum topology, and which could gravitate differently than conventional matter.


title: Piecewise Linear Activation Functions For More Efficient Deep Networks

abstract: This submission has been withdrawn by arXiv administrators because it is intentionally incomplete, which is in violation of our policies.


title: The Higgs Boson at LHC and the Vacuum Stability of the Standard Model

abstract: The main aim of this work is to study the conditions of absolute vacuum stability within the Standard Model (SM) by the knowledge of the behaviour of the Higgs quartic coupling up to high energy scales and using the new data on the Higgs mass given by ATLAS and CMS as an input parameter. The Higgs mass value observed by ATLAS and CMS leads to a negative value of the quartic coupling, making metastable the vacuum of the Standard Model, as it is seen by the renormalization group improved (RGI) effective potential. The stability status of SM crucially depends upon the precise values of the top and Higgs masses, a more precision determination of those masses and related uncertainties can modify drastically our conclusions about the SM stability properties. For this reason, we have computed the vacuum bubbles and the Higgs tadpole diagrams, at two-loop level of accuracy, in a renormalization scheme proposed by A. Sirlin and R. Zucchini, where the input parameters are obtained in terms of physical observables related with muon decay, and where the threshold effects are included. In particular, we focus on the analytic computation of the Higgs tadpoles contributions by two different methods. From one side we have computed the sum of the tadpoles as the first derivative of the two-loop Higgs 1PI effective potential in the Sirlin-Zucchini scheme, on the other side we have checked the obtained result toward the direct diagrammatic two- loop computation, by proposing a way of automatization of our procedure based on the generation of Feynman diagrams, reduction of their integrands and evaluation of scalar integrals and sum of all contributions in a non-redundant way. We used the code TARCER that implements the Tarasov method to reduce two-loop tensorial integrals.


title: Linear response to leadership, effective temperature and decision making in flocks

abstract: Large collections of autonomously moving agents, such as animals or micro-organisms, are able to 'flock' coherently in space even in the absence of a central control mechanism. While the direction of the flock resulting from this critical behavior is random, this can be controlled by a small subset of informed individuals acting as leaders of the group. In this article we use the Vicsek model to investigate how flocks respond to leadership and make decisions. Using a combination of numerical simulations and continuous modeling we demonstrate that flocks display a linear response to leadership that can be cast in the framework of the fluctuation-dissipation theorem, identifying an 'effective temperature' reflecting how promptly the flock reacts to the initiative of the leaders. The linear response to leadership also holds in the presence of two groups of informed individuals with competing interests, indicating that the flock's behavioral decision is determined by both the number of leaders and their degree of influence.


title: Quantum deformations of D=4 Euclidean, Lorentz, Kleinian and quaternionic o^*(4) symmetries in unified o(4;C) setting

abstract: We employ new calculational technique and present complete list of classical $r$-matrices for $D=4$ complex homogeneous orthogonal Lie algebra $\mathfrak{o}(4;\mathbb{C})$, the rotational symmetry of four-dimensional complex space-time. Further applying reality conditions we obtain the classical $r$-matrices for all possible real forms of $\mathfrak{o}(4;\mathbb{C})$: Euclidean $\mathfrak{o}(4)$, Lorentz $\mathfrak{o}(3,1)$, Kleinian $\mathfrak{o}(2,2)$ and quaternionic $\mathfrak{o}^{\star}(4)$ Lie algebras. For $\mathfrak{o}(3,1)$ we get known four classical $D=4$ Lorentz $r$-matrices, but for other real Lie algebras (Euclidean, Kleinian, quaternionic) we provide new results and mention some applications.


title: Hydrodynamical wind on vertically self-gravitating ADAFs in the presence of toroidal magnetic field

abstract: We present the effect of a hydrodynamical wind on the structure and the surface temperature of a vertically self-gravitating magnetized ADAFs using self-similar solutions. Also a model for an axisymmetric, steady-state, vertically self-gravitating hot accretion flow threaded by a toroidal magnetic field has been formulated. The model is based on $\alpha-$prescription for turbulence viscosity. It is found that the thickness and radial velocity of the disc are reduced significantly as wind gets stronger. In particular, the solutions indicated that the wind and advection have the same effects on the structure of the disc. We also find that in the optically thin ADAF becomes hotter by including the wind parameter and the self-gravity parameter.


title: Dynamical Instability of Cylindrically Symmetric Collapsing Star in Generalized Teleparallel Gravity

abstract: This paper is devoted to analyze the dynamical instability of a self-gravitating object undergoes to collapse process. We take the framework of generalized teleparallel gravity with cylindrically symmetric gravitating object. The matter distribution is represented by locally anisotropic energy-momentum tensor. We develop basic equations such as dynamical equations along with matching conditions and Harrison-Wheeler equation of state. By applying linear perturbation strategy, we construct collapse equation which is used to accomplish the instability ranges in Newtonian and post-Newtonian regimes. We find these ranges for isotropic pressure as well as reduce the results in general relativity. The unstable behavior depends on matter, metric, mass and torsion based terms.


title: A Deep Search for Prompt Radio Emission from the Short GRB 150424A With The Murchison Widefield Array

abstract: We present a search for prompt radio emission associated with the short-duration gamma-ray burst (GRB) 150424A using the Murchison Widefield Array (MWA) at frequencies from 80-133 MHz. Our observations span delays of 23 s-30 min after the GRB, corresponding to dispersion measures of 100-7700 pc/cm^3. We see no excess flux in images with timescales of 4 s, 2 min, or 30 min, and set a 3 sigma flux density limit of 3.0 Jy at 132 MHz on the shortest timescales: some of the most stringent limits to date on prompt radio emission from any type of GRB. We use these limits to constrain a number of proposed models for coherent emission from short-duration GRBs, although we show that our limits are not particularly constraining for fast radio bursts because of reduced sensitivity for this pointing. Finally, we discuss the prospects for using the MWA to search for prompt radio emission from gravitational wave transients and find that while the flux density and luminosity limits are likely to be very constraining, the latency of the gravitational wave alert may limit the robustness of any conclusions.


title: Sums of Random Matrices and the Potts Model on Random Planar Maps

abstract: We compute the partition function of the $q$-states Potts model on a random planar lattice with $p\leq q$ allowed, equally weighted colours on a connected boundary. To this end, we employ its matrix model representation in the planar limit, generalising a result by Voiculescu for the addition of random matrices to a situation beyond free probability theory. We show that the partition functions with $p$ and $q-p$ colours on the boundary are related algebraically. Finally, we investigate the phase diagram of the model when $0\leq q\leq 4$ and comment on the conformal field theory description of the critical points.

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


title: Around Sylvester's question in the plane

abstract: Pick $n$ points $Z_0,...,Z_{n-1}$ uniformly and independently at random in a compact convex set $H$ with non empty interior of the plane, and let $Q^n_H$ be the probability that the $Z_i$'s are the vertices of a convex polygon. Blaschke 1917 \cite{Bla} proved that $Q^4_T\leq Q^4_H\leq Q^4_D$, where $D$ is a disk and $T$ a triangle. In the present paper we prove $Q^5_T\leq Q^5_H\leq Q^5_D$. One of the main ingredients of our approach is a new formula for $Q^n_H$ which permits to prove that Steiner symmetrization does not decrease $Q^5_H$, and that shaking does not increases it (this is the method Blaschke used in the $n=4$ case). We conjecture that the new formula we provide will lead in the future to the complete proof that $Q^n_T\leq Q^n_H\leq Q^n_D$ , for any $n$.


title: He II Proximity Effect and the Lifetime of Quasars

abstract: The lifetime of quasars is fundamental for understanding the growth of supermassive black holes, and is an important ingredient in models of the reionization of the intergalactic medium. However, despite various attempts to determine quasar lifetimes, current estimates from a variety of methods are uncertain by orders of magnitude. This work combines cosmological hydrodynamical simulations and 1D radiative transfer to investigate the structure and evolution of the He II Ly$\alpha$ proximity zones around quasars at $z \simeq 3-4$. We show that the time evolution in the proximity zone can be described by a simple analytical model for the approach of the He II fraction $x_{\rm HeII}\left( t \right)$ to ionization equilibrium, and use this picture to illustrate how the transmission profile depends on the quasar lifetime, quasar UV luminosity, and the ionization state of helium in the ambient IGM (i.e. the average He II fraction, or equivalently the metagalactic He II ionizing background). A significant degeneracy exists between the lifetime and the average He II fraction, however the latter can be determined from measurements of the He II Ly$\alpha$ optical depth far from quasars, allowing the lifetime to be measured. We advocate stacking existing He II quasar spectra at $z\sim 3$, and show that the shape of this average proximity zone profile is sensitive to lifetimes as long as $\sim 30$ Myr. At higher redshift $z\sim 4$ where the He II fraction is poorly constrained, degeneracies will make it challenging to determine these parameters independently. Our analytical model for He II proximity zones should also provide a useful description of the properties of H I proximity zones around quasars at $z \simeq 6-7$.


title: The Phase Coherence of Interstellar Density Fluctuations

abstract: Studies of MHD turbulence often investigate the Fourier power spectrum to provide information on the nature of the turbulence cascade. However, the Fourier power spectrum only contains the Fourier amplitudes and rejects all information regarding the Fourier phases. Here we investigate the utility of two statistical diagnostics for recovering information on Fourier phases in ISM density data: the averaged amplitudes of the bispectrum and the phase coherence index (PCI), a new phase technique for the ISM. We create 3D density and 2D column density maps using a set of simulations of isothermal ideal MHD turbulence with a wide range of sonic and Alfv\'enic Mach numbers. We find that the bispectrum averaged along different angles with respect to either the $k_1$ or $k_2$ axis is primarily sensitive to the sonic Mach number while averaging the bispectral amplitudes over different annuli is sensitive to both the sonic and Alfv\'enic Mach numbers. The PCI of density suggests that the most correlated phases occur in supersonic sub-Alfv\'enic turbulence and also near the numerical dissipation regime. This suggests that non-linear interactions with correlated phases are strongest in shock dominated regions, in agreement with findings from the solar wind. Additionally, our results are particularly encouraging as they suggests the phase information contained in the bispectrum and PCI can be used to find parameters of turbulence in column density maps.


title: Auroral radio emission from late L and T dwarfs: A new constraint on dynamo theory in the substellar regime

abstract: We have observed 6 late-L and T dwarfs with the Karl G. Jansky Very Large Array (VLA) to investigate the presence of highly circularly polarized radio emission, associated with large-scale auroral currents. Previous surveys encompassing ~60 L6 or later targets in this spectral range have yielded only one detection. Our sample includes the previously detected T6.5 dwarf 2MASS 10475385+2124234 as well as 5 new targets selected for the presence of H-alpha emission or optical/infrared photometric variability, which are possible manifestations of auroral activity. We detect 2MASS 10475385+2124234, as well as 4 of the 5 targets in our biased sample, including the strong IR variable SIMP J01365662+0933473 and bright H-alpha emitter 2MASS 12373919+6526148, reinforcing the possibility that activity at these disparate wavelengths is related. The radio emission frequency corresponds to a precise determination of the lower-bound magnetic field strength near the surface of each dwarf and this new sample provides robust constraints on dynamo theory in the low mass brown dwarf regime. Magnetic fields >2.5 kG are confirmed for 5/6 targets. Our results provide tentative evidence that the dynamo operating in this mass regime may be inconsistent with predicted values from a recently proposed model. Further observations at higher radio frequencies are essential for verifying this assertion.


title: The mass-sheet degeneracy and time-delay cosmography: Analysis of the strong lens RXJ1131-1231

abstract: We present extended modeling of the strong lens system RXJ1131-1231 with archival data in two HST bands in combination with existing line-of-sight contribution and velocity dispersion estimates. Our focus is on source size and its influence on time-delay cosmography. We therefore examine the impact of mass-sheet degeneracy and especially the degeneracy pointed out by Schneider & Sluse (2013) using the source reconstruction scale. We also extend on previous work by further exploring the effects of priors on the kinematics of the lens and the external convergence in the environment of the lensing system. Our results coming from RXJ1131-1231 are given in a simple analytic form so that they can be easily combined with constraints coming from other cosmological probes. We find that the choice of priors on lens model parameters and source size are subdominant for the statistical errors for $H_0$ measurements of this systems. The choice of prior for the source is sub-dominant at present (2% uncertainty on $H_0$) but may be relevant for future studies. More importantly, we find that the priors on the kinematic anisotropy of the lens galaxy have a significant impact on our cosmological inference. When incorporating all the above modeling uncertainties, we find $H_0 = 86.6^{+6.8}_{-6.9}$ km s$^{-1}$Mpc$^{-1}$, when using kinematic priors similar to other studies. When we use a different kinematic prior motivated by Barnab\`e et al. (2012) but covering the same anisotropic range, we find $H_0 = 74.5^{+8.0}_{-7.8}$ km s$^{-1}$Mpc$^{-1}$. This means that the choice of kinematic modeling and priors have a significant impact on cosmographic inferences. The way forward is either to get better velocity dispersion measures which would down weight the impact of the priors or to construct physically motivated priors for the velocity dispersion model.


title: The Frequency of Giant impacts on Earth-like Worlds

abstract: The late stages of terrestrial planet formation are dominated by giant impacts that collectively influence the growth, composition and habitability of any planets that form. Hitherto, numerical models designed to explore these late stage collisions have been limited by assuming that all collisions lead to perfect accretion, and many of these studies lack the large number of realizations needed to account for the chaotic nature of N-body systems. We improve on these limitations by performing 280 simulations of planet formation around a Sun-like star, half of which used an N-body algorithm that has recently been modified to include fragmentation and hit-and-run (bouncing) collisions. We find that when fragmentation is included, the final planets formed are comparable in terms of mass and number, however their collision histories differ significantly and the accretion time approximately doubles. We explored impacts onto Earth-like planets which we parameterized in terms of their specific impact energies. Only 15 of our 164 Earth-analogs experienced an impact that was energetic enough to strip an entire atmosphere. To strip about half of an atmosphere requires energies comparable to recent models of the Moon-forming giant impact. Almost all Earth-analogs received at least one impact that met this criteria during the 2 Gyr simulations and the median was three giant impacts. The median time of the final giant impact was 43 Myr after the start of the simulations, leading us to conclude that the time-frame of the Moon-forming impact is typical amongst planetary systems around Sun-like stars.


title: The capabilities and performance of the Automated Planet Finder Telescope with the implementation of a dynamic scheduler

abstract: We report initial performance results emerging from 600 hours of observations with the Automated Planet Finder (APF) telescope and Levy Spectrometer located at UCO/Lick Observatory. We have obtained multiple spectra of 80 G, K and M-type stars, which comprise 4,954 individual Doppler radial velocity (RV) measurements with a median internal uncertainty of 1.35 ms$^{-1}$. We find a strong, expected correlation between the number of photons accumulated in the 5000-6200$\AA$ iodine region of the spectrum, and the resulting internal uncertainty estimates. Additionally, we find an offset between the population of G and K stars and the M stars within the data set when comparing these parameters. As a consequence of their increased spectral line densities, M-type stars permit the same level of internal uncertainty with 2x fewer photons than G-type and K-type stars. When observing M stars, we show that the APF/Levy has essentially the same speed-on-sky as Keck/HIRES for precision RVs. In the interest of using the APF for long-duration RV surveys, we have designed and implemented a dynamic scheduling algorithm. We discuss the operation of the scheduler, which monitors ambient conditions and combines on-sky information with a database of survey targets to make intelligent, real-time targeting decisions.


title: Bekenstein's generalized second law of thermodynamics: The role of the hoop conjecture

abstract: Bekenstein's generalized second law (GSL) of thermodynamics asserts that the sum of black-hole entropy, $S_{\text{BH}}=Ac^3/4\hbar G$ (here $A$ is the black-hole surface area), and the ordinary entropy of matter and radiation fields in the black-hole exterior region never decreases. We here re-analyze an intriguing gedanken experiment which was designed by Bekenstein to challenge the GSL. In this historical gedanken experiment an entropy-bearing box is lowered into a charged Reissner-Nordstr\"om black hole. For the GSL to work, the resulting increase in the black-hole surface area (entropy) must compensate for the loss of the box's entropy. We show that if the box can be lowered adiabatically all the way down to the black-hole horizon, as previously assumed in the literature, then for near-extremal black holes the resulting increase in black-hole surface-area (due to the assimilation of the box by the black hole) may become too small to compensate for the loss of the box's entropy. In order to resolve this apparent violation of the GSL, we here suggest to use a generalized version of the hoop conjecture. In particular, assuming that a physical system of mass $M$ and electric charge $Q$ forms a black hole if its circumference radius $r_{\text{c}}$ is equal to (or smaller than) the corresponding Reissner-Nordstr\"om black-hole radius $r_{\text{RN}}=M+\sqrt{M^2-Q^2}$, we prove that a new (and larger) horizon is already formed before the entropy-bearing box reaches the horizon of the original near-extremal black hole. This result, which seems to have been overlooked in previous analyzes of the composed black-hole-box system, ensures the validity of Bekenstein's GSL in this famous gedanken experiment.


title: Entanglement Entropy in a Holographic Kondo Model

abstract: We calculate entanglement and impurity entropies in a recent holographic model of a magnetic impurity interacting with a strongly coupled system. There is an RG flow to an IR fixed point where the impurity is screened, leading to a decrease in impurity degrees of freedom. This information loss corresponds to a volume decrease in our dual gravity model, which consists of a codimension one hypersurface embedded in a BTZ black hole background in three dimensions. There are matter fields defined on this hypersurface which are dual to Kondo field theory operators. In the large N limit, the formation of the Kondo cloud corresponds to the condensation of a scalar field. The entropy is calculated according to the Ryu-Takayanagi prescription. This requires to determine the backreaction of the hypersurface on the BTZ geometry, which is achieved by solving the Israel junction conditions. We find that the larger the scalar condensate gets, the more the volume of constant time slices in the bulk is reduced, shortening the bulk geodesics and reducing the impurity entropy. This provides a new non-trivial example of an RG flow satisfying the g-theorem. Moreover, we find explicit expressions for the impurity entropy which are in agreement with previous field theory results for free electrons. This demonstrates the universality of perturbing about an IR fixed point.


title: Stellar streams around the Magellanic Clouds

abstract: Using Blue Horizontal Branch stars identified in the Dark Energy Survey Year 1 data, we report the detection of an extended and lumpy stellar debris distribution around the Magellanic Clouds. At the heliocentric distance of the Clouds, overdensities of BHBs are seen to reach at least to ~30 degrees, and perhaps as far as ~50 degrees from the LMC. In 3D, the stellar halo is traceable to between 25 and 50 kpc from the LMC. We catalogue the most significant of the stellar sub-structures revealed, and announce the discovery of a number of narrow streams and diffuse debris clouds. Two narrow streams appear approximately aligned with the Magellanic Clouds' proper motion. Moreover, one of these overlaps with the gaseous Magellanic Stream on the sky. Curiously, two diffuse BHB agglomerations seem coincident with several of the recently discovered DES satellites. Given the enormous size and the conspicuous lumpiness of the LMC's stellar halo, we speculate that the dwarf could easily have been more massive than previously had been assumed.


title: Entanglement generation in periodically driven integrable systems: dynamical phase transitions and steady state

abstract: We study a class of periodically driven $d-$dimensional integrable models and show that after $n$ drive cycles with frequency $\omega$, pure states with non-area-law entanglement entropy $S_n(l) \sim l^{\alpha(n,\omega)}$ are generated, where $l$ is the linear dimension of the subsystem, and $d-1 \le \alpha(n,\omega) \le d$. We identify and analyze the crossover phenomenon from an area ($S \sim l^{ d-1}$ for $d\geq1$) to a volume ($S \sim l^{d}$) law and provide a criterion for their occurrence which constitutes a generalization of Hastings' theorem to driven integrable systems in one dimension. We also find that $S_n$ generically decays to $S_{\infty}$ as $(\omega/n)^{(d+2)/2}$ for fast and $(\omega/n)^{d/2}$ for slow periodic drives; these two dynamical phases are separated by a topological transition in the eigensprectrum of the Floquet Hamiltonian. This dynamical transition manifests itself in the temporal behavior of all local correlation functions and does not require a critical point crossing during the drive. We find that these dynamical phases show a rich re-entrant behavior as a function of $\omega$ for $d=1$ models, and also discuss the dynamical transition for $d>1$ models. Finally, we study entanglement properties of the steady state and show that singular features (cusps and kinks in $d=1$) appear in $S_{\infty}$ as a function of $\omega$ whenever there is a crossing of the Floquet bands. We discuss experiments which can test our theory.


title: Non-BPS multi-bubble microstate geometries

abstract: We construct the first smooth horizonless supergravity solutions that have two topologically-nontrivial three-cycles supported by flux, and that have the same mass and charges as a non-extremal D1-D5-P black hole. Our configurations are solutions to six-dimensional ungauged supergravity coupled to a tensor multiplet, and uplift to solutions of Type IIB supergravity. The solutions represent multi-center generalizations of the non-BPS solutions of Jejjala, Madden, Ross, and Titchener, which have over-rotating angular momenta. By adding an additional Gibbons-Hawking center, we succeed in lowering one of the two angular momenta below the cosmic censorship bound, and bringing the other very close to this bound. Our results demonstrate that it is possible to construct multi-center horizonless solutions corresponding to non-extremal black holes, and offer the prospect of ultimately establishing that finite-temperature black holes have nontrivial structure at the horizon.


title: Gravity as the main driver of non-thermal motions in massive star formation

abstract: The origin of non-thermal motions in massive star forming regions can be ascribed to turbulence acting against the gravitational collapse, or to the self-gravity itself driving the rapid global collapse. The dependence between velocity dispersion, radius and clouds surface density found by Heyer et al. (2009), $\sigma/R^{1/2}\propto \Sigma^{1/2}$, has been interpreted in terms of global collapse of clouds. In this work we demonstrate that this relation is an expression of a more general relation between accelerations. We introduce the gravo-turbulent acceleration, a$_k$, which describe the non-thermal motions in each region, and the acceleration generated by the gravitational field a$_G$, which is proportional to $\Sigma$. We also introduce a new coefficient, the force partition coefficient $\alpha_{for}$ which is equivalent to the virial parameter but does not distinguish between collapsing and non-collapsing regions. In this work we use the a$_k$ - a$_G$ formalism in the analysis of a new sample of 16 massive starless clumps (MSCls) combined with data from the literature. We show that a$_k$ and a$_G$ are not independent. The non-thermal motions in each region can originate from both local turbulence and self-gravity but overall the data in the a$_k$ vs. a$_G$ diagram demonstrate that the majority of the non-thermal motions originate from self-gravity. We further show that all the MSCls with $\Sigma\geq 0.1$ g cm$^{-2}$ show signs of infall motions, a strong indication that the denser regions are the first to collapse. Finally, we include in the formalism the contribution of an external pressure and the magnetic fields.


title: Strange and Charm Quark Spins from Anomalous Ward Identity

abstract: We present a calculation of the strange and charm quark contributions to the nucleon spin from the anomalous Ward identity (AWI). It is performed with overlap valence quarks on 2+1-flavor domain-wall fermion gauge configurations on a $24^3 \times 64$ lattice with the light sea mass at $m_{\pi} = 330$ MeV. To satisfy the AWI, the overlap fermion for the pseudoscalar density and the overlap Dirac operator for the topological density, which do not have multiplicative renormalization, are used to normalize the form factor of the local axial-vector current at finite $q^2$. For the charm quark, we find that the negative pseudoscalar term almost cancels the positive topological term. For the strange quark, the pseudoscalar term is less negative than that of the charm. By imposing the AWI, the strange $g_A(q^2)$ at $q^2 =0$ is obtained by a global fit of the pseudoscalar and the topological form factors, together with $g_A(q^2)$ and the induced pseudoscalar form factor $h_A(q^2)$ at finite $q^2$. The chiral extrapolation to the physical pion mass gives $\Delta s + \Delta {\bar{s}} = -0.0403(44)(78)$.


title: Estimating SI violation in CMB due to non-circular beam and complex scan in minutes

abstract: Mild, unavoidable deviations from circular-symmetry of instrumental beams along with scan strategy can give rise to measurable Statistical Isotropy (SI) violation in Cosmic Microwave Background (CMB) experiments. If not accounted properly, this spurious signal can complicate the extraction of other SI violation signals (if any) in the data. However, estimation of this effect through exact numerical simulation is computationally intensive and time consuming. A generalized analytical formalism not only provides a quick way of estimating this signal, but also gives a detailed understanding connecting the leading beam anisotropy components to a measurable BipoSH characterisation of SI violation. In this paper, we provide an approximate generic analytical method for estimating the SI violation generated due to a non-circular (NC) beam and arbitrary scan strategy, in terms of the Bipolar Spherical Harmonic (BipoSH) spectra. Our analytical method can predict almost all the features introduced by a NC beam in a complex scan and thus reduces the need for extensive numerical simulation worth tens of thousands of CPU hours into minutes long calculations. As an illustrative example, we use WMAP beams and scanning strategy to demonstrate the easability, usability and efficiency of our method. We test all our analytical results against that from exact numerical simulations.


title: PICS: Simulations of Strong Gravitational Lensing in Galaxy Clusters

abstract: Gravitational lensing has become one of the most powerful tools available for investigating the 'dark side' of the universe. Cosmological strong gravitational lensing, in particular, probes the properties of the dense cores of dark matter halos over decades in mass and offers the opportunity to study the distant universe at flux levels and spatial resolutions otherwise unavailable. Studies of strongly-lensed variable sources offer yet further scientific opportunities. One of the challenges in realizing the potential of strong lensing is to understand the statistical context of both the individual systems that receive extensive follow-up study, as well as that of the larger samples of strong lenses that are now emerging from survey efforts. Motivated by these challenges, we have developed an image-simulation pipeline, PICS (Pipeline for Images of Cosmological Strong lensing) to generate realistic strong gravitational lensing signals from group and cluster scale lenses. PICS uses a low-noise and unbiased density estimator based on (resampled) Delaunay Tessellations to calculate the density field; lensed images are produced by ray-tracing images of actual galaxies from deep Hubble Space Telescope observations. Other galaxies, similarly sampled, are added to fill in the light cone. The pipeline further adds cluster-member galaxies and foreground stars into the lensed images. The entire image ensemble is then observed using a realistic point spread function which includes appropriate detector artifacts for bright stars. Noise is further added, including such non-Gaussian elements as noise window-paning from mosaiced observations, residual bad pixels, and cosmic rays. The aim is to produced simulated images that appear identical - to the eye (expert or otherwise) - to real observations in various imaging surveys.


title: Strong tW Scattering at the LHC

abstract: Deviations of the top electroweak couplings from their Standard Model values imply that certain amplitudes for the scattering of third generation fermions and longitudinally polarized vector bosons or Higgses diverge quadratically with momenta. This high-energy growth is a genuine signal of models where the top quark is strongly coupled to the sector responsible for electroweak symmetry breaking. We propose to profit from the high energies accessible at the LHC to enhance the sensitivity to non-standard top-$Z$ couplings, which are currently very weakly constrained. To demonstrate the effectiveness of the approach, we perform a detailed analysis of $tW \to tW$ scattering, which can be probed at the LHC via $pp\to t\bar{t}Wj$. By recasting a CMS analysis at 8 TeV, we derive the strongest direct bounds to date on the $Ztt$ couplings. We also design a dedicated search at 13 TeV that exploits the distinctive features of the $t\bar{t}Wj$ signal. Finally, we present other scattering processes in the same class that could provide further tests of the top-Higgs sector.


title: Membership in moment polytopes is in NP and coNP

abstract: We show that the problem of deciding membership in the moment polytope associated with a finite-dimensional unitary representation of a compact, connected Lie group is in NP and coNP. This is the first non-trivial result on the computational complexity of this problem, which naively amounts to a quadratically-constrained program. Our result applies in particular to the Kronecker polytopes, and therefore to the problem of deciding positivity of the stretched Kronecker coefficients. In contrast, it has recently been shown that deciding positivity of a single Kronecker coefficient is NP-hard in general [Ikenmeyer, Mulmuley and Walter, arXiv:1507.02955]. We discuss the consequences of our work in the context of complexity theory and the quantum marginal problem.

categories: cs.CC math-ph math.MP math.RT math.SG quant-ph


title: Variable-range hopping through marginally localized phonons

abstract: We investigate the effect of coupling Anderson localized particles in one dimension to a system of marginally localized phonons having a symmetry protected delocalized mode at zero frequency. This situation is naturally realized for electrons coupled to phonons in a disordered nano-wire as well as for ultra-cold fermions coupled to phonons of a superfluid in a one dimensional disordered trap. To determine if the coupled system can be many-body localized we analyze the phonon-mediated hopping transport for both the weak and strong coupling regimes. We show that the usual variable-range hopping mechanism involving a low-order phonon processes is ineffective at low temperature due to discreteness of the bath at the required energy. Instead, the system thermalizes through a many-body process involving exchange of a diverging number $n\propto -\log T$ of phonons in the low temperature limit. This effect leads to a highly singular prefactor to Mott's well known formula and strongly suppresses the variable range hopping rate. Finally we comment on possible implications of this physics in higher dimensional electron-phonon coupled systems.


title: Learning to Diagnose with LSTM Recurrent Neural Networks

abstract: Clinical medical data, especially in the intensive care unit (ICU), consist of multivariate time series of observations. For each patient visit (or episode), sensor data and lab test results are recorded in the patient's Electronic Health Record (EHR). While potentially containing a wealth of insights, the data is difficult to mine effectively, owing to varying length, irregular sampling and missing data. Recurrent Neural Networks (RNNs), particularly those using Long Short-Term Memory (LSTM) hidden units, are powerful and increasingly popular models for learning from sequence data. They effectively model varying length sequences and capture long range dependencies. We present the first study to empirically evaluate the ability of LSTMs to recognize patterns in multivariate time series of clinical measurements. Specifically, we consider multilabel classification of diagnoses, training a model to classify 128 diagnoses given 13 frequently but irregularly sampled clinical measurements. First, we establish the effectiveness of a simple LSTM network for modeling clinical data. Then we demonstrate a straightforward and effective training strategy in which we replicate targets at each sequence step. Trained only on raw time series, our models outperform several strong baselines, including a multilayer perceptron trained on hand-engineered features.


title: Abelian Girth and Girth

abstract: We show that the abelian girth of a graph is at least three times its girth. We prove an analogue of the Moore bound for the abelian girth of regular graphs, where the degree of the graph is fixed and the number of vertices is large. We conclude that one could try to improve the Moore bound for graphs of fixed degree and many vertices by trying to improve its analogue concerning the abelian girth.


title: Invariance of the generalized oscillator under linear transformation of the related system of orthogonal polynomials

abstract: We consider two families of polynomials $\mathbb{P}=\polP$ and $\mathbb{Q}=\polQ$\footnote{Here and below we consider only monic polynomials.} orthogonal on the real line with respect to probability measures $\mu$ and $\nu$ respectively. Let $\polQ$ and $\polP$ connected by the linear relations $$ Q_n(x)=P_n(x)+a_1P_{n-1}(x)+...+a_kP_{n-k}(x).$$ Let us denote $\mathfrak{A}_P$ and $\mathfrak{A}_Q$ generalized oscillator algebras associated with the sequences $\mathbb{P}$ and $\mathbb{Q}$. In the case $k=2$ we describe all pairs ($\mathbb{P}$,$\mathbb{Q}$), for which the algebras $\mathfrak{A}_P$ and $\mathfrak{A}_Q$ are equal. In addition, we construct corresponding algebras of generalized oscillators for arbitrary $k\geq1$.


title: Cavity magnomechanics

abstract: A dielectric body couples with electromagnetic fields through radiation pressure and electrostrictive forces, which mediate phonon-photon coupling in cavity optomechanics. In a magnetic medium, according to Korteweg-Helmholtz formula, magnetostrictive forces should arise and lead to phonon-magnon interaction. Here we report such a coupled phonon-magnon system based on ferrimagnetic spheres, which we term as cavity magnomechanics, by analogy to cavity optomechanics. Coherent phonon-magnon interactions, including electromagnetically induced transparency and absorption, are demonstrated. Excitingly, due to strong hybridization of magnon and microwave photon modes and their high tunability, our platform exhibits new features including parametric amplification of magnons and phonons, triply resonant photon-magnon-phonon coupling and phonon lasing. Our work demonstrates the fundamental principle of cavity magnomechanics and its application as a new information transduction platform based on coherent coupling between photons, phonons and magnons.


title: Einstein-Vector Gravity, Emerging Gauge Symmetry and de Sitter Bounce

abstract: We construct a class of Einstein-vector theories where the vector field couples bilinearly to the curvature polynomials of arbitrary order in such a way that only Riemann tensor rather than its derivative enters the equations of motion. The theories can thus be ghost free. The U(1) gauge symmetry may emerge in the vacuum and also in some weak-field limit. We focus on the two-derivative theory and study a variety of applications. We find that in this theory, the energy-momentum tensor of dark matter provides a position-dependent gauge-violating term to the Maxwell field. We also use the vector as an inflaton and construct cosmological solutions. We find that the expansion can accelerate without a bared cosmological constant, indicating a new candidate for dark energy. Furthermore we obtain exact solutions of de Sitter bounce, generated by the vector which behaves like a Maxwell field in the later time. We also obtain a few new exact black holes that are asymptotic to flat and Lifshitz spacetimes. In addition, we construct exact wormholes, and Randall-Sundrum II domain walls.


title: Scrutinizing the Alignment Limit in Two-Higgs-Doublet Models. Part 2: $m_H=125$ GeV

abstract: In the alignment limit of a multi-doublet Higgs sector, one of the Higgs mass eigenstates aligns in field space with the direction of the scalar field vacuum expectation values, and its couplings approach those of the Standard Model (SM) Higgs boson. We consider CP-conserving Two-Higgs-Doublet Models (2HDMs) of Type I and Type II near the alignment limit in which the heavier of the two CP-even Higgs bosons, $H$, is the SM-like state observed with a mass of 125 GeV, and the couplings of $H$ to gauge bosons approach those of the SM. We review the theoretical structure and analyze the phenomenological implications of this particular realization of the alignment limit, where decoupling of the extra states cannot occur given that the lighter CP-even state $h$ must, by definition, have a mass below 125 GeV. For the numerical analysis, we perform scans of the 2HDM parameter space employing the software packages 2HDMC and Lilith, taking into account all relevant pre-LHC constraints, constraints from the measurements of the 125 GeV Higgs signal at the LHC, as well as the most recent limits coming from searches for heavy Higgs-like states. Implications for Run 2 at the LHC, including expectations for observing the other scalar states, are also discussed.


title: Generative Concatenative Nets Jointly Learn to Write and Classify Reviews

abstract: A recommender system's basic task is to estimate how users will respond to unseen items. This is typically modeled in terms of how a user might rate a product, but here we aim to extend such approaches to model how a user would write about the product. To do so, we design a character-level Recurrent Neural Network (RNN) that generates personalized product reviews. The network convincingly learns styles and opinions of nearly 1000 distinct authors, using a large corpus of reviews from BeerAdvocate.com. It also tailors reviews to describe specific items, categories, and star ratings. Using a simple input replication strategy, the Generative Concatenative Network (GCN) preserves the signal of static auxiliary inputs across wide sequence intervals. Without any additional training, the generative model can classify reviews, identifying the author of the review, the product category, and the sentiment (rating), with remarkable accuracy. Our evaluation shows the GCN captures complex dynamics in text, such as the effect of negation, misspellings, slang, and large vocabularies gracefully absent any machinery explicitly dedicated to the purpose.


title: Quantum Fingerprints in Higher Order Correlators of a Spin Qubit

abstract: The spin of an electron in a semiconductor quantum dot represents a natural nanoscale solid state qubit. Coupling to nuclear spins leads to decoherence that limits the number of allowed quantum logic operations for this qubit. Traditional approach to characterize decoherence is to explore spin Relaxation and the spin echo, which are equivalent to the studies of the spins 2nd order time-correlator at various external conditions. Here we develop an alternative technique by showing that considerable information about dynamics can be obtained from direct measurements of higher than the 2nd order correlators, which to date have been hindered in semiconductor quantum dots. We show that such correlators are sensitive to pure quantum effects that cannot be explained within the classical framework, and which allow direct determination of ensemble and quantum dephasing times with only repeated projective measurements without the need for coherent spin control. Our method enables studies of pure quantum effects, including tests of the Leggett-Garg type inequalities that rule out local hidden variable interpretation of the spin qubit dynamics.


title: Dynamics of wave fronts and filaments in anisotropic cardiac tissue

abstract: The heartbeat is mediated between cardiac cells by waves of electrical depolarisation. During cardiac arrhythmias, electrical activity was found to be organised in scroll waves which rotate around a dynamical filament curve. In this thesis, a curved-space approach is used to mathematically capture anisotropy of wave propagation. We derive for the first time the covariant laws of motion for traveling wave fronts and scroll wave filaments in anisotropic excitable media such as cardiac tissue. We show that locally varying anisotropy yields non-zero Riemann tensor components, which may alter the stability of scroll wave filaments. The instability of scroll wave filaments has been linked to transition from ventricular tachycardia to fibrillation.


title: On the Capacity of Block Multiantenna Channels

abstract: In this paper we consider point-to-point multiantenna channels with certain block distributional symmetries which do not require the entries of the channel matrix to be either Gaussian, or independent, or identically distributed. A main contribution is a capacity theorem for these channels, which might be regarded as a generalization of Telatar's theorem (1999), which reduces the numerical optimization domain in the capacity computation. With this information theoretic result and some free probability arguments, we prove an asymptotic capacity theorem that, in addition to reducing the optimization domain, does not depend on the dimension of the channel matrix. This theorem allows us to apply free probability techniques to numerically compute the asymptotic capacity of the channels under consideration. These theorems provide a very efficient method for numerically approximating both the capacity and a capacity achieving input covariance matrix of certain channels.


title: Variational Analysis of Convexly Generated Spectral Max Functions

abstract: The spectral abscissa is the largest real part of an eigenvalue of a matrix and the spectral radius is the largest modulus. Both are examples of spectral max functions---the maximum of a real-valued function over the spectrum of a matrix. These mappings arise in the control and stabilization of dynamical systems. In 2001, Burke and Overton characterized the regular subdifferential of the spectral abscissa and showed that the spectral abscissa is subdifferentially regular in the sense of Clarke when all active eigenvalues are nonderogatory. In this paper we develop new techniques to obtain these results for the more general class of convexly generated spectral max functions. In particular, we extend the Burke-Overton subdifferential regularity result to this class. These techniques allow us to obtain new variational results for the spectral radius.


title: Online Principal Component Analysis in High Dimension: Which Algorithm to Choose?

abstract: In the current context of data explosion, online techniques that do not require storing all data in memory are indispensable to routinely perform tasks like principal component analysis (PCA). Recursive algorithms that update the PCA with each new observation have been studied in various fields of research and found wide applications in industrial monitoring, computer vision, astronomy, and latent semantic indexing, among others. This work provides guidance for selecting an online PCA algorithm in practice. We present the main approaches to online PCA, namely, perturbation techniques, incremental methods, and stochastic optimization, and compare their statistical accuracy, computation time, and memory requirements using artificial and real data. Extensions to missing data and to functional data are discussed. All studied algorithms are available in the R package onlinePCA on CRAN.


title: Spatially resolved eastward winds and rotation of HD 189733b

abstract: We measure wind velocities on opposite sides of the hot Jupiter HD$\,$189733b by modeling sodium absorption in high-resolution HARPS transmission spectra. Our model implicitly accounts for the Rossiter-McLaughlin effect, which we show can explain the high wind velocities suggested by previous studies. Our results reveal a strong eastward motion of the atmosphere of HD$\,$189733b, with a redshift of $2.3^{+1.3}_{-1.5}$$\,$km$\,$s$^{-1}$ on the leading limb of the planet and a blueshift of $5.3^{+1.0}_{-1.4}$$\,$km$\,$s$^{-1}$ on the trailing limb. These velocities can be understood as a combination of tidally locked planetary rotation and an eastward equatorial jet; closely matching the predictions of atmospheric circulation models. Our results show that the sodium absorption of HD$\,$189733b is intrinsically velocity broadened and so previous studies of the average transmission spectrum are likely to have overestimated the role of pressure and thermal broadening.


title: Deep Multimodal Semantic Embeddings for Speech and Images

abstract: In this paper, we present a model which takes as input a corpus of images with relevant spoken captions and finds a correspondence between the two modalities. We employ a pair of convolutional neural networks to model visual objects and speech signals at the word level, and tie the networks together with an embedding and alignment model which learns a joint semantic space over both modalities. We evaluate our model using image search and annotation tasks on the Flickr8k dataset, which we augmented by collecting a corpus of 40,000 spoken captions using Amazon Mechanical Turk.


title: Polarization of Magnetic Dipole Emission and Spinning Dust Emission from Magnetic Nanoparticles

abstract: Magnetic dipole emission (MDE) from interstellar magnetic nanoparticles is an important Galactic foreground in the microwave frequencies, and its polarization level may pose great challenges for achieving reliable measurements of cosmic microwave background (CMB) B-mode signal. To obtain theoretical constraints on the polarization of MDE, we first compute the degree of alignment of big silicate grains incorporated with magnetic inclusions. We find that, in realistic conditions of the interstellar medium, thermally rotating big grains with magnetic inclusions are weakly aligned and achieve {\it alignment saturation} when the magnetic alignment rate becomes much faster than the rotational damping rate. We then compute the degree of alignment for free-flying magnetic nanoparticles, taking into account various interaction processes of grains with the ambient gas and radiation field, including neutral collisions, ion collisions, and infrared emission. We find that the rotational damping by infrared emission can significantly decrease the degree of alignment of small particles from the saturation level, whereas the excitation by ion collisions can enhance the alignment of ultrasmall particles. Using the computed degrees of alignment, we predict the polarization level of MDE from free-flying magnetic nanoparticles to be rather low. Such a polarization level is within the upper limits measured for anomalous microwave emission (AME), which indicates that MDE from free-flying iron particles may not be ruled out as a source of AME. We also quantify spinning dust emission from free-flying iron nanoparticles with permanent magnetic moments and find that its emissivity is one order of magnitude lower than that from spinning polycyclic aromatic hydrocarbons (PAHs). Finally, we compute the polarization spectra of spinning dust emission from PAHs for the different interstellar magnetic fields.


title: Optimal bounds for convergence of expected spectral distributions to the semi-circular law for the $4+\epsilon$ moment ensemble

abstract: This paper extends a previous bound of order $O(n^{-1})$ of the authors (arXiv:1405.7820[math.PR]), for the rate of convergence in Kolmogorov distance of the expected spectral distribution of a Wigner random matrix ensemble to the semicircular law. Here we relax the moment conditions for entries of the Wigner matrices from order $8$ to order $4+ \epsilon$ for an arbitrary small $\epsilon>0$.


title: Game characterizations and lower cones in the Weihrauch degrees

abstract: We introduce a parametrized version of the Wadge game for functions and show that each lower cone in the Weihrauch degrees is characterized by such a game. These parametrized Wadge games subsume the original Wadge game, the eraser and backtrack games as well as Semmes's tree games. In particular, we propose that the lower cones in the Weihrauch degrees are the answer to Andretta's question on which classes of functions admit game characterizations. We then discuss some applications of such parametrized Wadge games. Using machinery from Weihrauch reducibility theory, we introduce games characterizing every (transfinite) level of the Baire hierarchy via an iteration of a pruning derivative on countably branching trees.


title: Examination of the scientific review process: Ten best practice suggestions for an improved process

abstract: In this article we wish to provide a common set of best practice approaches that should be considered for all effective research grant proposal reviews. The federal government performs a critical role in American competitiveness and security by supporting basic research funded with taxpayer dollars. Effectively managing their allocation to scientists and researchers is a noble and crucial mission for advancing fundamental knowledge and deserves a heightened attention. Ensuring that proposals submitted are treated fairly and transparently is essential to both the health of any research program and also a duty to the public who ultimately funds the research. The paper describes the general requirements of a review process and at each step underlines the issues and suggests potential improvements and some fundamental requirements that should be included in any scientific review. We also included a series of tips geared to the scientific community. Our goals in this paper are 1) to demystify the process for everyone including policy makers who are sometimes flummoxed by the results of some scientific reviews, 2) to trigger some discussions about reviews and review process in the scientific community, 3) to inform scientists whose careers are directly impacted by review results about their own role in this process and 4) to suggest a road to more efficient, fairer and overall more transparent process. For experts in proposal reviews or for busy or impatient readers, the entire list of our recommendations is presented at the beginning. We describe in each section the context and rational of each recommendation.

categories: physics.soc-ph astro-ph.CO astro-ph.GA astro-ph.HE astro-ph.SR


title: Swift Observations of Gamma-Ray Burst Pulse Shapes: GRB Pulse Spectral Evolution Clarified

abstract: Isolated Swift gamma-ray burst (GRB) pulses, like their higher-energy BATSE counterparts, emit the bulk of their pulsed emission as a hard-to-soft component that can be fitted by the Norris et al. (2005) empirical pulse model. This signal is overlaid by a fainter, three-peaked signal that can be modeled by an empirical wave-like function (Hakkila and Preece, 2014): the two fits combine to reproduce GRB pulses with distinctive three-peaked shapes. The precursor peak appears on or before the pulse rise and is often the hardest component, the central peak is the brightest, and the decay peak converts exponentially decaying emission into a long, soft, power-law tail. Accounting for systematic instrumental differences, the general characteristics of the fitted pulses are remarkably similar. Isolated GRB pulses are dominated by hard-to-soft evolution; this is more pronounced for asymmetric pulses than for symmetric ones. Isolated GRB pulses can also exhibit intensity tracking behaviors that, when observed, are tied to the timing of the three peaks: pulses with the largest maximum hardnesses are hardest during the precursor, those with smaller maximum hardnesses are hardest during the central peak, and all pulses can re-harden during the central peak and/or during the decay peak. Since these behaviors are essentially seen in all isolated pulses, the distinction between "hard-to-soft" and "intensity-tracking" pulses really no longer applies. Additionally, the triple-peaked nature of isolated GRB pulses seems to indicate that energy is injected on three separate occasions during the pulse duration: theoretical pulse models need to account for this.


title: Grain Alignment: Role of Radiative Torques and Paramagnetic Relaxation

abstract: Polarization arising from aligned dust grains presents a unique opportunity to study magnetic fields in the diffuse interstellar medium and molecular clouds. Polarization from circumstellar regions, accretion disks and comet atmospheres can also be related to aligned dust.To reliably trace magnetic fields quantitative theory of grain alignment is required. Formulating the theory that would correspond to observations was one of the longstanding problems in astrophysics. Lately this problem has been successfully addressed, and in this review we summarize some of the most important theoretical advances in the theory of grain alignment by radiative torques (RATs) that act on realistic irregular dust grains. We discuss an analytical model of RATs and the ways to make RAT alignment more efficient, e.g. through paramagnetic relaxation when grains have inclusions with strong magnetic response. For very small grains for which RAT alignment is inefficient, we also discuss paramagnetic relaxation and a process termed resonance relaxation. We provide an extensive analysis of the observational tests of grain alignment theory.


title: Local Shtukas and Divisible Local Anderson Modules

abstract: We develop the analog of crystalline Dieudonn\'e theory for p-divisible groups in the arithmetic of function fields. In our theory p-divisible groups are replaced by divisible local Anderson modules, and Dieudonn\'e modules are replaced by local shtukas. We show that the categories of divisible local Anderson modules and of effective local shtukas are anti-equivalent over arbitrary base schemes. We also clarify their relation with formal Lie groups and with global objects like Drinfeld modules, Anderson's abelian t-modules and t-motives, and Drinfeld shtukas. Moreover, we discuss the existence of a Verschiebung map and apply it to deformations of local shtukas and divisible local Anderson modules. As a tool we use Faltings's and Abrashkin's theory of strict modules, which we review to some extent.


title: Cloud Offloading for Multi-Radio Enabled Mobile Devices

abstract: The advent of 5G networking technologies has increased the expectations from mobile devices, in that, more sophisticated, computationally intense applications are expected to be delivered on the mobile device which are themselves getting smaller and sleeker. This predicates a need for offloading computationally intense parts of the applications to a resource strong cloud. Parallely, in the wireless networking world, the trend has shifted to multi-radio (as opposed to multi-channel) enabled communications. In this paper, we provide a comprehensive computation offloading solution that uses the multiple radio links available for associated data transfer, optimally. Our contributions include: a comprehensive model for the energy consumption from the perspective of the mobile device; the formulation of the joint optimization problem to minimize the energy consumed as well as allocating the associated data transfer optimally through the available radio links and an iterative algorithm that converges to a locally optimal solution. Simulations on an HTC phone, running a 14-component application and using the Amazon EC2 as the cloud, show that the solution obtained through the iterative algorithm consumes only 3% more energy than the optimal solution (obtained via exhaustive search).


title: Truthful Multi-Parameter Auctions with Online Supply: an Impossible Combination

abstract: We study a basic auction design problem with online supply. There are two unit-demand bidders and two types of items. The first item type will arrive first for sure, and the second item type may or may not arrive. The auctioneer has to decide the allocation of an item immediately after each item arrives, but is allowed to compute payments after knowing how many items arrived. For this problem we show that there is no deterministic truthful and individually rational mechanism that, even with unbounded computational resources, gets any finite approximation factor to the optimal social welfare.


title: The Volume of the Trace-Nonnegative Polytope via the Irwin-Hall Distribution

abstract: In this work, we find an explicit expression for the volume of the trace nonnegative polytope via a generalization of the Irwin-Hall distribution. This volume is an upper bound for the volume of all projected, normalized realizable spectra. We provide ancillary results on realizable trace-zero spectra and pose several problems suitable for further inquiry.


title: The nature of the companion star in Circinus X-1

abstract: We present optical spectra and images of the X-ray binary Circinus X-1. The optical light curve of Cir X-1 is strongly variable, changing in brightness by 1.2 magnitudes in the space of four days. The shape of the light curve is consistent with that seen in the 1980s, when the X-ray and radio counterparts of the source were at least ten times as bright as they are currently. We detect strong, variable H-alpha emission lines, consisting of multiple components which vary with orbital phase. We estimate the extinction to the source from the strength of the diffuse interstellar bands and the Balmer decrement; the two methods give A_V = 7.6 +/- 0.6 mag and A_V > 9.1 mag respectively. The optical light curve can be modelled as arising from irradiation of the companion star by the central X-ray source, where a low-temperature star fills its Roche lobe in an orbit of moderate eccentricity (e ~ 0.4). We suggest that the companion star is over-luminous and under-dense, due to the impact of the supernova which occurred less than 5000 yr ago.


title: Light propagation in inhomogeneous and anisotropic cosmologies

abstract: The standard model of cosmology is based on the hypothesis that the Universe is spatially homogeneous and isotropic. When interpreting most observations, this cosmological principle is applied stricto sensu: the light emitted by distant sources is assumed to propagate through a Friedmann-Lema\^itre spacetime. The main goal of the present thesis was to evaluate how reliable this assumption is, especially when small scales are at stake. After having reviewed the laws of geometric optics in curved spacetime, and the standard interpretation of cosmological observables, the dissertation reports a comprehensive analysis of light propagation in Swiss-cheese models, designed to capture the clumpy character of the Universe. The resulting impact on the interpretation of the Hubble diagram is quantified, and shown to be relatively small, thanks to the cosmological constant. When applied to current supernova data, the associated corrections tend however to improve the agreement between the cosmological parameters inferred from the Hubble diagram and from the cosmic microwave background. This is a hint that the effect of small-scale structures on light propagation may become non-negligible in the era of precision cosmology. This motivated the development of a new theoretical framework, based on stochastic processes, which aims at describing small-scale gravitational lensing with a better accuracy. Regarding the isotropy side of the cosmological principle, this dissertation addresses, on the one hand, the potential effect of a large-scale anisotropy on light propagation, by solving all the equations of geometric optics in the Bianchi I spacetime. On the other hand, possible sources of such an anisotropy, namely scalar-vector models for inflation or dark energy, are analysed. Most of them turn out to be excluded as physically viable theories.


title: Embedded Ensemble Propagation for Improving Performance, Portability and Scalability of Uncertainty Quantification on Emerging Computational Architectures

abstract: Quantifying simulation uncertainties is a critical component of rigorous predictive simulation. A key component of this is forward propagation of uncertainties in simulation input data to output quantities of interest. Typical approaches involve repeated sampling of the simulation over the uncertain input data, and can require numerous samples when accurately propagating uncertainties from large numbers of sources. Often simulation processes from sample to sample are similar and much of the data generated from each sample evaluation could be reused. We explore a new method for implementing sampling methods that simultaneously propagates groups of samples together in an embedded fashion, which we call embedded ensemble propagation. We show how this approach takes advantage of properties of modern computer architectures to improve performance by enabling reuse between samples, reducing memory bandwidth requirements, improving memory access patterns, improving opportunities for fine-grained parallelization, and reducing communication costs. We describe a software technique for implementing embedded ensemble propagation based on the use of C++ templates and describe its integration with various scientific computing libraries within Trilinos. We demonstrate improved performance, portability and scalability for the approach applied to the simulation of partial differential equations on a variety of CPU, GPU, and accelerator architectures, including up to 131,072 cores on a Cray XK7 (Titan).