title: Origin of Hydrogen Ionization in the 1 pc Galactic Central Region

abstract: We study a possible connection between processes of gamma-ray emission and hydrogen ionization in a few pc of central region around Sgr A*. Previous investigations showed there is a discrepancy between interpretation of gamma-ray and ionization data if gamma-rays are generated by proton-proton collisions. Here we provided analysis of processes of ionization and emission basing on analytical and numerical calculations of kinetic equations which describe processes of particle propagation and their energy losses. The origin of gamma rays could be either due to collisions of relativistic protons with the dense gas of the surrounding circumnuclear disk (CND) or bremsstrahlung and inverse Compton scattering of relativistic electrons. The hydrogen ionization in this case is produced by a low energy component of the CR spectrum. We found that if ionization is produced by protons the expected ionization rate of hydrogen in the CND is of the same order as derived from IR observations. So we do not see any discrepancy between the gamma-ray and ionization data for the hadronic model. In the case of ionization by electrons we obtained the ionization rate one order of magnitude higher than follows from the IR data. In principle, a selection between the leptonic and hadronic interpretations can be performed basing on measurements of radio and X-ray fluxes from this region because the leptonic and hadronic models give different values of the fluxes from there. We do not exclude that gamma-ray production and hydrogen ionization in the CND are due to a past activity of Sgr A* which occurred about 100 year ago. Then we hypothesize that there may be connection between a past proton eruption and a flux of hard X-rays emitted by Sgr A* hundred years ago as follows from the observed time variability of the iron line seen in the direction of GC molecular clouds.


title: Thin film detection of High Energy Materials: Optical Pumping Approach

abstract: We present our work on High Energy Material detection based on thin film of Lithium using the phenomenon of Optical Pumping. The Li atoms present in the thin film are optically pumped to one of the ground hyperfine energy levels so that they can no more absorb light from the resonant light source. Now in presence of a RF signal, which quantifies the ambient magnetic field, this polarized atomic system is again randomized thus making it reabsorb the resonant light. This gives a quantified measurement of the magnetic field surrounding the thin film detector. This is then mapped to the presence of magnetic HEM and hence the HEM are detected. Our approach in this regard starts with verifying the stability of Lithium atoms in various solvents so as to get a suitable liquid medium to form a thin film. In this regard, various UV-visible characterization spectra are presented to finally approach a stable system for the detection. We have worked on around 10 polar and non- polar solvents to see the stability criteria. UV-Vis probe has been found to be quite effective in addressing issue of hosting Lithium atoms by various solvents. Further results on the atomic system size as dispersed in the various solvents are presented. Atomic system size gives us a good quantitative estimate of the density of Lithium dispersed system in the thin film. Finally a separate atomic system under Ultrahigh vacuum is presented to achieve optical pumping wherein the polarization and depolarization of atomic system is seen using a 670 nm diode laser system which is electronically locked using a double window Hollow cathode lamp. These both systems will finally be combined together to give the final form to the detection system. As this whole system measures relative change in magnetic field due to the presence of HEM,it gives a more reliable and accurate detection method compared to the absolute measurement detectors.


title: Wave Model and Wave Drift Caused by the Asymmetry of Crest and Trough

abstract: It follows from the review on classical wave models that the asymmetry of crest and trough is the direct cause for wave drift. Based on this, a new model of Lagrangian form is constructed. Relative to the Gerstner model, its improvement is reflected in the horizontal motion which includes an explicit drift term. On the one hand, the depth-decay factor for the new drift accords well with that of the particle's horizontal velocity. It is more rational than that of Stokes drift. On the other hand, the new formula needs no Taylor expansion as for Stokes drift and is applicable for the waves with big slopes. In addition, the new formula can also yield a more rational magnitude for the surface drift than that of Stokes.


title: Exponential sums over primes in short intervals and an application to the Waring--Goldbach problem

abstract: Let $\Lambda(n)$ be the von Mangoldt function, $x$ real and $2\leq y \leq x$. This paper improves the estimate on the exponential sum over primes in short intervals \[ S_k(x,y;\alpha) = \sum_{x< n \leq x+y} \Lambda(n) e\left( n^k \alpha \right) \] when $k\geq 3$ for $\alpha$ in the minor arcs. And then combined with the Hardy--Littlewood circle method, this enables us to investigate the Waring--Goldbach problem of representing a positive integer $n$ as the sum of $s$ $k$th powers of almost equal prime numbers, which improves the results in Wei and Wooley [12].


title: Some remarks on prediction models

abstract: This article is a discussion of some characteristic properties in connection with global models, particularly for the application of prediction, such as the approximation property, the interpolation property and the transmission property.


title: Deuterium Enrichment of the Interstellar Medium

abstract: Despite low elemental abundance of atomic deuterium in interstellar medium (ISM), observational evidences suggest that several species in gas-phase and in ices could be heavily fractionated. We explore various aspects of deuterium enrichment by constructing a chemical evolution model in gas and grain phases. Depending on various physical parameters, gas and grains are allowed to interact with each other through exchange of their chemical species. It is known that HCO+ and N2H+ are two abundant gas phase ions in ISM and their deuterium fractionation are generally used to predict degree of ionization in various regions of a molecular cloud. To have a more realistic estimation, we consider a density profile of a collapsing cloud. We present radial distributions of important interstellar molecules along with their deuterated isotopomers. We carry out quantum chemical simulation to study effects of isotopic substitution on spectral properties of these important interstellar species. We calculate vibrational (harmonic) frequency of the most important deuterated species (neutral & ions). Rotational and distortional constants of these molecules are also computed to predict rotational transitions of these species. We compare vibrational (harmonic) and rotational transitions as computed by us with existing observational, experimental and theoretical results. We hope that our results would assist observers in their quest of several hitherto unobserved deuterated species.


title: Taxation as an instrument of stimulation of innovation-active business entities

abstract: The analysis of the theoretical material revealed the lack of consensus on defini-tion of the tax stimulation of innovation-active business entities within the re-gional taxation. The definition tax stimulation of innovation-active business en-tities is specified.


title: A simple method to estimate fractal dimension of mountain surfaces

abstract: Fractal surfaces are ubiquitous in nature as well as in the sciences. The examples range from the cloud boundaries to the corroded surfaces. Fractal dimension gives a measure of the irregularity in the object under study. We present a simple method to estimate the fractal dimension of mountain surface. We propose to use easily available satellite images of lakes for this purpose. The fractal dimension of the boundary of a lake, which can be extracted using image analysis softwares, can be determined easily which gives the estimate of the fractal dimension of the mountain surface and hence a quantitative characterization of the irregularity of the topography of the mountain surface. This value will be useful in validating models of mountain formation


title: Local Group galaxies emerge from the dark

abstract: The "Lambda Cold Dark Matter" (LCDM) model of cosmic structure formation is eminently falsifiable: once its parameters are fixed on large scales, it becomes testable in the nearby Universe. Observations within our Local Group of galaxies, including the satellite populations of the Milky Way and Andromeda, appear to contradict LCDM predictions: there are far fewer satellite galaxies than dark matter halos (the "missing satellites" problem), galaxies seem to avoid the largest substructures (the "too big to fail" problem), and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). We present results from the first hydrodynamic simulations of the Local Group that match the observed abundance of galaxies. We find that when baryonic and dark matter are followed simultaneously in the context of a realistic galaxy formation model, all three "problems" are resolved within the LCDM paradigm.


title: The SILCC (SImulating the LifeCycle of molecular Clouds) project: I. Chemical evolution of the supernova-driven ISM

abstract: The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We simulate the evolution of the multi-phase ISM in a 500 pc x 500 pc x 10 kpc region of a galactic disc, with a gas surface density of $\Sigma_{_{\rm GAS}} = 10 \;{\rm M}_\odot/{\rm pc}^2$. The Flash 4.1 simulations include an external potential, self-gravity, magnetic fields, heating and radiative cooling, time-dependent chemistry of H$_2$ and CO considering (self-) shielding, and supernova (SN) feedback. We explore SN explosions at different (fixed) rates in high-density regions (peak), in random locations (random), in a combination of both (mixed), or clustered in space and time (clustered). Only random or clustered models with self-gravity (which evolve similarly) are in agreement with observations. Molecular hydrogen forms in dense filaments and clumps and contributes 20% - 40% to the total mass, whereas most of the mass (55% - 75%) is in atomic hydrogen. The ionised gas contributes <10%. For high SN rates (0.5 dex above Kennicutt-Schmidt) as well as for peak and mixed driving the formation of H$_2$ is strongly suppressed. Also without self-gravity the H$_2$ fraction is significantly lower ($\sim$ 5%). Most of the volume is filled with hot gas ($\sim$90% within $\pm$2 kpc). Only for random or clustered driving, a vertically expanding warm component of atomic hydrogen indicates a fountain flow. Magnetic fields have little impact on the final disc structure. However, they affect dense gas ($n\gtrsim 10\;{\rm cm}^{-3}$) and delay H$_2$ formation. We highlight that individual chemical species, in particular atomic hydrogen, populate different ISM phases and cannot be accurately accounted for by simple temperature-/density-based phase cut-offs.


title: Radially stabilized inflating cosmic strings

abstract: In General Relativity, cosmic strings are well known to produce a static, locally flat spacetime with a wedge removed. If the tension exceeds a critical value, the deficit angle becomes larger than $ 2\pi $, leading to a compact exterior that ends in a conical singularity. In this work, we investigate dynamical solutions for cosmic strings with super-critical tensions. To this end, we model the string as a cylindrical shell of finite and stabilized transverse width and show that there is a marginally super-critical regime in which the stabilization can be achieved by physically reasonable matter. We show numerically that the static deficit angle solution is unstable for super-critical string tensions. Instead, the geometry starts expanding in axial direction at an asymptotically constant rate, and a horizon is formed in the exterior, which has the shape of a growing cigar. We are able to find the analytic form of the attractor solution describing the interior of the cosmic string. In particular, this enables us to analytically derive the relation between the string tension and the axial expansion rate. Furthermore, we show that the exterior conical singularity can be avoided for dynamical solutions. Our results might be relevant for theories with two extra dimensions, modeling our universe as a cosmic string with a three-dimensional axis. We derive the corresponding Friedmann equation, relating the on-brane Hubble parameter to the string tension or, equivalently, brane cosmological constant.


title: Apodized Pupil Lyot Coronagraphs for Arbitrary Apertures. IV. Reduced Inner Working Angle and Increased Robustness to Low-Order Aberrations

abstract: The Apodized Pupil Lyot Coronagraph (APLC) is a diffraction suppression system installed in the recently deployed instruments Palomar/P1640, Gemini/GPI, and VLT/SPHERE to allow direct imaging and spectroscopy of circumstellar environments. Using a prolate apodization, the current implementations offer raw contrasts down to $10^{-7}$ at 0.2 arcsec from a star over a wide bandpass (20\%), in the presence of central obstruction and struts, enabling the study of young or massive gaseous planets. Observations of older or lighter companions at smaller separations would require improvements in terms of inner working angle (IWA) and contrast, but the methods originally used for these designs were not able to fully explore the parameter space. We here propose a novel approach to improve the APLC performance. Our method relies on the linear properties of the coronagraphic electric field with the apodization at any wavelength to develop numerical solutions producing coronagraphic star images with high-contrast region in broadband light. We explore the parameter space by considering different aperture geometries, contrast levels, dark-zone sizes, bandpasses, and focal plane mask sizes. We present an application of these solutions to the case of Gemini/GPI with a design delivering a $10^{-8}$ raw contrast at 0.19 arcsec and offering a significantly reduced sensitivity to low-order aberrations compared to the current implementation. Optimal solutions have also been found to reach $10^{-10}$ contrast in broadband light regardless of the telescope aperture shape (in particular the central obstruction size), with effective IWA in the $2-3.5\lambda/D$ range, therefore making the APLC a suitable option for the future exoplanet direct imagers on the ground or in space.


title: Detection of a distinct metal-poor stellar halo in the early-type galaxy NGC 3115

abstract: We present the resolved stellar populations in the inner and outer halo of the nearby lenticular galaxy NGC~3115. Using deep HST observations, we analyze stars two magnitudes fainter than the tip of the red giant branch (TRGB). We study three fields along the minor axis of this galaxy, 19, 37 and 54 kpc from its center -- corresponding to 7, 14, 21 effective radii (r_{e}). Even at these large galactocentric distances, all of the fields are dominated by a relatively enriched population, with the main peak in the metallicity distribution decreasing with radius from [Z/H] ~ -0.5 to -0.65. The fraction of metal-poor stars ([Z/H] < -0.95) increases from 17%, at 16-37 kpc, to 28%, at ~54 kpc. We observe a distinct low metallicity population (peaked at [Z/H] ~ -1.3 and with total mass 2*10^{10}M_{\odot} ~ 14% of the galaxy's stellar mass) and argue that this represents the detection of an underlying low metallicity stellar halo. Such halos are generally predicted by galaxy formation theories and have been observed in several late type galaxies including the Milky Way and M31. The metallicity and spatial distribution of the stellar halo of NGC~3115 are consistent with the galaxy's globular cluster system, which has a similar low metallicity population that becomes dominant at these large radii. This finding supports the use of globular clusters as bright chemo-dynamical tracers of galaxy halos. These data also allow us to make a precise measurement of the magnitude of the TRGB, from which we derive a distance modulus of NGC~3115 of 30.05\pm0.05\pm0.10_{sys} (10.2\pm0.2\pm0.5_{sys} Mpc).


title: A detection threshold in the amplitude spectra calculated from Kepler data obtained during K2 mission

abstract: We present our analysis of simulated data in order to derive a detection threshold which can be used in the pre-whitening process of amplitude spectra. In case of ground-based data of pulsating stars, this threshold is conventionally taken to be four times the mean noise level in an amplitude spectrum. This threshold is questionable when space-based data are analyzed. Our effort is aimed at revising this threshold in the case of continuous 90-day Kepler K2 phase observations. Our result clearly shows that a 95% confidence level, common for ground observations, can be reached at 5.4 times the mean noise level and is coverage dependent. In addition, this threshold varies between 4.8 and 5.7, if the number of cadences is changed. This conclusion should secure further pre-whitening and helps to avoid over-interpretation of spectra of pulsating stars observed with the Kepler spacecraft during K2 phase. We compare our results with the standard approach widely used in the literature.


title: Astrochemical Correlations in Molecular Clouds

abstract: We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spatial resolution on the emission distribution and perform a robust quantitative comparison of the species to each other. We use the Spectral Correlation Function (SCF), which quantifies the root mean squared difference between spectra separated by some length scale, to characterize the structure of the simulated cloud in position-position-velocity (PPV) space. We predict the observed SCF for a broad range of observational tracers, and thus, identify homologous species. In particular, we show that the pairs C and CO, C$^{+}$ and CN, NH$_3$ and H$_2$CS have very similar SCFs. We measure the SCF slope variation as a function of beam size for all species and demonstrate that the beam size has a distinct effect on different species emission. However, for beams of up to 10", placing the cloud at 1 kpc, the change is not large enough to move the SCF slopes into different regions of parameter space. The results from this study provide observational guidance for choosing the best tracer to probe various cloud length scales.


title: The equilibrium view on dust and metals in galaxies: Galactic outflows drive low dust-to-metal ratios in dwarf galaxies

abstract: Most galaxy evolution simulations as well as a variety of observational methods assume a linear scaling between the (galaxy-averaged) dust-to-gas ratio D and metallicity Z of the interstellar medium (ISM). Indeed, nearby galaxies with solar or moderately sub-solar metallicities clearly follow this trend albeit with significant scatter. However, a growing number of observations show that the linear scaling breaks down for metal-poor galaxies (Z<0.2 Z_sun), highlighting the need for a more sophisticated modeling of the dust-to-metal ratio of galaxies. Here we study the co-evolution of dust and metal abundances in galaxies with the help of a dynamical, one-zone model that incorporates dust formation and destruction processes in addition to gas inflows, outflows, and metal enrichment. The dynamical model is consistent with various observational constraints, including the stellar mass -- metallicity relation, the stellar mass -- halo mass relation, and the observed Z -- D relation for both metal-poor and metal-rich galaxies. The functional form of the Z -- D relation follows from a basic equilibrium ansatz, similar to the ideas used previously to model the stellar mass -- metallicity relation. Galactic outflows regulate the inflow rate of gas from the cosmic web for galaxies of a given star formation rate. The mass loading factor of outflows thus dictates the rate at which the dust and metal content of the ISM is diluted. The stellar mass dependence of the mass loading factor drives the evolution of metallicities, dust-to-gas ratios, and dust-to-metal ratios in galaxies.


title: The density structure and star formation rate of non-isothermal polytropic turbulence

abstract: The interstellar medium of galaxies is governed by supersonic turbulence, which likely controls the star formation rate (SFR) and the initial mass function (IMF). Interstellar turbulence is non-universal, with a wide range of Mach numbers, magnetic fields strengths, and driving mechanisms. Although some of these parameters were explored, most previous works assumed that the gas is isothermal. However, we know that cold molecular clouds form out of the warm atomic medium, with the gas passing through chemical and thermodynamic phases that are not isothermal. Here we determine the role of temperature variations by modelling non-isothermal turbulence with a polytropic equation of state (EOS), where pressure and temperature are functions of gas density, P~rho^Gamma, T~rho^(Gamma-1). We use grid resolutions of 2048^3 cells and compare polytropic exponents Gamma=0.7 (soft EOS), Gamma=1 (isothermal EOS), and Gamma=5/3 (stiff EOS). We find a complex network of non-isothermal filaments with more small-scale fragmentation occurring for Gamma<1, while Gamma>1 smoothes out density contrasts. The density probability distribution function (PDF) is significantly affected by temperature variations, with a power-law tail developing at low densities for Gamma>1. In contrast, the PDF becomes closer to a lognormal distribution for Gamma<=1. We derive and test a new density variance - Mach number relation that takes Gamma into account. This new relation is relevant for theoretical models of the SFR and IMF, because it determines the dense gas mass fraction of a cloud, from which stars form. We derive the SFR as a function of Gamma and find that it decreases by a factor of ~5 from Gamma=0.7 to Gamma=5/3.

categories: astro-ph.GA astro-ph.SR physics.comp-ph physics.flu-dyn


title: CO mass upper limits in the Fomalhaut ring - the importance of NLTE excitation in debris discs and future prospects with ALMA

abstract: In recent years, gas has been observed in an increasing number of debris discs, though its nature remains to be determined. Here, we analyse CO molecular excitation in optically thin debris discs, and search ALMA Cycle-0 data for CO J=3-2 emission in the Fomalhaut ring. No significant line emission is observed; we set a 3-$\sigma$ upper limit on the integrated line flux of 0.16 Jy km s$^{-1}$. We show a significant dependency of the CO excitation on the density of collisional partners $n$, on the gas kinetic temperature $T_k$ and on the ambient radiation field $J$, suggesting that assumptions widely used for protoplanetary discs (e.g. LTE) do not necessarily apply to their low density debris counterparts. When applied to the Fomalhaut ring, we consider a primordial origin scenario where H$_2$ dominates collisional excitation of CO, and a secondary origin scenario dominated by e$^-$ and H$_2$O. In either scenario, we obtain a strict upper limit on the CO mass of 4.9 $\times$ 10$^{-4}$ M$_{\oplus}$. This arises in the non-LTE regime, where the excitation of the molecule is determined solely by the well-known radiation field. In the secondary scenario, assuming any CO present to be in steady state allows us to set an upper limit of $\sim$55% on the CO/H$_2$O ice ratio in the parent planetesimals. This could drop to $\sim$3% if LTE applies, covering the range observed in Solar System comets (0.4-30%). Finally, in light of our analysis, we present prospects for CO detection and characterisation in debris discs with ALMA.


title: An empirical calibration to estimate cool dwarf fundamental parameters from H-band spectra

abstract: Interferometric radius measurements provide a direct probe of the fundamental parameters of M dwarfs, but is within reach for only a limited sample of nearby, bright stars. We use interferometrically-measured radii, bolometric luminosities, and effective temperatures to develop new empirical calibrations based on low-resolution, near-infrared spectra. We use H-band Mg and Al features to derive calibrations for effective temperature, radius and log luminosity; the standard deviations in the residuals of our best fits are, respectively, 73K, 0.027Rsun, and 0.049 dex (11% error on luminosity). These relationships are valid for mid K to mid M dwarf stars, roughly corresponding to temperatures between 3100 and 4800K. We apply our calibrations to M dwarfs targeted by the MEarth transiting planet survey and to the cool Kepler Objects of Interest (KOIs). We independently validate our calibrations by demonstrating a clear relationship between our inferred parameters and the absolute K magnitudes of the MEarth stars, and we identify objects with magnitudes too bright for their estimated luminosities as candidate multiple systems. We also use our inferred luminosities to address the applicability of near-infrared metallicity calibrations to mid and late M dwarfs. The temperatures we infer for the KOIs agree remarkably well with those from the literature; however, our stellar radii are systematically larger than those presented in previous works that derive radii from model isochrones. This results in a mean planet radius that is 15% larger than one would infer using the stellar properties from recent catalogs. Our results confirm those of previous in-depth studies of Kepler-42, Kepler-45, and Kepler-186.


title: Permutation orbifolds and holography

abstract: Two dimensional conformal field theories with large central charge and a sparse low-lying spectrum are expected to admit a classical string holographic dual. We construct a large class of such theories employing permutation orbifold technology. In particular, we describe the group theoretic constraints on permutation groups to ensure a (stringy) holographic CFT. The primary result we uncover is that in order for the degeneracy of states to be finite in the large central charge limit, the groups of interest are the so-called oligomorphic permutation groups. Further requiring that the low-lying spectrum be sparse enough puts a bound on the number of orbits of these groups (on finite element subsets). Along the way we also study familiar cyclic and symmetric orbifolds to build intuition. We also demonstrate how holographic spectral properties are tied to the geometry of covering spaces for permutation orbifolds.


title: Variations of the ISM Compactness Across the Main Sequence of Star-Forming Galaxies: Observations and Simulations

abstract: (abridged) The majority of star-forming galaxies follow a simple empirical correlation in the star formation rate (SFR) versus stellar mass ($M_*$) plane, usually referred to as the star formation Main Sequence (MS). Here we combine a set of hydro-dynamical simulations of interacting galactic disks with state-of-the-art radiative transfer codes to analyze how the evolution of mergers is reflected upon the properties of the MS. We present \textsc{Chiburst}, a Markov Chain Monte Carlo (MCMC) Spectral Energy Distribution (SED) code that fits the multi-wavelength, broad-band photometry of galaxies and derives stellar masses, star formation rates, and geometrical properties of the dust distribution. We apply this tool to the SEDs of simulated mergers and compare the derived results with the reference output from the simulations. Our results indicate that changes in the SEDs of mergers as they approach coalescence and depart from the MS are related to an evolution of dust geometry in scales larger than a few hundred parsecs. This is reflected in a correlation between the specific star formation rate (sSFR), and the compactness parameter $\mathcal{C}$, that parametrizes this geometry and hence the evolution of dust temperature ($T_{\rm{dust}}$) with time. As mergers approach coalescence, they depart from the MS and increase their compactness, which implies that moderate outliers of the MS are consistent with late-type mergers. By further applying our method to real observations of Luminous Infrared Galaxies (LIRGs), we show that the merger scenario is unable to explain these extreme outliers of the MS. Only by significantly increasing the gas fraction in the simulations are we able to reproduce the SEDs of LIRGs.


title: Supermassive Black Hole Formation at High Redshifts via Direct Collapse in a Cosmological Context

abstract: We study the early stage of the formation of seed supermassive black holes via direct collapse in dark matter (DM) halos, in the cosmological context. We perform high-resolution zoom-in simulations of such collapse at high-$z$. Using the adaptive mesh refinement code ENZO, we resolve the formation and growth of a DM halo, until its virial temperature reaches $\sim 10^4$K, atomic cooling turns on, and collapse ensues. We demonstrate that direct collapse proceeds in two stages, although they are not well separated. The first stage is triggered by the onset of atomic cooling, and leads to rapidly increasing accretion rate with radius, from $\dot M\sim 0.1\,M_\odot {\rm yr^{-1}}$ at the halo virial radius to few $M_\odot \,{\rm yr^{-1}}$, around the scale radius $R_{\rm s}\sim 30$pc of the NFW DM density profile. The second stage of the collapse commences when the gas density takes precedence over the DM density. This is associated with the gas decoupling from the DM gravitational potential. The ensuing collapse approximates that of an isothermal sphere with $\dot M ( r )\sim $const. We confirm that the gas loses its angular momentum through non-axisymmetric perturbations and gravitational torques, to overcome the centrifugal barrier. During the course of the collapse, this angular momentum transfer process happens on nearly all spatial scales, and the angular momentum vector of the gas varies with position and time. Collapsing gas also exhibits supersonic turbulent motions which suppress gas fragmentation, and are characterized by density PDF consisting of a lognormal part and a high-density power law tail.


title: Direct observation of Higgs mode oscillations in the pump-probe photoemission spectra of electron-phonon mediated superconductors

abstract: Using the non-equilibrium Keldysh formalism, we solve the equations of motion for electron-phonon superconductivity, including an ultrafast pump field. We present results for time-dependent photoemission spectra out of equilibrium which probes the dynamics of the superconducting gap edge. The partial melting of the order by the pump field leads to oscillations at twice the melted gap frequency, a hallmark of the Higgs or amplitude mode. Thus the Higgs mode can be directly excited through the nonlinear effects of an electromagnetic field and detected without any additional symmetry breaking.


title: Asymptotic Symmetries of Massless QED in Even Dimensions

abstract: We consider the scattering of massless particles coupled to an abelian gauge field in 2n-dimensional Minkowski spacetime. Weinberg's soft photon theorem is recast as Ward identities for infinitely many new nontrivial symmetries of the massless QED S-matrix, with one such identity arising for each propagation direction of the soft photon. These symmetries are identified as large gauge transformations with angle-dependent gauge parameters that are constant along the null generators of null infinity. Almost all of the symmetries are spontaneously broken in the standard vacuum and the soft photons are the corresponding Goldstone bosons. Our result establishes a relationship between soft theorems and asymptotic symmetry groups in any even dimension.


title: Higher-order massive neutrino perturbations in large-scale structure

abstract: We develop a higher-order perturbation theory for large-scale structure formation involving a free-streaming hot or warm dark matter species. We focus on the case of mixed cold dark matter and massive neutrinos, although our approach is applicable also to a single warm dark matter species. In order to capture the suppressed growth of neutrino density perturbations on small scales, we account for the full momentum dependence of the phase space distribution using the Vlaslov equation, and derive from it a formal closed-form nonlinear equation for the neutrino density. Using a systematic perturbative expansion of this equation we compute high-order corrections to the neutrino density contrast without the explicit need to track the perturbed neutrino momentum distribution. We calculate the leading-order total matter bispectrum for several neutrino masses. Using our result as a benchmark, we test the accuracy of the fluid approximation and a linear approximation used in perturbative and N-body analyses, as well as a new hybrid approach that combines the exact linear evolution with the nonlinear structure of the fluid equations. Aiming at $\lesssim1%$ accuracy, we find that the total matter bispectrum with a low neutrino mass m=0.046 eV can be reproduced by all but the fluid approximation, while for larger neutrino masses m=0.46-0.93 eV only the hybrid approach has the desired accuracy on a large range of scales. This result serves as a cautionary note that approximate nonlinear models of neutrino clustering that reproduce the gross features of some observables may not suffice for precision calculations, nor are they guaranteed to apply to other observables. All of the approximation schemes fail to reproduce the bispectrum of the neutrino density perturbations at better than 20% accuracy across all scales, indicating that an exact treatment of nonlinear neutrino perturbations is necessary.


title: A natural SM-like 126 GeV Higgs via non-decoupling D-terms

abstract: Accommodating both a 126 GeV mass and Standard Model (SM) like couplings for the Higgs has a fine tuning price in supersymmetric models. Examples are the MSSM, in which SM-like couplings are natural, but raising the Higgs mass up to 126 GeV requires a considerable tuning, or the NMSSM, in which the situation is reversed: the Higgs is naturally heavier, but being SM-like requires some tuning. We show that models with non-decoupling D-terms alleviate this tension - a 126 GeV SM-like Higgs comes out basically with no fine tuning cost. In addition, the analysis of the fine tuning of the extended gauge sector shows that naturalness requires the heavy gauge bosons to likely be within the LHC run II reach.


title: Interacting surface states of three-dimensional topological insulators

abstract: We numerically investigate the surface states of a strong topological insulator in the presence of strong electron-electron interactions. We choose a spherical topological insulator geometry to make the surface amenable to a finite size analysis. The single-particle problem maps to that of Landau orbitals on the sphere with a magnetic monopole at the center that has unit strength and opposite sign for electrons with opposite spin. Assuming density-density contact interactions, we find superconducting and anomalous (quantum) Hall phases for attractive and repulsive interactions, respectively, as well as chiral fermion and chiral Majorana fermion boundary modes between different phases. Our setup is preeminently adapted to the search for topologically ordered surface terminations that could be microscopically stabilized by tailored surface interaction profiles.


title: Z$_2$ topology and superconductivity from symmetry lowering of a 3D Dirac Metal Au$_2$Pb

abstract: 3D Dirac semi-metals (DSMs) are materials that have massless Dirac electrons and exhibit exotic physical properties It has been suggested that structurally distorting a DSM can create a Topological Insulator (TI), but this has not yet been experimentally verified. Furthermore, quasiparticle excitations known as Majorana Fermions have been theoretically proposed to exist in materials that exhibit superconductivity and topological surface states. Here we show that the cubic Laves phase Au$_2$Pb has a bulk Dirac cone above 100 K that gaps out upon cooling at a structural phase transition to create a topologically non trivial phase that superconducts below 1.2 K. The nontrivial Z$_2$ = -1 invariant in the low temperature phase indicates that Au$_2$Pb in its superconducting state must have topological surface states. These characteristics make Au$_2$Pb a unique platform for studying the transition between bulk Dirac electrons and topological surface states as well as studying the interaction of superconductivity with topological surface states.


title: Strings, Branes and the Self-dual Solutions of Exceptional Field Theory

abstract: It has been shown that membranes and fivebranes are wave-like or monopole-like solutions in some higher dimensional theory. Here the picture is completed by combining the wave and monopole solutions into a single solution of Exceptional Field Theory. This solution solves the twisted self-duality constraint. The 1/2 BPS brane spectrum, consisting of fundamental, solitonic and Dirichlet branes, in ten- and eleven-dimensional supergravity may all be extracted from this single solution of Exceptional Field Theory. The solution's properties such as the singularity structure and its asymptotic behaviour near the core and at infinity are investigated.


title: Softened Gravity and the Extension of the Standard Model up to Infinite Energy

abstract: Attempts to solve naturalness by having the weak scale as the only breaking of classical scale invariance have to deal with two severe difficulties: gravity and the absence of Landau poles. We show that solutions to the first problem require premature modifications of gravity at scales no larger than $10^{11}$ GeV, while the second problem calls for many new particles at the weak scale. To build models that fulfil these properties, we classify 4-dimensional Quantum Field Theories that satisfy Total Asymptotic Freedom (TAF): the theory holds up to infinite energy, where all coupling constants flow to zero. We develop a technique to identify such theories and determine their low-energy predictions. Since the Standard Model turns out to be asymptotically free only under the unphysical conditions $g_1 = 0$, $M_t = 186$ GeV, $M_\tau = 0$, $M_h = 163$ GeV, we explore some of its weak-scale extensions that satisfy the requirements for TAF.


title: Detection of a large Be circumstellar disk during X-ray quiescence of XTE J1946+274

abstract: We present a multiwavelength study of the Be/X-ray binary system XTE J1946+274 with the main goal of better characterizing its behavior during X-ray quiescence. We aim to shed light on the mechanism which triggers the X-ray activity for this source. XTE J1946+274 was observed by Chandra-ACIS during quiescence in 2013 March 12. In addition, this source has been monitored from the ground-based astronomical observatories of El Teide (Tenerife, Spain), Roque de los Muchachos (La Palma, Spain) and Sierra Nevada (Granada, Spain) since 2011 September, and from the TUBITAK National Observatory (Antalya, Turkey) since 2005 April. We have performed spectral and photometric temporal analyses in order to investigate the quiescent state and transient behavior of this binary system. In 2006, a long mass ejection event took place from the Be star, lasting for about seven years. We also found that a large Be circumstellar disk was present during quiescence, although major X-ray activity was not observed. We made an attempt to explain this scenario by assuming the permanently presence of a tilted and warped Be circumstellar disk. The 0.3-10 keV X-ray spectrum during quiescence was well fitted with either an absorbed blackbody or absorbed power-law models. The main parameters obtained for these models were kT=1.43+/-0.17 and Gamma=0.9+/-0.4 (with N_H~2-7E+22 cm-2). The 0.3-10 keV flux of the source was ~0.8E-12 erg-1 cm-2 s-1. Pulsations were found with P_pulse=15.757(1) s (MJD 56363.115) and an rms pulse fraction of 32.1(3)%. The observed X-ray luminosity during quiescent periods may be explained by the neutron star being in supersonic propeller regimen.


title: Soft-virtual corrections to Higgs production at N$^3$LO

abstract: In this paper, we compute the soft-virtual corrections to Higgs boson production in gluon fusion for infinite top quark mass at next-to-next-to-next-to-leading order in QCD. In addition, we present analogous soft-virtual terms for both Drell-Yan lepton production in QCD and scalar pair production in N = 4 super Yang-Mills theory. The result for Drell-Yan lepton production is derived from the result for Higgs boson production using Casimir scaling arguments together with well-known results available in the literature. For scalar pair production in the N = 4 model, we show by explicit calculation that the result is equal to the part of the Higgs boson soft-virtual term which is of maximal transcendentality weight.


title: Thermal and Residual Excited-State Population in a 3D Transmon Qubit

abstract: We present a systematic study of the first excited-state population in a 3D transmon qubit mounted in a dilution refrigerator with a variable temperature. Using a modified version of the protocol developed by Geerlings et al. [1], we observe the excited-state population to be consistent with a Maxwell-Boltzmann distribution, i.e., a qubit in thermal equilibrium with the refrigerator, over the temperature range 35-150 mK. Below 35 mK, the excited-state population saturates to 0.1%, near the resolution of our measurement. We verified this result using a flux qubit with ten-times stronger coupling to its readout resonator. We conclude that these qubits have effective temperature T_{eff} = 35 mK. Assuming T_{eff} is due solely to hot quasiparticles, the inferred qubit lifetime is 108 us and in plausible agreement with the measured 80 us.


title: Cooperative Change Detection for Online Power Quality Monitoring

abstract: This paper considers the real-time power quality monitoring in power grid systems. The goal is to detect the occurrence of disturbances in the nominal sinusoidal voltage/current signal as quickly as possible such that protection measures can be taken in time. Based on an autoregressive (AR) model for the disturbance, we propose a generalized local likelihood ratio (GLLR) detector which processes meter readings sequentially and alarms as soon as the test statistic exceeds a prescribed threshold. The proposed detector not only reacts to a wide range of disturbances, but also achieves lower detection delay compared to the conventional block processing method. Then we further propose to deploy multiple meters to monitor the power signal cooperatively. The distributed meters communicate wirelessly to a central meter, where the data fusion and detection are performed. In light of the limited bandwidth of wireless channels, we develop a level-triggered sampling scheme, where each meter transmits only one-bit each time asynchronously. The proposed multi-meter scheme features substantially low communication overhead, while its performance is close to that of the ideal case where distributed meter readings are perfectly available at the central meter.


title: Light Higgsino Decays as a Probe of the NMSSM

abstract: In the Next-to Minimal Supersymmetric Standard Model (NMSSM), a sizable coupling $\lambda$ between the singlet and Higgs fields can naturally accommodate the observed Higgs boson mass of 125 GeV. This large coupling also results in a large separation between the Higgsino and singlino mass scales in the neutralino sector to evade stringent constraints from direct detection experiments. Most of the natural parameter space in this setup therefore contains light Higgsinos that can decay into the singlino and the 125 GeV Higgs. If the Higgsinos are to be light enough to be produced at the LHC, the mass gap $m_{\tilde\chi^0_2} -(m_{\tilde\chi^0_1} +m_h)$ is generally fairly low, resulting in small missing energy. We study the collider phenomenology of this process and demonstrate that the 4$b+j+{E\!\!\!\!/_{\rm T}}$ signal arising from $pp\rightarrow\tilde\chi_i^0\tilde\chi_j^0jj$ process can be a viable channel to search for this low mass gap region at the 14 TeV LHC. In addition, we find that a potential signal from the LHC search, in tandem with direct detection experiments, can distinguish the NMSSM from the analogous process in the MSSM, where the bino plays the role of the singlino, for positive values of the $\mu$ parameter even when the additional singlet-dominated (pseudo)scalars are absent or indistinguishable from the Higgs in the decay.


title: The Phase Diagram in Electron-Doped LCCO

abstract: Superconductors are a striking example of a quantum phenomenon in which electrons move coherently over macroscopic distances without scattering. The high-temperature superconducting oxides(cuprates) are the most studied class of superconductors, composed of two-dimensional CuO2 planes separated by other layers which control the electron concentration in the planes. A key unresolved issue in cuprates is the relationship between superconductivity and magnetism. In this paper, we report a sharp phase boundary of static three-dimensional magnetic order in the electron-doped superconductor La2-xCexCuO4-d where small changes in doping or depth from the surface switch the material from superconducting to magnetic. Using low-energy spin polarized muons, we find static magnetism disappears close to where superconductivity begins and well below the doping where dramatic changes in the transport properties are reported. These results indicate a higher degree of symmetry between the electron and hole-doped cuprates than previously thought.


title: An Unobscured type II quasar candidate: SDSS J012032.19-005501.9

abstract: We report the finding of an unobscured type II Active Galactic Nuclei (AGN) candidate, SDSS J012032.19-005501.9 at a relatively high redshift of 0.601,which shows a number of unusual properties. It varies significantly on timescales of years as typical type I AGNs and marginally on timescales of weeks. The color-magnitude relation and the structure function are also consistent with that of type I AGNs, which imply that its variability likely originates from the black hole accretion system .However, no broad emission line is detected in the SDSS spectrum, and the upper limit of the equivalent width of the H$\rm \beta$ broad emission line is much less than that of type I AGNs. These properties suggest that SDSS J012032.19-005501.9 may be an unobscured quasar without broad emission lines intrinsically, namely an unobscured type II AGN or "true" type II AGN. Furthermore, its continuum luminosity is at least one order of magnitude fainter than the average value of thepast century expected from the [OIII] emission line. It indicates that SDSS J012032.19-005501.9 may be switching off. Additional possible scenarios to explain this intriguing source are also discussed. Future deep observations at multi-wavelengths are needed to reveal the nature of this peculiar and intriguing AGN.


title: Evidence for interacting dark energy from BOSS

abstract: The result presented by the BOSS-SDSS Collaboration measuring the baryon acoustic oscillations of the Lyman-$\alpha$ forest from high-redshift quasars indicates a $2.5\sigma$ departure from the standard $\Lambda$-cold-dark-matter model. This is the first time that the evolution of dark energy at high redshifts has been measured, and the current results cannot be explained by simple generalizations of the cosmological constant. We show here that a simple phenomenological interaction in the dark sector provides a good explanation for this deviation, naturally accommodating the Hubble parameter obtained by BOSS, $H(z=2.34)=222 \pm 7 ~\mathrm{km~s^{-1}~Mpc^{-1}}$. By performing a global fit of the parameters with the inclusion of this new data set together with the Planck data for the interacting model, we are able to show that some interacting models have constraints for $H(2.34)$ and $D_\mathrm{A}(2.34)$ that are compatible with the ones obtained by the BOSS Collaboration, showing a better concordance than $\Lambda$CDM. We also show that the interacting models that have a small positive coupling constant, which helps alleviate the coincidence problem, are compatible with the cosmological observations. Adding the likelihood of these new baryon acoustic oscillations data shows an improvement in the global fit, although it is not statistically significant. The coupling constant could not be fully constrained by the data sets used, but the dark energy equation of state shows a slight preference for a value different from a cosmological constant.


title: Massive Gravity: A Lorentz-Symmetric Aether

abstract: This is a heuristic introduction to massive gravity based on an analogy with perfect fluids. I will argue that massive gravity can be thought of as Einstein gravity in the presence of a medium with unusual properties.


title: Persistent homology analysis of protein structure, flexibility and folding

abstract: Proteins are the most important biomolecules for living organisms. The understanding of protein structure, function, dynamics and transport is one of most challenging tasks in biological science. In the present work, persistent homology is, for the first time, introduced for extracting molecular topological fingerprints (MTFs) based on the persistence of molecular topological invariants. MTFs are utilized for protein characterization, identification and classification. The method of slicing is proposed to track the geometric origin of protein topological invariants. Both all-atom and coarse-grained representations of MTFs are constructed. A new cutoff-like filtration is proposed to shed light on the optimal cutoff distance in elastic network models. Based on the correlation between protein compactness, rigidity and connectivity, we propose an accumulated bar length generated from persistent topological invariants for the quantitative modeling of protein flexibility. To this end, a correlation matrix based filtration is developed. This approach gives rise to an accurate prediction of the optimal characteristic distance used in protein B-factor analysis. Finally, MTFs are employed to characterize protein topological evolution during protein folding and quantitatively predict the protein folding stability. An excellent consistence between our persistent homology prediction and molecular dynamics simulation is found. This work reveals the topology-function relationship of proteins.


title: Repeated, Delayed Torque Variations Following X-ray Flux Enhancements in the Magnetar 1E$\;$1048.1$-$5937

abstract: We report on two years of flux and spin evolution monitoring of 1E 1048.1$-$5937, a 6.5-s X-ray pulsar identified as a magnetar. Using {\it Swift} XRT data, we observed an X-ray outburst consisting of an increase in the persistent 1--10 keV flux by a factor of 6.3$\pm$0.2, beginning on 2011 December 31 (MJD 55926). Following a delay of $\sim100$ days, the magnetar entered a period of large torque variability, with $\dot{\nu}$ reaching a factor of $4.55\pm0.05$ times the nominal value, before decaying in an oscillatory manner over a time scale of months. We show by comparing to previous outbursts from the source that this pattern of behavior may repeat itself with a quasi-period of $\sim1800$ days. We compare this phenomenology to periodic torque variations in radio pulsars, finding some similarities which suggest a magnetospheric origin for the behavior of 1E 1048.1$-$5937.


title: Surface defects and instanton partition functions

abstract: We study the superconformal index of five-dimensional SCFTs and the sphere partition function of four-dimensional gauge theories with eight supercharges in the presence of co-dimension two half-BPS defects. We derive a prescription which is valid for defects which can be given a "vortex construction", i.e. can be defined by RG flow from vortex configurations in a larger theory. We test the prescription against known results and expected dualities. We employ our prescription to develop a general computational strategy for defects defined by coupling the bulk degrees of freedom to a Gauged Linear Sigma Model living in co-dimension two.


title: A streamlined difference ring theory: Indefinite nested sums, the alternating sign and the parameterized telescoping problem

abstract: We present an algebraic framework to represent indefinite nested sums over hypergeometric expressions in difference rings. In order to accomplish this task, parts of Karr's difference field theory have been extended to a ring theory in which also the alternating sign can be expressed. The underlying machinery relies on algorithms that compute all solutions of a given parameterized telescoping equation. As a consequence, we can solve the telescoping and creative telescoping problem in such difference rings.


title: Coexistence of density wave and superfluid order in a dipolar Fermi gas

abstract: We analyse the coexistence of superfluid and density wave (stripe) order in a quasi-two-dimensional gas of dipolar fermions aligned by an external field. Remarkably, the anisotropic nature of the dipolar interaction allows for such a coexistence in a large region of the zero temperature phase diagram. In this region, the repulsive part of the interaction drives the stripe formation and the attractive part induces the pairing, resulting in a supersolid with $p$-wave Cooper pairs aligned along the stripes. From a momentum space perspective, the stability of the supersolid phase is due to the fact that the stripe order renders the Fermi surface only partially gapped, leaving gapless regions that are most important for $p$-wave pairing. We finally demonstrate how this supersolid phase can be detected in time-of-flight experiments.


title: SUSY searches with the ATLAS detector

abstract: Highlights of the search for SUSY at ATLAS are presented emphasising on the variety of searches that were performed in the LHC Run I in p-p collisions at a centre of mass of 7 TeV and 8TeV. Searches for squars, gluinos, third generation squarks and electroweakinos in R-parity conserving scenarios and an example of searches in R-parity violating scenarios are highlighted.


title: Synthetic CO, H2 and HI surveys of the Galactic 2nd Quadrant, and the properties of molecular gas

abstract: We present CO, H2, HI and HISA distributions from a set of simulations of grand design spirals including stellar feedback, self-gravity, heating and cooling. We replicate the emission of the 2nd Galactic Quadrant by placing the observer inside the modelled galaxies and post process the simulations using a radiative transfer code, so as to create synthetic observations. We compare the synthetic datacubes to observations of the 2nd Quadrant of the Milky Way to test the ability of the current models to reproduce the basic chemistry of the Galactic ISM, as well as to test how sensitive such galaxy models are to different recipes of chemistry and/or feedback. We find that models which include feedback and self-gravity can reproduce the production of CO with respect to H2 as observed in our Galaxy, as well as the distribution of the material perpendicular to the Galactic plane. While changes in the chemistry/feedback recipes do not have a huge impact on the statistical properties of the chemistry in the simulated galaxies, we find that the inclusion of both feedback and self-gravity are crucial ingredients, as our test without feedback failed to reproduce all of the observables. Finally, even though the transition from H2 to CO seems to be robust, we find that all models seem to underproduce molecular gas, and have a lower molecular to atomic gas fraction than is observed. Nevertheless, our fiducial model with feedback and self-gravity has shown to be robust in reproducing the statistical properties of the basic molecular gas components of the ISM in our Galaxy.


title: Fast and Anisotropic Flexibility-Rigidity Index

abstract: The flexibility-rigidity index (FRI) is a newly proposed method for the construction of atomic rigidity functions. The FRI method analyzes protein rigidity and flexibility and is capable of predicting protein B-factors without resorting to matrix diagonalization. A fundamental assumption used in the FRI is that protein structures are uniquely determined by various internal and external interactions, while the protein functions, such as stability and flexibility, are solely determined by the structure. As such, one can predict protein flexibility without resorting to the protein interaction Hamiltonian. Consequently, bypassing the matrix diagonalization, the original FRI has a computational complexity of O(N^2). This work introduces a fast FRI (fFRI) algorithm for the flexibility analysis of large macromolecules. The proposed fFRI further reduces the computational complexity to O(N). Additionally, we propose anisotropic FRI (aFRI) algorithms for the analysis of protein collective dynamics. The aFRI algorithms admit adaptive Hessian matrices, from a completely global 3N*3N matrix to completely local 3*3 matrices. However, these local 3*3 matrices have built in much non-local correlation. Furthermore, we compare the accuracy and efficiency of FRI with some {established} approaches to flexibility analysis, namely, normal mode analysis (NMA) and Gaussian network model (GNM). The accuracy of the FRI method is tested. The FRI, particularly the fFRI, is orders of magnitude more efficient and about 10% more accurate overall than some of the most popular methods in the field. The proposed fFRI is able to predict B-factors for alpha-carbons of the HIV virus capsid (313,236 residues) in less than 30 seconds on a single processor using only one core. Finally, we demonstrate the application of FRI and aFRI to protein domain analysis.


title: A Robust and Scalable Algorithm for the Steiner Problem in Graphs

abstract: We present an effective heuristic for the Steiner Problem in Graphs. Its main elements are a multistart algorithm coupled with aggressive combination of elite solutions, both leveraging recently-proposed fast local searches. We also propose a fast implementation of a well-known dual ascent algorithm that not only makes our heuristics more robust (by quickly dealing with easier cases), but can also be used as a building block of an exact (branch-and-bound) algorithm that is quite effective for some inputs. On all graph classes we consider, our heuristic is competitive with (and sometimes more effective than) any previous approach with similar running times. It is also scalable: with long runs, we could improve or match the best published results for most open instances in the literature.


title: Fitting and projecting HIV epidemics: Data, structure and parsimony

abstract: Understanding historical trends in the epidemic of HIV is important for assessing current and projecting future trends in prevalence, incidence and mortality and for evaluating the impact and cost-effectiveness of control measures. In generalized epidemics the available data are of variable quality among countries and limited mainly to trends in the prevalence of HIV among women attending ante-natal clinics. In concentrated epidemics one needs, at the very least, time trends in the prevalence of HIV among different risk groups, including intravenous drug users, men-who-have-sex-with-men, and commercial sex workers as well as the size of each group and the degree of overlap between them. Here we focus on the comparatively straight forward problems presented by generalized epidemics. We fit data from Kenya to a susceptible-infected model and then successively add structure to the model, drawing on our knowledge of the natural history of HIV, to explore the effect that different structural aspects of the model have on the fits and the projections. Both heterogeneity in risk and changes in behaviour over time are important but easily confounded. Using a Weibull rather than exponential survival function for people infected with HIV, in the absence of treatment, makes a significant difference to the estimated trends in incidence and mortality and to the projected trends. Allowing for population growth has a small effect on the fits and the projections but is easy to include. Including details of the demography adds substantially to the complexity of the model, increases the run time by several orders of magnitude, but changes the fits and projections only slightly and to an extent that is less than the uncertainty inherent in the data. We make specific recommendations for the kind of model that would be suitable for understanding and managing HIV epidemics in east and southern Africa.


title: Instanton Operators in Five-Dimensional Gauge Theories

abstract: We discuss instanton operators in five-dimensional gauge theories. These are defined as disorder operators which create a non-vanishing second Chern class on a four-sphere surrounding their insertion point. As such they may be thought of as higher-dimensional analogues of three-dimensional monopole (or `t Hooft) operators. We argue that they play an important role in the enhancement of the Lorentz symmetry for maximally supersymmetric Yang-Mills to SO(1,5) at strong coupling.


title: Inflation and Dark Energy with a Single Superfield

abstract: We discuss the possibility to construct supergravity models with a single superfield describing inflation as well as the tiny cosmological constant $V \sim 10^{{-120}}$. One could expect that the simplest way to do it is to study models with a supersymmetric Minkowski vacuum and then slightly uplift them. However, due to the recently proven no-go theorem, such a tiny uplifting cannot be achieved by a small modification of the parameters of the theory. We illustrate this general result by investigation of models with a single chiral superfield recently proposed by Ketov and Terada. We show that the addition of a small constant or a linear term to the superpotential of a model with a stable supersymmetric Minkowski vacuum converts it to an AdS vacuum, which results in a rapid cosmological collapse. One can avoid this problem and uplift a supersymmetric Minkowski vacuum to a dS vacuum with $V_{0}\sim 10^{-120}$ without violating the no-go theorem by making these extra terms large enough. However, we show that this leads to a strong supersymmetry breaking in the uplifted vacua.


title: Effect of Electron-Phonon Coupling on Thermal Transport across Metal-Nonmetal Interface - A Second Look

abstract: The effect of electron-phonon (e-ph) coupling on thermal transport across metal-nonmetal interfaces is yet to be completely understood. In this paper, we use a series of molecular dynamics (MD) simulations with e-ph coupling effect included by Langevin dynamics to calculate the thermal conductance at a model metal-nonmetal interface. It is found that while e-ph coupling can present additional thermal resistance on top of the phonon-phonon thermal resistance, it can also make the phonon-phonon thermal conductance larger than the pure phonon transport case. This is because the e-ph interaction can disturb the phonon subsystem and enhance the energy communication between different phonon modes inside the metal. This facilitates redistributing phonon energy into modes that can more easily transfer energy across the interfaces. Compared to the pure phonon thermal conduction, the total thermal conductance with e-ph coupling effect can become either smaller or larger depending on the coupling factor. This result helps clarify the role of e-ph coupling in thermal transport across metal-nonmetal interface.


title: Spin Transitions in Graphene Butterflies at an Integer Filling Factor

abstract: Recent experiments on the role of electron-electron interactions in fractal Dirac systems have revealed a host of interesting effects, in particular, the unique nature of the magnetic field dependence of butterfly gaps in graphene. The novel gap structure observed in the integer quantum Hall effect is quite intriguing [Nat. Phys. 10, 525 (2014)], where one observes a suppression of the ferromagnetic state at one value of the commensurable flux but a reentrant ferromagnetic state at another. Our present work that includes the interplay between the electron-electron interaction and the periodic potential, explains the underlying physical processes that can lead to such a unique behavior of the butterfly gaps in that system where spin flip transitions are involved in the ground state.


title: The Chandra Planetary Nebulae Survey (ChanPlaNS): III. X-ray Emission from the Central Stars of Planetary Nebulae

abstract: We present X-ray spectral analysis of 20 point-like X-ray sources detected in Chandra Planetary Nebula Survey (ChanPlaNS) observations of 59 planetary nebulae (PNe) in the solar neighborhood. Most of these 20 detections are associated with luminous central stars within relatively young, compact nebulae. The vast majority of these point-like X-ray-emitting sources at PN cores display relatively "hard" ($\geq0.5$~keV) X-ray emission components that are unlikely to be due to photospheric emission from the hot central stars (CSPN). Instead, we demonstrate that these sources are well modeled by optically-thin thermal plasmas. From the plasma properties, we identify two classes of CSPN X-ray emission: (1) high-temperature plasmas with X-ray luminosities, $L_{\rm X}$, that appear uncorrelated with the CSPN bolometric luminosity, $L_{\rm bol}$; and (2) lower-temperature plasmas with $L_{\rm X}/L_{\rm bol}\sim10^{-7}$. We suggest these two classes correspond to the physical processes of magnetically active binary companions and self-shocking stellar winds, respectively. In many cases this conclusion is supported by corroborative multiwavelength evidence for the wind and binary properties of the PN central stars. By thus honing in on the origins of X-ray emission from PN central stars, we enhance the ability of CSPN X-ray sources to constrain models of PN shaping that invoke wind interactions and binarity.


title: Multi-strange baryon production in pp, p-Pb and Pb-Pb collisions at LHC measured with ALICE

abstract: Transverse momentum spectra and yields of charged $\Xi$ and $\Omega$ at mid-rapidity in pp, p-Pb and Pb-Pb collisions at the LHC have been measured by the ALICE Collaboration. These baryons are identified by reconstruction of their weak decay topology, in modes with only charged decay products, using the excellent tracking and particle identification capabilities of the detector. The recent measurements of the multi-strange baryon production relative to non-strange particles in p-Pb collisions are presented: this would help to understand the change in relative strangeness production from pp collisions to Pb-Pb collisions. Results on the nuclear modification factors for the charged $\Xi$ and $\Omega$ particles, compared with those for other light particles, are also reported.


title: Linear Controller Design for Chance Constrained Systems

abstract: This paper is concerned with the design of a linear control law for linear systems with stationary additive disturbances. The objective is to find a state feedback gain that minimizes a quadratic stage cost function, while observing chance constraints on the input and/or the state. Unlike most of the previous literature, the chance constraints (and the stage cost) are not considered on each input/state of the transient response. Instead, they refer to the input/state of the closed-loop system in its stationary mode of operation. Hence the control is optimized for a long-run, rather than a finite-horizon operation. The controller synthesis can be cast as a convex semi-definite program (SDP). The chance constraints appear as linear matrix inequalities. Both single chance constraints (SCCs) and joint chance constraints (JCCs) on the input and/or the state can be included. If the disturbance is Gaussian, additionally to WSS, this information can be used to improve the controller design. The presented approach can also be extended to the case of output feedback. The entire design procedure is flexible and easy to implement, as demonstrated on a short illustrative example.


title: On a random search tree: asymptotic enumeration of vertices by distance from leaves

abstract: A random binary search tree grown from the uniformly random permutation of $[n]$ is studied. We analyze the exact and asymptotic counts of vertices by rank, the distance from the set of leaves. The asymptotic fraction $c_k$ of vertices of a fixed rank $k\ge 0$ is shown to decay exponentially with $k$. Notoriously hard to compute, the exact fractions $c_k$ had been determined for $k\le 3$ only. We computed $c_4$ and $c_5$ as well; both are ratios of enormous integers, denominator of $c_5$ being $274$ digits long. Prompted by the data, we proved that, in sharp contrast, the largest prime divisor of $c_k$'s denominator is $2^{k+1}+1$ at most. We conjecture that, in fact, the prime divisors of every denominator for $k>1$ form a single interval, from $2$ to the largest prime not exceeding $2^{k+1}+1$.


title: Correlation between crystalline order and vitrification in colloidal monolayers

abstract: We investigate experimentally the relationship between local structure and dynamical arrest in a quasi-2d colloidal model system which approximates hard discs. We introduce polydispersity to the system to suppress crystallisation. Upon compression, the increase in structural relaxation time is accompanied by the emergence of local hexagonal symmetry. Examining the dynamical heterogeneity of the system, we identify three types of motion : "zero-dimensional" corresponding to beta-relaxation, "one-dimensional" or stringlike motion and "two-dimensional" motion. The dynamic heterogeneity is correlated with the local order, that is to say locally hexagonal regions are more likely to be dynamically slow. However we find that lengthscales corresponding to dynamic heterogeneity and local structure do not appear to scale together approaching the glass transition.


title: Estimation of a non-stationary model for annual precipitation in southern Norway using replicates of the spatial field

abstract: Estimation of stationary dependence structure parameters using only a single realisation of the spatial process, typically leads to inaccurate estimates and poorly identified parameters. A common way to handle this is to fix some of the parameters, or within the Bayesian framework, impose prior knowledge. In many applied settings, stationary models are not flexible enough to model the process of interest, thus non-stationary spatial models are used. However, more flexible models usually means more parameters, and the identifiability problem becomes even more challenging. We investigate aspects of estimation of a Bayesian non-stationary spatial model for annual precipitation using observations from multiple years. The model contains replicates of the spatial field, which increases precision of the estimates and makes them less prior sensitive. Using R-INLA, we analyse precipitation data from southern Norway, and investigate statistical properties of the replicate model in a simulation study. The non-stationary spatial model we explore belongs to a recently introduced class of stochastic partial differential equation (SPDE) based spatial models. This model class allows for non-stationary models with explanatory variables in the dependence structure. We derive conditions to facilitate prior specification for these types of non-stationary spatial models.


title: Impact of User Pairing on 5G Non-Orthogonal Multiple Access

abstract: Non-orthogonal multiple access (NOMA) represents a paradigm shift from conventional orthogonal multiple access (MA) concepts, and has been recognized as one of the key enabling technologies for 5G systems. In this paper, the impact of user pairing on the performance of two NOMA systems, NOMA with fixed power allocation (F-NOMA) and cognitive radio inspired NOMA (CR-NOMA), is characterized. For FNOMA, both analytical and numerical results are provided to demonstrate that F-NOMA can offer a larger sum rate than orthogonal MA, and the performance gain of F-NOMA over conventional MA can be further enlarged by selecting users whose channel conditions are more distinctive. For CR-NOMA, the quality of service (QoS) for users with the poorer channel condition can be guaranteed since the transmit power allocated to other users is constrained following the concept of cognitive radio networks. Because of this constraint, CR-NOMA has different behavior compared to F-NOMA. For example, for the user with the best channel condition, CR-NOMA prefers to pair it with the user with the second best channel condition, whereas the user with the worst channel condition is preferred by F-NOMA.


title: A new non-Hermitian E2-quasi-exactly solvable model

abstract: We construct a previously unknown $E_2$-quasi-exactly solvable non-Hermitian model whose eigenfunctions involve weakly orthogonal polynomials obeying three-term recurrence relations that factorize beyond the quantization level. The model becomes Hermitian when one of its two parameters is fixed to a specific value. We analyze the double scaling limit of this model leading to the complex Mathieu equation. The norms, Stieltjes measures and moment functionals are evaluated for some concrete values of one of the two parameters.


title: Consimilarity and quaternion matrix equations $AX-\hat{X}B=C$, $X-A\hat{X}B=C$

abstract: L.Huang [Linear Algebra Appl. 331 (2001) 21-30] gave a canonical form of a quaternion matrix $A$ with respect to consimilarity transformations $\tilde{S}^{-1}AS$ in which $S$ is a nonsingular quaternion matrix and $\tilde{h}:=a-bi+cj-dk$ for each quaternion $h=a+bi+cj+dk$. We give an analogous canonical form of a quaternion matrix with respect to consimilarity transformations $\hat{S}^{-1}AS$ in which $h\mapsto\hat{h}$ is an arbitrary involutive automorphism of the skew field of quaternions. We apply the obtained canonical form to the quaternion matrix equations $AX-\hat{X}B=C$ and $X-A\hat{X}B=C$.


title: Structure of Lefschetz thimbles in simple fermionic systems

abstract: The Picard-Lefschetz theory offers a promising tool to solve the sign problem in QCD and other field theories with complex path-integral weight. In this paper the Lefschetz-thimble approach is examined in simple fermionic models which share some features with QCD. In zero-dimensional versions of the Gross-Neveu model and the Nambu-Jona-Lasinio model, we study the structure of Lefschetz thimbles and its variation across the chiral phase transition. We map out a phase diagram in the complex four-fermion coupling plane using a thimble decomposition of the path integral, and demonstrate an interesting link between anti-Stokes lines and Lee-Yang zeros. In the case of nonzero mass, it is shown that the approach to the chiral limit is singular because of intricate cancellation between competing thimbles, which implies the necessity to sum up multiple thimbles related by symmetry. We also consider a Chern-Simons theory with fermions in $0+1$-dimension and show how Lefschetz thimbles solve the complex phase problem caused by a topological term. These prototypical examples would aid future application of this framework to bona fide QCD.


title: KAM for the nonlinear beam equation 1: small-amplitude solutions

abstract: In this paper we prove a KAM result for the non linear beam equation on the d-dimensional torus $$u_{tt}+\Delta^2 u+m u + g(x,u)=0\ ,\quad t\in { \mathbb{R}} , \; x\in {\mathbb T}^d, \qquad \qquad (*) $$ where $g(x,u)=4u^3+ O(u^4)$. Namely, we show that, for generic $m$, most of the small amplitude invariant finite dimensional tori of the linear equation $(*)_{g=0}$, written as the system $$ u_t=-v,\quad v_t=\Delta^2 u+mu, $$, persist as invariant tori of the nonlinear equation $(*)$, re-written similarly. If $d\ge2$, then not all the persisted tori are linearly stable, and we construct explicit examples of partially hyperbolic invariant tori. The unstable invariant tori, situated in the vicinity of the origin, create around them some local instabilities, in agreement with the popular belief in nonlinear physics that small-amplitude solutions of space-multidimensonal hamiltonian PDEs behave in a chaotic way. The proof uses an abstract KAM theorem from another our publication.


title: On the elliptic genera of manifolds of Spin(7) holonomy

abstract: Superstring compactification on a manifold of Spin(7) holonomy gives rise to a 2d worldsheet conformal field theory with an extended supersymmetry algebra. The $\mathcal{N}=1$ superconformal algebra is extended by additional generators of spins 2 and 5/2, and instead of just superconformal symmetry one has a $c=12$ realization of the symmetry group $\mathcal{SW}(3/2,2)$. In this paper, we compute the characters of this supergroup, and decompose the elliptic genus of a general Spin(7) compactification in terms of these characters. We find suggestive relations to various sporadic groups, which are made more precise in a companion paper.


title: Surface properties of liquid mercury: a comparison of density-dependent and density-independent force fields

abstract: Motivated by an experimental interest we investigate by the means of atomistic Molecular Dynamics simulation the ability of density-independent, empiric density-dependent, and recently proposed embedded-atom force fields for liquid mercury to predict the surface tension of the free surface of liquid mercury at the temperature of 293~K. The effect of the density dependence of the studied models on the liquid-vapor coexistence and surface tension is discussed in detail. In view of computational efficiency of the density-independent model we optimize its functional form to obtain higher surface tension values in order to improve agreement with experiment. The results are also corroborated by Monte Carlo simulations and semi-analytic estimations of the liquid-vapor coexistence density.


title: Non-equilibrium scalar two point functions in AdS/CFT

abstract: In the first part of the paper, we discuss different versions of the AdS/CFT dictionary out of equilibrium. We show that the Skenderis - van Rees prescription and the "extrapolate" dictionary are equivalent at the level of "in-in" two point functions of free scalar fields in arbitrary asymptotically AdS spacetimes. In the second part of the paper, we calculate two point correlation functions in dynamical spacetimes using the "extrapolate" dictionary. These calculations are performed for conformally coupled scalar fields in examples of spacetimes undergoing gravitational collapse, the AdS$_2$-Vaidya spacetime and the AdS$_3$-Vaidya spacetime, which allow us to address the problem of thermalization following a quench in the boundary field theory. The computation of the correlators is formulated as an initial value problem in the bulk spacetime. Finally, we compare our results for AdS$_3$-Vaidya to results in the previous literature obtained using the geodesic approximation and we find qualitative agreement.


title: Conjugate Directions in Lattice Landau and Coulomb Gauge Fixing

abstract: We provide details expanding on our implementation of a non-linear conjugate gradient method with Fourier acceleration for lattice Landau and Coulomb gauge fixing. We find clear improvement over the Fourier accelerated steepest descent method, with the average time taken for the algorithm to converge to a fixed, high accuracy, being reduced by a factor of 2 to 4. We show such improvement for the logarithmic definition of the gauge fields here, having already shown this to be the case for a more common definition. We also discuss the implementation of an optimal Fourier accelerated steepest descent method.


title: The extension of distributions on manifolds, a microlocal approach

abstract: Let $M$ be a smooth manifold, $I\subset M$ a closed embedded submanifold of $M$ and $U$ an open subset of $M$. In this paper, we find conditions using a geometric notion of scaling for $t\in \mathcal{D}^\prime(U\setminus I)$ to admit an extension in $\mathcal{D}^\prime(U)$. We give microlocal conditions on $t$ which allow to control the wave front set of the extension generalizing a previous result of Brunetti--Fredenhagen. Furthermore, we show that there is a subspace of extendible distributions for which the wave front of the extension is minimal which has applications for the renormalization of quantum field theory on curved spacetimes.


title: On Distribution of Product of Stable Laws

abstract: We derive the probability distribution of product of two independent random variables, each distributed according the one-dimensional stable law. We represent the density by its power series and its asymptotic expansions. As Fox's H-functions with a particular choice of parameters well describe the densities of stable laws, we discuss the choice of parameters for Fox's H-function such that it matches the derived densities. As a consequence, we give representations of these particular Fox's H-functions in terms of its power series.


title: Vertical dipole above a dielectric or metallic half-space - energy flow considerations

abstract: The emission pattern from a classical dipole located above and oriented perpendicular to a metallic or dielectric half space is calculated for a dipole driven at constant amplitude. This is a problem considered originally by Sommerfeld and analyzed subsequently by numerous authors. In contrast to most previous treatments, however, we focus on the energy flow in the metal or dielectric. It is shown that the radial Poynting vector in the metal points inwards when the frequency of the dipole is below the surface plasmon resonance frequency. In this case, energy actually flows of the interface at small radii. The Joule heating in the metal is also calculated and it is shown explicitly that Poynting's theorem holds for a cylindrical surface in the metal. When the metal is replaced by a dielectric having permittivity less than that of the medium in which the dipole is immersed, it is found that energy flows out of the interface for sufficiently large radii. In all cases it is assumed that the imaginary part of the permittivity of the metal or dielectric is much less than unity.


title: Force-Extension for DNA in a Nanoslit: Using an Effective Dimensionality to Map between the 3D and 2D Limits

abstract: The force-extension relation for a semi-flexible polymer such as DNA confined in a nanoslit is investigated and it is found that both the effective persistence length and the form of the force-extension relation change as the chain goes from 3D (very large slit heights) to 2D (very tight confinement). Generalizations of the Marko-Siggia relation appropriate for polymers in nanoconfinement are presented. The forms for both strong and weak confinement regimes are characterized by an \textit{effective dimensionality}. At low forces, the effective dimensionality is given by the correlations along the polymer in the plane of the confining walls. At high forces, the theoretical force must account for reduced conformation space. Together the interpolations give good agreement for all slit heights at all forces. As DNA and other semi-flexible biopolymers are commonly confined \textit{in situ} to various degrees, both the idea of an effective dimensionality and the associated generalized Marko-Siggia interpolations are useful for qualitatively understanding and quantitatively modeling polymers in nanoconfinement.


title: Unsupervised Induction of Semantic Roles within a Reconstruction-Error Minimization Framework

abstract: We introduce a new approach to unsupervised estimation of feature-rich semantic role labeling models. Our model consists of two components: (1) an encoding component: a semantic role labeling model which predicts roles given a rich set of syntactic and lexical features; (2) a reconstruction component: a tensor factorization model which relies on roles to predict argument fillers. When the components are estimated jointly to minimize errors in argument reconstruction, the induced roles largely correspond to roles defined in annotated resources. Our method performs on par with most accurate role induction methods on English and German, even though, unlike these previous approaches, we do not incorporate any prior linguistic knowledge about the languages.


title: Joint Segmentation and Deconvolution of Ultrasound Images Using a Hierarchical Bayesian Model based on Generalized Gaussian Priors

abstract: This paper proposes a joint segmentation and deconvolution Bayesian method for medical ultrasound (US) images. Contrary to piecewise homogeneous images, US images exhibit heavy characteristic speckle patterns correlated with the tissue structures. The generalized Gaussian distribution (GGD) has been shown to be one of the most relevant distributions for characterizing the speckle in US images. Thus, we propose a GGD-Potts model defined by a label map coupling US image segmentation and deconvolution. The Bayesian estimators of the unknown model parameters, including the US image, the label map and all the hyperparameters are difficult to be expressed in closed form. Thus, we investigate a Gibbs sampler to generate samples distributed according to the posterior of interest. These generated samples are finally used to compute the Bayesian estimators of the unknown parameters. The performance of the proposed Bayesian model is compared with existing approaches via several experiments conducted on realistic synthetic data and in vivo US images.


title: Color Hom-Akivis algebras, Color Hom-Leibniz algebras and Modules over Color Hom-Leibniz algebras

abstract: In this paper we introduce color Hom-Akivis algebras and prove that the commutator of any color non-associative Hom-algebra structure map leads to a color Hom-akivis algebra. We give various constructions of color Hom-Akivis algebras. Next we study flexible and alternative color Hom-Akivis algebras. Likewise color Hom-Akivis algebras, we introduce non-commutative color Hom-Leibniz-Poisson algebras and presente several constructions. Moreover we give the relationship between Hom-dialgebras and Hom-Leibniz-Poisson algebras; i.e. a Hom-dialgebras give rise to a Hom-Leibniz-Poisson algebra. Finally we show that twisting a color Hom-Leibniz module structure map by a color Hom-Leibniz algebra endomorphism, we get another one.


title: Numerical modeling of the heterocycle intercalated proton-conducting polymers at various mole ratios

abstract: The kinetic Monte Carlo simulations are employed to study the proton conductivity for anhydrous heterocyclic based polymers. The proton transport is based on a two-step process called the Grotthuss mechanism. In the referring system the proton concentration depends on the relative molar ratio, $x$, of the benzimidazole and the polystyrene sulfonic acid. Available experimental data with contrasting behavior are fitted and interpreted in terms of our microscopic model. Moreover, it has been shown that the current behavior similar to the Vogel-Tamman-Fulcher law can be reproduced with high precision on the basis of the Grotthuss mechanism.


title: The proton conductivity in benzimidazolium azelate under moderate pressure

abstract: The kinetic Monte Carlo method is applied to examine effects of hydrostatic pressure on the benzimidazolium azelate (BenAze) proton conductivity. Following the experimental indications the recently proposed model has been modified to simulate the transport phenomena under moderate pressure, resulting in a very good agreement between numerical and experimental results. We demonstrate that the pressure-induced changes in the proton conductivity can be attributed to solely two parameters: the length of the hydrogen bond and the amplitude of lattice vibrations while other processes play a more minor role. It may provide an insight into tailoring new materials with improved proton-conducting properties. Furthermore, in high-pressure regime we anticipate the crossover from the increasing to decreasing temperature dependence of the proton conductivity resulting from the changes in the hydrogen-bond activation barrier with increased pressure.


title: Diffusion Estimation Over Cooperative Multi-Agent Networks With Missing Data

abstract: In many fields, and especially in the medical and social sciences and in recommender systems, data are gathered through clinical studies or targeted surveys. Participants are generally reluctant to respond to all questions in a survey or they may lack information to respond adequately to some questions. The data collected from these studies tend to lead to linear regression models where the regression vectors are only known partially: some of their entries are either missing completely or replaced randomly by noisy values. In this work, assuming missing positions are replaced by noisy values, we examine how a connected network of agents, with each one of them subjected to a stream of data with incomplete regression information, can cooperate with each other through local interactions to estimate the underlying model parameters in the presence of missing data. We explain how to adjust the distributed diffusion through (de)regularization in order to eliminate the bias introduced by the incomplete model. We also propose a technique to recursively estimate the (de)regularization parameter and examine the performance of the resulting strategy. We illustrate the results by considering two applications: one dealing with a mental health survey and the other dealing with a household consumption survey.


title: The hippocampal-striatal circuit for goal-directed and habitual choice

abstract: It is now widely accepted that one of the roles of the hippocampus is to maintain episodic spatial representations, while parallel striatal pathways contribute to both declarative and procedural value computations by encoding different input-specific outcome predictions. In this paper we investigate the use of these brain mechanisms for action selection, linking them to model-based and model-free controllers for decision making. To this aim we propose a biologically inspired computational model that embodies these theories and explains the functioning of the hippocampal-striatal circuit in a rat navigation task. Its main characteristic is to allow the cooperation of habitual and goal-directed behaviors, with the hippocampus primarily involved in encoding spatial information and simulating possible navigation paths, and the ventral and dorsal striatum involved in learning stimulus-response behaviors and evaluating the reward expectancies associated to predicted locations and sensed stimuli, respectively. The architecture we present employs an unsupervised reinforcement learning rule for the hippocampal-striatal network that is able to build a representation of the environment in which rewarding sites and informative landmarks produce value gradients that are used for planning and decision making. Additionally, it utilizes an arbitration mechanism that balances between exploitation, i.e. stimulus-response behaviors, and mental exploration, i.e. motor imagery processes, based on the intensity and the variability of the responses of striatal neurons. We interpret these results in light of recent experimental data that show anticipatory activations in hippocampal and striatal areas.


title: Dissipation, Energy Transfer, and Spindown Luminosity in 2.5D PIC Simulations of the Pulsar Magnetosphere

abstract: We perform 2.5D axisymmetric simulations of the pulsar magnetosphere (aligned dipole rotator) using the charge conservative, relativistic, electromagnetic particle in cell code PICsar. Particle in cell codes are a powerful tool to use for studying the pulsar magnetosphere, because they can handle the force-free and vacuum limits and provide a self-consistent treatment of magnetic reconnection. In the limit of dense plasma throughout the magnetosphere, our solutions are everywhere in the force-free regime except for dissipative regions at the polar caps, in the current layers, and at the Y-point. These dissipative regions arise self-consistently, since we do not have any explicit dissipation in the code. A minimum of $ \approx 15-20\%$ of the electromagnetic spindown luminosity is transferred to the particles inside 5 light cylinder radii. However, the particles can carry as much as $\gtrsim 50 \%$ of the spindown luminosity if there is insufficient plasma in the outer magnetosphere to screen the component of electric field parallel to the magnetic field. In reality, the component of the spindown luminosity carried by the particles could be radiated as gamma rays, but high-frequency synchrotron emission would need to be implemented as a sub-grid process in our simulations and is not present for the current suite of runs. The value of the spindown luminosity in our simulations is within $\approx 10\%$ of the force-free value, and the structure of the electromagnetic fields in the magnetosphere is on the whole consistent with the force-free model.


title: A Compact High Velocity Cloud near the Magellanic Stream: Metallicity and Small-Scale Structure

abstract: The Magellanic Stream (MS) is a well-resolved gaseous tail originating from the Magellanic Clouds. Studies of its physical properties and chemical composition are needed to understand its role in Galactic evolution. We investigate the properties of a compact HVC (CHVC 224.0-83.4-197) lying close on the sky to the MS to determine whether it is physically connected to the Stream and to examine its internal structure. Our study is based on analysis of HST/COS spectra of three QSOs (Ton S210, B0120-28, and B0117-2837) all of which pass through this single cloud at small angular separation (\lessim 0.72{\deg}), allowing us to compare physical conditions on small spatial scales. No significant variation is detected in the ionization structure from one part of the cloud to the other. Using Cloudy photoionization models, toward Ton S210 we derive elemental abundances of [C/H] = -1.21 +/- 0.11, [Si/H] = -1.16 +/- 0.11, [Al/H] = -1.19 +/- 0.17 and [O/H] = -1.12 +/- 0.22, which agree within 0.09 dex. The CHVC abundances match the 0.1 solar abundances measured along the main body of the Stream. This suggests that the CHVC (and by extension the extended network of filaments to which it belongs) has an origin in the MS. It may represent a fragment that has been removed from the Stream as it interacts with the gaseous Galactic halo.


title: Zipf's Law and the Frequency of Characters or Words of Oracles

abstract: The article discusses the frequency of characters of Oracle,concluding that the frequency and the rank of a word or character is fit to Zipf-Mandelboit Law or Zipf's law with three parameters,and figuring out the parameters based on the frequency,and pointing out that what some researchers of Oracle call the assembling on the two ends is just a description by their impression about the Oracle data.


title: The Algebraic Duality Resolution at $p=2$

abstract: The goal of this paper is to develop some of the machinery necessary for doing $K(2)$-local computations in the stable homotopy category using duality resolutions at the prime $p=2$. The Morava stabilizer group $\mathbb{S}_2$ admits a norm whose kernel we denote by $\mathbb{S}_2^1$. The algebraic duality resolution is a finite resolution of the trivial $\mathbb{Z}_2[[\mathbb{S}_2^1]]$-module $\mathbb{Z}_2$ by modules induced from representations of finite subgroups of $\mathbb{S}_2^1$. Its construction is due to Goerss, Henn, Mahowald and Rezk. It is an analogue of their finite resolution of the trivial $\mathbb{Z}_3[[\mathbb{G}_2^1]]$-module $\mathbb{Z}_3$ at the prime $p=3$. The construction was never published and it is the main result in this paper. In the process, we give a detailed description of the structure of Morava stabilizer group $\mathbb{S}_2$ at the prime $2$. We also describe the maps in the algebraic duality resolution with the precision necessary for explicit computations.


title: Decomposition into Propagating and Evanescent Modes of Graphene Ribbons

abstract: Bulk modes (BM) are basis solutions to the Schr\"{o}dinger equation and they are useful in a number of physical problems. In the present work, we establish a complete set of BMs for graphene ribbons at arbitrary energy. We derive analytical expressions for these modes and systematically classify them into propagating or evanescent mode. We also demonstrate their uses in efficient electronic transport simulations of graphene-based electronic devices within both the mode-matching method and the Green's function framework. Explicit constructions of Green's functions for infinite and semi-infinite graphene ribbons are presented.