title: Operator bases, $S$-matrices, and their partition functions

abstract: Relativistic quantum systems that admit scattering experiments are quantitatively described by effective field theories, where $S$-matrix kinematics and symmetry considerations are encoded in the operator spectrum of the EFT. In this paper we use the $S$-matrix to derive the structure of the EFT operator basis, providing complementary descriptions in (i) position space utilizing the conformal algebra and cohomology and (ii) momentum space via an algebraic formulation in terms of a ring of momenta with kinematics implemented as an ideal. These frameworks systematically handle redundancies associated with equations of motion (on-shell) and integration by parts (momentum conservation). We introduce a partition function, termed the Hilbert series, to enumerate the operator basis--correspondingly, the $S$-matrix--and derive a matrix integral expression to compute the Hilbert series. The expression is general, easily applied in any spacetime dimension, with arbitrary field content and (linearly realized) symmetries. In addition to counting, we discuss construction of the basis. Simple algorithms follow from the algebraic formulation in momentum space. We explicitly compute the basis for operators involving up to $n=5$ scalar fields. This construction universally applies to fields with spin, since the operator basis for scalars encodes the momentum dependence of $n$-point amplitudes. We discuss in detail the operator basis for non-linearly realized symmetries. In the presence of massless particles, there is freedom to impose additional structure on the $S$-matrix in the form of soft limits. The most na\"ive implementation for massless scalars leads to the operator basis for pions, which we confirm using the standard CCWZ formulation for non-linear realizations.


title: The Global Gauge Group Structure of F-theory Compactification with U(1)s

abstract: We show that F-theory compactifications with abelian gauge factors generally exhibit a non-trivial global gauge group structure. The geometric origin of this structure lies with the Shioda map of the Mordell--Weil generators. This results in constraints on the U(1) charges of non-abelian matter consistent with observations made throughout the literature. In particular, we find that F-theory models featuring the Standard Model algebra actually realise the precise gauge group [SU(3)xSU(2)xU(1)]/Z6. Furthermore, we explore the relationship between the gauge group structure and geometric (un-)higgsing. In an explicit class of models, we show that, depending on the global group structure, an SU(2)xU(1) gauge theory can either unhiggs into an SU(2)xSU(2) or an SU(3)xSU(2) theory. We also study implications of the charge constraints as a criterion for the F-theory 'swampland'.


title: Pure Natural Inflation

abstract: We point out that a simple inflationary model in which the axionic inflaton couples to a pure Yang-Mills theory may give the scalar spectral index (n_s) and tensor-to-scalar ratio (r) in complete agreement with the current observational data.


title: Perturbative and global anomalies in bosonic analogues of integer quantum Hall and topological insulator phases

abstract: We study perturbative and global anomalies at the boundaries of bosonic analogues of integer quantum Hall (BIQH) and topological insulator (BTI) phases using a description of the boundaries of these phases in terms of a nonlinear sigma model (NLSM) with Wess-Zumino term. One of the main results of the paper is that these anomalies are robust against arbitrary smooth deformations of the target space of the NLSM which describes the phase, provided that the deformations also respect the symmetry of the phase. In the first part of the paper we discuss the perturbative $U(1)$ anomaly at the boundary of BIQH states in all odd (spacetime) dimensions. In the second part we study global anomalies at the boundary of BTI states in even dimensions. In a previous work [Phys. Rev. B 95, 035149 (2017)] we argued that the boundary of the BTI phase exhibits a global anomaly which is an analogue of the parity anomaly of Dirac fermions in three dimensions. Here we elevate this argument to a proof for the boundary of the two-dimensional BTI state by explicitly computing the partition function of the gauged NLSM describing the boundary. We then use the powerful equivariant localization technique to show that this global anomaly is robust against all smooth deformations of the target space of the NLSM which preserve the $U(1)\rtimes\mathbb{Z}_2$ symmetry of the BTI state. We also comment on the difficulties of generalizing this latter proof to higher dimensions. Finally, we discuss the expected low energy behavior of the boundary theories studied in this paper when the coupling constants are allowed to flow under the renormalization group.


title: From the trees to the forest: a review of radiative neutrino mass models

abstract: A plausible explanation for the lightness of neutrino masses is that neutrinos are massless at tree level, with their mass (typically Majorana) being generated radiatively at one or more loops. The new couplings, together with the suppression coming from the loop factors, imply that the new degrees of freedom cannot be too heavy (they are typically at the TeV scale). Therefore, in these models there are no large mass hierarchies and they can be tested using different searches, making their detailed phenomenological study very appealing. In particular, the new particles can be searched for at colliders and generically induce signals in lepton-flavor and lepton-number violating processes (in the case of Majorana neutrinos), which are not independent from reproducing correctly the neutrino masses and mixings. The main focus of the review is on Majorana neutrinos. We order the allowed theory space from three different perspectives: (i) using an effective operator approach to lepton number violation, (ii) by the number of loops at which the Weinberg operator is generated, (iii) within a given loop order, by the possible irreducible topologies. We also discuss in more detail some popular radiative models which involve qualitatively different features, revisiting their most important phenomenological implications. Finally, we list some promising avenues to pursue.


title: High $\beta$ Effects on Cosmic Ray Streaming in Galaxy Clusters

abstract: Diffuse, extended radio emission in galaxy clusters, commonly referred to as radio halos, indicate the presence of high energy cosmic ray (CR) electrons and cluster-wide magnetic fields. We can predict from theory the expected surface brightness of a radio halo, given magnetic field and CR density profiles. Previous studies have shown that the nature of CR transport can radically effect the expected radio halo emission from clusters (Wiener et al. 2013). Reasonable levels of magnetohydrodynamic (MHD) wave damping can lead to significant CR streaming speeds. But a careful treatment of MHD waves in a high $\beta$ plasma, as expected in cluster environments, reveals damping rates may be enhanced by a factor of $\beta^{1/2}$. This leads to faster CR streaming and lower surface brightnesses than without this effect. In this work we re-examine the simplified, 1D Coma cluster simulations (with radial magnetic fields) of Wiener et al. (2013) and discuss observable consequences of this high $\beta$ damping. Future work is required to study this effect in more realistic simulations.


title: On the Light Massive Flavor Dependence of the Large Order Asymptotic Behavior and the Ambiguity of the Pole Mass

abstract: We provide a systematic renormalization group formalism for the mass effects in the relation of the pole mass $m_Q^{\rm pole}$ and short-distance masses such as the $\overline{\rm MS}$ mass $\overline{m}_Q$ of a heavy quark $Q$, coming from virtual loop insertions of massive quarks lighter than $Q$. The formalism reflects the constraints from heavy quark symmetry and entails a combined matching and evolution procedure that allows to disentangle and successively integrate out the corrections coming from the lighter massive quarks and the momentum regions between them and to precisely control the large order asymptotic behavior. With the formalism we systematically sum logarithms of ratios of the lighter quark masses and $m_Q$, relate the QCD corrections for different external heavy quarks to each other, predict the ${\cal O}(\alpha_s^4)$ virtual quark mass corrections in the pole-$\overline{\rm MS}$ mass relation, calculate the pole mass differences for the top, bottom and charm quarks with a precision of around $20$ MeV and analyze the decoupling of the lighter massive quark flavors at large orders. The summation of logarithms is most relevant for the top quark pole mass $m_t^{\rm pole}$, where the hierarchy to the bottom and charm quarks is large. We determine the ambiguity of the pole mass for top, bottom and charm quarks in different scenarios with massive or massless bottom and charm quarks in a way consistent with heavy quark symmetry, and we find that it is $250$ MeV. The ambiguity is larger than current projections for the precision of top quark mass measurements in the high-luminosity phase of the LHC.


title: Combinatorics and Topology of Kawai-Lewellen-Tye Relations

abstract: We revisit the relations between open and closed string scattering amplitudes discovered by Kawai, Lewellen, and Tye (KLT). We show that they emerge from the underlying algebro-topological identities known as the twisted period relations. In order to do so, we formulate tree-level string theory amplitudes in the language of twisted de Rham theory. There, open string amplitudes are understood as pairings between twisted cycles and cocycles. Similarly, closed string amplitudes are given as a pairing between two twisted cocycles. Finally, objects relating the two types of string amplitudes are the $\alpha'$-corrected bi-adjoint scalar amplitudes recently defined by the author [arXiv:1610.04230]. We show that they naturally arise as intersection numbers of twisted cycles. In this work we focus on the combinatorial and topological description of twisted cycles relevant for string theory amplitudes. In this setting, each twisted cycle is a polytope, known in combinatorics as the associahedron, together with an additional structure encoding monodromy properties of string integrals. In fact, this additional structure is given by higher-dimensional generalizations of the Pochhammer contour. An open string amplitude is then computed as an integral of a logarithmic form over an associahedron. We show that the inverse of the KLT kernel can be calculated from the knowledge of how pairs of associahedra intersect one another in the moduli space. In the field theory limit, contributions from these intersections localize to vertices of the associahedra, giving rise to the bi-adjoint scalar partial amplitudes.


title: Algebraic Cycles and Local Anomalies in F-Theory

abstract: We introduce a set of identities in the cohomology ring of elliptic fibrations which are equivalent to the cancellation of gauge and mixed gauge-gravitational anomalies in F-theory compactifications to four and six dimensions. The identities consist in (co)homological relations between complex codimension-two cycles. The same set of relations, once evaluated on elliptic Calabi-Yau three-folds and four-folds, is shown to universally govern the structure of anomalies and their Green-Schwarz cancellation in six- and four-dimensional F-theory vacua, respectively. We furthermore conjecture that these relations hold not only within the cohomology ring, but even at the level of the Chow ring, i.e. as relations among codimension-two cycles modulo rational equivalence. We verify this conjecture in non-trivial examples with Abelian and non-Abelian gauge groups factors. Apart from governing the structure of local anomalies, the identities in the Chow ring relate different types of gauge backgrounds on elliptically fibred Calabi-Yau four-folds.


title: Graph Theory Data for Topological Quantum Chemistry

abstract: Topological phases of noninteracting particles are distinguished by global properties of their band structure and eigenfunctions in momentum space. On the other hand, group theory as conventionally applied to solid-state physics focuses only on properties which are local (at high symmetry points, lines, and planes) in the Brillouin zone. To bridge this gap, we have previously [B. Bradlyn et al., Nature 547, 298--305 (2017)] mapped the problem of constructing global band structures out of local data to a graph construction problem. In this paper, we provide the explicit data and formulate the necessary algorithms to produce all topologically distinct graphs. Furthermore, we show how to apply these algorithms to certain "elementary" band structures highlighted in the aforementioned reference, and so identified and tabulated all orbital types and lattices that can give rise to topologically disconnected band structures. Finally, we show how to use the newly developed BANDREP program on the Bilbao Crystallographic Server to access the results of our computation.


title: Discerning Dark Energy Models with High-Redshift Standard Candles

abstract: Following the success of type Ia supernovae in constraining cosmologies at lower redshift $(z\lesssim2)$, effort has been spent determining if a similarly useful standardisable candle can be found at higher redshift. {In this work we determine the largest possible magnitude discrepancy between a constant dark energy $\Lambda$CDM cosmology and a cosmology in which the equation of state $w(z)$ of dark energy is a function of redshift for high redshift standard candles $(z\gtrsim2)$}. We discuss a number of popular parametrisations of $w(z)$ with two free parameters, $w_z$CDM cosmologies, including the Chevallier-Polarski-Linder and generalisation thereof, $n$CPL, as well as the Jassal-Bagla-Padmanabhan parametrisation. For each of these parametrisations we calculate and find extrema of $\Delta \mu$, the difference between the distance modulus of a $w_z$CDM cosmology and a fiducial $\Lambda$CDM cosmology as a function of redshift, given 68\% likelihood constraints on the parameters $P=(\Omega_{m,0}, w_0, w_a)$. The parameters are constrained using cosmic microwave background, baryon acoustic oscillations, and type Ia supernovae data using CosmoMC. We find that none of the tested cosmologies can deviate more than 0.05 mag from the fiducial $\Lambda$CDM cosmology at high redshift, implying that high redshift standard candles will not aid in discerning between a $w_z$CDM cosmology and the fiducial $\Lambda$CDM cosmology. Conversely, this implies that if high redshift standard candles are found to be in disagreement with $\Lambda$CDM at high redshift, then this is a problem not only for $\Lambda$CDM but for the entire family of $w_z$CDM cosmologies.


title: On the SFR-M$_*$ main sequence archetypal star-formation history and analytical models

abstract: We derive the SFH of MS galaxies showing how the SFH peak of a galaxy depends on its seed mass at e.g. z=5. Following the MS, galaxies undergo a drastic slow down of their stellar mass growth after reaching the peak of their SFH. According to abundance matching, these masses correspond to hot and massive DM halos which state could results in less efficient gas inflows on the galaxies and thus could be at the origin of the limited stellar mass growth. As a result, galaxies on the MS can enter the passive region of the UVJ diagram while still forming stars. The ability of the classical analytical SFHs to retrieve the SFR of galaxies from SED fitting is studied. Due to mathematical limitations, the exp-declining and delayed SFH struggle to model high SFR which starts to be problematic at z>2. The exp-rising and log-normal SFHs exhibit the opposite behavior with the ability to reach very high SFR, and thus model starburst galaxies, but not low values such as those expected at low redshift for massive galaxies. We show that these four analytical forms recover the SFR of MS galaxies with an error dependent on the model and the redshift. They are, however, sensitive enough to probe small variations of SFR within the MS but all the four fail to recover the SFR of rapidly quenched galaxies. However, these SFHs lead to an artificial gradient of age, parallel to the MS which is not exhibited by a simulated sample. This gradient is also produced on real data, using a sample of GOODS-South galaxies at 1.5<z<1.2. We propose a SFH composed of a delayed form to model the bulk of stellar population plus a flexibility in the recent SFH. This SFH provides very good estimates of the SFR of MS, starbursts, and rapidly quenched galaxies at all z. Furthermore, used on the GOODS-South sample, the age gradient disappears, showing its dependency on the SFH assumption made to perform the SED fitting.


title: Simulating the effect of high column density absorbers on the one-dimensional Lyman-alpha forest flux power spectrum

abstract: We measure the effect of high column density absorbing systems of neutral hydrogen (HI) on the one-dimensional (1D) Lyman-alpha forest flux power spectrum using cosmological hydrodynamical simulations from the Illustris project. High column density absorbers (which we define to be those with HI column densities $N(\mathrm{HI}) > 1.6 \times 10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$) cause broadened absorption lines with characteristic damping wings. These damping wings bias the 1D Lyman-alpha forest flux power spectrum by causing absorption in quasar spectra away from the location of the absorber itself. We investigate the effect of high column density absorbers on the Lyman-alpha forest using hydrodynamical simulations for the first time. We provide templates as a function of column density and redshift, allowing the flexibility to accurately model residual contamination, i.e., if an analysis selectively clips out the largest damping wings. This flexibility will improve cosmological parameter estimation, e.g., allowing more accurate measurement of the shape of the power spectrum, with implications for cosmological models containing massive neutrinos or a running of the spectral index. We provide fitting functions to reproduce these results so that they can be incorporated straightforwardly into a data analysis pipeline.


title: An extended cold gas absorber in a central cluster galaxy

abstract: We present the serendipitous discovery of an extended cold gas structure projected close to the brightest cluster galaxy (BCG) of the z=0.045 cluster Abell 3716, from archival integral field spectroscopy. The gas is revealed through narrow Na D line absorption, seen against the stellar light of the BCG, which can be traced for $\sim$25 kpc, with a width of 2-4 kpc. The gas is offset to higher velocity than the BCG (by $\sim$100 km/s), showing that it is infalling rather than outflowing; the intrinsic linewidth is $\sim$80 km/s (FWHM). Very weak H$\alpha$ line emission is detected from the structure, and a weak dust absorption feature is suggested from optical imaging, but no stellar counterpart has been identified. We discuss some possible interpretations for the absorber: as a projected low-surface-brightness galaxy, as a stream of gas that was stripped from an infalling cluster galaxy, or as a "retired" cool-core nebula filament.


title: CP-violation for Electroweak Baryogenesis from Dynamical CKM Matrix

abstract: We show that the CKM matrix can be the source of CP violation for electroweak baryogenesis if Yukawa couplings vary at the same time as the Higgs acquires its vacuum expectation value. This offers new avenues for explaining the baryon asymmetry of the universe. These ideas apply if the mechanism explaining the flavour structure of the Standard Model is connected to electroweak symmetry breaking. We compute the resulting baryon asymmetry for various low-scale flavour models and different configurations of the Yukawa coupling variation across the bubble wall, and show that it can naturally be of the right order.


title: SU(5) Unification without Proton Decay

abstract: We construct a four-dimensional SU(5) grand unified theory in which the proton is stable. The Standard Model leptons reside in the 5 and 10 irreps of SU(5), whereas the quarks live in the 40 and 50 irreps. The SU(5) gauge symmetry is broken by the vacuum expectation values of the scalar 24 and 75 irreps. All non-Standard Model fields are heavy. Stability of the proton requires three relations between the parameters of the model to hold. However, abandoning the requirement of absolute proton stability, the model fulfills current experimental constraints without fine-tuning.


title: Was the Universe Actually Radiation Dominated Prior to Nucleosynthesis?

abstract: Maybe not. String theory approaches to both beyond the Standard Model and Inflationary model building generically predict the existence of scalars (moduli) that are light compared to the scale of quantum gravity. These moduli become displaced from their low energy minima in the early universe and lead to a prolonged matter-dominated epoch prior to BBN. In this paper, we examine whether non-perturbative effects such as parametric resonance or tachyonic instabilities can shorten, or even eliminate, the moduli condensate and matter-dominated epoch. Such effects depend crucially on the strength of the couplings, and we find that unless the moduli become strongly coupled the matter-dominated epoch is unavoidable. In particular, we find that in string and M-theory compactifications where the lightest moduli are near the TeV-scale that a matter-dominated epoch will persist until the time of Big Bang Nucleosynthesis.


title: The Not-So-Well Tempered Neutralino

abstract: Light electroweakinos, the neutral and charged fermionic supersymmetric partners of the Standard Model $\mathrm{SU}(2)\times\mathrm{U}(1)$ gauge bosons and of the two $\mathrm{SU}(2)$ Higgs doublets, are an important target for searches for new physics with the Large Hadron Collider (LHC). However, if the lightest neutralino is the dark matter, constraints from direct dark matter detection experiments rule out large swaths of the parameter space accessible to the LHC, including in large part the so-called "well-tempered" neutralinos. We focus on the Minimal Supersymmetric Standard Model (MSSM) and explore in detail which regions of parameter space are not excluded by null results from direct dark matter detection, assuming exclusive thermal production of neutralinos in the early universe, and illustrate the complementarity with current and future LHC searches for electroweak gauginos. We consider both bino-Higgsino and bino-wino "not-so-well-tempered" neutralinos, i.e. we include models where the lightest neutralino constitutes only part of the cosmological dark matter, with the consequent suppression of the constraints from direct and indirect dark matter searches.


title: A deceleration search for magnetar pulsations in the X-ray plateaus of Short GRBs

abstract: A newly formed magnetar has been proposed as the central engine of short GRBs to explain on-going energy injection giving observed plateau phases in the X-ray lightcurves. These rapidly spinning magnetars may be capable of emitting pulsed emission comparable to known pulsars and magnetars. In this paper we show that, if present, a periodic signal would be detectable during the plateau phases observed using the Swift X-Ray Telescope recording data in Window Timing mode. We conduct a targeted deceleration search for a periodic signal from a newly formed magnetar in 2 Swift short GRBs and rule out any periodic signals in the frequency band 10-285 Hz to $\approx$15-30% rms. These results demonstrate that we would able to detect pulsations from the magnetar central engine of short GRBs if they contribute to 15-30% of the total emission. We consider these constraints in the context of the potential emission mechanisms. The non-detection is consistent with the emission being reprocessed in the surrounding environment or with the rotation axis being highly aligned with the observing angle. As the emission may be reprocessed, the expected periodic emission may only constitute a few percent of the total emission and be undetectable in our observations. Applying this strategy to future observations of the plateau phases with more sensitive X-ray telescopes may lead to the detection of the periodic signal.


title: The Galaxy in circular polarization: all-sky radio prediction, detection strategy, and the charge of the leptonic cosmic rays

abstract: The diffuse Galactic synchrotron emission should exhibit a low level of diffuse circular polarization (CP) due to the circular motions of the emitting relativistic electrons. This probes the Galactic magnetic field in a similar way as the product of total Galactic synchrotron intensity times Faraday depth. We use this to construct an all sky prediction of the so far unexplored Galactic CP from existing measurements. This map can be used to search for this CP signal in low frequency radio data even prior to imaging. If detected as predicted, it would confirm the expectation that relativistic electrons, and not positrons, are responsible for the Galactic radio emission. Furthermore, the strength of real to predicted circular polarization would provide statistical information on magnetic structures along the line-of-sights.


title: Cloud Scale ISM Structure and Star Formation in M51

abstract: We compare the structure of molecular gas at $40$ pc resolution to the ability of gas to form stars across the disk of the spiral galaxy M51. We break the PAWS survey into $370$ pc and $1.1$ kpc resolution elements, and within each we estimate the molecular gas depletion time ($\tau_{\rm Dep}^{\rm mol}$), the star formation efficiency per free fall time ($\epsilon_{\rm ff}$), and the mass-weighted cloud-scale (40 pc) properties of the molecular gas: surface density, $\Sigma$, line width, $\sigma$, and $b\equiv\Sigma/\sigma^2\propto\alpha_{\rm vir}^{-1}$, a parameter that traces the boundedness of the gas. We show that the cloud-scale surface density appears to be a reasonable proxy for mean volume density. Applying this, we find a typical star formation efficiency per free-fall time, $\epsilon_{ff} \left( \left< \Sigma_{40pc} \right> \right) \sim 0.3{-}0.36\%$, lower than adopted in many models and found for local clouds. More, the efficiency per free fall time anti-correlates with both $\Sigma$ and $\sigma$, in some tension with turbulent star formation models. The best predictor of the rate of star formation per unit gas mass in our analysis is $b \equiv \Sigma / \sigma^2$, tracing the strength of self gravity, with $\tau_{\rm Dep}^{\rm mol} \propto b^{-0.9}$. The sense of the correlation is that gas with stronger self-gravity (higher $b$) forms stars at a higher rate (low $\tau_{\rm Dep}^{\rm mol}$). The different regions of the galaxy mostly overlap in $\tau_{\rm Dep}^{\rm mol}$ as a function of $b$, so that low $b$ explains the surprisingly high $\tau_{\rm Dep}^{\rm mol}$ found towards the inner spiral arms found by by Meidt et al. (2013).


title: Exterior Distance Function

abstract: We introduce and study exterior distance function (EDF) and correspondent exterior point method (EPM) for convex optimization. The EDF is a classical Lagrangian for an equivalent problem obtained from the initial one by monotone transformation of both the objective function and the constraints. The constraints transformation is scaled by a positive scaling parameter. Thus, the EDF is a particular realization of the Nonlinear Rescaling (NR) principle. Along with the "center", the EDF has two extra tools: the barrier (scaling) parameter and the vector of Lagrange multipliers. We show that EPM generates primal - dual sequence, which converges to the primal - dual solution in value under minimum assumption on the input data. Moreover, the convergence is taking place under any fixed interior point as a "center" and any fixed positive scaling parameter, just due to the Lagrange multipliers update. If the second order sufficient optimality condition is satisfied, then the EPM converges with Q-linear rate under any fixed interior point as a "center" and any fixed, but large enough positive scaling parameter.


title: Are fossil groups early-forming galaxy systems?

abstract: Using the Illustris cosmological simulation, we investigate the origin of fossil groups in the $M_{200}=10^{13}-10^{13.5}M_{\odot}/h$ mass regime. We examine the formation of the two primary features of fossil groups: the large magnitude gap between their two brightest galaxies, and their exceptionally luminous brightest group galaxy (BGG). For fossils and non-fossils identified at $z=0$, we find no difference in their halo mass assembly at early times, departing from previous studies. However, we do find a significant difference in the recent accretion history of fossil and non-fossil halos; in particular, fossil groups show a lack of recent accretion and have in majority assembled 80\% of their $M_{200}(z=0)$ mass before $z\sim 0.4$. For fossils, massive satellite galaxies accreted during this period have enough time to merge with the BGG by the present day, producing a more massive central galaxy; and, the lack of recent group accretion prevents replenishment of the bright satellite population, allowing for a large magnitude gap to develop within the past few Gyr. We thus find that the origin of the magnitude gap and overmassive BGG of fossils in Illustris depends on the recent accretion history of the groups and merger history of the BGG after their collapse at $z\sim1$. This indicates that selecting galaxy groups by their magnitude gap does not guarantee obtaining neither early-forming galaxy systems nor undisturbed central galaxies.


title: Dark matter indirect signals with long-lived mediators

abstract: Dark matter particles could annihilate into light and metastable mediators subsequently decaying far away from where they are produced. In this scenario, the indirect signatures of dark matter are altered with respect to the conventional situation where standard model particles are directly injected where annihilation happens. We explore here this new phenomenology for which we devise the tools. We calculate the effective dark matter distribution resulting from the smearing by mediator propagation. We derive the fluxes of mediators and their decay particles. We study how the $J$-factor, which naturally appears in the calculation of the dark matter induced gamma ray signal, is modified in the presence of mediators. We also derive the anisotropy which the cosmic ray positron flux exhibits in this scenario. We finally comment upon a recent proposal based on long-lived mediators where the effective dark matter density at the Earth is increased such as to explain the cosmic ray positron anomaly. We conclude that this scenario is barely tenable as regards the very dense dark matter spike which it requires at the Galactic center. The associated positron anisotropy is very small and undetectable, except at high energies where it reaches a level of order $10^{-4}$ to $10^{-3}$.


title: Delay-coordinate maps and the spectra of Koopman operators

abstract: The Koopman operator induced by a dynamical system is inherently linear and provides an alternate method of studying many properties of the system, including attractor reconstruction and forecasting. Koopman eigenfunctions represent the non-mixing component of the dynamics. They factor the dynamics, which can be chaotic, into quasiperiodic rotations on tori. Here, we describe a method through which these eigenfunctions can be obtained from a kernel integral operator, which also annihilates the continuous spectrum. We show that incorporating a large number of delay coordinates in constructing the kernel of that operator results, in the limit of infinitely many delays, in the creation of a map into the discrete spectrum subspace of the Koopman operator. This enables efficient approximation of Koopman eigenfunctions from high-dimensional data in systems with pure point or mixed spectra.


title: Orthogonal Symmetric Chain Decompositions of Hypercubes

abstract: In 1979, Shearer and Kleitman conjectured that there exist $\lfloor n/2 \rfloor+1$ orthogonal chain decompositions of the hypercube $Q_n$, and constructed two orthogonal chain decompositions. In this paper, we make the first non-trivial progress on this conjecture since by constructing three orthogonal chain decompositions of $Q_n$ for $n$ large enough. To do this, we introduce the notion of "almost orthogonal symmetric chain decompositions". We explicitly describe three such decompositions of $Q_5$ and $Q_7$, and describe conditions which allow us to decompose products of hypercubes into $k$ almost orthogonal symmetric chain decompositions given such decompositions of the original hypercubes.


title: Symmetric Chain Decompositions of Products of Posets with Long Chains

abstract: We ask if there exists a symmetric chain decomposition of the cuboid $Q_k \times n$ such that no chain is "taut", i.e. no chain has a subchain of the form $(a_1,\ldots, a_k,0)\prec \ldots\prec (a_1,\ldots,a_k,n-1)$. In this paper, we show this is true precisely when $k \ge 5$ and $n\ge 3$. This question arises naturally when considering products of symmetric chain decompositions which induce orthogonal chain decompositions --- the existence of the decompositions provided in this paper unexpectedly resolves the most difficult case of previous work by the second author on almost orthogonal symmetric chain decompositions \cite{orth}, making progress on a conjecture of Shearer and Kleitman. In general, we show that for a finite graded poset $P$, there exists a canonical bijection between symmetric chain decompositions of $P \times m$ and $P \times n$ for $m, n\ge \text{rk}(P) + 1$, that preserves the existence of taut chains. If $P$ has a unique maximal and minimal element, then we also produce a canonical $(\text{rk}(P) +1)$ to $1$ surjection from symmetric chain decompositions of $P \times (\text{rk}(P) + 1)$ to symmetric chain decompositions of $P \times \text{rk}(P)$ which sends decompositions with taut chains to decompositions with taut chains.


title: 3D AMR hydrosimulations of a compact source scenario for the Galactic Centre cloud G2

abstract: The nature of the gaseous and dusty cloud G2 in the Galactic Centre is still under debate. We present three-dimensional hydrodynamical adaptive mesh refinement (AMR) simulations of G2, modeled as an outflow from a "compact source" moving on the observed orbit. The construction of mock position-velocity (PV) diagrams enables a direct comparison with observations and allow us to conclude that the observational properties of the gaseous component of G2 could be matched by a massive ($\dot{M}_\mathrm{w}=5\times 10^{-7} \;M_{\odot} \mathrm{yr^{-1}}$) and slow ($50 \;\mathrm{km \;s^{-1}}$) outflow, as observed for T Tauri stars. In order for this to be true, only the material at larger ($>100 \;\mathrm{AU}$) distances from the source must be actually emitting, otherwise G2 would appear too compact compared to the observed PV diagrams. On the other hand, the presence of a central dusty source might be able to explain the compactness of G2's dust component. In the present scenario, 5-10 years after pericentre the compact source should decouple from the previously ejected material, due to the hydrodynamic interaction of the latter with the surrounding hot and dense atmosphere. In this case, a new outflow should form, ahead of the previous one, which would be the smoking gun evidence for an outflow scenario.


title: Top Down Electroweak Dipole Operators

abstract: We derive present constraints on, and prospective sensitivity to, the electric dipole moment (EDM) of the top quark ($d_t$) implied by searches for the EDMs of the electron and nucleons. Above the electroweak scale $v$, the $d_t$ arises from two gauge invariant operators generated at a scale $\Lambda \gg v$ that also mix with the light fermion EDMs under renormalization group evolution at two-loop order. Bounds on the EDMs of first generation fermion systems thus imply bounds on $|d_t|$. Working in the leading log-squared approximation, we find that the present upper bound on $|d_t|$ is roughly $10^{-19}$ $e$ cm for $\Lambda = 1$ TeV, except in regions of finely tuned cancellations that allow for $|d_t|$ to be up to fifty times larger. Future $d_e$ and $d_n$ probes may yield an order of magnitude increase in $d_t$ sensitivity, while inclusion of a prospective proton EDM search may lead to an additional increase in reach.


title: Topology of (small) Lagrangian cobordisms

abstract: We study the following quantitative phenomenon in symplectic topology: In many situations, if a Lagrangian cobordism is sufficiently small (in a sense specified below) then its topology is to a large extend determined by its boundary. This principle allows us to derive several homological uniqueness results for small Lagrangian cobordisms. In particular, under the smallness assumption, we prove homological uniqueness of the class of Lagrangian cobordisms which, by Biran-Cornea's Lagrangian cobordism theory, induces operations on a version of the derived Fukaya category. We also establish a link between our results and Vassilyev's theory of Lagrange characteristic classes. Most currently known constructions of Lagrangian cobordisms yield small Lagrangian cobordisms in many examples.


title: Semidefinite Programming and Nash Equilibria in Bimatrix Games

abstract: We explore the power of semidefinite programming (SDP) for finding additive $epsilon$-approximate Nash equilibria in bimatrix games. We introduce an SDP relaxation for a quadratic programming formulation of the Nash equilibrium (NE) problem and provide a number of valid inequalities to improve the quality of the relaxation. If a rank-1 solution to this SDP is found, then an exact NE can be recovered. We show that for a strictly competitive game, our SDP is guaranteed to return a rank-1 solution. We propose two algorithms based on iterative linearization of smooth nonconvex objective functions whose global minima by design coincide with rank-1 solutions. Empirically, we demonstrate that these algorithms often recover solutions of rank at most two and $epsilon$ close to zero. Furthermore, we prove that if a rank-2 solution to our SDP is found, then a 5/11-NE can be recovered for any game, or a 1/3-NE for a symmetric game. We then show how our SDP approach can address two (NP-hard) problems of economic interest: finding the maximum welfare achievable under any NE, and testing whether there exists a NE where a particular set of strategies is not played. Finally, we show the connection between our SDP and the first level of the Lasserre/sum of squares hierarchy.


title: Universal phase transition and band structures for spinless nodal-line and Weyl semimetals

abstract: We study a general phase transition between spinless topological nodal-line semimetal and Weyl semimetal phases. We classify topological nodal lines into two types based on their positions and shapes, and their phase transitions depends on their types. We show that the topological nodal-line semimetal becomes the Weyl semimetal by breaking time-reversal symmetry when the nodal lines enclose time-reversal invariant momenta (type-A nodal lines). We also discuss an effect of crystallographic symmetries determining the band structure of the topological nodal-line semimetals. Thanks to protection by the symmetries, the topological nodal-line semimetals can transit into the spinless Weyl semimetals or maintain the nodal lines in many crystals after inversion symmetry is broken.


title: A spectral interpretation of zeros of certain functions

abstract: We prove that all the zeros of certain meromorphic functions are on the critical line $\text{Re}(s)=1/2$, and are simple (except possibly when $s=1/2$). We prove this by relating the zeros to the discrete spectrum of an unbounded self-adjoint operator. Specifically, we show for $h(s)$ a meromorphic function with no zeros in $\text{Re}(s)>1/2$ and no poles in $\text{Re}(s)<1/2$, real-valued on $\R$, $\frac{h(1-s)}{h(s)}\ll |s|^{1-\epsilon}$ in $\text{Re}(s)>1/2$ and $\frac{h(1-s)}{h(s)}\notin L^2(1/2+i\R)$, the only zeros of $h(s)\pm h(1-s)$ are on the critical line. One instance of such a function $h$ is $h(s)=\xi(2s)$, the completed zeta-function. We use spectral theory suggested by results of Lax-Phillips and Colin de Verdi\`{e}re. This simplifies ideas of W. M\"{u}ller, J. Lagarias, M. Suzuki, H. Ki, O. Vel\'{a}squez Casta\~{n}\'{o}n, D. Hejhal, L. de Branges and P.R. Taylor.


title: The role of vector potential coupling in hot electron cooling power in bilayer graphene at low temperature

abstract: We have studied, in bilayer graphene (BLG), the hot electron cooling power $F_{VP} (T, n_s)$ due to acoustic phonons via vector potential (VP) coupling. It is calculated as a function of electron concentration $ns$ and temperature $T$ and compared with $F_{DP} (T, n_s)$ that due to deformation potential coupling. For the $n_s$ around $1\times10^{12} cm^{-2}$, $F_{VP} (T, n_s)$ is much smaller than $F_{DP} (T, n_s)$ . With increase of $n_s$, $F_{DP} (T, n_s)$ decreases faster than $F_{VP} (T, n_s)$ does. A cross over is predicted and dominant contribution of $F_{VP} (T, n_s)$ can be observed at large $n_s$. In the Bloch- Gr\"uneisen (BG) regime $F_{VP} (T, n_s)~ n_s^{-1/2}$ and $F_{DP} (T, n_s)~ n_s^{-3/2}$. Both $F_{VP} (T, n_s)$ and $F_{DP} (T, n_s)$ have the same $T$ dependence with $T^4$ power law in BG regime. Behaviour of $F_{DP} (T, n_s)~ n_s^{-3/2}$ and $T^4$ is in agreement with the experimental results at moderate $n_s$. Besides, in BG regime, we have predicted, for both the VP and DP coupling, a relation between $F(T, n_s)$ and the acoustic phonon limited mobility $\mu_p$, opening a new door to determine $\mu_p$ from the measurements of $F(T, n_s)$.


title: Topological Equivalences of E-infinity Differential Graded Algebras

abstract: Two DGAs are called topologically equivalent if the corresponding Eilenberg-Mac Lane ring spectra are weakly equivalent as ring spectra. Quasi-isomorphic DGAs are topologically equivalent but the converse is not necessarily true. As a counter-example, Dugger and Shipley showed that there are DGAs that are non-trivially topologically equivalent, i.e. topologically equivalent but not quasi-isomorphic. In this work, we define $E_\infty$ topological equivalences and utilize the obstruction theories developed by Goerss, Hopkins and Miller to construct first examples of non-trivially $E_\infty$ topologically equivalent $E_\infty$ DGAs. Also, we show using these obstruction theories that for co-connective $E_\infty$ DGAs, $E_\infty$ topological equivalences and quasi-isomorphisms agree. For $E_\infty$ $\mathbb{F}_p$-DGAs with trivial first homology, we show that an $E_\infty$ topological equivalence induces an isomorphism in homology that preserves the Dyer-Lashof operations and therefore induces an $H_\infty $ $\mathbb{F}_p$-equivalence.


title: HI in Virgo's "Red and Dead" Dwarf Ellipticals - A Tidal Tail and Central Star Formation

abstract: We investigate a sample of 3 dwarf elliptical galaxies in the Virgo Cluster which have significant reservoirs of HI. We present deep optical imaging (from CFHT and KPNO), HI spectra (Arecibo) and resolved HI imaging (VLA) of this sample. These observations confirm their HI content and optical morphologies, and indicate that the gas is unlikely to be recently accreted. The sample has more in common with dwarf transitionals, although dwarf transitionals are generally lower in stellar mass and gas fraction. VCC 190 has an HI tidal tail from a recent encounter with the massive spiral galaxy NGC 4224. In VCC 611, blue star-forming features are observed which were unseen by shallower SDSS imaging.


title: Unsupervised Method for Correlated Noise Removal for Multi-wavelength Exoplanet Transit Observations

abstract: Exoplanetary atmospheric observations require an exquisite precision in the measurement of the relative flux among wavelengths. In this paper, we aim to provide a new adaptive method to treat light curves before fitting transit parameters in order to minimize systematic effects that affect, for instance, ground-based observations of exo-atmospheres. We propose a neural-network-based method that uses a reference built from the data itself with parameters that are chosen in an unsupervised fashion. To improve the performance of proposed method, K-means clustering and Silhouette criteria are used for identifying similar wavelengths in each cluster. We also constrain under which circumstances our method improves the measurement of planetary-to-stellar radius ratio without producing significant systematic offset. We tested our method in high quality data from WASP-19b and low quality data from GJ-1214. We succeed in providing smaller error bars for the former when using JKTEBOP, but GJ-1214 light curve was beyond the capabilities of this method to improve as it was expected from our validation tests.


title: New insights into non-central beta distributions

abstract: The beta family owes its privileged status within unit interval distributions to several relevant features such as, for example, easyness of interpretation and versatility in modeling different types of data. However, its flexibility at the unit interval endpoints is poor enough to prevent from properly modeling the portions of data having values next to zero and one. Such a drawback can be overcome by resorting to the class of the non-central beta distributions. Indeed, the latter allows the density to take on arbitrary positive and finite limits which have a really simple form. That said, new insights into such class are provided in this paper. In particular, new representations and moments expressions are derived. Moreover, its potential with respect to alternative models is highlighted through applications to real data.


title: Heat transport study of the spin liquid candidate 1T-TaS2

abstract: We present the ultra-low-temperature thermal conductivity measurements on single crystals of the prototypical charge-density-wave material 1$T$-TaS$_2$, which was recently argued to be a candidate for quantum spin liquid. Our experiments show that the residual linear term of thermal conductivity at zero field is essentially zero, within the experimental accuracy. Furthermore, the thermal conductivity is found to be insensitive to the magnetic field up to 9 T. These results clearly demonstrate the absence of itinerant magnetic excitations with fermionic statistics in bulk 1$T$-TaS$_2$ and, thus, put a strong constraint on the theories of the ground state of this material.


title: Polarization of Neural Rings

abstract: The "neural code" is the way the brain characterizes, stores, and processes information. Unraveling the neural code is a key goal of mathematical neuroscience. Topology, coding theory, and, recently, commutative algebra are some the mathematical areas that are involved in analyzing these codes. Neural rings and ideals are algebraic objects that create a bridge between mathematical neuroscience and commutative algebra. A neural ideal is an ideal in a polynomial ring that encodes the combinatorial firing data of a neural code. Using some algebraic techniques one hopes to understand more about the structure of a neural code via neural rings and ideals. In this paper, we introduce an operation, called "polarization," that allows us to relate neural ideals with squarefree monomial ideals, which are very well studied and known for their nice behavior in commutative algebra.


title: Skill Learning by Autonomous Robotic Playing using Active Learning and Creativity

abstract: We treat the problem of autonomous acquisition of manipulation skills where problem-solving strategies are initially available only for a narrow range of situations. We propose to extend the range of solvable situations by autonomous playing with the object. By applying previously-trained skills and behaviours, the robot learns how to prepare situations for which a successful strategy is already known. The information gathered during autonomous play is additionally used to learn an environment model. This model is exploited for active learning and the creative generation of novel preparatory behaviours. We apply our approach on a wide range of different manipulation tasks, e.g. book grasping, grasping of objects of different sizes by selecting different grasping strategies, placement on shelves, and tower disassembly. We show that the creative behaviour generation mechanism enables the robot to solve previously-unsolvable tasks, e.g. tower disassembly. We use success statistics gained during real-world experiments to simulate the convergence behaviour of our system. Experiments show that active improves the learning speed by around 9 percent in the book grasping scenario.


title: Group Synchronization on Grids

abstract: Group synchronization requires to estimate unknown elements $({\theta}_v)_{v\in V}$ of a compact group ${\mathfrak G}$ associated to the vertices of a graph $G=(V,E)$, using noisy observations of the group differences associated to the edges. This model is relevant to a variety of applications ranging from structure from motion in computer vision to graph localization and positioning, to certain families of community detection problems. We focus on the case in which the graph $G$ is the $d$-dimensional grid. Since the unknowns ${\boldsymbol \theta}_v$ are only determined up to a global action of the group, we consider the following weak recovery question. Can we determine the group difference ${\theta}_u^{-1}{\theta}_v$ between far apart vertices $u, v$ better than by random guessing? We prove that weak recovery is possible (provided the noise is small enough) for $d\ge 3$ and, for certain finite groups, for $d\ge 2$. Viceversa, for some continuous groups, we prove that weak recovery is impossible for $d=2$. Finally, for strong enough noise, weak recovery is always impossible.


title: Return Oriented Programming - Exploit Implementation using functions

abstract: Return Oriented programming was surfaced first a decade ago, and was built to overcome the buffer exploit defense mechanisms like ASLR, DEP (or W^ X) by method of reusing the system code in the form of gadgets which are stitched together to make a Turing complete attack. And to perform Turing complete attack would require greater efforts which are quite complex, and there is very little research available for performing a Turing complete attack. So, in this project, we are systematizing the knowledge of the existing research that can be used to perform a Turing complete ROP attack.


title: A Streaming Multi-GPU Implementation of Image Simulation Algorithms for Scanning Transmission Electron Microscopy

abstract: Simulation of atomic resolution image formation in scanning transmission electron microscopy can require significant computation times using traditional methods. A recently developed method, termed plane-wave reciprocal-space interpolated scattering matrix (PRISM), demonstrates potential for significant acceleration of such simulations with negligible loss of accuracy. Here we present a software package called Prismatic for parallelized simulation of image formation in scanning transmission electron microscopy (STEM) using both the PRISM and multislice methods. By distributing the workload between multiple CUDA-enabled GPUs and multicore processors, accelerations as high as 1000x for PRISM and 30x for multislice are achieved relative to traditional multislice implementations using a single 4-GPU machine. We demonstrate a potentially important application of Prismatic, using it to compute images for atomic electron tomography at sufficient speeds to include in the reconstruction pipeline. Prismatic is freely available both as an open-source CUDA/C++ package with a graphical user interface and as a Python package, PyPrismatic.


title: Illuminating Pedestrians via Simultaneous Detection & Segmentation

abstract: Pedestrian detection is a critical problem in computer vision with significant impact on safety in urban autonomous driving. In this work, we explore how semantic segmentation can be used to boost pedestrian detection accuracy while having little to no impact on network efficiency. We propose a segmentation infusion network to enable joint supervision on semantic segmentation and pedestrian detection. When placed properly, the additional supervision helps guide features in shared layers to become more sophisticated and helpful for the downstream pedestrian detector. Using this approach, we find weakly annotated boxes to be sufficient for considerable performance gains. We provide an in-depth analysis to demonstrate how shared layers are shaped by the segmentation supervision. In doing so, we show that the resulting feature maps become more semantically meaningful and robust to shape and occlusion. Overall, our simultaneous detection and segmentation framework achieves a considerable gain over the state-of-the-art on the Caltech pedestrian dataset, competitive performance on KITTI, and executes 2x faster than competitive methods.


title: Satellite conjunction analysis and the false confidence theorem

abstract: Satellite conjunction analysis is the assessment of collision risk during a close encounter between a satellite and another object in orbit. A counterintuitive phenomenon has emerged in the conjunction analysis literature, namely, probability dilution, in which lower quality data paradoxically appear to reduce the risk of collision. We show that probability dilution is a symptom of a fundamental deficiency in probabilistic representations of statistical inference, in which there are propositions that will consistently be assigned a high degree of belief, regardless of whether or not they are true. We call this deficiency false confidence. In satellite conjunction analysis, it results in a severe and persistent underestimate of collision risk exposure. We introduce the Martin--Liu validity criterion as a benchmark by which to identify statistical methods that are free from false confidence. Such inferences will necessarily be non-probabilistic. In satellite conjunction analysis, we show that uncertainty ellipsoids satisfy the validity criterion. Performing collision avoidance maneuvers based on ellipsoid overlap will ensure that collision risk is capped at the user-specified level. Further, this investigation into satellite conjunction analysis provides a template for recognizing and resolving false confidence issues as they occur in other problems of statistical inference.


title: SchNet: A continuous-filter convolutional neural network for modeling quantum interactions

abstract: Deep learning has the potential to revolutionize quantum chemistry as it is ideally suited to learn representations for structured data and speed up the exploration of chemical space. While convolutional neural networks have proven to be the first choice for images, audio and video data, the atoms in molecules are not restricted to a grid. Instead, their precise locations contain essential physical information, that would get lost if discretized. Thus, we propose to use continuous-filter convolutional layers to be able to model local correlations without requiring the data to lie on a grid. We apply those layers in SchNet: a novel deep learning architecture modeling quantum interactions in molecules. We obtain a joint model for the total energy and interatomic forces that follows fundamental quantum-chemical principles. This includes rotationally invariant energy predictions and a smooth, differentiable potential energy surface. Our architecture achieves state-of-the-art performance for benchmarks of equilibrium molecules and molecular dynamics trajectories. Finally, we introduce a more challenging benchmark with chemical and structural variations that suggests the path for further work.


title: Empirical priors and posterior concentration rates for a monotone density

abstract: In a Bayesian context, prior specification for inference on monotone densities is conceptually straightforward, but proving posterior convergence theorems is complicated by the fact that desirable prior concentration properties often are not satisfied. In this paper, I first develop a new prior designed specifically to satisfy an empirical version of the prior concentration property, and then I give sufficient conditions on the prior inputs such that the corresponding empirical Bayes posterior concentrates around the true monotone density at nearly the optimal minimax rate. Numerical illustrations also reveal the practical benefits of the proposed empirical Bayes approach compared to Dirichlet process mixtures.


title: Neural Question Answering at BioASQ 5B

abstract: This paper describes our submission to the 2017 BioASQ challenge. We participated in Task B, Phase B which is concerned with biomedical question answering (QA). We focus on factoid and list question, using an extractive QA model, that is, we restrict our system to output substrings of the provided text snippets. At the core of our system, we use FastQA, a state-of-the-art neural QA system. We extended it with biomedical word embeddings and changed its answer layer to be able to answer list questions in addition to factoid questions. We pre-trained the model on a large-scale open-domain QA dataset, SQuAD, and then fine-tuned the parameters on the BioASQ training set. With our approach, we achieve state-of-the-art results on factoid questions and competitive results on list questions.


title: Parareal Algorithm Implementation and Simulation in Julia

abstract: We present a full implementation of the parareal algorithm---an integration technique to solve differential equations in parallel---in the Julia programming language for a fully general, first-order, initial-value problem. We provide a brief overview of Julia---a concurrent programming language for scientific computing. Our implementation of the parareal algorithm accepts both coarse and fine integrators as functional arguments. We use Euler's method and another Runge-Kutta integration technique as the integrators in our experiments. We also present a simulation of the algorithm for purposes of pedagogy and as a tool for investigating the performance of the algorithm.


title: Distance-like functions and smooth approximations: a correction to "Logarithm laws for flows on homogeneous spaces"

abstract: One of the propositions in the paper [D. Kleinbock and G.A. Margulis, Logarithm laws for flows on homogeneous spaces, Invent. Math. 138 (1999), 451-494] related to approximating certain sets by smooth functions, was recently found to be incorrect. Here we correct the mistake.


title: Real-time measurements of dissipative solitons in a mode-locked fiber laser

abstract: Dissipative solitons are remarkable localized states of a physical system that arise from the dynamical balance between nonlinearity, dispersion and environmental energy exchange. They are the most universal form of soliton that can exist in nature, and are seen in far-from-equilibrium systems in many fields including chemistry, biology, and physics. There has been particular interest in studying their properties in mode-locked lasers producing ultrashort light pulses, but experiments have been limited by the lack of convenient measurement techniques able to track the soliton evolution in real-time. Here, we use dispersive Fourier transform and time lens measurements to simultaneously measure real-time spectral and temporal evolution of dissipative solitons in a fiber laser as the turn-on dynamics pass through a transient unstable regime with complex break-up and collision dynamics before stabilizing to a regular mode-locked pulse train. Our measurements enable reconstruction of the soliton amplitude and phase and calculation of the corresponding complex-valued eigenvalue spectrum to provide further physical insight. These findings are significant in showing how real-time measurements can provide new perspectives into the ultrafast transient dynamics of complex systems.


title: Analytic Moduli of Plane Branches and Holomorphic Flows

abstract: We study the behaviour (in the infinitesimal neighbourhood of the singularity) of a singular plane branch under the action of holomorphic flows. The techniques we develop provide a new elementary, geometric and dynamical solution to Zariski's moduli problem for singular branches in $({\mathbb C}^{2},0)$. Furthermore, we study whether elements of the same class of analytic conjugacy are conjugated by a holomorphic flow; in particular we show that there exists an analytic class that is not complete: meaning that there are two elements of the class that are not analytically conjugated by a local diffeomorphism embedded in a one-parameter flow.


title: Oxide-mediated self-limiting recovery of field effect mobility in plasma-treated MoS$_2$

abstract: Precise tunability of electronic properties of 2D nanomaterials is a key goal of current research in this field of materials science. Chemical modification of layered transition metal dichalcogenides leads to the creation of heterostructures of low-dimensional variants of these materials. In particular, the effect of oxygen-containing plasma treatment on molybdenum disulfide (MoS$_2$) has long been thought to be detrimental to the electrical performance of the material. Here we show that the mobility and conductivity of MoS$_2$ can be precisely controlled and improved by systematic exposure to oxygen:argon plasma, and characterise the material utilising advanced spectroscopy and microscopy. Through complementary theoretical modelling which confirms conductivity enhancement, we uncover the role of a two-dimensional phase of molybdenum trioxide (2D-MoO$_3$) in improving the electronic behaviour of the material. Deduction of the beneficial role of MoO$_3$ will serve to open the field to new approaches with regard to the tunability of 2D semiconductors by their low-dimensional oxides in nano-modified heterostructures.


title: Detecting Small Signs from Large Images

abstract: In the past decade, Convolutional Neural Networks (CNNs) have been demonstrated successful for object detections. However, the size of network input is limited by the amount of memory available on GPUs. Moreover, performance degrades when detecting small objects. To alleviate the memory usage and improve the performance of detecting small traffic signs, we proposed an approach for detecting small traffic signs from large images under real world conditions. In particular, large images are broken into small patches as input to a Small-Object-Sensitive-CNN (SOS-CNN) modified from a Single Shot Multibox Detector (SSD) framework with a VGG-16 network as the base network to produce patch-level object detection results. Scale invariance is achieved by applying the SOS-CNN on an image pyramid. Then, image-level object detection is obtained by projecting all the patch-level detection results to the image at the original scale. Experimental results on a real-world conditioned traffic sign dataset have demonstrated the effectiveness of the proposed method in terms of detection accuracy and recall, especially for those with small sizes.


title: Using Frame Theoretic Convolutional Gridding for Robust Synthetic Aperture Sonar Imaging

abstract: Recent progress in synthetic aperture sonar (SAS) technology and processing has led to significant advances in underwater imaging, outperforming previously common approaches in both accuracy and efficiency. There are, however, inherent limitations to current SAS reconstruction methodology. In particular, popular and efficient Fourier domain SAS methods require a 2D interpolation which is often ill conditioned and inaccurate, inevitably reducing robustness with regard to speckle and inaccurate sound-speed estimation. To overcome these issues, we propose using the frame theoretic convolution gridding (FTCG) algorithm to handle the non-uniform Fourier data. FTCG extends upon non-uniform fast Fourier transform (NUFFT) algorithms by casting the NUFFT as an approximation problem given Fourier frame data. The FTCG has been show to yield improved accuracy at little more computational cost. Using simulated data, we outline how the FTCG can be used to enhance current SAS processing.


title: Invariant Causal Prediction for Nonlinear Models

abstract: An important problem in many domains is to predict how a system will respond to interventions. This task is inherently linked to estimating the system's underlying causal structure. To this end, Invariant Causal Prediction (ICP) (Peters et al., 2016) has been proposed which learns a causal model exploiting the invariance of causal relations using data from different environments. When considering linear models, the implementation of ICP is relatively straightforward. However, the nonlinear case is more challenging due to the difficulty of performing nonparametric tests for conditional independence. In this work, we present and evaluate an array of methods for nonlinear and nonparametric versions of ICP for learning the causal parents of given target variables. We find that an approach which first fits a nonlinear model with data pooled over all environments and then tests for differences between the residual distributions across environments is quite robust across a large variety of simulation settings. We call this procedure "invariant residual distribution test". In general, we observe that the performance of all approaches is critically dependent on the true (unknown) causal structure and it becomes challenging to achieve high power if the parental set includes more than two variables. As a real-world example, we consider fertility rate modelling which is central to world population projections. We explore predicting the effect of hypothetical interventions using the accepted models from nonlinear ICP. The results reaffirm the previously observed central causal role of child mortality rates.


title: Incoherent qubit control using the quantum Zeno effect

abstract: The quantum Zeno effect is the suppression of Hamiltonian evolution by repeated observation, resulting in the pinning of the state to an eigenstate of the measurement observable. Using measurement only, control of the state can be achieved if the observable is slowly varied such that the state tracks the now time-dependent eigenstate. We demonstrate this using a circuit-QED readout technique that couples to a dynamically controllable observable of a qubit. Continuous monitoring of the measurement record allows us to detect an escape from the eigenstate, thus serving as a built-in form of error detection. We show this by post-selecting on realizations with arbitrarily high fidelity with respect to the target state. Our dynamical measurement operator technique offers a new tool for numerous forms of quantum feedback protocols, including adaptive measurements and rapid state purification.


title: Search for a heavy composite Majorana neutrino in the final state with two leptons and two quarks at sqrt(s) = 13 TeV

abstract: A search for physics beyond the standard model in the final state with two same-flavour leptons (electrons or muons) and two quarks produced in proton-proton collisions at sqrt(s) = 13 TeV is presented. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 2.3 inverse femtobarns. The observed data are in good agreement with the standard model background prediction. The results of the measurement are interpreted in the framework of a recently proposed model in which a heavy Majorana neutrino, N(l), stems from a composite-fermion scenario. Exclusion limits are set for the first time on the mass of the heavy composite Majorana neutrino, m[N(l)], and the compositeness scale Lambda. For the case m[N(l)] = Lambda, the existence of N(e) (N(mu)) is excluded for masses up to 4.60 (4.70) TeV at 95% confidence level.


title: Numerical Simulations of Collisional Cascades at the Roche Limits of White Dwarf Stars

abstract: We consider the long-term collisional and dynamical evolution of solid material orbiting in a narrow annulus near the Roche limit of a white dwarf. With orbital velocities of 300 km/sec, systems of solids with initial eccentricity $e \gtrsim 10^{-3}$ generate a collisional cascade where objects with radii $r \lesssim$ 100--300 km are ground to dust. This process converts 1-100 km asteroids into 1 $\mu$m particles in $10^2 - 10^6$ yr. Throughout this evolution, the swarm maintains an initially large vertical scale height $H$. Adding solids at a rate $\dot{M}$ enables the system to find an equilibrium where the mass in solids is roughly constant. This equilibrium depends on $\dot{M}$ and $r_0$, the radius of the largest solid added to the swarm. When $r_0 \lesssim$ 10 km, this equilibrium is stable. For larger $r_0$, the mass oscillates between high and low states; the fraction of time spent in high states ranges from 100% for large $\dot{M}$ to much less than 1% for small $\dot{M}$. During high states, the stellar luminosity reprocessed by the solids is comparable to the excess infrared emission observed in many metallic line white dwarfs.


title: Learning Local Feature Aggregation Functions with Backpropagation

abstract: This paper introduces a family of local feature aggregation functions and a novel method to estimate their parameters, such that they generate optimal representations for classification (or any task that can be expressed as a cost function minimization problem). To achieve that, we compose the local feature aggregation function with the classifier cost function and we backpropagate the gradient of this cost function in order to update the local feature aggregation function parameters. Experiments on synthetic datasets indicate that our method discovers parameters that model the class-relevant information in addition to the local feature space. Further experiments on a variety of motion and visual descriptors, both on image and video datasets, show that our method outperforms other state-of-the-art local feature aggregation functions, such as Bag of Words, Fisher Vectors and VLAD, by a large margin.


title: Treewidth Bounds for Planar Graphs Using Three-Sided Brambles

abstract: Square grids play a pivotal role in Robertson and Seymour's work on graph minors as planar obstructions to small treewidth. We introduce a three-sided bramble in a plane graph called a net, which generalizes the standard bramble of crosses in a square grid. We then characterize any minimal cover of a net as a tree drawn in the plane. We use nets in an $O(n^3)$ time algorithm that computes both upper and lower bounds on the bramble number (hence treewidth) of any planar graph. Let $G$ be a planar graph, $BN(G)$ be its bramble number and $\lambda(G)$ be the largest order of any net in a subgraph of $G$. Our algorithm outputs a constant, $KB$, so that $\lambda(G)/4 \leq KB \leq BN(G)\leq 4KB \leq 4\lambda(G)$. Let $s(G)$ be the size of a side of the largest square grid minor of $G$. Smith (2015) has shown that $\lambda(G) \geq s(G)$. Our upper bound improves that of Grigoriev (2011) when $\lambda(G)\leq (5/4)s(G)$. We correct a lower bound of Bodlaender, Grigoriev and Koster (2008) to $s(G)/5$ (instead of $s(G)/4$) and thus the lower bound of $\lambda(G)/4$ on our approximation is an improvement.


title: Similarity of operators on $l^{p}$

abstract: For $1<p<\infty$, we prove (i) a version of Voiculescu's absorption theorem for operators on $l^{p}$, (ii) that $\mathrm{Ext}_{\sim,s}(\mathcal{A},K(l^{p}))$ is a group for certain Banach algebra $\mathcal{A}$, and (iii) homotopy invariance of $\mathrm{Ext}_{\sim,s}(\mathcal{A},K(l^{p}))^{-1}$ in $\mathcal{A}$ for separable Banach algebra $\mathcal{A}$ that is isomorphic to a subalgebra of $B(l^{p})$.


title: Formation of a transient amorphous solid in low density aqueous charged sphere suspensions

abstract: Colloidal glasses form from hard spheres, nearly hard spheres, ellipsoids and platelets or their attractive variants have been studied in detail. Complementing and checking theoretical approaches and simulations, the many different types of model systems have significantly advanced our understanding of the glass transition in general. Despite their early prediction, however, no experimental charged sphere glasses have been found at low density, where the competing process of crystallization prevails. We here report the formation of a transient amorphous solid formed from charged polymer spheres suspended in thoroughly deionized water at volume fractions of 0.0002-0.01. From optical experiments, we observe the presence of short-range order and an enhanced shear rigidity as compared to the stable polycrystalline solid of body centred cubic structure. On a density dependent time scale of hours to days, the amorphous solid transforms into this stable structure. We further present preliminary dynamic light scattering data showing the evolution of a second slow relaxation process possibly pointing to a dynamic heterogeneity known from other colloidal glasses and gels.We compare our findings to the predicted phase behaviour of charged sphere suspensions and discuss possible mechanisms for the formation of this peculiar type of colloidal glass.


title: Learning to make external sensory stimulus predictions using internal correlations in populations of neurons

abstract: To compensate for sensory processing delays, the visual system must make predictions to ensure timely and appropriate behaviors. Recent work has found predictive information about the stimulus in neural populations early in vision processing, starting in the retina. However, to utilize this information, cells downstream must in turn be able to read out the predictive information from the spiking activity of retinal ganglion cells. Here we investigate whether a downstream cell could learn efficient encoding of predictive information in its inputs in the absence of other instructive signals, from the correlations in the inputs themselves. We simulate learning driven by spiking activity recorded in salamander retina. We model a downstream cell as a binary neuron receiving a small group of weighted inputs and quantify the predictive information between activity in the binary neuron and future input. Input weights change according to spike timing-dependent learning rules during a training period. We characterize the readouts learned under spike timing-dependent learning rules, finding that although the fixed points of learning dynamics are not associated with absolute optimal readouts, they convey nearly all the information conveyed by the optimal readout. Moreover, we find that learned perceptrons transmit position and velocity information of a moving bar stimulus nearly as efficiently as optimal perceptrons. We conclude that predictive information is, in principle, readable from the perspective of downstream neurons in the absence of other inputs, and consequently suggests that bottom-up prediction may play an important role in sensory processing.


title: Rapidity distribution of particle multiplicity in DIS at small x

abstract: Analytical study of the rapidity distribution of the final state particles in deep inelastic scattering at small x is presented. We separate and analyse three sources of particle production: fragmentation of the quark-antiquark pair, accompanying coherent soft gluon radiation due to octet color exchange in the t-channel, and fragmentation of gluons that form parton distribution functions. Connection to Catani-Ciafaloni-Fiorani-Marchesini (CCFM) equations and the role of gluon reggezation are also discussed.


title: Supersonic Cluster Beam Fabrication of a Stretchable Keyboard with Multilayer Electrical Interconnects Integrated with Paper-based Flexible Devices

abstract: The development of stretchable and flexible electronics based on soft polymeric substrates requires novel approaches in the design, fabrication, integration, and packaging of passive and active electronic components. Here we demonstrate the fabrication, based on the use of supersonic cluster beam implantation and deposition, of a system consisting of a stretchable keyboard made by PDMS connected with a touch sensor and a matrix of LEDs printed on paper, through paper-based interconnections. The stretchable keyboard has conductive vias to electrically connect the top and bottom layers. Gold conductive films printed on paper are used to connect the keyboard with the electronic components that drive the LED matrix.


title: On structural properties of trees with minimal atom-bond connectivity index IV: Solving a conjecture about the pendent paths of length three

abstract: The atom-bond connectivity (ABC) index is one of the most investigated degree-based molecular structure descriptors with a variety of chemical applications. It is known that among all connected graphs, the trees minimize the ABC index. However, a full characterization of trees with a minimal ABC index is still an open problem. By now, one of the proved properties is that a tree with a minimal ABC index may have, at most, one pendent path of length $3$, with the conjecture that it cannot be a case if the order of a tree is larger than $1178$. Here, we provide an affirmative answer of a strengthened version of that conjecture, showing that a tree with minimal ABC index cannot contain a pendent path of length $3$ if its order is larger than $415$.


title: Corrigendum for "Second-order reflected backward stochastic differential equations" and "Second-order BSDEs with general reflection and game options under uncertainty"

abstract: The aim of this short note is to fill in a gap in our earlier paper [16] on 2BSDEs with reflections, and to explain how to correct the subsequent results in the second paper [15]. We also provide more insight on the properties of 2RBSDEs, in the light of the recent contributions [13, 23] in the so--called $G-$framework.


title: Shallow water modeling of rolling pad instability in liquid metal batteries

abstract: Magnetohydrodynamically induced interface instability in liquid metal batteries is analyzed. The batteries are represented by a simplified system in the form of a rectangular cell, in which strong vertical electric current flows through three horizontal layers: the layer of a heavy metal at the bottom, the layer of a light metal at the top, and the layer of electrolyte in the middle. A new two-dimensional nonlinear model based on the conservative shallow water approximation is derived and utilized in a numerical study. It is found that in the case of small density difference between the electrolyte and one of the metals, the instability closely resembles the rolling pad instability observed earlier in the aluminum reduction cells. When the two electrolyte-metal density differences are comparable, the dynamics of unstable systems is more complex and characterized by interaction between two nearly symmetric or antisymmetric interfacial waves.


title: Robust Sonar ATR Through Bayesian Pose Corrected Sparse Classification

abstract: Sonar imaging has seen vast improvements over the last few decades due in part to advances in synthetic aperture Sonar (SAS). Sophisticated classification techniques can now be used in Sonar automatic target recognition (ATR) to locate mines and other threatening objects. Among the most promising of these methods is sparse reconstruction-based classification (SRC) which has shown an impressive resiliency to noise, blur, and occlusion. We present a coherent strategy for expanding upon SRC for Sonar ATR that retains SRC's robustness while also being able to handle targets with diverse geometric arrangements, bothersome Rayleigh noise, and unavoidable background clutter. Our method, pose corrected sparsity (PCS), incorporates a novel interpretation of a spike and slab probability distribution towards use as a Bayesian prior for class-specific discrimination in combination with a dictionary learning scheme for localized patch extractions. Additionally, PCS offers the potential for anomaly detection in order to avoid false identifications of tested objects from outside the training set with no additional training required. Compelling results are shown using a database provided by the United States Naval Surface Warfare Center.


title: Simulation optimization: A review of algorithms and applications

abstract: Simulation Optimization (SO) refers to the optimization of an objective function subject to constraints, both of which can be evaluated through a stochastic simulation. To address specific features of a particular simulation---discrete or continuous decisions, expensive or cheap simulations, single or multiple outputs, homogeneous or heterogeneous noise---various algorithms have been proposed in the literature. As one can imagine, there exist several competing algorithms for each of these classes of problems. This document emphasizes the difficulties in simulation optimization as compared to mathematical programming, makes reference to state-of-the-art algorithms in the field, examines and contrasts the different approaches used, reviews some of the diverse applications that have been tackled by these methods, and speculates on future directions in the field.


title: Profilometry with digital fringe-projection at the spatial and temporal Nyquist frequencies

abstract: A phase-demodulation method for digital fringe-projection profilometry using the spatial and temporal Nyquist frequencies is presented. It allows us to digitize tridimensional surfaces using the highest spatial frequency ({\pi} radians per pixel) and consequently with the highest sensitivity for a given digital fringe projector. Working with the highest temporal frequency ({\pi} radians per temporal sample), the proposed method rejects the DC component and all even-order distorting harmonics using a bare-minimum 2-step phase shift. The proposed method is suitable for digitization of piece-wise continuous surfaces because it does not require spatial low-pass filtering. Gamma calibration is unnecessary because the projected fringes are binary, and the harmonics produced by the binary profile can be easily attenuated with a slight defocusing on the digital projector. Viability of the proposed method is supported by experimental results showing complete agreement with the predicted behavior.


title: Two-dimensional Klein-Gordon Oscillator in the presence of a minimal length

abstract: Minimal length of a two-dimensional Klein-Gordon oscillator is investigated and illustrates the wave functions in the momentum space. The energy eigenvalues are found and the corresponding wave functions are calculated in terms of hyper-geometric functions.


title: On locally compact topological graph inverse semigroups

abstract: In this paper we characterise graph inverse semigroups which admit only discrete locally compact semigroup topology. This characterization provides a complete answer on the question of Z. Mesyan, J. D. Mitchell, M. Morayne and Y. H. P\'eresse posed in their paper "Topological graph inverse semigroups".


title: On general features of warm dark matter with reduced relativistic gas

abstract: Reduced Relativistic Gas (RRG) is a useful approach to describe the warm dark matter (WDM) or the warmness of baryonic matter in the approximation when the interaction between the particles is irrelevant. The use of Maxwell distribution leads to the complicated equation of state of the J\"{u}ttner model of relativistic ideal gas. The RRG enables one to reproduce the same physical situation but in a much simpler form. For this reason RRG can be a useful tool for the theories with some sort of a "new Physics". On the other hand, even without the qualitatively new physical implementations, the RRG can be useful to describe the general features of WDM in a model-independent way. In this sense one can see, in particular, to which extent the cosmological manifestations of WDM may be dependent on its Particle Physics background. In the present work RRG is used as a complementary approach to derive the main observational exponents for the WDM in a model-independent way. The only assumption concerns a non-negligible velocity $v$ for dark matter particles which is parameterized by the warmness parameter $b$. The relatively high values of $b$ ( $b^2\gtrsim 10^{-6}$) erase the radiation (photons and neutrinos) dominated epoch and cause an early warm matter domination after inflation. Furthermore, RRG approach enables one to quantify the lack of power in linear matter spectrum at small scales and in particular, reproduces the relative transfer function commonly used in context of WDM with accuracy of $\lesssim 1\%$. A warmness with $b^2\lesssim 10^{-6}$ (equivalent to $v\lesssim 300 km/s$) does not alter significantly the CMB power spectrum and is in agreement with the background observational tests.


title: Fractional Driven Damped Oscillator

abstract: The resonances associated with a fractional damped oscillator which is driven by an oscillatory external force are studied. It is shown that such resonances can be manipulated by tuning up either the coefficient of the fractional damping or the order of the corresponding fractional derivatives.


title: Democratizing Design for Future Computing Platforms

abstract: Information and communications technology can continue to change our world. These advances will partially depend upon designs that synergistically combine software with specialized hardware. Today open-source software incubates rapid software-only innovation. The government can unleash software-hardware innovation with programs to develop open hardware components, tools, and design flows that simplify and reduce the cost of hardware design. Such programs will speed development for startup companies, established industry leaders, education, scientific research, and for government intelligence and defense platforms.


title: Cyber Defense as a Complex Adaptive System: A model-based approach to strategic policy design

abstract: In a world of ever-increasing systems interdependence, effective cybersecurity policy design seems to be one of the most critically understudied elements of our national security strategy. Enterprise cyber technologies are often implemented without much regard to the interactions that occur between humans and the new technology. Furthermore, the interactions that occur between individuals can often have an impact on the newly employed technology as well. Without a rigorous, evidence-based approach to ground an employment strategy and elucidate the emergent organizational needs that will come with the fielding of new cyber capabilities, one is left to speculate on the impact that novel technologies will have on the aggregate functioning of the enterprise. In this paper, we will explore a scenario in which a hypothetical government agency applies a complexity science perspective, supported by agent-based modeling, to more fully understand the impacts of strategic policy decisions. We present a model to explore the socio-technical dynamics of these systems, discuss lessons using this platform, and suggest further research and development.


title: Assortment and Price Optimization Under the Two-Stage Luce model

abstract: This paper studies assortment and pricing optimization problems under the Two-Stage Luce model (2SLM), a discrete choice model introduced by Echenique and Saito (2018) that generalizes the multinomial logit model (MNL). The model employs an utility function as in the the MNL, and a dominance relation between products. When consumers are offered an assortmentS , they first discard all dominated products in S and then select one of the remaining products using the standard MNL. This model may violate the regularity condition, which states that the probability of choosing a product cannot increase if the offer set is enlarged. Therefore, the 2SLM falls outside the large family of discrete choice models based on random utility which contains almost all choice models studied in revenue management. We prove that the assortment problem under the 2SLM is polynomial-time solvable. Moreover, we show that the capacitated assortment optimization problem is NP-hard and but it admits polynomial-time algorithms for the relevant special cases cases where (1) the dominance relation is attractiveness-correlated and (2) its transitive reduction is a forest. The proofs exploit a strong connection between assortments under the 2SLM and independent sets in comparability graphs. Finally, we study the associated joint pricing and assortment problem under this model. First, we show that well known optimal pricing policy for the MNL can be arbitrarily bad. Our main result in this section is the development of an efficient algorithm for this pricing problem. The resulting optimal pricing strategy is simple to describe: it assigns the same price for all products, except for the one with the highest attractiveness and as well as for the one with the lowest attractiveness.


title: Multiband superconductivity in ${\rm BiS_2}$-based layered compounds

abstract: A mean-field treatment is presented of a square lattice two-orbital-model for ${\rm BiS_2}$ taking into account intra- and inter-orbital superconductivity. A rich phase diagram involving both types of superconductivity is presented as a function of the ratio between the couplings of electrons in the same and different orbitals (${\rm \eta= V_{XX}/V_{XY}}$) and electron doping $x$. With the help of a quantity we call orbital-mixing ratio, denoted as $R(\phi)$, the phase diagram is analyzed using a simple and intuitive picture based on how $R(\phi)$ varies as electron doping increases. The predictive power of $R(\phi)$ suggests that it could be a useful tool in qualitatively (or even semi-quantitatively) analyzing multiband superconductivity in BCS-like superconductors.


title: Domain reduction techniques for global NLP and MINLP optimization

abstract: Optimization solvers routinely utilize presolve techniques, including model simplification, reformulation and domain reduction techniques. Domain reduction techniques are especially important in speeding up convergence to the global optimum for challenging nonconvex nonlinear programming (NLP) and mixed-integer nonlinear programming (MINLP) optimization problems. In this work, we survey the various techniques used for domain reduction of NLP and MINLP optimization problems. We also present a computational analysis of the impact of these techniques on the performance of various widely available global solvers on a collection of 1740 test problems.


title: Second-Order Moment-Closure for Tighter Epidemic Thresholds

abstract: In this paper, we study the dynamics of contagious spreading processes taking place in complex contact networks. We specifically present a lower-bound on the decay rate of the number of nodes infected by a susceptible-infected-susceptible (SIS) stochastic spreading process. A precise quantification of this decay rate is crucial for designing efficient strategies to contain epidemic outbreaks. However, existing lower-bounds on the decay rate based on first-order mean-field approximations are often accompanied by a large error resulting in inefficient containment strategies. To overcome this deficiency, we derive a lower-bound based on a second-order moment-closure of the stochastic SIS processes. The proposed second-order bound is theoretically guaranteed to be tighter than existing first-order bounds. We also present various numerical simulations to illustrate how our lower-bound drastically improves the performance of existing first-order lower-bounds in practical scenarios, resulting in more efficient strategies for epidemic containment.


title: Tidal synchronization of an anelastic multi-layered body: Titan's synchronous rotation

abstract: This paper presents one analytical tidal theory for a viscoelastic multi-layered body with an arbitrary number of homogeneous layers. Starting with the static equilibrium figure, modified to include tide and differential rotation, and using the Newtonian creep approach, we find the dynamical equilibrium figure of the deformed body, which allows us to calculate the tidal potential and the forces acting on the tide generating body, as well as the rotation and orbital elements variations. In the particular case of the two-layer model, we study the tidal synchronization when the gravitational coupling and the friction in the interface between the layers is added. For high relaxation factors (low viscosity), the stationary solution of each layer is synchronous with the orbital mean motion (n) when the orbit is circular, but the spin rates increase if the orbital eccentricity increases. For low relaxation factors (high viscosity), as in planetary satellites, if friction remains low, each layer can be trapped in different spin-orbit resonances with frequencies n/2,n,3n/2,... . We apply the theory to Titan. The main results are: i) the rotational constraint does not allow us confirm or reject the existence of a subsurface ocean in Titan; and ii) the crust-atmosphere exchange of angular momentum can be neglected. Using the rotation estimate based on Cassini's observation, we limit the possible value of the shell relaxation factor, when a subsurface ocean is assumed, to 10^-9 Hz, which correspond to a shell's viscosity 10^18 Pa s, depending on the ocean's thickness and viscosity values. In the case in which the ocean does not exist, the maximum shell relaxation factor is one order of magnitude smaller and the corresponding minimum shell's viscosity is one order higher.


title: Metal-Insulator Transition in the Dimerized Organic Conductor $\kappa$-(BEDT-TTF)$_2$Hg(SCN)$_2$Br

abstract: The organic charge-transfer salt $\kappa$-(BEDT-TTF)$_{2}$Hg(SCN)$_{2}$Br is a quasi two-dimensional metal with a half-filled conduction band at ambient conditions. When cooled below $T=80$ K it undergoes a pronounced transition to an insulating phase where the resistivity increases many orders of magnitude. In order to elucidate the nature of this metal-insulator transition we have performed comprehensive transport, dielectric and optical investigations. The findings are compared with other dimerized $\kappa$-(BEDT-TTF) salts, in particular the Cl-analogue, where a charge-order transition takes place at $T_{\rm CO}=30$ K.


title: Developing Bug-Free Machine Learning Systems With Formal Mathematics

abstract: Noisy data, non-convex objectives, model misspecification, and numerical instability can all cause undesired behaviors in machine learning systems. As a result, detecting actual implementation errors can be extremely difficult. We demonstrate a methodology in which developers use an interactive proof assistant to both implement their system and to state a formal theorem defining what it means for their system to be correct. The process of proving this theorem interactively in the proof assistant exposes all implementation errors since any error in the program would cause the proof to fail. As a case study, we implement a new system, Certigrad, for optimizing over stochastic computation graphs, and we generate a formal (i.e. machine-checkable) proof that the gradients sampled by the system are unbiased estimates of the true mathematical gradients. We train a variational autoencoder using Certigrad and find the performance comparable to training the same model in TensorFlow.


title: Cognitive Psychology for Deep Neural Networks: A Shape Bias Case Study

abstract: Deep neural networks (DNNs) have achieved unprecedented performance on a wide range of complex tasks, rapidly outpacing our understanding of the nature of their solutions. This has caused a recent surge of interest in methods for rendering modern neural systems more interpretable. In this work, we propose to address the interpretability problem in modern DNNs using the rich history of problem descriptions, theories and experimental methods developed by cognitive psychologists to study the human mind. To explore the potential value of these tools, we chose a well-established analysis from developmental psychology that explains how children learn word labels for objects, and applied that analysis to DNNs. Using datasets of stimuli inspired by the original cognitive psychology experiments, we find that state-of-the-art one shot learning models trained on ImageNet exhibit a similar bias to that observed in humans: they prefer to categorize objects according to shape rather than color. The magnitude of this shape bias varies greatly among architecturally identical, but differently seeded models, and even fluctuates within seeds throughout training, despite nearly equivalent classification performance. These results demonstrate the capability of tools from cognitive psychology for exposing hidden computational properties of DNNs, while concurrently providing us with a computational model for human word learning.


title: Correlated electronic states at domain walls of a Mott-charge-density-wave insulator 1T-TaS2

abstract: Domain walls in interacting electronic systems can have distinct localized states, which often govern physical properties and may lead to unprecedented functionalities and novel devices. However, electronic states within domain walls themselves have not been clearly identified and understood for strongly correlated electron systems. Here, we resolve the electronic states localized on domain walls in a Mott-charge-density-wave(CDW) insulator 1T-TaS2 using scanning tunneling spectroscopy. We establish that the domain wall state decomposes into two nonconducting states located at the center of domain walls and edges of domains. Theoretical calculations reveal their atomistic origin as the local reconstruction of domain walls under the strong influence of electron correlation. Our results introduce a concept for the domain wall electronic property, the wall's own internal degrees of freedom, which is potentially related to the controllability of domain wall electronic properties.


title: Model-checking Counting Temporal Logics on Flat Structures

abstract: We study several extensions of linear-time and computation-tree temporal logics with quantifiers that allow for counting how often certain properties hold. For most of these extensions, the model-checking problem is undecidable, but we show that decidability can be recovered by considering flat Kripke structures where each state belongs to at most one simple loop. Most decision procedures are based on results on (flat) counter systems where counters are used to implement the evaluation of counting operators.