title: Quantum quenches in the many-body localized phase

abstract: Many-body localized (MBL) systems are characterized by the absence of transport and thermalization, and therefore cannot be described by conventional statistical mechanics. In this paper, using analytic arguments and numerical simulations, we study the behaviour of local observables in an isolated MBL system following a quantum quench. For the case of a global quench, we find that the local observables reach stationary, highly non-thermal values at long times as a result of slow dephasing characteristic of the MBL phase. These stationary values retain the local memory of the initial state due to the existence of local integrals of motion in the MBL phase. The temporal fluctuations around stationary values exhibit universal power-law decay in time, with an exponent set by the localization length and the diagonal entropy of the initial state. Such a power-law decay holds for any local observable and is related to the logarithmic in time growth of entanglement in the MBL phase. This behaviour distinguishes the MBL phase from both the Anderson insulator (where no stationary state is reached), and from the ergodic phase (where relaxation is expected to be exponential). For the case of a local quench, we also find a power-law approach of local observables to their stationary values when the system is prepared in a mixed state. Quench protocols considered in this paper can be naturally implemented in systems of ultra cold atoms in disordered optical lattices, and the behaviour of local observables provides a direct experimental signature of many-body localization.


title: Intersection theory and the Alesker product

abstract: Alesker has introduced the space $\mathcal V^\infty(M)$ of {\it smooth valuations} on a smooth manifold $M$, and shown that it admits a natural commutative multiplication. Although Alesker's original construction is highly technical, from a moral perspective this product is simply an artifact of the operation of intersection of two sets. Subsequently Alesker and Bernig gave an expression for the product in terms of differential forms. We show how the Alesker-Bernig formula arises naturally from the intersection interpretation, and apply this insight to give a new formula for the product of a general valuation with a valuation that is expressed in terms of intersections with a sufficiently rich family of smooth polyhedra.


title: Automorphism groups of countable algebraically closed graphs and endomorphisms of the random graph

abstract: We establish links between countable algebraically closed graphs and the endomorphisms of the countable universal graph $R$. As a consequence we show that, for any countable graph $\Gamma$, there are uncountably many maximal subgroups of the endomorphism monoid of $R$ isomorphic to the automorphism group of $\Gamma$. Further structural information about End $R$ is established including that Aut $\Gamma$ arises in uncountably many ways as a Sch\"{u}tzenberger group. Similar results are proved for the countable universal directed graph and the countable universal bipartite graph.


title: The impact of macroscopic epistasis on long-term evolutionary dynamics

abstract: Genetic interactions can strongly influence the fitness effects of individual mutations, yet the impact of these epistatic interactions on evolutionary dynamics remains poorly understood. Here we investigate the evolutionary role of epistasis over 50,000 generations in a well-studied laboratory evolution experiment in E. coli. The extensive duration of this experiment provides a unique window into the effects of epistasis during long-term adaptation to a constant environment. Guided by analytical results in the weak-mutation limit, we develop a computational framework to assess the compatibility of a given epistatic model with the observed patterns of fitness gain and mutation accumulation through time. We find that a decelerating fitness trajectory alone provides little power to distinguish between competing models, including those that lack any direct epistatic interactions between mutations. However, when combined with the mutation trajectory, these observables place strong constraints on the set of possible models of epistasis, ruling out many existing explanations of the data. Instead, we find that the data are consistent with "two-epoch" model of adaptation, in which an initial burst of diminishing returns epistasis is followed by a steady accumulation of mutations under a constant distribution of fitness effects. Our results highlight the need for additional DNA sequencing of these populations, as well as for more sophisticated models of epistasis that are compatible with all of the experimental data.


title: Limits of geometries

abstract: A geometric transition is a continuous path of geometric structures that changes type, meaning that the model geometry, i.e. the homogeneous space on which the structures are modeled, abruptly changes. In order to rigorously study transitions, one must define a notion of geometric limit at the level of homogeneous spaces, describing the basic process by which one homogeneous geometry may transform into another. We develop a general framework to describe transitions in the context that both geometries involved are represented as sub-geometries of a larger ambient geometry. Specializing to the setting of real projective geometry, we classify the geometric limits of any sub-geometry whose structure group is a symmetric subgroup of the projective general linear group. As an application, we classify all limits of three-dimensional hyperbolic geometry inside of projective geometry, finding Euclidean, Nil, and Sol geometry among the limits. We prove, however, that the other Thurston geometries, in particular $\mathbb{H}^2 \times \mathbb{R}$ and $\widetilde{\operatorname{SL}_2 \mathbb{R}}$, do not embed in any limit of hyperbolic geometry in this sense.


title: Augmentation Rank of Satellites with Braid Pattern

abstract: Given a knot $K$ in $S^3$, a question raised by Cappell and Shaneson asks if the meridional rank of $K$ equals the bridge number of $K$. Using augmentations in knot contact homology we consider the persistence of equality between these two invariants under satellite operations on $K$ with a braid pattern. In particular, we answer the question in the affirmative for a large class of iterated torus knots.


title: A New Approach to Customization of Collision Warning Systems to Individual Drivers

abstract: This paper discusses the need for individualizing safety systems and proposes an approach including the Real-Time estimation of the distribution of brake response times for an individual driver. While maintaining high level of safety, the collision warning system should send "tailored" responses to the driver. This method could be the first step to show that safety applications would potentially benefit from customizing to individual drivers' characteristics using VANET. Our simulation results show that, as one of the imminent and preliminary outcomes of the new improved system, the number of false alarms will be reduced by more than 40%. We think this tactic can reach to even beyond the safety applications for designing the future innovative systems.


title: Ameliorate Threshold Distributed Energy Efficient Clustering Algorithm for Heterogeneous Wireless Sensor Networks

abstract: Ameliorating the lifetime in heterogeneous wireless sensor network is an important task because the sensor nodes are limited in the resource energy. The best way to improve a WSN lifetime is the clustering based algorithms in which each cluster is managed by a leader called Cluster Head. Each other node must communicate with this CH to send the data sensing. The nearest base station nodes must also send their data to their leaders, this causes a loss of energy. In this paper, we propose a new approach to ameliorate a threshold distributed energy efficient clustering protocol for heterogeneous wireless sensor networks by excluding closest nodes to the base station in the clustering process. We show by simulation in MATLAB that the proposed approach increases obviously the number of the received packet messages and prolongs the lifetime of the network compared to TDEEC protocol.


title: Fast Shortest Path Routing in Transportation Networks with Time-Dependent Road Speeds

abstract: The current paper deals with the subject of shortest path routing in transportation networks (in terms of travelling time), where the speed in several of the network's roads is a function of the time interval. The main contribution of the paper is a procedure that is faster compared to the conventional approaches, that derives the road's traversal time according to the time instant of departure, for the case where the road's speed has a constant value inside each time interval (in general, different value for each time interval). Furthermore, the case where the road's speed is a linear function of time inside each time interval (in general, different linear function for each time interval) is investigated. A procedure that derives the road's traversal time according to the time instant of departure is proposed for this case as well. The proposed procedures are combined with Dijkstra's algorithm and the resulting algorithms, that are practically applicable and of low complexity, provide optimal shortest path routing in the networks under investigation.


title: An Energetic Variational Approach for ion transport

abstract: The transport and distribution of charged particles are crucial in the study of many physical and biological problems. In this paper, we employ an Energy Variational Approach to derive the coupled Poisson-Nernst-Planck-Navier-Stokes system. All physics is included in the choices of corresponding energy law and kinematic transport of particles. The variational derivations give the coupled force balance equations in a unique and deterministic fashion. We also discuss the situations with different types of boundary conditions. Finally, we show that the Onsager's relation holds for the electrokinetics, near the initial time of a step function applied field.


title: Non-Detection of X-Ray Emission From Sterile Neutrinos in Stacked Galaxy Spectra

abstract: We conduct a comprehensive search for X-ray emission lines from sterile neutrino dark matter, motivated by recent claims of unidentified emission lines in the stacked X-ray spectra of galaxy clusters and the centers of the Milky Way and M31. Since the claimed emission lines lie around 3.5 keV, we focus on galaxies and galaxy groups (masking the central regions), since these objects emit very little radiation above $\sim 2$ keV and offer a clean background against which to detect emission lines. We develop a formalism for maximizing the signal-to-noise of decaying dark matter emission lines by weighing each X-ray event according to the expected dark matter profile. In total, we examine 81 and 89 galaxies with Chandra and XMM-Newton respectively, totaling 15.0 and 14.6 Ms of integration time. We find no significant evidence of any emission lines, placing strong constraints on the mixing angle of sterile neutrinos with masses between 4.8-12.4 keV. In particular, if the 3.57 keV feature from Bulbul et al. (2014) were due to 7.1 keV sterile neutrino emission, we would have detected it at $4.4\sigma$ and $11.8\sigma$ in our two samples. The most conservative estimates of the systematic uncertainties reduce these constraints to $4.4\sigma$ and 7.8$\sigma$, or letting the line energy vary between 3.50 and 3.60 keV reduces these constraints to $2.7\sigma$ and $11.0\sigma$ respectively. Unlike previous constraints, our measurements do not depend on the model of the X-ray background or on the assumed logarithmic slope of the center of the dark matter profile.

categories: astro-ph.HE astro-ph.CO astro-ph.GA hep-ph


title: Do Two Temperature Debris Disks Have Multiple Belts?

abstract: We present a study of debris disks whose spectra are well modelled by dust emission at two different temperatures. These disks are typically assumed to be a sign of multiple belts, which in only a few cases have been confirmed via high resolution observations. We first compile a sample of two-temperature disks to derive their properties, summarised by the ratios of the warm and cool component temperatures and fractional luminosities. The ratio of warm to cool temperatures is constant in the range 2-4, and the temperatures of both warm and cool components increases with stellar mass. We then explore whether this emission can arise from dust in a single narrow belt, with the range of temperatures arising from the size variation of grain temperatures. This model can produce two-temperature spectra for Sun-like stars, but is not supported where it can be tested by observed disk sizes and far-IR/mm spectral slopes. Therefore, while some two-temperature disks arise from single belts, it is probable that most have multiple spatial components. These disks are plausibly similar to the outer Solar System's configuration of Asteroid and Edgeworth-Kuiper belts separated by giant planets. Alternatively, the inner component could arise from inward scattering of material from the outer belt, again due to intervening planets. In either case, we suggest that the ratio of warm/cool component temperatures is indicative of the scale of outer planetary systems, which typically span a factor of about ten in radius.


title: A low-loss, broadband antenna for efficient photon collection from a coherent spin in diamond

abstract: We report the creation of a low-loss, broadband optical antenna giving highly directed output from a coherent single spin in the solid-state. The device, the first solid-state realization of a dielectric antenna, is engineered for individual nitrogen vacancy (NV) electronic spins in diamond. We demonstrate a directionality close to 10. The photonic structure preserves the high spin coherence of single crystal diamond (T2>100us). The single photon count rate approaches a MHz facilitating efficient spin readout. We thus demonstrate a key enabling technology for quantum applications such as high-sensitivity magnetometry and long-distance spin entanglement.


title: Dynamics of Gauge Field Inflation

abstract: We analyze the existence and stability of dynamical attractor solutions for cosmological inflation driven by the coupling between fermions and a gauge field. Assuming a spatially homogeneous and isotropic gauge field and fermion current, the interacting fermion equation of motion reduces to that of a free fermion up to a phase shift. Consistency of the model is ensured via the Stuckelberg mechanism. We prove the existence of exactly one stable solution, and demonstrate the stability numerically. Inflation arises without fine tuning, and does not require postulating any effective potential or non-standard coupling.


title: Charged Higgs Search via $AW^\pm/HW^\pm$ Channel

abstract: Models of electroweak symmetry breaking with extended Higgs sectors are theoretically well motivated. In this study, we focus on models with a low energy spectrum containing a pair of charged scalars $H^\pm$, as well as a light scalar H and/or a pseudoscalar A. We study the $H^\pm tb$ associated production with $H^\pm \to AW/HW$, which could reach sizable branching fractions in certain parameter regions. With detailed collider analysis, we obtain the exclusion bounds as well as discovery reach at the 14 TeV LHC for the process $pp \to H^\pm tb \to AW^\pm tb/HW^\pm tb \to \tau\tau bbWW, bbbbWW$. We find that for a daughter particle mass of 70 GeV, the 95% C.L. exclusion reach in $\sigma$xBR varies from about 60 fb to 25 fb, for $m_{H^\pm}$ ranging from 150 GeV to 500 GeV with 300 fb$^{-1}$ integrated luminosity in the $\tau\tau$ mode. We further interpret these bounds in the context of Type II Two Higgs Doublet Model. The exclusion region in the $m_{H^\pm}-\tan\beta$ plane can be extended to $m_{H^\pm}=$ 600 GeV, while discovery is possible for $m_{H^\pm}\lesssim$ 400 GeV with 300 fb$^{-1}$ integrated luminosity. The exotic decay mode $H^\pm \to AW^\pm/HW^\pm$ offers a complementary channel to the conventional mode $H^\pm \to \tau\nu$ for charged Higgs searches.


title: Pseudo Three-dimensional Maps of the Diffuse Interstellar Band at 862 nm

abstract: The diffuse interstellar bands (DIBs) are absorption lines observed in visual and near infrared spectra of stars. Understanding their origin in the interstellar medium (ISM) is one of the oldest problems in astronomical spectroscopy, as DIBs ahave been known since 1922. In a completely new approach to understanding DIBs, we combined information from nearly 500,000 stellar spectra obtained by the massive spectroscopic survey RAVE (Radial Velocity Experiment) to produce the first pseudo three-dimensional map of the strength of the DIB at 8620~{\AA}ngstroms covering the nearest 3 kiloparsecs from the Sun, and show that it follows our independently constructed spatial distribution of extinction by interstellar dust along the Galactic plane. Despite having a similar distribution in the Galactic plane, the DIB 8620 carrier has a significantly larger vertical scale height than the dust. Even if one DIB may not represent the general DIB population, our observations outline the future direction of DIB research.


title: GRB 140206A: the most distant polarized Gamma-Ray Burst

abstract: The nature of the prompt gamma-ray emission of Gamma-Ray Bursts (GRBs) is still far from being completely elucidated. The measure of linear polarization is a powerful tool that can be used to put further constraints on the content and magnetization of the GRB relativistic outflows, as well as on the radiation processes at work. To date only a handful of polarization measurements are available for the prompt emission of GRBs. Here we present the analysis of the prompt emission of GRB 140206A, obtained with INTEGRAL/IBIS, Swift/BAT, and Fermi/GBM. Using INTEGRAL/IBIS as a Compton polarimeter we were able to constrain the linear polarization level of the second peak of this GRB as being larger than 28% at 90% c.l. We also present the GRB afterglow optical spectroscopy obtained at the Telescopio Nazionale Galileo (TNG), which allowed us the measure the distance of this GRB, z=2.739. This distance value together with the polarization measure obtained with IBIS, allowed us to derive the deepest and most reliable limit to date (xi <1x10-16) on the possibility of Lorentz Invariance Violation, measured through the vacuum birefringence effect on a cosmological source.


title: Interpreting high [O III]/Hbeta ratios with maturing starbursts

abstract: Star forming galaxies at high redshift show ubiquitously high ionization parameters, as measured by the ratio of optical emission lines. We demonstrate that local (z < 0.2) sources selected as Lyman break analogues also manifest high line ratios with a typical [O III]/H beta = 3.36(+0.14,-0.04) - comparable to all but the highest ratios seen in star forming galaxies at z ~ 2-4. We argue that the stellar population synthesis code BPASS can explain the high ionization parameters required through the ageing of rapidly formed star populations, without invoking any AGN contribution. Binary stellar evolution pathways prolong the age interval over which a starburst is likely to show elevated line ratios, relative to those predicted by single stellar evolution codes. As a result, model galaxies at near-Solar metallicities and with ages of up to ~100 Myr after a starburst typically have a line ratio [O III]/H beta~3, consistent with those seen in Lyman break galaxies and local sources with similar star formation densities. This emphasises the importance of including binary evolution pathways when simulating the nebular line emission of young or bursty stellar populations.


title: Charge noise, spin-orbit coupling, and coherence of single-spin qubits

abstract: Spin-orbit coupling is ubiquitous in quantum dot quantum computing architectures, and makes spin qubits susceptible to charge noise. We derive a Hamiltonian describing the effect of spin-orbit and noise on a single-spin qubit in a quantum dot. Relaxation is due to noise coupling different orbital levels and is dominated by screened whole charge defects near the dot. Dephasing stems from noise causing relative fluctuations between orbital levels, and is driven by screened whole charge defects and unscreened dipole defects in the substrate. Dephasing times are vastly different between common materials such as Si and GaAs. They can be enhanced by increasing gate fields, choosing materials with weak spin-orbit such as Si, making dots narrower, or using accumulation dots.


title: Transit Light Curves with Finite Integration Time: Fisher Information Analysis

abstract: Kepler has revolutionized the study of transiting planets with its unprecedented photometric precision on more than 150,000 target stars. Most of the transiting planet candidates detected by Kepler have been observed as long-cadence targets with 30 minute integration times, and the upcoming Transiting Exoplanet Survey Satellite (TESS) will record full frame images with a similar integration time. Integrations of 30 minutes affect the transit shape, particularly for small planets and in cases of low signal-to-noise. Using the Fisher information matrix technique, we derive analytic approximations for the variances and covariances on the transit parameters obtained from fitting light curve photometry collected with a finite integration time. We find that binning the light curve can significantly increase the uncertainties and covariances on the inferred parameters when comparing scenarios with constant total signal-to-noise (constant total integration time in the absence of read noise). Uncertainties on the transit ingress/egress time increase by a factor of 34 for Earth-size planets and 3.4 for Jupiter-size planets around Sun-like stars for integration times of 30 minutes compared to instantaneously-sampled light curves. Similarly, uncertainties on the mid-transit time for Earth and Jupiter-size planets increase by factors of 3.9 and 1.4. Uncertainties on the transit depth are largely unaffected by finite integration times. While correlations among the transit depth, ingress duration, and transit duration all increase in magnitude with longer integration times, the mid-transit time remains uncorrelated with the other parameters. We provide code in Python and Mathematica for predicting the variances and covariances at www.its.caltech.edu/~eprice.


title: Statistical diagnostics to identify Galactic foregrounds in B-mode maps

abstract: Recent developments in the search for inflationary gravitational waves in the cosmic microwave background (CMB) polarization motivate the search for new diagnostics to distinguish the Galactic foreground contribution to B modes from the cosmic signal. We show that B modes from these foregrounds should exhibit a local hexadecapolar departure in power from statistical isotropy (SI). We present a simple algorithm to search for a uniform SI violation of this sort, as may arise in a sufficiently small patch of sky. We then show how to search for these effects if the orientation of the SI violation varies across the survey region, as is more likely to occur in surveys with more sky coverage. If detected, these departures from Gaussianity would indicate some level of Galactic foreground contamination in the B-mode maps. Given uncertainties about foreground properties, though, caution should be exercised in attributing a null detection to an absence of foregrounds.


title: Constraining topological defects with temperature and polarization anisotropies

abstract: We analyse the possible contribution of topological defects to cosmic microwave anisotropies, both temperature and polarisation. We allow for the presence of both inflationary scalars and tensors, and of polarised dust foregrounds that may contribute to or dominate the B-mode polarisation signal. We confirm and quantify our previous statements that topological defects on their own are a poor fit to the B-mode signal. However, adding topological defects to a models with a tensor component or a dust component improves the fit around $\ell=200$. Fitting simultaneously to both temperature and polarisation data, we find that textures fit almost as well as tensors ($\Delta\chi^2 = 2.0$), while Abelian Higgs strings are ruled out as the sole source of the B-mode signal at low $\ell$. The 95% confidence upper limits on models combining defects and dust are $G\mu < 2.7\times 10^{-7}$ (Abelian Higgs strings), $G\mu < 9.8\times 10^{-7}$ (semilocal strings) and $G\mu < 7.3\times 10^{-7}$ (textures), a small reduction on the $Planck$ bounds. The most economical fit overall is obtained by the standard $\Lambda$CDM model with a polarised dust component.


title: Towards phase transitions between discrete and continuum quantum spacetime from the Renormalization Group

abstract: We establish the functional Renormalization Group as an exploratory tool to investigate a possible phase transition between a pre-geometric discrete phase and a geometric continuum phase in quantum gravity. In this paper, based on the analysis of [1], we study three new aspects of the double-scaling limit of matrix models as Renormalization Group fixed points: Firstly, we investigate multicritical fixed points, which are associated with quantum gravity coupled to conformal matter. Secondly, we discuss an approximation that reduces the scheme dependence of our results as well as computational effort while giving good numerical results. This is a consequence of the approximation being a solution to the unitary Ward-identity associated to the U(N) symmetry of the hermitian matrix model. Thirdly, we discuss a scenario that relates the double scaling limit to fixed points of continuum quantum gravity.


title: And yet it moves: The dangers of artificially fixing the Milky Way center of mass in the presence of a massive Large Magellanic Cloud

abstract: Motivated by recent studies suggesting that the Large Magellanic Cloud (LMC) could be significantly more massive than previously thought, we explore whether the approximation of an inertial Galactocentric reference frame is still valid in the presence of such a massive LMC. We find that previous estimates of the LMC's orbital period and apocentric distance derived assuming a fixed Milky Way are significantly shortened for models where the Milky Way is allowed to move freely in response to the gravitational pull of the LMC. Holding other parameters fixed, the fraction of models favoring first infall is reduced. Due to this interaction, the Milky Way center of mass within the inner 50 kpc can be significantly displaced in phase-space in a very short period of time that ranges from 0.3 to 0.5 Gyr by as much as 30 kpc and 75 km/s. Furthermore, we show that the gravitational pull of the LMC and response of the Milky Way are likely to significantly affect the orbit and phase space distribution of tidal debris from the Sagittarius dwarf galaxy (Sgr). Such effects are larger than previous estimates based on the torque of the LMC alone. As a result, Sgr deposits debris in regions of the sky that are not aligned with the present-day Sgr orbital plane. In addition, we find that properly accounting for the movement of the Milky Way around its common center of mass with the LMC significantly modifies the angular distance between apocenters and tilts its orbital pole, alleviating tensions between previous models and observations. While these models are preliminary in nature, they highlight the central importance of accounting for the mutual gravitational interaction between the MW and LMC when modeling the kinematics of objects in the Milky Way and Local Group.


title: Metastable spin textures and Nambu-Goldstone modes of a ferromagnetic spin-1 Bose-Einstein condensate confined in a ring trap

abstract: We investigate the metastability of a ferromagnetic spin-1 Bose-Einstein condensate confined in a quasi-one-dimensional rotating ring trap by solving the spin-1 Gross-Pitaevskii equation. We find analytical solutions that exhibit spin textures. By performing linear stability analysis, it is shown that the solutions can become metastable states. We also find that the number of Nambu-Goldstone modes changes at a certain rotation velocity without changing the continuous symmetry of the order parameter.


title: The Solar Twin Planet Search. I. Fundamental parameters of the stellar sample

abstract: We are carrying out a search for planets around a sample of solar twin stars using the HARPS spectrograph. The goal of this project is to exploit the advantage offered by solar twins to obtain chemical abundances of unmatched precision. This survey will enable new studies of the stellar composition -- planet connection. Here we used the MIKE spectrograph on the Magellan Clay Telescope to acquire high resolution, high signal-to-noise ratio spectra of our sample stars. We measured the equivalent widths of iron lines and used strict differential excitation/ionization balance analysis to determine atmospheric parameters of unprecedented internal precision (DTeff=7K, Dlogg=0.019, D[Fe/H]=0.006dex, Dvt=0.016km/s). Reliable relative ages and highly precise masses were then estimated using theoretical isochrones. The spectroscopic parameters we derived are in good agreement with those measured using other independent techniques. The root-mean-square scatter of the differences seen is fully compatible with the observational errors, demonstrating, as assumed thus far, that systematic uncertainties in the stellar parameters are negligible in the study of solar twins. We find a tight activity-age relation for our sample stars, which validates the internal precision of our dating method. Furthermore, we find that the solar cycle is perfectly consistent both with this trend and its star-to-star scatter. We present the largest sample of solar twins analyzed homogeneously using high quality spectra. The fundamental parameters derived from this work will be employed in subsequent work that aims to explore the connections between planet formation and stellar chemical composition.


title: A realistic assessment of the CTA sensitivity to dark matter annihilation

abstract: We estimate the sensitivity of the upcoming CTA gamma-ray telescope to DM annihilation at the Galactic centre, improving on previous analyses in a number of significant ways. First, we perform a detailed analyses of all backgrounds, including diffuse astrophysical emission for the first time in a study of this type. Second, we present a statistical framework for including systematic errors and estimate the consequent degradation in sensitivity. These errors may come from e.g. event reconstruction, Monte Carlo determination of the effective area or uncertainty in atmospheric conditions. Third, we show that performing the analysis on a set of suitably optimised regions of interest makes it possible to partially compensate for the degradation in sensitivity caused by systematics and diffuse emission. To probe dark matter with the canonical thermal annihilation cross-section, CTA systematics like non-uniform variations in acceptance over a single field of view must be kept below the 0.3% level, unless the dark matter density rises more steeply in the centre of the Galaxy than predicted by a typical Navarro-Frenk-White or Einasto profile. For a contracted $r^{-1.3}$ profile, and systematics at the 1% level, CTA can probe annihilation to $b\bar{b}$ at the canonical thermal level for dark matter masses between 100 GeV and 10 TeV.


title: Givental J-functions, Quantum integrable systems, AGT relation with surface operator

abstract: We study 4d $\mathcal{N}=2$ gauge theories with a co-dimension two full surface operator, which exhibit a fascinating interplay of supersymmetric gauge theories, equivariant Gromov-Witten theory and geometric representation theory. For pure Yang-Mills and $\mathcal{N}=2^*$ theory, we describe a full surface operator as the 4d gauge theory coupled to a 2d $\mathcal{N}=(2,2)$ gauge theory. By supersymmetric localizations, we present the exact partition functions of both 4d and 2d theories which satisfy integrable equations. In addition, the form of the structure constants with a semi-degenerate field in SL(N,R) WZNW model is predicted from one-loop determinants of 4d gauge theories with a full surface operator via the AGT relation.


title: Efficient geodesics and an effective algorithm for distance in the complex of curves

abstract: We give an algorithm for determining the distance between two vertices of the complex of curves. While there already exist such algorithms, for example by Leasure, Shackleton, and Webb, our approach is new, simple, and more effective for all distances accessible by computer. Our method gives a new preferred finite set of geodesics between any two vertices of the complex, called efficient geodesics, which are different from the tight geodesics introduced by Masur and Minsky.


title: MICC: A tool for computing short distances in the curve complex

abstract: The complex of curves $\mathcal{C}(S_g)$ of a closed orientable surface of genus $g \geq 2$ is the simplicial complex having its vertices, $\mathcal{C}^0(S_g)$, are isotopy classes of essential curves in $S_g$. Two vertices co-bound an edge of the $1$-skeleton, $\mathcal{C}^1(S_g)$, if there are disjoint representatives in $S_g$. A metric is obtained on $\mathcal{C}^0(S_g)$ by assigning unit length to each edge of $\mathcal{C}^1(S_g)$. Thus, the distance between two vertices, $d(v,w)$, corresponds to the length of a geodesic---a shortest edge-path between $v$ and $w$ in $\mathcal{C}^1 (S_g)$. Recently, Birman, Margalit and the second author introduced the concept of {\em initially efficient geodesics} in $\mathcal{C}^1(S_g)$ and used them to give a new algorithm for computing the distance between vertices. In this note we introduce the software package MICC ({\em Metric in the Curve Complex}), a partial implementation of the initially efficient geodesic algorithm. We discuss the mathematics underlying MICC and give applications. In particular, we give examples of distance four vertex pairs, for $g=2$ and 3. Previously, there was only one known example, in genus $2$, due to John Hempel.


title: How many nucleosynthesis processes exist at low metallicity?

abstract: Abundances of low-metallicity stars offer a unique opportunity to understand the contribution and conditions of the different processes that synthesize heavy elements. Many old, metal-poor stars show a robust abundance pattern for elements heavier than Ba, and a less robust pattern between Sr and Ag. Here we probe if two nucleosynthesis processes are sufficient to explain the stellar abundances at low metallicity, and we carry out a site independent approach to separate the contribution from these two processes or components to the total observationally derived abundances. Our approach provides a method to determine the contribution of each process to the production of elements such as Sr, Zr, Ba, and Eu. We explore the observed star-to-star abundance scatter as a function of metallicity that each process leads to. Moreover, we use the deduced abundance pattern of one of the nucleosynthesis components to constrain the astrophysical conditions of neutrino-driven winds from core-collapse supernovae.


title: Initial data for binary neutron stars with adjustable eccentricity

abstract: Binary neutron stars in circular orbits can be modeled as helically symmetric, i.e., stationary in a rotating frame. This symmetry gives rise to a first integral of the Euler equation, often employed for constructing equilibrium solutions via iteration. For eccentric orbits, however, the lack of helical symmetry has prevented the use of this method, and the numerical relativity community has often resorted to constructing initial data by superimposing boosted spherical stars without solving the Euler equation. The spuriously excited neutron star oscillations seen in evolutions of such data arise because such configurations lack the appropriate tidal deformations and are stationary in a linearly comoving---rather than rotating---frame. We consider eccentric configurations at apoapsis that are instantaneously stationary in a rotating frame. We extend the notion of helical symmetry to eccentric orbits, by approximating the elliptical orbit of each companion as instantaneously circular, using the ellipse's inscribed circle. The two inscribed helical symmetry vectors give rise to approximate instantaneous first integrals of the Euler equation throughout each companion. We use these integrals as the basis of a self-consistent iteration of the Einstein constraints to construct conformal thin-sandwich initial data for eccentric binaries. We find that the spurious stellar oscillations are reduced by at least an order of magnitude, compared with those found in evolutions of superposed initial data. The tidally induced oscillations, however, are physical and qualitatively similar to earlier evolutions. Finally, we show how to incorporate radial velocity due to radiation reaction in our inscribed helical symmetry vectors, which would allow one to obtain truly non-eccentric initial data when our eccentricity parameter $e$ is set to zero.


title: The WFCAM Multi-wavelength Variable Star Catalog

abstract: Stellar variability in the near-infrared (NIR) remains largely unexplored. The exploitation of public science archives with data-mining methods offers a perspective for the time-domain exploration of the NIR sky. We perform a comprehensive search for stellar variability using the optical-NIR multi-band photometric data in the public Calibration Database of the WFCAM Science Archive (WSA), with the aim of contributing to the general census of variable stars, and to extend the current scarce inventory of accurate NIR light curves for a number of variable star classes. We introduce new variability indices designed for multi-band data with correlated sampling, and apply them for pre-selecting variable star candidates, i.e., light curves that are dominated by correlated variations, from noise-dominated ones. Pre-selection criteria are established by robust numerical tests for evaluating the response of variability indices to colored noise characteristic to the data. We find 275 periodic variable stars and an additional 44 objects with suspected variability with uncertain periods or apparently aperiodic variation. Only 44 of these objects had been previously known, including 11 RR~Lyrae stars in the outskirts of the globular cluster M3 (NGC~5272). We provide a preliminary classification of the new variable stars that have well-measured light curves, but the variability types of a large number of objects remain ambiguous. We classify most of the new variables as contact binary stars, but we also find several pulsating stars, among which 34 are probably new field RR~Lyrae and 3 are likely Cepheids. We also identify 32 highly reddened variable objects close to previously known dark nebulae, suggesting that these are embedded young stellar objects. We publish our results and all light-curve data as the WFCAM Variable Star Catalog.


title: Recent developments of analysis for hydrodynamic flow of nematic liquid crystals

abstract: The study of hydrodynamics of liquid crystal leads to many fasci- nating mathematical problems, which has prompted various interesting works recently. This article reviews the static Oseen-Frank theory and surveys some recent progress on the existence, regularity, uniqueness, and large time asymp- totic of the hydrodynamic flow of nematic liquid crystals. We will also propose a few interesting questions for future investigations.


title: Shear Transformation Zones: State Determined or Protocol Dependent?

abstract: The concept of a Shear Transformation Zone (STZ) refers to a region in an amorphous solid that undergoes a plastic event when the material is put under an external mechanical load. An important question that had accompanied the development of the theory of plasticity in amorphous solids for many years now is whether an STZ is a {\em region} existing in the material (which can be predicted by analyzing the unloaded material), or is it an {\em event} that depends on the loading protocol (i.e., the event cannot be predicted without following the protocol itself). In this Letter we present strong evidence that the latter is the case. Infinitesimal changes of protocol result in macroscopically big jumps in the positions of plastic events, meaning that these can never be predicted from considering the unloaded material.


title: BET: Bayesian Ensemble Trees for Clustering and Prediction in Heterogeneous Data

abstract: We propose a novel "tree-averaging" model that utilizes the ensemble of classification and regression trees (CART). Each constituent tree is estimated with a subset of similar data. We treat this grouping of subsets as Bayesian ensemble trees (BET) and model them as an infinite mixture Dirichlet process. We show that BET adapts to data heterogeneity and accurately estimates each component. Compared with the bootstrap-aggregating approach, BET shows improved prediction performance with fewer trees. We develop an efficient estimating procedure with improved sampling strategies in both CART and mixture models. We demonstrate these advantages of BET with simulations, classification of breast cancer and regression of lung function measurement of cystic fibrosis patients. Keywords: Bayesian CART; Dirichlet Process; Ensemble Approach; Heterogeneity; Mixture of Trees.


title: Schwinger Effect in 4D de Sitter Space and Constraints on Magnetogenesis in the Early Universe

abstract: We investigate pair creation by an electric field in four-dimensional de Sitter space. The expectation value of the induced current is computed, using the method of adiabatic regularization. Under strong electric fields the behavior of the current is similar to that in flat space, while under weak electric fields the current becomes inversely proportional to the mass squared of the charged field. Thus we find that the de Sitter space obtains a large conductivity under weak electric fields in the presence of a charged field with a tiny mass. We then apply the results to constrain electromagnetic fields in the early universe. In particular, we study cosmological scenarios for generating large-scale magnetic fields during the inflationary era. Electric fields generated along with the magnetic fields can induce sufficiently large conductivity to terminate the phase of magnetogenesis. For inflationary magnetogenesis models with a modified Maxwell kinetic term, the generated magnetic fields cannot exceed 10^{-30} G on Mpc scales in the present epoch, when a charged field carrying an elementary charge with mass of order the Hubble scale or smaller exists in the Lagrangian. Similar constraints from the Schwinger effect apply for other magnetogenesis mechanisms.


title: Possible binary progenitors for the type Ib supernova iPTF13bvn

abstract: Cao et al. (2013) reported a possible progenitor detection for the type Ib supernovae iPTF13bvn for the first time. We find that the progenitor is in fact brighter than the magnitudes previously reported by approximately 0.7 to 0.2 mag with a larger error in the bluer filters. We compare our new magnitudes to our large set of binary evolution models and find that many binary models with initial masses in the range of 10 to 20M(Sun) match this new photometry and other constraints suggested from analysing the supernova. In addition these lower mass stars retain more helium at the end of the model evolution indicating that they are likely to be observed as type Ib supernovae rather than their more massive, Wolf-Rayet counter parts. We are able to rule out typical Wolf-Rayet models as the progenitor because their ejecta masses are too high and they do not fit the observed SED unless they have a massive companion which is the observed source at the supernova location. Therefore only late-time observations of the location will truly confirm if the progenitor was a helium giant and not a Wolf-Rayet star.


title: Self Organization Map based Texture Feature Extraction for Efficient Medical Image Categorization

abstract: Texture is one of the most important properties of visual surface that helps in discriminating one object from another or an object from background. The self-organizing map (SOM) is an excellent tool in exploratory phase of data mining. It projects its input space on prototypes of a low-dimensional regular grid that can be effectively utilized to visualize and explore properties of the data. This paper proposes an enhancement extraction method for accurate extracting features for efficient image representation it based on SOM neural network. In this approach, we apply three different partitioning approaches as a region of interested (ROI) selection methods for extracting different accurate textural features from medical image as a primary step of our extraction method. Fisherfaces feature selection is used, for selecting discriminated features form extracted textural features. Experimental result showed the high accuracy of medical image categorization with our proposed extraction method. Experiments held on Mammographic Image Analysis Society (MIAS) dataset.


title: On electromagnetic momentum of an electric dipole in a magnetic field

abstract: The total linear electromagnetic field momentum $\mathbf P_{\mathrm{em}}$ of a stationary electric dipole $\mathbf p$ in a static magnetic field $\mathbf B$ is considered. The expression $\mathbf P_{\mathrm{em}} = \frac12\mathbf B \times \mathbf p$, which has previously been implied to hold in all static magnetic field situations, is not valid in general. The contribution of the electromagnetic momentum of the fringing fields of the dipole is discussed. It is shown that when either the static magnetic field or the electric dipole moment is changed, the mechanical impulse on the system equals $-\Delta\mathbf P_{\mathrm{em}}$, and hidden momentum does not need to be invoked in order to conserve total momentum.


title: The Masses and Spins of Neutron Stars and Stellar-Mass Black Holes

abstract: Stellar-mass black holes and neutron stars represent extremes in gravity, density, and magnetic fields. They therefore serve as key objects in the study of multiple frontiers of physics. In addition, their origin (mainly in core-collapse supernovae) and evolution (via accretion or, for neutron stars, magnetic spindown and reconfiguration) touch upon multiple open issues in astrophysics. In this review, we discuss current mass and spin measurements and their reliability for neutron stars and stellar-mass black holes, as well as the overall importance of spins and masses for compact object astrophysics. Current masses are obtained primarily through electromagnetic observations of binaries, although future microlensing observations promise to enhance our understanding substantially. The spins of neutron stars are straightforward to measure for pulsars, but the birth spins of neutron stars are more difficult to determine. In contrast, even the current spins of stellar-mass black holes are challenging to measure. As we discuss, major inroads have been made in black hole spin estimates via analysis of iron lines and continuum emission, with reasonable agreement when both types of estimate are possible for individual objects, and future X-ray polarization measurements may provide additional independent information. We conclude by exploring the exciting prospects for mass and spin measurements from future gravitational wave detections, which are expected to revolutionize our understanding of strong gravity and compact objects.


title: Global existence of weak solutions of the nematic liquid crystal flow in dimensions three

abstract: For any bounded smooth domain $\Omega\subset\mathbb R^3$, we establish the global existence of a weak solution $u:\Omega\times (0,+\infty)\to\mathbb R^3\times\mathbb S^2$ of the initial-boundary value (or the Cauchy) problem of the simplified Ericksen-Leslie system (1.1) modeling the hydrodynamic flow of nematic liquid crystals for any initial and boundary (or Cauchy) data $(u_0. d_0)\in {\bf H}\times H^1(\Omega,\mathbb S^2$), with $d_0(\Omega)\subset\mathbb S^2_+$ (the upper hemisphere). Furthermore, ($u,d$) satisfies the global energy inequality (1.4).


title: Analysis of graded-index optical fibers by the spectral parameter power series method

abstract: Spectral parameter power series (SPPS) method is a recently introduced technique for solving linear differential equations and related spectral problems. In the present work we develop an approach based on the SPPS for analysis of graded-index optical fibers. The characteristic equation of the eigenvalue problem for calculation of guided modes is obtained in an analytical form in terms of SPPS. Truncation of the series and consideration in this way of the approximate characteristic equation gives us a simple and efficient numerical method for solving the problem. Comparison with the results obtained by other available techniques reveals clear advantages of the SPPS approach, in particular, with regards to accuracy. Based on the solution of the eigenvalue problem, parameters describing the dispersion are analyzed as well.


title: Lattices with many Borcherds products

abstract: We prove that there are only finitely many isometry classes of even lattices $L$ of signature $(2,n)$ for which the space of cusp forms of weight $1+n/2$ for the Weil representation of the discriminant group of $L$ is trivial. We compute the list of these lattices. They have the property that every Heegner divisor for the orthogonal group of $L$ can be realized as the divisor of a Borcherds product. We obtain similar classification results in greater generality for finite quadratic modules.


title: Global finite energy weak solutions to the compressible nematic liquid crystal flow in dimension three

abstract: In this paper, we consider the initial and boundary value problem of a simplified compressible nematic liquid crystal flow in $\Omega\subset\mathbb R^3$. We establish the existence of global weak solutions, provided the initial orientational director field $d_0$ lies in the hemisphere $\mathbb S^2_+$.


title: On the Inclination and Habitability of the HD 10180 System

abstract: There are numerous multi-planet systems that have now been detected via a variety of techniques. These systems exhibit a range of both planetary properties and orbital configurations. For those systems without detected planetary transits, a significant unknown factor is the orbital inclination. This produces an uncertainty in the mass of the planets and their related properties, such as atmospheric scale height. Here we investigate the HD~10180 system which was discovered using the radial velocity technique. We provide a new orbital solution for the system which allows for eccentric orbits for all planets. We show how the inclination of the system affects the mass/radius properties of the planets and how the detection of phase signatures may resolve the inclination ambiguity. We finally evaluate the Habitable Zone properties of the system and show that the g planet spends 100\% of an eccentric orbit within the Habitable Zone.


title: A Price Selective Centralized Algorithm for Resource Allocation with Carrier Aggregation in LTE Cellular Networks

abstract: In this paper, we consider a resource allocation with carrier aggregation optimization problem in long term evolution (LTE) cellular networks. In our proposed model, users are running elastic or inelastic traffic. Each user equipment (UE) is assigned an application utility function based on the type of its application. Our objective is to allocate multiple carriers resources optimally among users in their coverage area while giving the user the ability to select one of the carriers to be its primary carrier and the others to be its secondary carriers. The UE's decision is based on the carrier price per unit bandwidth. We present a price selective centralized resource allocation with carrier aggregation algorithm to allocate multiple carriers resources optimally among users while providing a minimum price for the allocated resources. In addition, we analyze the convergence of the algorithm with different carriers rates. Finally, we present simulation results for the performance of the proposed algorithm.


title: SNe Ia: Can Chandrasekhar Mass Explosions Reproduce the Observed Zoo?

abstract: The question of the nature of the progenitor of Type Ia supernovae (SNe Ia) is important both for our detailed understanding of stellar evolution and for their use as cosmological probes of the dark energy. Much of the basic features of SNe Ia can be understood directly from the nuclear physics, a fact which Gerry would have appreciated. We present an overview of the current observational and theoretical situation and show that it not incompatible with most SNe Ia being the results of thermonuclear explosions near the Chandrasekhar mass.


title: Central limit theorems, Lee-Yang zeros, and graph-counting polynomials

abstract: We consider the asymptotic normalcy of families of random variables $X$ which count the number of occupied sites in some large set. We write $Prob(X=m)=p_mz_0^m/P(z_0)$, where $P(z)$ is the generating function $P(z)=\sum_{j=0}^{N}p_jz^j$ and $z_0>0$. We give sufficient criteria, involving the location of the zeros of $P(z)$, for these families to satisfy a central limit theorem (CLT) and even a local CLT (LCLT); the theorems hold in the sense of estimates valid for large $N$ (we assume that $Var(X)$ is large when $N$ is). For example, if all the zeros lie in the closed left half plane then $X$ is asymptotically normal, and when the zeros satisfy some additional conditions then $X$ satisfies an LCLT. We apply these results to cases in which $X$ counts the number of edges in the (random) set of "occupied" edges in a graph, with constraints on the number of occupied edges attached to a given vertex. Our results also apply to systems of interacting particles, with $X$ counting the number of particles in a box $\Lambda$ whose size approaches infinity; $P(z)$ is then the grand canonical partition function and its zeros are the Lee-Yang zeros.

categories: math.CO cond-mat.stat-mech math-ph math.MP


title: A non-local rare mutations model for quasispecies and Prisoner's dilemma: numerical assessment of qualitative behaviour

abstract: An integro-differential model for evolutionary dynamics with mutations is investigated by improving the understanding of its behavior using numerical simulations. The proposed numerical approach can handle also density dependent fitness, and gives new insights about the role of mutation in the preservation of cooperation.


title: Harnack Estimates for Conjugate Heat Kernel on Evolving Manifolds

abstract: In this article we derive Harnack estimates for conjugate heat kernel in an abstract geometric flow. Our calculation involves a correction term D. When D is nonnegative, we are able to obtain a Harnack inequality. Our abstract formulation provides a unified framework for some known results, in particular including corresponding results of Ni, Perelman, and Tran as special cases. Moreover it leads to new results in the setting of Ricci-Harmonic flow and mean curvature flow in Lorentzian manifolds with nonnegative sectional curvature.


title: Efficient Online Strategies for Renting Servers in the Cloud

abstract: In Cloud systems, we often deal with jobs that arrive and depart in an online manner. Upon its arrival, a job should be assigned to a server. Each job has a size which defines the amount of resources that it needs. Servers have uniform capacity and, at all times, the total size of jobs assigned to a server should not exceed the capacity. This setting is closely related to the classic bin packing problem. The difference is that, in bin packing, the objective is to minimize the total number of used servers. In the Cloud, however, the charge for each server is proportional to the length of the time interval it is rented for, and the goal is to minimize the cost involved in renting all used servers. Recently, certain bin packing strategies were considered for renting servers in the Cloud [Li et al. SPAA'14]. There, it is proved that all Any-Fit bin packing strategy has a competitive ratio of at least $\mu$, where $\mu$ is the max/min interval length ratio of jobs. It is also shown that First Fit has a competitive ratio of $2\mu + 13$ while Best Fit is not competitive at all. We observe that the lower bound of $\mu$ extends to all online algorithms. We also prove that, surprisingly, Next Fit algorithm has competitive ratio of at most $2 \mu +1$. We also show that a variant of Next Fit achieves a competitive ratio of $K \times max\{1,\mu/(K-1)\}+1$, where $K$ is a parameter of the algorithm. In particular, if the value of $\mu$ is known, the algorithm has a competitive ratio of $\mu+2$; this improves upon the existing upper bound of $\mu+8$. Finally, we introduce a simple algorithm called Move To Front (MTF) which has a competitive ratio of at most $6\mu + 7$ and also promising average-case performance. We experimentally study the average-case performance of different algorithms and observe that the typical behaviour of MTF is distinctively better than other algorithms.


title: Compressive Sampling of Polynomial Chaos Expansions: Convergence Analysis and Sampling Strategies

abstract: Sampling orthogonal polynomial bases via Monte Carlo is of interest for uncertainty quantification of models with high-dimensional random inputs, using Polynomial Chaos (PC) expansions. It is known that bounding a probabilistic parameter, referred to as {\it coherence}, yields a bound on the number of samples necessary to identify coefficients in a sparse PC expansion via solution to an $\ell_1$-minimization problem. Utilizing asymptotic results for orthogonal polynomials, we bound the coherence parameter for polynomials of Hermite and Legendre type under the respective natural sampling distribution. In both polynomial bases we identify an importance sampling distribution which yields a bound with weaker dependence on the order of the approximation. For more general orthonormal bases, we propose the {\it coherence-optimal} sampling: a Markov Chain Monte Carlo sampling, which directly uses the basis functions under consideration to achieve a statistical optimality among all sampling schemes with identical support. We demonstrate these different sampling strategies numerically in both high-order and high-dimensional, manufactured PC expansions. In addition, the quality of each sampling method is compared in the identification of solutions to two differential equations, one with a high-dimensional random input and the other with a high-order PC expansion. In both cases the coherence-optimal sampling scheme leads to similar or considerably improved accuracy.


title: Sparse and compositionally robust inference of microbial ecological networks

abstract: 16S-ribosomal sequencing and other metagonomic techniques provide snapshots of microbial communities, revealing phylogeny and the abundances of microbial populations across diverse ecosystems. While changes in microbial community structure are demonstrably associated with certain environmental conditions, identification of underlying mechanisms requires new statistical tools, as these datasets present several technical challenges. First, the abundances of microbial operational taxonomic units (OTUs) from 16S datasets are compositional, and thus, microbial abundances are not independent. Secondly, microbial sequencing-based studies typically measure hundreds of OTUs on only tens to hundreds of samples; thus, inference of OTU-OTU interaction networks is severely under-powered, and additional assumptions are required for accurate inference. Here, we present SPIEC-EASI (SParse InversE Covariance Estimation for Ecological Association Inference), a statistical method for the inference of microbial ecological interactions from metagenomic datasets that addresses both of these issues. SPIEC-EASI combines data transformations developed for compositional data analysis with a graphical model inference framework that assumes the underlying ecological interaction network is sparse. To reconstruct the interaction network, SPIEC-EASI relies on algorithms for sparse neighborhood and inverse covariance selection. Because no large-scale microbial ecological networks have been experimentally validated, SPIEC-EASI comprises computational tools to generate realistic OTU count data from a set of diverse underlying network topologies. SPIEC-EASI outperforms state-of-the-art methods in terms of edge recovery and network properties on realistic synthetic data under a variety of scenarios. SPIEC-EASI also reproducibly predicts previously unknown microbial interactions using data from the American Gut project.


title: Building relativistic mean field models for finite nuclei and neutron stars

abstract: Background: Theoretical approaches based on density functional theory provide the only tractable method to incorporate the wide range of densities and isospin asymmetries required to describe finite nuclei, infinite nuclear matter, and neutron stars. Purpose: A relativistic energy density functional (EDF) is developed to address the complexity of such diverse nuclear systems. Moreover, a statistical perspective is adopted to describe the information content of various physical observables. Methods: We implement the model optimization by minimizing a suitably constructed chi-square objective function using various properties of finite nuclei and neutron stars. The minimization is then supplemented by a covariance analysis that includes both uncertainty estimates and correlation coefficients. Results: A new model, FSUGold2, is created that can well reproduce the ground-state properties of finite nuclei, their monopole response, and that accounts for the maximum neutron star mass observed up to date. In particular, the model predicts both a stiff symmetry energy and a soft equation of state for symmetric nuclear matter--suggesting a fairly large neutron-skin thickness in Pb208 and a moderate value of the nuclear incompressibility. Conclusions: We conclude that without any meaningful constraint on the isovector sector, relativistic EDFs will continue to predict significantly large neutron skins. However, the calibration scheme adopted here is flexible enough to create models with different assumptions on various observables. Such a scheme--properly supplemented by a covariance analysis--provides a powerful tool to identify the critical measurements required to place meaningful constraints on theoretical models.


title: The Deep Diffuse Extragalactic Radio Sky at 1.75 GHz

abstract: We present a study of diffuse extragalactic radio emission at $1.75\,$GHz from part of the ELAIS-S1 field using the Australia Telescope Compact Array. The resulting mosaic is $2.46\,$deg$^2$, with a roughly constant noise region of $0.61\,$deg$^2$ used for analysis. The image has a beam size of $150 \times60\,$arcsec and instrumental $\langle\sigma_{\rm n}\rangle= (52\pm5)\, \mu$Jy beam$^{-1}$. Using point-source models from the ATLAS survey, we subtract the discrete emission in this field for $S \ge 150\, \mu$Jy beam$^{-1}$. Comparison of the source-subtracted probability distribution, or \pd, with the predicted distribution from unsubtracted discrete emission and noise, yields an excess of $(76 \pm 23) \, \mu$Jy beam$^{-1}$. Taking this as an upper limit on any extended emission we constrain several models of extended source counts, assuming $\Omega_{\rm source} \le 2\,$arcmin. The best-fitting models yield temperatures of the radio background from extended emission of $T_{\rm b}=(10\pm7) \,$mK, giving an upper limit on the total temperature at $1.75\,$GHz of $(73\pm10)\,$mK. Further modelling shows that our data are inconsistent with the reported excess temperature of ARCADE2 to a source-count limit of $1\, \mu$Jy. Our new data close a loop-hole in the previous constraints, because of the possibility of extended emission being resolved out at higher resolution. Additionally, we look at a model of cluster halo emission and two WIMP dark matter annihilation source-count models, and discuss general constraints on any predicted counts from such sources. Finally, we report the derived integral count at $1.4\,$GHz using the deepest discrete count plus our new extended-emission limits, providing numbers that can be used for planning future ultra-deep surveys.


title: The Biermann Catastrophe in Numerical MHD

abstract: The Biermann Battery effect is frequently invoked in cosmic magnetogenesis and studied in High-Energy Density laboratory physics experiments. Generation of magnetic fields by the Biermann effect due to mis-aligned density and temperature gradients in smooth flow behind shocks is well known. We show that a Biermann-effect magnetic field is also generated within shocks. Direct implementation of the Biermann effect in MHD codes does not capture this physical process, and worse, produces unphysical magnetic fields at shocks whose value does not converge with resolution. We show that this convergence breakdown is due to naive discretization, which fails to account for the fact that discretized irrotational vector fields have spurious solenoidal components that grow without bound near a discontinuity. We show that careful consideration of the kinetics of ion viscous shocks leads to a formulation of the Biermann effect that gives rise to a convergent algorithm. We note two novel physical effects: a resistive magnetic precursor in which Biermann-generated field in the shock "leaks" resistively upstream; and a thermal magnetic precursor, in which field is generated by the Biermann effect ahead of the shock front due to gradients created by the shock's electron thermal conduction precursor. Both effects appear to be potentially observable in experiments at laser facilities. We re-examine published studies of magnetogenesis in galaxy cluster formation, and conclude that the simulations in question had inadequate resolution to reliably estimate the field generation rate. Corrected estimates suggest primordial field values in the range $B\sim 10^{-22}$G --- $10^{-19}$G by $z=3$.


title: Optimal factorizations of rational numbers using factorization trees

abstract: Let $m_t(\alpha)$ denote the $t$-metric Mahler measure of the algebraic number $\alpha$. Recent work of the first author established that the infimum in $m_t(\alpha)$ is attained by a single point $\bar\alpha = (\alpha_1,\ldots,\alpha_N)\in \overline{\mathbb Q}^N$ for all sufficiently large $t$. Nevertheless, no efficient method for locating $\bar \alpha$ is known. In this article, we define a new tree data structure, called a factorization tree, which enables us to find $\bar\alpha$ when $\alpha\in \mathbb Q$. We establish several basic properties of factorization trees, and use these properties to locate $\bar\alpha$ in previously unknown cases.


title: The effect of spatial resolution on optical and near-IR studies of stellar clusters: Implications for the origin of the red excess

abstract: Recent ground based near-IR studies of stellar clusters in nearby galaxies have suggested that young clusters remain embedded for 7-10Myr in their progenitor molecular cloud, in conflict with optical based studies which find that clusters are exposed after 1-3Myr. Here, we investigate the role that spatial resolution plays in this apparent conflict. We use a recent catalogue of young ($<10$~Myr) massive ($>5000$~\msun) clusters in the nearby spiral galaxy, M83, along with Hubble Space Telescope (HST) imaging in the optical and near-IR, and ground based near-IR imaging, to see how the colours (and hence estimated properties such as age and extinction) are affected by the aperture size employed, in order to simulate studies of differing resolution. We find that the near-IR is heavily affected by the resolution, and when aperture sizes $>40$~pc are used, all young/blue clusters move red-ward in colour space, which results in their appearance as heavily extincted clusters. However, this is due to contamination from nearby sources and nebular emission, and is not an extinction effect. Optical colours are much less affected by resolution. Due to the larger affect of contamination in the near-IR, we find that, in some cases, clusters will appear to show near-IR excess when large ($>20$~pc) apertures are used. Our results explain why few young ($<6$~Myr), low extinction ($\av < 1$~mag) clusters have been found in recent ground based near-IR studies of cluster populations, while many such clusters have been found in higher resolution HST based studies. Additionally, resolution effects appear to (at least partially) explain the origin of the near-IR excess that has been found in a number of extragalactic YMCs.


title: The generic Green-Lazarsfeld secant conjecture

abstract: Generalizing the well-known Green Conjecture on syzygies of canonical curves, Green and Lazarsfeld formulated in 1986 the Secant Conjecture predicting that a line bundle L of sufficiently high degree on a curve has a non-linear p-syzygy if and only if L fails to be (p+1)-very ample. Via lattice theory for special K3 surfaces, Voisin's solution of the classical Green Conjecture and calculations on moduli stacks of pointed curves, we prove: (1) The Green-Lazarsfeld Secant Conjecture in various degree of generality, including its strongest possible form in the divisorial case in the universal Jacobian. (2) The Prym-Green Conjecture on the naturality of the resolution of a general Prym-canonical curve of odd genus.


title: The infimum in the metric Mahler measure

abstract: Dubickas and Smyth defined the metric Mahler measure on the multiplicative group of non-zero algebraic numbers. The definition involves taking an infimum over representations of an algebraic number $\alpha$ by other algebraic numbers. We verify their conjecture that the infimum in its definition is always achieved as well as establish its analog for the ultrametric Mahler measure.


title: The $t$-metric Mahler measures of surds and rational numbers

abstract: A. Dubickas and C. Smyth introduced the metric Mahler measure $$ M_1(\alpha) = \inf\left\{\sum_{n=1}^N M(\alpha_n): N \in \mathbb N, \alpha_1 \cdots \alpha_N = \alpha\right\}, $$ where $M(\alpha)$ denotes the usual (logarithmic) Mahler measure of $\alpha \in \overline{\mathbb Q}$. This definition extends in a natural way to the $t$-metric Mahler measure by replacing the sum with the usual $L_t$ norm of the vector $(M(\alpha_1), \dots, M(\alpha_N))$ for any $t\geq 1$. For $\alpha \in \mathbb Q$, we prove that the infimum in $M_t(\alpha)$ may be attained using only rational points, establishing an earlier conjecture of the second author. We show that the natural analogue of this result fails for general $\alpha\in\overline{\mathbb Q}$ by giving an infinite family of quadratic counterexamples. As part of this construction, we provide an explicit formula to compute $M_t(D^{1/k})$ for a square-free $D \in \mathbb N$.


title: Estimating heights using auxiliary functions

abstract: Several recent papers construct auxiliary polynomials to bound the Weil height of certain classes of algebraic numbers from below. Following these techniques, the author gave a general method for introducing auxiliary polynomials to problems involving the Weil height. The height appears as a solution to a certain extremal problem involving polynomials. We further generalize the above techniques to acquire both the projective height and the height on subspaces in the same way. We further obtain lower bounds on the heights of points on some subvarieties of $\mathbb P^{N-1}(\overline{\mathbb Q})$.


title: To which densities is spin-polarized neutron matter a weakly interacting Fermi gas?

abstract: We study the properties of spin-polarized neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory, including two-, three-, and four-neutron interactions. The energy of spin-polarized neutrons is remarkably close to a non-interacting system at least up to saturation density, where interaction effects provide less than 10% corrections. This shows that the physics of neutron matter is similar to a unitary gas well beyond the scattering-length regime. Implications for energy-density functionals and for a possible ferromagnetic transition in neutron stars are discussed. Our predictions can be tested with lattice QCD, and we present results for varying pion mass.


title: High-precision quark masses and QCD coupling from $n_f=4$ lattice QCD

abstract: We present a new lattice QCD analysis of heavy-quark pseudoscalar-pseudoscalar correlators, using gluon configurations from the MILC collaboration that include vacuum polarization from $u$, $d$, $s$ and $c$~quarks($n_f=4$). We extract new values for the QCD coupling and for the $c$ quark's $\overline{\mathrm{MS}}$ mass: $\alpha_{\overline{\mathrm{MS}}}(M_Z,n_f=5) = 0.11822(74)$ and $m_c(3\mathrm{GeV}, n_f=4) = 0.9851(63)$GeV. These agree well with our earlier simulations using $n_f=3$ sea quarks, vindicating the perturbative treatment of $c$ quarks in that analysis. We also obtain a new nonperturbative result for the ratio of $c$~and $s$~quark masses: $m_c/m_s=11.652(65)$. This ratio implies $m_s(2\,\mathrm{GeV}, n_f=3)=93.6(8)$MeV when it is combined with our new~$c$~mass. Combining $m_c/m_s$ with our earlier $m_b/m_c$ gives $m_b/m_s=52.55(55)$, which is several standard deviations (but only 4%) away from the Georgi-Jarlskop prediction from certain GUTs. Finally we obtain an $n_f=4$ estimate for $m_b/m_c=4.528(54)$ which agrees well with our earlier $n_f=3$ result. The new ratio implies~$m_b(m_b,n_f=5)=4.162(48)$GeV.


title: Mixing time of A + B -> 0 in 1D

abstract: A mixing time density of $A + B \to 0$ on a finite one dimensional domain is defined for general initial and boundary conditions in which $A$ and $B$ diffuse at the same rate. The density is a measure of the number of $A$ and $B$ particles that mix through the center of the reaction zone. It also corresponds to the reaction density for the special case in which $A$ and $B$ annihilate upon contact. An exact expression is found for the generating function of the mixing time. The analysis is extended to multiple reaction fronts and finitely ramified fractals. The full method involves using the kernel of the Laplace transform integral operator to map and analyze a moving homogeneous Dirichlet interior point condition.


title: Solitary Waves in a Discrete Nonlinear Dirac equation

abstract: In the present work, we introduce a discrete formulation of the nonlinear Dirac equation in the form of a discretization of the Gross-Neveu model. The motivation for this discrete model proposal is both computational (near the continuum limit) and theoretical (using the understanding of the anti-continuum limit of vanishing coupling). Numerous unexpected features are identified including a staggered solitary pattern emerging from a single site excitation, as well as two- and three-site excitations playing a role analogous to one- and two-site, respectively, excitations of the discrete nonlinear Schr\"odinger analogue of the model. Stability exchanges between the two- and three-site states are identified, as well as instabilities that appear to be persistent over the coupling strength $\epsilon$, for a subcritical value of the propagation constant $\Lambda$. Variations of the propagation constant, coupling parameter and nonlinearity exponent are all examined in terms of their existence and stability implications and long dynamical simulations are used to unravel the evolutionary phenomenology of the system (when unstable).


title: Tunneling Conductance and Spin Transport in Clean Ferromagnet-Ferromagnet-Superconductor Heterostructures

abstract: We present a transfer matrix approach that combines the Blonder-Tinkham-Klapwijk (BTK) formalism and self-consistent solutions to the Bogolibuov-de Gennes (BdG) equations and use it to study the tunneling conductance and spin transport in ferromagnet (${\rm F}$)-superconductor (${\rm S}$) trilayers (${\rm F_1F_2 S}$) as functions of bias voltage. The self-consistency ensures that the spin and charge conservation laws are properly satisfied. We consider forward and angularly averaged conductances over a broad range of the strength of the exchange fields and ${\rm F}$ thicknesses, as the relative in-plane magnetization angle, $\phi$, between the two ferromagnets varies. The $\phi$-dependence of the self-consistent conductance curves in the trilayers can differ substantially from that obtained via a non-self-consistent approach. The zero bias forward conductance peak exhibits, as $\phi$ varies, resonance effects intricately associated with particular combinations of the geometrical and material parameters. We find, when the magnetizations are non-collinear, signatures of the anomalous Andreev reflections in the subgap regions of the angularly averaged conductances. When ${\rm F_1}$ is half-metallic, the angularly averaged subgap conductance chiefly arises from anomalous Andreev reflection. The in-plane components of the spin current are strongly bias dependent, while the out-of-plane spin current component is only weakly dependent upon voltage. The components of the spin current aligned with the local exchange field of one of the F layers are conserved in that layer and in the S region, while they oscillate in the other layer. We compute the spin transfer torques, in connection with the oscillatory behavior of spin currents, and verify that the spin continuity equation is strictly obeyed in our method.


title: Magnetohydrodynamics and charged currents in heavy ion collisions

abstract: The hot QCD matter produced in any heavy ion collision with a nonzero impact parameter is produced within a strong magnetic field. We study the imprint the magnetic fields produced in non-central heavy ion collisions leave on the azimuthal distributions and correlations of the produced charged hadrons. The magnetic field is time-dependent and the medium is expanding, which leads to the induction of charged currents due to the combination of Faraday and Hall effects. We find that these currents result in a charge-dependent directed flow $v_1$ that is odd in rapidity and odd under charge exchange. It can be detected by measuring correlations between the directed flow of charged hadrons at different rapidities, $\langle v_1^\pm (y_1) v_1^\pm (y_2) \rangle$.


title: Resonance Conversion as a Catalyser of Nuclear Reactions

abstract: It is shown that resonance interal conversion offers a feasible tool for mastering nuclear processes with laser or synchrotron radiation. Physics of the process is discussed in detail in historical aspect. Possible way of experimental applicaytion is shown in the case of the $M1$ 70.6-keV transition in nuclei of $^{169}$Yb. Nuclear transition rate in hydrogenlike ions of this nuclide can be enhanced by up to four orders of magnitude.


title: Improved V II log($gf$) Values, Hyperfine Structure Constants, and Abundance Determinations in the Photospheres of the Sun and Metal-poor Star HD 84937

abstract: New experimental absolute atomic transition probabilities are reported for 203 lines of V II. Branching fractions are measured from spectra recorded using a Fourier transform spectrometer and an echelle spectrometer. The branching fractions are normalized with radiative lifetime measurements to determine the new transition probabilities. Generally good agreement is found between this work and previously reported V II transition probabilities. Use of two spectrometers, independent radiometric calibration methods, and independent data analysis routines enables a reduction in systematic uncertainties, in particular those due to optical depth errors. In addition, new hyperfine structure constants are measured for selected levels by least squares fitting line profiles in the FTS spectra. The new V II data are applied to high resolution visible and UV spectra of the Sun and metal-poor star HD 84937 to determine new, more accurate V abundances. Lines covering a range of wavelength and excitation potential are used to search for non-LTE effects. Very good agreement is found between our new solar photospheric V abundance, log {\epsilon}(V) = 3.95 from 15 V II lines, and the solar-system meteoritic value. In HD 84937, we derive [V/H] = -2.08 from 68 lines, leading to a value of [V/Fe] = 0.24.


title: Linear and Nonlinear Modeling of a Traveling-Wave Thermoacoustic Heat Engine

abstract: We have carried out three-dimensional Navier-Stokes simulations, from quiescent conditions to the limit cycle, of a traveling-wave thermoacoustic heat engine (TAE) composed of a long variable-area resonator shrouding a smaller annular tube, which encloses the hot (HHX) and ambient (AHX) heat-exchangers, and the regenerator (REG). Simulations are wall-resolved, with no-slip and adiabatic conditions enforced at all boundaries, while the heat transfer and drag due to the REG and HXs are modeled. HHX temperatures have been investigated in the range 440K - 500K with AHX temperature fixed at 300K. The initial exponential growth of acoustic energy is due to a network of traveling waves amplified by looping around the REG/HX unit in the direction of the imposed temperature gradient. A simple analytical model demonstrates that such thermoacoustic instability is a Lagrangian thermodynamic process resembling a Stirling cycle. A system-wide linear stability model based on Rott's theory is able to accurately predict the frequency and the growth rate, exhibiting properties consistent with a supercritical Hopf bifurcation. The limit cycle is governed by acoustic streaming -- a rectified flow resulting from high-amplitude acoustics. Its key features are explained with an axially symmetric incompressible model driven by the wave-induced stresses extracted from the compressible calculations. These features include Gedeon streaming and mean recirculations due to flow separation. The first drives the mean advection of hot fluid from the HHX to a secondary heat exchanger (AHX2), in the thermal buffer tube (TBT), necessary to achieve saturation of the acoustic energy growth. The direct evaluation of the energy fluxes reveals that the efficiency of the device deteriorates with the drive ratio and that the acoustic power in the TBT is balanced primarily by the mean advection and thermoacoustic heat transport.


title: Radio continuum properties of luminous infrared galaxies. Identifying the presence of an AGN in the radio

abstract: Luminous infrared galaxies are systems enshrouded in dust, which absorbs most of their optical/UV emission and re-radiates it in the mid- and far-infrared. Radio observations are largely unaffected by dust obscuration, enabling us to study the central regions of LIRGs in an unbiased manner. The main goal of this project is to examine how the radio properties of local LIRGs relate to their infrared spectral characteristics. Here we present an analysis of the radio continuum properties of a subset of the Great Observatories All-sky LIRG Survey (GOALS), which consists of 202 nearby systems (z<0.088). Our radio sample consists of 35 systems, or 46 individual galaxies, that were observed at both 1.49 and 8.44 GHz with the VLA with a resolution of about 1 arcsec (FWHM). The aim of the project is to use the radio imagery to probe the central kpc of these LIRGs in search of active galactic nuclei. We used the archival data at 1.49 and 8.44 GHz to create radio-spectral-index maps using the standard relation between flux density Sv and frequency v, S~v^-a, where a is the radio spectral index. By studying the spatial variations in a, we classified the objects as radio-AGN, radio-SB, and AGN/SB (a mixture). We identified the presence of an active nucleus using the radio morphology, deviations from the radio/infrared correlation, and spatially resolved spectral index maps, and then correlated this to the usual mid-infrared ([NeV]/[NeII] and [OIV]/[NeII] line ratios and EQW of the 6.2 um PAH feature) and optical (BPT diagram) AGN diagnostics. We find that 21 out of the 46 objects in our sample are radio-AGN, 9 are classified as starbursts (SB), and 16 are AGN/SB. After comparing to other AGN diagnostics we find 3 objects out of the 46 that are identified as AGN based on the radio analysis, but are not classified as such based on the mid-infrared and optical AGN diagnostics presented in this study.


title: On the Two-user Multi-carrier Joint Channel Selection and Power Control Game

abstract: In this paper, we propose a hierarchical game approach to model the energy efficiency maximization problem where transmitters individually choose their channel assignment and power control. We conduct a thorough analysis of the existence, uniqueness and characterization of the Stackelberg equilibrium. Interestingly, we formally show that a spectrum orthogonalization naturally occurs when users decide sequentially about their transmitting carriers and powers, delivering a binary channel assignment. Both analytical and simulation results are provided for assessing and improving the performances in terms of energy efficiency and spectrum utilization between the simultaneous-move game (with synchronous decision makers), the social welfare (in a centralized manner) and the proposed Stackelberg (hierarchical) game. For the first time, we provide tight closed-form bounds on the spectral efficiency of such a model, including correlation across carriers and users. We show that the spectrum orthogonalization capability induced by the proposed hierarchical game model enables the wireless network to achieve the spectral efficiency improvement while still enjoying a high energy efficiency.


title: Note on Soft Graviton theorem by KLT Relation

abstract: Recently, new soft graviton theorem proposed by Cachazo and Strominger has inspired a lot of works. In this note, we use the KLT-formula to investigate the theorem. We have shown how the soft behavior of color ordered Yang-Mills amplitudes can be combined with KLT relation to give the soft behavior of gravity amplitudes. As a byproduct, we find two nontrivial identities of the KLT momentum kernel must hold.


title: Dissociable changes in functional network topology underlie early category learning and development of automaticity

abstract: Recent work has shown that multimodal association areas-including frontal, temporal and parietal cortex-are focal points of functional network reconfiguration during human learning and performance of cognitive tasks. On the other hand, neurocomputational theories of category learning suggest that the basal ganglia and related subcortical structures are focal points of functional network reconfiguration during early learning of some categorization tasks, but become less so with the development of automatic categorization performance. Using a combination of network science and multilevel regression, we explore how changes in the connectivity of small brain regions can predict behavioral changes during training in a visual categorization task. We find that initial category learning, as indexed by changes in accuracy, is predicted by increasingly efficient integrative processing in subcortical areas, with higher functional specialization, more efficient integration across modules, but a lower cost in terms of redundancy of information processing. The development of automaticity, as indexed by changes in the speed of correct responses, was predicted by lower clustering (particularly in subcortical areas), higher strength (highest in cortical areas) and higher betweenness centrality. By combining neurocomputational theories and network scientific methods, these results synthesize the dissociative roles of multimodal association areas and subcortical structures in the development of automaticity during category learning.


title: Accounting for inertia effects to access the high-frequency microrheology of viscoelastic fluids

abstract: We study the Brownian motion of microbeads immersed in water and in a viscoelastic wormlike micelles solution by optical trapping interferometry and diffusing wave spectroscopy. Through the mean-square displacement obtained from both techniques, we deduce the mechanical properties of the fluids at high frequencies by explicitly accounting for inertia effects of the particle and the surrounding fluid at short time scales. For wormlike micelle solutions, we recover the 3/4 scaling exponent for the loss modulus over two decades in frequency as predicted by the theory for semiflexible polymers.


title: The Szeg\"o kernel on a class of noncompact CR manifolds of high codimension

abstract: We generalize Nagel's formula for the Szeg\"o kernel and use it to compute the Szeg\"o kernel on a class of noncompact CR manifolds whose tangent space decomposes into one complex direction and several totally real directions. We also discuss the control metric on these manifolds and relate it to the size of the Szeg\"o kernel.


title: Identifying transitions in finite systems by means of partition function zeros and microcanonical inflection-point analysis: A comparison for elastic flexible polymers

abstract: For the estimation of transition points of finite elastic, flexible polymers with chain lengths from $13$ to $309$ monomers, we compare systematically transition temperatures obtained by the Fisher partition function zeros approach with recent results from microcanonical inflection-point analysis. These methods rely on accurate numerical estimates of the density of states, which have been obtained by advanced multicanonical Monte Carlo sampling techniques. Both the Fisher zeros method and microcanonical inflection-point analysis yield very similar results and enable the unique identification of transition points in finite systems, which is typically impossible in the conventional canonical analysis of thermodynamic quantities.

categories: physics.data-an cond-mat.soft cond-mat.stat-mech physics.comp-ph


title: Quadratic diameter bounds for dual network flow polyhedra

abstract: Both the combinatorial and the circuit diameters of polyhedra are of interest to the theory of linear programming for their intimate connection to a best-case performance of linear programming algorithms. We study the diameters of dual network flow polyhedra associated to $b$-flows on directed graphs $G=(V,E)$ and prove quadratic upper bounds for both of them: the minimum of $(|V|-1)\cdot |E|$ and $\frac{1}{6}|V|^3$ for the combinatorial diameter, and $\frac{|V|\cdot (|V|-1)}{2}$ for the circuit diameter. The latter strengthens the cubic bound implied by a result in [De Loera, Hemmecke, Lee; 2014]. Previously, bounds on these diameters have only been known for bipartite graphs. The situation is much more involved for general graphs. In particular, we construct a family of dual network flow polyhedra with members that violate the circuit diameter bound for bipartite graphs by an arbitrary additive constant. Further, it provides examples of circuit diameter $\frac{4}{3}|V| - 4$.


title: Bosonic Ghosts at $c=2$ as a Logarithmic CFT

abstract: Motivated by Wakimoto free field realisations, the bosonic ghost system of central charge $c=2$ is studied using a recently proposed formalism for logarithmic conformal field theories. This formalism addresses the modular properties of the theory with the aim being to determine the (Grothendieck) fusion coefficients from a variant of the Verlinde formula. The key insight, in the case of bosonic ghosts, is to introduce a family of parabolic Verma modules which dominate the spectrum of the theory. The results include S-transformation formulae for characters, non-negative integer Verlinde coefficients, and a family of modular invariant partition functions. The logarithmic nature of the corresponding ghost theories is explicitly verified using the Nahm-Gaberdiel-Kausch fusion algorithm.


title: Thermalization in closed quantum systems: semiclassical approach

abstract: Thermalization in closed quantum systems can be explained either by means of the eigenstate thermalization hypothesis or the concept of canonical typicality. Both concepts are based on quantum mechanical formalism such as spectral properties of the eigenstates or entanglement between subsystems respectively. Here we study the onset of thermalization of Bose particles in a two-band double well potential using the truncated Wigner approximation. This allows us to use the familiar classical formalism to explain quantum thermalization in this system. In particular, we demonstrate that sampling of an initial quantum state plays the role of a statistical mechanical ensemble, while subsequent chaotic classical evolution turns the initial quantum state into the thermal state.


title: Closed-Form Delay-Optimal Power Control for Energy Harvesting Wireless System with Finite Energy Storage

abstract: In this paper, we consider delay-optimal power control for an energy harvesting wireless system with finite energy storage. The wireless system is powered solely by a renewable energy source with bursty data arrivals, and is characterized by a data queue and an energy queue. We consider a delay-optimal power control problem and formulate an infinite horizon average cost Markov Decision Process (MDP). To deal with the curse of dimensionality, we introduce a virtual continuous time system and derive closed-form approximate priority functions for the discrete time MDP at various operating regimes. Based on the approximation, we obtain an online power control solution which is adaptive to the channel state information as well as the data and energy queue state information. The derived power control solution has a multi-level water-filling structure, where the water level is determined jointly by the data and energy queue lengths. We show through simulations that the proposed scheme has significant performance gain compared with various baselines.


title: Small universal families of graphs on $\aleph_{\omega+1}$

abstract: We prove that it is consistent that $\aleph_\omega$ is strong limit, $2^{\aleph_\omega}$ is large and the universality number for graphs on $\aleph_{\omega+1}$ is small. The proof uses Prikry forcing with interleaved collapsing.


title: Modular perverse sheaves on flag varieties III: positivity conditions

abstract: We further develop the general theory of the "mixed modular derived category" introduced by the authors in a previous paper in this series. We then use it to study positivity and Q-Koszulity phenomena on flag varieties.


title: Nonlinear Decelerator for Payloads in Aerial Delivery Systems. I: Design and Testing

abstract: We study the dynamics and the optimization of the shock deceleration supported by a payload when its airborne carrier impacts the ground. We build a nonlinear elastic model for a container prototype and an elastic suspension system for the payload. We model the dynamics of this system and extract information on maximum deceleration, energy transfer between the container and payload, and energy resonant damping. We designed the system and perform lab experiments for various terminal velocities and types of grounds (cement, grass, sand water, etc.). The results are compared with the theoretical model and results are commented, including predictions for deceleration at different types of ground impact. The results can be used for aerial delivery systems, splash-down of capsules, recoveries, weather balloons, coastal surveying systems, or the new introduced goal-line technology in sport competitions.


title: Impact of top-Higgs couplings on di-Higgs production at future colliders

abstract: Measuring the Higgs-self coupling is one of the most crucial goals of the future colliders, such as the LHC Run-II and the ILC-based photon collider. Since the new physics can affects the di-Higgs production not only from the Higgs self-coupling but also from the top-Higgs coupling, we investigate the di-Higgs production in the presence of the non-standard top-Higgs coupling at the LHC and ILC-based photon collider given the recent Higgs data. Due to the changed interference behaviors of the top quark loops with itself or $W$ boson loops, we find that the cross section of di-Higgs production at the LHC-14 TeV and ILC-500 GeV can be respectively enhanced up to nearly 3 and 2 times the SM predictions within 2$\sigma$ Higgs data allowed parameter region.


title: Heavy-traffic Asymptotics of Priority Polling System with Threshold Service Policy

abstract: In this paper, by the singular-perturbation technique, we investigate the heavy-traffic behavior of a priority polling system consisting of three M/M/1 queues with threshold policy. It turns out that the scaled queue-length of the critically loaded queue is exponentially distributed, independent of that of the stable queues. In addition, the queue lengths of stable queues possess the same distributions as a priority polling system with N-policy vacation. Based on this fact, we provide the exact tail asymptotics of the vacation polling system to approximate the tail distribution of the queue lengths of the stable queues, which shows that it has the same prefactors and decay rates as the classical M/M/1 preemptive priority queues. Finally, a stochastic simulation is taken to test the results aforementioned.


title: Functional Meyer-Tanaka Formula

abstract: The functional Ito formula, firstly introduced by Bruno Dupire for continuous semimartingales, might be extended in two directions: different dynamics for the underlying process and/or weaker assumptions on the regularity of the functional. In this paper, we pursue the former type by proving the functional version of the Meyer-Tanaka Formula. Following the idea of the proof of the classical time-dependent Meyer-Tanaka formula, we study the mollification of functionals and its convergence properties. As an example, we study the running maximum and the max-martingales of Yor and Obloj.