title: Global well-posedness of the one-dimensional cubic nonlinear Schr\"odinger equation in almost critical spaces

abstract: In this paper, we first introduce a new function space $MH^{\theta, p}$ whose norm is given by the $\ell^p$-sum of modulated $H^\theta$-norms of a given function. In particular, when $\theta < -\frac 12$, we show that the space $MH^{\theta, p}$ agrees with the modulation space $M^{2, p}(\mathbb R)$ on the real line and the Fourier-Lebesgue space $\mathcal F L^{p}(\mathbb T)$ on the circle. We use this equivalence of the norms and the Galilean symmetry to adapt the conserved quantities constructed by Killip-Vi\c{s}an-Zhang to the modulation space setting. By applying the scaling symmetry, we then prove global well-posedness of the one-dimensional cubic nonlinear Schr\"odinger equation (NLS) in almost critical spaces. More precisely, we show that the cubic NLS on $\mathbb R$ is globally well-posed in $M^{2, p}(\mathbb R)$ for any $p < \infty$, while the renormalized cubic NLS on $\mathbb T$ is globally well-posed in $\mathcal FL^p(\mathbb T)$ for any $p < \infty$. In Appendix, we also establish analogous global-in-time bounds for the modified KdV equation (mKdV) in the modulation spaces on the real line and in the Fourier-Lebesgue spaces on the circle. An additional key ingredient of the proof in this case is a Galilean transform which converts the mKdV to the mKdV-NLS equation.


title: Experimentally Probing the Algorithmic Randomness and Incomputability of Quantum Randomness

abstract: The advantages of quantum random number generators (QRNGs) over pseudo-random number generators (PRNGs) are normally attributed to the nature of quantum measurements. This is often seen as implying the superiority of the sequences of bits themselves generated by QRNGs, despite the absence of empirical tests supporting this. Nonetheless, one may expect sequences of bits generated by QRNGs to have properties that pseudo-random sequences do not; indeed, pseudo-random sequences are necessarily computable, a highly nontypical property of sequences. In this paper, we discuss the differences between QRNGs and PRNGs and the challenges involved in certifying the quality of QRNGs theoretically and testing their output experimentally. While QRNGs are often tested with standard suites of statistical tests, such tests are designed for PRNGs and only verify statistical properties of a QRNG, but are insensitive to many supposed advantages of QRNGs. We discuss the ability to test the incomputability and algorithmic complexity of QRNGs. While such properties cannot be directly verified with certainty, we show how one can construct indirect tests that may provide evidence for the incomputability of QRNGs. We use these tests to compare various PRNGs to a QRNG, based on superconducting transmon qutrits and certified by the Kochen-Specker Theorem, to see whether such evidence can be found in practice. While our tests fail to observe a strong advantage of the quantum random sequences due to algorithmic properties, the results are nonetheless informative: some of the test results are ambiguous and require further study, while others highlight difficulties that can guide the development of future tests of algorithmic randomness and incomputability.


title: Election Score Can Be Harder Than Winner

abstract: Election systems based on scores generally determine the winner by computing the score of each candidate and the winner is the candidate with the best score. It would be natural to expect that computing the winner of an election is at least as hard as computing the score of a candidate. We show that this is not always the case. In particular, we show that for Young elections for dichotomous preferences the winner problem is easy, while determining the score of a candidate is hard. This complexity behavior has not been seen before and is unusual. The easiness of the winner problem for dichotomous Young crucially uses the fact that dichotomous preferences guarantee the transitivity of the majority relation. In addition to dichotomous preferences we also look at single-peaked preferences, the most well-studied domain restriction that guarantees the transitivity of the majority relation. We show that for the three major hard election systems and their natural variants, dichotomous Young is the only case where winner is easy and score is hard. This also solves an open question from Lackner and Peters (AAAI 2017), by providing a polynomial-time algorithm for Dodgson score for single-peaked electorates.


title: Learning Traffic Flow Dynamics using Random Fields

abstract: This paper presents a mesoscopic traffic flow model that explicitly describes the spatio-temporal evolution of the probability distributions of vehicle trajectories. The dynamics are represented by a sequence of factor graphs, which enable learning of traffic dynamics from limited Lagrangian measurements using an efficient message passing technique. The approach ensures that estimated speeds and traffic densities are non-negative with probability one. The estimation technique is tested using vehicle trajectory datasets generated using an independent microscopic traffic simulator and is shown to efficiently reproduce traffic conditions with probe vehicle penetration levels as little as 10\%. The proposed algorithm is also compared with state-of-the-art traffic state estimation techniques developed for the same purpose and it is shown that the proposed approach can outperform the state-of-the-art techniques in terms reconstruction accuracy.


title: Gravitational self-force corrections to tidal invariants for particles on circular orbits in a Kerr spacetime

abstract: We generalize to the Kerr spacetime existing self-force results on tidal invariants for particles moving along circular orbits around a Schwarzschild black hole. We obtain linear-in-mass-ratio corrections to the quadratic and cubic electric-type invariants and the quadratic magnetic-type invariant in series of the rotation parameter up to the fourth order. We then construct the eigenvalues of both electric and magnetic tidal tensors and analytically compute them through high post-Newtonian orders.


title: The $A_\infty$-structure of the index map

abstract: Let $F$ be a local field with residue field $k$. The classifying space of $GL_n(F)$ comes canonically equipped with a map to the delooping of the $K$-theory space of $k$. Passing to loop spaces, such a map abstractly encodes a homotopy coherently associative map of A-infinity-spaces $GL_n(F)\to K_k$. Using a generalized Waldhausen construction, we construct an explicit model built for the $A_\infty$-structure of this map, built from nested systems of lattices in $F^n$. More generally, we construct this model in the framework of Tate objects in exact categories, with finite dimensional vector spaces over local fields as a motivating example.


title: Filament flexibility enhances power transduction of F-actin bundles

abstract: The dynamic behavior of bundles of actin filaments growing against a loaded obstacle is investigated through a generalized version of the standard multi filaments Brownian Ratchet model in which the (de)polymerizing filaments are treated not as rigid rods but as semi-flexible discrete wormlike chains with a realistic value of the persistence length. By stochastic dynamic simulations we study the relaxation of a bundle of $N_f$ filaments with staggered seed arrangement against a harmonic trap load in supercritical conditions. Thanks to the time scale separation between the wall motion and the filament size relaxation, mimiking realistic conditions, this set-up allows us to extract a full load-velocity curve from a single experiment over the trap force/size range explored. We observe a systematic evolution of steady non-equilibrium states over three regimes of bundle lengths $L$. A first threshold length $\Lambda$ marks the transition between the rigid dynamic regime ($L<\Lambda$), characterized by the usual rigid filament load-velocity relationship $V(F)$, and the flexible dynamic regime ($L>\Lambda$), where the velocity $V(F,L)$ is an increasing function of the bundle length $L$ at fixed load $F$, the enhancement being the result of an improved level of work sharing among the filaments induced by flexibility. A second critical length corresponds to the beginning of an unstable regime characterized by a high probability to develop escaping filaments which start growing laterally and thus do not participate anymore to the generation of the polymerization force. This phenomenon prevents the bundle from reaching at this critical length the limit behavior corresponding to Perfect Load Sharing.


title: Coherency and incoherency in neutrino-nucleus elastic and inelastic scattering

abstract: Neutrino-nucleus scattering $\nu A\to \nu A$, in which the nucleus conserves its integrity, is considered. We show that elastic interactions keeping the nucleus in the same quantum state lead to a quadratic enhancement of the corresponding cross-section in terms of the number of nucleons. Meanwhile, the cross-section of inelastic processes in which the quantum state of the nucleus is changed, essentially has a linear dependence on the number of nucleons. These two classes of processes are referred to as coherent and incoherent, respectively. The coherent and incoherent cross-sections are driven by factors $|F_{p/n}|^2$ and $(1-|F_{p/n}|^2)$, where $|F_{p/n}|^2$ is a proton/neutron form-factor of the nucleus, averaged over its initial states. The coherent cross-section formula used in the literature is revised and corrections depending on kinematics are estimated. As an illustration of the importance of the incoherent channel we considered three experimental setups with different nuclei. Experiments attempting to measure coherent neutrino scattering by solely detecting the recoiling nucleus, as is typical, might be including an incoherent background that is indistinguishable from the signal if the excitation gamma eludes its detection. However, as is shown here, the incoherent component can be measured directly by searching for photons released by the excited nuclei inherent to the incoherent channel. For a beam experiment these gammas should be correlated in time with the beam, and their higher energies make the corresponding signal easily detectable at a rate governed by the ratio of incoherent to coherent cross-sections. The detection of signals due to the nuclear recoil and excitation gammas provides a more sensitive instrument in studies of nuclear structure and possible signs of new physics.


title: Gauge Field and Fermion Production during Axion Inflation

abstract: We study the dual production of helical Abelian gauge fields and chiral fermions through the Chern-Simons (CS) coupling with a pseudo-scalar inflaton in the presence of a chiral anomaly. Through the CS term, the motion of the inflaton induces a tachyonic instability for one of the two helicities of the gauge field. We show that the resulting helical gauge field necessarily leads to the production of chiral fermions by deforming their Fermi sphere into discrete Landau levels. The population of the lowest Landau level leads to a chiral asymmetry as inferred from the chiral anomaly, while the higher levels are populated symmetrically through pair production. From the backreaction of the fermions on the gauge field production we derive a conservative but stringent upper bound on the magnitude of the gauge fields. Consequently, we find that the scalar perturbations sourced by these helical gauge fields, responsible for enhanced structure formation on small scales, get reduced significantly. We also discuss the fate of the primordial chiral asymmetry and of the helical gauge fields after inflation, and show that the instability in the chiral plasma tends to erase these primordial asymmetries. This result may impact scenarios where the baryon asymmetry of the Universe is connected to primordial magnetic fields.


title: Uniform 2D-Monotone Minimum Spanning Graphs

abstract: A geometric graph $G$ is $xy-$monotone if each pair of vertices of $G$ is connected by a $xy-$monotone path. We study the problem of producing the $xy-$monotone spanning geometric graph of a point set $P$ that (i) has the minimum cost, where the cost of a geometric graph is the sum of the Euclidean lengths of its edges, and (ii) has the least number of edges, in the cases that the Cartesian System $xy$ is specified or freely selected. Building upon previous results, we easily obtain that the two solutions coincide when the Cartesian System is specified and are both equal to the rectangle of influence graph of $P$. The rectangle of influence graph of $P$ is the geometric graph with vertex set $P$ such that two points $p,q \in P$ are adjacent if and only if the rectangle with corners $p$ and $q$ does not include any other point of $P$. When the Cartesian System can be freely chosen, we note that the two solutions do not necessarily coincide, however we show that they can both be obtained in $O(|P|^3)$ time. We also give a simple $2-$approximation algorithm for the problem of computing the spanning geometric graph of a $k-$rooted point set $P$, in which each root is connected to all the other points (including the other roots) of $P$ by $y-$monotone paths, that has the minimum cost.


title: What Does This Notation Mean Anyway?

abstract: Following the introduction of BNF notation by Backus for the Algol 60 report and subsequent notational variants, a metalanguage involving formal "grammars" has developed for discussing structured objects in Computer Science and Mathematical Logic. We refer to this offspring of BNF as Math-BNF or MBNF, to the original BNF and its notational variants just as BNF, and to aspects common to both as BNF-style. What all BNF-style notations share is the use of production rules roughly of this form: $$\bullet \mathrel{::=} \circ_1 \mid \cdots \mid \circ_n $$ Normally, such a rule says "every instance of $\circ_i$ for $i \in \{1, \ldots, n\}$ is also an in stance of $\bullet$". MBNF is distinct from BNF in the entities and operations it allows. Instead of strings, MBNF builds arrangements of symbols that we call math-text. Sometimes "syntax" is defined by interleaving MBNF production rules and other mathematical definitions that can contain chunks of math-text. There is no clear definition of MBNF. Readers do not have a document which tells them how MBNF is to be read and must learn MBNF through a process of cultural initiation. To the extent that MBNF is defined, it is largely through examples scattered throughout the literature and which require readers to guess the mathematical structures underpinning them. This paper gives MBNF examples illustrating some of the differences between MBNF and BNF. We propose a definition of syntactic math text (SMT) which handles many (but far from all) uses of math-text and MBNF in the wild. We aim to balance the goal of being accessible and not requiring too much prerequisite knowledge with the conflicting goal of providing a rich mathematical structure that already supports many uses and has possibilities to be extended to support more challenging cases.


title: Inverse problems for Maxwell's equations in a slab with partial boundary data

abstract: We consider two inverse boundary value problems for the time-harmonic Maxwell equations in an infinite slab. Assuming that tangential boundary data for the electric and magnetic fields at a fixed frequency is available either on subsets of one boundary hyperplane, or on subsets of different boundary hyperplanes, we show that the electromagnetic material parameters, the conductivity, electric permittivity, and magnetic permeability, are uniquely determined by these partial measurements.


title: Problems, Questions, and Conjectures about Mapping Class Groups

abstract: We discuss a number of open problems about mapping class groups of surfaces. In particular, we discuss problems related to linearity, congruence subgroups, cohomology, pseudo-Anosov stretch factors, Torelli subgroups, and normal subgroups.


title: Impact of detuning and dephasing on a laser-corrected subnatural-linewidth single-photon source

abstract: The elastic scattering peak of a resonantly driven two-level system has been argued to provide narrow-linewidth antibunched photons. Although independent measurements of spectral width on the one hand and antibunching on the other hand do seem to show that this is the case, a joint measurement reveals that only one or the other of these attributes can be realised in the direct emission. We discuss a scheme which interferes the emission with a laser to produce simultaneously single photons of subnatural linewidth. In particular, we consider the effect of dephasing and of the detuning between the driving laser and/or the detector with the emitter. We find that also in presence of dephasing, our scheme brings considerable improvement as compared to the standard scheme.


title: CVC4 at the SMT Competition 2018

abstract: This paper is a description of the CVC4 SMT solver as entered into the 2018 SMT Competition. We only list important differences from the 2017 SMT Competition version of CVC4. For further and more detailed information about CVC4, please refer to the original paper, the CVC4 website, or the source code on GitHub.


title: Age of Information and Throughput in a Shared Access Network with Heterogeneous Traffic

abstract: We consider a cognitive shared access scheme consisting of a high priority primary node and a low priority network with $N$ secondary nodes accessing the spectrum. Assuming bursty traffic at the primary node, saturated queues at the secondary nodes, and multipacket reception capabilities at the receivers, we derive analytical expressions of the time average age of information of the primary node and the throughput of the secondary nodes. We formulate two optimization problems, the first aiming to minimize the time average age of information of the primary node subject to an aggregate secondary throughput requirement. The second problem aims to maximize the aggregate secondary throughput of the network subject to a maximum time average staleness constraint. Our results provide guidelines for the design of a multiple access system with multipacket reception capabilities that fulfils both timeliness and throughput requirements.


title: Wireless Channel Dynamics and Robustness for Ultra-Reliable Low-Latency Communications

abstract: Interactive, immersive and critical applications demand ultra-reliable low-latency communication (URLLC). To build wireless communication systems that can support these applications, understanding the characteristics of the wireless medium is paramount. Although wireless channel characteristics and dynamics have been extensively studied, it is important to revisit these concepts in the context of the strict demands of low latency and ultra-reliability. In this paper, we bring a modeling approach from robust control to wireless communication -- the wireless channel characteristics are given a nominal model around which we allow for some quantified uncertainty. We propose certain key "directions" along which to bound model uncertainty that are relevant to URLLC. For the nominal model, we take an in-depth look at wireless channel characteristics such as spatial and temporal correlations based on Jakes' model. Contrary to what has been claimed in the literature, we find that standard Rayleigh fading processes are not bandlimited. This has significant implications on the predictability of channels. We also find that under reasonable conditions the spatial correlation of channels provide a fading distribution that is not too far off from an independent spatial fading model. Additionally, we look at the impact of these channel models on cooperative communication based systems. We find that while spatial-diversity-based techniques are necessary to combat the adverse effects of fading, time-diversity-based techniques are necessary to be robust against unmodeled errors. Robust URLLC systems need to operate with both an SNR margin and a time/repetition margin.


title: Reduced Scaling Hilbert Space Variational Monte Carlo

abstract: We show that for both single-Slater-Jastrow and Jastrow geminal power wave functions, the formal cost scaling of Hilbert space variational Monte Carlo can be reduced from fifth to fourth order in the system size, thus bringing it in line with the long-standing scaling of its real space counterpart. While traditional quantum chemistry methods can reduce costs related to the two-electron integral tensor through resolution of the identity and Cholesky decomposition approaches, we show that such approaches are ineffective in the presence of Hilbert space Jastrow factors. Instead, we develop a simple semi-stochastic approach that can take similar advantage of the near-sparsity of this four-index tensor. Through demonstrations on alkanes of increasing length, we show that accuracy and overall statistical uncertainty are not meaningfully affected and that a total cost crossover is reached as early as 50 electrons.


title: Quantum current in dissipative systems

abstract: Describing current in open quantum systems can be problematic due to the subtle interplay of quantum coherence and environmental noise. Probing the noise-induced current can be detrimental to the tunneling-induced current and vice versa. We derive a general theory for the probability current in quantum systems arbitrarily interacting with their environment that overcomes this difficulty. We show that the current can be experimentally measured by performing a sequence of weak and standard quantum measurements. We exemplify our theory by analyzing a simple Smoluchowski-Feynman-type ratchet consisting of two particles, operating deep in the quantum regime. Fully incorporating both thermal and quantum effects, the current generated in the model can be used to detect the onset of "genuine quantumness" in the form of quantum contextuality. The model can also be used to generate steady-state entanglement in the presence of arbitrarily hot environment.

categories: quant-ph cond-mat.mes-hall cond-mat.other cond-mat.stat-mech


title: Subsystem symmetries, quantum cellular automata, and computational phases of quantum matter

abstract: Quantum phases of matter are resources for notions of quantum computation. In this work, we establish a new link between concepts of quantum information theory and condensed matter physics by presenting a unified understanding of symmetry-protected topological (SPT) order protected by subsystem symmetries and its relation to measurement-based quantum computation (MBQC). The key unifying ingredient is the concept of quantum cellular automata (QCA) which we use to define subsystem symmetries acting on rigid lower-dimensional lines or fractals on a 2D lattice. Notably, both types of symmetries are treated equivalently in our framework. We show that states within a non-trivial SPT phase protected by these symmetries are indicated by the presence of the same QCA in a tensor network representation of the state, thereby characterizing the structure of entanglement that is uniformly present throughout these phases. By also formulating schemes of MBQC based on these QCA, we are able to prove that most of the phases we construct are computationally universal phases of matter, in which every state is a resource for universal MBQC. Interestingly, our approach allows us to construct computational phases which have practical advantages over previous examples, including a computational speedup. The significance of the approach stems from constructing novel computationally universal phases of matter and showcasing the power of tensor networks and quantum information theory in classifying subsystem SPT order.


title: Quantum Codes from Neural Networks

abstract: We examine the usefulness of applying neural networks as a variational state ansatz for many-body quantum systems in the context of quantum information-processing tasks. In the neural network state ansatz, the complex amplitude function of a quantum state is computed by a neural network. The resulting multipartite entanglement structure captured by this ansatz has proven rich enough to describe the ground states and unitary dynamics of various physical systems of interest. In the present paper, we initiate the study of neural network states in quantum information-processing tasks. We demonstrate that neural network states are capable of efficiently representing quantum codes for quantum information transmission and quantum error correction, supplying further evidence for the usefulness of neural network states to describe multipartite entanglement. In particular, we show the following main results: a) Neural network states yield quantum codes with a high coherent information for two important quantum channels, the generalized amplitude damping channel and the dephrasure channel. These codes outperform all other known codes for these channels, and cannot be found using a direct parametrization of the quantum state. b) For the depolarizing channel, the neural network state ansatz reliably finds the best known codes given by repetition codes. c) Neural network states can be used to represent absolutely maximally entangled states, a special type of quantum error-correcting codes. In all three cases, the neural network state ansatz provides an efficient and versatile means as a variational parametrization of these highly entangled states.


title: Finding Local Minima via Stochastic Nested Variance Reduction

abstract: We propose two algorithms that can find local minima faster than the state-of-the-art algorithms in both finite-sum and general stochastic nonconvex optimization. At the core of the proposed algorithms is $\text{One-epoch-SNVRG}^+$ using stochastic nested variance reduction (Zhou et al., 2018a), which outperforms the state-of-the-art variance reduction algorithms such as SCSG (Lei et al., 2017). In particular, for finite-sum optimization problems, the proposed $\text{SNVRG}^{+}+\text{Neon2}^{\text{finite}}$ algorithm achieves $\tilde{O}(n^{1/2}\epsilon^{-2}+n\epsilon_H^{-3}+n^{3/4}\epsilon_H^{-7/2})$ gradient complexity to converge to an $(\epsilon, \epsilon_H)$-second-order stationary point, which outperforms $\text{SVRG}+\text{Neon2}^{\text{finite}}$ (Allen-Zhu and Li, 2017) , the best existing algorithm, in a wide regime. For general stochastic optimization problems, the proposed $\text{SNVRG}^{+}+\text{Neon2}^{\text{online}}$ achieves $\tilde{O}(\epsilon^{-3}+\epsilon_H^{-5}+\epsilon^{-2}\epsilon_H^{-3})$ gradient complexity, which is better than both $\text{SVRG}+\text{Neon2}^{\text{online}}$ (Allen-Zhu and Li, 2017) and Natasha2 (Allen-Zhu, 2017) in certain regimes. Furthermore, we explore the acceleration brought by third-order smoothness of the objective function.


title: Exploring $B\to\pi\pi, \pi K$ Decays at the High-Precision Frontier

abstract: The $B\to\pi\pi,\pi K$ system offers a powerful laboratory to probe strong and weak interactions. Using the isospin symmetry, we determine hadronic $B\to\pi\pi$ parameters from data where new measurements of direct CP violation in $B^0_d\to\pi^0\pi^0$ resolve a discrete ambiguity. With the help of the $SU(3)$ flavour symmetry, the $B\to\pi\pi$ parameters can be converted into their $B\to\pi K$ counterparts, thereby allowing us to make predictions of observables. A particularly interesting decay is $B^0_d\to\pi^0 K_{\rm S}$ as it exhibits mixing-induced CP violation. Using an isospin relation, complemented with a robust $SU(3)$ input, we calculate correlations between the direct and mixing-induced CP asymmetries of $B^0_d\to\pi^0 K_{\rm S}$, which are the theoretically cleanest $B\to\pi K$ probes. Interestingly, they show tensions with respect to the Standard Model. Should this $B\to\pi K$ puzzle originate from New Physics, electroweak penguins offer an attractive scenario for new particles to enter. We present a strategy to determine the parameters characterising these topologies and obtain the state-of-the-art picture from current data. In the future, this method will allow us to reveal the $B\to\pi K$ dynamics and to obtain insights into the electroweak penguin sector with unprecedented precision.


title: Local unambiguous discrimination of symmetric ternary states

abstract: We investigate unambiguous discrimination between given quantum states with a sequential measurement, which is restricted to local measurements and one-way classical communication. If the given states are binary or those each of whose individual systems is two-dimensional, then it is in some cases known whether a sequential measurement achieves a globally optimal unambiguous measurement. In contrast, for more than two states each of whose individual systems is more than two-dimensional, the problem becomes extremely complicated. This paper focuses on symmetric ternary pure states each of whose individual systems is three-dimensional, which include phase shift keyed (PSK) optical coherent states and a lifted version of "double trine" states. We provide a necessary and sufficient condition for an optimal sequential measurement to be globally optimal for the bipartite case. A sufficient condition of global optimality for multipartite states is also presented. One can easily judge whether these conditions hold for given states. Some examples are given, which demonstrate that, despite the restriction to local measurements and one-way classical communication, a sequential measurement can be globally optimal in quite a few cases.


title: CoHSI I; Detailed properties of the Canonical Distribution for Discrete Systems such as the Proteome

abstract: The CoHSI (Conservation of Hartley-Shannon Information) distribution is at the heart of a wide-class of discrete systems, defining the length distribution of their components amongst other global properties. Discrete systems such as the known proteome where components are proteins, computer software, where components are functions and texts where components are books, are all known to fit this distribution accurately. In this short paper, we explore its solution and its resulting properties and lay the foundation for a series of papers which will demonstrate amongst other things, why the average length of components is so highly conserved and why long components occur so frequently in these systems. These properties are not amenable to local arguments such as natural selection in the case of the proteome or human volition in the case of computer software, and indeed turn out to be inevitable global properties of discrete systems devolving directly from CoHSI and shared by all. We will illustrate this using examples from the Uniprot protein database as a prelude to subsequent studies.


title: Projective geometry in the Poincar\'e disk of a $C^*$-algebra

abstract: We study the Poincar\'e disk ${\cal D}=\{a\in {\cal A}: \|a\|<1\}$ of a C$^*$-algebra ${\cal A}$ from a projective point of view: ${\cal D}$ is regarded as an open subset of the projective line $\mathbb{P}_1{\cal A}$, the space of complemented rank one submodules of ${\cal A}^2$. We introduce the concept of cross ratio of four points in $\mathbb{P}_1{\cal A}$. Our main result establishes the relation between the exponential map $Exp_{z_0}(z_1)$ of ${\cal D}$ ($z_0,z_1\in {\cal D}$) and the cross ratio of the four-tuple $$ \delta(-\infty), \delta(0)=z_0, \delta(1)=z_1 , \delta(+\infty), $$ where $\delta$ is the unique geodesic of ${\cal D}$ joining $z_0$ and $z_1$ at times $t=0$ and $t=1$, respectively.


title: A proposed solution for analysis management in high energy physics

abstract: This paper presents an architecture for the analysis management in high energy physics experiments. Some new concepts on data analysis are introduced. A protocol for organizing and operating an analysis is raised. A toolkit following this architecture is developed, which provides a solution of analysis management with both flexibility and reproducibility. A foreseen development of this toolkit is discussed.


title: On the Structure and Function of Scientific Perspectivism in Quantum Mechanics

abstract: Contemporary scientific perspectivism is re-evaluated and extended to a comprehensive perspectivist methodology and 'mediated' realistic epistemology, especially, with reference to quantum mechanics. In the present study, this is realized by representing categorically the global structure of a quantum algebra of events in terms of structured multitudes of interrelated local Boolean frames, realized as suitable perspectives or contexts for measuring physical quantities. The philosophical meaning of the proposed approach implies that the quantum world can be consistently approached and comprehended through a multilevel structure of locally variable perspectives, which interlock, in a category-theoretical environment, to form a coherent picture of the whole in a nontrivial way. Thus, in contrast to a panoptical "view from nowhere" of the classical paradigm, quantum theory acknowledges in an essential way a perspectival/contextual character of scientific knowledge.


title: Physics of muscle contraction

abstract: In this paper we report, clarify and broaden various recent efforts to complement the chemistry-centered models of force generation in muscles by mechanics-centered models. The physical mechanisms of interest can be grouped into two classes: passive and active. The main passive effect is the fast force recovery which does not require the detachment of myosin cross-bridges from actin filaments and can operate without a specialized supply ATP. In mechanical terms, it can be viewed as a collective folding-unfolding phenomenon in the system of interacting bistable units and modeled by near equilibrium Langevin dynamics. The parallel active force generation mechanism operates at slow timescales, requires detachment and is crucially dependent on ATP hydrolysis. The underlying mechanical processes take place far from equilibrium and are represented by stochastic models with broken time reversal symmetry implying non-potentiality, correlated noise or multiple reservoirs. The modeling approaches reviewed in this paper deal with both active and passive processes and support from the mechanical perspective the biological point of view that phenomena involved in slow (active) and fast (passive) force generation are tightly intertwined. They reveal, however, that biochemical studies in solution, macroscopic physiological measurements and structural analysis do not provide by themselves all the necessary insights into the functioning of the contractile system. In particular, the reviewed body of work emphasizes the important role of long-range interactions and criticality in securing the targeted mechanical response in the physiological regime of isometric contractions. The importance of the purely mechanical microscale modeling is accentuated at the end of the paper where we address the puzzling issue of the stability of muscle response on the so called descending limb of the isometric tetanus.


title: The Physics of baking good Pizza

abstract: Physical principles are involved in almost any aspect of cooking. Here we analyze the specific process of baking pizzas, deriving in simple terms the baking times for two different situations: For a brick oven in a pizzeria and a modern metallic oven at home. Our study is based on basic thermodynamic principles relevant to the cooking process and is accessible to undergraduate students. We start with a historical overview of the development and art of pizza baking, illustrate the underlying physics by some simple common examples, and then apply them in detail to the example of baking pizza.


title: The Hierarchion, a Relaxion Addressing the Standard Model's Hierarchies

abstract: We present a mechanism that addresses the electroweak, the strong CP, and the flavor hierarchies of the Standard Model (including neutrino masses) in a unified way. The naturalness of the electroweak scale is solved together with the strong CP problem by the Nelson-Barr relaxion: the relaxion field is identified with the pseudo-Nambu-Goldstone boson of an abelian symmetry with no QCD anomaly. The Nelson-Barr sector generates the "rolling" potential and the relaxion vacuum expectation value at the stopping point is mapped to the Cabibbo-Kobayashi-Maskawa phase. An abelian symmetry accounts for the Standard Model's mass hierarchies and flavor textures through the Froggatt-Nielsen mechanism. We show how the "backreaction" potential of the relaxion can be induced by a sterile neutrino sector, without any extra state with electroweak quantum numbers. The same construction successfully explains neutrino masses and mixings. The only light field in our construction is the relaxion, which we call the hierarchion because it is essentially linked to our construct that accounts for all the Standard Model hierarchies. Given its interplay with flavor symmetries, the hierarchion can be probed in flavor-violating decays of the Standard Model fermions, motivating a further experimental effort in looking for new physics in rare decays of leptons and mesons.


title: Relative Binning and Fast Likelihood Evaluation for Gravitational Wave Parameter Estimation

abstract: We present a method to accelerate the evaluation of the likelihood in gravitational wave parameter estimation. Parameter estimation codes compute likelihoods of similar waveforms, whose phases and amplitudes differ smoothly with frequency. We exploit this by precomputing frequency-binned overlaps of the best-fit waveform with the data. We show how these summary data can be used to approximate the likelihood of any waveform that is sufficiently probable within the required accuracy. We demonstrate that $\simeq 60$ bins suffice to accurately compute likelihoods for strain data at a sampling rate of $4096\,$Hz and duration of $T=2048\,$s around the binary neutron star merger GW170817. Relative binning speeds up parameter estimation for frequency domain waveform models by a factor of $\sim 10^4$ compared to naive matched filtering and $\sim 10$ compared to reduced order quadrature.


title: Parameter Estimation for GW170817 using Relative Binning

abstract: Relative binning is a new method for fast and accurate evaluation of the likelihood of gravitational wave strain data. This technique can be used to produce reliable posterior distributions for compact object mergers with very moderate computational resources. We use a fast likelihood evaluation code based on this technique to estimate the parameters of the double neutron-star merger event GW170817 using publicly available LIGO data. We obtain statistically similar posteriors using either Markov-chain Monte-Carlo or nested sampling. The results do not favor non-zero aligned spins at a statistically significant level. There is no significant sign of non-zero tidal deformability (as quantified by the Bayesian evidence), whether or not high-spin or low-spin priors are adopted. Our posterior samples are publicly available, and we also provide a tutorial Python code to implement fast likelihood evaluation using the relative binning method.


title: Optical and Infrared Photometry of the Nearby SN 2017cbv

abstract: On 2017 March 11, the DLT40 Transient Discovery Survey discovered SN 2017cbv in NGC5643, a Type 2 Seyfert Galaxy in the Lupus Constellation. SN 2017cbv went on to become a bright Type Ia supernova, with a $V_{max}$ of 11.51 $\pm$ 0.05 mag. We present early time optical and infrared photometry of SN 2017cbv covering the rise and fall of over 68 days. We find that SN 2017cbv has a broad light curve $\Delta m_{15}(B)$ = 0.88 $\pm$ 0.07, a $B$-band maximum at 2457840.97 $\pm$ 0.43, a negligible host galaxy reddening where $E(B-V)_{host}$ $\approx$ 0, and a distance modulus of 30.49 $\pm$ 0.32 to the SN, corresponding to a distance of $12.58_{-1.71}^{+1.98}$ Mpc. We also present the results of two different numerical models we used for analysis in this paper: SALT2, an empirical model for Type Ia supernova optical light curves that accounts for variability components; and SNooPy, the CSP-II light-curve model that covers both optical and near-infrared wavelengths and is used for distance estimates.


title: Topological superconductivity in Dirac honeycomb systems

abstract: We predict two topological superconducting phases in microscopic models arising from the Berry phase associated with the valley degree of freedom in gapped Dirac honeycomb systems. The first one is a topological helical spin-triplet superconductor with a nonzero center-of-mass momentum that does not break time-reversal symmetry. We also find a topological chiral-triplet superconductor with Chern number $\pm 1$ with equal-spin-pairing in one valley and opposite-spin-triplet pairing in the other valley. Our results are obtained for the Kane-Mele model in which we have explored the effect of three different interactions, onsite attraction $U$, nearest-neighbor density-density attraction $V$, and nearest-neighbor antiferromagnetic exchange $J$, within self-consistent Bogoliubov--de Gennes theory. Transition metal dichalcogenides and cold atom experiments are promising platforms to explore these phases.


title: In-medium loop corrections and longitudinally polarized gauge bosons in high-energy showers

abstract: The splitting processes of bremsstrahlung and pair production in a medium are coherent over large distances in the very high energy limit, which leads to a suppression known as the Landau-Pomeranchuk-Migdal (LPM) effect. We continue study of the case when the coherence lengths of two consecutive splitting processes overlap (which is important for understanding corrections to standard treatments of the LPM effect in QCD), avoiding soft-emission approximations. In this particular paper, we show (i) how the "instantaneous" interactions of Light-Cone Perturbation Theory must be included in the calculation to account for effects of longitudinally-polarized gauge bosons in intermediate states, and (ii) how to compute virtual corrections to LPM emission rates, which will be necessary in order to make infrared-safe calculations of the characteristics of in-medium QCD showering of high-energy partons. In order to develop these topics in as simple a context as possible, we will focus in the current paper not on QCD but on large-$N_f$ QED, where $N_f$ is the number of electron flavors.


title: Gauss's Law, Duality, and the Hamiltonian Formulation of U(1) Lattice Gauge Theory

abstract: Quantum computers have the potential to explore the vast Hilbert space of entangled states that play an important role in the behavior of strongly interacting matter. This opportunity motivates reconsidering the Hamiltonian formulation of gauge theories, with a suitable truncation scheme to render the Hilbert space finite-dimensional. Conventional formulations lead to a Hilbert space largely spanned by unphysical states; given the current inability to perform large scale quantum computations, we examine here how one might restrict wave function evolution entirely or mostly to the physical subspace. We consider such constructions for the simplest of these theories containing dynamical gauge bosons -- U(1) lattice gauge theory without matter in $d=2,3$ spatial dimensions -- and find that electric-magnetic duality naturally plays an important role. We conclude that this approach is likely to significantly reduce computational overhead in $d=2$ by a reduction of variables and by allowing one to regulate magnetic fluctuations instead of electric. The former advantage does not exist in $d=3$, but the latter might be important for asymptotically-free gauge theories.


title: z ~ 2.5 - 3 Ionizers in the GOODS-N Field

abstract: We use deep F275W imaging from the $Hubble$ Deep UV Legacy Survey (HDUV) and G280 grism spectroscopy from $HST$/WFC3, along with new and archival optical spectra from Keck/DEIMOS, to search for candidate ionizing sources in the GOODS-N field at z ~ 2.5 - 3. Spectroscopic identification of our UV-selected sources are 99% complete to F275W = 25.5 in the region of the UV imaging, and we identify 6 potential ionizing galaxies or AGNs at z ~ 3. By far the brightest of these is a z = 2.583 AGN that totally dominates the ionizing flux in the region, with a specific ionizing volume emissivity at $912~ \AA$ of $\epsilon_{912}$ = $8.3^{27}_{1.4} \times 10^{24}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$. Based on our spectroscopic data, we find four candidates are contaminated by foreground galaxies at z ~ 0.5 - 0.7. At $\epsilon_{912}$ = $2.2^{7.2}_{0.4} \times 10^{23}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$, the remaining candidate galaxy's contribution to the ionizing background lies well below the flux required to ionize the intergalactic medium at z ~ 2.5 - 3, consistent with previous observations that show AGNs provide the bulk of the ionizing background at these redshifts.


title: An Excess of Jupiter Analogs in Super-Earth Systems

abstract: We use radial velocity observations to search for long-period gas giant companions in systems hosting inner super-Earth (1-4 R_Earth, 1-10 M_Earth) planets to constrain formation and migration scenarios for this population. We consistently re-fit published RV datasets for 65 stars and find 9 systems with statistically significant trends indicating the presence of an outer companion. We combine these RV data with AO images to constrain the masses and semi-major axes of these companions. We quantify our sensitivity to the presence of long-period companions by fitting the sample with a power law distribution and find an occurrence rate of 39+/-7% for companions 0.5-20 M_Jup and 1-20 AU. Half of our systems were discovered by the transit method and half were discovered by the RV method. While differences in RV baselines and number of data points between the two samples lead to different sensitivities to distant companions, we find that occurrence rates of gas giant companions in each sample are consistent at the 0.5$\sigma$ level. We compare the frequency of Jupiter analogs in these systems to the equivalent rate from field star surveys and find that Jupiter analogs are more common around stars hosting super-Earths. We conclude that the presence of outer gas giants does not suppress the formation of inner super-Earths, and that these two populations of planets instead appear to be correlated. We also find that the stellar metallicities of systems with gas giant companions are higher than those without companions, in agreement with the well-established metallicity correlation from RV surveys of field stars.


title: The Disk Wind in the Neutron Star Low-mass X-Ray Binary GX 13+1

abstract: We present the analysis of seven \emph{Chandra} High Energy Transmission Grating Spectrometer and six simultaneous \emph{RXTE} Proportional Counter Array observations of the persistent neutron star (NS) low-mass X-ray binary GX 13+1 on its normal and horizontal branches. Across nearly 10 years, GX 13+1 is consistently found to be accreting at $50-70$\% Eddington, and all observations exhibit multiple narrow, blueshifted absorption features, the signature of a disk wind, despite the association of normal and horizontal branches with jet activity. A single absorber with standard abundances cannot account for all seven major disk wind features, indicating that multiple absorption zones may be present. Two or three absorbers can produce all of the absorption features at their observed broadened widths and reveal that multiple kinematic components produce the accretion disk wind signature. Assuming the most ionized absorber reflects the physical conditions closest to the NS, we estimate a wind launching radius of $7\times10^{10}$ cm, for an electron density of $10^{12}$ cm$^{-3}$. This is consistent with the Compton radius and also with a thermally driven wind. Because of the source's high Eddington fraction, radiation pressure likely facilitates the wind launching.


title: Galaxies Probing Galaxies in PRIMUS - II. The Coherence Scale of the Cool Circumgalactic Medium

abstract: The circumgalactic medium (CGM) close to ~L* star-forming galaxies hosts strong MgII 2796 absorption (with equivalent width W_2796>0.1 Ang) with a near-unity covering fraction. To characterize the spatial coherence of this absorption, we analyze the W_2796 distribution in the CGM of 27 star-forming galaxies detected in deep spectroscopy of bright background (b/g) galaxies first presented in Rubin et al. (2018). The sample foreground (f/g) systems have redshifts 0.35<z<0.8 and stellar masses 9.1<log M_*/M_sun<11.1, and the b/g galaxies provide spatially-extended probes with half-light radii 1.0 kpc<R_eff<7.9 kpc at projected distances R_perp<50 kpc. Our analysis also draws on literature W_2796 values measured in b/g QSO spectroscopy probing the halos of f/g galaxies with a similar range in M_* at z ~ 0.25. By making the assumptions that (1) samples of like galaxies exhibit similar circumgalactic W_2796 distributions; and that (2) the quantity log W_2796 has a Gaussian distribution with a dispersion that is constant with M_* and R_perp, we use this QSO-galaxy pair sample to construct a model for the log W_2796 distribution in the CGM. We then demonstrate the dependence of this distribution on the ratio of the surface area of the b/g probe to the projected absorber surface area (x_A=A_G/A_A), finding that distributions which assume x_A>=15 are statistically inconsistent with that observed toward our b/g galaxies at a 95% confidence level. This limit, in combination with the b/g galaxy sizes, requires that the length scale over which W_2796 does not vary (the "coherence scale" of MgII absorption) is l_A>1.9 kpc. This novel constraint on the morphology of cool, photoionized structures in the inner CGM suggests that either these structures each extend over kiloparsec scales, or that the numbers and velocity dispersion of these structures are spatially correlated over the same scales.


title: Empty gaps? Depleting annular regions in debris discs by secular resonance with a two-planet system

abstract: We investigate the evolution on secular time-scales of a radially extended debris disc under the influence of a system of two coplanar planets interior to the disc, showing that the secular resonances of the system can produce a depleted region in the disc by exciting the eccentricities of planetesimals. Using Laplace-Lagrange theory, we consider how the two exterior secular resonance locations, time-scales and widths depend on the masses, semi-major axes and eccentricities of the planets. In particular, we find that unless the resonances are very close to each other, one of them is very narrow and therefore unimportant for determining the observable structure of the disc. We apply these considerations to the debris disc of HD 107146, identifying combinations of the parameters of a possible unobserved two-planet system that could configure the secular resonances appropriately to reproduce the depletion observed in the disc. By performing N-body simulations of such systems, we find that planetesimal eccentricities do indeed become large near the theoretical secular resonance locations. The N-body output is post-processed to set the initial surface density profile of the disc, and to include the possible effects of collisional depletion. We find that it is possible to obtain a double-ringed disc in these simulations but not an axisymmetric one, with the inner ring having an offset whose magnitude depends on the eccentricities of the planets, and the outer ring showing spiral structure.


title: Higher-order stationary dispersive equations on bounded intervals: a relation between the order of an equation and the growth of its convective term

abstract: A boundary value problem for a stationary nonlinear dispersive equation of order $2l+1\;\; l\in \mathbb{N}$ with a convective term in the form $u^ku_x\;\; k\in \mathbb{N}$ was considered on an interval $(0,L)$. The existence, uniqueness and continuous dependence of a regular solution as well as a relation between $l$ and critical values of $k$ have been established.


title: Hierarchical Graph Representation Learning with Differentiable Pooling

abstract: Recently, graph neural networks (GNNs) have revolutionized the field of graph representation learning through effectively learned node embeddings, and achieved state-of-the-art results in tasks such as node classification and link prediction. However, current GNN methods are inherently flat and do not learn hierarchical representations of graphs---a limitation that is especially problematic for the task of graph classification, where the goal is to predict the label associated with an entire graph. Here we propose DiffPool, a differentiable graph pooling module that can generate hierarchical representations of graphs and can be combined with various graph neural network architectures in an end-to-end fashion. DiffPool learns a differentiable soft cluster assignment for nodes at each layer of a deep GNN, mapping nodes to a set of clusters, which then form the coarsened input for the next GNN layer. Our experimental results show that combining existing GNN methods with DiffPool yields an average improvement of 5-10% accuracy on graph classification benchmarks, compared to all existing pooling approaches, achieving a new state-of-the-art on four out of five benchmark data sets.


title: PCA of high dimensional random walks with comparison to neural network training

abstract: One technique to visualize the training of neural networks is to perform PCA on the parameters over the course of training and to project to the subspace spanned by the first few PCA components. In this paper we compare this technique to the PCA of a high dimensional random walk. We compute the eigenvalues and eigenvectors of the covariance of the trajectory and prove that in the long trajectory and high dimensional limit most of the variance is in the first few PCA components, and that the projection of the trajectory onto any subspace spanned by PCA components is a Lissajous curve. We generalize these results to a random walk with momentum and to an Ornstein-Uhlenbeck processes (i.e., a random walk in a quadratic potential) and show that in high dimensions the walk is not mean reverting, but will instead be trapped at a fixed distance from the minimum. We finally compare the distribution of PCA variances and the PCA projected training trajectories of a linear model trained on CIFAR-10 and ResNet-50-v2 trained on Imagenet and find that the distribution of PCA variances resembles a random walk with drift.


title: Domination and regularity

abstract: We discuss the close relationship between structural theorems in generalized stability theory and graph regularity theorems.


title: Estimating the treatment effect in a subgroup defined by an early post-baseline biomarker measurement in randomized clinical trials with time-to-event endpoint

abstract: Biomarker measurements can be relatively easy and quick to obtain and they are useful to investigate whether a compound works as intended on a mechanistic, pharmacological level. In some situations, it is realistic to assume that patients, whose post-baseline biomarker levels indicate that they do not sufficiently respond to the drug, are also unlikely to respond on clinically relevant long term outcomes (such as time-to-event). However the determination of the treatment effect in the subgroup of patients that sufficiently respond to the drug according to their biomarker levels is not straightforward: It is unclear which patients on placebo would have responded had they been given the treatment, so that naive comparisons between treatment and placebo will not estimate the treatment effect of interest. The purpose of this paper is to investigate assumptions necessary to obtain causal conclusions in such a setting, utilizing the formalism of causal inference. Three approaches for estimation of subgroup effects will be developed and illustrated using simulations and a case-study.


title: Complexity Matching and Requisite Variety

abstract: Complexity matching characterizes the role of information in interactions between systems and can be traced back to the 1957 Introduction to Cybernetics by Ross Ashby. We argue that complexity can be expressed in terms of crucial events, which are generated by the processes of spontaneous self-organization. Complex processes, ranging from biological to sociological, must satisfy the homeodynamic condition and host crucial events that have recently been shown to drive the information transport between complex systems. We adopt a phenomenological approach, based on the subordination to periodicity that makes it possible to combine homeodynamics and self-organization induced crucial events. The complexity of crucial events is defined by the waiting-time probability density function (PDF) of the intervals between consecutive crucial events, which have an inverse power law (IPL) PDF $\psi (\tau )\propto 1/(\tau )^{\mu }$ with $1<\mu <3$. We show that the action of crucial events has an effect compatible with the shared notion of complexity-induced entropy reduction, while making the synchronization between systems sharing the same complexity different from chaos synchronization. We establish the coupling between two temporally complex systems using a phenomenological approach inspired by models of swarm cognition and prove that complexity matching, namely sharing the same IPL index $\mu $, facilitates the transport of information, generating perfect synchronization. This new form of complexity matching is expected to contribute significantly to progress in understanding and improving biofeedback therapies.


title: Packing and covering directed triangles

abstract: We prove that if a directed multigraph $D$ has at most $t$ pairwise arc-disjoint directed triangles, then there exists a set of less than $2t$ arcs in $D$ which meets all directed triangles in $D$, except in the trivial case $t=0$. This answers affirmatively a question of Tuza from 1990.


title: Self-Driving Vehicle Verification Towards a Benchmark

abstract: Industrial cyber-physical systems are hybrid systems with strict safety requirements. Despite not having a formal semantics, most of these systems are modeled using Stateflow/Simulink for mainly two reasons: (1) it is easier to model, test, and simulate using these tools, and (2) dynamics of these systems are not supported by most other tools. Furthermore, with the ever growing complexity of cyber-physical systems, grows the gap between what can be modeled using an automatic formal verification tool and models of industrial cyber-physical systems. In this paper, we present a simple formal model for self-deriving cars. While after some simplification, safety of this system has already been proven manually, to the best of our knowledge, no automatic formal verification tool supports its dynamics. We hope this serves as a challenge problem for formal verification tools targeting industrial applications.


title: Structure property relationship in $(TiZrNbCu)_{1-x}Ni_x$ metallic glasses

abstract: The atomic structure, electronic structure and physical properties of $(TiZrNbCu)_{(1-x)}Ni_x$ (x $\leq$ 0.5) metallic glasses (MG) were studied in both the high-entropy (0<x<0.35) and the higher Ni concentration range (x $\geq$ 0.35). Atomic structure studies performed with X-ray diffraction and synchrotron powder diffraction provided average atomic volumes, structure factors, radial distribution functions, coordination numbers and packing densities. Electronic structure studies performed using photoemission spectroscopy and low-temperature specific heat provided information about the electronic density of states within the valence band and at the Fermi level and also about interatomic bonding and atomic vibrations (from the Debye temperature and the boson peak). Variations of both atomic structure and electronic structure with x showed a clear change for x $\geq$ 0.35, which corresponds to a valence electron number $\geq$ 7.4. All physical properties, namely thermal stability parameters, Debye temperatures, boson peaks, magnetic, elastic and electronic transport properties change their concentration-dependence for x$\geq$ 0.35. The results are compared with those for binary and ternary MGs of the same elements.


title: Performance of stochastic clocks in the Alternate Ticks Game

abstract: In standard quantum theory, time is not an observable. It enters as a parameter in the Schr\"odinger equation, but there is no measurement operator associated to it. Nevertheless, one may take an operational viewpoint and regard time as the information one can read from clocks. The Alternate Ticks Game, introduced in arXiv:1506.01373, is a completely operational means to quantify the accuracy of time scales generated by clocks. The idea is to count the number of ticks that two copies of a clock can produce until they run out of synchronisation. Here we investigate the performance of stochastic clocks in this game. These are clocks which are classical in the sense that they do not exploit quantum coherence. Our results support earlier conjectures that their accuracy grows linearly in the size of the clockwork, measured in terms of the dimension of the associated Hilbert space. In particular, we derive explicit bounds on the accuracy of a natural class of stochastic clocks, the stochastic ladder clocks.


title: An Annealed Sequential Monte Carlo Method for Bayesian Phylogenetics

abstract: We describe an "embarrassingly parallel" method for Bayesian phylogenetic inference, annealed Sequential Monte Carlo, based on recent advances in the Sequential Monte Carlo literature such as adaptive determination of annealing parameters. The algorithm provides an approximate posterior distribution over trees and evolutionary parameters as well as an unbiased estimator for the marginal likelihood. This unbiasedness property can be used for the purpose of testing the correctness of posterior simulation software. We evaluate the performance of phylogenetic annealed Sequential Monte Carlo by reviewing and comparing with other computational Bayesian phylogenetic methods, in particular, different marginal likelihood estimation methods. Unlike previous Sequential Monte Carlo methods in phylogenetics, our annealed method can utilize standard Markov chain Monte Carlo tree moves and hence benefit from the large inventory of such moves available in the literature. Consequently, the annealed Sequential Monte Carlo method should be relatively easy to incorporate into existing phylogenetic software packages based on Markov chain Monte Carlo algorithms. We illustrate our method using simulation studies and real data analysis.


title: Augmented Reality-based Feedback for Technician-in-the-loop C-arm Repositioning

abstract: Interventional C-arm imaging is crucial to percutaneous orthopedic procedures as it enables the surgeon to monitor the progress of surgery on the anatomy level. Minimally invasive interventions require repeated acquisition of X-ray images from different anatomical views to verify tool placement. Achieving and reproducing these views often comes at the cost of increased surgical time and radiation dose to both patient and staff. This work proposes a marker-free "technician-in-the-loop" Augmented Reality (AR) solution for C-arm repositioning. The X-ray technician operating the C-arm interventionally is equipped with a head-mounted display capable of recording desired C-arm poses in 3D via an integrated infrared sensor. For C-arm repositioning to a particular target view, the recorded C-arm pose is restored as a virtual object and visualized in an AR environment, serving as a perceptual reference for the technician. We conduct experiments in a setting simulating orthopedic trauma surgery. Our proof-of-principle findings indicate that the proposed system can decrease the 2.76 X-ray images required per desired view down to zero, suggesting substantial reductions of radiation dose during C-arm repositioning. The proposed AR solution is a first step towards facilitating communication between the surgeon and the surgical staff, improving the quality of surgical image acquisition, and enabling context-aware guidance for surgery rooms of the future. The concept of technician-in-the-loop design will become relevant to various interventions considering the expected advancements of sensing and wearable computing in the near future.


title: Physics-Inspired Optimization for Quadratic Unconstrained Problems Using a Digital Annealer

abstract: The Fujitsu Digital Annealer (DA) is designed to solve fully connected quadratic unconstrained binary optimization (QUBO) problems. It is implemented on application-specific CMOS hardware and currently solves problems of up to 1024 variables. The DA's algorithm is currently based on simulated annealing; however, it differs from it in its utilization of an efficient parallel-trial scheme and a dynamic escape mechanism. In addition, the DA exploits the massive parallelization that custom application-specific CMOS hardware allows. We compare the performance of the DA to simulated annealing and parallel tempering with isoenergetic cluster moves on two-dimensional and fully connected spin-glass problems with bimodal and Gaussian couplings. These represent the respective limits of sparse versus dense problems, as well as high-degeneracy versus low-degeneracy problems. Our results show that the DA currently exhibits a time-to-solution speedup of roughly two orders of magnitude for fully connected spin-glass problems with bimodal or Gaussian couplings, over the single-core implementations of simulated annealing and parallel tempering Monte Carlo used in this study. The DA does not appear to exhibit a speedup for sparse two-dimensional spin-glass problems, which we explain on theoretical grounds. We also benchmarked an early implementation of the Parallel Tempering DA. Our results suggest an improved scaling over the other algorithms for fully connected problems of average difficulty with bimodal disorder. The next generation of the DA is expected to be able to solve fully connected problems up to 8192 variables in size. This would enable the study of fundamental physics problems and industrial applications that were previously inaccessible using standard computing hardware or special-purpose quantum annealing machines.


title: Existence of infinitely many minimal hypersurfaces in closed manifolds

abstract: Using min-max theory, we show that in any closed Riemannian manifold of dimension at least 3 and at most 7, there exist infinitely many smoothly embedded closed minimal hypersurfaces. It proves a conjecture of S.-T. Yau. This paper builds on the methods developed by F. C. Marques and A. Neves.


title: Aggregation-Based Certificate Transparency Gossip

abstract: Certificate Transparency (CT) requires that every CA-issued TLS certificate must be publicly logged. While a CT log need not be trusted in theory, it relies on the assumption that every client observes and cryptographically verifies the same log. As such, some form of gossip mechanism is needed in practice. Despite CT being adopted by several major browser vendors, no gossip mechanism is widely deployed. We suggest an aggregation-based gossip mechanism that passively observes cryptographic material that CT logs emit in plaintext, aggregating at packet processors (such as routers and switches) to periodically verify log consistency off-path. In other words, gossip is provided as-a-service by the network. Based on 20 days of RIPE Atlas measurements that represent clients from 3500 autonomous systems and 40% of the IPv4 space, our proposal can be deployed incrementally for a realistic threat model with significant protection against split-viewing CT logs. We also show that aggregation-based gossip can be implemented for a variety of packet processors using P4 and XDP, running at 10 Gbps line-speed.


title: Robust Resource Allocation for MISO Cognitive Radio Networks Under Two Practical Non-Linear Energy Harvesting Models

abstract: To achieve a good tradeoff between the consumed power and the harvested power, a robust resource optimization problem is studied in a multiple-input single-output cognitive radio network with simultaneous wireless information and power transfer under imperfect channel state information. Unlike most of the existing works that assume an ideal linear energy harvesting (EH) model, we assume two practical non-linear EH models. In order to solve the resulting challenging non-convex problem, we propose an algorithm based on successive convex approximation and prove that a rank-one solution is obtained in each iteration of the proposed algorithm. Our simulation results quantify the effect of the sensitivity threshold of the EH circuit on the harvested power.


title: Forecasting Internally Displaced Population Migration Patterns in Syria and Yemen

abstract: Armed conflict has led to an unprecedented number of internally displaced persons (IDPs) - individuals who are forced out of their homes but remain within their country. IDPs often urgently require shelter, food, and healthcare, yet prediction of when large fluxes of IDPs will cross into an area remains a major challenge for aid delivery organizations. Accurate forecasting of IDP migration would empower humanitarian aid groups to more effectively allocate resources during conflicts. We show that monthly flow of IDPs from province to province in both Syria and Yemen can be accurately forecasted one month in advance, using publicly available data. We model monthly IDP flow using data on food price, fuel price, wage, geospatial, and news data. We find that machine learning approaches can more accurately forecast migration trends than baseline persistence models. Our findings thus potentially enable proactive aid allocation for IDPs in anticipation of forecasted arrivals.


title: Warped product submanifolds in metallic Riemannian manifolds

abstract: In this paper, we study the existence of proper warped product submanifolds in metallic (or Golden) Riemannian manifolds and we discuss about semi-invariant, semi-slant and, respectively, hemi-slant warped product submanifolds in metallic and Golden Riemannian manifolds. Also, we provide some examples of warped product submanifolds in Euclidean spaces.


title: Homogenization of obstacle problems in Orlicz-Sobolev spaces

abstract: We study the homogenization of obstacle problems in Orlicz-Sobolev spaces for a wide class of monotone operators (possibly degenerate or singular) of the $p(\cdot)$-Laplacian type. Our approach is based on the Lewy-Stampacchia inequalities, which then give access to a compactness argument. We also prove the convergence of the coincidence sets under non-degeneracy conditions.


title: Cooling+Heating Flows in Galaxy clusters: Turbulent heating, spectral modelling, and cooling efficiency

abstract: The discrepancy between expected and observed cooling rates of X-ray emitting gas has led to the {\it cooling flow problem} at the cores of clusters of galaxies. A variety of models have been proposed to model the observed X-ray spectra and resolve the cooling flow problem, which involves heating the cold gas through different mechanisms. As a result, realistic models of X-ray spectra of galaxy clusters need to involve both heating {\it and} cooling mechanisms. In this paper, we argue that the heating time-scale is set by the magnetohydrodynamic (MHD) turbulent viscous heating for the Intracluster plasma, parametrised by the Shakura-Sunyaev viscosity parameter, $\alpha$. Using a cooling+heating flow model, we show that a value of $\alpha\simeq 0.05$ (with 10\% scatter) provides improved fits to the X-ray spectra of cooling flow, while at the same time, predicting reasonable cooling efficiency, $\epsilon_{cool} = 0.33^{+0.63}_{-0.15}$. Our inferred values for $\alpha$ based on X-ray spectra are also in line with direct measurements of turbulent pressure in simulations and observations of galaxy clusters. This simple picture unifies astrophysical accretion, as a balance of MHD turbulent heating and cooling, across more than 16 orders of magnitudes in scale, from neutron stars to galaxy clusters.


title: Null boundary terms for Lanczos-Lovelock gravity

abstract: We derive boundary terms appropriate for the general Lanczos-Lovelock action on a null boundary, when Dirichlet boundary conditions are imposed. We believe that these boundary terms have been derived for the first time in the literature. In this derivation, we rely only on the structure of the boundary variation of the action for Lanczos-Lovelock gravity. We also provide the null boundary term for Gauss-Bonnet gravity separately.


title: On Banach Structure of Multivariate BV Spaces I

abstract: We introduce and study multivariate generalizations of the classical BV spaces of Jordan, F. Riesz and Wiener. The family of the introduced spaces contains or is intimately related to a considerable class of function spaces of modern analysis including BMO, BV, Morrey spaces and those of Sobolev of arbitrary smoothness, Besov and Triebel-Lizorkin spaces. We prove under mild restrictions that the BV spaces of this family are dual and present constructive characterizations of their preduals via atomic decompositions. Moreover, we show that under additional restrictions such a predual space is isometrically isomorphic to the dual space of the separable subspace of the related BV space generated by $C^\infty$ functions. As a corollary we obtain the "two stars theorem" asserting that the second dual of this separable subspace is isometrically isomorphic to the BV space. An essential role in the proofs play approximation properties of the BV spaces under consideration, in particular, weak$^*$ denseness of their subspaces of $C^\infty$ functions. Our results imply the similar ones (old and new) for the classical function spaces listed above obtained by the unified approach.


title: Staircase-PIR: Universally Robust Private Information Retrieval

abstract: We consider the problem of designing private information retrieval (PIR) schemes on data of $m$ files replicated on $n$ servers that can possibly collude. We focus on devising robust PIR schemes that can tolerate stragglers, i.e., slow or unresponsive servers. In many settings, the number of stragglers is not known a priori or may change with time. We define universally robust PIR as schemes that achieve PIR capacity asymptotically in $m$ and simultaneously for any number of stragglers up to a given threshold. We introduce Staircase-PIR schemes and prove that they are universally robust. Towards that end, we establish an equivalence between robust PIR and communication efficient secret sharing.


title: Measuring the Variation in Nuclear Charge Radius of Xe Isotopes by EUV Spectroscopy of Highly-Charged Na-like Ions

abstract: The variation in mean-square nuclear charge radius of xenon isotopes was measured utilizing a new method based on extreme ultraviolet spectroscopy of highly charged Na-like ions. The isotope shift of the Na-like D1 (3s $^{2}$S$_{1/2}$ - 3p $^2$P$_{1/2}$) transition between the $^{124}$Xe and $^{136}$Xe isotopes was experimentally determined using the electron beam ion trap facility at the National Institute of Standards and Technology. The mass shift and the field shift coefficients were calculated with enhanced precision by relativistic many-body perturbation theory and multi-configuration Dirac-Hartree-Fock method. The mean-square nuclear charge radius difference was found to be $\delta<r^2>^{136, 124}$ = 0.269(0.042) fm$^2$. Our result has smaller uncertainty than previous experimental results and agrees with the recommended value by Angeli and Marinova [I. Angeli and K. P. Marinova, At. Data and Nucl. Data Tables {\bf 99}, 69-95 (2013)].


title: The Relation Between Classical and Quantum Mechanics

abstract: Familiar formulations of classical and quantum mechanics are shown to follow from a general theory of mechanics based on pure states with an intrinsic probability structure. This theory is developed to the stage where theorems from quantum logic enable expression of the state geometry in Hilbert space. Quantum and classical mechanics are then elaborated and applied to subsystems and the measurement process. Consideration is also given to space-time geometry and the constraints this places on the dynamics. Physics and Mathematics, it is argued, are growing apart; the inadequate treatment of approximations in general and localisation in quantum mechanics in particular are seen as contributing factors. In the description of systems, the link between localisation and lack of knowledge shows that quantum mechanics should reflect the domain of applicability. Restricting the class of states provides a means of achieving this goal. Localisation is then shown to have a mathematical expression in terms of compactness, which in turn is applied to yield a topological theory of bound and scattering states. Finally, the thesis questions the validity of classical limits and quantisations in intertheoretic reduction, and demonstrates that a widely accepted classical limit does not constitute a proof of reduction. It proposes a procedure for determining whether classical and quantum mechanics are weakly equivalent over a domain of applicability, and concludes that, in this restricted sense, classical mechanics reduces to quantum mechanics.

categories: quant-ph math-ph math.MP physics.hist-ph


title: Astrophysical factors of ${}^{12}{\rm C}+{}^{12}{\rm C}$ fusion from Trojan horse method

abstract: Carbon-carbon burning plays an important role in many stellar environments. Recently, using the indirect Trojan horse method A. Tumino {\em et al.} reported [Nature {\bf 557} 687 (2018)] a strong rise of the astrophysical factor for the carbon-carbon fusion at low resonance energies. In this paper, we demonstrate that this rise is the artifact of using an invalid plane-wave approximation. It is shown that the calculated renormalization factor decreases the astrophysical factor from [A. Tumino {\em et al.}, Nature {\bf 557} 687 (2018)] at the resonance energies of $E=0.8-0.9$ MeV by as much as $\approx 10^{3}$ times.


title: Diffusion Scattering Transforms on Graphs

abstract: Stability is a key aspect of data analysis. In many applications, the natural notion of stability is geometric, as illustrated for example in computer vision. Scattering transforms construct deep convolutional representations which are certified stable to input deformations. This stability to deformations can be interpreted as stability with respect to changes in the metric structure of the domain. In this work, we show that scattering transforms can be generalized to non-Euclidean domains using diffusion wavelets, while preserving a notion of stability with respect to metric changes in the domain, measured with diffusion maps. The resulting representation is stable to metric perturbations of the domain while being able to capture "high-frequency" information, akin to the Euclidean Scattering.


title: Differential Geometry of Rotation Minimizing Frames, Spherical Curves, and Quantum Mechanics of a Constrained Particle

abstract: This thesis is devoted to the Differential Geometry of curves and surfaces along with applications in Quantum Mechanics. In the 1st part we introduce the well known Frenet frame. Later, we show that the curvature function is a lower bound for the scalar angular velocity of any other orthonormal moving frame, from which one defines Rotation Minimizing (RM) frames as the ones that achieve this minimum. Remarkably, RM frames are ideal to study spherical curves and allow us to characterize them through a simple linear equation. We also apply these ideas to curves that lie on level surfaces, by reinterpreting the problem in the context of a metric induced by a Hessian, which may fail to be positive or non-degenerate and naturally leads us to a Lorentz-Minkowski $\mathbb{E}_1^3$ or isotropic $\mathbb{I}^3$ space. Here we develop a systematic approach to construct RM frames and characterize spherical curves in $\mathbb{E}_1^3$ and $\mathbb{I}^3$ and furnish a criterion for a curve to lie on a level surface. Finally, we extend these tools to characterize geodesic spherical curves in hyperbolic and spherical spaces. In the 2nd part we apply the previous ideas in the quantum dynamics of a constrained particle. After describing the confining potential formalism, from which emerges a geometry-induced potential (GIP), we devote our attention to tubular surfaces to model curved nanotubes. The use of RM frames offers a simpler description for the constrained dynamics and allows us to show that the torsion of the centerline of a tube gives rise to a geometric phase. Later, we study the problem of prescribed GIP for surfaces. Here we explore the GIP for surfaces invariant by a 1-parameter group of isometries, i.e., cylindrical, revolution, and helicoidal surfaces. Finally, we devote a special attention to helicoidal minimal surfaces and prove the existence of geometry-induced bound and localized states.


title: Efficient Graph Compression Using Huffman Coding Based Techniques

abstract: Graphs have been extensively used to represent data from various domains. In the era of Big Data, information is being generated at a fast pace, and analyzing the same is a challenge. Various methods have been proposed to speed up the analysis of the data and also mining it for information. All of this often involves using a massive array of compute nodes, and transmitting the data over the network. Of course, with the huge quantity of data, this poses a major issue to the task of gathering intelligence from data. Therefore, in order to address such issues with Big Data, using data compression techniques is a viable option. Since graphs represent most real world data, methods to compress graphs have been in the forefront of such endeavors. In this paper we propose techniques to compress graphs by finding specific patterns and replacing those with identifiers that are of variable length, an idea inspired by Huffman Coding. Specifically, given a graph G = (V, E), where V is the set of vertices and E is the set of edges, and |V| = n, we propose methods to reduce the space requirements of the graph by compressing the adjacency representation of the same. The proposed methods show up to 80% reduction is the space required to store the graphs as compared to using the adjacency matrix. The methods can also be applied to other representations as well. The proposed techniques help solve the issues related to computing on the graphs on resources limited compute nodes, as well as reduce the latency for transfer of data over the network in case of distributed computing.


title: Higgs boson production at large transverse momentum within SMEFT: analytical results

abstract: We consider Higgs boson production through gluon fusion at large transverse momentum in hadronic collisions. We present the analytic expressions of the relevant one-loop QCD amplitudes including the effects of the complete set of dimension-six operators. The latter correspond to modifications of the top and bottom Yukawa couplings, to an effective point-like Higgs coupling to gluons and to the chromomagnetic operator of the top quark. The quantitative impact of the chromomagnetic operator is also studied. Our results confirm previous findings that the effect of the chromomagnetic operator at high $p_T$ can be large and should not be neglected.


title: Coherency-broken Bragg filters: surpassing on-chip rejection limitations

abstract: Selective on-chip optical filters with high rejection levels are key components for a wide range of advanced photonic circuits. However, maximum achievable rejection in state-of-the-art on-chip devices is seriously limited by phase errors arising from fabrication imperfections. Due to coherent interactions, unwanted phase-shifts result in detrimental destructive interferences that distort the filter response, whatever the chosen strategy (resonators, interferometers, Bragg filters, etc.). Here we propose and experimentally demonstrate a radically different approach to overcome this fundamental limitation, based on coherency-broken Bragg filters. We exploit non-coherent interaction among modal-engineered waveguide Bragg gratings separated by single-mode waveguides to yield effective cascading, even in the presence of fabrication errors. This technologically independent approach allows seamless combination of filter stages with moderate performance, providing a dramatic increase of on-chip rejection. Based on this concept, we experimentally demonstrate on-chip non-coherent cascading of Si Bragg filters with a record light rejection exceeding 80 dB in the C-band.


title: Smart Inverter Grid Probing for Learning Loads: Part I - Identifiability Analysis

abstract: Distribution grids currently lack comprehensive real-time metering. Nevertheless, grid operators require precise knowledge of loads and renewable generation to accomplish any feeder optimization task. At the same time, new grid technologies, such as solar photovoltaics and energy storage units are interfaced via inverters with advanced sensing and actuation capabilities. In this context, this two-part work puts forth the idea of engaging power electronics to probe an electric grid and record its voltage response at actuated and metered buses, to infer non-metered loads. Probing can be accomplished by commanding inverters to momentarily perturb their power injections. Multiple probing actions can be induced within a few tens of seconds. In Part I, load inference via grid probing is formulated as an implicit nonlinear system identification task, which is shown to be topologically observable under certain conditions. The conditions can be readily checked upon solving a max-flow problem on a bipartite graph derived from the feeder topology and the placement of actuated and non-metered buses. The analysis holds for single- and multi-phase grids, radial or meshed, and applies to phasor or magnitude-only voltage data. The topological observability of distribution systems using smart meter or phasor data is cast and analyzed a special case.


title: Domain Adaptation for Infection Prediction from Symptoms Based on Data from Different Study Designs and Contexts

abstract: Acute respiratory infections have epidemic and pandemic potential and thus are being studied worldwide, albeit in many different contexts and study formats. Predicting infection from symptom data is critical, though using symptom data from varied studies in aggregate is challenging because the data is collected in different ways. Accordingly, different symptom profiles could be more predictive in certain studies, or even symptoms of the same name could have different meanings in different contexts. We assess state-of-the-art transfer learning methods for improving prediction of infection from symptom data in multiple types of health care data ranging from clinical, to home-visit as well as crowdsourced studies. We show interesting characteristics regarding six different study types and their feature domains. Further, we demonstrate that it is possible to use data collected from one study to predict infection in another, at close to or better than using a single dataset for prediction on itself. We also investigate in which conditions specific transfer learning and domain adaptation methods may perform better on symptom data. This work has the potential for broad applicability as we show how it is possible to transfer learning from one public health study design to another, and data collected from one study may be used for prediction of labels for another, even collected through different study designs, populations and contexts.


title: Smart Inverter Grid Probing for Learning Loads: Part II - Probing Injection Design

abstract: This two-part work puts forth the idea of engaging power electronics to probe an electric grid to infer non-metered loads. Probing can be accomplished by commanding inverters to perturb their power injections and record the induced voltage response. Once a probing setup is deemed topologically observable by the tests of Part I, Part II provides a methodology for designing probing injections abiding by inverter and network constraints to improve load estimates. The task is challenging since system estimates depend on both probing injections and unknown loads in an implicit nonlinear fashion. The methodology first constructs a library of candidate probing vectors by sampling over the feasible set of inverter injections. Leveraging a linearized grid model and a robust approach, the candidate probing vectors violating voltage constraints for any anticipated load value are subsequently rejected. Among the qualified candidates, the design finally identifies the probing vectors yielding the most diverse system states. The probing task under noisy phasor and non-phasor data is tackled using a semidefinite-program (SDP) relaxation. Numerical tests using synthetic and real-world data on a benchmark feeder validate the conditions of Part I; the SDP-based solver; the importance of probing design; and the effects of probing duration and noise.


title: Rearrangement and Prekopa-Leindler type inequalities

abstract: We investigate the interactions of functional rearrangements with Prekopa-Leindler type inequalities. It is shown that that a general class of integral inequalities tighten on rearrangement to "isoperimetric" sets with respect to a relevant measure. Applications to the Borell-Brascamp-Lieb, Borell-Ehrhart, and the recent polar Prekopa-Leindler inequalities are demonstrated. It is also proven that an integrated form of the Gaussian log-Sobolev inequality decreases on half-space rearrangement.


title: Bayesian Optimization of Combinatorial Structures

abstract: The optimization of expensive-to-evaluate black-box functions over combinatorial structures is an ubiquitous task in machine learning, engineering and the natural sciences. The combinatorial explosion of the search space and costly evaluations pose challenges for current techniques in discrete optimization and machine learning, and critically require new algorithmic ideas. This article proposes, to the best of our knowledge, the first algorithm to overcome these challenges, based on an adaptive, scalable model that identifies useful combinatorial structure even when data is scarce. Our acquisition function pioneers the use of semidefinite programming to achieve efficiency and scalability. Experimental evaluations demonstrate that this algorithm consistently outperforms other methods from combinatorial and Bayesian optimization.


title: Quantitative principles for precise engineering of sensitivity in carbon-based electrochemical sensors

abstract: A major practical barrier for implementing carbon-based electrode arrays with high device-packing density is to ensure large, predictable, and homogeneous sensitivities across the array. Overcoming this barrier depends on quantitative models to predict electrode sensitivity from its material structure. However, such models are currently lacking. Here, we show that the sensitivity of multilayer graphene electrodes increases linearly with the average point defect density, whereas it is unaffected by line defects or oxygen-containing groups. These quantitative relationships persist until the electrode material transitions to a fully disordered sp2 carbon, where sensitivity declines sharply. We show that our results generalize to a variety of graphene production methods and use them to derive a predictive model that guides nano-engineering graphene structure for optimum sensitivity. Our approach achieves reproducible fabrication of miniaturized sensors with extraordinarily higher sensitivity than conventional material. These results lay the foundation for new integrated electrochemical sensor arrays based on nano-engineered graphene.


title: Deep SNP: An End-to-end Deep Neural Network with Attention-based Localization for Break-point Detection in SNP Array Genomic data

abstract: Diagnosis and risk stratification of cancer and many other diseases require the detection of genomic breakpoints as a prerequisite of calling copy number alterations (CNA). This, however, is still challenging and requires time-consuming manual curation. As deep-learning methods outperformed classical state-of-the-art algorithms in various domains and have also been successfully applied to life science problems including medicine and biology, we here propose Deep SNP, a novel Deep Neural Network to learn from genomic data. Specifically, we used a manually curated dataset from 12 genomic single nucleotide polymorphism array (SNPa) profiles as truth-set and aimed at predicting the presence or absence of genomic breakpoints, an indicator of structural chromosomal variations, in windows of 40,000 probes. We compare our results with well-known neural network models as well as Rawcopy though this tool is designed to predict breakpoints and in addition genomic segments with high sensitivity. We show, that Deep SNP is capable of successfully predicting the presence or absence of a breakpoint in large genomic windows and outperforms state-of-the-art neural network models. Qualitative examples suggest that integration of a localization unit may enable breakpoint detection and prediction of genomic segments, even if the breakpoint coordinates were not provided for network training. These results warrant further evaluation of DeepSNP for breakpoint localization and subsequent calling of genomic segments.


title: Reaction dynamics of cyanohydrins with hydrosulfide in water

abstract: We studied the reaction dynamics of a proposed prebiotic reaction theoretically. The chemical process involves acetone cyanohydrin or formalcyanohydrin reacting with hydrosulfide in an aqueous environment. Rate constants and populations of reactant and product bimolecular geometric orientations for the reactions were obtained by using density functional theory for the energies, transition state theory for the rates and matrix exponentiation as well as the hybrid tau-leaping algorithm for the population dynamics. The role of including the solvent explicitly versus implicitly was also investigated. We found that adding water or hydrogen sulfide molecules explicitly lowers the activation energy barrier and leads to a more efficient reaction pathway. In particular, hydrogen sulfide was a better proton donor. Further, when adding the solvent explicitly a synchronous reaction mechanism was observed while a concerted reaction mechanism was seen when using an implicit solvent model. Finally, we studied the role of including more than one reactant and product bimolecular orientation geometry in the dynamics. Including all pathways, reactant to reactant, product to product, reactant to product and product to reactant reduced the reaction half-time by two orders of magnitude and was therefore determined to be fundamental for an accurate description of the dynamics. Overall, we found that most reactions which involve formalcyanohydrin occur more rapidly or at the same speed as reactions which involve acetone cyanohydrin at room temperature.


title: The Variable Wavelength Dependence of the Dipping event of KIC 8462852

abstract: First observed with the Kepler mission, KIC 8462852 undergoes unexplained dimming events, "dips," on the timescale of days which were again observed from the ground from May to December 2017. Monitored with multi-band photometry by the Los Cumbres Observatory, all four dips of the "Elsie dip family" display clear wavelength dependence. We measure how the wavelength dependence changes over the whole dimming event, including the dimming between the dips and the brightening event (the `blip') which occurs after the dips. We find that a single wavelength dependence does not fit the entire light curve and the dimming occurring between the dips is non-gray and varies in time. Because of the non-gray dimming between the dips, we measure the wavelength dependence of the dips separately and without the extra depth from this dimming. Such measurements yield a different estimate of the wavelength dependence the wavelength dependence of the dips but remains consistent with the previous measurement except for Elsie (the first dip), which is surrounded by dimming with strong wavelength dependence. We find the range of the wavelength dependence variation of the entire 2017 light curve is consistent with optically-thin dust with an average radius of $r<1\,\mu$m and the dust causing just the dips being $r<0.5\,\mu$m. Since the dependence is time-dependent, the dust occulting the star must be heterogeneous in size, composition, or both and the distributions of these properties along the line of sight must change over time.


title: The Meeting Time of Multiple Random Walks

abstract: This article rigorously analyzes the meeting time between pursuers and evaders performing random walks on digraphs. There exist several bounds on the expected meeting time between random walkers on graphs in the literature, however, closed-form expressions are limited in scope. By utilizing the notion that multiple random walks on a common graph can be understood as a single random walk on the Kronecker product graph, we are able to provide the first analytic expression for the meeting time in terms of the transition matrices of the random walkers when modeled by either discrete-time Markov chains or continuous-time Markov chains. We further extend the results to the case of multiple pursuers and multiple evaders performing independent random walks. We present various sufficient conditions for pairs (or tuples) of transition matrices that satisfy certain conditions on the absorbing classes for which finite meeting times are guaranteed to exist.


title: A linear state feedback switching rule for global stabilization of switched nonlinear systems about a nonequilibrium point

abstract: A switched equilibrium of a switched system of two subsystems is a such a point where the vector fields of the two subsystems point strictly towards one another. Using the concept of stable convex combination that was developed by Wicks-Peleties-DeCarlo (1998) for linear systems, Bolzern-Spinelli (2004) offered a design of a state feedback switching rule that is capable to stabilize an affine switched system to any switched equilibrium. The state feedback switching rule of Bolzern-Spinelli gives a nonlinear (quadratic) switching threshold passing through the switched equilibrium. In this paper we prove that the switching threshold (i.e. the associated switching rule) can be chosen linear, if each of the subsystems of the switched system under consideration are stable.