title: Gradient methods for variational optimization of projected entangled-pair states

abstract: We present a conjugate-gradient method for the ground-state optimization of projected entangled-pair states (PEPS) in the thermodynamic limit, as a direct implementation of the variational principle within the PEPS manifold. Our optimization is based on an efficient and accurate evaluation of the gradient of the global energy functional by using effective corner environments, and is robust with respect to the initial starting points. It has the additional advantage that physical and virtual symmetries can be straightforwardly implemented. We provide the tools to compute static structure factors directly in momentum space, as well as the variance of the Hamiltonian. We benchmark our method on Ising and Heisenberg models, and show a significant improvement on the energies and order parameters as compared to algorithms based on imaginary-time evolution.


title: The dynamics of a spinning particle in a linear in spin Hamiltonian approximation

abstract: We investigate for order and chaos the dynamical system of a spinning test particle of mass $m$ moving in the spacetime background of a Kerr black hole of mass M. This system is approximated in our investigation by the linear in spin Hamiltonian function provided in [E. Barausse, and A. Buonanno, Phys.Rev. D 81, 084024 (2010)]. We study the corresponding phase space by using 2D projections on a surface of section and the method of color and rotation on a 4D Poincar\'e section. Various topological structures coming from the non-integrability of the linear in spin Hamiltonian are found and discussed. Moreover, an interesting result is that from the value of the dimensionless spin $S/(m M)=10^{-4}$ of the particle and below, the impact of the non-integrability of the system on the motion of the particle seems to be negligible.


title: Lightning climatology of exoplanets and brown dwarfs guided by Solar System data

abstract: Clouds form on extrasolar planets and brown dwarfs where lightning could occur. Lightning is a tracer of atmospheric convection, cloud formation and ionization processes as known from the Solar System, and may be significant for the formation of prebiotic molecules. We study lightning climatology for the different atmospheric environments of Earth, Venus, Jupiter and Saturn. We present lightning distribution maps for Earth, Jupiter and Saturn, and flash densities for these planets and Venus, based on optical and/or radio measurements from the WWLLN and STARNET radio networks, the LIS/OTD satellite instruments, the Galileo, Cassini, New Horizons and Venus Express spacecraft. We also present flash densities calculated for several phases of two volcano eruptions, Eyjafjallaj\"okull's (2010) and Mt Redoubt's (2009). We estimate lightning rates for sample, transiting and directly imaged extrasolar planets and brown dwarfs. Based on the large variety of exoplanets, six categories are suggested for which we use the lightning occurrence information from the Solar System. We examine lightning energy distributions for Earth, Jupiter and Saturn. We discuss how strong stellar activity may support lightning activity. We provide a lower limit of the total number of flashes that might occur on transiting planets during their full transit as input for future studies. We find that volcanically very active planets might show the largest lightning flash densities. When applying flash densities of the large Saturnian storm from 2010/11, we find that the exoplanet HD 189733b would produce high lightning occurrence even during its short transit.

categories: astro-ph.EP astro-ph.SR physics.ao-ph physics.geo-ph


title: Approximation methods for the calculation of eigenvalues in ODE with periodic or anti periodic boundary conditions: Application to nanotubes

abstract: We compare three different methods to obtain solutions of Sturm-Liouville problems: a successive approximation method and two other iterative methods. We look for solutions with periodic or anti periodic boundary conditions. With some numerical test over the Mathieu equation, we compare the efficiency of these three methods. As an application, we make a numerical analysis on a model for carbon nanotubes.


title: On the Formation and Chemical Composition of Super Earths

abstract: Super Earths are the largest population of exoplanets and are seen to exhibit a rich diversity of compositions as inferred through their mean densities. Here we present a model that combines equilibrium chemistry in evolving disks with core accretion that tracks materials accreted onto planets during their formation. In doing so, we aim to explain why super Earths form so frequently and how they acquire such a diverse range of compositions. A key feature of our model is disk inhomogeneities, or planet traps, that act as barriers to rapid type-I migration. The traps we include are the dead zone, which can be caused by either cosmic ray or X-ray ionization, the ice line, and the heat transition. We find that in disks with sufficiently long lifetimes ($\gtrsim$ 4 Myr), all traps produce Jovian planets. In these disks, planet formation in the heat transition and X-ray dead zone produces hot Jupiters while the ice line and cosmic ray dead zones produce Jupiters at roughly 1 AU. Super Earth formation takes place within short-lived disks ($\lesssim$ 2 Myr), whereby the disks are photoevaporated while planets are in a slow phase of gas accretion. We find that super Earth compositions range from dry and rocky ($<$ 6 % ice by mass) to those with substantial water contents ($>$ 30 % ice by mass). The traps play a crucial role in our results, as they dictate where in the disk particular planets can accrete from, and what compositions they are able to acquire.


title: Comments related to infinite wedge representations

abstract: We study the infinite wedge representation and show how it is related to the universal extension of $g[t,t^{-1}]$ the loop algebra of a complex semi-simple Lie algebra $g$. We also give an elementary proof of the boson-fermion correspondence. Our approach to proving this result is based on a combinatorial construction with partitions combined with an application of the Murnaghan-Nakayama rule.


title: Improved Lower Bounds on Mutual Information Accounting for Nonlinear Signal-Noise Interaction

abstract: In fiber-optic communications, evaluation of mutual information (MI) is still an open issue due to the unavailability of an exact and mathematically tractable channel model. Traditionally, lower bounds on MI are computed by approximating the (original) channel with an auxiliary forward channel. In this paper, lower bounds are computed using an auxiliary backward channel, which has not been previously considered in the context of fiber-optic communications. Distributions obtained through two variations of the stochastic digital backpropagation (SDBP) algorithm are used as auxiliary backward channels and these bounds are compared with bounds obtained through the conventional digital backpropagation (DBP). Through simulations, higher information rates were achieved with SDBP, {which can be explained by the ability of SDBP to account for nonlinear signal--noise interactions


title: Measuring rare and exclusive Higgs boson decays into light resonances

abstract: We evaluate the LHC's potential of observing Higgs boson decays into light elementary or composite resonances through their hadronic decay channels. We focus on the Higgs boson production processes with the largest cross sections, $pp\to h$ and $pp\to h+\mathrm{jet}$, with subsequent decays $h \to ZA$ or $h\to Z\,\eta_c$, and comment on the production process $pp\to hZ$. By exploiting track-based jet substructure observables and extrapolating to $3000~\mathrm{fb}^{-1}$ we find ${\cal BR}(h \to ZA) \simeq {\cal BR}(h \to Z \eta_c) \lesssim 0.02$ at 95% CL. We interpret this limit in terms of the 2HDM Type 1. We find that searches for $h\to ZA$ are complementary to existing measurements and can constrain large parts of the currently allowed parameter space.


title: High-frequency asymptotic compression of dense BEM matrices for general geometries without ray tracing

abstract: Wave propagation and scattering problems in acoustics are often solved with boundary element methods. They lead to a discretization matrix that is typically dense and large: its size and condition number grow with increasing frequency. Yet, high frequency scattering problems are intrinsically local in nature, which is well represented by highly localized rays bouncing around. Asymptotic methods can be used to reduce the size of the linear system, even making it frequency independent, by explicitly extracting the oscillatory properties from the solution using ray tracing or analogous techniques. However, ray tracing becomes expensive or even intractable in the presence of (multiple) scattering obstacles with complicated geometries. In this paper, we start from the same discretization that constructs the fully resolved large and dense matrix, and achieve asymptotic compression by explicitly localizing the Green's function instead. This results in a large but sparse matrix, with a faster associated matrix-vector product and, as numerical experiments indicate, a much improved condition number. Though an appropriate localisation of the Green's function also depends on asymptotic information unavailable for general geometries, we can construct it adaptively in a frequency sweep from small to large frequencies in a way which automatically takes into account a general incident wave. We show that the approach is robust with respect to non-convex, multiple and even near-trapping domains, though the compression rate is clearly lower in the latter case. Furthermore, in spite of its asymptotic nature, the method is robust with respect to low-order discretizations such as piecewise constants, linears or cubics, commonly used in applications. On the other hand, we do not decrease the total number of degrees of freedom compared to a conventional classical discretization. The combination of the ...


title: New debris disks in nearby young moving groups

abstract: A significant fraction of nearby young moving group members harbor circumstellar debris dust disks. Due to their proximity and youth, these disks are attractive targets for studying the early evolution of debris dust and planetesimal belts. Here we present 70 and 160$\mu$m observations of 31 systems in the $\beta$ Pic moving group, and in the Tucana-Horologium, Columba, Carina and Argus associations, using the Herschel Space Observatory. None of these stars were observed at far-infrared wavelengths before. Our Herschel measurements were complemented by photometry from the WISE satellite for the whole sample, and by submillimeter/millimeter continuum data for one source, HD 48370. We identified six stars with infrared excess, four of them are new discoveries. By combining our new findings with results from the literature, we examined the incidence and general characteristics of debris disks around Sun-like members of the selected groups. With their dust temperatures of <45 K the newly identified disks around HD 38397, HD 48370, HD 160305, and BD-20 951 represent the coldest population within this sample. For HD 38397 and HD 48370, the emission is resolved in the 70$\mu$m PACS images, the estimated radius of these disks is ~90 au. Together with the well-known disk around HD 61005, these three systems represent the highest mass end of the known debris disk population around young G-type members of the selected groups. In terms of dust content, they resemble the hypothesized debris disk of the ancient Solar System.


title: All correlations must die: Assessing the significance of a stochastic gravitational-wave background in pulsar-timing arrays

abstract: We present two methods for determining the significance of a stochastic gravitational-wave (GW) background affecting a pulsar-timing array, where detection is based on evidence for quadrupolar spatial correlations between pulsars. Rather than constructing noise simulations, we eliminate the GWB spatial correlations in the true datasets to assess detection significance with all real data features intact. In our first method, we perform random phase shifts in the signal-model basis functions. This phase shifting eliminates signal phase coherence between pulsars, while keeping the statistical properties of the pulsar timing residuals intact. We then explore a method to null correlations between pulsars by using a "scrambled" overlap-reduction function in the signal model for the array. This scrambled function is orthogonal to what we expect of a real GW background signal. We demonstrate the efficacy of these methods using Bayesian model selection on a set of simulated datasets that contain a stochastic GW signal, timing noise, undiagnosed glitches, and uncertainties in the Solar system ephemeris. Finally, we introduce an overarching formalism under which these two techniques are naturally linked. These methods are immediately applicable to all current pulsar-timing array datasets, and should become standard tools for future analyses.


title: Frequency-dependent dielectric function of semiconductors with application to physisorption

abstract: The dielectric function is one of the most important quantities that describes the electrical and optical properties of solids. Accurate modeling of the frequency-dependent dielectric function has great significance in the study of the long-range van der Waals (vdW) interaction for solids and adsorption. In this work, we calculate the frequency-dependent dielectric functions of semiconductors and insulators using the $GW$ method with and without exciton effects, as well as efficient semilocal density functional theory (DFT), and compare these calculations with a model frequency-dependent dielectric function. We find that for semiconductors with moderate band gaps, the model dielectric functions, $GW$ values, and DFT calculations all agree well with each other. However, for insulators with strong exciton effects, the model dielectric functions have a better agreement with accurate $GW$ values than the DFT calculations, particularly in high-frequency region. To understand this, we repeat the DFT calculations with scissors correction, by shifting DFT Kohn-Sham energy gap to match the experimental band gap. We find that scissors correction only moderately improves the DFT dielectric function in low-frequency region. Based on the dielectric functions calculated with different methods, we make a comparative study by applying these dielectric functions to calculate the vdW coefficients ($C_3$ and $C_5$) for adsorption of rare-gas atoms on a variety of surfaces. We find that the vdW coefficients obtained with the nearly-free electron gas-based model dielectric function agree quite well with those obtained from the $GW$ dielectric function, in particular for adsorption on semiconductors, leading to an overall error of less than 7% for $C_3$ and 5% for $C_5$. This demonstrates the reliability of the model dielectric function for the study of physisorption.


title: Density Scaling and Decoupling in o-Terphenyl, Salol, and Dibutyphthalate

abstract: We present new viscosity and equation of state (EoS) results extending to high pressures for o-terphenyl, salol, and dibutylphthalate. Using these and data from the literature, we show that the three liquids all conform to density scaling; that is, their reduce viscosities and reorientational relaxation times are a function of the ratio of temperature and density with the latter raised to a constant. Moreover, the functional form of the dependence on this ratio is independent of the experimental probe of the dynamics. This means that there is no decoupling of the viscosities and relaxation times over the measured range of conditions. Previous literature at odds with these results were based on erroneous extrapolations of the EoS or problematic diamond anvil viscosity data. Thus, there are no exceptions to the experimental fact that every non-associated liquid complies with density scaling with an invariant scaling exponent.


title: Graphlike families of multiweights

abstract: Let ${\cal G}=(G,w)$ be a weighted graph , that is, a graph $G$ endowed with a function $w$ from the edge set of $G$ to the set of real numbers; for any subset $S$ of the vertex set of $G$, we define $D_S({\cal G})$ to be the minimum of the weights of the subgraphs of $G$ whose vertex set contains $S$; we call $D_S({\cal G})$ a multiweight of ${\cal G}$. Let $X$ be a finite set and let $\{D_S\}_{S \subset X, \; \sharp S \geq 2} $ be a family of positive real numbers. We find necessary and sufficient conditions for the family to be the family of multiweights of a positive-weighted graph with vertex set $X$. Moreover we study the analogous problem for trees. Finally, we find a criterion to say if there exists a nonnegative-weighted tree ${\cal T}$ with leaf set $X$ and such that $D_S ({\cal T})=D_S $ for any $S \subset X$.


title: Disease Trajectory Maps

abstract: Medical researchers are coming to appreciate that many diseases are in fact complex, heterogeneous syndromes composed of subpopulations that express different variants of a related complication. Time series data extracted from individual electronic health records (EHR) offer an exciting new way to study subtle differences in the way these diseases progress over time. In this paper, we focus on answering two questions that can be asked using these databases of time series. First, we want to understand whether there are individuals with similar disease trajectories and whether there are a small number of degrees of freedom that account for differences in trajectories across the population. Second, we want to understand how important clinical outcomes are associated with disease trajectories. To answer these questions, we propose the Disease Trajectory Map (DTM), a novel probabilistic model that learns low-dimensional representations of sparse and irregularly sampled time series. We propose a stochastic variational inference algorithm for learning the DTM that allows the model to scale to large modern medical datasets. To demonstrate the DTM, we analyze data collected on patients with the complex autoimmune disease, scleroderma. We find that DTM learns meaningful representations of disease trajectories and that the representations are significantly associated with important clinical outcomes.


title: Cognitive Control in the Controllable Connectome

abstract: Cognition is supported by neurophysiological processes that occur both in local anatomical neighborhoods and in distributed large-scale circuits. Recent evidence from network control theory suggests that white matter pathways linking large-scale brain regions provide a critical substrate constraining the ability of single areas to affect control on those processes. Yet, no direct evidence exists for a relationship between brain network controllability and cognitive control performance. Here, we address this gap by constructing structural brain networks from diffusion tensor imaging data acquired in 125 healthy adult individuals. We define a simplified model of brain dynamics and simulate network control to quantify modal and boundary controllability, which together describe complementary features of a region's theoretically predicted preference to drive the brain into different cognitive states. We observe that individual differences in these control features derived from structural connectivity are significantly correlated with individual differences in cognitive control performance, as measured by a continuous performance attention test, a color/shape switching task, the Stroop inhibition task, and a spatial n-back working memory task. Indeed, control hubs like anterior cingulate are distinguished from default mode and frontal association areas in terms of the relationship between their control properties and individual differences in cognitive function. These results provide the first empirical evidence that network control forms a fundamental mechanism of cognitive control.


title: Rainbow metric from quantum gravity: anisotropic cosmology

abstract: In this paper we present a construction of effective cosmological models which describe the propagation of a massive quantum scalar field on a quantum anisotropic cosmological spacetime. Each obtained effective model is represented by a rainbow metric in which particles of distinct momenta propagate on different classical geometries. Our analysis shows that upon certain assumptions and conditions on the parameters determining such anisotropic models, we surprisingly obtain a unique deformation parameter $\beta$ in the modified dispersion relation of the modes. Hence inducing an isotropic deformation despite the general starting considerations. We then ensure the recovery of the dispersion relation realized in the isotropic case, studied in [arXiv:1412.6000], when some proper symmetry constraints are imposed, and we estimate the value of the deformation parameter for this case in loop quantum cosmology context.


title: Object Boundary Detection and Classification with Image-level Labels

abstract: Semantic boundary and edge detection aims at simultaneously detecting object edge pixels in images and assigning class labels to them. Systematic training of predictors for this task requires the labeling of edges in images which is a particularly tedious task. We propose a novel strategy for solving this task, when pixel-level annotations are not available, performing it in an almost zero-shot manner by relying on conventional whole image neural net classifiers that were trained using large bounding boxes. Our method performs the following two steps at test time. Firstly it predicts the class labels by applying the trained whole image network to the test images. Secondly, it computes pixel-wise scores from the obtained predictions by applying backprop gradients as well as recent visualization algorithms such as deconvolution and layer-wise relevance propagation. We show that high pixel-wise scores are indicative for the location of semantic boundaries, which suggests that the semantic boundary problem can be approached without using edge labels during the training phase.


title: Three dimensional graph drawing with fixed vertices and one bend per edge

abstract: We prove that for every graph $G$, given fixed locations for the vertices of $G$ in $\mathbb{Z}^3$, there is a three-dimensional grid-drawing of $G$ with one bend per edge. The best previous bound was three bends per edge.


title: Multiple mixing and parabolic divergence in smooth area-preserving flows on higher genus surfaces

abstract: We consider typical area preserving flows on higher genus surfaces and prove that the flow restricted to mixing minimal components is mixing of all orders, thus answering affimatively to Rohlin's multiple mixing question in this context. The main tool is a variation of the Ratner property (a property originally proved by Ratner for the horocycle flow), i.e. the switchable Ratner property introduced by Fayad and Kanigowski for special flows over rotations. This property, which is of independent interest, provides a quantitative description of the parabolic behaviour of these flows and has implications to joinings classification. The main result is formulated in the language of special flows over interval exchange transformations with asymmetric logarithmic singularities. We also prove a strengthening of one of Fayad and Kanigowski's main results, by showing that Arnold's flows are mixing of all orders for almost every location of the singularities.


title: A Semi-Definite Programming approach to low dimensional embedding for unsupervised clustering

abstract: This paper proposes a variant of the method of Gu\'edon and Verhynin for estimating the cluster matrix in the Mixture of Gaussians framework via Semi-Definite Programming. A clustering oriented embedding is deduced from this estimate. The procedure is suitable for very high dimensional data because it is based on pairwise distances only. Theoretical garantees are provided and an eigenvalue optimisation approach is proposed for computing the embedding. The performance of the method is illustrated via Monte Carlo experiements and comparisons with other embeddings from the literature.


title: Non-universal Z' models with protected flavour-changing interactions

abstract: We define a new class of Z' models with neutral flavour-changing interactions at tree level in the down-quark sector. They are related in an exact way to elements of the quark mixing matrix due to an underlying flavoured U(1)' gauge symmetry, rendering these models particularly predictive. The same symmetry implies lepton-flavour non-universal couplings, fully determined by the gauge structure of the model. Our models allow to address several presently observed deviations from the SM and specific correlations among the new physics contributions to the Wilson coefficients C9,10(l) can be tested in b->sll transitions. We furthermore predict lepton-universality violations in Z' decays, testable at the LHC.


title: Reorientation of the bicollinear antiferromagnetic structure at the surface of Fe$_{1+y}$Te bulk and thin films

abstract: Establishing the relation between the ubiquitous antiferromagnetism in the non-superconducting parent compounds of unconventional superconductors and their superconducting phase is believed to be important for the understanding of the complex physics in these materials. Going from the bulk systems to thin films strongly affects the phase diagram of unconventional superconductors. For Fe$_{1+y}$Te, the parent compound of the Fe$_{1+y}$Se$_{1-x}$Te$_x$ superconductors, bulk sensitive neutron diffraction has revealed an in-plane oriented bicollinear antiferromagnetic structure. Here, we show by spin-resolved scanning tunneling microscopy that on the surfaces of bulk Fe$_{1+y}$Te, as well as on thin films grown on the topological insulator Bi$_2$Te$_3$, the spin direction is canted both away from the surface plane and from the high-symmetry directions of the surface unit cell, while keeping the bicollinear magnetic structure. Our results demonstrate that the magnetism at the Fe-chalcogenide surface markedly deviates from a simple in-plane oriented bicollinear antiferromagnetic structure, which implies that the pairing at the surface of the related superconducting compounds might be different from that in the bulk.


title: Small coherence implies the weak Null Space Property

abstract: In the Compressed Sensing community, it is well known that given a matrix $X \in \mathbb R^{n\times p}$ with $\ell_2$ normalized columns, the Restricted Isometry Property (RIP) implies the Null Space Property (NSP). It is also well known that a small Coherence $\mu$ implies a weak RIP, i.e. the singular values of $X_T$ lie between $1-\delta$ and $1+\delta$ for "most" index subsets $T \subset \{1,\ldots,p\}$ with size governed by $\mu$ and $\delta$. In this short note, we show that a small Coherence implies a weak Null Space Property, i.e. $\Vert h_T\Vert_2 \le C \ \Vert h_{T^c}\Vert_1/\sqrt{s}$ for most $T \subset \{1,\ldots,p\}$ with cardinality $|T|\le s$. We moreover prove some singular value perturbation bounds that may also prove useful for other applications.


title: A multilayer approach for price dynamics in financial markets

abstract: We introduce a new Self-Organized Criticality (SOC) model for simulating price evolution in an artificial financial market, based on a multilayer network of traders. The model also implements, in a quite realistic way with respect to previous studies, the order book dy- namics, by considering two assets with variable fundamental prices. Fat tails in the probability distributions of normalized returns are observed, together with other features of real financial markets.


title: Toward unbiased estimations of the statefinder parameters

abstract: With the use of simulated supernova catalogs, we show that the statefinder parameters turn out to be poorly and biased estimated by standard cosmography. To this end, we compute their standard deviations and several bias statistics on cosmologies near the concordance model, demonstrating that these are very large, making standard cosmography unsuitable for future and wider compilations of data. To overcome this issue, we propose a new method that consists in introducing the series of the Hubble function into the luminosity distance, instead of considering the usual direct Taylor expansions of the luminosity distance. Moreover, in order to speed up the numerical computations, we estimate the coefficients of our expansions in a hierarchical manner, in which the order of the expansion depends on the redshift of every single piece of data. In addition, we propose two hybrids methods that incorporates standard cosmography at low redshifts. The methods presented here perform better than the standard approach of cosmography both in the errors and bias of the estimated statefinders. We further propose a one-parameter diagnostic to reject non-viable methods in cosmography.


title: The Role of Magnetic Helicity in the Structure and Heating of the Sun's Corona

abstract: Two of the most important features of the solar atmosphere are its hot, smooth coronal loops and the concentrations of magnetic shear, known as filament channels, that reside above photospheric polarity inversion lines (PILs). The shear observed in filament channels represents magnetic helicity, while the smoothness of the coronal loops indicates an apparent lack of magnetic helicity in the rest of the corona. At the same time, models that attempt to explain the high temperatures observed in these coronal loops require magnetic energy, in the form of twist, to be injected at the photosphere. In addition to magnetic energy, this twist also represents magnetic helicity. Unlike magnetic energy, magnetic helicity is conserved under reconnection, and is consequently expected to accumulate and be observed in the corona. However, filament channels, rather than the coronal loops, are the locations in the corona where magnetic helicity is observed, and it manifests itself in the form of shear, rather than twist. This naturally raises the question: if magnetic helicity needs to be injected to heat coronal loops, why is it only observed in filament channels, while coronal loops are observed to be laminar and smooth? This thesis addresses this question using a series of numerical simulations that demonstrate that magnetic helicity is transported throughout the solar corona by magnetic reconnection in such a way that it accumulates above PILs, forming filament channels, and leaving the rest of the corona generally smooth. In the process, it converts magnetic energy into heat, accounting for the large observed temperatures. This thesis presents a model for the formation of filament channels in the solar corona and the presence of smooth, hot coronal loops, and shows how the transport of magnetic helicity throughout the solar corona by magnetic reconnection is responsible for both of these phenomena.


title: Model-Free Trajectory-based Policy Optimization with Monotonic Improvement

abstract: Many of the recent trajectory optimization algorithms alternate between linear approximation of the system dynamics around the mean trajectory and conservative policy update. One way of constraining the policy change is by bounding the Kullback-Leibler (KL) divergence between successive policies. These approaches already demonstrated great experimental success in challenging problems such as end-to-end control of physical systems. However, the linear approximation of the system dynamics can introduce a bias in the policy update and prevent convergence to the optimal policy. In this article, we propose a new model-free trajectory-based policy optimization algorithm with guaranteed monotonic improvement. The algorithm backpropagates a local, quadratic and time-dependent \qfunc~learned from trajectory data instead of a model of the system dynamics. Our policy update ensures exact KL-constraint satisfaction without simplifying assumptions on the system dynamics. We experimentally demonstrate on highly non-linear control tasks the improvement in performance of our algorithm in comparison to approaches linearizing the system dynamics. In order to show the monotonic improvement of our algorithm, we additionally conduct a theoretical analysis of our policy update scheme to derive a lower bound of the change in policy return between successive iterations.


title: On the integrability of the isotropic almost complex structures and harmonic unit vector fields

abstract: Aguilar introduced isotropic almost complex structures $J_{\delta , \sigma}$ on the tangent bundle of a Riemannian manifold $(M, g)$. In this paper, some results will be obtained on the integrability of these structures. These structures with the Liouville 1-form define a class of Riemannian metrics $g_{\delta , \sigma}$ on $T M$ which are a generalization of the Sasaki metric. Moreover, the notion of a harmonic unit vector field is introduced with respect to these metrics like as the Sasaki metric and the necessary and sufficient conditions for a unit vector field to be a harmonic unit vector field are obtained.


title: Inflation Model Selection meets Dark Radiation

abstract: We investigate how inflation model selection is affected by the presence of additional free-streaming relativistic degrees of freedom, i.e. dark radiation. We perform a full Bayesian analysis of both inflation parameters and cosmological parameters taking reheating into account self-consistently. We compute the Bayesian evidence for a few representative inflation scenarios in both the standard $\Lambda$CDM model and an extension including dark radiation parametrised by its effective number of relativistic species $N_\mathrm{eff}$. Using a minimal dataset (Planck low-$\ell$ polarisation, temperature power spectrum and lensing reconstruction), we find that the observational status of most inflationary models is unchanged. The exceptions are potentials such as power-law inflation that predict large values for the scalar spectral index that can only be realised when $N_\mathrm{eff}$ is allowed to vary. Adding baryon acoustic oscillations data and the B-mode data from BICEP2/Keck makes power-law inflation disfavoured, while adding local measurements of the Hubble constant $H_0$ makes power-law inflation slightly favoured compared to the best single-field plateau potentials. This illustrates how the dark radiation solution to the $H_0$ tension would have deep consequences for inflation model selection.


title: Multiparty Delegated Quantum Computing

abstract: Quantum computing has seen tremendous progress in the past years. However, due to limitations in scalability of quantum technologies, it seems that we are far from constructing universal quantum computers for everyday users. A more feasible solution is the delegation of computation to powerful quantum servers on the network. This solution was proposed in previous studies of Blind Quantum Computation, with guarantees for both the secrecy of the input and of the computation being performed. In this work, we further develop this idea of computing over encrypted data, to propose a multiparty delegated quantum computing protocol in the measurement-based quantum computing framework.


title: Reconnection-Driven Coronal-Hole Jets with Gravity and Solar Wind

abstract: Coronal-hole jets occur ubiquitously in solar coronal holes, at EUV and X-ray bright points associated with intrusions of minority magnetic polarity. The embedded-bipole model for these jets posits that they are driven by explosive, fast reconnection between the stressed closed field of the embedded bipole and the open field of the surrounding coronal hole. Previous numerical studies in Cartesian geometry, assuming uniform ambient magnetic field and plasma while neglecting gravity and solar wind, demonstrated that the model is robust and can produce jet-like events in simple configurations. We have extended these investigations by including spherical geometry, gravity, and solar wind in a nonuniform, coronal hole-like ambient atmosphere. Our simulations confirm that the jet is initiated by the onset of a kink-like instability of the internal closed field, which induces a burst of reconnection between the closed and external open field, launching a helical jet. Our new results demonstrate that the jet propagation is sustained through the outer corona, in the form of a traveling nonlinear Alfven wave front trailed by slower-moving plasma density enhancements that are compressed and accelerated by the wave. This finding agrees well with observations of white-light coronal-hole jets, and can explain microstreams and torsional Alfven waves detected in situ in the solar wind. We also use our numerical results to deduce scaling relationships between properties of the coronal source region and the characteristics of the resulting jet, which can be tested against observations.


title: Tighter bounds lead to improved classifiers

abstract: The standard approach to supervised classification involves the minimization of a log-loss as an upper bound to the classification error. While this is a tight bound early on in the optimization, it overemphasizes the influence of incorrectly classified examples far from the decision boundary. Updating the upper bound during the optimization leads to improved classification rates while transforming the learning into a sequence of minimization problems. In addition, in the context where the classifier is part of a larger system, this modification makes it possible to link the performance of the classifier to that of the whole system, allowing the seamless introduction of external constraints.


title: Bias-dependent local structure of water molecule at a metallic interface

abstract: Understanding the local structure of water at the interfaces of metallic electrodes is a key problem in aqueous-based electrochemistry. Nevertheless, a realistic simulation of such setup is challenging, particularly when the electrodes are maintained at different potentials. To correctly compute the effect of an external bias potential applied to truly semi-infinite surfaces, we combine Density Functional Theory (DFT) and Non-Equilibrium Green's Functions (NEGF) methods. This framework allows for the out-of-equilibrium calculation of forces and dynamics, and directly correlates to the chemical potential of the electrodes, which is the one introduced experimentally. In this work, we apply this methodology to study the electronic properties and atomic forces of one water molecule at the interface of gold surface. We find that the water molecule tends to align its dipole moment with the electric field, and it is either repelled or attracted to the metal depending on the sign and magnitude of the applied bias, in an asymmetric fashion.


title: Extinction Law at a Distance up to 25 kpc toward the Galactic Poles

abstract: Photometry from the Tycho-2, 2MASS, and WISE catalogues for clump and branch giants at a distance up to 25 kpc toward the Galactic poles has allowed the variations of various characteristics of the infrared interstellar extinction law with distance to be analyzed. The results obtained by the extinction law extrapolation method are consistent for different classes of stars and different characteristics as well as with previous studies. The conventional extinction law with a low infrared extinction is characteristic of only a thin layer no farther than 100 pc from the Galactic plane and of two thin layers near $Z=-600$ and $+500$ pc. Far from the Galactic plane, in the Galactic halo, the infrared extinction law is different: the extinction in the $Ks$, $W1$, $W2$, $W3$, and $W4$ bands is, respectively, $0.17$, $0.16$, $0.16$, $0.07$, and $0.03$ of the extinction in the $V$ band. The accuracy of these coefficients is $0.03$. If the extinction law reflects primarily the grain size distribution, then the fraction of large dust grains far from the Galactic plane is greater than that in the circumsolar interstellar medium.


title: Next Generation Robotics

abstract: The National Robotics Initiative (NRI) was launched 2011 and is about to celebrate its 5 year anniversary. In parallel with the NRI, the robotics community, with support from the Computing Community Consortium, engaged in a series of road mapping exercises. The first version of the roadmap appeared in September 2009; a second updated version appeared in 2013. While not directly aligned with the NRI, these road-mapping documents have provided both a useful charting of the robotics research space, as well as a metric by which to measure progress. This report sets forth a perspective of progress in robotics over the past five years, and provides a set of recommendations for the future. The NRI has in its formulation a strong emphasis on co-robot, i.e., robots that work directly with people. An obvious question is if this should continue to be the focus going forward? To try to assess what are the main trends, what has happened the last 5 years and what may be promising directions for the future a small CCC sponsored study was launched to have two workshops, one in Washington DC (March 5th, 2016) and another in San Francisco, CA (March 11th, 2016). In this report we brief summarize some of the main discussions and observations from those workshops. We will present a variety of background information in Section 2, and outline various issues related to progress over the last 5 years in Section 3. In Section 4 we will outline a number of opportunities for moving forward. Finally, we will summarize the main points in Section 5.


title: Performance of spatial Multi-LRU caching under traffic with temporal locality

abstract: In this work a novel family of decentralised caching policies for wireless networks is introduced, referred to as spatial multi-LRU. These improve cache-hit probability by exploiting multi-coverage. Two variations are proposed, the multi-LRU-One and -All, which differ in the number of replicas inserted in the covering edge-caches. The evaluation is done under spatial traffic that exhibits temporal locality, with varying content catalogue and dependent demands. The performance metric is hit probability and the policies are compared to (1) the single-LRU and (2) an upper bound for all centralised policies with periodic popularity updates. Numerical results show the multi-LRU policies outperform both comparison policies. The reason is their passive adaptability to popularity changes. Between the -One and -All variation, which one is preferable strongly depends on the available storage space and on traffic characteristics. The performance also depends on the popularity shape.


title: Spherical collapse of dark matter haloes in tidal gravitational fields

abstract: We study the spherical collapse model in the presence of external gravitational tidal shear fields for different dark energy scenarios and investigate the impact on the mass function and cluster number counts. While previous studies of the influence of shear and rotation on $\delta_\mathrm{c}$ have been performed with heuristically motivated models, we try to avoid this model dependence and sample the external tidal shear values directly from the statistics of the underlying linearly evolved density field based on first order Lagrangian perturbation theory. Within this self-consistent approach, in the sense that we restrict our treatment to scales where linear theory is still applicable, only fluctuations larger than the scale of the considered objects are included into the sampling process which naturally introduces a mass dependence of $\delta_\mathrm{c}$. We find that shear effects are predominant for smaller objects and at lower redshifts, i. e. the effect on $\delta_\mathrm{c}$ is at or below the percent level for the $\Lambda$CDM model. For dark energy models we also find small but noticeable differences, similar to $\Lambda$CDM. The virial overdensity $\Delta_\mathrm{V}$ is nearly unaffected by the external shear. The now mass dependent $\delta_c$ is used to evaluate the mass function for different dark energy scenarios and afterwards to predict cluster number counts, which indicate that ignoring the shear contribution can lead to biases of the order of $1\sigma$ in the estimation of cosmological parameters like $\Omega_\mathrm{m}$, $\sigma_8$ or $w$.


title: The maximum size of a partial spread II: Upper bounds

abstract: Let $n$ and $t$ be positive integers with $t<n$, and let $q$ be a prime power. A partial $(t-1)$-spread of ${\rm PG}(n-1,q)$ is a set of $(t-1)$-dimensional subspaces of ${\rm PG}(n-1,q)$ that are pairwise disjoint. Let $r\equiv n\pmod{t}$ with $0\leq r<t$, and let ${\theta}_i=(q^i-1)/(q-1)$. We essentially prove that if $2\leq r<t\leq {\theta}_r$, then the maximum size of a partial $(t-1)$-spread of ${\rm PG}(n-1,q)$ is bounded from above by $({\theta}_n-{\theta}_{t+r})/{\theta}_t+q^r-(q-1)(t-3)+1$. We actually give tighter bounds when certain divisibility conditions are satisfied. These bounds improve on the previously known upper bound for the maximum size partial ($t-1$)-spreads of ${\rm PG}(n-1,q)$; for instance, when $\lceil\frac{{\theta}_r}{2}\rceil+4\leq t\leq {\theta}_r$ and $q>2$. The exact value of the maximum size partial $(t-1)$-spread has been recently determined for $t>{\theta}_r$ by the authors of this paper (see N\u{a}stase-Sissokho [21]).


title: Weyl node with random vector potential

abstract: We study Weyl semimetals in the presence of generic disorder, consisting of a random vector potential as well as a random scalar potential. We derive renormalization group flow equations to second order in the disorder strength. These flow equations predict a disorder-induced phase transition between a pseudo-ballistic weak-disorder phase and a diffusive strong-disorder phase for sufficiently strong random scalar potential or for a pure three-component random vector potential. We verify these predictions using a numerical study of the density of states near the Weyl point and of quantum transport properties at the Weyl point. In contrast, for a pure single-component random vector potential the diffusive strong-disorder phase is absent.


title: Renormalization group functions of $\phi^4$ theory in the MS-scheme to six loops

abstract: Subdivergences constitute a major obstacle to the evaluation of Feynman integrals and an expression in terms of finite quantities can be a considerable advantage for both analytic and numeric calculations. We report on our implementation of the suggestion by F. Brown and D. Kreimer, who proposed to use a modified BPHZ scheme where all counterterms are single-scale integrals. Paired with parametric integration via hyperlogarithms, this method is particularly well suited for the computation of renormalization group functions and easily automated. As an application of this approach we compute the 6-loop beta function and anomalous dimensions of the $\phi^4$ model.


title: Synchronization of electrically coupled stochastic magnetic oscillators induced by thermal and electrical noise

abstract: Superparamagnetic tunnel junctions are nanostructures that auto-oscillate stochastically under the effect of thermal noise. Recent works showed that despite their stochasticity, such junctions possess a capability to synchronize to subthreshold voltage drives, in a way that can be enhanced or controlled by adding noise. In this work, we investigate a system composed of two electrically coupled junctions, connected in series to a periodic voltage source. We make use of numerical simulations and of an analytical model to demonstrate that both junctions can be phase-locked to the drive, in phase or in anti-phase. This synchronization phenomenon can be controlled by both thermal and electrical noises, although the two types of noises induce qualitatively different behaviors. Namely, thermal noise can stabilize a regime where one junction is phase-locked to the drive voltage while the other is blocked in one state. On the contrary, electrical noise causes the junctions to have highly correlated behaviors and thus cannot induce the latter. These results open the way for the design of superparamagnetic tunnel junctions that can perform computation through synchronization, and which harvest the largest part of their energy consumption from thermal noise.


title: Unique Measure for Time-Dependent Random Dynamical Systems

abstract: This paper proves the uniqueness of measure for the two-dimensional Navier-Stokes equations under a random kick-force and a time-dependent deterministic force. By extending a result for uniqueness of measure for time-homogeneous Markov processes to the time-inhomogeneous case, it is shown that the measures are exponentially mixing for the 2D Navier-Stokes equations on the sphere.


title: Leading two-loop corrections to the Higgs boson masses in SUSY models with Dirac gauginos

abstract: We compute the two-loop O(as*at) corrections to the Higgs boson masses in supersymmetric extensions of the Standard Model with Dirac gaugino masses. We rely on the effective-potential technique, allow for both Dirac and Majorana mass terms for the gluinos, and compute the corrections in both the DRbar and on-shell renormalisation schemes. We give detailed results for the MDGSSM and the MRSSM, and simple approximate formulae valid in the decoupling limit for all currently-studied variants of supersymmetric models with Dirac gluinos. These results represent the first explicit two-loop calculation of Higgs boson masses in supersymmetric models beyond the MSSM and the NMSSM.


title: HiPERCAM: A high-speed, quintuple-beam CCD camera for the study of rapid variability in the Universe

abstract: HiPERCAM is a high-speed camera for the study of rapid variability in the Universe. The project is funded by a 3.5MEuro European Research Council Advanced Grant. HiPERCAM builds on the success of our previous instrument, ULTRACAM, with very significant improvements in performance thanks to the use of the latest technologies. HiPERCAM will use 4 dichroic beamsplitters to image simultaneously in 5 optical channels covering the u'g'r'i'z' bands. Frame rates of over 1000 per second will be achievable using an ESO CCD controller (NGC), with every frame GPS timestamped. The detectors are custom-made, frame-transfer CCDs from e2v, with 4 low-noise (2.5e-) outputs, mounted in small thermoelectrically-cooled heads operated at 180 K, resulting in virtually no dark current. The two reddest CCDs will be deep-depletion devices with anti-etaloning, providing high quantum efficiencies across the red part of the spectrum with no fringing. The instrument will also incorporate scintillation noise correction via the conjugate-plane photometry technique. The opto-mechanical chassis will make use of additive manufacturing techniques in metal to make a light-weight, rigid and temperature-invariant structure. First light is expected on the 4.2m William Herschel Telescope on La Palma in 2017 (on which the field of view will be 10' with a 0.3"/pixel scale), with subsequent use planned on the 10.4m Gran Telescopio Canarias on La Palma (on which the field of view will be 4' with a 0.11"/pixel scale) and the 3.5m New Technology Telescope in Chile.


title: Observation and analysis of the new W-type W UMa eclipsing binary VSX J053024.8+842243

abstract: Using multicolour photometry we have confirmed the binary nature of the new W-type W UMa eclipsing binary VSX J053024.8+842243 and established its primary eclipse ephemeris to be HJD = 2455924.38150(26) + 0.4322929(1) * E. Using the light curve modelling code PHOEBE and published data on the evolution of W-type contact binaries we found the primary and secondary components to have masses 0.50 Msun and 1.44 Msun, radii 0.87 Rsun and 1.42 Rsun, luminosities 0.98 Lsun and 1.91 Lsun, temperatures 6145 K and 5702 K and binary orbit inclination 59.4{\deg}. We found the distance to the binary to be 511 parsec, its E(B-V) colour excess 0.04 and its intrinsic (B-V) colour index 0.62. A low resolution spectrum corrected for interstellar reddening confirmed its spectral type as G2V.


title: True Error Control for the Localized Reduced Basis Method for Parabolic Problems

abstract: We present an abstract framework for a posteriori error estimation for approximations of scalar parabolic evolution equations, based on elliptic reconstruction techniques [10, 9, 3, 5]. In addition to its original application (to derive error estimates on the discretization error), we extend the scope of this framework to derive offline/online decomposable a posteriori estimates on the model reduction error in the context of Reduced Basis (RB) methods. In addition, we present offline/online decomposable a posteriori error estimates on the full approximation error (including discretization as well as model reduction error) in the context of the localized RB method [14]. Hence, this work generalizes the localized RB method with true error certification to parabolic problems. Numerical experiments are given to demonstrate the applicability of the approach.


title: Orbital angular momentum of scalar field generated by gravitational scatterings

abstract: It has been expected that astronomical observations to detect the orbital angular momenta of electromagnetic waves may give us a new insight into astrophysics. Previous works pointed out the possibility that a rotating black hole can produce orbital angular momenta of electromagnetic waves through gravitational scattering, and the spin parameter of the black hole can be measured by observing them. However, the mechanism how the orbital angular momentum of the electromagnetic wave is generated by the gravitational scattering has not been clarified sufficiently. In this paper, in order to understand it from a point of view of gravitational lensing effects, we consider an emitter which radiates a spherical wave of the real massless scalar field and study the deformation of the scalar wave by the gravitational scattering due to a black hole by invoking the geometrical optics approximation. We show that the frame dragging caused by the rotating black hole is not a necessary condition for generating the orbital angular momentum of the scalar wave. However, its components parallel to the direction cosines of images appear only if the black hole is rotating.


title: Range-separated tensor formats for numerical modeling of many-particle interaction potentials

abstract: We introduce and analyze the new range-separated (RS) canonical/Tucker tensor format which aims for numerical modeling of the 3D long-range interaction potentials in multi-particle systems. The main idea of the RS tensor format is the independent grid-based low-rank representation of the localized and global parts in the target tensor which allows the efficient numerical approximation of $N$-particle interaction potentials. The single-particle reference potential like $1/\|x\|$ is split into a sum of localized and long-range low-rank canonical tensors represented on a fine 3D $n\times n\times n$ Cartesian grid. The smoothed long-range contribution to the total potential sum is represented on the 3D grid in $O(n)$ storage via the low-rank canonical/Tucker tensor. We prove that the Tucker rank parameters depend only logarithmically on the number of particles $N$ and the grid-size $n$. Agglomeration of the short range part in the sum is reduced to an independent treatment of $N$ localized terms with almost disjoint effective supports, calculated in $O(N)$ operations. Thus, the cumulated sum of short range clusters is parametrized by a single low-rank canonical reference tensor with a local support, accomplished by a list of particle coordinates and their charges. The RS canonical/Tucker tensor representations reduce the cost of multi-linear algebraic operations on the 3D potential sums arising in modeling of multi-dimensional data by radial basis functions, say, in computation of the electrostatic potential of a protein, in 3D integration and convolution transforms, computation of gradients, forces and the interaction energy of a many-particle systems, and in low parametric fitting of multi-dimensional scattered data by reducing all of them to 1D calculations.


title: Simulation Results on Selector Adaptation in Behavior Trees

abstract: Behavior trees (BTs) emerged from video game development as a graphical language for modeling intelligent agent behavior. However as initially implemented, behavior trees are static plans. This paper adds to recent literature exploring the ability of BTs to adapt to their success or failure in achieving tasks. The "Selector" node of a BT tries alternative strategies (its children) and returns success only if all of its children return failure. This paper studies several means by which Selector nodes can learn from experience, in particular, learn conditional probabilities of success based on sensor information, and modify the execution order based on the learned iformation. Furthermore, a "Greedy Selector" is studied which only tries the child having the highest success probability. Simulation results indicate significantly increased task performance, especially when frequentist probability estimate is conditioned on sensor information. The Greedy selector was ineffective unless it was preceded by a period of training in which all children were exercised.


title: Scalar Gravitational Waves in the Effective Theory of Gravity

abstract: As a low energy effective field theory, classical General Relativity receives an infrared relevant modification from the conformal trace anomaly of the energy-momentum tensor of massless, or nearly massless, quantum fields. The local form of the effective action associated with the trace anomaly is expressed in terms of a dynamical scalar field that couples to the conformal factor of the spacetime metric, allowing it to propagate over macroscopic distances. Linearized around flat spacetime, this semi-classical EFT admits scalar gravitational wave solutions in addition to the transversely polarized tensor waves of the classical Einstein theory. The amplitude of the scalar wave modes, as well as their energy and energy flux which are positive and contain a monopole moment, are computed. Astrophysical sources for scalar gravitational waves are considered, with the excited gluonic condensates in the interiors of neutron stars in merger events with other compact objects likely to provide the strongest burst signals.


title: The new challenges of multiplex networks: measures and models

abstract: What do societies, the Internet, and the human brain have in common? They are all examples of complex relational systems, whose emerging behaviours are largely determined by the non-trivial networks of interactions among their constituents, namely individuals, computers, or neurons, rather than the properties of the units themselves. In the last two decades, network scientists have proposed models of increasing complexity to better understand real-world systems. Only recently we have realised that multiplexity, i.e. the coexistence of several types of interactions among the constituents of a complex system, is responsible for substantial qualitative and quantitative differences in the type and variety of behaviours that a complex system can exhibit. As a consequence, multilayer and multiplex networks have become a hot topic in complexity science. Here we provide an overview of some of the measures proposed so far to characterise the structure of multiplex networks, and a selection of models aiming at reproducing those structural properties and quantifying their statistical significance. Focusing on a subset of relevant topics, this brief review is a quite comprehensive introduction to the most basic tools for the analysis of multiplex networks observed in the real-world. The wide applicability of multiplex networks as a framework to model complex systems in different fields, from biology to social sciences, and the colloquial tone of the paper will make it an interesting read for researchers working on both theoretical and experimental analysis of networked systems.

categories: physics.soc-ph cond-mat.dis-nn cond-mat.stat-mech cs.SI


title: Penambahan emosi menggunakan metode manipulasi prosodi untuk sistem text to speech bahasa Indonesia

abstract: Adding an emotions using prosody manipulation method for Indonesian text to speech system. Text To Speech (TTS) is a system that can convert text in one language into speech, accordance with the reading of the text in the language used. The focus of this research is a natural sounding concept, the make "humanize" for the pronunciation of voice synthesis system Text To Speech. Humans have emotions / intonation that may affect the sound produced. The main requirement for the system used Text To Speech in this research is eSpeak, the database MBROLA using id1, Human Speech Corpus database from a website that summarizes the words with the highest frequency (Most Common Words) used in a country. And there are 3 types of emotional / intonation designed base. There is a happy, angry and sad emotion. Method for develop the emotional filter is manipulate the relevant features of prosody (especially pitch and duration value) using a predetermined rate factor that has been established by analyzing the differences between the standard output Text To Speech and voice recording with emotional prosody / a particular intonation. The test results for the perception tests of Human Speech Corpus for happy emotion is 95 %, 96.25 % for angry emotion and 98.75 % for sad emotions. For perception test system carried by intelligibility and naturalness test. Intelligibility test for the accuracy of sound with the original sentence is 93.3%, and for clarity rate for each sentence is 62.8%. For naturalness, accuracy emotional election amounted to 75.6 % for happy emotion, 73.3 % for angry emotion, and 60 % for sad emotions. ----- Text To Speech (TTS) merupakan suatu sistem yang dapat mengonversi teks dalam format suatu bahasa menjadi ucapan sesuai dengan pembacaan teks dalam bahasa yang digunakan.


title: The magnetic precursor of the pressure-induced superconductivity in Fe-ladder compound

abstract: The pressure effects on the antiferromagentic orders in iron-based ladder compounds CsFe$_2$Se$_3$ and BaFe$_2$S$_3$ have been studied using neutron diffraction. With identical crystal structure and similar magnetic structures, the two compounds exhibit highly contrasting magnetic behaviors under moderate external pressures. In CsFe$_2$Se$_3$ the ladders are brought much closer to each other by pressure, but the stripe-type magnetic order shows no observable change. In contrast, the stripe order in BaFe$_2$S$_3$, undergoes a quantum phase transition where an abrupt increase of N$\acute{e}$el temperature by more than 50$\%$ occurs at about 1 GPa, accompanied by a jump in the ordered moment. With its spin structure unchanged, BaFe$_2$S$_3$ enters an enhanced magnetic phase that bears the characteristics of an orbital selective Mott phase, which is the true neighbor of superconductivity emerging at higher pressures.


title: The habitability of the Universe through 13 billion years of cosmic time

abstract: The field of astrobiology has made tremendous progress in modelling galactic-scale habitable zones which offer a stable environment for life to form and evolve in complexity. Recently, this idea has been extended to cosmological scales by studies modelling the habitability of the local Universe in its entirety (e.g. Dayal et al. 2015; Li & Zhang 2015). However, all of these studies have solely focused on estimating the potentially detrimental effects of either Type II supernovae (SNII) or Gamma Ray Bursts (GRBs), ignoring the contributions from Type Ia supernovae (SNIa) and active galactic nuclei (AGN). In this study we follow two different approaches, based on (i) the amplitude of deleterious radiation and (ii) the total planet-hosting volume irradiated by deleterious radiation. We simultaneously track the contributions from the key astrophysical sources (SNII, SNIa, AGN and GRBs) for the entire Universe, for both scenarios, to determine its habitability through 13.8 billion years of cosmic time. We find that SNII dominate the total radiation budget and the volume irradiated by deleterious radiation at any cosmic epoch closely followed by SNIa (that contribute half as much as SNII), with GRBs and AGN making up a negligible portion (<1%). Secondly, as a result of the total mass in stars (or the total number of planets) slowly building-up with time and the total deleterious radiation density, and volume affected, falling-off after the first 3 billion years, we find that the Universe has steadily increased in habitability through cosmic time. We find that, depending on the exact model assumptions, the Universe is 2.5 to 20 times more habitable today compared to when life first appeared on the Earth 4 billion years ago. We find that this increase in habitability will persist until the final stars die out over the next hundreds of billions of years.


title: Quintuple: a Python 5-qubit quantum computer simulator to facilitate cloud quantum computing

abstract: In May 2016 IBM released access to its 5-qubit quantum computer to the scientific community, its "IBM Quantum Experience" since acquiring over 25,000 users from students, educators and researchers around the globe. In the short time since the "IBM Quantum Experience" became available, a flurry of research results on 5-qubit systems have been published derived from the platform hardware. Quintuple is an open-source object-oriented Python module implementing the simulation of the "IBM Quantum Experience" hardware. Quintuple quantum algorithms can be programmed and run via a custom language fully compatible with the "IBM Quantum Experience" or in pure Python. Over 40 example programs are provided with expected results, including Grover's Algorithm and the Deutsch-Jozsa algorithm. Quintuple contributes to the study of 5-qubit systems and the development and debugging of quantum algorithms for deployment on the "IBM Quantum Experience" hardware.


title: All-strain based valley filter in graphene nanoribbons using snake states

abstract: A pseudo-magnetic field kink can be realized along a graphene nanoribbon using strain engineering. Electron transport along this kink is governed by snake states that are characterized by a single propagation direction. Those pseudo-magnetic fields point towards opposite directions in the K and K' valleys, leading to valley polarized snake states. In a graphene nanoribbon with armchair edges this effect results in a valley filter that is based only on strain engineering. We discuss how to maximize this valley filtering by adjusting the parameters that define the stress distribution along the graphene ribbon.


title: Density scaling relation in Orion A: effects of region selection

abstract: Recently Stanchev et al. (2015) proposed a technique to derive density scaling relations in a star-forming region from analysis of the probability distribution function of column density. We address the possible dependence of the outcome on the selection of probe zones, applying the method to Planck dust-opacity data on Orion A. The derived steep scaling relation of mean density with index -1.6 in the molecular cloud (so called `Central filament') points to its self-gravitating nature. The result is reproduced also for large parts of the clouds' vicinity which indicates major role of gravity in the energy balance of the entire star-forming region.


title: How flat is flat in random interface growth?

abstract: Domains of attraction are identified for the universality classes of one-point asymptotic fluctuations for the Kardar-Parisi-Zhang (KPZ) equation with general initial data. The criterion is based on a large deviation rate function for the rescaled initial data, which arises naturally from the Hopf-Cole transformation. This allows us, in particular, to distinguish the domains of attraction of curved, flat, and Brownian initial data, and to identify the boundary between the curved and flat domains of attraction, which turns out to correspond to square root initial data. The distribution of the asymptotic one-point fluctuations is characterized by means of a variational formula written in terms of certain limiting processes (arising as subsequential limits of the spatial fluctuations of KPZ equation with narrow wedge initial data, as shown in [CH16]) which are widely believed to coincide with the Airy$_2$ process. In order to identify these distributions for general initial data, we extend earlier results on continuum statistics of the Airy$_2$ process to probabilities involving the process on the entire line. In particular, this allows us to write an explicit Fredholm determinant formula for the case of square root initial data.


title: Shock waves in dispersive hydrodynamics with non-convex dispersion

abstract: Dissipationless hydrodynamics regularized by dispersion describe a number of physical media including water waves, nonlinear optics, and Bose-Einstein condensates. As in the classical theory of hyperbolic equations where a non-convex flux leads to non-classical solution structures, a non-convex linear dispersion relation provides an intriguing dispersive hydrodynamic analogue. Here, the fifth order Korteweg-de Vries (KdV) equation, also known as the Kawahara equation, a classical model for shallow water waves, is shown to be a universal model of Eulerian hydrodynamics with higher order dispersive effects. Utilizing asymptotic methods and numerical computations, this work classifies the long-time behavior of solutions for step-like initial data. For convex dispersion, the result is a dispersive shock wave (DSW), qualitatively and quantitatively bearing close resemblance to the KdV DSW. For non-convex dispersion, three distinct dynamic regimes are observed. For small jumps, a perturbed KdV DSW with positive polarity and orientation is generated, accompanied by small amplitude radiation from an embedded solitary wave leading edge, termed a radiating DSW or RDSW. For moderate jumps, a crossover regime is observed with waves propagating forward and backward from the sharp transition region. For jumps exceeding a critical threshold, a new type of DSW is observed we term a translating DSW or TDSW. The TDSW consists of a traveling wave that connects a partial, non-monotonic, negative solitary wave at the trailing edge to an interior nonlinear periodic wave. Its speed, a generalized Rankine-Hugoniot jump condition, is determined by the far-field structure of the traveling wave. The TDSW is resolved at the leading edge by a harmonic wavepacket moving with the linear group velocity. The non-classical TDSW exhibits features common to both dissipative and dispersive shock waves.


title: Feedback stabilization of the Cahn-Hilliard type system for phase separation

abstract: This article is concerned with the internal feedback stabilization of the phase field system of Cahn-Hilliard type, modeling the phase separation in a binary mixture. Under suitable assumptions on an arbitrarily fixed stationary solution, we construct via spectral separation arguments a feedback controller having the support in an arbitrary open subset of the space domain, such that the closed loop nonlinear system exponentially reach the prescribed stationary solution. This feedback controller has a finite dimensional structure in the state space of solutions. In particular, every constant stationary solution is admissible.


title: Aharonov-Bohm effect for a fermion field in a planar black hole "spacetime"

abstract: In this paper we consider the dynamics of a massive spinor field in the background of the acoustic black hole spacetime. Although this effective metric is acoustic and describes the propagation of sound waves, it can be considered as a toy model for the gravitational black hole. In this manner, we study the properties of the dynamics of the fermion field in this "gravitational" rotating black hole as well as the vortex background. We compute the differential cross section through the use of the partial wave approach and show that an effect similar to the gravitational Aharonov-Bohm effect occurs for the massive fermion field moving in this effective metric. We discuss the limiting cases and compare the results with the massless scalar field case.


title: Neutron Star Mass-Radius Constraints using Evolutionary Optimization

abstract: The equation of state of cold supra-nuclear-density matter, such as in neutron stars, is an open question in astrophysics. A promising method for constraining the neutron star equation of state is modelling pulse profiles of thermonuclear X-ray burst oscillations from hotspots on accreting neutron stars. The pulse profiles, constructed using spherical and oblate neutron star models, are comparable to what would be observed by a next-generation X-ray timing instrument like ASTROSAT, NICER, or LOFT. In this paper, we showcase the use of an evolutionary optimization algorithm to fit pulse profiles to determine the best-fit masses and radii. By fitting synthetic data, we assess how well the optimization algorithm can recover the input parameters. Multiple Poisson realizations of the synthetic pulse profiles, constructed with 1.6 million counts and no background, were fitted with the Ferret algorithm to analyze both statistical and degeneracy-related uncertainty, and to explore how the goodness-of-fit depends on the input parameters. For the regions of parameter space sampled by our tests, the best-determined parameter is the projected velocity of the spot along the observer's line-of-sight, with an accuracy of $\le3$% compared to the true value and with $\le5$\% statistical uncertainty. The next best-determined are the mass and radius; for a neutron star with a spin frequency of 600\,Hz, the best-fit mass and radius are accurate to $\le5$%, with respective uncertainties of $\le7$% and $\le10$%. The accuracy and precision depend on the observer inclination and spot co-latitude, with values of $\sim1$ % achievable in mass and radius if both the inclination and co-latitude are > 60 degrees.


title: Unequal Error Protection Querying Policies for the Noisy 20 Questions Problem

abstract: In this paper, we propose an open-loop unequal-error-protection querying policy based on superposition coding for the noisy 20 questions problem. In this problem, a player wishes to successively refine an estimate of the value of a continuous random variable by posing binary queries and receiving noisy responses. When the queries are designed non-adaptively as a single block and the noisy responses are modeled as the output of a binary symmetric channel the 20 questions problem can be mapped to an equivalent problem of channel coding with unequal error protection (UEP). A new non-adaptive querying strategy based on UEP superposition coding is introduced whose estimation error decreases with an exponential rate of convergence that is significantly better than that of the UEP repetition coding introduced by Variani et al. (2015). With the proposed querying strategy, the rate of exponential decrease in the number of queries matches the rate of a closed-loop adaptive scheme where queries are sequentially designed with the benefit of feedback. Furthermore, the achievable error exponent is significantly better than that of random block codes employing equal error protection.


title: Intermodulation in Nonlinear SQUID Metamaterials: Experiment and Theory

abstract: The response of nonlinear metamaterials and superconducting electronics to two-tone excitation is critical for understanding their use as low-noise amplifiers and tunable filters. A new setting for such studies is that of metamaterials made of radio frequency superconducting quantum interference devices (rf-SQUIDs). The two-tone response of self-resonant rf-SQUID meta-atoms and metamaterials is studied here via intermodulation (IM) measurement over a broad range of tone frequencies and tone powers. A sharp onset followed by a surprising strongly suppressed IM region near the resonance is observed. Using a two time scale analysis technique, we present an analytical theory that successfully explains our experimental observations. The theory predicts that the IM can be manipulated with tone power, center frequency, frequency difference between the two tones, and temperature. This quantitative understanding potentially allows for the design of rf-SQUID metamaterials with either very low or very high IM response.


title: Graphone (one-side hydrogenated graphene) formation on different substrates

abstract: In this work we present a fully atomistic reactive (ReaxFF force field) molecular dynamics study of the structural and dynamical aspects of the one-side hydrogenation of graphene membranes, leading to the formation of the so-called graphone structure. We have considered different substrates: graphene, few-layers graphene, graphite and platinum at different temperatures. Our results showed that the hydrogenation rates are very dependent on the substrate and thermal effects. Our results also showed that, similarly to graphane, large hydrogenated domains are unlikely to be formed. These hydrogenation processes occur through the formation of uncorrelated cluster domains.


title: The Future of Computing Research: Industry-Academic Collaborations

abstract: IT-driven innovation is an enormous factor in the worldwide economic leadership of the United States. It is larger than finance, construction, or transportation, and it employs nearly 6% of the US workforce. The top three companies, as measured by market capitalization, are IT companies - Apple, Google (now Alphabet), and Microsoft. Facebook, a relatively recent entry in the top 10 list by market capitalization has surpassed Walmart, the nation's largest retailer, and the largest employer in the world. The net income of just the top three exceeds $80 billion - roughly 100 times the total budget of the NSF CISE directorate which funds 87% of computing research. In short, the direct return on federal research investments in IT research has been enormously profitable to the nation. The IT industry ecosystem is also evolving. The time from conception to market of successful products has been cut from years to months. Product life cycles are increasingly a year or less. This change has pressured companies to focus industrial R&D on a pipeline or portfolio of technologies that bring immediate, or almost immediate, value to the companies. To defeat the competition and stay ahead of the pack, a company must devote resources to realizing gains that are shorter term, and must remain agile to respond quickly to market changes driven by new technologies, new startups, evolving user experience expectations, and the continuous consumer demand for new and exciting products. Amidst this landscape, the Computing Community Consortium convened a round-table of industry and academic participants to better understand the landscape of industry-academic interaction, and to discuss possible actions that might be taken to enhance those interactions. We close with some recommendations for actions that could expand the lively conversation we experienced at the round-table to a national scale.


title: Kaehler structures on spin 6-manifolds

abstract: We show that many spin 6-manifolds have the homotopy type but not the homeomorphism type of a Kaehler manifold. Moreover, for given Betti numbers, there are only finitely many deformation types and hence topological types of smooth complex projectve spin threefolds of general type. Finally, on a fixed spin 6-manifold, the Chern numbers take on only finitely many values on all possible Kaehler structures.


title: A Chebotarev variant of the Brun-Titchmarsh theorem and bounds for the Lang-Trotter conjectures

abstract: We improve the Chebotarev variant of the Brun-Titchmarsh theorem proven by Lagarias, Montgomery, and Odlyzko using the log-free zero density estimate and zero repulsion phenomenon for Hecke L-functions that were recently proved by the authors. Our result produces an improvement for the best unconditional bounds toward two conjectures of Lang and Trotter regarding the distribution of traces of Frobenius for elliptic curves and holomorphic cuspidal modular forms. We also obtain new results on the distribution of primes represented by positive-definite integral binary quadratic forms.


title: Learning Concept Taxonomies from Multi-modal Data

abstract: We study the problem of automatically building hypernym taxonomies from textual and visual data. Previous works in taxonomy induction generally ignore the increasingly prominent visual data, which encode important perceptual semantics. Instead, we propose a probabilistic model for taxonomy induction by jointly leveraging text and images. To avoid hand-crafted feature engineering, we design end-to-end features based on distributed representations of images and words. The model is discriminatively trained given a small set of existing ontologies and is capable of building full taxonomies from scratch for a collection of unseen conceptual label items with associated images. We evaluate our model and features on the WordNet hierarchies, where our system outperforms previous approaches by a large gap.


title: Abelian $n$-division fields of elliptic curves and Brauer groups of product Kummer & abelian surfaces

abstract: Let $Y$ be a principal homogeneous space of an abelian surface, or a K3 surface, over a finitely generated extension of $\mathbb{Q}$. In 2008, Skorobogatov and Zarhin showed that the Brauer group modulo algebraic classes $\text{Br}\, Y/ \text{Br}_1\, Y$ is finite. We study this quotient for the family of surfaces that are geometrically isomorphic to a product of isogenous non-CM elliptic curves, as well as the related family of geometrically Kummer surfaces; both families can be characterized by their geometric N\'eron-Severi lattices. Over a field of characteristic $0$, we prove that the existence of a strong uniform bound on the size of the odd-torsion of $\text{Br}\, Y / \text{Br}_1\, Y$ is equivalent to the existence of a strong uniform bound on integers $n$ for which there exist non-CM elliptic curves with abelian $n$-division fields. Using the same methods we show that, for a fixed prime $p$, a number field $k$ of fixed degree $r$, and a fixed discriminant of the geometric N\'eron-Severi lattice, $(\text{Br}\, Y / \text{Br}_1\, Y)[p^\infty]$ is bounded by a constant that depends only on $p$, $r$, and the discriminant.


title: Quantum superpositions of "common-cause" and "direct-cause" causal structures

abstract: The constraints arising for a general set of causal relations, both classically and quantumly, are still poorly understood. As a step in exploring this question, we consider a coherently controlled superposition of "direct-cause" and "common-cause" relationships between two events. We propose an implementation involving the spatial superposition of a mass and general relativistic time dilation. Finally, we develop a computationally efficient method to distinguish such genuinely quantum causal structures from classical (incoherent) mixtures of causal structures and show how to design experimental verifications of the nonclassicality of a causal structure.


title: Swift: Compiled Inference for Probabilistic Programming Languages

abstract: A probabilistic program defines a probability measure over its semantic structures. One common goal of probabilistic programming languages (PPLs) is to compute posterior probabilities for arbitrary models and queries, given observed evidence, using a generic inference engine. Most PPL inference engines---even the compiled ones---incur significant runtime interpretation overhead, especially for contingent and open-universe models. This paper describes Swift, a compiler for the BLOG PPL. Swift-generated code incorporates optimizations that eliminate interpretation overhead, maintain dynamic dependencies efficiently, and handle memory management for possible worlds of varying sizes. Experiments comparing Swift with other PPL engines on a variety of inference problems demonstrate speedups ranging from 12x to 326x.


title: Supernova neutrino physics with xenon dark matter detectors: A timely perspective

abstract: Dark matter detectors that utilize liquid xenon have now achieved tonne-scale targets, giving them sensitivity to all flavours of supernova neutrinos via coherent elastic neutrino-nucleus scattering. Considering for the first time a realistic detector model, we simulate the expected supernova neutrino signal for different progenitor masses and nuclear equations of state in existing and upcoming dual-phase liquid xenon experiments. We show that the proportional scintillation signal (S2) of a dual-phase detector allows for a clear observation of the neutrino signal and guarantees a particularly low energy threshold, while the backgrounds are rendered negligible during the supernova burst. XENON1T (XENONnT and LZ; DARWIN) experiments will be sensitive to a supernova burst up to 25 (35; 65) kpc from Earth at a significance of more than 5 sigma, observing approximately 35 (123; 704) events from a 27 Msun supernova progenitor at 10 kpc. Moreover, it will be possible to measure the average neutrino energy of all flavours, to constrain the total explosion energy, and to reconstruct the supernova neutrino light curve. Our results suggest that a large xenon detector such as DARWIN will be competitive with dedicated neutrino telescopes, while providing complementary information that is not otherwise accessible.


title: A Conic Section Problem Involving the Maximum Generalized Golden Right Triangle

abstract: An interesting conic section problem involving the maximum generalized golden right triangle $T_2$ is solved, and two simple constructions of $T_2$ are shown.


title: String Theory Origin of Constrained Multiplets

abstract: We study the non-linearly realized spontaneously broken supersymmetry of the (anti-)D3-brane action in type IIB string theory. The worldvolume fields are one vector $A_\mu$, three complex scalars $\phi^i$ and four 4d fermions $\lambda^0$, $\lambda^i$. These transform, in addition to the more familiar N=4 linear supersymmetry, also under 16 spontaneously broken, non-linearly realized supersymmetries. We argue that the worldvolume fields can be packaged into the following constrained 4d non-linear N=1 multiplets: four chiral multiplets $S$, $Y^i$ that satisfy $S^2=SY^i=0$ and contain the worldvolume fermions $\lambda^0$ and $\lambda^i$; and four chiral multiplets $W_\alpha$, $H^i$ that satisfy $S W_\alpha=0$ and $S \bar D_{\dot \alpha} \bar H^{\bar \imath}=0$ and contain the vector $A_\mu$ and the scalars $\phi^i$. We also discuss how placing an anti-D3-brane on top of intersecting O7-planes can lead to an orthogonal multiplet $\Phi$ that satisfies $S(\Phi-\bar \Phi)=0$, which is particularly interesting for inflationary cosmology.


title: No Timing Variations Observed in Third Transit of Snow-Line Exoplanet Kepler-421b

abstract: We observed Kepler-421 during the anticipated third transit of the snow-line exoplanet Kepler-421b in order to constrain the existence and extent of transit timing variations (TTVs). Previously, the Kepler Spacecraft only observed two transits of Kepler-421b leaving the planet's transit ephemeris unconstrained. Our visible light, time-series observations from the 4.3-meter Discovery Channel Telescope were designed to capture pre-transit baseline and the partial transit of Kepler-421b barring significant TTVs. We use the light curves to assess the probabilities of various transit models using both the posterior odds ratio and the Bayesian Information Criterion (BIC) and find that a transit model with no TTVs is favored to 3.6-sigma confidence. These observations suggest that Kepler-421b is either alone in its system or is only experiencing minor dynamic interactions with an unseen companion. With the Kepler-421b ephemeris constrained, we calculate future transit times and discuss the opportunity to characterize the atmosphere of this cold, long-period exoplanet via transmission spectroscopy. Our investigation emphasizes the difficulties associated with observing long-period exoplanet transits and the consequences that arise from failing to refine transit ephemerides.


title: The properties of hypervelocity stars and S-stars originating from an eccentric disc around a supermassive black hole

abstract: Hypervelocity stars (HVSs) that are observed in the Galactic halo, are believed to be accelerated to large velocities by a process of tidal disruption of binary stars passing close to a supermassive black hole (SMBH) which resides in the center of the Galaxy. It is, however, still unclear, where these relatively young stars were born and which dynamical process pushed them to nearly radial orbits around the SMBH. In this paper we investigate the possibility that the young binaries originated from a thin eccentric disc, similar to the one observed in the Galactic center nowadays. By means of direct N-body simulations, we follow the dynamical evolution of an initially thin and eccentric disc of stars with a 100% binary fraction orbiting around the SMBH. Such a configuration leads to Kozai-Lidov oscillations of orbital elements, bringing considerable amount of binaries to close vicinity of the black hole. Subsequent tidal disruption of these binaries accelerates one of their component to velocities well above the escape velocity from the SMBH while the second component becomes tightly bound to the SMBH. We describe the main kinematic properties of the escaping and tightly bound stars within our model and compare them qualitatively to the properties of the observed HVSs and S-stars, respectively. The most prominent feature is a strong anisotropy in the directions of the escaping stars which is observed for the Galactic HVSs but not explained yet.


title: SMPDF Web: a web-based application for specialized minimal parton distribution functions

abstract: We present SMPDF Web, a web interface for the construction of parton distribution functions (PDFs) with a minimal number of error sets needed to represent the PDF uncertainty of specific processes (SMPDF).


title: Damped Topological Magnons in the Kagom\'{e}-Lattice Ferromagnets

abstract: We demonstrate that interactions can substantially undermine the free-particle description of magnons in ferromagnets on geometrically frustrated lattices. The anharmonic coupling, facilitated by the Dzyaloshinskii-Moriya interaction, and a highly-degenerate two-magnon continuum yield a strong, non-perturbative damping of the high-energy magnon modes. We provide a detailed account of the effect for the $S=1/2$ ferromagnet on the kagom\'e lattice and propose further experiments.


title: The Gamma-Ray Pulsar Population of Globular Clusters: Implications for the GeV Excess

abstract: It has been suggested that the GeV excess, observed from the region surrounding the Galactic Center, might originate from a population of millisecond pulsars that formed in globular clusters. With this in mind, we employ the publicly available Fermi data to study the gamma-ray emission from 157 globular clusters, identifying a statistically significant signal from 25 of these sources (ten of which are not found in existing gamma-ray catalogs). We combine these observations with the predicted pulsar formation rate based on the stellar encounter rate of each globular cluster to constrain the gamma-ray luminosity function of millisecond pulsars in the Milky Way's globular cluster system. We find that this pulsar population exhibits a luminosity function that is quite similar to those millisecond pulsars observed in the field of the Milky Way (i.e. the thick disk). After pulsars are expelled from a globular cluster, however, they continue to lose rotational kinetic energy and become less luminous, causing their luminosity function to depart from the steady-state distribution. Using this luminosity function and a model for the globular cluster disruption rate, we show that millisecond pulsars born in globular clusters can account for only a few percent or less of the observed GeV excess. Among other challenges, scenarios in which the entire GeV excess is generated from such pulsars are in conflict with the observed mass of the Milky Way's Central Stellar Cluster.

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


title: SPARC: Mass Models for 175 Disk Galaxies with Spitzer Photometry and Accurate Rotation Curves

abstract: We introduce SPARC (Spitzer Photometry & Accurate Rotation Curves): a sample of 175 nearby galaxies with new surface photometry at 3.6 um and high-quality rotation curves from previous HI/Halpha studies. SPARC spans a broad range of morphologies (S0 to Irr), luminosities (~5 dex), and surface brightnesses (~4 dex). We derive [3.6] surface photometry and study structural relations of stellar and gas disks. We find that both the stellar mass-HI mass relation and the stellar radius-HI radius relation have significant intrinsic scatter, while the HI mass-radius relation is extremely tight. We build detailed mass models and quantify the ratio of baryonic-to-observed velocity (Vbar/Vobs) for different characteristic radii and values of the stellar mass-to-light ratio (M/L) at [3.6]. Assuming M/L=0.5 Msun/Lsun (as suggested by stellar population models) we find that (i) the gas fraction linearly correlates with total luminosity, (ii) the transition from star-dominated to gas-dominated galaxies roughly corresponds to the transition from spiral galaxies to dwarf irregulars in line with density wave theory; and (iii) Vbar/Vobs varies with luminosity and surface brightness: high-mass, high-surface-brightness galaxies are nearly maximal, while low-mass, low-surface-brightness galaxies are submaximal. These basic properties are lost for low values of M/L=0.2 Msun/Lsun as suggested by the DiskMass survey. The mean maximum-disk limit in bright galaxies is M/L=0.7 Msun/Lsun at [3.6]. The SPARC data are publicly available and represent an ideal test-bed for models of galaxy formation.


title: Metal flows of the circumgalactic medium, and the metal budget in galactic halos

abstract: We present an analysis of the flow of metals through the circumgalactic medium (CGM) in the Feedback in Realistic Environments (FIRE) simulations of galaxy formation, ranging from isolated dwarfs to $L*$ galaxies. We find that nearly all metals produced in high-redshift galaxies are carried out in winds that reach $0.25 R_{\rm vir}$. When measured at $0.25 R_{\rm vir}$ the metallicity of outflows is slightly higher than the interstellar medium (ISM) metallicity. Many metals thus reside in the CGM. Cooling and recycling from this reservoir determine the metal budget in the ISM. The outflowing metal flux decreases by a factor of $\sim2-5$ between $0.25 R_{\rm vir}$ and $R_{\rm vir}$. Furthermore, outflow metallicity is typically lower at $R_{\rm vir}$ owing to dilution of the remaining outflow by metal-poor material swept up from the CGM. The inflow metallicity at $R_{\rm vir}$ is generally low, but outflow and inflow metallicities are similar in the inner halo. At low redshift, massive galaxies no longer generate outflows that reach the CGM, causing a divergence in CGM and ISM metallicity. Dwarf galaxies continue to generate outflows, although they preferentially retain metal ejecta. In all but the least massive galaxy considered, a majority of the metals are within the halo at $z=0$. We measure the fraction of metals in CGM, ISM and stars, and quantify the thermal state of CGM metals in each halo. The total amount of metals in the low-redshift CGM of two simulated $L*$ galaxies is consistent with estimates from the COS halos survey, while for the other two it appears to be lower.


title: Constraining Light-Quark Yukawa Couplings from Higgs Distributions

abstract: We propose a novel strategy to constrain the bottom and charm Yukawa couplings by exploiting LHC measurements of transverse momentum distributions in Higgs production. Our method does not rely on the reconstruction of exclusive final states or heavy-flavour tagging. Compared to other proposals it leads to an enhanced sensitivity to the Yukawa couplings due to distortions of the differential Higgs spectra from emissions which either probe quark loops or are associated to quark-initiated production. We derive constraints using data from LHC Run I, and we explore the prospects of our method at future LHC runs. Finally, we comment on the possibility of bounding the strange Yukawa coupling.


title: The Herschel-ATLAS Data Release 1 Paper II: Multi-wavelength counterparts to submillimetre sources

abstract: This paper is the second in a pair of articles presenting data release 1 (DR1) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), the largest single open-time key project carried out with the Herschel Space Observatory. The H-ATLAS is a wide-area imaging survey carried out in five photometric bands at 100, 160, 250, 350 and 500$\mu$m covering a total area of 600deg$^2$. In this paper we describe the identification of optical counterparts to submillimetre sources in DR1, comprising an area of 161 deg$^2$ over three equatorial fields of roughly 12$^\circ$x4.5$^\circ$ centred at 9$^h$, 12$^h$ and 14.5$^h$ respectively. Of all the H-ATLAS fields, the equatorial regions benefit from the greatest overlap with current multi-wavelength surveys spanning ultraviolet (UV) to mid-infrared regimes, as well as extensive spectroscopic coverage. We use a likelihood-ratio technique to identify SDSS counterparts at r<22.4 for 250-$\mu$m-selected sources detected at $\geq$ 4$\sigma$ ($\approx$28mJy). We find `reliable' counterparts (reliability R$\geq$0.8) for 44,835 sources (39 per cent), with an estimated completeness of 73.0 per cent and contamination rate of 4.7 per cent. Using redshifts and multi-wavelength photometry from GAMA and other public catalogues, we show that H-ATLAS-selected galaxies at $z<0.5$ span a wide range of optical colours, total infrared (IR) luminosities, and IR/UV ratios, with no strong disposition towards mid-IR-classified AGN in comparison with optical selection. The data described herein, together with all maps and catalogues described in the companion paper (Valiante et al. 2016), are available from the H-ATLAS website at www.h-atlas.org.