title: Estimation of BBR shift due to Stark effect for the Microwave $^{113}$Cd$^{+}$ Ion Clock

abstract: The microwave clock frequency of the $|5s~^2S_{1/2}, F=0,m_F=0 \rangle \leftrightarrow |5s~^2S_{1/2}, F=1,m_F=0 \rangle$ transition in the $^{113}$Cd$^+$ ion has been reported as 15199862855.0192(10) Hz [Opt. Lett. {\bf 40}, 4249 (2015)]. Fractional systematic due to the black-body radiation (BBR) shift ($\beta$) arising from the Stark effect in the above clock transition was used as $-1.1 \times 10^{-16}$ from our unpublished preliminary estimation. We present here a precise value as $\beta=-1.815(77) \times 10^{-16}$ by carrying out rigorous calculations of third-order polarizabilities of the hyperfine levels associated with the clock transition. This is determined by evaluating matrix elements of the magnetic dipole hyperfine interaction Hamiltonian, electric dipole operator and energies between many low-lying states of $^{113}$Cd$^+$. We employ all-order relativistic many-body theories in the frameworks of Fock-space coupled-cluster and relativistic multi-configuration Dirac-Fock methods.


title: Construction of algebraic covers

abstract: Let $Y$ be an algebraic variety, $\mathcal{F}$ a locally free sheaf of $\mathcal{O}_Y$-modules, and $\mathcal{R}(\mathcal{F})$ the $\mathcal{O}_Y$-algebra $\operatorname{Sym}^\bullet \mathcal{F}$. In this paper we study local properties of sheaves of $\mathcal{O}_{\mathcal{R}(\mathcal{F})}$-ideals $\mathcal{I}$ such that $\mathcal{R}(\mathcal{F}))/\mathcal{I}$ is an algebraic cover of $Y$. Following the work of Miranda for triple covers, for $\mathcal{Q}$ a direct summand of $\mathcal{R}(\mathcal{F})$, we say that a morphism $\Phi\colon \mathcal{Q}\rightarrow\mathcal{R}(\mathcal{F})/\langle\mathcal{Q}\rangle$ is a covering homomorphism if it induces such an ideal. As an application we study in detail the case of Gorenstein covering maps of degree $6$ for which the direct image of $\varphi_*\mathcal{O}_X$ admits an orthogonal decomposition. These are deformation of $S_3$-Galois branch covers.


title: Optimal non-asymptotic bound of the Ruppert-Polyak averaging without strong convexity

abstract: This paper is devoted to the non-asymptotic control of the mean-squared error for the Ruppert-Polyak stochastic averaged gradient descent introduced in the seminal contributions of [Rup88] and [PJ92]. In our main results, we establish non-asymptotic tight bounds (optimal with respect to the Cramer-Rao lower bound) in a very general framework that includes the uniformly strongly convex case as well as the one where the function f to be minimized satisfies a weaker Kurdyka-Lojiasewicz-type condition [Loj63, Kur98]. In particular, it makes it possible to recover some pathological examples such as on-line learning for logistic regression (see [Bac14]) and recursive quan- tile estimation (an even non-convex situation).


title: Manipulating magnetism by ultrafast control of the exchange interaction

abstract: In recent years, the optical control of exchange interactions has emerged as an exciting new direction in the study of the ultrafast optical control of magnetic order. Here we review recent theoretical works on antiferromagnetic systems, devoted to i) simulating the ultrafast control of exchange interactions, ii) modeling the strongly nonequilibrium response of the magnetic order and iii) the relation with relevant experimental works developed in parallel. In addition to the excitation of spin precession, we discuss examples of rapid cooling and the control of ultrafast coherent longitudinal spin dynamics in response to femtosecond optically induced perturbations of exchange interactions. These elucidate the potential for exploiting the control of exchange interactions to find new scenarios for both faster and more energy-efficient manipulation of magnetism.


title: A damped forward EMI model for a horizontally stratified earth

abstract: If a magnetic dipole is placed above the surface of the earth, the Electromagnetic Induction (EMI) effect, encoded in Maxwell's equations, causes eddy currents in the soil which, on their turn, induce response electromagnetic fields. The magnetic field can be measured in geophysical surveys to determine the conductivity profile of the ground in a non-destructive manner. The forward model used in the inversion of experimental data usually consists of a set of horizontal homogeneous layers. A frequently used analytical model, proposed by McNeill, does not include the interaction between the eddy currents, and therefore fails for larger conductivities. In this paper we construct a new forward, analytical, model to estimate the magnetic field caused by a horizontally stratified earth but which approximates the interaction between eddy currents. This makes it valid for a broader range of parameters than the current state of the art. Furthermore, the error with the (numerically obtainable) exact result is substantially decreased. We also calculate the vertical sensitivity ("depth of exploration") of the model and observe that it is in good agreement with the values obtained from the exact model.


title: The origin and properties of massive prolate galaxies in the Illustris simulation

abstract: We study galaxy shapes in the Illustris cosmological hydrodynamic simulation. We find that massive galaxies have a higher probability of being prolate. For galaxies with stellar mass larger than $10^{11}\rm M_{\odot}$, 35 out of total 839 galaxies are prolate. For 21 galaxies with stellar mass larger than $10^{12}\rm M_{\odot}$, 9 are prolate, 4 are triaxial while the others are close to being oblate. There are almost no prolate galaxies with stellar mass smaller than $3\times10^{11}\rm M_{\odot}$. We check the merger history of the prolate galaxies, and find that they are formed by major dry mergers. All the prolate galaxies have at least one such merger, with most having mass ratios between $1:1$ and $1:3$. The gas fraction (gas mass to total baryon mass) of the progenitors is 0-3 percent for nearly all these mergers, except for one whose second progenitor contains $\sim 15\%$ gas mass, while its main progenitor still contains less than $5\%$. For the 35 massive prolate galaxies that we find, 18 of them have minor axis rotation, and their angular momenta mostly come from the spin angular momenta of the progenitors (usually that of the main progenitor). We analyse the merger orbits of these prolate galaxies and find that most of them experienced a nearly radial merger orbit. Oblate galaxies with major dry mergers can have either radial or circular merger orbits. We further discuss various properties of these prolate galaxies, such as spin parameter $\lambda_{\rm R}$, spherical anisotropy parameter $\beta$, dark matter fraction, as well as inner density slopes for the stellar, dark matter and total mass distributions.


title: Simultaneous Dimension Reduction and Clustering via the NMF-EM Algorithm

abstract: Mixture models are among the most popular tools for clustering. However, when the dimension and the number of clusters is large, the estimation of the clusters become challenging, as well as their interpretation. Restriction on the parameters can be used to reduce the dimension. An example is given by mixture of factor analyzers for Gaussian mixtures. The extension of MFA to non-Gaussian mixtures is not straightforward. We propose a new constraint for parameters in non-Gaussian mixture model: the $K$ components parameters are combinations of elements from a small dictionary, say $H$ elements, with $H \ll K$. Including a nonnegative matrix factorization (NMF) in the EM algorithm allows us to simultaneously estimate the dictionary and the parameters of the mixture. We propose the acronym NMF-EM for this algorithm, implemented in the R package {\tt nmfem}. This original approach is motivated by passengers clustering from ticketing data: we apply NMF-EM to data from two Transdev public transport networks. In this case, the words are easily interpreted as typical slots in a timetable.


title: Comparison of damping mechanisms for transverse waves in solar coronal loops

abstract: We present a method to assess the plausibility of alternative mechanisms to explain the damping of magnetohydrodynamic (MHD) transverse waves in solar coronal loops. The considered mechanisms are resonant absorption of kink waves in the Alfv\'en continuum, phase-mixing of Alfv\'en waves, and wave leakage. Our methods make use of Bayesian inference and model comparison techniques. We first infer the values for the physical parameters that control the wave damping, under the assumption of a particular mechanism, for typically observed damping time-scales. Then, the computation of marginal likelihoods and Bayes factors enable us to quantify the relative plausibility between the alternative mechanisms. We find that, in general, the evidence is not large enough to support a single particular damping mechanism as the most plausible one. Resonant absorption and wave leakage offer the most probable explanations in strong damping regimes, while phase mixing is the best candidate for weak/moderate damping. When applied to a selection of 89 observed transverse loop oscillations, with their corresponding measurements of damping times scales and taking into account data uncertainties, we find that only in a few cases positive evidence for a given damping mechanism is available.


title: Relativity, Anomalies and Objectivity Loophole in Recent Tests of Local Realism

abstract: Local realism is in conflict with special quantum Bell-type models. Recently, several experiments have demonstrated violation of local realism if we trust their setup assuming special relativity valid. In this paper we question the assumption of relativity, point out not commented anomalies and show that the experiments have not closed objectivity loophole because clonability of the result has not been demonstrated. We propose several improvements in further experimental tests of local realism make the violation more convincing.


title: Topological properties of chains of magnetic impurities on a superconducting substrate: Interplay between the Shiba band and ferromagnetic wire limits

abstract: We consider a one-dimensional system combining local magnetic moments and a delocalized metallic band on top of a superconducting substrate. This system can describe a chain of magnetic impurities with both localized polarized orbitals and delocalized s-like orbitals or a conducting wire with embedded magnetic impurities. We study the interplay between the one-dimensional Shiba band physics arising from the interplay between magnetic moments and the substrate and the delocalized wire-like conduction band on top of the superconductor. We derive an effective low-energy Hamiltonian in terms of two coupled asymmetric Kitaev-like Hamiltonians and analyze its topological properties. We have found that this system can host multiple Majorana bound states at its extremities provided a magnetic mirror symmetry is present. We compute the phase diagram of the system depending on the magnetic exchange interactions, the impurity distance and especially the coupling between both bands. In presence of inhomogeneities which typically break this magnetic mirror symmetry, we show that the coexistence of a Shiba and wire delocalized topological band can drive the system into a non-topological regime with a splitting of Majorana bound states.


title: Strong convergence of the Euler--Maruyama approximation for a class of L\'evy-driven SDEs

abstract: Consider the following stochastic differential equation (SDE) $$dX_t = b(t,X_{t-}) \, dt+ dL_t, \quad X_0 = x,$$ driven by a $d$-dimensional L\'evy process $(L_t)_{t \geq 0}$. We establish conditions on the L\'evy process and the drift coefficient $b$ such that the Euler--Maruyama approximation converges strongly to a solution of the SDE with an explicitly given rate. The convergence rate depends on the regularity of $b$ and the behaviour of the L\'evy measure at the origin. As a by-product of the proof, we obtain that the SDE has a pathwise unique solution. Our result covers many important examples of L\'evy processes, e.g. isotropic stable, relativistic stable, tempered stable and layered stable.


abstract: A Fourier transform from momentum space to twistor space is introduced in twistor string theory, for the first time, for the case where the twistor space is a three-dimensional real projective space, corresponding to ultra-hyperbolic spacetime. In this case, the Fourier transform is the same as in the standard analysis. However, when the twistor space is a three-dimensional complex projective space, corresponding to the complexified Minkowski space, some aspects of the Fourier transform have yet to be clarified. For example, no concrete method is known for calculating the complex integral, and nor is it known which functions it should be applied to or what results are obtained. In this paper, we define a Fourier transform for momentum-space functions in terms of a certain complex integral, assuming that the functions can be expanded as power series. We also define the inverse Fourier transform for twistor-space functions in terms of a certain complex integral. Then, we show that the momentum-space functions to which the Fourier transform can be applied and the twistor-space functions that are obtained can be represented in terms of cohomology groups. Furthermore, we show that the twistor operator representations popular in the twistor theory literature are valid for the Fourier transform and its inverse transform. We also show that the functions over which the application of the operators is closed can be represented in terms of cohomology groups.


title: The Ramsey-Tur\'{a}n problem for cliques

abstract: An important question in extremal graph theory raised by Vera T. S\'os asks to determine for a given integer $t\ge 3$ and a given positive real number $\delta$ the asymptotically supremal edge density $f_t(\delta)$ that an $n$-vertex graph can have provided it contains neither a complete graph $K_t$ nor an independent set of size $\delta n$. Building upon recent work of Fox, Loh, and Zhao [The critical window for the classical Ramsey-Tur\'an problem, Combinatorica 35 (2015), 435-476], we prove that if $\delta$ is sufficiently small (in a sense depending on $t$), then \[ f_t(\delta)= \begin{cases} \frac{3t-10}{3t-4}+\delta-\delta^2 & \text{ if $t$ is even,} \cr \frac{t-3}{t-1}+\delta & \text{ if $t$ is odd.} \end{cases} \]


title: Optimal verification of entangled states with local measurements

abstract: Consider the task of verifying that a given quantum device, designed to produce a particular entangled state, does indeed produce that state. One natural approach would be to characterise the output state by quantum state tomography; or alternatively to perform some kind of Bell test, tailored to the state of interest. We show here that neither approach is optimal amongst local verification strategies for two qubit states. We find the optimal strategy in this case and show that quadratically fewer total measurements are needed to verify to within a given fidelity than in published results for quantum state tomography, Bell test, or fidelity estimation protocols. We also give efficient verification protocols for any stabilizer state. Additionally, we show that requiring that the strategy be constructed from local, non-adaptive and non-collective measurements only incurs a constant-factor penalty over a strategy without these restrictions.


title: Chromatic symmetric functions and H-free graphs

abstract: Two celebrated conjectures in chromatic symmetric function theory concern the positivity chromatics symmetric functions of claw-free graphs. Here we extend the claw-free idea to general graphs and consider the e-positivity question for H-free graphs where H = {claw, F} and H={claw, F, co-F}, where F is a four-vertex graph. We settle the question for all cases except H={claw, co-diamond}, and we provide some partial results in that case.


title: Evidence for Z=6 `magic number' in neutron-rich carbon isotopes

abstract: The nuclear shell structure, which originates in the nearly independent motion of nucleons in an average potential, provides an important guide for our understanding of nuclear structure and the underlying nuclear forces. Its most remarkable fingerprint is the existence of the so-called `magic numbers' of protons and neutrons associated with extra stability. Although the introduction of a phenomenological spin-orbit (SO) coupling force in 1949 helped explain the nuclear magic numbers, its origins are still open questions. Here, we present experimental evidence for the smallest SO-originated magic number (subshell closure) at the proton number 6 in 13-20C obtained from systematic analysis of point-proton distribution radii, electromagnetic transition rates and atomic masses of light nuclei. Performing ab initio calculations on 14,15C, we show that the observed proton distribution radii and subshell closure can be explained by the state-of-the-art nuclear theory with chiral nucleon-nucleon and three-nucleon forces, which are rooted in the quantum chromodynamics.


title: Local invariants in non-ideal flows of neutral fluids and two-fluid plasmas

abstract: Local structures, beyond the well-known `frozen-in' to the barotropic flows of the generalized vorticities, of the two-fluid model of plasma flows are presented. More general non-barotropic situations are also considered. A modified Euler equation [T. Tao, Ann. PDE \textbf{2}, 9 (2016)] is also accordingly analyzed and remarked from the angle of view of two-fluid model, with emphasis on the local structures. And, the local constraints of high-order differential forms such as helicity, among others, find simple formulation for possible applications in modelling the dynamics. Thus, the Cauchy invariants equation [N. Besse and U. Frisch, J. Fluid Mech. \textbf{825}, 412 (2017)] may find practical application in non-ideal flows. Some formal examples are offered to outline the calculations, and particularly interestingly the two-dimensional-three-component (2D3C) or the 2D passive scalar problem presents that a locally invariant $\Theta = 2\theta \zeta$, with $\theta$ and $\zeta$ being respectively the scalar value of the `vertical velocity' (or the passive scalar) and the `vertical vorticity', may be used as if it were the spatial density of the globally invariant helicity, providing a Lagrangian prescription to control the latter in some situations of studying its physical effects in rapidly rotating flows (ubiquitous in atmosphere of astrophysical objects) with marked 2D3C vortical modes or in purely 2D passive scalars.

categories: physics.flu-dyn astro-ph.SR math-ph math.MP nlin.CD


title: Characterizations of centrality by local convexity of certain functions on $C^*$-algebras

abstract: We provide a function class which is useful to distinguish central and non-central elements of a $C^*$-algebra in the following sense: for each element $f$ of this function class, a self-adjoint element $a$ of a $C^*$-algebra is central if and only if $f$ is locally convex at $a.$


title: qDSA: Small and Secure Digital Signatures with Curve-based Diffie--Hellman Key Pairs

abstract: qDSA is a high-speed, high-security signature scheme that facilitates implementations with a very small memory footprint, a crucial requirement for embedded systems and IoT devices, and that uses the same public keys as modern Diffie--Hellman schemes based on Montgomery curves (such as Curve25519) or Kummer surfaces. qDSA resembles an adaptation of EdDSA to the world of Kummer varieties, which are quotients of algebraic groups by $\pm$1. Interestingly, qDSA does not require any full group operations or point recovery: all computations, including signature verification, occur on the quotient where there is no group law. We include details on four implementations of qDSA, using Montgomery and fast Kummer surface arithmetic on the 8-bit AVR ATmega and 32-bit ARM Cortex M0 platforms. We find that qDSA significantly outperforms state-of-the-art signature implementations in terms of stack usage and code size. We also include an efficient compression algorithm for points on fast Kummer surfaces, reducing them to the same size as compressed elliptic curve points for the same security level.


title: Unboundedness of some higher Euler classes

abstract: In this paper, we study Euler classes in groups of homeomorphisms of Seifert fibered 3-manifolds. We show that, in contrast to the familiar Euler class for $\mathrm{Homeo}_0(S^1)^\delta$, these Euler classes for $\mathrm{Homeo}_0(M^3)^\delta$ are unbounded classes. In fact, we give examples of flat topological M bundles over a genus 3 surface whose Euler class takes arbitrary values.


title: A Multilayer and Time-varying Structural Analysis of the Brazilian Air Transportation Network

abstract: This paper provides a multilayer and time-varying structural analysis of one air transportation network, having the Brazilian air transportation network as a case study. Using a single mathematical object called MultiAspect Graph (MAG) for this analysis, the multilayer perspective enables the unveiling of the particular strategies of each airline to both establish and adapt in a moment of crisis its specific flight network. Similarly, the time-varying perspective allows multi-scale analysis considering different time periods, and thus assessing the impact of the economic crisis on how the different airlines establish their routes as well as the flights that use these routes. Altogether, besides the multilayer and time-varying structural analysis of the Brazilian air transportation network, this paper also acts as a proof-of-concept for the MAG potential for the modeling and analysis of high-order networks.


title: Crossover from impurity-controlled to granular superconductivity in (TMTSF)2ClO4

abstract: Using a proper cooling procedure, a controllable amount of non-magnetic structural disorder can be introduced at low temperature in (TMTSF)2ClO4. Here we performed simultaneous measurements of transport and magnetic properties of (TMTSF)2ClO4 in its normal and superconducting states, while finely covering three orders of magnitude of the cooling rate around the anion ordering temperature. Our result reveals, with increasing density of disorder, the existence of a crossover between homogeneous defect-controlled d-wave superconductivity and granular superconductivity. At slow cooling rates, with small amount of disorder, the evolution of superconducting properties is well described with the Abrikosov-Gorkov theory, providing further confirmation of non-s-wave pairing in this compound. In contrast, at fast cooling rates, zero resistance and diamagnetic shielding are achieved through a randomly distributed network of superconducting puddles embedded in an normal conducting background and interconnected by proximity effect coupling. The temperature dependence of the AC complex susceptibility reveals features typical for a network of granular superconductors. This makes (TMTSF)2ClO4 a model system for granular superconductivity where the grain size and their concentration are tunable within the same sample.


title: Cartan Invariants and Event Horizon Detection, Extended Version

abstract: We show that it is possible to locate the event horizons of a black hole (in arbitrary dimensions) as the zeros of certain Cartan invariants. This approach accounts for the recent results on the detection of stationary horizons using scalar polynomial curvature invariants, and improves upon them since the proposed method is computationally less expensive. As an application, we produce Cartan invariants that locate the event horizons for various exact four-dimensional and five-dimensional stationary, asymptotically flat (or (anti) de Sitter) black hole solutions and compare the Cartan invariants with the corresponding scalar curvature invariants that detect the event horizon. In particular, for each of the four-dimensional examples we express the scalar polynomial curvature invariants introduced by Abdelqader and Lake in terms of the Cartan invariants and show a direct relationship between the scalar polynomial curvature invariants and the Cartan invariants that detect the horizon.


title: A Planning Approach to Monitoring Behavior of Computer Programs

abstract: We describe a novel approach to monitoring high level behaviors using concepts from AI planning. Our goal is to understand what a program is doing based on its system call trace. This ability is particularly important for detecting malware. We approach this problem by building an abstract model of the operating system using the STRIPS planning language, casting system calls as planning operators. Given a system call trace, we simulate the corresponding operators on our model and by observing the properties of the state reached, we learn about the nature of the original program and its behavior. Thus, unlike most statistical detection methods that focus on syntactic features, our approach is semantic in nature. Therefore, it is more robust against obfuscation techniques used by malware that change the outward appearance of the trace but not its effect. We demonstrate the efficacy of our approach by evaluating it on actual system call traces.


title: Detecting localized eigenstates of linear operators

abstract: We describe a way of detecting the location of localized eigenvectors of a linear system $Ax = \lambda x$ for eigenvalues $\lambda$ with $|\lambda|$ comparatively large. We define the family of functions $f_{\alpha}: \left\{1.2. \dots, n\right\} \rightarrow \mathbb{R}_{}$ $$ f_{\alpha}(k) = \log \left( \| A^{\alpha} e_k \|_{\ell^2} \right),$$ where $\alpha \geq 0$ is a parameter and $e_k = (0,0,\dots, 0,1,0, \dots, 0)$ is the $k-$th standard basis vector. We prove that eigenvectors associated to eigenvalues with large absolute value localize around local maxima of $f_{\alpha}$: the metastable states in the power iteration method (slowing down its convergence) can be used to predict localization. We present a fast randomized algorithm and discuss different examples: a random band matrix, discretizations of the local operator $-\Delta + V$ and the nonlocal operator $(-\Delta)^{3/4} + V$.


title: Holomorphic Scalar Portals and the OPE

abstract: Visible-sector SUSY-breaking effects are computed in terms of hidden-sector correlation functions for generic holomorphic scalar portals. The solutions, which are valid irrespective of the hidden-sector dynamics, are approximated with the help of the operator product expansion (OPE). Indeed, for theories with superconformal symmetry at high energy, the superconformal OPE formalism can be used to disentangle the high-energy dynamics, encoded in the OPE coefficients, from the low-energy dynamics of the SUSY-breaking vacuum expectation values. A systematic method is proposed to compute the OPE coefficients, using relations between correlation functions of superfields and correlation functions of their quasi-primary component fields. The method, which is quite general, could be useful in building models of gauge- or gravity-mediated SUSY breaking and in analysing the viability of such models in a systematic way.


title: On the harmonic scattering amplitudes from a nonlinear elastic spherical inclusion

abstract: In this paper, the interaction of an incident finite amplitude longitudinal wave with a localized region of nonlinearity is considered. This interaction produces a secondary field represented by a superposition of first-, second-, and third-harmonic components. The secondary field is solely a result of the quadratic and cubic elastic nonlinearity present within the region of the inclusion. The second-harmonic scattering amplitude depends on the quadratic nonlinearity parameter $\beta$, while the first- and third-harmonic amplitudes depend on the cubic nonlinearity parameter $\gamma$. The special cases of forward and backward scattering amplitudes were analyzed. For each harmonic, the forward-scattering amplitude is always greater than or equal to the backward scattering amplitude in which the equality is only realized in the Rayleigh scattering limit. Lastly, the amplitudes of the scattered harmonic waves are compared to predicted harmonic amplitudes derived from a plane wave model.


title: The life cycle of starbursting circumnuclear gas discs

abstract: High-resolution observations from the sub-mm to the optical wavelength regime resolve the central few 100pc region of nearby galaxies in great detail. They reveal a large diversity of features: thick gas and stellar discs, nuclear starbursts, in- and outflows, central activity, jet interaction, etc. Concentrating on the role circumnuclear discs play in the life cycles of galactic nuclei, we employ 3D adaptive mesh refinement hydrodynamical simulations with the RAMSES code to self-consistently trace the evolution from a quasi-stable gas disc, undergoing gravitational (Toomre) instability, the formation of clumps and stars and the disc's subsequent, partial dispersal via stellar feedback. Our approach builds upon the observational finding that many nearby Seyfert galaxies have undergone intense nuclear starbursts in their recent past and in many nearby sources star formation is concentrated in a handful of clumps on a few 100pc distant from the galactic centre. We show that such observations can be understood as the result of gravitational instabilities in dense circumnuclear discs. By comparing these simulations to available integral field unit observations of a sample of nearby galactic nuclei, we find consistent gas and stellar masses, kinematics, star formation and outflow properties. Important ingredients in the simulations are the self-consistent treatment of star formation and the dynamical evolution of the stellar distribution as well as the modelling of a delay time distribution for the supernova feedback. The knowledge of the resulting simulated density structure and kinematics on pc scale is vital for understanding inflow and feedback processes towards galactic scales.


title: Spin ensemble-based AC magnetometry using concatenated dynamical decoupling at low temperatures

abstract: Ensembles of nitrogen-vacancy (NV) centers in diamond are widely used as AC magnetometers. While such measurements are usually performed using standard (XY) dynamical decoupling (DD) protocols at room temperature, we study the sensitivities achieved by utilizing various DD protocols, for measuring magnetic AC fields at frequencies in the 10-250 kHz range, at room temperature and 77 K. By performing measurements on an isotopically pure $^{12}$C sample, we find that the Carr-Purcell-Meiboom-Gill (CPMG) protocol, which is not robust against pulse imperfections, is less efficient for magnetometry than robust XY-based sequences. The concatenation of a standard XY-based protocol may enhance the sensitivities only for measuring high-frequency fields, for which many ($> 500$) DD pulses are necessary and the robustness against pulse imperfections is critical. Moreover, we show that cooling is effective only for measuring low-frequency fields (~10 kHz), for which the experiment time apporaches $T_1$ at a small number of applied DD pulses.


title: Vortex dynamics in type II superconductors under strong pinning conditions

abstract: We study effects of pinning on the dynamics of a vortex lattice in a type II superconductor in the strong-pinning situation and determine the force--velocity (or current--voltage) characteristic combining analytical and numerical methods. Our analysis deals with a small density $n_p$ of defects that act with a large force $f_p$ on the vortices, thereby inducing bistable configurations that are a characteristic feature of strong pinning theory. We determine the velocity-dependent average pinning-force density $\langle F_p(v)\rangle$ and find that it changes on the velocity scale $v_p \sim f_p/\eta a_0^3$, where $\eta$ is the viscosity of vortex motion and $a_0$ the distance between vortices. In the small pin-density limit, this velocity is much larger than the typical flow velocity $v_c \sim F_c/\eta$ of the free vortex system at drives near the critical force-density $F_c = \langle F_p(v=0)\rangle \propto n_p f_p$. As a result, we find a generic excess-force characteristic, a nearly linear force--velocity characteristic shifted by the critical force-density $F_c$; the linear flux-flow regime is approached only at large drives. Our analysis provides a derivation of Coulomb's law of dry friction for the case of strong vortex pinning.


title: Identifying and decoupling many-body interactions in spin ensembles in diamond

abstract: We simulate the dynamics of varying density quasi-two-dimensional spin ensembles in solid-state systems, focusing on the nitrogen-vacancy centers in diamond. We consider the effects of various control sequences on the averaged dynamics of large ensembles of spins, under a realistic "spin-bath" environment. We reveal that spin locking is efficient for decoupling spins initialized along the driving axis, both from coherent dipolar interactions and from the external spin-bath environment, when the driving is two orders of magnitude stronger than the relevant coupling energies. Since the application of standard pulsed dynamical decoupling sequences leads to strong decoupling from the environment, while other specialized pulse sequences can decouple coherent dipolar interactions, such sequences can be used to identify the dominant interaction type. Moreover, a proper combination of pulsed decoupling sequences could lead to the suppression of both interaction types, allowing additional spin manipulations. Finally, we consider the effect of finite-width pulses on these control protocols and identify improved decoupling efficiency with increased pulse duration, resulting from the interplay of dephasing and coherent dynamics.


title: Boundary Regularity for the Free Boundary in the One-phase Problem

abstract: We consider the Bernoulli one-phase free boundary problem in a domain $\Omega$ and show that the free boundary $F$ is $C^{1,1/2}$ regular in a neighborhood of the fixed boundary $\partial \Omega$. We achieve this by relating the behavior of $F$ near $\partial \Omega$ to a Signorini-type obstacle problem.


title: Approximate formula for the macroscopic polarization including quantum fluctuations

abstract: The many-body Berry phase formula for the macroscopic polarization is approximated by a sum of natural orbital geometric phases with fractional occupation numbers accounting for the dominant correlation effects. This reduced formula accurately reproduces the exact polarization in the Rice-Mele-Hubbard model across the band insulator-Mott insulator transition. A similar formula based on a one-body reduced Berry curvature accurately predicts the interaction-induced quenching of Thouless topological charge pumping.


title: On Zariski's multiplicity problem at infinity

abstract: We address a metric version of Zariski's multiplicity conjecture at infinity that says that two complex algebraic affine sets which are bi-Lipschitz homeomorphic at infinity must have the same degree. More specifically, we prove that the degree is a bi-Lipschitz invariant at infinity when the bi-Lipschitz homeomorphism has Lipschitz constants close to 1. In particular, we have that a family of complex algebraic sets bi-Lipschitz equisingular at infinity has constant degree. Moreover, we prove that if two polynomials are weakly rugose equivalent at infinity, then they have the same degree. In particular, we obtain that if two polynomials are rugose equivalent at infinity or bi-Lipschitz contact equivalent at infinity or bi-Lipschitz right-left equivalent at infinity, then they have the same degree.


title: When to arrive at a queue with earliness, tardiness and waiting costs

abstract: We consider a queueing facility where customers decide when to arrive. All customers have the same desired arrival time (w.l.o.g.\ time zero). There is one server, and the service times are independent and exponentially distributed. The total number of customers that demand service is random, and follows the Poisson distribution. Each customer wishes to minimize the sum of three costs: earliness, tardiness and waiting. We assume that all three costs are linear with time and are defined as follows. Earliness is the time between arrival and time zero, if there is any. Tardiness is simply the time of entering service, if it is after time zero. Waiting time is the time from arrival until entering service. We focus on customers' rational behaviour, assuming that each customer wants to minimize his total cost, and in particular, we seek a symmetric Nash equilibrium strategy. We show that such a strategy is mixed, unless trivialities occur. We construct a set of equations that its solution provides the symmetric Nash equilibrium. The solution is a continuous distribution on the real line. We also compare the socially optimal solution (that is, the one that minimizes total cost across all customers) to the overall cost resulting from the Nash equilibrium.


title: Excitation energy shift and size difference of low-energy levels in $p$-shell $\Lambda$ hypernuclei

abstract: We investigated structures of low-lying $0s$-orbit $\Lambda$ states in $p$-shell $\Lambda$ hypernuclei ($^A_\Lambda Z$) by applying microscopic cluster models for nuclear structure and a single-channel folding potential model for a $\Lambda$ particle. For $A>10$ systems, the size reduction of core nuclei is small, and the core polarization effect is regarded as a higher-order perturbation in the $\Lambda$ binding. The present calculation qualitatively describes the systematic trend of experimental data for excitation energy change from $^{A-1} Z$ to $^A_\Lambda Z$, in $A>10$ systems. The energy change shows a clear correlation with the nuclear size difference between the ground and excited states. In $^7_\Lambda \textrm{Li}$ and $^9_\Lambda \textrm{Be}$, the significant shrinkage of cluster structures occurs consistently with the prediction of other calculations.


title: Stack-Captioning: Coarse-to-Fine Learning for Image Captioning

abstract: The existing image captioning approaches typically train a one-stage sentence decoder, which is difficult to generate rich fine-grained descriptions. On the other hand, multi-stage image caption model is hard to train due to the vanishing gradient problem. In this paper, we propose a coarse-to-fine multi-stage prediction framework for image captioning, composed of multiple decoders each of which operates on the output of the previous stage, producing increasingly refined image descriptions. Our proposed learning approach addresses the difficulty of vanishing gradients during training by providing a learning objective function that enforces intermediate supervisions. Particularly, we optimize our model with a reinforcement learning approach which utilizes the output of each intermediate decoder's test-time inference algorithm as well as the output of its preceding decoder to normalize the rewards, which simultaneously solves the well-known exposure bias problem and the loss-evaluation mismatch problem. We extensively evaluate the proposed approach on MSCOCO and show that our approach can achieve the state-of-the-art performance.


title: Notes on Harmonic Analysis Part I: The Fourier Transform

abstract: Fourier Transforms is a first in a series of monographs we present on harmonic analysis. Harmonic analysis is one of the most fascinating areas of research in mathematics. Its centrality in the development of many areas of mathematics such as partial differential equations and integration theory and its many and diverse applications in sciences and engineering fields makes it an attractive field of study and research. The purpose of these notes is to introduce the basic ideas and theorems of the subject to students of mathematics, physics or engineering sciences. Our goal is to illustrate the topics with utmost clarity and accuracy, readily understandable by the students or interested readers. Rather than providing just the outlines or sketches of the proofs, we have actually provided the complete proofs of all theorems. This will illuminate the necessary steps taken and the machinery used to complete each proof. The prerequisite for understanding the topics presented is the knowledge of Lebesgue measure and integral. This will provide ample mathematical background for an advanced undergraduate or a graduate student in mathematics.


title: Analytic evaluation of Coulomb integrals for one, two and three-electron distance operators, $R_{C1}^{-n}R_{D1}^{-m}$, $R_{C1}^{-n}r_{12}^{-m}$ and $r_{12}^{-n}r_{13}^{-m}$ with $n, m=0,1,2$

abstract: The state of the art for integral evaluation is that analytical solutions to integrals are far more useful than numerical solutions. We evaluate certain integrals analytically that are necessary in some approaches in quantum chemistry. In the title, where R stands for nucleus-electron and r for electron-electron distances, the $(n,m)=(0,0)$ case is trivial, the $(n,m)=(1,0)$ and (0,1) cases are well known, fundamental milestone in integration and widely used in computation chemistry, as well as based on Laplace transformation with integrand exp(-$a^2t^2$). The rest of the cases are new and need the other Laplace transformation with integrand exp(-$a^2t$) also, as well as the necessity of a two dimensional version of Boys function comes up in case. These analytic expressions (up to Gaussian function integrand) are useful for manipulation with higher moments of inter-electronic distances, for example in correlation calculations.


title: Automated Generation of Anomalous Cancellations

abstract: If you cancel out the digit 6 from the ratio 16/64, you get the right answer by the wrong method. In 1979 R.P. Boas made an extensive study of such "anomalous cancellations" and generated many examples, in many bases. We continue his pioneering work, and generate many more examples, using the Extended Euclidean Algorithm.


title: Construction of divisible formal weight enumerators and extremal polynomials not satisfying the Riemann hypothesis

abstract: The formal weight enumerators were first introduced by M. Ozeki. They form a ring of invariant polynomials which is similar to that of the weight enumerators of Type II codes. Later, the zeta functions for linear codes were discovered and their theory was developed by I. Duursma. It was generalized to certain invariant polynomials including Ozeki's formal weight enumerators by the present author. One of the famous and important problems is whether extremal weight enumerators satisfy the Riemann hypothesis. In this paper, first we formulate the notion of divisible formal weight enumerators and propose an algorithm for the efficient search of the formal weight enumerators divisible by two. The main tools are the binomial moments. It leads to the discovery of several new families of formal weight enumerators. Then, as a result, we find examples of extremal formal weight enumerators which do not satisfy the Riemann hypothesis.


title: Self-propulsion against a moving membrane: enhanced accumulation and drag force

abstract: Self-propulsion (SP) is a main feature of active particles (AP), such as bacteria or biological micromotors, distinguishing them from passive colloids. A renowned consequence of SP is accumulation at static interfaces, even in the absence of hydrodynamic interactions. Here we address the role of SP in the interaction between AP and a moving semipermeable membrane. In particular, we implement a model of noninteracting AP in a channel crossed by a partially penetrable wall, moving at a constant velocity $c$. With respect to both the cases of passive colloids with $c>0$ and AP with $c=0$, the AP with finite $c$ show enhancement of accumulation in front of the obstacle and experience a largely increased drag force. This effect is understood in terms of an effective potential localised at the interface between particles and membrane, of height proportional to $c\tau/\xi$, where $\tau$ is the AP's re-orientation time and $\xi$ the width characterising the surface's smoothness ($\xi\to 0$ for hard core obstacles). An approximate analytical scheme is able to reproduce the observed density profiles and the measured drag force, in very good agreement with numerical simulations. The effects discussed here can be exploited for automatic selection and filtering of AP with desired parameters.


title: A Short and Elementary Proof of the Two-sidedness of the Matrix-Inverse

abstract: An elementary proof of the two-sidedness of the matrix-inverse is given using only linear independence and the reduced row-echelon form of a matrix. In addition, it is shown that a matrix is invertible if and only if it is row-equivalent to the identity matrix without appealing to elementary matrices. This proof underscores the importance of a basis and provides a proof of the invertible matrix theorem.


title: Derived Picard groups of preprojective algebras of Dynkin type

abstract: In this paper, we study two-sided tilting complexes of preprojective algebras of Dynkin type. We construct the most fundamental class of two-sided tilting complexes, which has a group structure by derived tensor products and induces a group of auto-equivalences of the derived category. We show that the group structure of the two-sided tilting complexes is isomorphic to the braid group of the corresponding folded graph. Moreover we show that these two-sided tilting complexes induce tilting mutation and any tilting complex is given as the derived tensor products of them. Using these results, we determine the derived Picard group of preprojective algebras for type $A$ and $D$.


title: Ghost Penalties in Nonconvex Constrained Optimization: Diminishing Stepsizes and Iteration Complexity

abstract: We consider nonconvex constrained optimization problems and propose a new approach to the convergence analysis based on penalty functions. We make use of classical penalty functions in an unconventional way, in that penalty functions only enter in the theoretical analysis of convergence while the algorithm itself is penalty-free. Based on this idea, we are able to establish several new results, including the first general analysis for diminishing stepsize methods in nonconvex, constrained optimization, showing convergence to generalized stationary points, and a complexity study for SQP-type algorithms.


title: The Recursive Stopping Time Structure of the $3x+1$ Function

abstract: The $3x + 1$ problem concerns iteration of the map $T:\mathbb{Z}\rightarrow\mathbb{Z}$ given by \begin{align*} T(x)=\left\{\begin{array}{lcr}\;\;\;\;\displaystyle{\frac{x}{2}} & \mbox{if $x\equiv 0\ (\text{mod}\ 2)$},\\ \\\displaystyle{\frac{3x+1}{2}} & \mbox{if $x\equiv 1\ (\text{mod}\ 2)$}.\end{array}\right. \end{align*} The $3x+1$ Conjecture states that every $x\geq1$ has some iterate $T^s(x)=1$. The least $s\in\mathbb{N}$ such that $T^s(x)<x$ is called the stopping time of $x$. It is shown that the congruence classes $(\text{mod}\ 2^k)$ of the integers having finite stopping time are given by a recursive algorithm producing a directed rooted tree.


title: Multiple scattering limit in optical microscopy

abstract: Optical microscopy offers a unique insight of biological structures with a sub-micrometer resolution and a minimum invasiveness. However, the inhomogeneities of the specimen itself can induce multiple scattering of light and optical aberrations which limit the observation to depths close to the surface. To predict quantitatively the penetration depth in microscopy, we theoretically derive the single-to-multiple scattering ratio in reflection. From this key quantity, the multiple scattering limit is deduced for various microscopic imaging techniques such as confocal microscopy, optical coherence tomography and related methods.


title: With the Brexit referendum Britain has spoken or has it

abstract: The Brexit referendum on the United Kingdom membership of the European Union took place on 23 June 2016. On the basis of a 52-48 split, statements such as: the majority of the UK chose to leave the EU, or the British people have voted to leave the European Union, or the will of the British people is to leave the EU, have pervaded political discourse and newspaper articles. In a recent paper, the authors present scientific arguments to contend that all those statements are without basis. They point out the fundamental issues with the Brexit referendum, which lie in how the data was summarized, ignoring fundamental statistical principles, and in the lack of terms of reference fixing a threshold for the result and for the minimum turnout. If carefully planned, the actual division of the United Kingdom would have been properly captured and the result---either way, for Leave or Remain---would have been readily accepted. In the present paper, the authors summarize the discussion and add some context in non-technical terms.


title: Statistical Physics of Design

abstract: A key challenge in complex design problems that permeate science and engineering is the need to balance design objectives for specific design elements or subsystems with global system objectives. Global objectives give rise to competing design pressures, whose effects can be difficult to trace in subsystem design. Here, using examples from arrangement problems, we show that the systems-level application of statistical physics principles, which we term "systems physics", provides a detailed characterization of subsystem design in terms of the concepts of stress and strain from materials physics. We analyze instances of routing problems in naval architectures, and show that systems physics provides a direct means of classifying architecture types, and quantifying trade-offs between subsystem- and overall performance. Our approach generalizes straightforwardly to design problems in a wide range of other disciplines that require concrete understanding of how the pressure to meet overall design objectives drives the outcomes for component subsystems.


title: Extremal invariant polynomials not satisfying the Riemann hypothesis

abstract: Zeta functions for linear codes were defined by Iwan Duursma in 1999. They were generalized to the case of some invariant polynomials by the preset author. One of the most important problems is whether extremal weight enumerators satisfy the Riemann hypothesis. In this article, we show there exist extremal polynomials of the weight enumerator type which are invariant under the MacWilliams transform and do not satisfy the Riemann hypothesis.


title: Deploy-As-You-Go Wireless Relay Placement: An Optimal Sequential Decision Approach using the Multi-Relay Channel Model

abstract: We use information theoretic achievable rate formulas for the multi-relay channel to study the problem of as-you-go deployment of relay nodes. The achievable rate formulas are for full-duplex radios at the relays and for decode-and-forward relaying. Deployment is done along the straight line joining a source node and a sink node at an unknown distance from the source. The problem is for a deployment agent to walk from the source to the sink, deploying relays as he walks, given that the distance to the sink is exponentially distributed with known mean. As a precursor, we apply the multi-relay channel achievable rate formula to obtain the optimal power allocation to relays placed along a line, at fixed locations. This permits us to obtain the optimal placement of a given number of nodes when the distance between the source and sink is given. Numerical work suggests that, at low attenuation, the relays are mostly clustered near the source in order to be able to cooperate, whereas at high attenuation they are uniformly placed and work as repeaters. We also prove that the effect of path-loss can be entirely mitigated if a large enough number of relays are placed uniformly between the source and the sink. The structure of the optimal power allocation for a given placement of the nodes, then motivates us to formulate the problem of as-you-go placement of relays along a line of exponentially distributed length, and with the exponential path-loss model, so as to minimize a cost function that is additive over hops. The hop cost trades off a capacity limiting term, motivated from the optimal power allocation solution, against the cost of adding a relay node. We formulate the problem as a total cost Markov decision process, establish results for the value function, and provide insights into the placement policy and the performance of the deployed network via numerical exploration.


title: The driven oscillator, with friction

abstract: This paper develops further the semi-classical theory of an harmonic oscillator acted on by a Gaussian white noise force discussed in (arXiv:1508.02379). Here I add to that theory the effects of Brownian damping (friction). Albeit semi-classical, the theory can be used to model quantum expectations and probabilities. I consider several examples.


title: A game theoretic model of wealth distribution

abstract: In this work we consider an agent based model in order to study the wealth distribution problem where the interchange is determined with a symmetric zero sum game. Simultaneously, the agents update their way of play trying to learn the optimal one. Here, the agents use mixed strategies. We study this model using both simulations and theoretical tools. We derive the equations for the learning mechanism, and we show that the mean strategy of the population satisfies an equation close to the classical replicator equation. Concerning the wealth distribution, there are two interesting situations depending on the equilibrium of the game. If the equilibrium is a pure strategy, the wealth distribution is fixed after some transient time, and those players which are close to optimal strategy are richer. When the game has an equilibrium in mixed strategies, the stationary wealth distribution is close to a Gamma distribution with variance depending on the coefficients of the game matrix. We compute theoretically their second moment in this case.


title: Optimal prediction in the linearly transformed spiked model

abstract: We consider the linearly transformed spiked model, where observations $Y_i$ are noisy linear transforms of unobserved signals of interest $X_i$: \begin{align*} Y_i = A_i X_i + \varepsilon_i, \end{align*} for $i=1,\ldots,n$. The transform matrices $A_i$ are also observed. We model $X_i$ as random vectors lying on an unknown low-dimensional space. How should we predict the unobserved signals (regression coefficients) $X_i$? The naive approach of performing regression for each observation separately is inaccurate due to the large noise. Instead, we develop optimal linear empirical Bayes methods for predicting $X_i$ by "borrowing strength" across the different samples. Our methods are applicable to large datasets and rely on weak moment assumptions. The analysis is based on random matrix theory. We discuss applications to signal processing, deconvolution, cryo-electron microscopy, and missing data in the high-noise regime. For missing data, we show in simulations that our methods are faster, more robust to noise and to unequal sampling than well-known matrix completion methods.


title: Massive fermions with low mobility in antiferromagnet orthorhombic CuMnAs single crystals

abstract: We report the physical properties of orthorhombic o-CuMnAs single crystal, which is predicted to be a topological Dirac semimetal with magnetic ground state and inversion symmetry broken. o-CuMnAs exhibits an antiferromagnetic transition with TN ~ 312 K. Further characterizations of magnetic properties suggest that the AFM order may be canted with the spin orientation in the bc plane. Small isotropic MR and linearly field-dependent Hall resistivity with positive slope indicate that single hole-type carries with high density and low mobility dominate the transport properties of o-CuMnAs. Furthermore, the result of low-temperature heat capacity shows that the effective mass of carriers is much larger than those in typical topological semimetals. These results imply that the carriers in o-CuMnAs exhibit remarkably different features from those of Dirac fermions predicted in theory.


title: Low-memory GEMM-based convolution algorithms for deep neural networks

abstract: Deep neural networks (DNNs) require very large amounts of computation both for training and for inference when deployed in the field. A common approach to implementing DNNs is to recast the most computationally expensive operations as general matrix multiplication (GEMM). However, as we demonstrate in this paper, there are a great many different ways to express DNN convolution operations using GEMM. Although different approaches all perform the same number of operations, the size of temporary data structures differs significantly. Convolution of an input matrix with dimensions $C \times H \times W$, requires $O(K^2CHW)$ additional space using the classical im2col approach. More recently memory-efficient approaches requiring just $O(KCHW)$ auxiliary space have been proposed. We present two novel GEMM-based algorithms that require just $O(MHW)$ and $O(KW)$ additional space respectively, where $M$ is the number of channels in the result of the convolution. These algorithms dramatically reduce the space overhead of DNN convolution, making it much more suitable for memory-limited embedded systems. Experimental evaluation shows that our low-memory algorithms are just as fast as the best patch-building approaches despite requiring just a fraction of the amount of additional memory. Our low-memory algorithms have excellent data locality which gives them a further edge over patch-building algorithms when multiple cores are used. As a result, our low memory algorithms often outperform the best patch-building algorithms using multiple threads.


title: On some families of divisible formal weight enumerators and their zeta functions

abstract: The formal weight enumerators were first introduced by M. Ozeki, and it was shown in the author's previous paper that there are various families of divisible formal weight enumerators. Among them, three families are dealt with in this paper and their properties are investigated: they are analogs of the Mallows-Sloane bound, the extremal property, the Riemann hypothesis, etc. In the course of the investigation, some generalizations of the theory of invariant differential operators developed by I. Duursma and T. Okuda are deduced.


title: Lepton number violation phenomenology of d=7 neutrino mass models

abstract: We study the phenomenology of d=7 1-loop neutrino mass models. All models in this particular class require the existence of several new $SU(2)_L$ multiplets, both scalar and fermionic, and thus predict a rich phenomenology at the LHC. The observed neutrino masses and mixings can easily be fitted in these models. Interestingly, despite the smallness of the observed neutrino masses, some particular lepton number violating (LNV) final states can arise with observable branching ratios. These LNV final states consists of leptons and gauge bosons with high multiplicities, such as 4l+4W, 6l+2W, etc. We study current constraints on these models from upper bounds on charged lepton flavour violating decays, existing lepton number conserving searches at the LHC and discuss possible future LNV searches.


title: More Infinite Products: Thue-Morse and the Gamma function

abstract: Letting $(t_n)$ denote the Thue-Morse sequence with values $0, 1$, we note that the Woods-Robbins product $$ \prod_{n \geq 0} \left(\frac{2n+1}{2n+2}\right)^{(-1)^{t_n}} = 2^{-1/2} $$ involves a rational function in $n$ and the $\pm 1$ Thue-Morse sequence $((-1)^{t_n})_{n \geq 0}$. The purpose of this paper is twofold. On the one hand, we try to find other rational functions for which similar infinite products involving the $\pm 1$ Thue-Morse sequence have an expression in terms of known constants. On the other hand, we also try to find (possibly different) rational functions $R$ for which the infinite product $\prod R(n)^{t_n}$ also has an expression in terms of known constants.


title: Automated Identification of Trampoline Skills Using Computer Vision Extracted Pose Estimation

abstract: A novel method to identify trampoline skills using a single video camera is proposed herein. Conventional computer vision techniques are used for identification, estimation, and tracking of the gymnast's body in a video recording of the routine. For each frame, an open source convolutional neural network is used to estimate the pose of the athlete's body. Body orientation and joint angle estimates are extracted from these pose estimates. The trajectories of these angle estimates over time are compared with those of labelled reference skills. A nearest neighbour classifier utilising a mean squared error distance metric is used to identify the skill performed. A dataset containing 714 skill examples with 20 distinct skills performed by adult male and female gymnasts was recorded and used for evaluation of the system. The system was found to achieve a skill identification accuracy of 80.7% for the dataset.


title: Fragmentation functions of charged hadrons

abstract: I present a preliminary determination of the Fragmentation Functions (FFs) of unidentified charged hadrons at next-to-leading order in quantum chromodynamics. The analysis is based on hadron production cross section data in single-inclusive electron-positron annihilation. It extends a recent determination of the FFs of identified charged pions, charged kaons and protons/antiprotons performed by the NNPDF Collaboration. I illustrate the quality of the FFs determined in this analysis and show how they describe the charged hadron spectra measured in proton-(anti)proton collisions.


title: EndoSensorFusion: Particle Filtering-Based Multi-sensory Data Fusion with Switching State-Space Model for Endoscopic Capsule Robots

abstract: A reliable, real time multi-sensor fusion functionality is crucial for localization of actively controlled capsule endoscopy robots, which are an emerging, minimally invasive diagnostic and therapeutic technology for the gastrointestinal (GI) tract. In this study, we propose a novel multi-sensor fusion approach based on a particle filter that incorporates an online estimation of sensor reliability and a non-linear kinematic model learned by a recurrent neural network. Our method sequentially estimates the true robot pose from noisy pose observations delivered by multiple sensors. We experimentally test the method using 5 degree-of-freedom (5-DoF) absolute pose measurement by a magnetic localization system and a 6-DoF relative pose measurement by visual odometry. In addition, the proposed method is capable of detecting and handling sensor failures by ignoring corrupted data, providing the robustness expected of a medical device. Detailed analyses and evaluations are presented using ex-vivo experiments on a porcine stomach model prove that our system achieves high translational and rotational accuracies for different types of endoscopic capsule robot trajectories.


title: Ultrafast carbon monoxide photolysis and heme spin-crossover in myoglobin via nonadiabatic quantum dynamics

abstract: Light absorption of myoglobin triggers diatomic ligand photolysis and spin crossover transition of iron(II) that initiate protein conformational change. The photolysis and spin crossover reactions happen concurrently on a femtosecond timescale. The microscopic origin of these reactions remains controversial. Here, we apply quantum wavepacket dynamics to elucidate the ultrafast photochemical mechanism for a heme--carbon monoxide (heme--CO) complex. We observe coherent oscillations of the Fe-CO bond distance with a period of 42 fs and an amplitude of $\sim$1 \AA{}. These nuclear motions induce pronounced geometric reorganization, which makes the CO dissociation irreversible. The reaction is initially dominated by symmetry breaking vibrations inducing an electron transfer from porphyrin to iron. Subsequently, the wavepacket relaxes to the triplet manifold in $\sim$75 fs and to the quintet manifold in $\sim$430 fs. Our results highlight the central role of nuclear vibrations at the origin of the ultrafast photodynamics of organometallic complexes.


title: On some families of invariant polynomials divisible by three and their zeta functions

abstract: In this note, we establish an analog of the Mallows-Sloane bound for Type III formal weight enumerators. This completes the bounds for all types (Types I through IV) in synthesis of our previous results. Next we show by using the binomial moments that there exists a family of polynomials divisible by three, which are not related to linear codes but are invariant under the MacWilliams transform for the value 3/2. We also discuss some properties of the zeta functions for such polynomials.


title: A Domain-specific Language for High-reliability Software used in the JUICE SWI Instrument - The hO Language Manual

abstract: hO is a custom restricted dialect of Oberon, developed at the Max-Planck Institute for Solar System Research in G\"ottingen and used in the SWI flight software for the JUICE mission. hO is applied to reduce the possibility of syntactically valid but incorrect code, provide better means of statically analyzing source code, is more readable than C and gives syntactic support for the software architecture used in the SWI instrument software. By using a higher-level, application-specific notation a whole range of possible errors is eliminated and source code size is reduced, while making the code itself easier to understand, review and analyze.


title: The Castalia Mission to Main Belt Comet 133P/Elst-Pizarro

abstract: We describe Castalia, a proposed mission to rendezvous with a Main Belt Comet (MBC), 133P/Elst-Pizarro. MBCs are a recently discovered population of apparently icy bodies within the main asteroid belt between Mars and Jupiter, which may represent the remnants of the population which supplied the early Earth with water. Castalia will perform the first exploration of this population by characterising 133P in detail, solving the puzzle of the MBC's activity, and making the first in situ measurements of water in the asteroid belt. In many ways a successor to ESA's highly successful Rosetta mission, Castalia will allow direct comparison between very different classes of comet, including measuring critical isotope ratios, plasma and dust properties. It will also feature the first radar system to visit a minor body, mapping the ice in the interior. Castalia was proposed, in slightly different versions, to the ESA M4 and M5 calls within the Cosmic Vision programme. We describe the science motivation for the mission, the measurements required to achieve the scientific goals, and the proposed instrument payload and spacecraft to achieve these.


title: Social Media Text Processing and Semantic Analysis for Smart Cities

abstract: With the rise of Social Media, people obtain and share information almost instantly on a 24/7 basis. Many research areas have tried to gain valuable insights from these large volumes of freely available user generated content. With the goal of extracting knowledge from social media streams that might be useful in the context of intelligent transportation systems and smart cities, we designed and developed a framework that provides functionalities for parallel collection of geo-located tweets from multiple pre-defined bounding boxes (cities or regions), including filtering of non-complying tweets, text pre-processing for Portuguese and English language, topic modeling, and transportation-specific text classifiers, as well as, aggregation and data visualization. We performed an exploratory data analysis of geo-located tweets in 5 different cities: Rio de Janeiro, S\~ao Paulo, New York City, London and Melbourne, comprising a total of more than 43 million tweets in a period of 3 months. Furthermore, we performed a large scale topic modelling comparison between Rio de Janeiro and S\~ao Paulo. Interestingly, most of the topics are shared between both cities which despite being in the same country are considered very different regarding population, economy and lifestyle. We take advantage of recent developments in word embeddings and train such representations from the collections of geo-located tweets. We then use a combination of bag-of-embeddings and traditional bag-of-words to train travel-related classifiers in both Portuguese and English to filter travel-related content from non-related. We created specific gold-standard data to perform empirical evaluation of the resulting classifiers. Results are in line with research work in other application areas by showing the robustness of using word embeddings to learn word similarities that bag-of-words is not able to capture.


title: Asymptotic normality of Laplacian coefficients of graphs

abstract: Let $G$ be a simple graph with $n$ vertices and let $$C(G;x)=\sum_{k=0}^n(-1)^{n-k}c(G,k)x^k$$ denote the Laplacian characteristic polynomial of $G$. Then if the size $|E(G)|$ is large compared to the maximum degree $\Delta(G)$, Laplacian coefficients $c(G,k)$ are approximately normally distributed (by central and local limit theorems). We show that Laplacian coefficients of the paths, the cycles, the stars, the wheels and regular graphs of degree $d$ are approximately normally distributed respectively. We also point out that Laplacian coefficients of the complete graphs and the complete bipartite graphs are approximately Poisson distributed respectively.


title: Plant survival and keystone pollinator species in stochastic coextinction models: role of intrinsic dependence on animal-pollination

abstract: Coextinction models are useful to understand community robustness to species loss and resilience to disturbances. We simulated pollinator extinctions in pollination networks by using a hybrid model that combined a recently developed stochastic coextinction model (SCM) for plant extinctions and a topological model (TCM) for animal extinctions. Our model accounted for variation in interaction strengths and included empirical estimates of plant dependence on pollinators to set seeds. The stochastic nature of such model allowed us determining plant survival to single (and multiple) extinction events, and identifying which pollinators (keystone species) were more likely to trigger secondary extinctions. Consistently across three different pollinator removal sequences, plant robustness was lower than in a pure TCM, and plant survival was more determined by dependence on the mutualism than by interaction strength. As expected, highly connected and dependent plants were the most sensitive to pollinator loss and collapsed faster in extinction cascades. We predict that the relationship between dependence and plant connectivity is crucial to determine network robustness to interaction loss. Finally, we showed that honeybees and several beetles were keystone species in our communities. This information is of great value to foresee consequences of pollinator losses facing current global change and to identify target species for effective conservation.


title: Deep Shape Matching

abstract: We cast shape matching as metric learning with convolutional networks. We break the end-to-end process of image representation into two parts. Firstly, well established efficient methods are chosen to turn the images into edge maps. Secondly, the network is trained with edge maps of landmark images, which are automatically obtained by a structure-from-motion pipeline. The learned representation is evaluated on a range of different tasks, providing improvements on challenging cases of domain generalization, generic sketch-based image retrieval or its fine-grained counterpart. In contrast to other methods that learn a different model per task, object category, or domain, we use the same network throughout all our experiments, achieving state-of-the-art results in multiple benchmarks.


title: One-Shot Learning for Semantic Segmentation

abstract: Low-shot learning methods for image classification support learning from sparse data. We extend these techniques to support dense semantic image segmentation. Specifically, we train a network that, given a small set of annotated images, produces parameters for a Fully Convolutional Network (FCN). We use this FCN to perform dense pixel-level prediction on a test image for the new semantic class. Our architecture shows a 25% relative meanIoU improvement compared to the best baseline methods for one-shot segmentation on unseen classes in the PASCAL VOC 2012 dataset and is at least 3 times faster.


title: Acute sets of exponentially optimal size

abstract: We present a simple construction of an acute set of size $2^{d-1}+1$ in $\mathbb{R}^d$ for any dimension $d$. That is, we explicitly give $2^{d-1}+1$ points in the $d$-dimensional Euclidean space with the property that any three points form an acute triangle. It is known that the maximal number of such points is less than $2^d$. Our result significantly improves upon a recent construction, due to Dmitriy Zakharov, with size of order $\varphi^d$ where $\varphi = (1+\sqrt{5})/2 \approx 1.618$ is the golden ratio.


title: Generalized Poisson integral and sharp estimates for harmonic and biharmonic functions in the half-space

abstract: A representation for the sharp coefficient in a pointwise estimate for the gradient of a generalized Poisson integral of a function $f$ on ${\mathbb R}^{n-1}$ is obtained under the assumption that $f$ belongs to $L^p$. It is assumed that the kernel of the integral depends on the parameters $\alpha$ and $\beta$. The explicit formulas for the sharp coefficients are found for the cases $p=1$, $p=2$ and for some values of $\alpha , \beta$ in the case $p=\infty$. Conditions ensuring the validity of some analogues of the Khavinson's conjecture for the generalized Poisson integral are obtained. The sharp estimates are applied to harmonic and biharmonic functions in the half-space.


title: Gigamachine: incremental machine learning on desktop computers

abstract: We present a concrete design for Solomonoff's incremental machine learning system suitable for desktop computers. We use R5RS Scheme and its standard library with a few omissions as the reference machine. We introduce a Levin Search variant based on a stochastic Context Free Grammar together with new update algorithms that use the same grammar as a guiding probability distribution for incremental machine learning. The updates include adjusting production probabilities, re-using previous solutions, learning programming idioms and discovery of frequent subprograms. The issues of extending the a priori probability distribution and bootstrapping are discussed. We have implemented a good portion of the proposed algorithms. Experiments with toy problems show that the update algorithms work as expected.


title: X-Shooter observations of low-mass stars in the Eta Chamaeleontis association

abstract: The nearby Eta Chamaeleontis association is a collection of 4-10 Myr old stars with a disk fraction of 35-45%. In this study, the broad wavelength coverage of VLT/X-Shooter is used to measure the stellar and mass accretion properties of 15 low mass stars in the Eta Chamaeleontis association. For each star, the observed spectrum is fitted with a non-accreting stellar template and an accretion spectrum obtained from assuming a plane-parallel hydrogen slab. Five of the eight stars with an IR disk excess show excess UV emission, indicating ongoing accretion. The accretion rates measured here are similar to those obtained from previous measurements of excess UV emission, but tend to be higher than past measurements from H-alpha modeling. The mass accretion rates are consistent with those of other young star forming regions.


title: What Quark-Gluon Plasma in small systems might tell us about nucleons

abstract: The origin of flow-like effects in small systems, such as those produced in ultra-relativistic proton-proton and proton-lead collisions, is still widely debated. In this paper the goal is to look at possible consequences if indeed a mini-Quark-Gluon Plasma is formed in these collisions. It is argued that this could indicate a duality between the QGP phase and the color fields in hadrons. A qualitative dense field picture is presented for this duality and discussed.


title: High Impedance Detector Arrays for Magnetic Resonance

abstract: Resonant inductive coupling is commonly seen as an undesired fundamental phenomenon emergent in densely packed resonant structures, such as nuclear magnetic resonance phased array detectors. The need to mitigate coupling imposes rigid constraints on the detector design, impeding performance and limiting the scope of magnetic resonance experiments. Here we introduce a high impedance detector design, which can cloak itself from electrodynamic interactions with neighboring elements. We verify experimentally that the high impedance detectors do not suffer from signal-to-noise degradation mechanisms observed with traditional low impedance elements. Using this new-found robustness, we demonstrate an adaptive wearable detector array for magnetic resonance imaging of the hand. The unique properties of the detector glove reveal new pathways to study the biomechanics of soft tissues, and exemplify the enabling potential of high-impedance detectors for a wide range of demanding applications that are not well suited to traditional coil designs.


title: Determination of Trace Levels of Uranium and Thorium in High Purity Gadolinium Sulfate Using ICP-MS with Solid-Phase Chromatographic Extraction Resin

abstract: The new Super-Kamiokande-Gadolinium (SK-Gd) project is an upgrade of the Super-Kamiokande (SK) detector. In the SK-Gd project, 0.2\% Gd$_2$(SO$_4$)$_3$ is loaded into the 50 kton water tank of the SK. One of the main purposes of the project is to discover Supernova Relic Neutrinos. Neutrino measurements and proton decay searches will also be performed in the SK-Gd. In order to measure solar neutrinos with a low energy threshold of $\sim$3.5 MeV in the SK-Gd, the main radioactive contaminations, $^{238}$U and $^{232}$Th, in Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O, should be minimized before loading. Our maximum levels for U and Th are 5 mBq (U)/kg (Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O) and 0.05 mBq (Th)/kg (Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O). In order to measure such low concentrations of U and Th in Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O, we developed the solid-phase extraction technique. Using this method, about 90\% or more U and Th could be efficiently extracted while Gd was reduced by a factor of about $10^{4}$. This allowed these radioactivity contaminations to be measured precisely as 0.04 mBq/kg (Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O) for U and 0.01 mBq/kg (Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O) for Th. We measured three pure Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O samples using this method and estimated that the purest one contained $<0.04$ mBq (U)/kg (Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O) and 0.06 $\pm$ 0.01 mBq (Th)/kg (Gd$_2$(SO$_4$)$_3{\cdot}$8H$_2$O) by ICP-MS.


title: Simulation of Integro-Differential Equation and Application in Estimation of Ruin Probability with Mixed Fractional Brownian Motion

abstract: In this paper, we are concerned with the numerical solution of one type integro-differential equation by a probability method based on the fundamental martingale of mixed Gaussian processes. As an application, we will try to simulate the estimation of ruin probability with an unknown parameter driven not by the classical L\'evy process but by the mixed fractional Brownian motion.


title: Planetary Radio Interferometry and Doppler Experiment (PRIDE) Technique: a Test Case of the Mars Express Phobos Fly-by. 2. Doppler tracking: Formulation of observed and computed values, and noise budget

abstract: Context. Closed-loop Doppler data obtained by deep space tracking networks (e.g., NASA's DSN and ESA's Estrack) are routinely used for navigation and science applications. By "shadow tracking" the spacecraft signal, Earth-based radio telescopes involved in Planetary Radio Interferometry and Doppler Experiment (PRIDE) can provide open-loop Doppler tracking data when the dedicated deep space tracking facilities are operating in closed-loop mode only. Aims. We explain in detail the data processing pipeline, discuss the capabilities of the technique and its potential applications in planetary science. Methods. We provide the formulation of the observed and computed values of the Doppler data in PRIDE tracking of spacecraft, and demonstrate the quality of the results using as a test case an experiment with ESA's Mars Express spacecraft. Results. We find that the Doppler residuals and the corresponding noise budget of the open-loop Doppler detections obtained with the PRIDE stations are comparable to the closed-loop Doppler detections obtained with the dedicated deep space tracking facilities.


title: Clustering dark energy and halo abundances

abstract: Within the standard paradigm, dark energy is taken as a homogeneous fluid that drives the accelerated expansion of the universe and does not contribute to the mass of collapsed objects such as galaxies and galaxy clusters. The abundance of galaxy clusters -- measured through a variety of channels -- has been extensively used to constrain the normalization of the power spectrum: it is an important probe as it allows us to test if the standard $\Lambda$CDM model can indeed accurately describe the evolution of structures across billions of years. It is then quite significant that the Planck satellite has detected, via the Sunyaev-Zel'dovich effect, less clusters than expected according to the primary CMB anisotropies. One of the simplest generalizations that could reconcile these observations is to consider models in which dark energy is allowed to cluster, i.e., allowing its sound speed to vary. In this case, however, the standard methods to compute the abundance of galaxy clusters need to be adapted to account for the contributions of dark energy. In particular, we examine the case of clustering dark energy -- a dark energy fluid with negligible sound speed -- with a redshift-dependent equation of state. We carefully study how the halo mass function is modified in this scenario, highlighting corrections that have not been considered before in the literature. We address modifications in the growth function, collapse threshold, virialization densities and also changes in the comoving scale of collapse and mass function normalization. Our results show that clustering dark energy can impact halo abundances at the level of 10\%--30\%, depending on the halo mass, and that cluster counts are modified by about 30\% at a redshift of unity.


title: Modeling and identification of uncertain-input systems

abstract: In this work, we present a new class of models, called uncertain-input models, that allows us to treat system-identification problems in which a linear system is subject to a partially unknown input signal. To encode prior information about the input or the linear system, we use Gaussian-process models. We estimate the model from data using the empirical Bayes approach: the input and the impulse responses of the linear system are estimated using the posterior means of the Gaussian-process models given the data, and the hyperparameters that characterize the Gaussian-process models are estimated from the marginal likelihood of the data. We propose an iterative algorithm to find the hyperparameters that relies on the EM method and results in simple update steps. In the most general formulation, neither the marginal likelihood nor the posterior distribution of the unknowns is tractable. Therefore, we propose two approximation approaches, one based on Markov-chain Monte Carlo techniques and one based on variational Bayes approximation. We also show special model structures for which the distributions are treatable exactly. Through numerical simulations, we study the application of the uncertain-input model to the identification of Hammerstein systems and cascaded linear systems. As part of the contribution of the paper, we show that this model structure encompasses many classical problems in system identification such as classical PEM, Hammerstein models, errors-in-variables problems, blind system identification, and cascaded linear systems. This allows us to build a systematic procedure to apply the algorithms proposed in this work to a wide class of classical problems.


title: On the de-Rham cohomology of hyperelliptic curves

abstract: For any hyperelliptic curve X, we give an explicit basis of the first de-Rham cohomology of X in terms of \v{C}ech cohomology. We use this to produce a family of curves in characteristic p>2 for which the Hodge-de-Rham short exact sequence does not split equivariantly; this generalises a result of Hortsch. Further, we use our basis to show that the hyperelliptic involution acts on the first de-Rham cohomology by multiplication by -1, i.e., acts as the identity when p=2.


title: Ensemble Methods as a Defense to Adversarial Perturbations Against Deep Neural Networks

abstract: Deep learning has become the state of the art approach in many machine learning problems such as classification. It has recently been shown that deep learning is highly vulnerable to adversarial perturbations. Taking the camera systems of self-driving cars as an example, small adversarial perturbations can cause the system to make errors in important tasks, such as classifying traffic signs or detecting pedestrians. Hence, in order to use deep learning without safety concerns a proper defense strategy is required. We propose to use ensemble methods as a defense strategy against adversarial perturbations. We find that an attack leading one model to misclassify does not imply the same for other networks performing the same task. This makes ensemble methods an attractive defense strategy against adversarial attacks. We empirically show for the MNIST and the CIFAR-10 data sets that ensemble methods not only improve the accuracy of neural networks on test data but also increase their robustness against adversarial perturbations.


title: Constant-Weight Array Codes

abstract: Binary constant-weight codes have been extensively studied, due to both their numerous applications and to their theoretical significance. In particular, constant-weight codes have been proposed for error correction in store and forward. In this paper, we introduce constant-weight array codes (CWACs), which offer a tradeoff between the rate gain of general constant-weight codes and the low decoding complexity of liftings. CWACs can either be used in the on-shot setting introduced earlier or in a multi-shot approach, where one codeword consists of several messages. The multi-shot approach generalizes the one-shot approach and hence allows for higher rate gains. We first give a construction of CWACs based on concatenation, which generalizes the traditional erasure codes, and also provide a decoding algorithm for these codes. Since CWACs can be viewed as a generalization of both binary constant-weight codes and nonrestricted Hamming metric codes, CWACs thus provide an additional degree of freedom to both problems of determining the maximum cardinality of constant-weight codes and nonrestricted Hamming metric codes. We then investigate their theoretical significance. We first generalize many classical bounds derived for Hamming metric codes or constant-weight codes in the CWAC framework. We finally relate the maximum cardinality of a CWAC to that of a constant-weight code, of a nonrestricted Hamming metric code, and of a spherical code.


title: On the Hartogs extension theorem for unbounded domains in $\mathbb{C}^n$

abstract: Let $\Omega\subset\mathbb{C}^n$, $n\geq 2$, be a domain with smooth connected boundary. If $\Omega$ is relatively compact, the Hartogs-Bochner theorem ensures that every CR distribution on $\partial\Omega$ has a holomorphic extension to $\Omega$. For unbounded domains this extension property may fail, for example if $\Omega$ contains a complex hypersurface. The main result in this paper tells that the extension property holds if and only if the envelope of holomorphy of $\mathbb{C}^n\backslash\overline{\Omega}$ is $\mathbb{C}^n$. It seems that it is a first result in the literature which gives a geometric characterization of unbounded domains in $\mathbb C^n$ for which the Hartogs phenomenon holds. Comparing this to earlier work by the first two authors and Z.~S{\l}odkowski, one observes that the extension problem sensitively depends on a finer geometry of the contact of a complex hypersurface and the boundary of the domain.


title: Trajectory Synthesis for Fisher Information Maximization

abstract: Estimation of model parameters in a dynamic system can be significantly improved with the choice of experimental trajectory. For general, nonlinear dynamic systems, finding globally "best" trajectories is typically not feasible; however, given an initial estimate of the model parameters and an initial trajectory, we present a continuous-time optimization method that produces a locally optimal trajectory for parameter estimation in the presence of measurement noise. The optimization algorithm is formulated to find system trajectories that improve a norm on the Fisher information matrix. A double-pendulum cart apparatus is used to numerically and experimentally validate this technique. In simulation, the optimized trajectory increases the minimum eigenvalue of the Fisher information matrix by three orders of magnitude compared to the initial trajectory. Experimental results show that this optimized trajectory translates to an order of magnitude improvement in the parameter estimate error in practice.


title: Electronic and magnetic excitations in the "half-stuffed" Cu--O planes of Ba$_2$Cu$_3$O$_4$Cl$_2$ measured by resonant inelastic x-ray scattering

abstract: We use resonant inelastic x-ray scattering (RIXS) at the Cu L$_3$ edge to measure the charge and spin excitations in the "half-stuffed" Cu--O planes of the cuprate antiferromagnet Ba$_2$Cu$_3$O$_4$Cl$_2$. The RIXS line shape reveals distinct contributions to the $dd$ excitations from the two structurally inequivalent Cu sites, which have different out-of-plane coordinations. The low-energy response exhibits magnetic excitations. We find a spin-wave branch whose dispersion follows the symmetry of a CuO$_2$ sublattice, similar to the case of the "fully-stuffed" planes of tetragonal CuO (T-CuO). Its bandwidth is closer to that of a typical cuprate material, such as Sr$_2$CuO$_2$Cl$_2$, than it is to that of T-CuO. We interpret this result as arising from the absence of the effective four-spin inter-sublattice interactions that act to reduce the bandwidth in T-CuO.


title: Coherent metamaterial absorption of two-photon states with 40% efficiency

abstract: Multi-photon absorption processes have a nonlinear dependence on the amplitude of the incident optical field i.e. the number of photons. However, multi-photon absorption is generally weak and multi-photon events occur with extremely low probability. Consequently, it is extremely challenging to engineer quantum nonlinear devices that operate at the single photon level and the majority of quantum technologies have to rely on single photon interactions. Here, we demonstrate experimentally and theoretically that exploiting coherent absorption of N = 2 N00N states makes it possible to enhance the number of two-photon states that are absorbed. An absorbing metasurface placed inside a Sagnac-style interferometer into which we inject an N = 2 N00N state, exhibits two-photon absorption with 40.5% efficiency, close to the theoretical maximum. This high probability of simultaneous absorption of two photons holds the promise for applications in fields that require multi-photon upconversion but are hindered by high peak intensities.


title: String chopping and time-ordered products of linear string-localized quantum fields

abstract: For a renormalizability proof of perturbative models in the Epstein--Glaser scheme with string-localized quantum fields, one needs to know what freedom one has in the definition of time-ordered products of the interaction Lagrangian. This paper provides a first step in that direction. The basic issue is the presence of an open set of $n$-tuples of strings which cannot be chronologically ordered. We resolve it by showing that almost all such string configurations can be dissected into finitely many pieces which can indeed be chronologically ordered. This fixes the time-ordered products of linear field factors outside a nullset of string configurations. (The extension across the nullset, as well as the definition of time-ordered products of Wick monomials, will be discussed elsewhere.)


title: A Tannakian interpretation of the elliptic infinitesimal braid Lie algebras

abstract: Let $n\geq 1$. The pro-unipotent completion of the pure braid group of $n$ points on a genus 1 surface has been shown to be isomorphic to an explicit pro-unipotent group with graded Lie algebra using two types of tools: (a) minimal models (Bezrukavnikov), (b) the choice of a complex structure on the genus 1 surface, making it into an elliptic curve $E$, and an appropriate flat connection on the configuration space of $n$ points in $E$ (joint work of the authors with D. Calaque). Following a suggestion by P. Deligne, we give an interpretation of this isomorphism in the framework of the Riemann-Hilbert correspondence, using the total space $E^\#$ of an affine line bundle over $E$, which identifies with the moduli space of line bundles over $E$ equipped with a flat connection.


title: A Longitudinal Higher-Order Diagnostic Classification Model

abstract: Providing diagnostic feedback about growth is crucial to formative decisions such as targeted remedial instructions or interventions. This paper proposed a longitudinal higher-order diagnostic classification modeling approach for measuring growth. The new modeling approach is able to provide quantitative values of overall and individual growth by constructing a multidimensional higher-order latent structure to take into account the correlations among multiple latent attributes that are examined across different occasions. In addition, potential local item dependence among anchor (or repeated) items can also be taken into account. Model parameter estimation is explored in a simulation study. An empirical example is analyzed to illustrate the applications and advantages of the proposed modeling approach.


title: Stability of the Kasner Universe in f(T) Gravity

abstract: $f(T)$ gravity offers an alternative context in which to consider gravitational interactions where torsion, rather than curvature, is the mechanism by which gravitation is communicated. We investigate the stability of the Kasner solution with several forms of the arbitrary lagrangian function examined within the $f(T)$ context. This is a Bianchi type--I vacuum solution with anisotropic expansion factors. In the $f(T)$ gravity setting, the solution must conform to a set of conditions in order to continue to be a vacuum solution of the generalized field equations. With this solution in hand, the perturbed field equations are determined for power-law and exponential forms of the $f(T)$ function. We find that the point which describes the Kasner solution is a saddle point which means that the singular solution is unstable. However, we find the de Sitter universe is a late-time attractor. In general relativity, the cosmological constant drives the the isotropization of the spacetime while in this setting the extra $f(T)$ contributions now provide this impetus.


title: Asymptotic Geometry of the Hitchin Metric

abstract: We study the asymptotics of the natural $L^2$ metric on the Hitchin moduli space with group $G = \mathrm{SU}(2)$. Our main result, which addresses a detailed conjectural picture made by Gaiotto, Neitzke and Moore \cite{gmn13}, is that on the regular part of the Hitchin system, this metric is well-approximated by the semiflat metric from \cite{gmn13}. We prove that the asymptotic rate of convergence for gauged tangent vectors to the moduli space has a precise polynomial expansion, and hence that the the difference between the two sets of metric coefficients in a certain natural coordinate system also has polynomial decay. Very recent work by Dumas and Neitzke indicates that the convergence rate for the metric is exponential, at least in certain directions.


title: Neutrino Interferometry for High-Precision Tests of Lorentz Symmetry with IceCube

abstract: Lorentz symmetry is a fundamental space-time symmetry underlying the Standard Model of particle physics and gravity. However, unified theories, such as string theory, allow for violation of this symmetry. Thus, the discovery of Lorentz symmetry violation could be the first hint of these theories. Here, we use high-energy atmospheric neutrinos observed at the IceCube Neutrino Observatory to search for anomalous neutrino oscillations as signals of Lorentz violation. The large range of neutrino energies and propagation baselines, together with high statistics, let us perform the most precise test of space-time symmetry in the neutrino sector to date. We find no evidence for Lorentz violation. This allows us to constrain the size of the dimension-four operator in the Standard-Model Extension for Lorentz violation to the $10^{-28}$ level and to set limits on higher dimensional operators of that theory. These are among the most stringent limits on Lorentz violation across all fields of physics.


title: Positive polynomials on unbounded domains

abstract: Certificates of non-negativity such as Putinar's Positivstellensatz have been used to obtain powerful numerical techniques to solve polynomial optimization (PO) problems. Putinar's certificate uses sum-of-squares (sos) polynomials to certify the non-negativity of a given polynomial over a domain defined by polynomial inequalities. This certificate assumes the Archimedean property of the associated quadratic module, which in particular implies compactness of the domain. In this paper we characterize the existence of a certificate of non-negativity for polynomials over a possibly unbounded domain, without the use of the associated quadratic module. Next, we show that the certificate can be used to convergent linear matrix inequality (LMI) hierarchies for PO problems with unbounded feasible sets. Furthermore, by using copositive polynomials to certify non-negativity, instead of sos polynomials, the certificate allows the use of a very rich class of convergent LMI hierarchies to approximate the solution of general PO problems. Throughout the article we illustrate our results with various examples certifying the non-negativity of polynomials over possibly unbounded sets defined by polynomial equalities or inequalities.