title: Koopman Performance Analysis of Nonlinear Consensus Networks

abstract: Spectral decomposition of dynamical systems is a popular methodology to investigate the fundamental qualitative and quantitative properties of these systems and their solutions. In this chapter, we consider a class of nonlinear cooperative protocols, which consist of multiple agents that are coupled together via an undirected state-dependent graph. We develop a representation of the system solution by decomposing the nonlinear system utilizing ideas from the Koopman operator theory and its spectral analysis. We use recent results on the extensions of the well-known Hartman theorem for hyperbolic systems to establish a connection between the original nonlinear dynamics and the linearized dynamics in terms of Koopman spectral properties. The expected value of the output energy of the nonlinear protocol, which is related to the notions of coherence and robustness in dynamical networks, is evaluated and characterized in terms of Koopman eigenvalues, eigenfunctions, and modes. Spectral representation of the performance measure enables us to develop algorithmic methods to assess the performance of this class of nonlinear dynamical networks as a function of their graph topology. Finally, we propose a scalable computational method for approximation of the components of the Koopman mode decomposition, which is necessary to evaluate the systemic performance measure of the nonlinear dynamic network.


title: Spectra of Hadamard matrices

abstract: A Butson Hadamard matrix $H$ has entries in the kth roots of unity, and satisfies the matrix equation $HH^{\ast} = nI_{n}$. We write $\mathrm{BH}(n, k)$ for the set of such matrices. A complete morphism of Butson matrices is a map $\mathrm{BH}(n, k) \rightarrow \mathrm{BH}(m, \ell)$. In this paper, we develop a technique for controlling the spectra of certain Hadamard matrices. For each integer $t$, we construct a real Hadamard matrix $H_{t}$ of order $n_{t} = 2^{2^{t-1}-1}$ such that the minimal polynomial of $\frac{1}{\sqrt{n_{t}}}H_{t}$ is the cyclotomic polynomial $\Phi_{2^{t+1}}(x)$. Such matrices yield new examples of complete morphisms \[ \mathrm{BH}(n, 2^{t}) \rightarrow \mathrm{BH}(2^{2^{t-1}-1}n, 2)\,, \] for each $t \geq 2$, generalising a well-known result of Turyn.


title: Morse Code Datasets for Machine Learning

abstract: We present an algorithm to generate synthetic datasets of tunable difficulty on classification of Morse code symbols for supervised machine learning problems, in particular, neural networks. The datasets are spatially one-dimensional and have a small number of input features, leading to high density of input information content. This makes them particularly challenging when implementing network complexity reduction methods. We explore how network performance is affected by deliberately adding various forms of noise and expanding the feature set and dataset size. Finally, we establish several metrics to indicate the difficulty of a dataset, and evaluate their merits. The algorithm and datasets are open-source.


title: Collider signatures for dark matter and long-lived particles with Pythia 8

abstract: We describe the implementation of four models for production of dark matter or as- sociated particles at the LHC based on the simplest extensions of the Standard Model. The first kind of models include dark matter production via s-channel mediators. This includes production in association with a jet for a vector boson (Z') or scalar (S) mediator as well as mono-higgs production via associated hZ' production. We implement the simplest t-channel mediator in the form of a scalar with leptonic quantum numbers and completely generic Yukawa couplings to the dark matter fermion and a right-handed SM lepton. Finally, we implement a generalised model of mixed dark matter where the dark matter is a mixture of an SU(2) singlet and N-plet. We find that the last two models are also ideally suited to study the production of a range of long-lived particle signatures. We illustrate this by reinterpreting the limits from CMS search for heavy charged particles for the t-channel case.


title: DeepMove: Learning Place Representations through Large Scale Movement Data

abstract: Understanding and reasoning about places and their relationships are critical for many applications. Places are traditionally curated by a small group of people as place gazetteers and are represented by an ID with spatial extent, category, and other descriptions. However, a place context is described to a large extent by movements made from/to other places. Places are linked and related to each other by these movements. This important context is missing from the traditional representation. We present DeepMove, a novel approach for learning latent representations of places. DeepMove advances the current deep learning based place representations by directly model movements between places. We demonstrate DeepMove's latent representations on place categorization and clustering tasks on large place and movement datasets with respect to important parameters. Our results show that DeepMove outperforms state-of-the-art baselines. DeepMove's representations can provide up to 15% higher than competing methods in matching rate of place category and result in up to 39% higher silhouette coefficient value for place clusters. DeepMove is spatial and temporal context aware. It is scalable. It outperforms competing models using much smaller training dataset (a month or 1/12 of data). These qualities make it suitable for a broad class of real-world applications.


title: Multiparameter perturbation theory of matrices and linear operators

abstract: We show that a normal matrix $A$ with coefficient in $\mathbb C[[X]]$, $X=(X_1, \ldots, X_n)$, can be diagonalized, provided the discriminant $\Delta_A $ of its characteristic polynomial is a monomial times a unit. The proof is an adaptation of the algorithm of proof of Abhyankar-Jung Theorem. As a corollary we obtain the singular value decomposition for an arbitrary matrix $A$ with coefficient in $\mathbb C[[X]]$ under a similar assumption on $\Delta_{AA^*} $ and $\Delta_{A^*A} $. We also show real versions of these results, i.e. for coefficients in $\mathbb R[[X]]$, and deduce several results on multiparameter perturbation theory for normal matrices with real analytic, quasi-analytic, or Nash coefficients.


title: Transforming a large introductory lab course: impacts on views about experimental physics

abstract: Laboratory courses are key components of most undergraduate physics programs. Lab courses often aim to achieve the following learning outcomes: developing students' experimental skills, engaging students in authentic scientific practices, reinforcing concepts, and inspiring students' interest and engagement in physics. Some of these outcomes can be measured by the Colorado Learning Attitudes about Science Survey for Experimental Physics (E-CLASS), a research-based assessment that measures students' views about experimental physics. We used E-CLASS at the University of Colorado Boulder to measure learning outcomes during a course transformation process in which views about experimental physics were reflected in the learning goals. We collected over 600 student responses per semester from the large introductory laboratory course, both before and after implementing the course transformation. While we observed no statistically significant difference in overall post-instruction E-CLASS scores before and after the transformation, in the transformed course, student responses to three E-CLASS items that related to the goals of the transformation were more favorable than in the original course.


title: The three- and four-Higgs couplings in the general two-Higgs-doublet model

abstract: We apply the unitarity bounds and the bounded-from-below (BFB) bounds to the most general scalar potential of the two-Higgs-doublet model (2HDM). We do this in the Higgs basis, i.e. in the basis for the scalar doublets where only one doublet has vacuum expectation value. In this way we obtain bounds on the scalar masses and couplings that are valid for all 2HDMs. We compare those bounds to the analogous bounds that we have obtained for other simple extensions of the Standard Model (SM), namely the 2HDM extended by one scalar singlet and the extension of the SM through two scalar singlets.


title: Estimating Cellular Goals from High-Dimensional Biological Data

abstract: Optimization-based models have been used to predict cellular behavior for over 25 years. The constraints in these models are derived from genome annotations, measured macro-molecular composition of cells, and by measuring the cell's growth rate and metabolism in different conditions. The cellular goal (the optimization problem that the cell is trying to solve) can be challenging to derive experimentally for many organisms, including human or mammalian cells, which have complex metabolic capabilities and are not well understood. Existing approaches to learning goals from data include (a) estimating a linear objective function, or (b) estimating linear constraints that model complex biochemical reactions and constrain the cell's operation. The latter approach is important because often the known/observed biochemical reactions are not enough to explain observations, and hence there is a need to extend automatically the model complexity by learning new chemical reactions. However, this leads to nonconvex optimization problems, and existing tools cannot scale to realistically large metabolic models. Hence, constraint estimation is still used sparingly despite its benefits for modeling cell metabolism, which is important for developing novel antimicrobials against pathogens, discovering cancer drug targets, and producing value-added chemicals. Here, we develop the first approach to estimating constraint reactions from data that can scale to realistically large metabolic models. Previous tools have been used on problems having less than 75 biochemical reactions and 60 metabolites, which limits real-life-size applications. We perform extensive experiments using 75 large-scale metabolic network models for different organisms (including bacteria, yeasts, and mammals) and show that our algorithm can recover cellular constraint reactions, even when some measurements are missing.


title: Notes on pulse-width modulation appropriate for a sophomore-level electronics or programming class

abstract: Pulse-width modulation is a nice way to create pseudo analog outputs from a microcontroller (eg, an Arduino), which is limited to digital output signals. The quality of the analog output signal depends on the filter one uses to average out the digital outputs, and this filter can be analyzed with sophomore-level knowledge of resistor-capacitor circuits in charging and discharging state. These notes were used as a text supplement in a sophomore/junior level Microcontroller programming class for second and third-year Physics majors. I don't know where this material is normally treated in the Physics curriculum and figured others might find the resource useful. Comments welcome!


title: Self-Regulation in Infinite Populations with Fission-Death Dynamics

abstract: The evolution of an infinite population of interacting point entities placed in $\mathbb{R}^d$ is studied. The elementary evolutionary acts are death of an entity with rate that includes a competition term and independent fission into two entities. The population states are probability measures on the corresponding configuration space and the result is the construction of the evolution of states in the class of sub-Poissonian measures, that corresponds to the lack of clusters in such states. This is considered as a self-regulation in the population due to competition.


title: Optical Pumping of TeH+: Implications for the Search for Varying mp/me

abstract: Molecular overtone transitions provide optical frequency transitions sensitive to variation in the proton-to-electron mass ratio ($\mu\equiv m_p/m_e$). However, robust molecular state preparation presents a challenge critical for achieving high precision. Here, we characterize infrared and optical-frequency broadband laser cooling schemes for TeH$^+$, a species with multiple electronic transitions amenable to sustained laser control. Using rate equations to simulate laser cooling population dynamics, we estimate the fractional sensitivity to $\mu$ attainable using TeH$^+$. We find that laser cooling of TeH$^+$ can lead to significant improvements on current $\mu$ variation limits.


title: PSR J0636+5128: A Heated Companion in a Tight Orbit

abstract: We present an analysis of archival Gemini g', r', K and Keck H, K_s imaging of this nearby short period binary (P_B=95.8min) 2.87ms pulsar. The heated companion is clearly detected. Direct pulsar heating provides an acceptable model at the revised >700pc parallax distance. The relatively shallow light curve modulation prefers an inclination i<40deg; this high latitude view provides a likely explanation for the lack of radio signatures of wind dispersion or eclipse. It also explains the low minimum companion mass and, possibly, the faintness of the source in X- and gamma-rays.


title: Stellar chemo-kinematics of the Cetus dwarf spheroidal galaxy

abstract: In order to minimize environmental effects and gain an insight into the internal mechanisms that shape the properties of the early-type dwarf systems, we study one of the few isolated dwarf spheroidal galaxies (dSphs) of the Local Group (LG): Cetus. We obtained VLT/FORS2 spectra ($R\sim2600$) in the region of the nIR CaII triplet lines for 80 candidate red giant branch stars. The analysis yielded line-of-sight velocities and metallicities ([Fe/H]) for 54 bona fide member stars. The kinematic analysis shows that Cetus is a mainly pressure-supported ($\sigma_v = 11.0_{-1.3}^{+1.6}$ km/s), dark-matter-dominated system ($M_{1/2}/L_V = 23.9_{-8.9}^{+9.7} M_\odot/L_\odot$) with no significant signs of internal rotation. We find Cetus to be a metal-poor system with a significant [Fe/H] spread (median [Fe/H] = -1.71 dex, median-absolute-deviation = 0.49 dex), as expected for its stellar mass. We report the presence of a mild metallicity gradient compatible with those found in other dSphs of the same luminosity; we trace the presence of a stellar population gradient also in the spatial distribution of stars in different evolutionary phases in ancillary SuprimeCam photometric data. There are tentative indications of two chemo-kinematically distinct sub-populations, with the more metal-poor stars showing a hotter kinematics than the metal-richer ones. Furthermore, the photometric dataset reveals the presence of a foreground population that most likely belongs to the Sagittarius stream. This study represents a first comprehensive analysis of Cetus chemo-kinematic properties. Our results add Cetus to the growing scatter in stellar-dark matter halo properties in low-mass galactic systems. The presence of a metallicity gradient akin to those found in similar systems inhabiting different environments may hint at metallicity gradients in LG early-type dwarfs being driven by internal mechanisms.


title: A study of Schr\"oder's method for the matrix $p$th root using power series expansions

abstract: When $A$ is a matrix with all eigenvalues in the disk $|z-1|<1$, the principal $p$th root of $A$ can be computed by Schr\"oder's method, among many other methods. In this paper we present a further study of Schr\"oder's method for the matrix $p$th root, through an examination of power series expansions of some sequences of scalar functions. Specifically, we obtain a new and informative error estimate for the matrix sequence generated by the Schr\"oder's method, a monotonic convergence result when $A$ is a nonsingular $M$-matrix, and a structure preserving result when $A$ is a nonsingular $M$-matrix or a real nonsingular $H$-matrix with positive diagonal entries.


title: Last-Iterate Convergence: Zero-Sum Games and Constrained Min-Max Optimization

abstract: Motivated by applications in Game Theory, Optimization, and Generative Adversarial Networks, recent work of Daskalakis et al \cite{DISZ17} and follow-up work of Liang and Stokes \cite{LiangS18} have established that a variant of the widely used Gradient Descent/Ascent procedure, called "Optimistic Gradient Descent/Ascent (OGDA)", exhibits last-iterate convergence to saddle points in {\em unconstrained} convex-concave min-max optimization problems. We show that the same holds true in the more general problem of {\em constrained} min-max optimization under a variant of the no-regret Multiplicative-Weights-Update method called "Optimistic Multiplicative-Weights Update (OMWU)". This answers an open question of Syrgkanis et al \cite{SALS15}. The proof of our result requires fundamentally different techniques from those that exist in no-regret learning literature and the aforementioned papers. We show that OMWU monotonically improves the Kullback-Leibler divergence of the current iterate to the (appropriately normalized) min-max solution until it enters a neighborhood of the solution. Inside that neighborhood we show that OMWU is locally (asymptotically) stable converging to the exact solution. We believe that our techniques will be useful in the analysis of the last iterate of other learning algorithms.


title: The localization theorem for framed motivic spaces

abstract: We prove the analog of the Morel-Voevodsky localization theorem for framed motivic spaces. We deduce that framed motivic spectra are equivalent to motivic spectra over arbitrary schemes, and we give a new construction of the motivic cohomology of arbitrary schemes.


title: On persistence of superoscillations for Schr\"{o}dinger equation with time-dependent quadratic Hamiltonians

abstract: In this work we study the persistence in time of superoscillations for the Schr\"{o}dinger equation with quadratic time-dependent Hamiltonians. We have solved explicitly the Cauchy initial value problem with three different kind of oscillatory initial data. In order to prove the persistence of superoscillations we have defined explicitly an operator in terms of solutions of a Riccati system associated with the variable coefficients of the Hamiltonian. The operator is defined on a space of entire functions. Particular examples include Caldirola-Kanai and degenerate parametric harmonic oscillator Hamiltonians and more. For these examples we have illustrated numerically the convergence on real and imaginary parts.


title: On the Fundamental Limits of Coded Data Shuffling for Distributed Machine Learning

abstract: We consider the data shuffling problem in a distributed learning system, in which a master node is connected to a set of worker nodes, via a shared link, in order to communicate a set of files to the worker nodes. The master node has access to a database of files. In every shuffling iteration, each worker node processes a new subset of files, and has excess storage to partially cache the remaining files, assuming the cached files are uncoded. The caches of the worker nodes are updated every iteration, and they should be designed to satisfy any possible unknown permutation of the files in subsequent iterations. For this problem, we characterize the exact load-memory trade-off for worst-case shuffling by deriving the minimum communication load for a given storage capacity per worker node. As a byproduct, the exact load-memory trade-off for any shuffling is characterized when the number of files is equal to the number of worker nodes. We propose a novel deterministic coded shuffling scheme, which improves the state of the art, by exploiting the cache memories to create coded functions that can be decoded by several worker nodes. Then, we prove the optimality of our proposed scheme by deriving a matching lower bound and showing that the placement phase of the proposed coded shuffling scheme is optimal over all shuffles.


title: Electron-phonon coupling in undoped cuprate $YBa_2Cu_3O_6$ estimated from Raman and optical conductivity spectra

abstract: We study experimentally the Raman response of the undoped high-Tc parent compound $YBa_2Cu_3O_6$, and give a unified theory of the two-magnon Raman peak and optical conductivity based on the Hubbard-Holstein model with electron-phonon coupling (EPC). The Hubbard model without EPC can qualitatively account for the experimentally observed resonance of the Raman response, but only the Hubbard-Holstein model (i) reproduces asymmetry of the Raman spectrum, (ii) validates experimental visibility of the two-magnon peak, and (iii) predicts the correct shape and energy of the lower edge of the charge transfer gap in optical conductivity. Comparison of experiments with the theory gives the EPC strength $\lambda$ = 0.6. This result convincingly indicates the vital role of EPC in high-Tc cuprates providing a clue to the mechanism of high-Tc.


title: Fundamental solution for super-critical non-symmetric L\'evy-type operators

abstract: We prove the existence and give estimates of the fundamental solution (the heat kernel) for the equation $\partial_t =\mathcal{L}^{\kappa}$ for non-symmetric non-local operators $$ \mathcal{L}^{\kappa}f(x):= \int_{\mathbb{R}^d}( f(x+z)-f(x)- 1_{|z|<1} \left<z,\nabla f(x)\right>)\kappa(x,z)J(z)\, dz\,, $$ under broad assumptions on $\kappa$ and $J$. Of special interest is the case when the order of the operator $\mathcal{L}^{\kappa}$ is smaller than or equal to 1. Our approach rests on imposing suitable cancellation conditions on the internal drift coefficient $$ \int_{r\leq |z|<1} z \kappa(x,z)J(z)dz\,,\qquad 0<r\leq 1\,, $$ which allows us to handle the non-symmetry of $z\mapsto \kappa(x,z)J(z)$. The results are new even for the $1$-stable L\'evy measure $J(z)=|z|^{-d-1}$.


title: Mechanisms of DNA Hybridization: Transition Path Analysis of a Simulation-Informed Markov Model

abstract: Complementary DNA strands in solution reliably hybridize to form stable duplexes. We study the kinetics of the hybridization process and the mechanisms by which two initially isolated strands come together to form a stable double helix. We adopt a multi-step computational approach. First, we perform a large number of Brownian dynamics simulations of the hybridization process using the coarse-grained oxDNA2 model. Second, we use these simulations to construct a Markov State Model of DNA dynamics that uses a state decomposition based on the inter-strand hydrogen bonding pattern. Third, we take advantage of Transition Path Theory to obtain quantitative information about the thermodynamic and dynamic properties of the hybridization process. We find that while there is a large ensemble of possible hybridization pathways there is a single dominant mechanism in which an initial base pair forms close to either end of the nascent double helix, and the remaining bases pair sequentially in a zipper-like fashion. We also show that the number of formed base pairs by itself is insufficient to describe the transition state of the hybridization process.


title: The fundamentals of quantum machine learning

abstract: Within the past few years, we have witnessed the rising of quantum machine learning (QML) models which infer electronic properties of molecules and materials, rather than solving approximations to the electronic Schrodinger equation. The increasing availability of large quantum mechanics reference data sets have enabled these developments. We review the basic theories and key ingredients of popular QML models such as choice of regressor, data of varying trustworthiness, the role of the representation, and the effect of training set selection. Throughout we emphasize the indispensable role of learning curves when it comes to the comparative assessment of different QML models.


title: Phase structures emerging from holography with Einstein gravity -- dilaton models at finite temperature

abstract: Asymptotic AdS Riemann space-times in five dimensions with a black brane (horizon) sourced by a fully back-reacted scalar field (dilaton) offer -- via the holographic dictionary -- various options for the thermodynamics of the flat four-dimensional boundary theory, uncovering Hawking-Page, first-order and second-order phase transitions up to a cross-over or featureless behavior. The relation of these phase structures to the dilaton potential is clarified and illustrating examples are presented. Having in mind applications to QCD we study probe vector mesons with the goal to figure out conditions for forming Regge type series of radial excitations and address the issue of meson melting.


title: Phase Retrieval Under a Generative Prior

abstract: The phase retrieval problem asks to recover a natural signal $y_0 \in \mathbb{R}^n$ from $m$ quadratic observations, where $m$ is to be minimized. As is common in many imaging problems, natural signals are considered sparse with respect to a known basis, and the generic sparsity prior is enforced via $\ell_1$ regularization. While successful in the realm of linear inverse problems, such $\ell_1$ methods have encountered possibly fundamental limitations, as no computationally efficient algorithm for phase retrieval of a $k$-sparse signal has been proven to succeed with fewer than $O(k^2\log n)$ generic measurements, exceeding the theoretical optimum of $O(k \log n)$. In this paper, we propose a novel framework for phase retrieval by 1) modeling natural signals as being in the range of a deep generative neural network $G : \mathbb{R}^k \rightarrow \mathbb{R}^n$ and 2) enforcing this prior directly by optimizing an empirical risk objective over the domain of the generator. Our formulation has provably favorable global geometry for gradient methods, as soon as $m = O(kd^2\log n)$, where $d$ is the depth of the network. Specifically, when suitable deterministic conditions on the generator and measurement matrix are met, we construct a descent direction for any point outside of a small neighborhood around the unique global minimizer and its negative multiple, and show that such conditions hold with high probability under Gaussian ensembles of multilayer fully-connected generator networks and measurement matrices. This formulation for structured phase retrieval thus has two advantages over sparsity based methods: 1) deep generative priors can more tightly represent natural signals and 2) information theoretically optimal sample complexity. We corroborate these results with experiments showing that exploiting generative models in phase retrieval tasks outperforms sparse phase retrieval methods.


title: The 2018 Census of Luminous Blue Variables

abstract: We present a census of all known and candidate luminous blue variables in the Local Group.


title: Knowledge Compilation, Width and Quantification

abstract: We generalize many results concerning the tractability of SAT and #SAT on bounded treewidth CNF-formula in the context of Quantified Boolean Formulas (QBF). To this end, we start by studying the notion of width for OBDD and observe that the blow up in size while existentially or universally projecting a block of variables in an OBDD only affects its width. We then generalize this notion of width to the more general representation of structured (deterministic) DNNF and give a similar algorithm to existentially or universally project a block of variables. Using a well-known algorithm transforming bounded treewidth CNF formula into deterministic DNNF, we are able to generalize this connection to quantified CNF which gives us as a byproduct that one can count the number of models of a bounded treewidth and bounded quantifier alternation quantified CNF in FPT time. We also give an extensive study of bounded width d-DNNF and proves the optimality of several of our results.


title: The Unification of Nonassociative Structures

abstract: The main non-associative algebras are Lie algebras and Jordan algebras. There are several ways to unify these non-associative algebras and associative algebras.


title: Photon-mediated interactions between quantum emitters in a diamond nanocavity

abstract: Photon-mediated interactions between quantum systems are essential for realizing quantum networks and scalable quantum information processing. We demonstrate such interactions between pairs of silicon-vacancy (SiV) color centers strongly coupled to a diamond nanophotonic cavity. When the optical transitions of the two color centers are tuned into resonance, the coupling to the common cavity mode results in a coherent interaction between them, leading to spectrally-resolved superradiant and subradiant states. We use the electronic spin degrees of freedom of the SiV centers to control these optically-mediated interactions. Our experiments pave the way for implementation of cavity-mediated quantum gates between spin qubits and for realization of scalable quantum network nodes.


title: Accurate cosmic shear errors: do we need ensembles of simulations?

abstract: Accurate inference of cosmology from weak lensing shear requires an accurate shear power spectrum covariance matrix. Here, we investigate this accuracy requirement and quantify the relative importance of the Gaussian (G), super-sample covariance (SSC) and connected non-Gaussian (cNG) contributions to the covariance. Specifically, we forecast cosmological parameter constraints for future wide-field surveys and study how different covariance matrix components affect parameter bounds. Our main result is that the cNG term represents only a small and potentially negligible contribution to statistical parameter errors: the errors obtained using the G+SSC subset are within $\lesssim 5\%$ of those obtained with the full G+SSC+cNG matrix for a Euclid-like survey. This result also holds for the shear two-point correlation function, variations in survey specifications and for different analytical prescriptions of the cNG term. The cNG term is that which is often tackled using numerically expensive ensembles of survey realizations. Our results suggest however that the accuracy of analytical or approximate numerical methods to compute the cNG term is likely to be sufficient for cosmic shear inference from the next generation of surveys.


title: Fault-tolerant quantum metrology

abstract: We introduce the notion of fault-tolerant quantum metrology to overcome noise beyond our control -- associated with sensing the parameter, by reducing the noise in operations under our control -- associated with preparing and measuring probes and ancillae. To that end, we introduce noise thresholds to quantify the noise resilience of parameter estimation schemes. We demonstrate improved noise thresholds over the non-fault-tolerant schemes. We use quantum Reed-Muller codes to retrieve more information about a single phase parameter being estimated in the presence of full-rank Pauli noise. Using only error detection, as opposed to error correction, allows us to retrieve higher thresholds. We show that better devices, which can be engineered, can enable us to counter larger noise in the field beyond our control. Further improvements in fault-tolerant quantum metrology could be achieved by optimising in tandem parameter-specific estimation schemes and transversal quantum error correcting codes.


title: A proof of the elliptical range theorem via Kippenhahn's theorem

abstract: The elliptical range theorem asserts that the field of values (or numerical range) of a two-by-two matrix with complex entries is an elliptical disk, the foci of which are the eigenvalues of the given matrix. Many proofs of this result are available in the literature, but most, with one exception, are computational and quite involved. In this note, it is shown that the elliptical range theorem follows from the properties of plane algebraic curves and a straightforward application of a well-known result due to Kippenhahn.


title: Applying a Dynamical Systems Model and Network Theory to Major Depressive Disorder

abstract: Mental disorders like major depressive disorder can be seen as complex dynamical systems. In this study we investigate the dynamic behaviour of individuals to see whether or not we can expect a transition to another mood state. We introduce a mean field model to a binomial process, where we reduce a dynamic multidimensional system (stochastic cellular automaton) to a one-dimensional system to analyse the dynamics. Using maximum likelihood estimation, we can estimate the parameter of interest which, in combination with a bifurcation diagram, reflects the expectancy that someone has to transition to another mood state. After validating the proposed method with simulated data, we apply this method to two empirical examples, where we show its use in a clinical sample consisting of patients diagnosed with major depressive disorder, and a general population sample. Results showed that the majority of the clinical sample was categorized as having an expectancy for a transition, while the majority of the general population sample did not have this expectancy. We conclude that the mean field model has great potential in assessing the expectancy for a transition between mood states. With some extensions it could, in the future, aid clinical therapists in the treatment of depressed patients.


title: Attack and defence in cellular decision-making: lessons from machine learning

abstract: Machine learning algorithms can be fooled by small well-designed adversarial perturbations. This is reminiscent of cellular decision-making where ligands (called antagonists) prevent correct signalling, like in early immune recognition. We draw a formal analogy between neural networks used in machine learning and models of cellular decision-making (adaptive proofreading). We apply attacks from machine learning to simple decision-making models, and show explicitly the correspondence to antagonism by weakly bound ligands. Such antagonism is absent in more nonlinear models, which inspired us to implement a biomimetic defence in neural networks filtering out adversarial perturbations. We then apply a gradient-descent approach from machine learning to different cellular decision-making models, and we reveal the existence of two regimes characterized by the presence or absence of a critical point for the gradient. This critical point causes the strongest antagonists to lie close to the decision boundary. This is validated in the loss landscapes of robust neural networks and cellular decision-making models, and observed experimentally for immune cells. For both regimes, we explain how associated defence mechanisms shape the geometry of the loss landscape, and why different adversarial attacks are effective in different regimes. Our work connects evolved cellular decision-making to machine learning, and motivates the design of a general theory of adversarial perturbations, both for in vivo and in silico systems.


title: A quantum-inspired classical algorithm for recommendation systems

abstract: We give a classical analogue to Kerenidis and Prakash's quantum recommendation system, previously believed to be one of the strongest candidates for provably exponential speedups in quantum machine learning. Our main result is an algorithm that, given an $m \times n$ matrix in a data structure supporting certain $\ell^2$-norm sampling operations, outputs an $\ell^2$-norm sample from a rank-$k$ approximation of that matrix in time $O(\text{poly}(k)\log(mn))$, only polynomially slower than the quantum algorithm. As a consequence, Kerenidis and Prakash's algorithm does not in fact give an exponential speedup over classical algorithms. Further, under strong input assumptions, the classical recommendation system resulting from our algorithm produces recommendations exponentially faster than previous classical systems, which run in time linear in $m$ and $n$. The main insight of this work is the use of simple routines to manipulate $\ell^2$-norm sampling distributions, which play the role of quantum superpositions in the classical setting. This correspondence indicates a potentially fruitful framework for formally comparing quantum machine learning algorithms to classical machine learning algorithms.


title: Towards a Complete Picture of Stationary Covariance Functions on Spheres Cross Time

abstract: With the advent of wide-spread global and continental-scale spatiotemporal datasets, increased attention has been given to covariance functions on spheres over time. This paper provides results for stationary covariance functions of random fields defined over $d$-dimensional spheres cross time. Specifically, we provide a bridge between the characterization in \cite{berg-porcu} for covariance functions on spheres cross time and Gneiting's lemma \citep{gneiting2002} that deals with planar surfaces. We then prove that there is a valid class of covariance functions similar in form to the Gneiting class of space-time covariance functions \citep{gneiting2002} that replaces the squared Euclidean distance with the great circle distance. Notably, the provided class is shown to be positive definite on every $d$-dimensional sphere cross time, while the Gneiting class is positive definite over $\R^d \times \R$ for fixed $d$ only. In this context, we illustrate the value of our adapted Gneiting class by comparing examples from this class to currently established nonseparable covariance classes using out-of-sample predictive criteria. These comparisons are carried out on two climate reanalysis datasets from the National Centers for Environmental Prediction and National Center for Atmospheric Research. For these datasets, we show that examples from our covariance class have better predictive performance than competing models.


title: A single photon detection system for visible and infrared spectrum range

abstract: We demonstrate niobium nitride based superconducting single-photon detectors sensitive in the spectral range 452 nm - 2300 nm. The system performance was tested in a real-life experiment with correlated photons generated by means of spontaneous parametric down conversion, where one of photon was in the visible range and the other was in the infrared range. We measured a signal to noise ratio as high as $4\times 10^4$ in our detection setting. A photon detection efficiency as high as 64% at 1550 nm and 15% at 2300 nm was observed.

categories: physics.ins-det cond-mat.supr-con physics.optics quant-ph


title: Gluing I: Integrals and Symmetries

abstract: We review some aspects of the cutting and gluing law in local quantum field theory. In particular, we emphasize the description of gluing by a path integral over a space of polarized boundary conditions, which are given by leaves of some Lagrangian foliation in the phase space. We think of this path integral as a non-local $(d-1)$-dimensional gluing theory associated to the parent local $d$-dimensional theory. We describe various properties of this procedure and spell out conditions under which symmetries of the parent theory lead to symmetries of the gluing theory. The purpose of this paper is to set up a playground for the companion paper where these techniques are applied to obtain new results in supersymmetric theories.


title: Modelling the coincident observation of a high-energy neutrino and a bright blazar flare

abstract: In September 2017, the IceCube Neutrino Observatory recorded a very-high-energy neutrino in directional coincidence with a blazar in an unusually bright gamma-ray state, TXS0506+056. Blazars are prominent photon sources in the universe because they harbor a relativistic jet whose radiation is strongly collimated and amplified. High-energy atomic nuclei known as cosmic rays can produce neutrinos; thus the recent detection may help identifying the sources of the diffuse neutrino flux and the energetic cosmic rays. Here we report on a self-consistent analysis of the physical relation between the observed neutrino and the blazar, in particular the time evolution and spectral behavior of neutrino and photon emission. We demonstrate that a moderate enhancement in the number of cosmic rays during the flare can yield a very strong increase of the neutrino flux which is limited by co-produced hard X-rays and TeV gamma rays. We also test typical radiation models for compatibility and identify several model classes as incompatible with the observations. We investigate to what degree the findings can be generalized to the entire population of blazars, to determine the relation between their output in photons, neutrinos, and cosmic rays, and suggest how to optimize the strategy of future observations.


title: Entanglement Holonomies

abstract: We introduce a quantum notion of parallel transport between subsystems of a quantum state whose holonomies characterize the structure of entanglement. In AdS/CFT, entanglement holonomies are reflected in the bulk spacetime connection. When the subsystems are a pair of holographic CFTs in an entangled state, our quantum transport measures Wilson lines threading the dual wormhole. For subregions of a single CFT it is generated by the modular Berry connection and computes the effect of the AdS curvature on the transport of minimal surfaces. Our observation reveals a new aspect of the spacetime-entanglement duality and yet another concept shared between gravity and quantum mechanics.


title: Gravitational lensing detection of an extremely dense environment around a galaxy cluster

abstract: Galaxy clusters form at the highest density nodes of the cosmic web. The clustering of massive halos is enhanced relative to the general mass distribution and matter beyond the virial region is strongly correlated to the halo mass (halo bias). Clustering can be further enhanced depending on halo properties other than mass (secondary bias). The questions of how much and why the regions surrounding rich clusters are over-dense are still unanswered. Here, we report the analysis of the environment bias in a sample of very massive clusters, selected through the Sunyaev-Zel'dovich effect by the Planck mission. We present the first detection of the correlated dark matter associated to a single cluster, PSZ2 G099.86+58.45. The system is extremely rare in the current paradigm of structure formation. The gravitational lensing signal was traced up to 30 megaparsecs with high signal-to-noise ratio ~3.4. The measured shear is very large and points at environment matter density in notable excess of the cosmological mean. The boosting of the correlated dark matter density around high mass halos can be very effective. Together with ensemble studies of the large scale structure, lensing surveys can picture the surroundings of single haloes.


title: Gluing II: Boundary Localization and Gluing Formulas

abstract: We describe applications of the gluing formalism discussed in the companion paper. When a $d$-dimensional local theory $\text{QFT}_d$ is supersymmetric, and if we can find a supersymmetric polarization for $\text{QFT}_d$ quantized on a $(d-1)$-manifold $W$, gluing along $W$ is described by a non-local $\text{QFT}_{d-1}$ that has an induced supersymmetry. Applying supersymmetric localization to $\text{QFT}_{d-1}$, which we refer to as the boundary localization, allows in some cases to represent gluing by finite-dimensional integrals over appropriate spaces of supersymmetric boundary conditions. We follow this strategy to derive a number of `gluing formulas' in various dimensions, some of which are new and some of which have been previously conjectured. First we show how gluing in supersymmetric quantum mechanics can reduce to a sum over a finite set of boundary conditions. Then we derive two gluing formulas for 3D $\mathcal{N}=4$ theories on spheres: one providing the Coulomb branch representation of gluing, and another providing the Higgs branch representation. This allows to study various properties of their $(2,2)$-preserving boundary conditions in relation to Mirror Symmetry. After that we derive a gluing formula in 4D $\mathcal{N}=2$ theories on spheres, both squashed and round. First we apply it to predict the hemisphere partition function, then we apply it to the study of boundary conditions and domain walls in these theories. Finally, we mention how to glue half-indices of 4D $\mathcal{N}=2$ theories.


title: The Bearable Compositeness of Leptons

abstract: Partial compositeness as a theory of flavor in the lepton sector is assessed. We begin presenting the first systematic analysis of neutrino mass generation in this context, and identifying the distinctive mass textures. We then update the bounds from charged lepton flavor and CP violating observables. We put forward a $U(1)^3\times CP$ symmetry of the composite sector, in order to allow the new physics to be not far above the TeV scale. This hypothesis effectively suppresses the new contributions to the electron EDM and $\mu\to e\gamma$, by far the most constraining observables, and results in a novel pattern of flavor violation and neutrino masses. The CP violation in the elementary-composite mixing is shown to induce a CKM phase of the correct size, as well as order-one phases in the PMNS matrix. We compare with the alternative possibility of introducing multiple scales of compositeness for leptons, that also allow to evade flavor and CP constraints. Finally, we examine violations of lepton flavor universality in $B$-meson semi-leptonic decays. The neutral-current anomalies can be accommodated, predicting strong correlations among different lepton flavors, with a few channels close to the experimental sensitivity.


title: SU Lyn: diagnosing the boundary layer with UV and hard X-ray data

abstract: Symbiotic stars in which the symbiotic phenomenon is powered solely by accretion, often at an average rate that is higher than in cataclysmic variable stars, provide an important opportunity to diagnose boundary layers around disk-accreting white dwarfs. Here we investigate SU Lyncis, a recently discovered example of a purely accretion-powered symbiotic star, using the first reliable X-ray spectroscopy, obtained with NuSTAR, and UV photometry obtained with Swift. SU Lyn has hard, thermal, X-ray emission that is strongly affected by a variable local absorber - that has little impact on the UV emission. Its X-ray spectrum is described well using a plasma cooling from $k$T $\approx$ 21 keV, with a 3 to 30 keV luminosity of approximately 4.9$\times$10$^{32}$ ergs s$^{-1}$. The spectrum is also consistent with the presence of reflection with an amplitude of 1.0, although in that case, the best-fit plasma temperature is 20-25% lower. The UV to X-ray luminosity ratio of SU Lyn changed significantly between 2015 and 2016. We interpret this as a consequence of a drop by almost 90% in the accretion rate. Whereas the UV luminosity of the disk responded linearly, the luminosity of the optically thin (hard X-ray) emission from the boundary layer remained roughly constant because the boundary layer changed from partially optically thick to almost completely optically thin. Under this interpretation, we place a lower limit on the white dwarf mass of 0.7 M$_{\odot}$ (0.8 M$_{\odot}$ if we neglect reflection).


title: The luminous host galaxy, faint supernova and rapid afterglow rebrightening of GRB 100418A

abstract: Long gamma-ray bursts give us the chance to study both their extreme physics and the star-forming galaxies in which they form. GRB 100418A, at a z = 0.6239, had a bright optical and radio afterglow, and a luminous star-forming host galaxy. This allowed us to study the radiation of the explosion as well as the interstellar medium of the host both in absorption and emission. We collected photometric data from radio to X-ray wavelengths to study the evolution of the afterglow and the contribution of a possible supernova and three X-shooter spectra obtained during the first 60 hr. The light curve shows a very fast optical rebrightening, with an amplitude of 3 magnitudes, starting 2.4 hr after the GRB onset. This cannot be explained by a standard external shock model and requires other contributions, such as late central-engine activity. Two weeks after the burst we detect an excess in the light curve consistent with a SN with peak absolute magnitude M_V = -18.5 mag, among the faintest GRB-SNe detected to date. The host galaxy shows two components in emission, with velocities differing by 130 km s^-1, but otherwise having similar properties. While some absorption and emission components coincide, the absorbing gas spans much higher velocities, indicating the presence of gas beyond the star-forming regions. The host has a star-formation rate of 12.2 M_sol yr^-1, a metallicity of 12 + log(O/H) = 8.55 and a mass of 1.6x10^9 M_sol. GRB 100418A is a member of a class of afterglow light curves which show a steep rebrightening in the optical during the first day, which cannot be explained by traditional models. Its very faint associated SN shows that GRB-SNe can have a larger dispersion in luminosities than previously seen. Furthermore, we have obtained a complete view of the host of GRB 100418A owing to its spectrum, which contains a remarkable number of both emission and absorption lines.


title: A Thermal Neutrino Portal to Sub-MeV Dark Matter

abstract: Thermal relics lighter than an MeV contribute to the energy density of the universe at the time of nucleosynthesis and recombination. Constraints on extra radiation degrees of freedom typically exclude even the simplest of such dark sectors. We explore the possibility that a sub-MeV dark sector entered equilibrium with the Standard Model after neutrino-photon decoupling, which significantly weakens these constraints and naturally arises in the context of neutrino mass generation through the spontaneous breaking of lepton number. Acquiring an adequate dark matter abundance independently motivates the MeV-scale in these models through the coincidence of gravitational, matter-radiation equality, and neutrino mass scales, $(m_\text{Pl} / T^\text{MRE})^{1/4} \, m_\nu \sim \text{MeV}$. This class of scenarios will be decisively tested by future measurements of the cosmic microwave background and matter structure of the universe. While the dark sector dominantly interacts with Standard Model neutrinos, large couplings to nucleons are possible in principle, leading to observable signals at proposed low-threshold direct detection experiments.


title: Majorana Neutrino Masses in the RGEs for Lepton Flavour Violation

abstract: We suppose that the observed neutrino masses can be parametrised by a lepton number violating dimension-five operator, and calculate the mixing of double insertions of this operator into lepton flavour changing dimension-six operators of the standard model effective theory. This allows to predict the log-enhanced, but $m_\nu^2$-suppressed lepton flavour violation that is generic to high-scale Majorana neutrino mass models. We also consider the Two Higgs Doublet Model, where the second Higgs allows the construction of three additional dimension-five operators, and evaluate the corresponding anomalous dimensions. The sensitivity of current searches for lepton flavour violation to these additional Wilson coefficients is then examined.


title: Signs of heavy Higgs bosons at CLIC: An $e^+ e^-$ road to the Electroweak Phase Transition

abstract: We analyse the sensitivity of the proposed Compact Linear Collider (CLIC) to the existence of beyond the Standard Model (SM) Higgs bosons through their decays into pairs of massive gauge bosons $H \to VV$ and SM-like Higgses $H \to hh$, considering CLIC centre of mass energies $\sqrt{s} = 1.4$ TeV and $3$ TeV. We find that resonant di-Higgs searches at CLIC would allow for up to two orders of magnitude improvement w.r.t. the sensitivity achievable by HL-LHC in the mass range $m_H \in [250\,\mathrm{GeV},\, 1 \,\mathrm{TeV}]$. Focusing then on a real singlet extension of the SM, we explore the prospects of heavy Higgs searches at CLIC for probing the regions of parameter space yielding a strongly first order electroweak phase transition that could generate the observed matter-antimatter asymmetry of the Universe. Our study illustrates the complementarity between CLIC and other possible future colliders like FCC-ee in probing singlet extensions of the SM, and shows that high-energy $e^+ e^-$ colliders provide a powerful means to unravel the nature of electroweak symmetry breaking in the early Universe.


title: Analytically solvable renormalization group for the many-body localization transition

abstract: We introduce a simple, exactly solvable strong-randomness renormalization group (RG) model for the many-body localization (MBL) transition in one dimension. Our approach relies on a family of RG flows parametrized by the asymmetry between thermal and localized phases. We identify the physical MBL transition in the limit of maximal asymmetry, reflecting the instability of MBL against rare thermal inclusions. We find a critical point that is localized with power-law distributed thermal inclusions. The typical size of critical inclusions remains finite at the transition, while the average size is logarithmically diverging. We propose a two-parameter scaling theory for the many-body localization transition that falls into the Kosterlitz-Thouless universality class, with the MBL phase corresponding to a stable line of fixed points with multifractal behavior.

categories: cond-mat.dis-nn cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el


title: The HP2 Survey - IV. The Pipe nebula: Effective dust temperatures in dense cores

abstract: Multi-wavelength observations in the sub-mm regime provide information on the distribution of both the dust column density and the effective dust temperature in molecular clouds. In this study, we created high-resolution and high-dynamic-range maps of the Pipe nebula region and explored the value of dust-temperature measurements in particular towards the dense cores embedded in the cloud. The maps are based on data from the Herschel and Planck satellites, and calibrated with a near-infrared extinction map based on 2MASS observations. We have considered a sample of previously defined cores and found that the majority of core regions contain at least one local temperature minimum. Moreover, we observed an anti-correlation between column density and temperature. The slope of this anti-correlation is dependent on the region boundaries and can be used as a metric to distinguish dense from diffuse areas in the cloud if systematic effects are addressed appropriately. Employing dust-temperature data thus allows us to draw conclusions on the thermodynamically dominant processes in this sample of cores: external heating by the interstellar radiation field and shielding by the surrounding medium. In addition, we have taken a first step towards a physically motivated core definition by recognising that the column-density-temperature anti-correlation is sensitive to the core boundaries. Dust-temperature maps therefore clearly contain valuable information about the physical state of the observed medium.


title: Hacking Alice's box in CV-QKD

abstract: Security analyses of quantum cryptographic protocols typically rely on certain conditions; one such condition is that the sender (Alice) and receiver (Bob) have isolated devices inaccessible to third parties. If an eavesdropper (Eve) has a side-channel into one of the devices, then the key rate may be sensibly reduced. In this paper, we consider an attack on a coherent-state protocol, where Eve not only taps the main communication channel but also hacks Alice's device. This is done by introducing a Trojan horse mode with low mean number of photons $\bar{n}$ which is then modulated in a similar way to the signal state. First we show that this strategy can be reduced to an attack without side channels but with higher loss and noise in the main channel. Then we show how the key rate rapidly deteriorates for increasing photons $\bar{n}$, being halved at long distances each time $\bar{n}+1$ doubles. Our work suggests that Alice's device should also be equipped with sensing systems that are able to detect and estimate the total number of incoming and outgoing photons.

categories: quant-ph cond-mat.other physics.ins-det physics.optics


title: Cosmological simulations of dwarfs: the need for ISM physics beyond SN feedback alone

abstract: The dominant feedback mechanism in low mass haloes is usually assumed to take the form of massive stars exploding as supernovae (SNe). We perform very high resolution cosmological zoom-in simulations of five dwarf galaxies to z = 4 with our mechanical SN feedback model. This delivers the correct amount of momentum corresponding to the stage of the SN remnant evolution resolved, and has been shown to lead to realistic dwarf properties in isolated simulations. We find that in 4 out of our 5 simulated cosmological dwarfs, SN feedback has insufficient impact resulting in excessive stellar masses, extremely compact sizes and central super-solar stellar metallicities. The failure of the SN feedback in our dwarfs is physical in nature within our model and is the result of the build up of very dense gas in the early universe due to mergers and cosmic inflows prior to the first SN occurring. We demonstrate that our results are insensitive to resolution (provided that it is high enough), details of the (spatially uniform) UV background and reasonable alterations within our star formation prescription. We therefore conclude that the ability of SNe to regulate dwarf galaxy properties is dependent on other physical processes, such as turbulent pressure support, clustering and runaway of SN progenitors and other sources of stellar feedback.


title: Searching for the Lowest Metallicity Galaxies

abstract: We report a method of identifying candidate low-metallicity blue compact dwarf galaxies (BCDs) from the Sloan Digital Sky Survey (SDSS) imaging data and present 3-m Lick Observatory and 10-m W.M. Keck Observatory optical spectroscopic observations of 94 new systems that have been discovered with this method. The candidate BCDs are selected from Data Release 12 (DR12) of SDSS based on their photometric colors and morphologies. Using the Kast spectrometer on the 3-m telescope, we confirm that the candidate low-metallicity BCDs are emission-line galaxies and we make metallicity estimates using the empirical R and S calibration methods. Follow-up observations on a subset of the lowest-metallicity systems are made at Keck using the Low Resolution Imaging Spectrometer (LRIS), which allow for a direct measurement of the oxygen abundance. We determine that 45 of the reported BCDs are low-metallicity candidates with 12 + log(O/H) <= 7.65, including six systems which are either confirmed or projected to be among the lowest-metallicity galaxies known, at 1/30 of the solar oxygen abundance, or 12 + log(O/H) ~ 7.20.


title: The Halo's Ancient Metal-Rich Progenitor Revealed with BHB Stars

abstract: Using the data from the Sloan Digital Sky Survey and the Gaia satellite, we assemble a pure sample of $\sim$3000 Blue Horizontal Branch (BHB) stars with 7-D information, including positions, velocities and metallicities. We demonstrate that, as traced with BHBs, the Milky Way's stellar halo is largely unmixed and can not be well represented with a conventional Gaussian velocity distribution. A single-component model fails because the inner portions of the halo are swamped with metal-rich tidal debris from an ancient, head-on collision, known as the "Gaia Sausage". Motivated by the data, we build a flexible mixture model which allows us to track the evolution of the halo make-up across a wide range of radii. It is built from two components, one representing the radially anisotropic Sausage stars with their lobed velocity distribution, the other representing a more metal-poor and more isotropic component built up from minor mergers. We show that inside 25 kpc the "Sausage" contributes at least 50 % of the Galactic halo. The fraction of "Sausage" stars diminishes sharply beyond 30 kpc, which is the long-established break radius of the classical stellar halo.


title: Flat Rotation Curves found in Merging Dusty Starbursts at $z=2.3$ through Tilted-Ring Modeling

abstract: The brightest 500$\,\mu$m source in the XMM field, HXMM01, is a rare merger of luminous starburst galaxies at $z=2.3$ with a dust-obscured star-formation rate of 2,000$\,M_{\odot}\,{\rm yr}^{-1}$. Here we present high-resolution spectroscopic observations of HXMM01 with the Atacama Large Millimeter/submillimeter Array (ALMA). We detect line emission from ${\rm CO\,{\it J}=7\to6}$, [C I]$\,{{^3P_2}\to{^3P_1}}$, and p-${\rm H_{2}O}\,{2_{11}}\to{2_{02}}$ and continuum emission at $230\,$GHz. At a spatial resolution of 0.2" and a spectral resolution of 40$\,\rm km\,s^{-1}$, the source is resolved into three distinct components, which are spatially and dynamically associated within a projected radius of 20$\,$kpc and a radial velocity range of 2,000$\,\rm km\,s^{-1}$. For two major components, our Bayesian-based tilted-ring modeling of the ALMA spectral cubes shows almost flat rotation curves peaking at $\sim500\,\rm km\,s^{-1}$ at galactocentric distances between 2 and 5$\,$kpc. Each of them has a dynamical mass of $\sim10^{11}\,M_\odot$. The combination of the dynamical masses and the archival ${\rm CO\,{\it J}=1\to0}$ data places strong upper limits on the CO$\to$H$_2$ conversion factor of $\alpha_{\rm CO}\lesssim1.4-2.0\,{M_{\odot}}\,\rm (K\,km\,s^{-1}\,pc^{2})^{-1}$. These limits are significantly below the Galactic inner disk $\alpha_{\rm CO}$ value of $4.3\,{M_{\odot}}\,\rm (K\,km\,s^{-1}\,pc^{2})^{-1}$ but are consistent with those of local starbursts. Therefore, the previously estimated short gas depletion timescale of $\sim200\,$Myr remains unchanged.


title: The Pristine Survey IV: Approaching the Galactic metallicity floor with the discovery of an ultra metal-poor star

abstract: The early Universe presented a star formation environment that was almost devoid of heavy elements. The lowest metallicity stars thus provide a unique window into the earliest Galactic stages, but are exceedingly rare and difficult to find. Here we present the discovery of an ultra-metal-poor star, Pristine_221.8781+9.7844, using narrow-band Ca H&K photometry from the Pristine survey. Follow-up medium and high-resolution spectroscopy confirms the ultra-metal-poor nature of Pristine_221.8781+9.7844 ([Fe/H] = -4.66 +/- 0.13 in 1D LTE) with an enhancement of 0.3-0.4 dex in alpha-elements relative to Fe, and an unusually low carbon abundance. We derive an upper limit of A(C) = 5.6, well below typical A(C) values for such ultra metal-poor stars. This makes Pristine_221.8781+9.7844 one of the most metal-poor stars; in fact, it is very similar to the most metal-poor star known (SDSS J102915+172927). The existence of a class of ultra metal-poor stars with low(er) carbon abundances suggest that there must have been several formation channels in the early Universe through which long-lived, low-mass stars were formed.


title: The Galactic WC and WO stars: The impact of revised distances from Gaia DR2 and their role as massive black hole progenitors

abstract: Wolf-Rayet stars of the carbon sequence (WC stars) are an important cornerstone in the late evolution of massive stars before their core collapse. As core-helium burning, hydrogen-free objects with huge mass-loss, they are likely the last observable stage before collapse and thus promising progenitor candidates for type Ib/c supernovae. Their strong mass-loss furthermore provides challenges and constraints to the theory of radiatively driven winds. Thus, the determination of the WC star parameters is of major importance for several astrophysical fields. With Gaia DR2, for the first time parallaxes for a large sample of Galactic WC stars are available, removing major uncertainties inherent to earlier studies. In this work, we re-examine the sample from Sander et al. (2012) to derive key properties of the Galactic WC population. All quantities depending on the distance are updated, while the underlying spectral analyses remain untouched. Contrasting earlier assumptions, our study yields that WC stars of the same subtype can significantly vary in absolute magnitude. With Gaia DR2, the picture of the Galactic WC population becomes more complex: We obtain luminosities ranging from log L = 4.9 to 6.0 with one outlier having log L = 4.7. This indicates that the WC stars are likely formed from a broader initial mass range than previously assumed. We obtain mass-loss rates ranging between log Mdot = -5.1 and -4.1, with Mdot propto L^0.68 and a linear scaling of the modified wind momentum with luminosity. We discuss the implications for stellar evolution, including unsolved issues regarding the need of envelope inflation to address the WR radius problem, and the open questions in regard to the connection of WR stars with Gamma-ray bursts. WC and WO stars are progenitors of massive black holes, collapsing either silently or in a supernova that most-likely has to be preceded by a WO stage.


title: Bit Threads in Higher Curvature Gravity

abstract: We generalize holographic bit threads to bulk theories with a gravitational action containing higher-curvature terms. Bit threads are a reformulation of holographic entanglement entropy, where the entropy is given by the maximum number of threads emanating from a boundary region into the bulk. We show that the addition of higher-curvature terms adds corrections to the bit thread thickness that depend on the local geometry and thread orientation. Two different methods are given: determination of the density bound by requiring the maximum number of threads through a given surface to reproduce the entanglement entropy functional on that surface, and application of Lagrange dualization. The results of the two methods are applied to Gauss-Bonnet gravity as the simplest non-trivial example.


title: Evidence of a past disc-disc encounter: HV and DO Tau

abstract: Theory and observations suggest that star formation occurs hierarchically due to the fragmentation of giant molecular clouds. In this case we would expect substructure and enhanced stellar multiplicity in the primordial cluster. This substructure is expected to decay quickly in most environments, however historic stellar encounters might leave imprints in a protoplanetary disc (PPD) population. In a low density environment such as Taurus, tidal tails from violent star-disc or disc-disc encounters might be preserved over time-scales sufficient to be observed. In this work, we investigate the possibility that just such an event occured between HV Tau C (itself a component of a triple system) and DO Tau $\sim 0.1$ Myr ago, as evidenced by an apparent `bridge' structure evident in the $160$ $\mu$m emission. By modelling the encounter using smoothed particle hydrodynamics (SPH) we reproduce the main features of the observed extended structure (`V'-shaped emission pointing west of HV Tau and a tail-like structure extending east of DO Tau). We suggest that HV Tau and DO Tau formed together in a quadruple system on a scale of $\sim 5000$ au ($0.025$ pc).


title: Geometry of Higgs-branch superconformal primary bundles

abstract: It is known that the two- and three-point functions of Higgs-branch superconformal primaries in 4d N=2 superconformal field theories obey non-renormalization theorems on N=2 superconformal manifolds. In this paper we prove a stronger statement, that the bundles of Higgs-branch superconformal primaries over N=2 superconformal manifolds are endowed with a flat connection, or equivalently that Higgs-branch superconformal primaries have vanishing Berry phases under N=2 exactly marginal deformations. This statement fits well with the proposed correspondence between the rigid structures of 2d chiral algebras and the sector of Schur operators in 4d N=2 theories. We also discuss the general interplay between non-renormalization theorems and the curvature of bundles of protected operators and provide a new simpler proof of the vanishing curvature of 1/2-BPS operators in 4d N=4 SYM theory that does not require the use of the 4d tt* equations.


title: An Experimental Quantum Bernoulli Factory

abstract: There has been a concerted effort to identify problems computable with quantum technology which are intractable with classical technology or require far fewer resources to compute. Recently, randomness processing in a Bernoulli factory has been identified as one such task. Here, we report two quantum photonic implementations of a Bernoulli factory, one utilising quantum coherence and single-qubit measurements and the other which uses quantum coherence and entangling measurements of two qubits. We show that the former consumes three orders of magnitude fewer resources than the best known classical method, while entanglement offers a further five-fold reduction. These concepts may provide a means for quantum enhanced-performance in the simulation of stochastic processes and sampling tasks.


title: Prevalence of neutral gas in centres of merging galaxies

abstract: We present Giant Metrewave Radio Telescope and Very Large Array observations of HI 21-cm absorption in ten z<=0.2 galaxy mergers that host strong radio sources. Seven of these mergers show strong absorption [with N(HI) ~10^{21-22} cm^{-2}, for spin temperature of 100 K and unit covering factor], leading to a total detection rate of 83 +/- 17% in low-z radio-loud galaxy mergers. This is ~3-4 times higher than that found for intrinsic HI 21-cm absorption in low-z radio sources in general. The fraction of intrinsic absorbers that are associated with mergers increases with increasing N(HI) threshold, i.e. ~40% and 100% of the absorbers with N(HI) >10^{21} cm^{-2} and >10^{22} cm^{-2} arise from mergers, respectively. The distribution of N(HI) among the mergers is significantly different from that found in isolated systems, with mergers giving rise to six times stronger absorption on average. The fraction of redshifted absorption components (with respect to the systemic velocity of the radio source obtained from optical emission lines) among mergers is found to be higher by two to three times compared to that found for non-interacting systems. Follow-up spatially-resolved multi-wavelength spectroscopy is essential to understand the exact connection between the presence of circumnuclear neutral gas and the AGN activity in these mergers.


title: A multiwavelength view of BL Lacs neutrino candidates

abstract: The discovery of high-energy astrophysical neutrinos by IceCube kicked off a new line of research to identify the electromagnetic counterparts producing these neutrinos. Among the extragalactic sources, active galactic nuclei (AGN), and in particular Blazars, are promising candidate neutrino emitters. Their structure, with a relativistic jet pointing to the Earth, offers a natural accelerator of particles and for this reason a perfect birthplace of high energy neutrinos. A good characterisation of the spectral energy distribution (SED) of these sources can improve the understanding of the physical composition of the source and the emission processes involved. Starting from our previous works in which we assumed a correlation between the $\gamma$-ray and the neutrino flux of the BL Lacs of the 2FHL catalogue (detected by Fermi above 50GeV), we select those BL Lac in spatial correlation with the IceCube events. We obtain a sample of 7 sources and we start an observational campaign to have a better characterisation of the synchrotron peak. During the analysis of the data a new source has been added because of its position inside the angular uncertainty of a muon track event detected by IceCube. This source, namely TXS0506+056, was in a high-state during the neutrino event and we will consider it as benchmark to check the proprieties of the other sources of the sample during the related neutrino detection. We obtain a better characterisation of the SED for the sources of our sample. A prospective extreme Blazar, a very peculiar low synchrotron peak (LSP) source with a large separation of the two peaks and a \textit{twin} of TXS0506+056 come up. We also provide the $\gamma$-ray light curve to check the trend of the sources around the neutrino detection but no clears patterns are in common among the sources.


title: The blazar TXS 0506+056 associated with a high-energy neutrino: insights into extragalactic jets and cosmic ray acceleration

abstract: A neutrino with energy of $\sim$290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS~0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at $\sim 3\sigma$ level. We monitored the object in the very-high-energy (VHE) band with the MAGIC telescopes for $\sim$41 hours from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability is clearly resolved. We interpret the quasi-simultaneous neutrino and broadband electromagnetic observations with a novel one-zone lepto-hadronic model, based on interactions of electrons and protons co-accelerated in the jet with external photons originating from a slow-moving plasma sheath surrounding the faster jet spine. We can reproduce the multiwavelength spectra of TXS 0506+056 with neutrino rate and energy compatible with IceCube-170922A, and with plausible values for the jet power of $\sim 10^{45} - 4 \times 10^{46} {\rm erg \ s^{-1}}$. The steep spectrum observed by MAGIC is concordant with internal $\gamma\gamma$ absorption above a few tens of GeV entailed by photohadronic production of a $\sim$290 TeV neutrino, corroborating a genuine connection between the multi-messenger signals. In contrast to previous predictions of predominantly hadronic emission from neutrino sources, the gamma-rays can be mostly ascribed to inverse Compton up-scattering of external photons by accelerated electrons. The X-ray and VHE bands provide crucial constraints on the emission from both accelerated electrons and protons. We infer that the maximum energy of protons in the jet co-moving frame can be in the range $\sim 10^{14}$ to $10^{18}$ eV.


title: The 55 Cnc system reassessed

abstract: Orbiting a bright, nearby star the 55 Cnc system offers a rare opportunity to study a multiplanet system that has a wide range of planetary masses and orbital distances. Using two decades of photometry and spectroscopy data, we have measured the rotation of the host star and its solar-like magnetic cycle. Accounting for this cycle in our velocimetric analysis of the system allows us to revise the properties of the outermost giant planet and its four planetary companions. The innermost planet 55 Cnc e is an unusually close-in super-Earth, whose transits have allowed for detailed follow-up studies. Recent observations favor the presence of a substantial atmosphere yet its composition, and the nature of the planet, remain unknown. We combined our derived planet mass (Mp = 8.0+-0.3 Mearth) with refined measurement of its optical radius derived from HST/STIS observations (Rp = 1.88+-0.03 Rearth over 530-750nm) to revise the density of 55 Cnc e (rho = 6.7+-0.4 g cm-3). Based on these revised properties we have characterized possible interiors of 55 Cnc e using a generalized Bayesian model. We confirm that the planet is likely surrounded by a heavyweight atmosphere, contributing a few percents of the planet radius. While we cannot exclude the presence of a water layer underneath the atmosphere, this scenario is unlikely given the observations of the planet across the entire spectrum and its strong irradiation. Follow-up observations of the system in photometry and in spectroscopy over different time-scales are needed to further investigate the nature and origin of this iconic super-Earth.


title: Structured Bayesian Gaussian process latent variable model: applications to data-driven dimensionality reduction and high-dimensional inversion

abstract: We introduce a methodology for nonlinear inverse problems using a variational Bayesian approach where the unknown quantity is a spatial field. A structured Bayesian Gaussian process latent variable model is used both to construct a low-dimensional generative model of the sample-based stochastic prior as well as a surrogate for the forward evaluation. Its Bayesian formulation captures epistemic uncertainty introduced by the limited number of input and output examples, automatically selects an appropriate dimensionality for the learned latent representation of the data, and rigorously propagates the uncertainty of the data-driven dimensionality reduction of the stochastic space through the forward model surrogate. The structured Gaussian process model explicitly leverages spatial information for an informative generative prior to improve sample efficiency while achieving computational tractability through Kronecker product decompositions of the relevant kernel matrices. Importantly, the Bayesian inversion is carried out by solving a variational optimization problem, replacing traditional computationally-expensive Monte Carlo sampling. The methodology is demonstrated on an elliptic PDE and is shown to return well-calibrated posteriors and is tractable with latent spaces with over 100 dimensions.


title: The VVV Survey RR Lyrae Population in the Galactic Centre Region

abstract: Deep near-IR images from the VVV Survey were used to search for RR Lyrae type ab (RRab) stars within 100' from the Galactic Centre (GC). A sample of 960 RRab stars were discovered. We use the reddening-corrected magnitudes in order to isolate RRab belonging to the GC. The mean period for our RRab sample is $P=0.5446$ days, yielding a mean metallicity of $[Fe/H] = -1.30$ dex and a median distance from the Sun of $D=8.05$. We measure the RRab surface density using the less reddened region sampled here, finding $1000$ RRab/sq deg at a projected Galactocentric distance $R_G=1.6$ deg. This implies a large total mass ($M>10^9 M_\odot$) for the old and metal-poor population contained inside $R_G$. We measure accurate relative proper motions, from which we derive tangential velocity dispersions of $\sigma V_l = 125.0$ and $\sigma V_b = 124.1$ km/s along the Galactic longitude and latitude coordinates, respectively. The fact that these quantities are similar indicate that the bulk rotation of the RRab population is negligible, and implies that this population is supported by velocity dispersion. There are two main conclusions of this study. First, the population as a whole is no different from the outer bulge RRab, predominantly a metal-poor component that is shifted respect the Oosterhoff type I population defined by the globular clusters in the halo. Second, the RRab sample, as representative of the old and metal-poor stellar population in the region, have high velocity dispersions and zero rotation, suggesting a formation via dissipational collapse.


title: On the Trajectory Dependence of Atmospheric Boundary Layer Turbulence Sensing using Small Unmanned Vehicles

abstract: Atmospheric turbulence, especially in the near-surface boundary layer is known to be under-sampled due to the need to capture a wide separation in length and time-scales and limitation in the number of sensors. Over the past decade, the use if Unmanned Aircraft Systems (UAS) technology is approaching ubiquitous proportions for a wide variety of applications, especially with the recent FAA relaxation of flying restrictions. From a geophysical sciences perspective, such technology would allow for sensing of large-scale atmospheric flows, particularly, atmospheric boundary layer (ABL) turbulence, air quality monitoring in urban settings where multitude of small, minimally-invasive and mobile sensors can drastically alter our ability to study such complex phenomena. Currently available observational data of atmospheric boundary layer physics is so sparse and infrequent which significantly limits analysis. With the quantity and resolution of the data that can be measured using a swarm of UAS, three-dimensional reconstruction and deduction of coherent structures in ABL turbulence may be feasible. However, key challenges remain in the form of identifying optimal trajectories to fly the UAS to obtain the relevant quantifications of the turbulence, interpretation of the sensor data from mobile sensors and understanding how representative are the sparse measurements of the overall turbulent boundary layer. This leads to many fundamentally interesting questions if which a few are explored in this article.


title: The twisted tensor product of dg categories and a contractible 2-operad

abstract: It is well-known that the "pre-2-category" $\mathscr{C}at_\mathrm{dg}^\mathrm{coh}(k)$ of small dg categories over a field $k$, with 1-morphisms defined as dg functors, and with 2-morphisms defined as the complexes of coherent natural transformations, fails to be a strict 2-category. In [T2], D.Tamarkin constructed a contractible 2-operad in the sense of M.Batanin [Ba3], acting on $\mathscr{C}at_\mathrm{dg}^\mathrm{coh}(k)$. According to Batanin loc.cit., it is a possible way to define a "weak 2-category". In this paper, we provide a construction of {\it another} contractible 2-operad $\mathcal{O}$, acting on $\mathscr{C}at_\mathrm{dg}^\mathrm{coh}(k)$. Our main tool is the {\it twisted tensor product} of small dg categories, introduced in [Sh3]. We establish a one-side associativity for the twisted tensor product, making $(\mathscr{C}at_\mathrm{dg}^\mathrm{coh}(k),\overset{\sim}{\otimes})$ a skew monoidal category in the sense of [LS], and construct a {\it twisted composition} $\mathscr{C}oh_\mathrm{dg}(D,E)\overset{\sim}{\otimes}\mathscr{C}oh_\mathrm{dg}(C,D)\to\mathscr{C}oh_\mathrm{dg}(C,E)$, and prove some compatibility between these two structures. Taken together, the two structures give rise to a 2-operad $\mathcal{O}$, acting on $\mathscr{C}at_\mathrm{dg}^\mathrm{coh}(k)$. Its contractibility is a consequence of a general result of [Sh3].


title: Probing Higgs boson couplings in $t\bar{t}b\bar{b}$ production at the LHC

abstract: In this paper, we investigate the effects of dimension six operators in the framework of the standard model (SM) effective field theory on using $t\bar{t}b\bar{b}$ channel at the LHC with the center-of-mass energy of 13 TeV. In this analysis, we only assume leptonic top-quarks. Once the events were simulated, the CMS-measured cross-section was used to constrain the corresponding effective Wilson coefficients at $95\%$ confidence level. Considering the SM effective Lagrangian, the $t\bar{t}b\bar{b}$ events in proton-proton collision were very sensitive to the effective Wilson coefficients $ \bar{C}_{g},\tilde{C_{g}},\bar{C}_{uG},\bar{C}_{uW}$. We show that using $t\bar{t}b\bar{b}$ events, we can provide stronger constraints on these Wilson coefficients than the current ones obtained from other channels.


title: Manifold regularization with GANs for semi-supervised learning

abstract: Generative Adversarial Networks are powerful generative models that are able to model the manifold of natural images. We leverage this property to perform manifold regularization by approximating a variant of the Laplacian norm using a Monte Carlo approximation that is easily computed with the GAN. When incorporated into the semi-supervised feature-matching GAN we achieve state-of-the-art results for GAN-based semi-supervised learning on CIFAR-10 and SVHN benchmarks, with a method that is significantly easier to implement than competing methods. We also find that manifold regularization improves the quality of generated images, and is affected by the quality of the GAN used to approximate the regularizer.


title: A Framework for Vehicle Routing Approximation Schemes in Trees

abstract: We develop a general framework for designing polynomial-time approximation schemes (PTASs) for various vehicle routing problems in trees. In these problems, the goal is to optimally route a fleet of vehicles, originating at a depot, to serve a set of clients, subject to various constraints. For example, in Minimum Makespan Vehicle Routing, the number of vehicles is fixed, and the objective is to minimize the longest distance traveled by a single vehicle. Our main insight is that we can often greatly restrict the set of potential solutions without adding too much to the optimal solution cost. This simplification relies on partitioning the tree into clusters such that there exists a near-optimal solution in which every vehicle that visits a given cluster takes on one of a few forms. In particular, only a small number of vehicles serve clients in any given cluster. By using these coarser building blocks, a dynamic programming algorithm can find a near-optimal solution in polynomial time. We show that the framework is flexible enough to give PTASs for many problems, including Minimum Makespan Vehicle Routing, Distance-Constrained Vehicle Routing, Capacitated Vehicle Routing, and School Bus Routing, and can be extended to the multiple depot setting.


title: Search for neutrinos from TXS 0506+056 with the ANTARES telescope

abstract: The results of three different searches for neutrino candidates, associated with the IceCube-170922A event or from the direction of TXS 0506+056, by the ANTARES neutrino telescope are presented. The first search refers to the online follow-up of the IceCube alert; the second is based on the standard time-integrated method employed by the Collaboration to search for point-like neutrino sources; the third uses the information from the IceCube time-dependent analysis reporting a bursting activity centered on December 13, 2014, as input for an ANTARES time-dependent analysis. The online follow-up and the time-dependent analysis yield no events related to the source. The time-integrated study performed over a period from 2007 to 2017 fits 1.03 signal events, which corresponds to a p-value of 3.4% (not considering trial factors). Only for two other astrophysical objects in our candidate list, a smaller p-value had been found. When considering that 107 sources have been investigated, the post-trial p-value for TXS 0506+056 corresponds to 87%.


title: On the zeros of a class of modular functions

abstract: We generalize a number of works on the zeros of certain level 1 modular forms to a class of weakly holomorphic modular functions whose $q$-expansions satisfy \[ f_k(A, \tau) \colon = q^{-k}(1+a(1)q+a(2)q^2+...) + O(q),\] where $a(n)$ are numbers satisfying a certain analytic condition. We show that the zeros of such $f_k(\tau)$ in the fundamental domain of $SL_2(\mathbb{Z})$ lie on $|\tau|=1$ and are transcendental. We recover as a special case earlier work of Witten on extremal "partition" functions $Z_k(\tau)$. These functions were originally conceived as possible generalizations of constructions in three-dimensional quantum gravity.


title: Optical and mechanical design of the extreme AO coronagraphic instrument MagAO-X

abstract: Here we review the current optical mechanical design of MagAO-X. The project is post-PDR and has finished the design phase. The design presented here is the baseline to which all the optics and mechanics have been fabricated. The optical/mechanical performance of this novel extreme AO design will be presented here for the first time. Some highlights of the design are: 1) a floating, but height stabilized, optical table; 2) a Woofer tweeter (2040 actuator BMC MEMS DM) design where the Woofer can be the current f/16 MagAO ASM or, more likely, fed by the facility f/11 static secondary to an ALPAO DM97 woofer; 3) 22 very compact optical mounts that have a novel locking clamp for additional thermal and vibrational stability; 4) A series of four pairs of super-polished off-axis parabolic (OAP) mirrors with a relatively wide FOV by matched OAP clocking; 5) an advanced very broadband (0.5-1.7micron) ADC design; 6) A Pyramid (PWFS), and post-coronagraphic LOWFS NCP wavefront sensor; 7) a vAPP coronagraph for starlight suppression. Currently all the OAPs have just been delivered, and all the rest of the optics are in the lab. Most of the major mechanical parts are in the lab or instrument, and alignment of the optics has occurred for some of the optics (like the PWFS) and most of the mounts. First light should be in 2019A.


title: Micro-Bunched Beam Production at FAST for Narrow Band THz Generation Using a Slit-Mask

abstract: We discuss simulations and experiments on creating micro-bunch beams for generating narrow band THz radiation at the Fermilab Accelerator Science and Technology (FAST) facility. The low-energy electron beamline at FAST consists of a photoinjector-based RF gun, two L-band superconducting accelerating cavities, a chicane, and a beam dump. The electron bunches are lengthened with cavity phases set off-crest for better longitudinal separation and then micro-bunched with a slit-mask installed in the chicane. We carried out the experiments with 30 MeV electron beams and detected signals of the micro-bunching using a skew quadrupole magnet in the chicane. In this paper, the details of micro-bunch beam production, the detection of micro-bunching and comparison with simulations are described.


title: $\eta^{\prime}\to\eta\pi\pi$ decays in unitarized resonance chiral theory

abstract: We study the hadronic $\eta^{\prime}\to\eta\pi\pi$ decays within the framework of $U(3)_{L}\otimes U(3)_{R}$ Chiral Perturbation Theory including resonance states and the complete one-loop corrections. The amplitude is projected in partial waves and unitarized by means of the $N/D$ method resumming both the important $S$-and $D$-wave $\pi\pi$ and the subleading $S$-wave $\pi\eta$ final-state interactions. The participating scalar multiplet mass and coupling strengths are determined from fits to the Dalitz plot experimental data recently released by the A2 collaboration. As a byproduct of our analysis, the associated Dalitz-plot slope parameters are found to be $a=-0.072(7)_{\rm{stat}}(8)_{\rm{syst}}\,, b=-0.052(1)_{\rm{stat}}(2)_{\rm{syst}}\,, d=-0.051(8)_{\rm{stat}}(6)_{\rm{syst}}$, which lie in the ballpark of the current experimental and theoretical determinations.


title: Problems with variable Hilbert space in quantum mechanics

abstract: The general problem is studied on a simple example. A quantum particle in an infinite one-dimensional well potential is considered. Let the boundaries of well changes in a finite time $T$. The standard methods for calculating probability of transition from an initial to the final state are in general inapplicable since the states of different wells belong to different Hilbert spaces. If the final well covers only a part of the initial well (and, possibly, some outer part of the configuration space), the total probability of the transition from any stationary state of the initial well into {\bf all} possible states of the final well is less than 1 at $T\to 0$. If the problem is regularized with a finite-height potential well, this missing probability can be understood as a non-zero probability of transitions into the continuous spectrum, despite the fact that this spectrum disappears at the removal of regularization. This phenomenon ("transition to nowhere'') can result new phenomena in some fundamental problems, in particular at description of earlier Universe. We discuss also how to calculate the probabilities of discussed transitions at final $T$ for some ranges of parameters.


title: MagAO-X: project status and first laboratory results

abstract: MagAO-X is an entirely new "extreme" adaptive optics system for the Magellan Clay 6.5 m telescope, funded by the NSF MRI program starting in Sep 2016. The key science goal of MagAO-X is high-contrast imaging of accreting protoplanets at H$\alpha$. With 2040 actuators operating at up to 3630 Hz, MagAO-X will deliver high Strehls (>70%), high resolution (19 mas), and high contrast ($< 1\times10^{-4}$) at H$\alpha$ (656 nm). We present an overview of the MagAO-X system, review the system design, and discuss the current project status.


title: Design of the MagAO-X Pyramid Wavefront Sensor

abstract: Adaptive optics systems correct atmospheric turbulence in real time. Most adaptive optics systems used routinely correct in the near infrared, at wavelengths greater than 1 micron. MagAO- X is a new extreme adaptive optics (ExAO) instrument that will offer corrections at visible-to- near-IR wavelengths. MagAO-X will achieve Strehl ratios greater than 70% at H-alpha when running the 2040 actuator deformable mirror at 3.6 kHz. A visible pyramid wavefront sensor (PWFS) optimized for sensing at 600-1000 nm wavelengths will provide the high-order wavefront sensing on MagAO- X. We present the optical design and predicted performance of the MagAO-X pyramid wavefront sensor.


title: A Formal Account of Effectiveness Evaluation and Ranking Fusion

abstract: This paper proposes a theoretical framework which models the information provided by retrieval systems in terms of Information Theory. The proposed framework allows to formalize: (i) system effectiveness as an information theoretic similarity between system outputs and human assessments, and (ii) ranking fusion as an information quantity measure. As a result, the proposed effectiveness metric improves popular metrics in terms of formal constraints. In addition, our empirical experiments suggest that it captures quality aspects from traditional metrics, while the reverse is not true. Our work also advances the understanding of theoretical foundations of the empirically known phenomenon of effectiveness increase when combining retrieval system outputs in an unsupervised manner.


title: On the Discrepancy of Random Matrices with Many Columns

abstract: Motivated by the Koml\'os conjecture in combinatorial discrepancy, we study the discrepancy of random matrices with $m$ rows and $n$ independent columns drawn from a bounded lattice random variable. It is known that for $n$ tending to infinity and $m$ fixed, with high probability the $\ell_\infty$-discrepancy is at most twice the $\ell_\infty$-covering radius of the integer span of the support of the random variable. However, the easy argument for the above fact gives no concrete bounds on the failure probability in terms of $n$. We prove that the failure probability is inverse polynomial in $m, n$ and some well-motivated parameters of the random variable. We also obtain the analogous bounds for the discrepancy in arbitrary norms. We apply these results to two random models of interest. For random $t$-sparse matrices, i.e. uniformly random matrices with $t$ ones and $m-t$ zeroes in each column, we show that the $\ell_\infty$-discrepancy is at most $2$ with probability $1 - O(\sqrt{ \log n/n})$ for $n = \Omega(m^3 \log^2 m)$. This improves on a bound proved by Ezra and Lovett (Ezra and Lovett, Approx+Random, 2016) showing that the same is true for $n$ at least $m^t$. For matrices with random unit vector columns, we show that the $\ell_\infty$-discrepancy is $O(\exp(\sqrt{ n/m^3}))$ with probability $1 - O(\sqrt{ \log n/n})$ for $n = \Omega(m^3 \log^2 m).$ Our approach, in the spirit of Kuperberg, Lovett and Peled (G. Kuperberg, S. Lovett and R. Peled, STOC 2012), uses Fourier analysis to prove that for $m \times n$ matrices $M$ with i.i.d. columns, and $m$ sufficiently large, the distribution of $My$ for random $y \in \{-1,1\}^n$ obeys a local limit theorem.


title: GUIMesh: a tool to import STEP geometries into Geant4 via GDML

abstract: Detailed radiation analysis of instruments flown in space is critical to ensure mission safety, often requiring the use of state-of-the-art particle transport simulation tools. Geant4 is one of the most powerful toolkits to simulate the interaction and the passage of particles through matter. This framework however, is not prepared to receive Standard for The Exchange of Product data (STEP) files, the most versatile Computer-Aided Design (CAD) format, as inputs requiring previous conversion to other formats. This, especially when the instruments have complex shapes, and/or a large number of volumes may lead to loss of detail and under or overestimation of the quantities under study. Though several solutions have been proposed to import complex geometries into Geant4, so far, only commercial options are available. In this paper we present a new tool, GUIMesh, that embeds FreeCAD libraries, an open-source CAD editor, to tessellate volumes, and convert them to Geometry Description Markup Language (GDML), a Geant4 readable format, in a straightforward way. Several degrees of freedom are given to the user regarding the mesh and choice of material. Different geometries were tested for p material definition, geometry and navigation errors, successfully validating the method used.


title: Automated Vulnerability Detection in Source Code Using Deep Representation Learning

abstract: Increasing numbers of software vulnerabilities are discovered every year whether they are reported publicly or discovered internally in proprietary code. These vulnerabilities can pose serious risk of exploit and result in system compromise, information leaks, or denial of service. We leveraged the wealth of C and C++ open-source code available to develop a large-scale function-level vulnerability detection system using machine learning. To supplement existing labeled vulnerability datasets, we compiled a vast dataset of millions of open-source functions and labeled it with carefully-selected findings from three different static analyzers that indicate potential exploits. The labeled dataset is available at: https://osf.io/d45bw/. Using these datasets, we developed a fast and scalable vulnerability detection tool based on deep feature representation learning that directly interprets lexed source code. We evaluated our tool on code from both real software packages and the NIST SATE IV benchmark dataset. Our results demonstrate that deep feature representation learning on source code is a promising approach for automated software vulnerability detection.