title: A Single Photon Can Simultaneously Excite Two or More Atoms

abstract: We consider two separate atoms interacting with a single-mode optical resonator. When the frequency of the resonator field is twice the atomic transition frequency, we show that there exists a resonant coupling between \textit{one} photon and \textit{two} atoms, via intermediate virtual states connected by counter-rotating processes. If the resonator is prepared in its one-photon state, the photon can be jointly absorbed by the two atoms in their ground state which will both reach their excited state with probability close to one. Like ordinary quantum Rabi oscillations, this process is coherent and reversible, so that two atoms in their excited state will undergo a downward transition jointly emitting a single cavity photon. This joint absorption and emission processes can also occur with three atoms. The parameters used to investigate this process correspond to experimentally demonstrated values in circuit quantum electrodynamics systems.


title: Charge-dependent correlations from event-by-event anomalous hydrodynamics

abstract: We report on our recent attempt of quantitative modeling of the Chiral Magnetic Effect (CME) in heavy-ion collisions. We perform 3+1 dimensional anomalous hydrodynamic simulations on an event-by-event basis, with constitutive equations that contain the anomaly-induced effects. We also develop a model of the initial condition for the axial charge density that captures the statistical nature of random chirality imbalances created by the color flux tubes. Basing on the event-by-event hydrodynamic simulations for hundreds of thousands of collisions, we calculate the correlation functions that are measured in experiments, and discuss how the anomalous transport affects these observables.


title: Hyper innerproduct spaces II

abstract: In this paper, we extend the definition of hyperinner product defined on weak hypervector spaces with a hyperoperation scalar product to weak hypervector spaces with the hyperoperations sum and scalar products.


title: Scaling limits for sub-ballistic biased random walks in random conductances

abstract: We consider biased random walks in positive random conductances on the d-dimensional lattice in the zero-speed regime and study their scaling limits. We obtain a functional Law of Large Numbers for the position of the walker, properly rescaled. Moreover, we state a functional Central Limit Theorem where an atypical process, related to the Fractional Kinetics, appears in the limit.


title: Intrinsic Anomalous Hall Effect in Type-II Weyl Semimetals

abstract: Recently, a new type of Weyl semimetal called type-II Weyl semimetal has been proposed. Unlike the usual (type-I) Weyl semimetal, which has a point-like Fermi surface, this new type of Weyl semimetal has a tilted conical spectrum around the Weyl point. Here we calculate the anomalous Hall conductivity of a Weyl semimetal with a tilted conical spectrum for a pair of Weyl points, using the Kubo formula. We find that the Hall conductivity is not universal and can change sign as a function of the parameters quantifying the tilts. Our results suggest that even for the case where the separation between the Weyl points vanishes, tilting of the conical spectrum could give rise to a finite anomalous Hall effect, if the tilts of the two cones are not identical.


title: An expanded evaluation of protein function prediction methods shows an improvement in accuracy

abstract: Background: The increasing volume and variety of genotypic and phenotypic data is a major defining characteristic of modern biomedical sciences. At the same time, the limitations in technology for generating data and the inherently stochastic nature of biomolecular events have led to the discrepancy between the volume of data and the amount of knowledge gleaned from it. A major bottleneck in our ability to understand the molecular underpinnings of life is the assignment of function to biological macromolecules, especially proteins. While molecular experiments provide the most reliable annotation of proteins, their relatively low throughput and restricted purview have led to an increasing role for computational function prediction. However, accurately assessing methods for protein function prediction and tracking progress in the field remain challenging. Methodology: We have conducted the second Critical Assessment of Functional Annotation (CAFA), a timed challenge to assess computational methods that automatically assign protein function. One hundred twenty-six methods from 56 research groups were evaluated for their ability to predict biological functions using the Gene Ontology and gene-disease associations using the Human Phenotype Ontology on a set of 3,681 proteins from 18 species. CAFA2 featured significantly expanded analysis compared with CAFA1, with regards to data set size, variety, and assessment metrics. To review progress in the field, the analysis also compared the best methods participating in CAFA1 to those of CAFA2. Conclusions: The top performing methods in CAFA2 outperformed the best methods from CAFA1, demonstrating that computational function prediction is improving. This increased accuracy can be attributed to the combined effect of the growing number of experimental annotations and improved methods for function prediction.


title: Weyl-Einstein structures on K-contact manifolds

abstract: We show that a compact K-contact manifold $(M,g,\xi)$ has a closed Weyl-Einstein connection compatible with the conformal structure $[g]$ if and only if it is Sasaki-Einstein.


title: The Role of Context Types and Dimensionality in Learning Word Embeddings

abstract: We provide the first extensive evaluation of how using different types of context to learn skip-gram word embeddings affects performance on a wide range of intrinsic and extrinsic NLP tasks. Our results suggest that while intrinsic tasks tend to exhibit a clear preference to particular types of contexts and higher dimensionality, more careful tuning is required for finding the optimal settings for most of the extrinsic tasks that we considered. Furthermore, for these extrinsic tasks, we find that once the benefit from increasing the embedding dimensionality is mostly exhausted, simple concatenation of word embeddings, learned with different context types, can yield further performance gains. As an additional contribution, we propose a new variant of the skip-gram model that learns word embeddings from weighted contexts of substitute words.


title: Configurable memory systems for embedded many-core processors

abstract: The memory system of a modern embedded processor consumes a large fraction of total system energy. We explore a range of different configuration options and show that a reconfigurable design can make better use of the resources available to it than any fixed implementation, and provide large improvements in both performance and energy consumption. Reconfigurability becomes increasingly useful as resources become more constrained, so is particularly relevant in the embedded space. For an optimised architectural configuration, we show that a configurable cache system performs an average of 20% (maximum 70%) better than the best fixed implementation when two programs are competing for the same resources, and reduces cache miss rate by an average of 70% (maximum 90%). We then present a case study of AES encryption and decryption, and find that a custom memory configuration can almost double performance, with further benefits being achieved by specialising the task of each core when parallelising the program.


title: Anisotropic confinement effects in a two-dimensional plasma crystal

abstract: The spectral asymmetry of the wave energy distribution of dust particles during mode-coupling induced melting, observed for the first time in plasma crystals by Cou\"edel et al. [Phys. Rev. E 89, 053108 (2014)], is studied theoretically and by molecular-dynamics simulations. It is shown that an anisotropy of the well confining the microparticles selects the directions of preferred particle motion. The observed differences in intensity of waves of opposed directions is explained by a nonvanishing phonon flux. Anisotropic phonon scattering by defects and Umklapp scattering are proposed as possible reasons for the mean phonon flux.


title: Rogue waves for a system of coupled derivative nonlinear Schr\"odinger equations

abstract: Rogue waves (RWs) are unexpectedly strong excitations emerging from an otherwise tranquil background. The nonlinear Schr\"odinger equation (NLSE), a ubiquitous model with wide applications to fluid mechanics, optics and plasmas, exhibits RWs only in the regime of modulation instability (MI) of the background. For system of multiple waveguides, the governing coupled NLSEs can produce regimes of MI and RWs, even if each component has dispersion and cubic nonlinearity of opposite signs. A similar effect will be demonstrated for a system of coupled derivative NLSEs (DNLSEs), where the special feature is the nonlinear self-steepening of narrow pulses. More precisely, these additional regimes of MI and RWs for coupled DNLSEs will depend on the mismatch in group velocities between the components, as well as the parameters for cubic nonlinearity and self-steepening. RWs considered in this work differ from those of the NLSEs in terms of the amplification ratio and criteria of existence. Applications to optics and plasma physics are discussed.


title: Quantum subgroups of simple twisted quantum groups at roots of one

abstract: Let $G$ be a connected, simply connected simple complex algebraic group and let $\epsilon$ be a primitive $\ell$th root of unity with $\ell$ odd and coprime with $3$ if $G$ is of type $G_{2}$. We determine all Hopf algebra quotients of the twisted multiparameter quantum function algebra $\mathcal{O}_{\epsilon}^{\varphi}(G)$ introduced by Costantini and Varagnolo. This extends the results of Andruskiewitsch and the first author, where the untwisted case is treated.


title: Kinetic modelling of runaway electrons in dynamic scenarios

abstract: Improved understanding of runaway-electron formation and decay processes are of prime interest for the safe operation of large tokamaks, and the dynamics of the runaway electrons during dynamical scenarios such as disruptions are of particular concern. In this paper, we present kinetic modelling of scenarios with time-dependent plasma parameters; in particular, we investigate hot-tail runaway generation during a rapid drop in plasma temperature. With the goal of studying runaway-electron generation with a self-consistent electric-field evolution, we also discuss the implementation of a collision operator that conserves momentum and energy and demonstrate its properties. An operator for avalanche runaway-electron generation, which takes the energy dependence of the scattering cross section and the runaway distribution into account, is investigated. We show that the simplified avalanche model of Rosenbluth & Putvinskii [Nucl. Fusion 1997 37 1355] can give inaccurate results for the avalanche growth rate (either lower or higher) for many parameters, especially when the average runaway energy is modest, such as during the initial phase of the avalanche multiplication. The developments presented pave the way for improved modelling of runaway-electron dynamics during disruptions or other dynamic events.


title: Secret Key Generation with Limited Interaction

abstract: A basic two-terminal secret key generation model is considered, where the interactive communication rate between the terminals may be limited, and in particular may not be enough to achieve the maximum key rate. We first prove a multi-letter characterization of the key-communication rate region (where the number of auxiliary random variables depend on the number of rounds of the communication), and then provide an equivalent but simpler characterization in terms of concave envelopes in the case of unlimited number of rounds. Two extreme cases are given special attention. First, in the regime of very low communication rates, the \emph{key bits per interaction bit} (KBIB) is expressed with a new "symmetric strong data processing constant", which has a concave envelope characterization analogous to that of the conventional strong data processing constant. The symmetric strong data processing constant can be upper bounded by the supremum of the maximal correlation coefficient over a set of distributions, which allows us to determine the KBIB for binary symmetric sources, and conclude, in particular, that the interactive scheme is not more efficient than the one-way scheme at least in the low communication-rate regime. Second, a new characterization of the \emph{minimum interaction rate needed for achieving the maximum key rate} (MIMK) is given, and we resolve a conjecture by Tyagi regarding the MIMK for (possibly nonsymmetric) binary sources. We also propose a new conjecture for binary symmetric sources that the interactive scheme is not more efficient than the one-way scheme at any communication rate.


title: Too good to be true: when overwhelming evidence fails to convince

abstract: Is it possible for a large sequence of measurements or observations, which support a hypothesis, to counterintuitively decrease our confidence? Can unanimous support be too good to be true? The assumption of independence is often made in good faith, however rarely is consideration given to whether a systemic failure has occurred. Taking this into account can cause certainty in a hypothesis to decrease as the evidence for it becomes apparently stronger. We perform a probabilistic Bayesian analysis of this effect with examples based on (i) archaeological evidence, (ii) weighing of legal evidence, and (iii) cryptographic primality testing. We find that even with surprisingly low systemic failure rates high confidence is very difficult to achieve and in particular we find that certain analyses of cryptographically-important numerical tests are highly optimistic, underestimating their false-negative rate by as much as a factor of $2^{80}$.


title: Joint learning of ontology and semantic parser from text

abstract: Semantic parsing methods are used for capturing and representing semantic meaning of text. Meaning representation capturing all the concepts in the text may not always be available or may not be sufficiently complete. Ontologies provide a structured and reasoning-capable way to model the content of a collection of texts. In this work, we present a novel approach to joint learning of ontology and semantic parser from text. The method is based on semi-automatic induction of a context-free grammar from semantically annotated text. The grammar parses the text into semantic trees. Both, the grammar and the semantic trees are used to learn the ontology on several levels -- classes, instances, taxonomic and non-taxonomic relations. The approach was evaluated on the first sentences of Wikipedia pages describing people.


title: Some studies on quantum equivalents of non-commutative operators via commutating eigenvalue relation: PT-symmetry

abstract: We study quantum equivalents of non-commutative operators in quantum mechanics. Any matrix "$B$" satisfying the non-commuting relation $[A,B]\neq 0$ with "$A$", can be used via $B^{-1} AB$ to reproduce eigenvalues of "$A$". This universality relation is also equally valid for any matrix in any branch of physical or social science and also any operator involving co-ordinate$(x)$ or momentum$(p)$. Pictorially this is represented in fig. 1. Many interesting models including logarithmic potential have been considered.


title: Long memory and multifractality: A joint test

abstract: The properties of statistical tests for hypotheses concerning the parameters of the multifractal model of asset returns (MMAR) are investigated, using Monte Carlo techniques. We show that, in the presence of multifractality, conventional tests of long memory tend to over-reject the null hypothesis of no long memory. Our test addresses this issue by jointly estimating long memory and multifractality. The estimation and test procedures are applied to exchange rate data for 12 currencies. In 11 cases, the exchange rate returns are accurately described by compounding a NIID series with a multifractal time-deformation process. There is no evidence of long memory.


title: Transverse Energy Measurements from the Beam Energy Scan in PHENIX

abstract: Transverse energy distributions at midrapidity have been measured by the PHENIX experiment at the BNL Relativistic Heavy Ion Collider (RHIC) for Au+Au, U+U, Cu+Au, Cu+Cu, He+Au, d+Au, and p+p collisions over a wide energy range from $\sqrt{s_{NN}}$ = 7.7 GeV to $\sqrt{s_{NN}}$ = 200 GeV as a function of centrality. For central Au+Au collisions, it is observed that the midrapidity Bjorken energy density demonstrates a power law behavior from $\sqrt{s_{NN}}$ = 7.7 GeV to $\sqrt{s_{NN}}$ = 2.76 TeV. At a given collision energy, the data presented as a function of $N_{part}$ are independent of the size of the collision system. For Au+Au, Cu+Au, and Cu+Cu collisions, the centrality-dependent data are better described by scaling with the number of constituent quark participants than scaling with the number of nucleon participants.


title: Secular evolution of asteroid families: the role of Ceres

abstract: We consider the role of the dwarf planet Ceres on the secular dynamics of the asteroid main belt. Specifically, we examine the post impact evolution of asteroid families due to the interaction of their members with the linear nodal secular resonance with Ceres. First, we find the location of this resonance and identify which asteroid families are crossed by its path. Next, we summarize our results for three asteroid families, namely (1726) Hoffmeister, (1128) Astrid and (1521) Seinajoki which have irregular distributions of their members in the proper elements space, indicative of the effect of the resonance. We confirm this by performing a set of numerical simulations, showcasing that the perturbing action of Ceres through its linear nodal secular resonance is essential to reproduce the actual shape of the families.


title: CS emission near MIR-bubbles

abstract: We survey 44 young stellar objects located near the edges of mid-IR-identified bubbles in CS (1-0) using the Green Bank Telescope. We detect emission in 18 sources, indicating young protostars that are good candidates for being triggered by the expansion of the bubble. We calculate CS column densities and abundances. Three sources show evidence of infall through non-Gaussian line-shapes. Two of these sources are associated with dark clouds and are promising candidates for further exploration of potential triggered star formation. We obtained on-the-fly maps in CS (1-0) of three sources, showing evidence of significant interactions between the sources and the surrounding environment.


title: Dust production in debris discs: constraints on the smallest grains

abstract: The surface energy constraint puts a limit on the smallest fragment $s_{surf}$ that can be produced after a collision. Based on analytical considerations, this mechanism has been recently identified as been potentially able to prevent the production of small dust grains in debris discs and cut off their size distribution at sizes larger than the blow-out size. We numerically investigate the importance of this effect to find under which conditions it can leave a signature in the small-size end of a disc's particle size distribution (PSD). An important part of this work is to map out, in a disc at steady-state, what is the most likely collisional origin for micron-sized grains, in terms of the sizes of their collisional progenitors. We implement, for the first time, the surface energy constraint into a collisional evolution code. We consider a debris disc extending from 50 to 100AU and 2 different stellar types. We also consider two levels of stirring in the disc: dynamically "hot" (e=0.075) and "cold" (e=0.01). For all cases, we derive $s_{surf}$ maps as a function of target and projectile sizes, $s_t$ and $s_p$, and compare them to equivalent maps for the dust-production rate. We then compute disc-integrated PSDs and estimate the imprint of the surface energy constraint. We find that the ($s_p$,$s_t$) regions of high $s_{surf}$ values do not coincide with those of high dust production rate. As a consequence, the surface energy constraint has generally a weak effect on the system's PSD. The maximum $s_{surf}$-induced depletion of $\mu$m-sized grains is $\sim 30$% and is obtained for a sun-like star and a dynamically hot case. For the e=0.01 cases, the surface energy effect is negligible compared to the massive small grain depletion induced by another mechanism: the natural imbalance between dust production and destruction rates in low-stirring discs identified by Thebault&Wu(2008).


title: PVC Extrusion Development and Production for the NOvA Neutrino Experiment

abstract: We have produced large and highly-reflective open-cell PVC extrusions for the NOvA neutrino oscillation experiment. The extrusions were sealed, instrumented, assembled into self-supporting detector blocks, and filled with liquid scintillator. Each Far Detector block stands 15.7 m high, is 15.7 m wide and 2.1 m thick. More than 22,000 extrusions were produced with high dimensional tolerance and robust mechanical strength. This paper provides an overview of the NOvA Far Detector, describes the preparation of the custom PVC powder, and the making of the extrusions. Quality control was a key element in the production and is described in detail.


title: The high-conductance state enables neural sampling in networks of LIF neurons

abstract: The apparent stochasticity of in-vivo neural circuits has long been hypothesized to represent a signature of ongoing stochastic inference in the brain. More recently, a theoretical framework for neural sampling has been proposed, which explains how sample-based inference can be performed by networks of spiking neurons. One particular requirement of this approach is that the neural response function closely follows a logistic curve. Analytical approaches to calculating neural response functions have been the subject of many theoretical studies. In order to make the problem tractable, particular assumptions regarding the neural or synaptic parameters are usually made. However, biologically significant activity regimes exist which are not covered by these approaches: Under strong synaptic bombardment, as is often the case in cortex, the neuron is shifted into a high-conductance state (HCS) characterized by a small membrane time constant. In this regime, synaptic time constants and refractory periods dominate membrane dynamics. The core idea of our approach is to separately consider two different "modes" of spiking dynamics: burst spiking and transient quiescence, in which the neuron does not spike for longer periods. We treat the former by propagating the PDF of the effective membrane potential from spike to spike within a burst, while using a diffusion approximation for the latter. We find that our prediction of the neural response function closely matches simulation data. Moreover, in the HCS scenario, we show that the neural response function becomes symmetric and can be well approximated by a logistic function, thereby providing the correct dynamics in order to perform neural sampling. We hereby provide not only a normative framework for Bayesian inference in cortex, but also powerful applications of low-power, accelerated neuromorphic systems to relevant machine learning tasks.

categories: q-bio.NC cond-mat.dis-nn cs.NE physics.bio-ph stat.ML


title: New non-equilibrium matrix imbibition equation for Kondaurov's double porosity model

abstract: The paper deals with the global Kondaurov double porosity model describing a non-equilibrium two-phase immiscible flow in fractured-porous reservoirs when non-equilibrium phenomena occur in the matrix blocks, only. It is shown that the homogenized model can be represented as usual equations of two-phase incompressible immiscible flow, except for the addition of two source terms calculated by a solution to a local problem which is a boundary value problem for a non-equilibrium imbibition equation given in terms of the real saturation and a non-equilibrium parameter.


title: $Fermi$ GBM Observations of V404 Cyg During its 2015 Outburst

abstract: V404 Cygni was discovered in 1989 by the $Ginga$ X-ray satellite during its only previously observed X-ray outburst and soon after confirmed as a black hole binary. On June 15, 2015, the Gamma Ray Burst Monitor (GBM) triggered on a new outburst of V404 Cygni. We present 13 days of GBM observations of this outburst including Earth occultation flux measurements, spectral and temporal analysis. The Earth occultation fluxes reached 30 Crab with detected emission to 100 keV and determined, via hardness ratios, that the source was in a hard state. At high luminosity, spectral analysis between 8 and 300 keV showed that the electron temperature decreased with increasing luminosity. This is expected if the protons and electrons are in thermal equilibrium during an outburst with the electrons cooled by the Compton scattering of softer seed photons from the disk. However, the implied seed photon temperatures are unusually high, suggesting a contribution from another source, such as the jet. No evidence of state transitions is seen during this time period. The temporal analysis reveals power spectra that can be modeled with two or three strong, broad Lorentzians, similar to the power spectra of black hole binaries in their hard state.


title: Assessing Mission Impact of Cyberattacks: Report of the NATO IST-128 Workshop

abstract: This report presents the results of a workshop conducted by the North Atlantic Treaty Organization (NATO) Information Systems Technology (IST) Panel in Istanbul, Turkey, in June 2015 to explore science and technology for characterizing the impact of cyber-attacks on missions. Military mission success is highly dependent on the communications and information systems (CISs) that support the mission and their use in the cyber battlespace. The inexorably growing dependency on computational information processing for weapons, intelligence, communication, and logistics systems continues to increase the vulnerability of missions to various cyber threats. Attacks on CISs or other cyber incidents degrade or disrupt the usage of CISs, and the resulting mission capability, performance, and completion. These incidents are expected to increase in frequency and sophistication. The workshop participants concluded that the key to solving the mission impact assessment problem was in adopting and developing a new model-driven paradigm that creates and validates mechanisms of modeling the mission organization, the mission(s), and the cyber-vulnerable systems that support the mission(s). Such models then simulate or portray the impacts of the cyber-attacks. In addition, such model-based analysis could explore multiple alternative mitigation and work-around strategies - an essential part of coping with mission impact - and select the optimal course of mitigating actions. Only such a paradigm can be expected to provide meaningful, actionable information about mission impacts that have not been seen before or do not match prior experiences and patterns. The papers presented at this workshop are available in an accompanying volume, Proceedings of the NATO Workshop IST-128, Assessing Mission Impact of Cyber Attacks.


title: Reproducing sterile neutrinos and the behavior of flavor oscillations with superconducting-magnetic proximity effects

abstract: The physics of a superconductor subjected to a magnetic field is known to be equivalent to neutrino oscillations. Examining the properties of singlet-triplet oscillations in the magnetic field, a sterile neutrino--shown to be a Majorana fermion--is suggested to be represented by singlet Cooper pairs and moderates flavor oscillations between three flavor neutrinos (triplet Cooper pairs). A superconductor-exchange spring system's rotating magnetization profile is used to simulate the mass-flavor oscillations in the neutrino case and the physics of neutrino oscillations are discussed. Symmetry protected triplet components are presented as weak process states. Phases acquired due to the Fulde-Ferrell-Larkin-Ovchinnikov effect produce a complex phase that may be responsible for charge-parity violation in flavor oscillations.


title: Riemann Hypothesis and Random Walks: the Zeta case

abstract: In previous work it was shown that if certain series based on sums over primes of non-principal Dirichlet characters have a conjectured random walk behavior, then the Euler product formula for its $L$-function is valid to the right of the critical line $\Re (s) > \tfrac{1}{2}$, and the Riemann Hypothesis for this class of $L$-functions follows. Building on this work, here we propose how to extend this line of reasoning to the Riemann zeta function and other principal Dirichlet $L$-functions. We apply these results to the study of the argument of the zeta function. In another application, we define and study a 1-point correlation function of the Riemann zeros, which leads to the construction of a probabilistic model for them. Based on these results we describe a new algorithm for computing very high Riemann zeros, and we calculate the googol-th zero, namely $10^{100}$-th zero to over 100 digits, far beyond what is currently known.


title: Superlattice gain in positive differential conductivity region

abstract: We analyze theoretically a superlattice structure proposed by A. Andronov et al. [JETP Lett 102, 207 (2015)] to give Terahertz gain for an operation point with positive differential conductivity. Here we confirm the existence of gain and show that an optimized structure displays gain above 20 cm$^{-1}$ at low temperatures, so that lasing may be observable. Comparing a variety of simulations, this gain is found to be strongly affected by elastic scattering. It is shown that the dephasing modifies the nature of the relevant states, so that the common analysis based on Wannier-Stark states is not reliable for a quantitative description of the gain in structures with extremely diagonal transitions.


title: The annihilating-submodule graph of modules over commutative rings

abstract: Let M be a module over a commutative ring R. In this paper, we continue our study of annihilating-submodule graph AG(M) which was introduced in (The Zariski topology-graph of modules over commutative rings, Comm. Algebra., 42 (2014), 3283{3296). AG(M) is a (undirected) graph in which a nonzero submodule N of M is a vertex if and only if there exists a nonzero proper submodule K of M such that NK = (0), where NK, the product of N and K, is defined by (N : M)(K : M)M and two distinct vertices N and K are adjacent if and only if NK = (0). We obtain useful characterizations for those modules M for which either AG(M) is a complete (or star) graph or every vertex of AG(M) is a prime (or maximal) submodule of M. Moreover, we study coloring of annihilating-submodule graphs.


title: DrMAD: Distilling Reverse-Mode Automatic Differentiation for Optimizing Hyperparameters of Deep Neural Networks

abstract: The performance of deep neural networks is well-known to be sensitive to the setting of their hyperparameters. Recent advances in reverse-mode automatic differentiation allow for optimizing hyperparameters with gradients. The standard way of computing these gradients involves a forward and backward pass of computations. However, the backward pass usually needs to consume unaffordable memory to store all the intermediate variables to exactly reverse the forward training procedure. In this work we propose a simple but effective method, DrMAD, to distill the knowledge of the forward pass into a shortcut path, through which we approximately reverse the training trajectory. Experiments on several image benchmark datasets show that DrMAD is at least 45 times faster and consumes 100 times less memory compared to state-of-the-art methods for optimizing hyperparameters with minimal compromise to its effectiveness. To the best of our knowledge, DrMAD is the first research attempt to make it practical to automatically tune thousands of hyperparameters of deep neural networks. The code can be downloaded from https://github.com/bigaidream-projects/drmad


title: Simulation of residual oil displacement in a sinusoidal channel with the lattice Boltzmann method

abstract: We simulate oil slug displacement in a sinusoidal channel in order to validate computational models and algorithms for multi-component flow. This case fits in the gap between fully realistic cases characterized by complicated geometry and academic cases with simplistic geometry. Our computational model is based on the lattice Boltzmann method and allows for variation of physical parameters such as wettability and viscosity. The effect of variation of model parameters is analyzed, in particular via comparison with analytical solutions. We discuss the requirements for accurate solution of the oil slug displacement problem.


title: Exponential integrability properties of Euler discretization schemes for the Cox-Ingersoll-Ross process

abstract: We analyze exponential integrability properties of the Cox-Ingersoll-Ross (CIR) process and its Euler discretizations with various types of truncation and reflection at 0. These properties play a key role in establishing the finiteness of moments and the strong convergence of numerical approximations for a class of stochastic differential equations arising in finance. We prove that both implicit and explicit Euler-Maruyama discretizations for the CIR process preserve the exponential integrability of the exact solution for a wide range of parameters, and find lower bounds on the explosion time.


title: The homology core and invariant measures

abstract: Here we shall consider the topology and dynamics associated to a wide class of matchbox manifolds, including a large selection of tiling spaces and all minimal matchbox manifolds of dimension one. For such spaces we introduce topological invariants related to their expansions as an inverse sequence of simplicial complexes. These invariants are related to corresponding inverse sequences of groups arising from applying the top--dimension homology to these sequences. In many cases this leads to a computable invariant based on an inverse sequence of matrices. Significantly, we show that when the space is obtained by suspending a topologically transitive action of the fundamental group $\G$ of a closed orientable on a zero--dimensional compact space this invariant at the same time corresponds to the space of Borel measures on the Cantor set which are invariant under the action of $\G$. This leads to connections between the rank of homology groups we consider and the number of invariant, ergodic Borel probability measures for such actions. We illustrate with several examples how these invariants can be calculated and used for topological classification and how it leads to an understanding of the invariant measures.


title: Soft Hair on Black Holes

abstract: It has recently been shown that BMS supertranslation symmetries imply an infinite number of conservation laws for all gravitational theories in asymptotically Minkowskian spacetimes. These laws require black holes to carry a large amount of soft ($i.e.$ zero-energy) supertranslation hair. The presence of a Maxwell field similarly implies soft electric hair. This paper gives an explicit description of soft hair in terms of soft gravitons or photons on the black hole horizon, and shows that complete information about their quantum state is stored on a holographic plate at the future boundary of the horizon. Charge conservation is used to give an infinite number of exact relations between the evaporation products of black holes which have different soft hair but are otherwise identical. It is further argued that soft hair which is spatially localized to much less than a Planck length cannot be excited in a physically realizable process, giving an effective number of soft degrees of freedom proportional to the horizon area in Planck units.


title: Exact diagonalization of Heisenberg $SU(N)$ chains in the fully symmetric and antisymmetric representations

abstract: Motivated by recent experimental progress in the context of ultra-cold multi-color fermionic atoms in optical lattices, we have developed a method to exactly diagonalize the Heisenberg $SU(N)$ Hamiltonian with several particles per site living in a fully symmetric or antisymmetric representation of $SU(N)$. The method, based on the use of standard Young tableaux, takes advantage of the full $SU(N)$ symmetry, allowing one to work directly in each irreducible representations of the global $SU(N)$ group. Since the $SU(N)$ singlet sector is often much smaller than the full Hilbert space, this enables one to reach much larger system sizes than with conventional exact diagonalizations. The method is applied to the study of Heisenberg chains in the symmetric representation with two and three particles per site up to $N=10$ and up to 20 sites. For the length scales accessible to this approach, all systems except the Haldane chain ($SU(2)$ with two particles per site) appear to be gapless, and the central charge and scaling dimensions extracted from the results are consistent with a critical behaviour in the $SU(N)$ level $k$ Wess-Zumino-Witten universality class, where $k$ is the number of particles per site. These results point to the existence of a cross-over between this universality class and the asymptotic low-energy behavior with a gapped spectrum or a critical behavior in the $SU(N)$ level $1$ WZW universality class.


title: Observations and Analysis of the GK Persei Nova Shell and its "Jet-like" Feature

abstract: GK Persei (1901, the 'Firework Nebula') is an old but bright nova remnant that offers a chance to probe the physics and kinematics of nova shells. The kinematics in new and archival longslit optical echelle spectra were analysed using the shape software. New imaging from the Aristarchos telescope continues to track the proper motion, extinction and structural evolution of the knots, which have been observed intermittently over several decades. We present for the first time, kinematical constraints on a large faint 'jet' feature, that was previously detected beyond the shell boundary. These observational constraints allow for the generation of models for individual knots, interactions within knot complexes, and the 'jet' feature. Put together, and taking into account dwarf-nova accelerated winds emanating from the central source, these data and models give a deeper insight into the GK Per nova remnant as a whole.


title: Existence and stability of Lagrangian points in the relativistic restricted three body problem

abstract: It is well known that objects can oscillate around the Lagrangian point L4. In this manuscript we compute the period of these oscillations by computing the exact expression of the characteristic polynomial of the matrix that determined the stability of the Lagrangian point. When we use relativity theory we are supposed to get a small improvement of this period. This paper compares the values of the the characteristic polynomials computed in the following way: 1. Using relativity with an error smaller than 10^{-30}. 2. The characteristic polynomial coming from a paper written in 2002, 3. The characteristic polynomial coming from a paper written in 2006 and 4. The characteristic polynomial obtain without using relativity. We show that the polynomial 1 is closer to the polynomial 4 than it is to the polynomials 2 and 3. This is, the error in both polynomials that dealt with relativity coming from previous papers is bigger than the error coming from the characteristic polynomial that do not use relativity at all. We also prove the existence of the Lagrangian point L4 in the Sun-Earth system using the Poincare Miranda theorem.


title: Complex Decomposition of the Negative Distance kernel

abstract: A Support Vector Machine (SVM) has become a very popular machine learning method for text classification. One reason for this relates to the range of existing kernels which allow for classifying data that is not linearly separable. The linear, polynomial and RBF (Gaussian Radial Basis Function) kernel are commonly used and serve as a basis of comparison in our study. We show how to derive the primal form of the quadratic Power Kernel (PK) -- also called the Negative Euclidean Distance Kernel (NDK) -- by means of complex numbers. We exemplify the NDK in the framework of text categorization using the Dewey Document Classification (DDC) as the target scheme. Our evaluation shows that the power kernel produces F-scores that are comparable to the reference kernels, but is -- except for the linear kernel -- faster to compute. Finally, we show how to extend the NDK-approach by including the Mahalanobis distance.


title: The part-frequency matrices of a partition

abstract: A new combinatorial object is introduced, the part-frequency matrix sequence of a partition, which is elementary to describe and is naturally motivated by Glaisher's bijection. We prove results that suggest surprising usefulness for such a simple tool, including the existence of a related statistic that realizes every possible Ramanujan-type congruence for the partition function. To further exhibit its research utility, we give an easy generalization of a theorem of Andrews, Dixit and Yee on the mock theta functions. Throughout, we state a number of observations and questions that can motivate an array of investigations.


title: Degeneracies in long-baseline neutrino experiments from nonstandard interactions

abstract: We study parameter degeneracies that can occur in long-baseline neutrino appearance experiments due to nonstandard interactions (NSI) in neutrino propagation. For a single off-diagonal NSI parameter, and neutrino and antineutrino measurements at a single L/E, there exists a continuous four-fold degeneracy (related to the mass hierarchy and $\theta_{23}$ octant) that renders the mass hierarchy, octant, and CP phase unknowable. Even with a combination of NO$\nu$A and T2K data, which in principle can resolve the degeneracy, both NSI and the CP phase remain unconstrained because of experimental uncertainties. A wide-band beam experiment like DUNE will resolve this degeneracy if the nonzero off-diagonal NSI parameter is $\epsilon_{e\mu}$. If $\epsilon_{e\tau}$ is nonzero, or the diagonal NSI parameter $\epsilon_{ee}$ is O(1), a wrong determination of the mass hierarchy and of CP violation can occur at DUNE. The octant degeneracy can be further complicated by $\epsilon_{e\tau}$, but is not affected by $\epsilon_{ee}$.


title: A greedy algorithm for $B_h[g]$ sequences

abstract: For any positive integers $h\ge 2$ and $g\ge 1$, we present a greedy algorithm that provides an infinite $B_h[g]$ sequence with $a_n\le 2gn^{h+(h-1)/g}.$


title: The Mumford-Tate conjecture for the product of an abelian surface and a K3 surface

abstract: In this paper we prove the Mumford-Tate conjecture in degree 2 for the product of an abelian surface $A$ and a K3 surface $X$ over a finitely generated field $K \subset \mathbb{C}$. The Mumford-Tate conjecture is a precise way of saying that the Hodge structure on singular cohomology conveys the same information as the Galois representation on $\ell$-adic \'{e}tale cohomology. To make this precise, let $G_{\mathrm{B}}$ be the Mumford-Tate group of the Hodge structure $H^{2}_{\text{sing}}(A(\mathbb{C}) \times X(\mathbb{C}), \mathbb{Q})$. Let $G_{\ell}^{\circ}$ be the connected component of the identity of the Zariski closure of the image of the Galois group $\textrm{Gal}(\bar{K}/K)$ in $\mathrm{GL}(H^{2}_{\text{\'{e}t}}(A_{\bar{K}} \times X_{\bar{K}}, \mathbb{Q}_{\ell}))$. The Mumford-Tate conjecture asserts that $G_{\mathrm{B}} \otimes \mathbb{Q}_{\ell} \cong G_{\ell}^{\circ}$. The proof presented in this paper uses input from number theory (Chebotaryov's density theorem), Lie theory, and some facts about K3 surfaces over finite fields.


title: Cohomology of finite modules over short Gorenstein rings

abstract: Let $R$ be a Gorenstein local ring with maximal ideal $\mathfrak{m}$ satisfying $\mathfrak{m}^3=0\ne\mathfrak{m}^2$. Set $k=R/\mathfrak{m}$ and $e=\text{rank}_{k}(\mathfrak{m}/\mathfrak{m}^2)$. If $e>2$ and $M$, $N$ are finitely generated $R$-modules, we show that the formal power series $\sum_{i=0}^\infty\text{rank}_{k}\left(\text{Ext}^i_R(M,N)\otimes_R k \right)t^i$ and $\sum_{i=0}^\infty\text{rank}_{k}\left(\text{Tor}_i^R(M,N)\otimes_R k \right)t^i$ are rational, with denominator $1-et+t^2$.


title: A direct Eulerian GRP scheme for spherically symmetric general relativistic hydrodynamics

abstract: The paper proposes a second-order accurate direct Eulerian generalized Riemann problem (GRP) scheme for the spherically symmetric general relativistic hydrodynamical (RHD) equations and a second-order accurate discretization for the spherically symmetric Einstein (SSE) equations. The former is directly using the Riemann invariants and the Runkine-Hugoniot jump conditions to analytically resolve the left and right nonlinear waves of the local GRP in the Eulerian formulation together with the local change of the metrics to obtain the limiting values of the time derivatives of the conservative variables along the cell interface and the numerical flux for the GRP scheme. While the latter utilizes the energy-momentum tensor obtained in the GRP solver to evaluate the fluid variables in the SSE equations and keeps the continuity of the metrics at the cell interfaces. Several numerical experiments show that the GRP scheme can achieve second-order accuracy and high resolution, and is effective for spherically symmetric general RHD problems.


title: Fourier integrals operators on lie groupoids

abstract: As announced in [12], we develop a calculus of Fourier integral G-operators on any Lie groupoid G. For that purpose, we study convolability and invertibility of Lagrangian conic submanifolds of the symplectic groupoid T * G. We also identify those Lagrangian which correspond to equivariant families parametrized by the unit space G (0) of homogeneous canonical relations in (T * Gx \ 0) x (T * G x \ 0). This allows us to select a subclass of Lagrangian distributions on any Lie groupoid G that deserve the name of Fourier integral G-operators (G-FIO). By construction, the class of G-FIO contains the class of equivariant families of ordinary Fourier integral operators on the manifolds Gx, x $\in$ G (0). We then develop for G-FIO the first stages of the calculus in the spirit of Hormander's work. Finally, we work out an example proving the efficiency of the present approach for studying Fourier integral operators on singular manifolds.


title: Density functional theory study of skyrmion pinning by atomic defects in MnSi

abstract: A magnetic skyrmion observed experimentally in chiral magnets is a topologically protected spin texture. For their unique properties, such as high mobility under current drive, skyrmions have huge potential for applications in next-generation spintronic devices. Defects naturally occurring in magnets have profound effects on the static and dynamical properties of skyrmions. In this work, we study the effect of an atomic defect on a skyrmion using the first-principles calculations within the density functional theory, taking MnSi as an example. By substituting one site of Mn or Si with different elements, we can tune the pinning energy. The effects of pinning by an atomic defect can be understood qualitatively within a phenomenological model.


title: Confidence Intervals for Projections of Partially Identified Parameters

abstract: We propose a bootstrap-based calibrated projection procedure to build confidence intervals for single components and for smooth functions of a partially identified parameter vector in moment (in)equality models. The method controls asymptotic coverage uniformly over a large class of data generating processes. The extreme points of the calibrated projection confidence interval are obtained by extremizing the value of the function of interest subject to a proper relaxation of studentized sample analogs of the moment (in)equality conditions. The degree of relaxation, or critical level, is calibrated so that the function of theta, not theta itself, is uniformly asymptotically covered with prespecified probability. This calibration is based on repeatedly checking feasibility of linear programming problems, rendering it computationally attractive. Nonetheless, the program defining an extreme point of the confidence interval is generally nonlinear and potentially intricate. We provide an algorithm, based on the response surface method for global optimization, that approximates the solution rapidly and accurately, and we establish its rate of convergence. The algorithm is of independent interest for optimization problems with simple objectives and complicated constraints. An empirical application estimating an entry game illustrates the usefulness of the method. Monte Carlo simulations confirm the accuracy of the solution algorithm, the good statistical as well as computational performance of calibrated projection (including in comparison to other methods), and the algorithm's potential to greatly accelerate computation of other confidence intervals.


title: Generic properties of magnetic flows

abstract: We obtain Kupka-Smale theorem and Franks lemma for magnetic flows on manifolds with any dimension. This improves Miranda result on surfaces. However our methods relies on geometric control theory, like in Rifford and Ruggiero articles.


title: Detectable singularities from dynamic Radon data

abstract: In this paper, we use microlocal analysis to understand what X-ray tomographic data acquisition does to singularities of an object which changes during the measuring process. Depending on the motion model, we study which singularities are detected by the measured data. In particular, this analysis shows that, due to the dynamic behavior, not all singularities might be detected, even if the radiation source performs a complete turn around the object. Thus, they cannot be expected to be (stably) visible in any reconstruction. On the other hand, singularities could be added (or masked) as well. To understand this precisely, we provide a characterization of visible and added singularities by analyzing the microlocal properties of the forward and reconstruction operators. We illustrate the characterization using numerical examples.


title: Distribution and mass of diffuse and dense CO gas in the Milky Way

abstract: Emission from carbon monoxide (CO) is ubiquitously used as a tracer of dense star forming molecular clouds. There is, however, growing evidence that a significant fraction of CO emission originates from diffuse molecular gas. Quantifying the contribution of diffuse CO-emitting gas is vital for understanding the relation between molecular gas and star formation. We examine the Galactic distribution of two CO-emitting gas components, a high column density component detected in 13CO and 12CO, and a low column density component detected in 12CO, but not in 13CO. The "diffuse" and "dense" components are identified using a combination of smoothing, masking, and erosion/dilation procedures, making use of three large-scale 12CO and 13CO surveys of the Inner and Outer Milky Way. The diffuse component, which globally represents 25 (1.5x1e8 Mo) of the total molecular gas mass (6.5x1e8 Mo), is more extended perpendicular to the Galactic plane. The fraction of diffuse gas increases from 15% at a galactocentric radius of 3 kpc to 50% at 15 kpc, and increases with decreasing surface density. In the Inner Galaxy, a yet denser component traced by CS emission represents 14% of the total molecular gas mass traced by 12CO emission. Only 14% of the molecular gas mass traced by 12CO emission is identified as part of molecular clouds in 12CO surveys by cloud identification algorithms. This study indicates that CO emission not only traces star forming clouds, but also a significant diffuse molecular ISM component.


title: Sharp inequalities for one-sided Muckenhoupt weights

abstract: Let $A_\infty ^+$ denote the class of one-sided Muckenhoupt weights, namely all the weights $w$ for which $\mathsf M^+:L^p(w)\to L^{p,\infty}(w)$ for some $p>1$, where $\mathsf M^+$ is the forward Hardy-Littlewood maximal operator. We show that $w\in A_\infty ^+$ if and only if there exist numerical constants $\gamma\in(0,1)$ and $c>0$ such that $$ w(\{x \in \mathbb{R} : \, \mathsf M ^+\mathbf 1_E (x)>\gamma\})\leq c w(E) $$ for all measurable sets $E\subset \mathbb R$. Furthermore, letting $$ \mathsf C_w ^+(\alpha):= \sup_{0<w(E)<+\infty} \frac{1}{w(E)} w(\{x\in\mathbb R:\,\mathsf M^+\mathbf 1_E (x)>\alpha\}) $$ we show that for all $w\in A_\infty ^+$ we have the asymptotic estimate $\mathsf C_w ^+ (\alpha)-1\lesssim (1-\alpha)^\frac{1}{c[w]_{A_\infty ^+}}$ for $\alpha$ sufficiently close to $1$ and $c>0$ a numerical constant, and that this estimate is best possible. We also show that the reverse H\"older inequality for one-sided Muckenhoupt weights, previously proved by Mart\'in-Reyes and de la Torre, is sharp, thus providing a quantitative equivalent definition of $A_\infty ^+$. Our methods also allow us to show that a weight $w\in A_\infty ^+$ satisfies $w\in A_p ^+$ for all $p>e^{c[w]_{A_\infty ^+}}$.


title: Global analysis of nuclear parton distribution functions and their uncertainties at next-to-next-to-leading order

abstract: We perform a next-to-next-to-leading order (NNLO) analysis of nuclear parton distribution functions (nPDFs) using neutral current charged-lepton ($\ell ^\pm$ + nucleus) deeply inelastic scattering (DIS) data and Drell-Yan (DY) cross-section ratios $\sigma_{DY}^{A}/\sigma_{DY}^{A^\prime}$ for several nuclear targets. We study in detail the parametrizations and the atomic mass (A) dependence of the nuclear PDFs at this order. The present nuclear PDFs global analysis provides us a complete set of nuclear PDFs, $f_i^{(A,Z)}(x,Q^2)$, with a full functional dependence on $x$, A, Q$^2$. The uncertainties of the obtained nuclear modification factors for each parton flavour are estimated using the well-known Hessian method. The nuclear charm quark distributions are also added into the analysis. We compare the parametrization results with the available data and the results of other nuclear PDFs groups. We found our nuclear PDFs to be in reasonably good agreement with them. The estimates of errors provided by our global analysis are rather smaller than those of other groups. In general, a very good agreement is achieved. We also briefly review the recent heavy-ion collisions data including the first experimental data from the LHC proton+lead and lead+lead run which can be used in the global fits of nuclear PDFs. We highlight different aspects of the high luminosity Pb--Pb and p--Pb data which have been recorded by the CMS Collaboration.


title: Pricing barrier options with discrete dividends

abstract: The presence of discrete dividends complicates the derivation and form of pricing formulas even for vanilla options. Existing analytic, numerical, and theoretical approximations provide results of varying quality and performance. Here, we compare the analytic approach, developed and effective for European puts and calls, of Buryak and Guo with the formulas, designed in the context of barrier option pricing, of Dai and Chiu.


title: The Chandra COSMOS Legacy survey: overview and point source catalog

abstract: The COSMOS-Legacy survey is a 4.6 Ms Chandra program that has imaged 2.2 deg$^2$ of the COSMOS field with an effective exposure of $\simeq$160 ks over the central 1.5 deg$^2$ and of $\simeq$80 ks in the remaining area. The survey is the combination of 56 new observations, obtained as an X-ray Visionary Project, with the previous C-COSMOS survey. We describe the reduction and analysis of the new observations and the properties of 2273 point sources detected above a spurious probability of 2$\times 10^{-5}$. We also present the updated properties of the C-COSMOS sources detected in the new data. The whole survey includes 4016 point sources (3814, 2920 and 2440 in the full, soft and hard band). The limiting depths are 2.2 $\times$ 10$^{-16}$, 1.5 $\times$ 10$^{-15}$ and 8.9$\times$ 10$^{-16}$ ${\rm erg~cm}^{-2}~{\rm s}^{-1}$ in the 0.5-2, 2-10 and 0.5-10 keV bands, respectively. The observed fraction of obscured AGN with column density $> 10^{22}$ cm$^{-2}$ from the hardness ratio (HR) is $\sim$50$^{+17}_{-16}$%. Given the large sample, we compute source number counts in the hard and soft bands, significantly reducing the uncertainties of 5-10%. For the first time, we compute number counts for obscured (HR$>$-0.2) and unobscured (HR$<$-0.2) sources and find significant differences between the two populations in the soft band. Due to the un-precedent large exposure, COSMOS-Legacy area is 3 times larger than surveys at similar depth and its depth is 3 times fainter than surveys covering similar area. The area-flux region occupied by COSMOS-Legacy is likely to remain unsurpassed for years to come.


title: Self-consistent chaotic transport in a high-dimensional mean-field Hamiltonian map model

abstract: Self-consistent chaotic transport is studied in a Hamiltonian mean-field model. The model provides a simplified description of transport in marginally stable systems including vorticity mixing in strong shear flows and electron dynamics in plasmas. Self-consistency is incorporated through a mean-field that couples all the degrees-of-freedom. The model is formulated as a large set of $N$ coupled standard-like area-preserving twist maps in which the amplitude and phase of the perturbation, rather than being constant like in the standard map, are dynamical variables. Of particular interest is the study of the impact of periodic orbits on the chaotic transport and coherent structures. Numerical simulations show that self-consistency leads to the formation of a coherent macro-particle trapped around the elliptic fixed point of the system that appears together with an asymptotic periodic behavior of the mean field. To model this asymptotic state, we introduced a non-autonomous map that allows a detailed study of the onset of global transport. A turnstile-type transport mechanism that allows transport across instantaneous KAM invariant circles in non-autonomous systems is discussed. As a first step to understand transport, we study a special type of orbits referred to as sequential periodic orbits. Using symmetry properties we show that, through replication, high-dimensional sequential periodic orbits can be generated starting from low-dimensional periodic orbits. We show that sequential periodic orbits in the self-consistent map can be continued from trivial (uncoupled) periodic orbits of standard-like maps using numerical and asymptotic methods. Normal forms are used to describe these orbits and to find the values of the map parameters that guarantee their existence. Numerical simulations are used to verify the prediction from the asymptotic methods.

categories: math.DS cond-mat.stat-mech nlin.CD physics.plasm-ph


title: Representation of large matchings in bipartite graphs

abstract: Let $f(n)$ be the smallest number such that every collection of $n$ matchings, each of size at least $f(n)$, in a bipartite graph, has a full rainbow matching. Generalizing famous conjectures of Ryser, Brualdi and Stein, Aharoni and Berger conjectured that $f(n)=n+1$ for every $n>1$. Clemens and Ehrenm{\"u}ller proved that $f(n) \le \frac{3}{2}n +o(n)$. We show that the $o(n)$ term can be reduced to a constant, namely $f(n) \le \lceil \frac{3}{2}n \rceil+1$.


title: On the number of nonisomorphic subtrees of a tree

abstract: We show that a tree of order $n$ has at most $O(5^{n/4})$ nonisomorphic subtrees, and that this bound is best possible. We also prove an analogous result for the number of nonisomorphic rooted subtrees of a rooted tree.


title: On the inward drift of runaway electrons during the plateau phase of runaway current

abstract: The well observed inward drift of current carrying runaway electrons during runaway plateau regime after disruption is studied by considering the phase space dynamic of runaways in a large aspect ratio toroidal system. We consider the case where the toroidal field is unperturbed and the toroidal symmetry of the system is preserved. The balance between the change in canonical angular momentum and the input of mechanical angular momentum in such system requires runaways to drift horizontally in configuration space for any given change in momentum space. The dynamic of this drift can be obtained by taking the variation of canonical angular momentum. It is then found that runaway electrons will always drift inward as long as they are decelerating. This drift motion is essentially non-linear, since the current is carried by runaways themselves, and any runaway drift relative to the magnetic axis will cause further displacement of the axis itself. A simplified analytical model is constructed to describe such inward drift both in ideal wall case and no wall case, and the runaway current center displacement as a function of parallel momentum variation is obtained. The time scale of such displacement is estimated by considering effective radiation drag, which shows reasonable agreement with observed displacement time scale. This indicates that the phase space dynamic studied here plays a major role in the horizontal displacement of runaway electrons during plateau regime.


title: Groups quasi-isometric to RAAG's

abstract: We characterize groups quasi-isometric to a right-angled Artin group $G$ with finite outer automorphism group. In particular all such groups admit a geometric action on a $CAT(0)$ cube complex that has an equivariant "fibering" over the Davis building of $G$.


title: Tests of GR with INPOP15a planetary ephemerides: estimations of possible supplementary advances of perihelia for Mercury and Saturn

abstract: Planetary ephemerides are a good tool for studying general relativity at the scale of our solar system. We present here new evaluations of advances of perihelia for Mercury and Saturn.


title: Restricted invertibility revisited

abstract: Suppose that $m,n\in \mathbb{N}$ and that $A:\mathbb{R}^m\to \mathbb{R}^n$ is a linear operator. It is shown here that if $k,r\in \mathbb{N}$ satisfy $k<r\le \mathrm{\bf rank(A)}$ then there exists a subset $\sigma\subseteq \{1,\ldots,m\}$ with $|\sigma|=k$ such that the restriction of $A$ to $\mathbb{R}^{\sigma}\subseteq \mathbb{R}^m$ is invertible, and moreover the operator norm of the inverse $A^{-1}:A(\mathbb{R}^{\sigma})\to \mathbb{R}^m$ is at most a constant multiple of the quantity $\sqrt{mr/((r-k)\sum_{i=r}^m \mathsf{s}_i(A)^2)}$, where $\mathsf{s}_1(A)\geqslant\ldots\geqslant \mathsf{s}_m(A)$ are the singular values of $A$. This improves over a series of works, starting from the seminal Bourgain--Tzafriri Restricted Invertibility Principle, through the works of Vershynin, Spielman--Srivastava and Marcus--Spielman--Srivastava. In particular, this directly implies an improved restricted invertibility principle in terms of Schatten--von Neumann norms.


title: On the Diversity of Compact Objects within Supernova Remnants. I: A Parametric Model for Magnetic Field Evolution

abstract: A wealth of X-ray and radio observations has revealed in the past decade a growing diversity of neutron stars (NSs) with properties spanning orders of magnitude in magnetic field strength and ages, and with emission processes explained by a range of mechanisms dictating their radiation properties. However, serious difficulties exist with the magneto-dipole model of isolated neutron star fields and their inferred ages, such as a large range of observed braking indices ($n$, with values often $<$3) and a mismatch between the neutron star and associated supernova remnant (SNR) ages. This problem arises primarily from the assumptions of a constant magnetic field with $n$=3, and an initial spin period that is much smaller than the observed current period. It has been suggested that a solution to this problem involves magnetic field evolution, with some NSs having magnetic fields buried within the crust by accretion of fall-back supernova material following their birth. In this work we explore a parametric phenomenological model for magnetic field growth that generalizes previous suggested field evolution functions, and apply it to a variety of NSs with both secure SNR associations and known ages. We explore the flexibility of the model by recovering the results of previous work on buried magnetic fields in young neutron stars. Our model fits suggest that apparently disparate classes of NSs may be related to one another through the time-evolution of the magnetic field.


title: Odd zeta motive and linear forms in odd zeta values

abstract: We study a family of mixed Tate motives over $\mathbb{Z}$ whose periods are linear forms in the zeta values $\zeta(n)$. They naturally include the Beukers-Rhin-Viola integrals for $\zeta(2)$ and the Ball-Rivoal linear forms in odd zeta values. We give a general integral formula for the coefficients of the linear forms and a geometric interpretation of the vanishing of the coefficients of a given parity. The main underlying result is a geometric construction of a minimal ind-object in the category of mixed Tate motives over $\mathbb{Z}$ which contains all the non-trivial extensions between simple objects. In a joint appendix with Don Zagier, we prove the compatibility between the structure of the motives considered here and the representations of their periods as sums of series.


title: Higher Moments of Net-Particle Multiplicity Distributions

abstract: Studying fluctuations of conserved quantities, such as baryon number, strangeness, and charge, provides insights into the properties of matter created in high-energy nuclear collisions. Lattice QCD calculations suggest that higher moments of these quantities are sensitive to the phase structure of the hot and dense nuclear matter created in such collisions. In this paper, we present first experimental results of volume and temperature independent cumulant ratios of net-charge and net-proton distributions in Au+Au collisions at $\sqrt{s_{\rm NN}}\!=\!14.5 \,{\rm GeV}$ completing the first RHIC Beam Energy Scan (BES-I) program for $\sqrt{s_{\rm NN}}\!=\!7.7 \,{\rm to}\,200\,{\rm GeV}$, together with the first measurement of fully corrected net-kaon results, measured with the STAR detector at RHIC at mid-rapidity and a transverse momentum up to $p_{\rm T}\!=\!2\,{\rm GeV/{\it c}}$. The pseudo-rapidity dependence of the $\sqrt{s_{\rm NN}}\!=\!14.5 \,{\rm GeV}$ net-charge cumulant ratios is discussed. The estimated uncertainties on the ratio $c_4/c_2$, the most statistics-hungry of the present observables, at $\sqrt{s_{\rm NN}}\!=\!7.7 \,{\rm GeV}$ in the upcoming RHIC BES-II program will also be presented.


title: Implications of the diphoton excess on Left-Right models and gauge unification

abstract: The recent diphoton excess signal at an invariant mass of 750 GeV can be interpreted in the framework of left-right symmetric models with additional scalar singlets and vector-like fermions. We propose a minimal scenario for such a purpose. Extending the LRSM framework to include these new vector-like fermionic fields, on the other hand, results in interesting phenomenological implications for the LRSM fermion masses and mixing. Furthermore, existence of such vector-like fermions can also have interesting implications for baryogenesis and the dark matter sector. The introduction of a real bi-triplet scalar which contains a potential DM candidate will allow the gauge couplings to unify at $\approx 10^{17.7}$ GeV.


title: Large Magellanic Cloud Near-Infrared Synoptic Survey. III. A Statistical Study of Non-Linearity in the Leavitt Laws

abstract: We present a detailed statistical analysis of possible non-linearities in the Period-Luminosity (P-L), Period-Wesenheit (P-W) and Period-Color (P-C) relations for Cepheid variables in the LMC at optical ($VI$) and near-infrared ($JHK_{s}$) wavelengths. We test for the presence of possible non-linearities and determine their statistical significance by applying a variety of robust statistical tests ($F$-test, Random-Walk, Testimator and the Davies test) to optical data from OGLE III and near-infrared data from LMCNISS. For fundamental-mode Cepheids, we find that the optical P-L, P-W and P-C relations are non-linear at 10 days. The near-infrared P-L and the $W^H_{V,I}$ relations are non-linear around 18 days; this break is attributed to a distinct variation in mean Fourier amplitude parameters near this period for longer wavelengths as compared to optical bands. The near-infrared P-W relations are also non-linear except for the $W_{H,K_s}$ relation. For first-overtone mode Cepheids, a significant change in the slope of P-L, P-W and P-C relations is found around 2.5 days only at optical wavelengths. We determine a global slope of $\textrm{-}3.212\pm0.013$ for the $W^H_{V,I}$ relation by combining our LMC data with observations of Cepheids in Supernovae host galaxies \citep{riess11}. We find this slope to be consistent with the corresponding LMC relation at short periods, and significantly different to the long-period value. We do not find any significant difference in the slope of the global-fit solution using a linear or non-linear LMC P-L relation as calibrator, but the linear version provides a $2\times$ better constraint on the slope and metallicity coefficient.


title: Jamming and Attraction of Interacting Run-and-Tumble Random Walkers

abstract: We study a model of bacterial dynamics where two interacting random walkers perform run-and-tumble motion on a one-dimensional lattice under mutual exclusion and find an exact expression for the probability distribution in the steady state. This stationary distribution has a rich structure comprising three components: a jammed component, where the particles are adjacent and block each other; an attractive component, where the probability distribution for the distance between particles decays exponentially; and an extended component in which the distance between particles is uniformly distributed. The attraction between the particles is sufficiently strong that even in the limit where continuous space is recovered for a finite system, the two walkers spend a finite fraction of time in a jammed configuration. Our results potentially provide a route to understanding the motility-induced phase separation characteristic of active matter from a microscopic perspective.


title: Optimally Pruning Decision Tree Ensembles With Feature Cost

abstract: We consider the problem of learning decision rules for prediction with feature budget constraint. In particular, we are interested in pruning an ensemble of decision trees to reduce expected feature cost while maintaining high prediction accuracy for any test example. We propose a novel 0-1 integer program formulation for ensemble pruning. Our pruning formulation is general - it takes any ensemble of decision trees as input. By explicitly accounting for feature-sharing across trees together with accuracy/cost trade-off, our method is able to significantly reduce feature cost by pruning subtrees that introduce more loss in terms of feature cost than benefit in terms of prediction accuracy gain. Theoretically, we prove that a linear programming relaxation produces the exact solution of the original integer program. This allows us to use efficient convex optimization tools to obtain an optimally pruned ensemble for any given budget. Empirically, we see that our pruning algorithm significantly improves the performance of the state of the art ensemble method BudgetRF.


title: A solution of Gromov's H\"older equivalence problem for the Heisenberg group

abstract: We show that a map with H\"older exponent bigger than $1/2$ from a quasi-convex metric space with vanishing first Lipschitz homology into the Sub-Riemannian Heisenberg group factors through a tree. In particular, if the domain contains a disk, such a map can't be injective. This gives an answer to a question of Gromov for the simplest nontrivial case. The same tools allow to improve on a result of Borisov and it is shown that an isometric immersion of class $C^{1,\alpha}$ of a Riemannian surface with positive Gauss curvature into $\mathbb{R}^3$ has bounded extrinsic curvature if $\alpha > 1/2$.


title: Infinity Laplacian equation with strong absorptions

abstract: We study regularity properties of solutions to reaction-diffusion equations ruled by the infinity laplacian operator. We focus our analysis in models presenting plateaus, i.e. regions where a non-negative solution vanishes identically. We obtain sharp geometric regularity estimates for solutions along the boundary of plateaus sets. In particular we show that the $(n-\epsilon)$-Hausdorff measure of the plateaus boundary is finite, for a universal number $\epsilon>0$.


title: Chiral spin liquids in triangular lattice SU(N) fermionic Mott insulators with artificial gauge fields

abstract: We show that, in the presence of a $\pi/2$ artificial gauge field per plaquette, Mott insulating phases of ultra-cold fermions with $SU(N)$ symmetry and one particle per site generically possess an extended chiral phase with intrinsic topological order characterized by a multiplet of $N$ low-lying singlet excitations for periodic boundary conditions, and by chiral edge states described by the $SU(N)_1$ Wess-Zumino-Novikov-Witten conformal field theory for open boundary conditions. This has been achieved by extensive exact diagonalizations for $N$ between $3$ and $9$, and by a parton construction based on a set of $N$ Gutzwiller projected fermionic wave-functions with flux $\pi/N$ per triangular plaquette. Experimental implications are briefly discussed.


title: SU(6) Heisenberg model on the honeycomb lattice: competition between plaquette and chiral order

abstract: We revisit the SU(6) Heisenberg model on the honeycomb lattice, which has been predicted to be a chiral spin liquid by mean-field theory [G. Szirmai et al., Phys. Rev. A 84, 011611 (2011)]. Using exact diagonalizations of finite clusters, infinite projected entangled pair states simulations, and variational Monte Carlo simulations based on Gutzwiller projected wave functions, we provide strong evidence in favour of the competing plaquette state, which was reported to be higher but close by in energy according to mean-field theory. This is further confirmed by the investigation of the model with a ring exchange term, which shows that there is a transition between the plaquette state and the chiral state at a finite value of the ring exchange term.


title: High Frequency Remote Monitoring of Parkinson's Disease via Smartphone: Platform Overview and Medication Response Detection

abstract: Objective: The aim of this study is to develop a smartphone-based high-frequency remote monitoring platform, assess its feasibility for remote monitoring of symptoms in Parkinson's disease, and demonstrate the value of data collected using the platform by detecting dopaminergic medication response. Methods: We have developed HopkinsPD, a novel smartphone-based monitoring platform, which measures symptoms actively (i.e. data are collected when a suite of tests is initiated by the individual at specific times during the day), and passively (i.e. data are collected continuously in the background). After data collection, we extract features to assess measures of five key behaviors related to PD symptoms -- voice, balance, gait, dexterity, and reaction time. A random forest classifier is used to discriminate measurements taken after a dose of medication (treatment) versus before the medication dose (baseline). Results: A worldwide study for remote PD monitoring was established using HopkinsPD in July, 2014. This study used entirely remote, online recruitment and installation, demonstrating highly cost-effective scalability. In six months, 226 individuals (121 PD and 105 controls) contributed over 46,000 hours of passive monitoring data and approximately 8,000 instances of structured tests of voice, balance, gait, reaction, and dexterity. To the best of our knowledge, this is the first study to have collected data at such a scale for remote PD monitoring. Moreover, we demonstrate the initial ability to discriminate treatment from baseline with 71.0(+-0.4)% accuracy, which suggests medication response can be monitored remotely via smartphone-based measures.


title: Tunable and direction-dependent group velocities in topologically protected edge states

abstract: Topological effects continue to fascinate physicists since more than three decades. One of their main applications are high-precision measurements of the resistivity. We propose to make also use of the spatially separated edge states. It is possible to realize strongly direction-dependent group velocities. They can also be tuned over orders of magnitude so that robust delay lines and interference devices are within reach.


title: Einstein-Podolsky-Rosen correlations and Bell correlations in the simplest scenario

abstract: Einstein-Podolsky-Rosen (EPR) steering is an intermediate type of quantum nonlocality which sits between entanglement and Bell nonlocality. A set of correlations is Bell nonlocal if it does not admit a local hidden variable (LHV) model, while it is EPR nonlocal if it does not admit a local hidden variable-local hidden state (LHV-LHS) model. It is interesting to know what states can generate EPR-nonlocal correlations in the simplest nontrivial scenario, that is, two projective measurements for each party sharing a two-qubit state. Here we show that a two-qubit state can generate EPR-nonlocal full correlations (excluding marginal statistics) in this scenario if and only if it can generate Bell-nonlocal correlations. If full statistics (including marginal statistics) is taken into account, surprisingly, the same scenario can manifest the simplest one-way steering and the strongest hierarchy between steering and Bell nonlocality. To illustrate these intriguing phenomena in simple setups, several concrete examples are discussed in detail, which facilitates experimental demonstration. In the course of study, we introduce the concept of restricted LHS models and thereby derive a necessary and sufficient semidefinite-programming criterion to determine the steerability of any bipartite state under given measurements. Analytical criteria are further derived in several scenarios of strong theoretical and experimental interest.


title: Identifying System-Wide Early Warning Signs of Instability in Stochastic Power Systems

abstract: Prior research has shown that spectral decomposition of the reduced power flow Jacobian (RPFJ) can yield participation factors that describe the extent to which particular buses contribute to particular spectral components of a power system. Research has also shown that both variance and autocorrelation of time series voltage data tend to increase as a power system with stochastically fluctuating loads approaches certain critical transitions. This paper presents evidence suggesting that a system's participation factors predict the relative bus voltage variance values for all nodes in a system. As a result, these participation factors can be used to filter, weight, and combine real time PMU data from various locations dispersed throughout a power network in order to develop coherent measures of global voltage stability. This paper first describes the method of computing the participation factors. Next, two potential uses of the participation factors are given: (1) predicting the relative bus voltage variance magnitudes, and (2) locating generators at which the autocorrelation of voltage measurements clearly indicate proximity to critical transitions. The methods are tested using both analytical and numerical results from a dynamic model of a 2383-bus test case.


title: Extremely Isolated Galaxies I. Sample and Simulation Analysis

abstract: We have selected a sample of extremely isolated galaxies (EIGs) from the local Universe ($\mbox{z} < 0.024$), using a simple isolation criterion: having no known neighbours closer than $300\,{\rm km\,s}^{-1}$ ($3\,h^{-1}\,\mbox{Mpc}$) in the three-dimensional redshift space $(\alpha,\delta,\mbox{z})$. The sample is unique both in its level of isolation and in the fact that it utilizes HI redshifts from the Arecibo Legacy Fast ALFA survey (ALFALFA). We analysed the EIG sample using cosmological simulations and found that it contains extremely isolated galaxies with normal mass haloes which have evolved gradually with little or no "major events" (major mergers, or major mass-loss events) in the last $3\,\mbox{Gyr}$. The fraction of EIGs which deviate from this definition (false positives) is 5%-10%. For the general population of dark matter haloes it was further found that the mass accretion (relative to the current halo mass) is affected by the halo environment mainly through strong interactions with its neighbours. As long as a halo does not experience major events, its Mass Accretion History (MAH) does not depend significantly on its environment. "Major events" seem to be the main mechanism that creates low-mass subhaloes ($M_{halo} < 10^{10}\,h^{-1}\,M_\odot$) that host galaxies (with $\mbox{M}_{g} \lesssim -14$).


title: Astrometry with MCAO at Gemini and at ELTs

abstract: We present in this study a first analysis of the astrometric error budget of absolute astrometry relative to background galaxies using adaptive optics. We use for this analysis multi-conjugated adaptive optics (MCAO) images obtained with GeMS/GSAOI at Gemini South. We find that it is possible to obtain 0.3 mas reference precision in a random field with 1 hour on source using faint background galaxies. Systematic errors are correctable below that level, such that the overall error is approximately 0.4 mas. Because the reference sources are extended, we find it necessary to correct for the dependency of the PSF centroid on the used aperture size, which would otherwise cause an important bias. This effect needs also to be considered for Extremely Large Telescopes (ELTs). When this effect is corrected, ELTs have the potential to measure proper motions of dwarfs galaxies around M31 with 10 km/s accuracy over a baseline of 5 years.


title: Capacities of repeater-assisted quantum communications

abstract: We consider quantum and private communications assisted by repeaters, from the basic scenario of a single repeater chain to the general case of an arbitrarily-complex quantum network, where systems may be routed through single or multiple paths. In this context, we investigate the ultimate rates at which two end-parties may transmit quantum information, distribute entanglement, or generate secret keys. These end-to-end capacities are defined by optimizing over the most general adaptive protocols that are allowed by quantum mechanics. Combining techniques from quantum information and classical network theory, we derive single-letter upper bounds for the end-to-end capacities in repeater chains and quantum networks connected by arbitrary quantum channels, establishing exact formulas under basic decoherence models, including bosonic lossy channels, quantum-limited amplifiers, dephasing and erasure channels. For the converse part, we adopt a teleportation-inspired simulation of a quantum network which leads to upper bounds in terms of the relative entropy of entanglement. For the lower bounds we combine point-to-point quantum protocols with classical network algorithms. Depending on the type of routing (single or multiple), optimal strategies corresponds to finding the widest path or the maximum flow in the quantum network. Our theory can also be extended to simultaneous quantum communication between multiple senders and receivers.

categories: quant-ph cond-mat.other math-ph math.MP physics.optics


title: NIHAO VII: Predictions for the galactic baryon budget in dwarf to Milky Way mass haloes

abstract: We use the NIHAO galaxy formation simulations to make predictions for the baryonic budget in present day galaxies ranging from dwarf to Milky Way masses. The sample is made of 88 independent high resolution cosmological zoom-in simulations. NIHAO galaxies reproduce key properties of observed galaxies, such as the stellar mass vs halo mass and cold gas vs stellar mass relations. Thus they make plausible predictions for the baryon budget. We present the mass fractions of stars, cold gas ($T<10^4$K), cool gas ($10^4 < T < 10^5$K), warm-hot gas ($10^5 < T < 5\times10^6$K), and hot gas (T$> 5\times10^6$K) inside the virial radius, $R_{200}$. Compared to the predicted baryon mass, using the dark halo mass and the universal baryon fraction, $f_{\rm b}\equiv \Omega_{\rm b}/\Omega_{\rm m}=0.15$, we find that all of our haloes are missing baryons. The missing mass has been relocated past 2 virial radii, and cool gas dominates the corona at low mass (M$_{200} < 3 \times 10^{11} M_\odot$) while the warm-hot gas dominates at high mass (M$_{200} > 3 \times 10^{11} M_\odot$). Haloes of mass $\sim 10^{10} M_\odot$ are missing $\sim 90\%$ of their baryons. More massive haloes ($\sim 10^{12} M_\odot$) retain a higher fraction of their baryons, with $\sim 30\%$ missing, consistent with recent observational estimates. Moreover, these more massive haloes reproduce the observed fraction of cold, warm-hot and hot gas. The fraction of cool gas we predict ($0.11\pm0.06$) is significantly lower than the observation from COS-HALOs (0.3-0.47), but agrees with the alternative analysis of Stern et al. 2016.


title: Chiral spin liquid from magnetic Wannier states

abstract: We present a mapping of a two-dimensional system of interacting bosons in a strong perpendicular magnetic field to an equivalent system of interacting bosons on the square lattice in the absence of the field. The mapping utilizes a magnetic Bloch and the corresponding magnetic Wannier single-particle basis in the lowest Landau level. By construction, the ground states of the resulting model of interacting bosons on the square lattice are gapped fractionalized liquids or gapless Bose metal states with broken time reversal symmetry at specific rational filling fractions.


title: Homomorphisms of Strongly Regular Graphs

abstract: We prove that if $G$ and $H$ are primitive strongly regular graphs with the same parameters and $\varphi$ is a homomorphism from $G$ to $H$, then $\varphi$ is either an isomorphism or a coloring (homomorphism to a complete subgraph). Therefore, the only endomorphisms of a primitive strongly regular graph are automorphisms or colorings. This confirms and strengthens a conjecture of Cameron and Kazanidis that all strongly regular graphs are cores or have complete cores. The proof of the result is elementary, mainly relying on linear algebraic techniques. In the second half of the paper we discuss implications of the result and the idea underlying the proof. We also show that essentially the same proof can be used to obtain a more general statement.


title: Lepton flavor (universality) violation in rare kaon decays

abstract: Recent anomalies in the decays of $B$ mesons and the Higgs boson provide hints towards lepton flavor (universality) violating physics beyond the Standard Model. We observe that four-fermion operators which can explain the $B$-physics anomalies have corresponding analogs in the kaon sector, and we analyze their impact on $K\to\pi\ell \ell'$ and $K\to\ell \ell'$ decays $(\ell=\mu,e)$. For these processes, we note the corresponding physics opportunities at the NA62 experiment. In particular, assuming minimal flavor violation, we comment on the required improvements in sensitivity necessary to test the $B$-physics anomalies in the kaon sector.


title: Mapping kiloparsec-scale structures in the extended HI disc of the galaxy UGC 00439 by HI 21-cm absorption

abstract: We study the properties of HI gas in the outer regions (~2r_25) of a spiral galaxy, UGC 00439 (z = 0.01769), using HI 21-cm absorption towards different components of an extended background radio source, J0041$-$0043 (z = 1.679). The radio source exhibits a compact core coincident with the optical quasar and two lobes separated by ~7 kpc, all at an impact parameter ~25 kpc. The HI 21-cm absorption detected towards the southern lobe is found to extend over ~2 kpc^2. The absorbing gas shows sub-kpc-scale structures with the line-of-sight velocities dominated by turbulent motions. Much larger optical depth variations over 4-7 kpc-scale are revealed by the non-detection of HI 21-cm absorption towards the radio core and the northern lobe, and the detection of NaI and CaII absorption towards the quasar. This could reflect a patchy distribution of cold gas in the extended HI disc. We also detect HI 21-cm emission from UGC 00439 and two other galaxies within ~150 kpc to it, that probably form an interacting group. However, no HI 21-cm emission from the absorbing gas is detected. Assuming a linear extent of ~4 kpc, as required to cover both the core and the southern lobe, we constrain the spin temperature <~300 K for the absorbing gas. The kinematics of the gas and the lack of signatures of any ongoing in situ star formation are consistent with the absorbing gas being at the kinematical minor axis and corotating with the galaxy. Deeper HI 21-cm observations would help to map in greater detail both the large- and small-scale structures in the HI gas associated with UGC 00439.


title: Updated determination of the solar neutrino fluxes from solar neutrino data

abstract: We present an update of the determination of the solar neutrino fluxes from a global analysis of the solar and terrestrial neutrino data in the framework of three-neutrino mixing. Using a Bayesian analysis we reconstruct the posterior probability distribution function for the eight normalization parameters of the solar neutrino fluxes plus the relevant masses and mixing, with and without imposing the luminosity constraint. We then use these results to compare the description provided by different Standard Solar Models. Our results show that, at present, both models with low and high metallicity can describe the data with equivalent statistical agreement. We also argue that even with the present experimental precision the solar neutrino data have the potential to improve the accuracy of the solar model predictions.


title: Elastic and Transition Form Factors in DSEs

abstract: A symmetry preserving framework for the study of continuum Quantum Chromodynamics (QCD) is obtained from a truncated solution of the QCD equations of motion or QCD's Dyson-Schwinger equations (DSEs). A nonperturbative solution of the DSEs enables the study of, e.g., hadrons as composites of dressed-quarks and dressed-gluons, the phenomena of confinement and dynamical chiral symmetry breaking (DCSB), and therefrom an articulation of any connection between them. It is within this context that we present a unified study of Nucleon, Delta and Roper elastic and transition form factors, and compare predictions made using a framework built upon a Faddeev equation kernel and interaction vertices that possess QCD-like momentum dependence with results obtained using a symmetry-preserving treatment of a vector$\,\otimes\,$vector contact-interaction.


title: The SINFONI Black Hole Survey: The Black Hole Fundamental Plane revisited and the paths of (co-) evolution of supermassive black holes and bulges

abstract: We investigate the correlations between the black hole mass $M_{BH}$, the velocity dispersion $\sigma$, the bulge mass $M_{Bu}$, the bulge average spherical density $\rho_h$ and its spherical half mass radius $r_h$, constructing a database of 97 galaxies (31 core ellipticals, 17 power-law ellipticals, 30 classical bulges, 19 pseudo bulges) by joining 72 galaxies from the literature to 25 galaxies observed during our recent SINFONI black hole survey. For the first time we discuss the full error covariance matrix. We analyse the well known $M_{BH}-\sigma$ and $M_{BH}-M_{Bu}$ relations and establish the existence of statistically significant correlations between $M_{Bu}$ and $r_h$ and anti-correlations between $M_{Bu}$ and $\rho_h$. We establish five significant bivariate correlations ($M_{BH}-\sigma-\rho_h$, $M_{BH}-\sigma-r_h$, $M_{BH}-M_{Bu}-\sigma$, $M_{BH}-M_{Bu}-\rho_h$, $M_{BH}-M_{Bu}-r_h$) that predict $M_{BH}$ of 77 core and power-law ellipticals and classical bulges with measured and intrinsic scatter as small as $\approx 0.36$ dex and $\approx 0.33$ dex respectively, or 0.26 dex when the subsample of 45 galaxies defined by Kormendy and Ho (2013) is considered. In contrast, pseudo bulges have systematically lower $M_{BH}$, but approach the predictions of all the above relations at spherical densities $\rho_h\ge 10^{10} M_\odot/kpc^3$ or scale lengths $r_h\le 1$ kpc. These findings fit in a scenario of co-evolution of BH and classical-bulge masses, where core ellipticals are the product of dry mergers of power-law bulges and power-law Es and bulges the result of (early) gas-rich mergers and of disk galaxies. In contrast, the (secular) growth of BHs is decoupled from the growth of their pseudo bulge hosts, except when (gas) densities are high enough to trigger the feedback mechanism responsible for the existence of the correlations between $M_{BH}$ and galaxy structural parameters.


title: Topological Phase Transitions in Line-nodal Superconductors

abstract: Fathoming interplay between symmetry and topology of many-electron wave-functions has deepened understanding of quantum many body systems, especially after the discovery of topological insulators. Topology of electron wave-functions enforces and protects emergent gapless excitations, and symmetry is intrinsically tied to the topological protection in a certain class. Namely, unless the symmetry is broken, the topological nature is intact. We show novel interplay phenomena between symmetry and topology in topological phase transitions associated with line-nodal superconductors. The interplay may induce an exotic universality class in sharp contrast to that of the phenomenological Landau-Ginzburg theory. Hyper-scaling violation and emergent relativistic scaling are main characteristics, and the interplay even induces unusually large quantum critical region. We propose characteristic experimental signatures around the phase transitions in three spatial dimensions, for example, a linear phase boundary in a temperature-tuning parameter phase-diagram.


title: How plausible are the proposed formation scenarios of CEMP-r/s stars?

abstract: CEMP-$r/s$ stars are metal-poor stars with enhanced abundances of carbon and heavy elements associated with the slow ($s$-) and rapid ($r$-) neutron-capture process. It is believed that carbon and $s$-elements were accreted from the wind of an AGB primary star, a scenario that is generally accepted to explain the formation of CEMP stars that are only enhanced in $s$-elements (CEMP-$s$ stars). The origin of $r$-element-enrichment in CEMP-$r/s$ stars is debated and many formation scenarios have been put forward. We aim to determine the likelihood of the scenarios proposed to explain the formation of CEMP-$r/s$ stars. We calculate the frequency of CEMP-$r/s$ stars among CEMP-$s$ stars for a variety of scenarios, and we compare it with that determined from an observed sample of CEMP-$r/s$ stars collected from the literature. The theoretical frequency of CEMP-$r/s$ stars predicted in most scenarios underestimates the observed ratio by at least a factor of 5. If the enrichments in $s$- and $r$-elements are independent, the model ratio of CEMP-$r/s$ to CEMP-$s$ stars is about 22%, that is approximately consistent with the lowest estimate of the observed ratio. However, this model predicts that about one third of all carbon-normal stars have [Ba/Fe] and [Eu/Fe] higher than 1, and that 40% of all CEMP stars have [Ba/Eu]$\le0$. Stars with these properties are at least ten times rarer in our observed sample. The $intermediate$ or $i$-process, which is supposedly active in some circumstances during the AGB phase, could provide an explanation of the origin of CEMP-$r/s$ stars, similar to that of CEMP-$s$ stars, in the context of wind mass accretion in binary systems. Further calculations of the nucleosynthesis of the $i$-process and of the detailed evolution of late AGB stars are needed to investigate if this scenario predicts a CEMP-$r/s$ star frequency consistent with the observations.


title: A new and trivial CP symmetry for extended $A_4$ flavor

abstract: The combination of $\nu_\mu$-$\nu_\tau$ exchange together with CP conjugation in the neutrino sector (known as $\mathrm{CP}^{\mu\tau}$ symmetry or $\mu\tau$-reflection) is known to predict the viable pattern: $\theta_{23}=45^\circ$, maximal Dirac CP phase and trivial Majorana phases. We implement such a CP symmetry as a new CP symmetry in theories with $A_4$ flavor. The implementation in a complete renormalizable model leads to a new form for the neutrino mass matrix that leads to further predictions: normal hierarchical spectrum with lightest mass and $m_{\beta\beta}$ ($0\nu\beta\beta$) of only few meV, and either $\nu_1$ or $\nu_2$ has opposite CP parity. An approximate $L_\mu-L_\tau$ symmetry arises naturally and controls the flavor structure of the model. The light neutrino masses are generated by the extended seesaw mechanism with 6 right-handed neutrinos (RHNs). The requirement of negligible one-loop corrections to light neutrino masses, validity of the extended seesaw approximation and not too long-lived BSM states to comply with BBN essentially restricts the parameters of the model to a small region: three relatively light right-handed neutrinos at the GeV scale, heavier neutrinos at the electroweak scale and Yukawa couplings smaller than the electron Yukawa. Such a small Yukawa couplings render these RHNs unobservable in terrestrial experiments.