title: Fine structure of an exciton coupled to a single Fe2+ ion in a CdSe/ZnSe quantum dot

abstract: We present a polarization-resolved photoluminescence study of the exchange interaction effects in a prototype system consisting of an individual Fe2+ ion and a single neutral exciton confined in a CdSe/ZnSe quantum dot. Maximal possible number of eight fully linearly-polarized lines in the bright exciton emission spectrum is observed, evidencing complete degeneracy lifting in the investigated system. We discuss conditions required for such a scenario to take place: anisotropy of the electron-hole interaction and the zero-field splitting of the Fe2+ ion spin states. Neglecting of either of these components is shown to restore partial degeneracy of the transitions, making the excitonic spectrum similar to those previously reported for all other systems of quantum dots with single magnetic dopants.


title: An optimal XP algorithm for Hamiltonian cycle on graphs of bounded clique-width

abstract: In this paper, we prove that, given a clique-width $k$-expression of an $n$-vertex graph, \textsc{Hamiltonian Cycle} can be solved in time $n^{\mathcal{O}(k)}$. This improves the naive algorithm that runs in time $n^{\mathcal{O}(k^2)}$ by Espelage et al. (WG 2001), and it also matches with the lower bound result by Fomin et al. that, unless the Exponential Time Hypothesis fails, there is no algorithm running in time $n^{o(k)}$ (SIAM. J. Computing 2014). We present a technique of representative sets using two-edge colored multigraphs on $k$ vertices. The essential idea is that, for a two-edge colored multigraph, the existence of an Eulerian trail that uses edges with different colors alternately can be determined by two information: the number of colored edges incident with each vertex, and the connectedness of the multigraph. With this idea, we avoid the bottleneck of the naive algorithm, which stores all the possible multigraphs on $k$ vertices with at most $n$ edges.


title: Stable quotients and the holomorphic anomaly equation

abstract: We study the fundamental relationship between stable quotient invariants and the B-model for local CP2 in all genera. Our main result is a direct geometric proof of the holomorphic anomaly equation in the precise form predicted by B-model physics. The method yields new holomorphic anomaly equations for an infinite class of twisted theories on projective spaces. An example of such a twisted theory is the formal quintic defined by a hyperplane section of CP4 in all genera via the Euler class of a complex. The formal quintic theory is found to satisfy the holomorphic anomaly equations conjectured for the true quintic theory. Therefore, the formal quintic theory and the true quintic theory should be related by transformations which respect the holomorphic anomaly equations.


title: Scanning the Earth with solar neutrinos and DUNE

abstract: We explore oscillations of the solar $^8$B neutrinos in the Earth in detail. The relative excess of night $\nu_e$ events (the Night-Day asymmetry) is computed as function of the neutrino energy and the nadir angle $\eta$ of its trajectory. The finite energy resolution of the detector causes an important attenuation effect, while the layer-like structure of the Earth density leads to an interesting parametric suppression of the oscillations. Different features of the $\eta-$ dependence encode information about the structure (such as density jumps) of the Earth density profile; thus measuring the $\eta$ distribution allows the scanning of the interior of the Earth. We estimate the sensitivity of the DUNE experiment to such measurements. About 75 neutrino events are expected per day in 40 kt. For high values of $\Delta m^2_{21}$ and $E_\nu > $11 MeV, the corresponding D-N asymmetry is about 4\% and can be measured with $15\%$ accuracy after 5 years of data taking. The difference of the D-N asymmetry between high and low values of $\Delta m^2_{21}$ can be measured at the $4\sigma$ level. The relative excess of the $\nu_e$ signal varies with the nadir angle up to 50\%. DUNE may establish the existence of the dip in the $\eta-$ distribution at the $(2 - 3) \sigma$ level.


title: Temporal-multiplexing interferometry applied to co-phased profilometry

abstract: Fringe-projection profilometry with 1 camera and 1 fringe-projector is a well-known and widely used technique in optical metrology. Spatial-frequency multiplexing interferometry with several spatial-carriers having non-overlapping spatial-spectra is also well known and productive in optical metrology. In this paper we propose temporal-multiplexing phase-shifting interferometry applied to profilometry. That is, instead of having fringe-patterns with well separated spatial-spectra, we propose instead to separate the fringe information in the temporal-spectra. In other words, we may have overlapping spatial-spectra, but separated in the temporal-spectra by frequency multiplexing. Using 1-camera and several fringe-projectors one minimizes the object shadows and specular reflections from the digitizing solid. Temporal multiplexing profilometry allows us to illuminate the object from several projectors turned-on simultaneously. In previous phase-shifting co-phased profilometry, the projectors were turned-on and off sequentially. As seen in this work temporal-multiplexing allow us to demodulate the several fringe-patterns without crosstalk from other simultaneously projected fringes. This is entirely analogous to having several television stations broadcasting simultaneously, each TV-transmitter having its own broadcasting frequency. A given TV-receiver tunes into a single TV-station and filter-out all other broadcasters. Following this analogy, each fringe-projector must have its own temporal broadcasting frequency to remain well separated from all other projectors in the time-spectra domain. In addition to the general theory presented, we assess its feasibility with experimental results.


title: Preemptive Online Partitioning of Sequences

abstract: Online algorithms process their inputs piece by piece, taking irrevocable decisions for each data item. This model is too restrictive for most partitioning problems, since data that is yet to arrive may render it impossible to extend partial partitionings to the entire data set reasonably well. In this work, we show that preemption might be a potential remedy. We consider the problem of partitioning online sequences, where $p-1$ separators need to be inserted into a sequence of integers that arrives online so as to create $p$ contiguous partitions of similar weight. While without preemption no algorithm with non-trivial competitive ratio is possible, if preemption is allowed, i.e., inserted partition separators may be removed but not reinserted again, then we show that constant competitive algorithms can be obtained. Our contributions include: We first give a simple deterministic $2$-competitive preemptive algorithm for arbitrary $p$ and arbitrary sequences. Our main contribution is the design of a highly non-trivial partitioning scheme, which, under some natural conditions and $p$ being a power of two, allows us to improve the competitiveness to $1.68$. We also show that the competitiveness of deterministic (randomized) algorithms is at least $\frac{4}{3}$ (resp. $\frac{6}{5}$). For $p=2$, the problem corresponds to the interesting special case of preemptively guessing the center of a weighted request sequence. While deterministic algorithms fail here, we provide a randomized $1.345$-competitive algorithm for all-ones sequences and prove that this is optimal. For weighted sequences, we give a $1.628$-competitive algorithm and a lower bound of $1.5$.


title: Off-shell persistence of composite pions and kaons

abstract: In order for a Sullivan-like process to provide reliable access to a meson target as $t$ becomes spacelike, the pole associated with that meson should remain the dominant feature of the quark-antiquark scattering matrix and the wave function describing the related correlation must evolve slowly and smoothly. Using continuum methods for the strong-interaction bound-state problem, we explore and delineate the circumstances under which these conditions are satisfied: for the pion, this requires $-t \lesssim 0.6\,$GeV$^2$, whereas $-t\lesssim 0.9\,$GeV$^2$ will suffice for the kaon. These results should prove useful in evaluating the potential of numerous experiments at existing and proposed facilities.


title: Getting Super-Excited with Modified Dispersion Relations

abstract: We demonstrate that in some regions of parameter space, modified dispersion relations can lead to highly populated excited states, which we dub as "super-excited" states. In order to prepare such super-excited states, we invoke dispersion relations that have negative slope in an interim sub-horizon phase at high momenta. This behaviour of quantum fluctuations can lead to large corrections relative to the Bunch-Davies power spectrum, which mimics highly excited initial conditions. We identify the Bogolyubov coefficients that can yield these power spectra. In the course of this computation, we also point out the shortcomings of the gluing method for evaluating the power spectrum and the Bogolyubov coefficients. As we discuss, there are other regions of parameter space, where the power spectrum does not get modified. Therefore, modified dispersion relations can also lead to so-called "calm excited states" as well. We conclude by commenting on the possibility of obtaining these modified dispersion relations within the Effective Field Theory of Inflation.


title: On positive local combinatorial dividing-lines in model theory

abstract: We introduce the notion of positive local combinatorial dividing-lines in model theory. We show these are equivalently characterized by indecomposable algebraically trivial Fraisse classes and by complete prime filter classes. We exhibit the relationship between this and collapse-of-indiscernibles dividing-lines. We examine several test cases, including those arising from various classes of hypergraphs.


title: An Improved Parametrization and Analysis of the EXP3++ Algorithm for Stochastic and Adversarial Bandits

abstract: We present a new strategy for gap estimation in randomized algorithms for multiarmed bandits and combine it with the EXP3++ algorithm of Seldin and Slivkins (2014). In the stochastic regime the strategy reduces dependence of regret on a time horizon from $(\ln t)^3$ to $(\ln t)^2$ and eliminates an additive factor of order $\Delta e^{1/\Delta^2}$, where $\Delta$ is the minimal gap of a problem instance. In the adversarial regime regret guarantee remains unchanged.


title: Temperature Dependence of Magnetic Excitations: Terahertz Magnons above the Curie Temperature

abstract: When an ordered spin system of a given dimensionality undergoes a second order phase transition the dependence of the order parameter i.e. magnetization on temperature can be well-described by thermal excitations of elementary collective spin excitations (magnons). However, the behavior of magnons themselves, as a function of temperature and across the transition temperature TC, is an unknown issue. Utilizing spin-polarized high resolution electron energy loss spectroscopy we monitor the high-energy (terahertz) magnons, excited in an ultrathin ferromagnet, as a function of temperature. We show that the magnons' energy and lifetime decrease with temperature. The temperature-induced renormalization of the magnons' energy and lifetime depends on the wave vector. We provide quantitative results on the temperature-induced damping and discuss the possible mechanism e.g., multi-magnon scattering. A careful investigation of physical quantities determining the magnons' propagation indicates that terahertz magnons sustain their propagating character even at temperatures far above TC.


title: On feebly compact shift-continuous topologies on the semilattice $\exp_n\lambda$

abstract: We study feebly compact topologies $\tau$ on the semilattice $\left(\exp_n\lambda,\cap\right)$ such that $\left(\exp_n\lambda,\tau\right)$ is a semitopological semilattice and prove that for any shift-continuous $T_1$-topology $\tau$ on $\exp_n\lambda$ the following conditions are equivalent: $(i)$~$\tau$ is countably pracompact; $(ii)$ $\tau$ is feebly compact; $(iii)$ $\tau$ is $d$-feebly compact; $(iv)$ $\left(\exp_n\lambda,\tau\right)$ is an $H$-closed space.


title: An Attention-Based Deep Net for Learning to Rank

abstract: In information retrieval, learning to rank constructs a machine-based ranking model which given a query, sorts the search results by their degree of relevance or importance to the query. Neural networks have been successfully applied to this problem, and in this paper, we propose an attention-based deep neural network which better incorporates different embeddings of the queries and search results with an attention-based mechanism. This model also applies a decoder mechanism to learn the ranks of the search results in a listwise fashion. The embeddings are trained with convolutional neural networks or the word2vec model. We demonstrate the performance of this model with image retrieval and text querying data sets.


title: Moduli of weighted stable elliptic surfaces and invariance of log plurigenera

abstract: This is the third paper in a series of work on weighted stable elliptic surfaces - elliptic fibrations with section and marked fibers weighted between zero and one. Motivated by Hassett's weighted pointed stable curves, we use the log minimal model program to construct compact moduli spaces parameterizing these objects. Moreoever, we show that the domain of weights admits a wall and chamber structure, we describe the induced wall crossing morphisms on the moduli spaces as the weight vector varies, and describe the surfaces that appear on the boundary of the moduli space. The main technical result is a proof of invariance of log plurigenera for slc elliptic surface pairs with arbitrary weights.


title: General Lossless Spatial Polarization Transformations

abstract: Liquid crystals allow for the real-time control of the polarization of light. We describe and provide some experimental examples of the types of general polarization transformations, including universal polarization transformations, that can be accomplished with liquid crystals in tandem with fixed waveplates. Implementing these transformations with an array of liquid crystals, e.g., a spatial light modulator, allows for the manipulation of the polarization across a beam's transverse plane. We outline applications of such general spatial polarization transformations in the generation of exotic types of vector polarized beams, a polarization magnifier, and the correction of polarization aberrations in light fields.


title: How Anomalous Resistivity Accelerates Magnetic Reconnection

abstract: Whether turbulence induced anomalous resistivity (AR) can facilitate a fast magnetic reconnection in collisionless plasma is a subject of active debate for decades. Recent space observations suggest that the reconnection rate can be higher than the Hall-reconnection rate and turbulent dissipation is required. In this paper, using particle-in-cell simulations, we present a case study of how AR produced by Buneman instability accelerates magnetic reconnection. We first show that the AR/drag produced by Buneman instability in a thin electron current layer 1) can dissipate magnetic energy stored in the current layer through dissipation of the kinetic energy of electron beams; 2) The inhomogeneous drag caused by wave couplings spontaneously breaks the magnetic field lines and causes impulsive fast non-Hall magnetic reconnection on electron-scales with a mean rate reaching 0.6 $V_A$. We then show that a Buneman instability driven by intense electron beams around the x-point in a 3D magnetic reconnection significantly enhances the dissipation of the magnetic energy. Electron-scale magnetic reconnections driven by the inhomogeneous drag around the x-line enhances the reconnection electric field and the in-plane perpendicular magnetic field. About 40\% of the released magnetic energy is converted into the electron thermal energy by AR while 50\% is converted into the kinetic energy of the electron beams through the acceleration by the reconnection electric field. The enhanced magnetic energy dissipation is balanced by a net Poynting flux in-flow. About 10\% of the released magnetic energy is brought out by an enhanced Poynting flux out-flow. These results suggest that AR with sufficient intensity and electron-scale inhomogeneity can significantly accelerate magnetic reconnection.

categories: physics.space-ph astro-ph.EP astro-ph.SR physics.plasm-ph


title: Density Independent Algorithms for Sparsifying $k$-Step Random Walks

abstract: We give faster algorithms for producing sparse approximations of the transition matrices of $k$-step random walks on undirected, weighted graphs. These transition matrices also form graphs, and arise as intermediate objects in a variety of graph algorithms. Our improvements are based on a better understanding of processes that sample such walks, as well as tighter bounds on key weights underlying these sampling processes. On a graph with $n$ vertices and $m$ edges, our algorithm produces a graph with about $n\log{n}$ edges that approximates the $k$-step random walk graph in about $m + n \log^4{n}$ time. In order to obtain this runtime bound, we also revisit "density independent" algorithms for sparsifying graphs whose runtime overhead is expressed only in terms of the number of vertices.


title: Antenna Count for Massive MIMO: 1.9 GHz versus 60 GHz

abstract: If we assume line-of-sight propagation and perfect channel state information at the base station -- consistent with slow moving terminals -- then a direct performance comparison between single-cell Massive MIMO at PCS and mmWave frequency bands is straightforward and highly illuminating. Line-of-sight propagation is considered favorable for mmWave because of minimal attenuation, and its facilitation of hybrid beamforming to reduce the required number of active transceivers. We quantify the number of mmWave (60 GHz) service antennas that are needed to duplicate the performance of a specified number of PCS (1.9 GHz) service antennas. As a baseline we consider a modest PCS deployment of 128 antennas serving 18 terminals. We find that, to achieve the same per-terminal max-min 95%-likely downlink throughput, 10000 mmWave antennas are needed. To match the total antenna area of the PCS array would require 128000 half-wavelength mmWave antennas, but a much reduced number is adequate because the large number of antennas also confers greater channel orthogonality. The principal alleged benefit of mmWave technology--vast amounts of inexpensive spectrum--is at least partially offset by the complexity of possibly unwieldy amounts of hardware.


title: A general framework for path convexities

abstract: In this work we deal with the so-called path convexities, defined over special collections of paths. For example, the collection of the shortest paths in a graph is associated with the well-known geodesic convexity, while the collection of the induced paths is associated with the monophonic convexity; and there are many other examples. Besides reviewing the path convexities in the literature, we propose a general path convexity framework, of which most existing path convexities can be viewed as particular cases. Some benefits of the proposed framework are the systematization of the algorithmic study of related problems and the possibility of defining new convexities not yet investigated.


title: Modeling of power distribution systems with solar generation: A case study

abstract: This document presents the data for a single-phase distribution bus (based on IEEE 37 bus) together with the model and data for a PV inverter and active and reactive power loads.


title: Quark-hadron Phase Transition in Proto-Neutron Stars Cores based on a Non-local NJL Model

abstract: We study the QCD phase diagram using a non-local SU(3) NJL model with vector interactions among quarks. We analyze several thermodynamic quantities such as entropy and specific heat, and study the influence of vector interactions on the thermodynamic properties of quark matter. Upon imposing electric charge neutrality and baryon number conservation on the field equations, we compute models for the equation of state of the inner cores of proto-neutron stars providing a non-local treatment of quark matter for astrophysics.


title: Critical exponents of normal subgroups, the spectrum of group extended transfer operators, and Kazhdan distance

abstract: For a pinched Hadamard manifold $X$ and a discrete group of isometries $\Gamma$ of $X$, the critical exponent $\delta_\Gamma$ is the exponential growth rate of the orbit of a point in $X$ under the action of $\Gamma$. We show that the critical exponent for any family $\mathcal{N}$ of normal subgroups of $\Gamma_0$ has the same coarse behaviour as the Kazhdan distances for the right regular representations of the quotients $\Gamma_0/\Gamma$. The key tool is to analyse the spectrum of transfer operators associated to subshifts of finite type, for which we obtain a result of independent interest.


title: Galaxies grow their bulges and black holes in diverse ways

abstract: Galaxies with Milky Way-like stellar masses have a wide range of bulge and black hole masses; in turn, these correlate with other properties such as star formation history. While many processes may drive bulge formation, major and minor mergers are expected to play a crucial role. Stellar halos offer a novel and robust measurement of galactic merger history; cosmologically-motivated models predict that mergers with larger satellites produce more massive, higher metallicity stellar halos, reproducing the recently-observed stellar halo metallicity-mass relation. We quantify the relationship between stellar halo mass and bulge or black hole prominence using a sample of eighteen Milky Way-mass galaxies with newly-available measurements of (or limits on) stellar halo properties. There is an order of magnitude range in bulge mass, and two orders of magnitude in black hole mass, at a given stellar halo mass (or, equivalently, merger history). Galaxies with low mass bulges show a wide range of quiet merger histories, implying formation mechanisms that do not require intense merging activity. Galaxies with massive 'classical' bulges and central black holes also show a wide range of merger histories. While three of these galaxies have massive stellar halos consistent with a merger origin, two do not - merging appears to have had little impact in making these two massive 'classical' bulges. Such galaxies may be ideal laboratories to study massive bulge formation through pathways such as early gas-rich accretion, violent disk instabilities or misaligned infall of gas throughout cosmic time.


title: Evaluating Reactor Antineutrino Signals for WATCHMAN

abstract: Increasing the distance from which an antineutrino detector is capable of monitoring the operation of a registered reactor, or discovering a clandestine reactor, strengthens the Non-Proliferation of Nuclear Weapons Treaty. This paper presents calculations of reactor antineutrino interactions from quasi-elastic neutrino-proton scattering and elastic neutrino-electron scattering in a water-based detector operated $\gtrsim10$ km from a commercial power reactor. It separately calculates signal from the proximal reactor and background from all other registered reactors. The main results are differential and integral interaction rates from the quasi-elastic and elastic processes. There are two underground facilities capable of hosting a detector ($\sim1$ kT H$_2$O) project nearby ($L\sim20$ km) an operating commercial reactor ($P_{th}\sim3$ GW). These reactor-site combinations, which are under consideration for project WATCHMAN, are Perry-Morton on the southern shore of Lake Erie in the United States and Hartlepool-Boulby on the western shore of the North Sea in England. The signal rate from the proximal reactor is about five times greater at the Morton site than at the Boulby site due to shorter reactor-site separation distance, larger reactor thermal power, and greater neutrino oscillation survival probability. Although the background rate from all other reactors is larger at Morton than at Boulby, it is a smaller fraction of the signal rate from the proximal reactor at Morton than at Boulby. Moreover, the Hartlepool power plant has two cores whereas the Perry plant has a single core. The Boulby site, therefore, offers an opportunity for remotely monitoring the on/off cycle of a reactor core under more stringent conditions than does the Morton site.


title: Derivative Based Focal Plane Array Nonuniformity Correction

abstract: This paper presents a fast and robust method for fixed pattern noise nonuniformity correction of infrared focal plane arrays. The proposed method requires neither shutter nor elaborate calibrations and therefore enables a real time correction with no interruptions. Based on derivative estimation of the fixed pattern noise from pixel sized translations of the focal plane array, the proposed method has the advantages of being invariant to the noise magnitude and robust to unknown camera and inter-scene movements while being virtually transparent to the end-user.


title: Shannon-inspired Statistical Computing to Enable Spintronics

abstract: Modern computing systems based on the von Neumann architecture are built from silicon complementary metal oxide semiconductor (CMOS) transistors that need to operate under practically error free conditions with 1 error in $10^{15}$ switching events. The physical dimensions of CMOS transistors have scaled down over the past five decades leading to exponential increases in functional density and energy consumption. Today, the energy and delay reductions from scaling have stagnated, motivating the search for a CMOS replacement. Of these, spintronics offers a path for enhancing the functional density and scaling the energy down to fundamental thermodynamic limits of 100kT to 1000kT. However, spintronic devices exhibit high error rates of 1 in 10 or more when operating at these limits, rendering them incompatible with deterministic nature of the von Neumann architecture. We show that a Shannon-inspired statistical computing framework can be leveraged to design a computer made from such stochastic spintronic logic gates to provide a computational accuracy close to that of a deterministic computer. This extraordinary result allowing a $10^{13}$ fold relaxation in acceptable error rates is obtained by engineering the error distribution coupled with statistical error compensation.


title: On the Consistency of $k$-means++ algorithm

abstract: We prove in this paper that the expected value of the objective function of the $k$-means++ algorithm for samples converges to population expected value. As $k$-means++, for samples, provides with constant factor approximation for $k$-means objectives, such an approximation can be achieved for the population with increase of the sample size. This result is of potential practical relevance when one is considering using subsampling when clustering large data sets (large data bases).


title: Reconstructing binary matrices under window constraints from their row and column sums

abstract: The present paper deals with the discrete inverse problem of reconstructing binary matrices from their row and column sums under additional constraints on the number and pattern of entries in specified minors. While the classical consistency and reconstruction problems for two directions in discrete tomography can be solved in polynomial time, it turns out that these window constraints cause various unexpected complexity jumps back and forth from polynomial-time solvability to $\mathbb{N}\mathbb{P}$-hardness.


title: METAQ: Bundle Supercomputing Tasks

abstract: We describe a light-weight system of bash scripts for efficiently bundling supercomputing tasks into large jobs, so that one can take advantage of incentives or discounts for requesting large allocations. The software can backfill computational tasks, avoiding wasted cycles, and can streamline collaboration between different users. It is simple to use, functioning similarly to batch systems like PBS, MOAB, and SLURM.


title: Physical properties of the first quasars

abstract: Since the beginning of the new millennium, more than 100 $z\sim 6$ quasars have been discovered through several surveys and followed-up with multi-wavelength observations. These data provided a large amount of information on the growth of supermassive black holes at the early epochs, the properties of quasar host galaxies and the joint formation and evolution of these massive systems. We review the properties of the highest-$z$ quasars known so far, especially focusing on some of the most recent results obtained in (sub-)millimeter bands. We discuss key observational challenges and open issues in theoretical models and highlight possible new strategies to improve our understanding of the galaxy-black hole formation and evolution in the early Universe.


title: Gravitational wave, collider and dark matter signals from a scalar singlet electroweak baryogenesis

abstract: We analyse a simple extension of the SM with just an additional scalar singlet coupled to the Higgs boson. We discuss the possible probes for electroweak baryogenesis in this model including collider searches, gravitational wave and direct dark matter detection signals. We show that a large portion of the model parameter space exists where the observation of gravitational waves would allow detection while the indirect collider searches would not.


title: Anisotropic Andreev Reflection and Josephson Effect in Ballistic Phosphorene

abstract: We study Andreev reflection and the Josephson effect in a ballistic monolayer of black phosphorous, known as phosphorene. Due to the anisotropic band structure of this system, the supercurrent changes with an order of magnitude when comparing tunneling along two perpendicular directions in the monolayer. We show that the main reason for this effect is a large difference in the number of transverse modes in Andreev bound states. The oscillatory behavior of the supercurrent as a function of the length and chemical potential of the junction also differs substantially depending on the orientation of the superconducting electrodes deposited on the phosphorene sheet. For Andreev reflection, we show that a gate voltage controls the probability of this process and that the anisotropic behavior found in the supercurrent case is also present for conductance spectra.


title: The mass, colour, and structural evolution of today's massive galaxies since z~5

abstract: In this paper, we use stacking analysis to trace the mass-growth, colour evolution, and structural evolution of present-day massive galaxies ($\log(M_{*}/M_{\odot})=11.5$) out to $z=5$. We utilize the exceptional depth and area of the latest UltraVISTA data release, combined with the depth and unparalleled seeing of CANDELS to gather a large, mass-selected sample of galaxies in the NIR (rest-frame optical to UV). Progenitors of present-day massive galaxies are identified via an evolving cumulative number density selection, which accounts for the effects of merging to correct for the systematic biases introduced using a fixed cumulative number density selection, and find progenitors grow in stellar mass by $\approx1.5~\mathrm{dex}$ since $z=5$. Using stacking, we analyze the structural parameters of the progenitors and find that most of the stellar mass content in the central regions was in place by $z\sim2$, and while galaxies continue to assemble mass at all radii, the outskirts experience the largest fractional increase in stellar mass. However, we find evidence of significant stellar mass build up at $r<3~\mathrm{kpc}$ beyond $z>4$ probing an era of significant mass assembly in the interiors of present day massive galaxies. We also compare mass assembly from progenitors in this study to the EAGLE simulation and find qualitatively similar assembly with $z$ at $r<3~\mathrm{kpc}$. We identify $z\sim1.5$ as a distinct epoch in the evolution of massive galaxies where progenitors transitioned from growing in mass and size primarily through in-situ star formation in disks to a period of efficient growth in $r_{e}$ consistent with the minor merger scenario.


title: Rainbow scattering in the gravitational field of a compact object

abstract: We study the elastic scattering of a planar wave in the curved spacetime of a compact object such as a neutron star, via a heuristic model: a scalar field impinging upon a spherically-symmetric uniform density star of radius $R$ and mass $M$. For $R < r_c$, there is a divergence in the deflection function at the light-ring radius $r_c = 3GM/c^2$, which leads to spiral scattering (orbiting) and a backward glory; whereas for $R > r_c$ there instead arises a stationary point in the deflection function which creates a caustic and rainbow scattering. As in nuclear rainbow scattering, there is an Airy-type oscillation on a Rutherford-like cross section, followed by a shadow zone. We show that, for $R \sim 3.5 GM/c^2$, the rainbow angle lies close to $180^\circ$, and thus there arises enhanced back-scattering and glory. We explore possible implications for gravitational wave astronomy, and dark matter models.


title: Geodesic Diagrams, Gravitational Interactions & OPE Structures

abstract: We give a systematic procedure to evaluate conformal partial waves involving symmetric tensors for an arbitrary CFT$_d$ using geodesic Witten diagrams in AdS$_{d+1}$. Using this procedure we discuss how to draw a line between the tensor structures in the CFT and cubic interactions in AdS. We contrast this map to known results using three-point Witten diagrams: the maps obtained via volume versus geodesic integrals differ. Despite these differences, we show how to decompose four-point exchange Witten diagrams in terms of geodesic diagrams, and we discuss the product expansion of local bulk fields in AdS.


title: Local Group Ultra-Faint Dwarf Galaxies in the Reionization Era

abstract: Motivated by the stellar fossil record of Local Group (LG) dwarf galaxies, we show that the star-forming ancestors of the faintest ultra-faint dwarf galaxies (UFDs; ${\rm M}_{\rm V}$ $\sim -2$ or ${\rm M}_{\star}$ $\sim 10^{2}$ at $z=0$) had ultra-violet (UV) luminosities of ${\rm M}_{\rm UV}$ $\sim -3$ to $-6$ during reionization ($z\sim6-10$). The existence of such faint galaxies has substantial implications for early epochs of galaxy formation and reionization. If the faint-end slopes of the UV luminosity functions (UVLFs) during reionization are steep ($\alpha\lesssim-2$) to ${\rm M}_{\rm UV}$ $\sim -3$, then: (i) the ancestors of UFDs produced $>50$% of UV flux from galaxies; (ii) galaxies can maintain reionization with escape fractions that are $>$2 times lower than currently-adopted values; (iii) direct HST and JWST observations may detect only $\sim10-50$% of the UV light from galaxies; (iv) the cosmic star formation history increases by $\gtrsim4-6$ at $z\gtrsim6$. Significant flux from UFDs, and resultant tensions with LG dwarf galaxy counts, are reduced if the high-redshift UVLF turns over. Independent of the UVLF shape, the existence of a large population of UFDs requires a non-zero luminosity function to ${\rm M}_{\rm UV}$ $\sim -3$ during reionization.


title: The Weak Gravity Conjecture and the Axionic Black Hole Paradox

abstract: In theories with a perturbatively massless 2-form (dual to an axion), a paradox may arise in the process of black hole evaporation. Schwarzschild black holes can support a non-trivial Wilson-line-type field, the integral of the 2-form around their horizon. After such an 'axionic black hole' evaporates, the Wilson line must be supported by the corresponding 3-form field strength in the region formerly occupied by the black hole. In the limit of small axion decay-constant f, the energy required for this field configuration is too large. The natural resolution is through the presence of light strings, which allow the black hole to "shed" its axionic hair sufficiently early. This gives rise to a new Weak-Gravity-type argument in the 2-form context: Small coupling, in this case f, enforces the presence of light strings or a low cutoff. We also discuss how this argument may be modified in situations where the weak coupling regime is achieved in the low-energy effective theory through an appropriate gauging of a model with a vector field and two 2-forms.


title: Multiwavelength follow-up of a rare IceCube neutrino multiplet

abstract: On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.


title: Fundamental Constant Observational Bounds on the Variability of the QCD Scale

abstract: Many physical theories beyond the Standard Model predict time variations of basic physics parameters. Direct measurement of the time variations of these parameters is very difficult or impossible to achieve. By contrast, measurements of fundamental constants are relatively easy to achieve, both in the laboratory and by astronomical spectra of atoms and molecules in the early universe. In this work measurements of the proton to electron mass ratio $\mu$ and the fine structure constant $\alpha$ are combined to place mildly model dependent limits on the fractional variation of the Quantum Chromodynamic Scale and the sum of the fractional variations of the Higgs Vacuum Expectation Value and the Yukawa couplings on time scales of more than half the age of the universe. The addition of another model parameter allows the fractional variation of the Higgs VEV and the Yukawa couplings to be computed separately. Limits on their variation are found at the level of less than $5 \times 10^{-5}$ over the past seven gigayears. A model dependent relation between the expected fractional variation of $\alpha$ relative to $\mu$ tightens the limits to $10^{-7}$ over the same time span. Limits on the present day rate of change of the constants and parameters are then calculated using slow roll quintessence. A primary result of this work is that studies of the dimensionless fundamental constants such as $\alpha$ and $\mu$, whose values depend on the values of the physics parameters, are excellent monitors of the limits on the time variation of these parameters.


title: Moving mesh simulations of star forming cores in magneto-gravo-turbulence

abstract: Star formation in our Galaxy occurs in molecular clouds that are self-gravitating, highly turbulent, and magnetized. We study the conditions under which cloud cores inherit large-scale magnetic field morphologies and how the field is governed by cloud turbulence. We present four moving-mesh simulations of supersonic, turbulent, isothermal, self-gravitating gas with a range of magnetic mean-field strengths characterized by the Alfv\'enic Mach number $\mathcal{M}_{{\rm A}, 0}$, resolving pre-stellar core formation from parsec to a few AU scales. In our simulations with the turbulent kinetic energy density dominating over magnetic pressure ($\mathcal{M}_{{\rm A}, 0}>1$), we find that the collapse is approximately isotropic with $B\propto\rho^{2/3}$, core properties are similar regardless of initial mean-field strength, and the field direction on $100$ AU scales is uncorrelated with the mean field. However, in the case of a dominant large-scale magnetic field ($\mathcal{M}_{{\rm A}, 0}=0.35$), the collapse is anisotropic with $B\propto\rho^{1/2}$. This transition at $\mathcal{M}_{{\rm A}, 0}\sim1$ is not expected to be sharp, but clearly signifies two different paths for magnetic field evolution in star formation. Based on observations of different star forming regions, we conclude that star formation in the interstellar medium may occur in both regimes. Magnetic field correlation with the mean-field extends to smaller scales as $\mathcal{M}_{{\rm A}, 0}$ decreases, making future ALMA observations useful for constraining $\mathcal{M}_{{\rm A}, 0}$ of the interstellar medium.


title: Positivity Bounds for Scalar Theories

abstract: Assuming the existence of a local, analytic, unitary UV completion in a Poincar\'{e} invariant scalar field theory with a mass gap, we derive an infinite number of positivity requirements using the known properties of the amplitude at and away from the forward scattering limit. These take the form of bounds on combinations of the pole subtracted scattering amplitude and its derivatives. In turn, these positivity requirements act as constraints on the operator coefficients in the low energy effective theory. For certain theories these constraints can be used to place an upper bound on the mass of the next lightest state that must lie beyond the low energy effective theory if such a UV completion is to ever exist.


title: Enabling Multi-Source Neural Machine Translation By Concatenating Source Sentences In Multiple Languages

abstract: In this paper, we explore a simple solution to "Multi-Source Neural Machine Translation" (MSNMT) which only relies on preprocessing a N-way multilingual corpus without modifying the Neural Machine Translation (NMT) architecture or training procedure. We simply concatenate the source sentences to form a single long multi-source input sentence while keeping the target side sentence as it is and train an NMT system using this preprocessed corpus. We evaluate our method in resource poor as well as resource rich settings and show its effectiveness (up to 4 BLEU using 2 source languages and up to 6 BLEU using 5 source languages). We also compare against existing methods for MSNMT and show that our solution gives competitive results despite its simplicity. We also provide some insights on how the NMT system leverages multilingual information in such a scenario by visualizing attention.


title: The atoms of spacetime and the cosmological constant

abstract: I describe an approach which relates classical gravity to the quantum microstructure of spacetime. In this approach, the field equations arise from maximizing the density of states of the matter plus geometry. The former is identified using the thermodynamics of null surfaces. The latter arises from the existence of a zero-point length in the spacetime which associates an internal degree of freedom with each event in the spacetime, in the form of a fluctuating vector of constant norm. The density of states, as well as the resulting gravitational field equations, remain invariant when a constant is added to the matter Lagrangian, which is a symmetry of the matter sector. The cosmological constant arises as an integration constant and renders the amount of cosmic information accessible to an eternal observer finite. The relation between the cosmological constant and the cosmic information allows us to determine the numerical value of cosmological constant from the information content of the quantum spacetime, within the context of cosmology.


title: T-Branes at the Limits of Geometry

abstract: Singular limits of 6D F-theory compactifications are often captured by T-branes, namely a non-abelian configuration of intersecting 7-branes with a nilpotent matrix of normal deformations. The long distance approximation of such 7-branes is a Hitchin-like system in which simple and irregular poles emerge at marked points of the geometry. When multiple matter fields localize at the same point in the geometry, the associated Higgs field can exhibit irregular behavior, namely poles of order greater than one. This provides a geometric mechanism to engineer wild Higgs bundles. Physical constraints such as anomaly cancellation and consistent coupling to gravity also limit the order of such poles. Using this geometric formulation, we unify seemingly different wild Hitchin systems in a single framework in which orders of poles become adjustable parameters dictated by tuning gauge singlet moduli of the F-theory model.


title: The Co-Evolution of a Magnetized Intracluster Medium and Hot Galactic Coronae: Magnetic Field Amplification and Turbulence Generation

abstract: We use adaptive-mesh magnetohydrodynamic simulations to study the effect of magnetic fields on ram pressure stripping of galaxies in the intracluster medium (ICM). Although the magnetic pressure in typical clusters is not strong enough to affect the gas mass loss rate from galaxies, magnetic fields can affect the morphology of stripped galaxies. ICM magnetic fields are draped around orbiting galaxies and aligned with their stripped tails. Magnetic fields suppress shear instabilities at the galaxy-ICM interface, and magnetized tails are smoother and narrower than tails in comparable hydrodynamic simulations in Vijayaraghavan & Ricker (2015). Orbiting galaxies stretch and amplify ICM magnetic fields, amplifying magnetic power spectra on $10 - 100$ kpc scales. Galaxies inject turbulent kinetic energy into the ICM via their turbulent wakes and $g$-waves. The magnetic energy and kinetic energy in the ICM increase up to $1.5 - 2$ Gyr of evolution, after which galaxies are stripped of most of their gas, and do not have sufficiently large gaseous cross sections to further amplify magnetic fields and inject turbulent kinetic energy. The increase in turbulent pressure due to galaxy stripping and generation of $g$-waves results in an increase in the turbulent volume fraction of the ICM. This turbulent kinetic energy is not a significant contributor to the overall ICM energy budget, but greatly impacts the evolution of the ICM magnetic field. Additionally, the effect of galaxies on magnetic fields can potentially be observed in high resolution Faraday rotation measure (RM) maps as small scale fluctuations in the RM structure.


title: Spinning Geodesic Witten Diagrams

abstract: We present an expression for the four-point conformal blocks of symmetric traceless operators of arbitrary spin as an integral over a pair of geodesics in Anti-de Sitter space, generalizing the geodesic Witten diagram formalism of Hijano et al [arXiv:1508.00501] to arbitrary spin. As an intermediate step in the derivation, we identify a convenient basis of bulk three-point interaction vertices which give rise to all possible boundary three point structures. We highlight a direct connection between the representation of the conformal block as a geodesic Witten diagram and the shadow operator formalism.


title: Hourglass semimetals with nonsymmorphic symmetries in three dimensions

abstract: It was recently shown that nonsymmorphic space group symmetries can protect novel surface states with hourglass-like dispersions. In this paper, we show that such dispersions can also appear in the bulk of three-dimensional (3D) systems which respect nonsymmorphic symmetries. Specifically, we construct 3D lattice models featuring hourglass-like dispersions in the bulk, which are protected by nonsymmorphic and time-reversal symmetries. We call such systems hourglass semimetals, as they have point or line nodes associated with hourglass-like dispersions. Hourglass nodal lines appear in glide-invariant planes, while hourglass Weyl points can occur on screw-invariant axes. The Weyl points and surface Fermi arcs in hourglass Weyl semimetals are stable against weak perturbations breaking those nonsymmorphic symmetries. Our results may shed light on searching for exotic Weyl semimetals in nonsymmorphic materials.


title: Low-control and robust quantum refrigerator and applications with electronic spins in diamond

abstract: We propose a general protocol for low-control refrigeration and thermometry of thermal qubits, which can be implemented using electronic spins in diamond. The refrigeration is implemented by a probe, consisting of a network of interacting spins. The protocol involves two operations: (i) free evolution of the probe; and (ii) a swap gate between one spin in the probe and the thermal qubit we wish to cool. We show that if the initial state of the probe falls within a suitable range, and the free evolution of the probe is both unital and conserves the excitation in the $z$-direction, then the cooling protocol will always succeed, with an efficiency that depends on the rate of spin dephasing and the swap gate fidelity. Furthermore, measuring the probe after it has cooled many qubits provides an estimate of their temperature. We provide a specific example where the probe is a Heisenberg spin chain, and suggest a physical implementation using electronic spins in diamond. Here the probe is constituted of nitrogen vacancy (NV) centers, while the thermal qubits are dark spins. By using a novel pulse sequence, a chain of NV centers can be made to evolve according to a Heisenberg Hamiltonian. This proposal allows for a range of applications, such as NV-based nuclear magnetic resonance of photosensitive molecules kept in a dark spot on a sample, and it opens up possibilities for the study of quantum thermodynamics, environment-assisted sensing, and many-body physics.


title: Locality from the Spectrum

abstract: Essential to the description of a quantum system are its local degrees of freedom, which enable the interpretation of subsystems and dynamics in the Hilbert space. While a choice of local tensor factorization of the Hilbert space is often implicit in the writing of a Hamiltonian or Lagrangian, the identification of local tensor factors is not intrinsic to the Hilbert space itself. Instead, the only basis-invariant data of a Hamiltonian is its spectrum, which does not manifestly determine the local structure. This ambiguity is highlighted by the existence of dualities, in which the same energy spectrum may describe two systems with very different local degrees of freedom. We argue that in fact, the energy spectrum alone almost always encodes a unique description of local degrees of freedom when such a description exists, allowing one to explicitly identify local subsystems and how they interact. As a consequence, we can almost always write a Hamiltonian in its local presentation given only its spectrum. In special cases, multiple dual local descriptions can be extracted from a given spectrum, but generically the local description is unique.


title: Considerations on how to investigate planes of satellite galaxies

abstract: The existence of a spatially thin, kinematically coherent Disk of Satellites (DoS) around the Milky Way (MW), is a problem that often garners vivacious debate in the literature or at scientific meetings. One of the most recent incarnations of this wrangle occurred with two papers by Maji et al, who argued that these structures "maybe a misinterpretation of the data". These claims are in stark contrast to previous works. Motivated by this and other recent publications on this problem, we discuss necessary considerations to make, observational effects to consider, and pitfalls to avoid when investigating satellite galaxy planes such as the MW's DoS. In particular, we emphasize that conclusions need to have a statistical basis including a determination of the significance of satellite alignments, observational biases must not be ignored, and measurement errors (e.g. for proper motions) need to be considered. We discuss general problems faced by attempts to determine the dynamical stability of the DoS via orbit integrations of MW satellite galaxies, and demonstrate that to interpret simulations, it is helpful to compare them with a null case of isotropically distributed satellite positions and velocities. Based on these criteria, we find that the conclusions of Maji et al. do not hold up to scrutiny, and that their hydrodynamic cosmological simulation of a single host shows no evidence for a significant kinematic coherence among the simulated satellite galaxies, in contrast to the observed MW system.


title: An Inequality for the Correlation of Two Functions Operating on Symmetric Bivariate Normal Variables

abstract: An inequality is derived for the correlation of two univariate functions operating on symmetric bivariate normal random variables. The inequality is a simple consequence of the Cauchy-Schwarz inequality.


title: Search for right-handed neutrinos from dark matter annihilation with gamma-rays

abstract: Several extensions of the Standard Model contain right-handed (sterile) neutrinos in the GeV-TeV mass range. Due to their mixing with the active neutrinos, they may give rise to novel effects in cosmology, neutrino physics, and collider searches. In addition, right-handed neutrinos can also appear as final states from dark matter annihilations, with important implications for dark matter indirect detection searches. In this paper, we use current data from the Fermi Large Area Telescope (6-year observation of dwarf spheroidal galaxies) and H.E.S.S. (10-year observation of the Galactic center) to constrain the annihilation of dark matter into right-handed neutrinos. We consider right-handed neutrino with masses between 10 GeV and 1 TeV, including both two-body and three-body decays, to derive bounds on the dark matter annihilation rate, $\langle \sigma v\rangle$, as a function of the dark matter mass. Our results show, in particular, that the thermal dark matter annihilation cross section, $3\times 10^{-26}\, {\rm cm^3 s^{-1}} $, into right-handed neutrinos is excluded for dark matter masses smaller than $200$ GeV.


title: Bethe-Boltzmann Hydrodynamics and Spin Transport in the XXZ Chain

abstract: Quantum integrable systems, such as the interacting Bose gas in one dimension and the XXZ quantum spin chain, have an extensive number of local conserved quantities that endow them with exotic thermalization and transport properties. We discuss recently introduced hydrodynamic approaches for such integrable systems from the viewpoint of kinetic theory and extend the previous works by proposing a numerical scheme to solve the hydrodynamic equations for finite times and arbitrary locally equilibrated initial conditions. We then discuss how such methods can be applied to describe non-equilibrium steady states involving ballistic heat and spin currents. In particular, we show that the spin Drude weight in the XXZ chain, previously accessible only by rigorous techniques of limited scope or controversial thermodynamic Bethe ansatz arguments, may be evaluated from hydrodynamics in very good agreement with density-matrix renormalization group calculations.

categories: cond-mat.stat-mech cond-mat.str-el math-ph math.MP quant-ph


title: A Chern-Simons Pandemic

abstract: In this paper we study the consistency of generalized global symmetries in theories of quantum gravity, in particular string theory. Such global symmetries arise in theories with $(p+1)$-form gauge fields, and for spacetime dimension $d\leq p+3$ there are obstructions to their breaking even by quantum effects of charged objects. In 4d theories with a 2-form gauge field (or with an axion scalar), these fields endow Schwarzschild black holes with quantum hair, a global charge leading to usual trouble with remnants. We describe precise mechanisms, and examples from string compactifications and holographic pairs, in which these problems are evaded by either gauging or breaking the global symmetry, via (suitable versions of) Stuckelberg or Kaloper-Sorbo couplings. We argue that even in the absence of such couplings, the generic solution in string theory is the breaking of the global symmetries by cubic Chern-Simons terms involving different antisymmetric tensor fields. We conjecture that any theory with (standard or higher-degree antisymmetric tensor) gauge fields is in the Swampland unless its effective action includes such Chern-Simons terms. This conjecture implies that many familiar theories, like QED (even including the charged particles required by the Weak Gravity Conjecture) or $\mathcal{N}=8$ supergravity in four dimensions, are inconsistent in quantum gravity unless they are completed by these Chern-Simons terms.


title: FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation

abstract: The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.


title: The Invisible AGN Catalogue: A Mid-Infrared - Radio Selection Method for Optically-Faint Active Galactic Nuclei

abstract: A large fraction of active galactic nuclei (AGN) are "invisible" in extant optical surveys due to either distance or dust-obscuration. The existence of this large population of dust-obscured, infrared-bright AGN is predicted by models of galaxy - supermassive black hole coevolution and is required to explain the observed X-ray and infrared backgrounds. Recently, infrared colour-cuts with WISE have identified a portion of this missing population. However, as the host galaxy brightness relative to that of the AGN increases, it becomes increasingly difficult to differentiate between infrared emission originating from the AGN and from its host galaxy. As a solution, we have developed a new method to select obscured AGN using their 20 cm continuum emission to identify the objects as AGN. We created the resulting Invisible AGN Catalogue by selecting objects that are detected in AllWISE (mid-IR) and FIRST (20 cm), but are not detected in SDSS (optical) or 2MASS (near-IR), producing a final catalogue of 46,258 objects. 30 per cent of the objects are selected by existing selection methods, while the remaining 70 per cent represent a potential previously-unidentified population of candidate AGN that are missed by mid-infrared colour cuts. Additionally, by relying on a radio continuum detection, this technique is efficient at detecting radio-loud AGN at z > 0.29, regardless of their level of dust obscuration or their host galaxy's relative brightness.


title: Bijections for Dyck paths with all peak heights of the same parity

abstract: We show bijectively that Dyck paths with all peaks at odd height are counted by the Motzkin numbers and Dyck paths with all peaks at even height are counted by the Riordan numbers.


title: Developing a comprehensive framework for multimodal feature extraction

abstract: Feature extraction is a critical component of many applied data science workflows. In recent years, rapid advances in artificial intelligence and machine learning have led to an explosion of feature extraction tools and services that allow data scientists to cheaply and effectively annotate their data along a vast array of dimensions---ranging from detecting faces in images to analyzing the sentiment expressed in coherent text. Unfortunately, the proliferation of powerful feature extraction services has been mirrored by a corresponding expansion in the number of distinct interfaces to feature extraction services. In a world where nearly every new service has its own API, documentation, and/or client library, data scientists who need to combine diverse features obtained from multiple sources are often forced to write and maintain ever more elaborate feature extraction pipelines. To address this challenge, we introduce a new open-source framework for comprehensive multimodal feature extraction. Pliers is an open-source Python package that supports standardized annotation of diverse data types (video, images, audio, and text), and is expressly with both ease-of-use and extensibility in mind. Users can apply a wide range of pre-existing feature extraction tools to their data in just a few lines of Python code, and can also easily add their own custom extractors by writing modular classes. A graph-based API enables rapid development of complex feature extraction pipelines that output results in a single, standardized format. We describe the package's architecture, detail its major advantages over previous feature extraction toolboxes, and use a sample application to a large functional MRI dataset to illustrate how pliers can significantly reduce the time and effort required to construct sophisticated feature extraction workflows while increasing code clarity and maintainability.


title: Throughput Optimal Beam Alignment in Millimeter Wave Networks

abstract: Millimeter wave communications rely on narrow-beam transmissions to cope with the strong signal attenuation at these frequencies, thus demanding precise beam alignment between transmitter and receiver. The communication overhead incurred to achieve beam alignment may become a severe impairment in mobile networks. This paper addresses the problem of optimizing beam alignment acquisition, with the goal of maximizing throughput. Specifically, the algorithm jointly determines the portion of time devoted to beam alignment acquisition, as well as, within this portion of time, the optimal beam search parameters, using the framework of Markov decision processes. It is proved that a bisection search algorithm is optimal, and that it outperforms exhaustive and iterative search algorithms proposed in the literature. The duration of the beam alignment phase is optimized so as to maximize the overall throughput. The numerical results show that the throughput, optimized with respect to the duration of the beam alignment phase, achievable under the exhaustive algorithm is 88.3% lower than that achievable under the bisection algorithm. Similarly, the throughput achievable by the iterative search algorithm for a division factor of 4 and 8 is, respectively, 12.8% and 36.4% lower than that achievable by the bisection algorithm.


title: Community Detection with Colored Edges

abstract: In this paper, we prove a sharp limit on the community detection problem with colored edges. We assume two equal-sized communities and there are $m$ different types of edges. If two vertices are in the same community, the distribution of edges follows $p_i=\alpha_i\log{n}/n$ for $1\leq i \leq m$, otherwise the distribution of edges is $q_i=\beta_i\log{n}/n$ for $1\leq i \leq m$, where $\alpha_i$ and $\beta_i$ are positive constants and $n$ is the total number of vertices. Under these assumptions, a fundamental limit on community detection is characterized using the Hellinger distance between the two distributions. If $\sum_{i=1}^{m} {(\sqrt{\alpha_i} - \sqrt{\beta_i})}^2 >2$, then the community detection via maximum likelihood (ML) estimator is possible with high probability. If $\sum_{i=1}^m {(\sqrt{\alpha_i} - \sqrt{\beta_i})}^2 < 2$, the probability that the ML estimator fails to detect the communities does not go to zero.


title: Role model detection using low rank similarity matrix

abstract: Computing meaningful clusters of nodes is crucial to analyse large networks. In this paper, we apply new clustering methods to improve the computational time. We use the properties of the adjacency matrix to obtain better role extraction. We also define a new non-recursive similarity measure and compare its results with the ones obtained with Browet's similarity measure. We will show the extraction of the different roles with a linear time complexity. Finally, we test our algorithm with real data structures and analyse the limit of our algorithm.


title: The Meissner effect for weakly isolated horizons

abstract: Black holes are important astrophysical objects describing an end state of stellar evolution, which are observed frequently. There are theoretical predictions that Kerr black holes with high spins expel magnetic fields. However, Kerr black holes are pure vacuum solutions, which do not include accretion disks, and additionally previous investigations are mainly limited to weak magnetic fields. We prove for the first time in full general relativity that generic rapidly spinning black holes including those deformed by accretion disks still expel even strong magnetic fields. Analogously to a similar property of superconductors, this is called Meissner effect.


title: How Much Urban Traffic is Searching for Parking? Simulating Curbside Parking as a Network of Finite Capacity Queues

abstract: With the increasing availability of transaction data collected by digital parking meters, paid curbside parking can be advantageously modeled as a network of interdependent queues. In this article we introduce methods for analyzing a special class of networks of finite capacity queues, where tasks arrive from an exogenous source, join the queue if there is an available server or are rejected and move to another queue in search of service according to the network topology. Such networks can be useful for modeling curbside parking since queues in the network perform the same function and drivers searching for an available server are under combinatorial constraints and jockeying is not instantaneous. Further, we provide a motivating example for such networks of finite capacity queues in the context of drivers searching for parking in the neighborhood of Belltown in Seattle, Washington, USA. Lastly, since the stationary distribution of such networks used to model parking are difficult to satisfactorily characterize, we also introduce a simulation tool for the purpose of testing the assumptions made to estimate interesting performance metrics. Our results suggest that a Markovian relaxation of the problem when solving for the mean rate metrics is comparable to deterministic service times reflective of a driver's tendency to park for the maximum allowable time.


title: Robust Phase Retrieval via ADMM with Outliers

abstract: An outlier-resistance phase retrieval algorithm based on alternating direction method of multipliers (ADMM) is devised in this letter. Instead of the widely used least squares criterion that is only optimal for Gaussian noise environment, we adopt the least absolute deviation criterion to enhance the robustness against outliers. Considering both intensity- and amplitude-based observation models, the framework of ADMM is developed to solve the resulting non-differentiable optimization problems. It is demonstrated that the core subproblem of ADMM is the proximity operator of the L1-norm, which can be computed efficiently by soft-thresholding in each iteration. Simulation results are provided to validate the accuracy and efficiency of the proposed approach compared to the existing schemes.


title: Con los juegos tambi\'en se educa: un enfoque educativo de los juegos de la Oca y el Parch\'is

abstract: In this work, we have designed and implemented, based on traditional board games such as the game of the Goose or Parchis, an educational software that aim to reinforce the learning of children. The idea that we are going to develop to do this is very simple: the children play with the same rules as in the traditional game but we add the functionality that after throwing the dice the system asks a question randomly chosen from a predefined database that can be easily modified, and the game piece only moves in the case the question is answered correctly.


title: DEEProtect: Enabling Inference-based Access Control on Mobile Sensing Applications

abstract: Personal sensory data is used by context-aware mobile applications to provide utility. However, the same data can also be used by an adversary to make sensitive inferences about a user thereby violating her privacy. We present DEEProtect, a framework that enables a novel form of inference control, in which mobile apps with access to sensor data are limited (provably) in their ability to make inferences about user's sensitive data and behavior. DEEProtect adopts a two-layered privacy strategy. First, it leverages novel autoencoder techniques to perform data minimization and limits the amount of information being shared; the learning network is used to derive a compact representation of sensor data consisting only of features relevant to authorized utility-providing inferences. Second, DEEProtect obfuscates the previously learnt features, thereby providing an additional layer of protection against sensitive inferences. Our framework supports both conventional as well as a novel relaxed notion of local differential privacy that enhances utility. Through theoretical analysis and extensive experiments using real-world datasets, we demonstrate that when compared to existing approaches DEEProtect provides provable privacy guarantees with up to 8x improvement in utility. Finally, DEEProtect shares obfuscated but raw sensor data reconstructed from the perturbed features, thus requiring no changes to the existing app interfaces.


title: Mass Transfer and Disc Formation in AGB Binary Systems

abstract: We investigate mass transfer and the formation of disc in binary systems using a combination of numerical simulations and theory. We consider six models distinguished by binary separation, secondary mass and outflow mechanisms. Each system consists of an asymptotic-giant-branch (AGB) star and an accreting secondary. The AGB star loses its mass via a wind. In one of our six models, the AGB star incurs a short period of outburst. In all cases, the secondary accretes part of the ejected mass and also influences the mass-loss rate of the AGB star. The ejected mass may remain gravitationally bound to the binary system and form a circumbinary disk, or contribute to an accretion disk around the secondary. In other cases, the ejecta will escape the binary system. The accretion rate on to the secondary changes non-linearly with binary separation. In our closest binary simulations, our models exemplify the wind Roche lobe overflow while in our wide binary cases, the mass transfer exhibits Bondi-Hoyle accretion. The morphologies of the outflows in the binary systems are varied. The variety may provide clues to how the late AGB phase influences planetary nebulae shaping. We employ the adaptive-mesh-refinement code ASTROBEAR for our simulations and include ray-tracing, radiation transfer, cooling and dust formation. To attain the highest computational efficiency and the most stable results, all simulations are run in the corotating frame.


title: Water Radiolysis by Low-Energy Carbon Projectiles from First-Principles Molecular Dynamics

abstract: Water radiolysis by low-energy carbon projectiles is studied by first-principles molecular dynamics. Carbon projectiles of kinetic energies between 175 eV and 2.8 keV are shot across liquid water. Apart from translational, rotational and vibrational excitation, they produce water dissociation. The most abundant products are H and OH fragments. We find that the maximum spatial production of radiolysis products, not only occurs at low velocities, but also well below the maximum of energy deposition, reaching one H every 5 Ang at the lowest speed studied (1 Bohr/fs), dissociative collisions being more significant at low velocity while the amount of energy required to dissociate water is constant and much smaller than the projectile's energy. A substantial fraction of the energy transferred to fragments, especially for high velocity projectiles, is in the form of kinetic energy, such fragments becoming secondary projectiles themselves. High velocity projectiles give rise to well-defined binary collisions, which should be amenable to binary approximations. This is not the case for lower velocities, where multiple collision events are observed. H secondary projectiles tend to move as radicals at high velocity, as cations when slower. We observe the generation of new species such as hydrogen peroxide and formic acid. The former occurs when an O radical created in the collision process attacks a water molecule at the O site. The latter when the C projectile is completely stopped and reacts with two water molecules.


title: Survey of Automated Vulnerability Detection and Exploit Generation Techniques in Cyber Reasoning Systems

abstract: Software is everywhere, from mission critical systems such as industrial power stations, pacemakers and even household appliances. This growing dependence on technology and the increasing complexity software has serious security implications as it means we are potentially surrounded by software that contain exploitable vulnerabilities. These challenges have made binary analysis an important area of research in computer science and has emphasized the need for building automated analysis systems that can operate at scale, speed and efficacy; all while performing with the skill of a human expert. Though great progress has been made in this area of research, there remains limitations and open challenges to be addressed. Recognizing this need, DARPA sponsored the Cyber Grand Challenge (CGC), a competition to showcase the current state of the art in systems that perform; automated vulnerability detection, exploit generation and software patching. This paper is a survey of the vulnerability detection and exploit generation techniques, underlying technologies and related works of two of the winning systems Mayhem and Mechanical Phish.


title: 1-Fan-Bundle-Planar Drawings of Graphs

abstract: Edge bundling is an important concept heavily used for graph visualization purposes. To enable the comparison with other established near-planarity models in graph drawing, we formulate a new edge-bundling model which is inspired by the recently introduced fan-planar graphs. In particular, we restrict the bundling to the end segments of the edges. Similarly to 1-planarity, we call our model 1-fan-bundle-planarity, as we allow at most one crossing per bundle. For the two variants where we allow either one or, more naturally, both end segments of each edge to be part of bundles, we present edge density results and consider various recognition questions, not only for general graphs, but also for the outer and 2-layer variants. We conclude with a series of challenging questions.


title: Search for standard model production of four top quarks in proton-proton collisions at sqrt(s) = 13 TeV

abstract: A search for events containing four top quarks (ttbar-ttbar) is reported from proton-proton collisions recorded by the CMS experiment at sqrt(s) = 13 TeV and corresponding to an integrated luminosity of 2.6 inverse femtobarns. The analysis considers the single-lepton (e or mu)+jets and the opposite-sign dilepton (mu+mu-, mu+/- e-/+, or e+e-)+jets channels. It uses boosted decision trees to combine information on the global event and jet properties to distinguish between ttbar-ttbar and ttbar production. The number of events observed after all selection requirements is consistent with expectations from background and standard model signal predictions, and an upper limit is set on the cross section for ttbar-ttbar production in the standard model of 94 fb at 95% confidence level (10.2 times the prediction), with an expected limit of 118 fb. This is combined with the results from the published CMS search in the same-sign dilepton channel, resulting in an improved limit of 69 fb at 95% confidence level (7.4 times the prediction), with an expected limit of 71 fb. These are the strongest constraints on the rate of ttbar-ttbar production to date.


title: Deep spectroscopy of nearby galaxy clusters: II. The Hercules cluster

abstract: We carried out the deep spectroscopic observations of the nearby cluster A2151 with AF2/WYFFOS@WHT. The caustic technique enables us to identify 360 members brighter than $M_r = -16$ and within 1.3$R_{200}$. We separated the members into subsamples according to photometrical and dynamical properties such as colour, local environment and infall time. The completeness of the catalogue and our large sample allow us to analyse the velocity dispersion and the luminosity functions of the identified populations. We found evidence of a cluster still in its collapsing phase. The LF of the red population of A2151 shows a deficit of dwarf red galaxies. Moreover, the normalized LFs of the red and blue populations of A2151 are comparable to the red and blue LFs of the field, even if the blue galaxies start dominating one magnitude fainter and the red LF is well represented by a single Schechter function rather than a double Schechter function. We discuss how the evolution of cluster galaxies depends on their mass: bright and intermediate galaxies are mainly affected by dynamical friction and internal/mass quenching, while the evolution of dwarfs is driven by environmental processes which need time and a hostile cluster environment to remove the gas reservoirs and halt the star formation.


title: Bayesian Boolean Matrix Factorisation

abstract: Boolean matrix factorisation aims to decompose a binary data matrix into an approximate Boolean product of two low rank, binary matrices: one containing meaningful patterns, the other quantifying how the observations can be expressed as a combination of these patterns. We introduce the OrMachine, a probabilistic generative model for Boolean matrix factorisation and derive a Metropolised Gibbs sampler that facilitates efficient parallel posterior inference. On real world and simulated data, our method outperforms all currently existing approaches for Boolean matrix factorisation and completion. This is the first method to provide full posterior inference for Boolean Matrix factorisation which is relevant in applications, e.g. for controlling false positive rates in collaborative filtering and, crucially, improves the interpretability of the inferred patterns. The proposed algorithm scales to large datasets as we demonstrate by analysing single cell gene expression data in 1.3 million mouse brain cells across 11 thousand genes on commodity hardware.

categories: stat.ML cs.LG cs.NA q-bio.GN q-bio.QM stat.ME


title: A Rollback in the History of Communication-Induced Checkpointing

abstract: The literature on communication-induced checkpointing presents a family of protocols that use logical clocks to control whether forced checkpoints must be taken. Efficiency of these protocols is measured by how many forced checkpoints are needed to ensure no checkpoint will be useless to the application; the fewer forced checkpoints the better. For many years, HMNR, also called Fully Informed (FI), was the most complex and efficient protocol of this family. The Lazy-FI protocol applies a lazy strategy that defers the increase of logical clocks, resulting in a protocol with better efficiency for distributed systems where processes can take basic checkpoints at different rates. Recently, the Fully Informed aNd Efficient (FINE) protocol was proposed using the same control structures as FI, but with a stronger and, presumably better, checkpoint-inducing condition. FINE and its lazy version, called Lazy-FINE, would now be the most efficient checkpointing protocols based on logical clocks. This paper reviews this family of protocols, proves a theorem on a condition that must be enforced by all stronger versions of FI, and proves that both FINE and Lazy-FINE do not guarantee the absence of useless checkpoints. As a consequence, FI and Lazy-FI can be rolled back to the position of most efficient protocols of this family of index-based checkpointing protocols.


title: Banach Algebra of Complex Bounded Radon Measures on Homogeneous Space

abstract: Let $ H $ be a compact subgroup of a locally compact group $G$. In this paper we define a convolution on $ M(G/H) $, the space of all complex bounded Radon measures on the homogeneous space G/H. Then we prove that the measure space $ M(G/H, *) $ is a non-unital Banach algebra that possesses an approximate identity. Finally, it is shown that the Banach algebra $ M(G/H, *) $ is not involutive and also $ L^1(G/H, *) $ is a two-sided ideal of it.


title: Critical points of master functions and mKdV hierarchy of type $A^{(2)}_{2n}$

abstract: We consider the population of critical points generated from the critical point of the master function with no variables, which is associated with the trivial representation of the twisted affine Lie algebra $A^{(2)}_{2n}$. The population is naturally partitioned into an infinite collection of complex cells $\mathbb{C}^m$, where $m$ are some positive integers. For each cell we define an injective rational map $\mathbb{C}^m \to M(A^{(2)}_{2n})$ of the cell to the space $M(A^{(2)}_{2n})$ of Miura opers of type $A^{(2)}_{2n}$. We show that the image of the map is invariant with respect to all mKdV flows on $M(A^{(2)}_{2n})$ and the image is point-wise fixed by all mKdV flows $\frac\partial{\partial t_r}$ with index $r$ greater than $4m$.


title: Limit multiplicities for ${\rm SL}_2(\mathcal{O}_F)$ in ${\rm SL}_2(\mathbb{R}^{r_1}\oplus\mathbb{C}^{r_2})$

abstract: We prove that the family of lattices ${\rm SL}_2(\mathcal{O}_F)$, $F$ running over number fields with fixed archimedean signature $(r_1, r_2)$, in ${\rm SL}_2(\mathbb{R}^{r_1}\oplus\mathbb{C}^{r_2})$ has the limit multiplicity property.


title: Continuous deformations of harmonic maps and their unitons

abstract: We consider harmonic maps on simply connected Riemann surfaces into the group $\mathrm{U}(n)$ of unitary matrices of order $n$. It is known that a harmonic map with an associated algebraic extended solution can be deformed into a new harmonic map that has an $S^1$-invariant associated extended solution. We study this deformation in detail and show that the corresponding unitons are smooth functions of the deformation parameter and real analytic along any line through the origin.


title: Exponential B-spline Collocation Solutions to the Gardner Equation

abstract: The exponential B-spline basis function set is used to develop a collocation method for some initial boundary value problems (IBVPs) to the Gardner equation. The Gardner equation has two nonlinear terms, namely quadratic and cubic ones. The order reduction of the equation is resulted in a coupled system of PDEs that enables the exponential B-splines to be implemented. The system is integrated in time by Crank-Nicolson implicit method. The validity of the method is investigated by calculating the discrete maximum error norm and observing the absolute relative changes of the conservation laws at the end of the simulations.


title: Asymmetric metallicity patterns in the stellar velocity space with RAVE

abstract: We explore the correlations between velocity and metallicity and the possible distinct chemical signatures of the velocity over-densities of the local Galactic neighbourhood. We use the large spectroscopic survey RAVE and the Geneva Copenhagen Survey. We compare the metallicity distribution of regions in the velocity plane ($v_R,v_\phi$) with that of their symmetric counterparts ($-v_R,v_\phi$). We expect similar metallicity distributions if there are no tracers of a sub-population (e.g., a dispersed cluster, accreted stars), if the disk of the Galaxy is axisymmetric, and if the orbital effects of the spiral arms and the bar are weak. We find that the metallicity-velocity space of the solar neighbourhood is highly patterned. A large fraction of the velocity plane shows differences in the metallicity distribution when comparing symmetric $v_R$ regions. The typical differences in the median metallicity are of $0.05$ dex with a statistical significance of at least $95\%$, and with values up to $0.6$ dex. For low azimuthal velocity $v_\phi$, stars moving outwards in the Galaxy have on average higher metallicity than those moving inwards. These include stars in the Hercules and Hyades moving groups and other velocity branch-like structures. For higher $v_\phi$, the stars moving inwards have higher metallicity than those moving outwards. The most likely interpretation of the metallicity asymmetry is that it is due to the orbital effects of the bar and the radial metallicity gradient of the disk. We present a simulation that supports this idea. We have also discovered a positive gradient in $v_\phi$ with respect to metallicity at high metallicities, apart from the two known positive and negative gradients for the thick and thin disks, respectively.


title: Electron-ion coupling in semiconductors beyond Fermi's golden rule

abstract: In the present work, a theoretical study of electron-phonon (electron-ion) coupling rates in semiconductors driven out of equilibrium is performed. Transient change of optical coefficients reflects the band gap shrinkage in covalently bonded materials, and thus, the heating of atomic lattice. Utilizing this dependence, we test various models of electron-ion coupling. The simulation technique is based on tight-binding molecular dynamics. Our simulations with the dedicated hybrid approach (XTANT) indicate that the widely used Fermi's golden rule can break down describing material excitation on femtosecond time scales. In contrast, dynamical coupling proposed in this work yields a reasonably good agreement of simulation results with available experimental data.


title: Structured signal recovery from quadratic measurements: Breaking sample complexity barriers via nonconvex optimization

abstract: This paper concerns the problem of recovering an unknown but structured signal $x \in R^n$ from $m$ quadratic measurements of the form $y_r=|<a_r,x>|^2$ for $r=1,2,...,m$. We focus on the under-determined setting where the number of measurements is significantly smaller than the dimension of the signal ($m<<n$). We formulate the recovery problem as a nonconvex optimization problem where prior structural information about the signal is enforced through constrains on the optimization variables. We prove that projected gradient descent, when initialized in a neighborhood of the desired signal, converges to the unknown signal at a linear rate. These results hold for any constraint set (convex or nonconvex) providing convergence guarantees to the global optimum even when the objective function and constraint set is nonconvex. Furthermore, these results hold with a number of measurements that is only a constant factor away from the minimal number of measurements required to uniquely identify the unknown signal. Our results provide the first provably tractable algorithm for this data-poor regime, breaking local sample complexity barriers that have emerged in recent literature. In a companion paper we demonstrate favorable properties for the optimization problem that may enable similar results to continue to hold more globally (over the entire ambient space). Collectively these two papers utilize and develop powerful tools for uniform convergence of empirical processes that may have broader implications for rigorous understanding of constrained nonconvex optimization heuristics. The mathematical results in this paper also pave the way for a new generation of data-driven phase-less imaging systems that can utilize prior information to significantly reduce acquisition time and enhance image reconstruction, enabling nano-scale imaging at unprecedented speeds and resolutions.

categories: cs.LG cs.IT math.FA math.IT math.OC stat.ML


title: MOLIERE: Automatic Biomedical Hypothesis Generation System

abstract: Hypothesis generation is becoming a crucial time-saving technique which allows biomedical researchers to quickly discover implicit connections between important concepts. Typically, these systems operate on domain-specific fractions of public medical data. MOLIERE, in contrast, utilizes information from over 24.5 million documents. At the heart of our approach lies a multi-modal and multi-relational network of biomedical objects extracted from several heterogeneous datasets from the National Center for Biotechnology Information (NCBI). These objects include but are not limited to scientific papers, keywords, genes, proteins, diseases, and diagnoses. We model hypotheses using Latent Dirichlet Allocation applied on abstracts found near shortest paths discovered within this network, and demonstrate the effectiveness of MOLIERE by performing hypothesis generation on historical data. Our network, implementation, and resulting data are all publicly available for the broad scientific community.


title: Coinfection in a stochastic model for bacteriophage systems

abstract: A system modeling bacteriophage treatments with coinfections in a noisy context is analyzed. We prove that in a small noise regime, the system converges in the long term to a bacteria free equilibrium. Moreover, we compare the treatment with coinfection with the treatment without coinfection, showing how the coinfection affects the dose of bacteriophages that is needed to eliminate the bacteria and the velocity of convergence to the free bacteria equilibrium.


title: The Twist of the Draped Interstellar Magnetic Field Ahead of the Heliopause: A Magnetic Reconnection Driven Rotational Discontinuity

abstract: Based on the difference between the orientation of the interstellar $B_{ISM}$ and the solar magnetic fields, there was an expectation that the magnetic field direction would rotate dramatically across the heliopause (HP). However, the Voyager 1 spacecraft measured very little rotation across the HP. Previously we showed that the $B_{ISM}$ twists as it approaches the HP and acquires a strong T component (East-West). Here we establish that reconnection in the eastern flank of the heliosphere is responsible for the twist. On the eastern flank the solar magnetic field has twisted into the positive N direction and reconnects with the Southward pointing component of the $B_{ISM}$. Reconnection drives a rotational discontinuity (RD) that twists the $B_{ISM}$ into the -T direction and propagates upstream in the interstellar medium towards the nose. The consequence is that the N component of $B_{ISM}$ is reduced in a finite width band upstream of the HP. Voyager 1 currently measures angles ($\delta=sin^{-1}(B_{N}/B)$) close to solar values. We present MHD simulations to support this scenario, suppressing reconnection in the nose region while allowing it in the flanks, consistent with recent ideas about reconnection suppression from diamagnetic drifts. The jump in plasma $\beta$ (the plasma to magnetic pressure) across the nose of HP is much greater than in the flanks because the heliosheath $\beta$ is greater there than in the flanks. Large-scale reconnection is therefore suppressed in the nose but not at the flanks. Simulation data suggest that $B_{ISM}$ will return to its pristine value $10-15~AU$ past the HP.


title: Square root of gerbe holonomy and invariants of time-reversal-symmetric topological insulators

abstract: The Feynman amplitudes with the two-dimensional Wess-Zumino action functional have a geometric interpretation as bundle gerbe holonomy. We present details of the construction of a distinguished square root of such holonomy and of a related 3d-index and briefly recall the application of those to the building of topological invariants for time-reversal-symmetric two- and three-dimensional crystals, both static and periodically forced.


title: Stochastic epidemic SEIRS models with a constant latency period

abstract: In this paper we consider the stability of a class of deterministic and stochastic SEIRS epidemic models with delay. Indeed, we assume that the transmission rate could be stochastic and the presence of a latency period of $r$ consecutive days, where $r$ is a fixed positive integer, in the "exposed" individuals class E. Studying the eigenvalues of the linearized system, we obtain conditions for the stability of the free disease equilibrium, in both the cases of the deterministic model with and without delay. In this latter case, we also get conditions for the stability of the coexistence equilibrium. In the stochastic case we are able to derive a concentration result for the random fluctuations and then, using the Lyapunov method, that under suitable assumptions the free disease equilibrium is still stable.


title: Lie symmetry and the Bethe ansatz solution of a new quasi-exactly solvable double-well potential

abstract: In this paper, we study the Schr\"odinger equation with a new quasi-exactly solvable double-well potential. Exact expressions for the energies, the corresponding wave functions and the allowed values of the potential parameters are obtained using two different methods, the Bethe ansatz method and the Lie algebraic approach. Some numerical results are reported and it is shown that the results are in good agreement with each other and with those obtained previously via a different method.


title: Evolution of Dust-Obscured Star Formation and Gas to z=2.2 from HiZELS

abstract: We investigate the far-infrared properties of galaxies selected via deep, narrow-band imaging of the H$\alpha$ emission line in four redshift slices from $z=0.40$--$2.23$ over $\sim 1$deg$^2$ as part of the High-redshift Emission Line Survey (HiZELS). We use a stacking approach in the Herschel PACS/SPIRE bands, along with $850\,\mu$m imaging from SCUBA-2 to study the evolution of the dust properties of H$\alpha$-emitters selected above an evolving characteristic luminosity threshold, $0.2L^\star_{{\rm H}\alpha}(z)$. We investigate the relationship between the dust temperatures and the far-infrared luminosities of our stacked samples, finding that H$\alpha$-selection identifies cold, low-$L_{\rm IR}$ galaxies ($T_{\rm dust}\sim 14$k; $\log[L_{\rm IR}/{\rm L}_\odot]\sim 9.9$) at $z=0.40$, and more luminous, warmer systems ($T_{\rm dust}\sim 34$k; $\log[L_{\rm IR}/{\rm L}_\odot]\sim 11.5$) at $z=2.23$. Using a modified greybody model, we estimate "characteristic sizes" for the dust-emitting regions of HiZELS galaxies of $\sim 0.5$kpc, nearly an order of magnitude smaller than their stellar continuum sizes, which may provide indirect evidence of clumpy ISM structure. Lastly, we measure the dust masses from our far-IR SEDs along with metallicity-dependent gas-to-dust ratios ($\delta_{\rm GDR}$) to measure typical molecular gas masses of $\sim 10^{10}$M$_\odot$ for these bright H$\alpha$-emitters. The gas depletion timescales are shorter than the Hubble time at each redshift, suggesting probable replenishment of their gas reservoirs from the intergalactic medium. Based on the number density of H$\alpha$-selected galaxies, we find that typical star-forming galaxies brighter than $0.2L^{\star}_{{\rm H}\alpha}(z)$ host a significant fraction ($35\pm10$%) of the total gas content of the Universe, consistent with the predictions of the latest cosmological simulations.


title: New reconstruction of the sunspot group number since 1739 using the direct calibration and "backbone" methods

abstract: Group sunspot number (GSN) series constitute the longest instrumental astronomical database providing information on solar activity. It is a compilation of observations by many individual observers, and their inter-calibration has usually been performed using linear rescaling. There are multiple published series that show different long-term trends for solar activity. We aim at producing a GSN series, with a non-linear non-parametric calibration. The only underlying assumptions are that the differences between the various series are due to different acuity thresholds of the observers, and that the threshold of each observer remains constant throughout the observing period. We use a daisy chain process with backbone (BB) observers and calibrate all overlapping observers to them. The calibration of each individual observer is performed with a probability distribution function (PDF) matrix, constructed considering all daily values for the overlapping period with the BB. The calibration of the BBs is done in a similar manner. Propagation of errors is modelled with Monte Carlo simulations. The final series extends back to 1739 and includes data from 314 observers. It suggests moderate activity during the 18th and 19th century, which is significantly lower than the high level of solar activity predicted by other recent reconstructions applying linear regressions. The new series provides a robust reconstruction, based on modern and non-parametric methods, of sunspot group numbers since 1739, and it confirms the existence of the Modern grand maximum of solar activity in the second half of the 20th century.