title: Jordan chains of elliptic partial differential operators and Dirichlet-to-Neumann maps

abstract: Let $\Omega \subset {\bf R}^d$ be a bounded open set with Lipschitz boundary $\Gamma$. It will be shown that the Jordan chains of m-sectorial second-order elliptic partial differential operators with measurable coefficients and (local or non-local) Robin boundary conditions in $L_2(\Omega)$ can be characterized with the help of Jordan chains of the Dirichlet-to-Neumann map and the boundary operator from $H^{1/2}(\Gamma)$ into $H^{-1/2}(\Gamma)$. This result extends the Birman--Schwinger principle in the framework of elliptic operators for the characterization of eigenvalues, eigenfunctions and geometric eigenspaces to the complete set of all generalized eigenfunctions and algebraic eigenspaces.


title: Blocksworld Revisited: Learning and Reasoning to Generate Event-Sequences from Image Pairs

abstract: The process of identifying changes or transformations in a scene along with the ability of reasoning about their causes and effects, is a key aspect of intelligence. In this work we go beyond recent advances in computational perception, and introduce a more challenging task, Image-based Event-Sequencing (IES). In IES, the task is to predict a sequence of actions required to rearrange objects from the configuration in an input source image to the one in the target image. IES also requires systems to possess inductive generalizability. Motivated from evidence in cognitive development, we compile the first IES dataset, the Blocksworld Image Reasoning Dataset (BIRD) which contains images of wooden blocks in different configurations, and the sequence of moves to rearrange one configuration to the other. We first explore the use of existing deep learning architectures and show that these end-to-end methods under-perform in inferring temporal event-sequences and fail at inductive generalization. We then propose a modular two-step approach: Visual Perception followed by Event-Sequencing, and demonstrate improved performance by combining learning and reasoning. Finally, by showing an extension of our approach on natural images, we seek to pave the way for future research on event sequencing for real world scenes.


title: Video-to-Video Translation for Visual Speech Synthesis

abstract: Despite remarkable success in image-to-image translation that celebrates the advancements of generative adversarial networks (GANs), very limited attempts are known for video domain translation. We study the task of video-to-video translation in the context of visual speech generation, where the goal is to transform an input video of any spoken word to an output video of a different word. This is a multi-domain translation, where each word forms a domain of videos uttering this word. Adaptation of the state-of-the-art image-to-image translation model (StarGAN) to this setting falls short with a large vocabulary size. Instead we propose to use character encodings of the words and design a novel character-based GANs architecture for video-to-video translation called Visual Speech GAN (ViSpGAN). We are the first to demonstrate video-to-video translation with a vocabulary of 500 words.


title: Generation of Policy-Level Explanations for Reinforcement Learning

abstract: Though reinforcement learning has greatly benefited from the incorporation of neural networks, the inability to verify the correctness of such systems limits their use. Current work in explainable deep learning focuses on explaining only a single decision in terms of input features, making it unsuitable for explaining a sequence of decisions. To address this need, we introduce Abstracted Policy Graphs, which are Markov chains of abstract states. This representation concisely summarizes a policy so that individual decisions can be explained in the context of expected future transitions. Additionally, we propose a method to generate these Abstracted Policy Graphs for deterministic policies given a learned value function and a set of observed transitions, potentially off-policy transitions used during training. Since no restrictions are placed on how the value function is generated, our method is compatible with many existing reinforcement learning methods. We prove that the worst-case time complexity of our method is quadratic in the number of features and linear in the number of provided transitions, $O(|F|^2 |tr\_samples|)$. By applying our method to a family of domains, we show that our method scales well in practice and produces Abstracted Policy Graphs which reliably capture relationships within these domains.


title: Dirac electron in graphene with magnetic fields arising from first-order intertwining operators

abstract: The behaviour of a Dirac electron in graphene, under magnetic fields which are orthogonal to the layer, is studied. The initial problem is reduced to an equivalent one, where two one-dimensional Schr\"{o}dinger Hamiltonians $H^{\pm}$ are intertwined by a first order differential operator. Special magnetic field are initially chosen, in order that $V^{\pm}$ will be shape invariant exactly solvable potentials. When looking for more general first order operators, intertwining $H^-$ with a non-necessarily shape invariant Hamiltonian, new magnetic fields associated also to analytic solutions will be generated. The iteration of this procedure is as well discussed.


title: Blocks in ASEP with step-Bernoulli initial condition

abstract: This paper extends work by Tracy and Widom on blocks in the asymmetric simple exclusion process (ASEP) to the case of step-Bernoulli initial condition. We consider the probability that a particle at site $x$ is the beginning of a block of $L$ consecutive particles at time $t$ in ASEP with step-Bernoulli initial condition. A Fredholm determinant representation for this probability is derived, and the asymptotics are computed for the KPZ regime.


title: A model of quantum collapse induced by gravity

abstract: We discuss a model where a spontaneous quantum collapse is induced by the gravitational interaction, treated classically. Its dynamics couples the standard wave function of a system with the Bohmian positions of its particles, which are considered as the only source of the gravitational attraction. The collapse is obtained by adding a small imaginary component to the gravitational coupling. It predicts extremely small perturbations of microscopic systems, but very fast collapse of QSMDS (quantum superpositions of macroscopically distinct quantum states) of a solid object, varying as the fifth power of its size. The model does not require adding any dimensional constant to those of standard physics.


title: Examining the Radius Valley: a Machine Learning Approach

abstract: The ''radius valley" is a relative dearth of planets between two potential populations of exoplanets, super-Earths and mini-Neptunes. This feature appears in examining the distribution of planetary radii, but has only ever been characterized on small samples. The valley could be a result of photoevaporation, which has been predicted in numerous theoretical models, or a result of other processes. Here, we investigate the relationship between planetary radius and orbital period through 2-dimensional kernel density estimator and various clustering methods, using all known super-Earths ($R<4.0R_E$). With our larger sample, we confirm the radius valley and characterize it as a power law. Using a variety of methods, we find a range of slopes that are consistent with each other and distinctly negative. We average over these results and find the slope to be $m=-0.319^{+0.088}_{-0.116}$. We repeat our analysis on samples from previous studies. For all methods we use, the resulting line has a negative slope, which is consistent with models of photoevaporation and core-powered mass loss but inconsistent with planets forming in a gas-poor disk.


title: Linear response theory in asymmetric nuclear matter for Skyrme functionals including spin-orbit and tensor terms II: Charge Exchange

abstract: We present the formalism of linear response theory both at zero and finite temperature in the case of asymmetric nuclear matter excited by an isospin flip probe. The particle-hole interaction is derived from a general Skyrme functional that includes spin-orbit and tensor terms. Response functions are obtained by solving a closed algebraic system of equations. Spin strength functions are analyzed for typical values of density, momentum transfer, asymmetry and temperature.


title: Coupled-cluster impurity solvers for dynamical mean-field theory

abstract: We describe the use of coupled-cluster theory as an impurity solver in dynamical mean-field theory (DMFT) and its cluster extensions. We present numerical results at the level of coupled-cluster theory with single and double excitations (CCSD) for the density of states and self-energies of cluster impurity problems in the one- and two-dimensional Hubbard models. Comparison to exact diagonalization shows that CCSD produces accurate density of states and self-energies at a variety of values of $U/t$ and filling fractions. However, the low cost allows for the use of many bath sites, which we define by a discretization of the hybridization directly on the real frequency axis. We observe convergence of dynamical quantities using approximately 30 bath sites per impurity site, with our largest 4-site cluster DMFT calculation using 120 bath sites. We suggest coupled cluster impurity solvers will be attractive in ab initio formulations of dynamical mean-field theory.


title: The compact multiple system HIP 41431

abstract: The nearby (50 pc) K7V dwarf HIP~41431 (EPIC 212096658) is a compact 3-tier hierarchy. Three K7V stars with similar masses, from 0.61 to 0.63 solar, make a triple-lined spectroscopic system where the inner binary with a period of 2.9 days is eclipsing, and the outer companion on a 59-day orbit exerts strong dynamical influence revealed by the eclipse time variation in the {\em Kepler} photometry. Moreover, the centre-of-mass of the triple system moves on a 3.9-year orbit, modulating the proper motion. The mass of the 4-th star is 0.35 solar. The Kepler and ground-based photometry and radial velocities from four different spectrographs are used to adjust the spectro-photodynamical model that accounts for dynamical interaction. The mutual inclination between the two inner orbits is 2\fdg16$\pm$0\fdg11, while the outer orbit is inclined to their common plane by 21\degr$\pm$16\degr. The inner orbit precesses under the influence of both outer orbits, causing observable variation of the eclipse depth. Moreover, the phase of the inner binary is strongly modulated with a 59-day period and its line of apsides precesses. The middle orbit with eccentricity $e=0.28$ also precesses, causing the observed variation of its radial velocity curve. Masses and other parameters of stars in this unique hierarchy are determined. This system is dynamically stable and likely old.


title: A New Distribution on the Simplex with Auto-Encoding Applications

abstract: We construct a new distribution for the simplex using the Kumaraswamy distribution and an ordered stick-breaking process. We explore and develop the theoretical properties of this new distribution and prove that it exhibits symmetry under the same conditions as the well-known Dirichlet. Like the Dirichlet, the new distribution is adept at capturing sparsity but, unlike the Dirichlet, has an exact and closed form reparameterization--making it well suited for deep variational Bayesian modeling. We demonstrate the distribution's utility in a variety of semi-supervised auto-encoding tasks. In all cases, the resulting models achieve competitive performance commensurate with their simplicity, use of explicit probability models, and abstinence from adversarial training.


title: Unitary designs from statistical mechanics in random quantum circuits

abstract: Random quantum circuits are proficient information scramblers and efficient generators of randomness, rapidly approximating moments of the unitary group. We study the convergence of local random quantum circuits to unitary $k$-designs. Employing a statistical mechanical mapping, we give an exact expression of the distance to forming an approximate design as a lattice partition function. In the statistical mechanics model, the approach to randomness has a simple interpretation in terms of domain walls extending through the circuit. We analytically compute the second moment, showing that random circuits acting on $n$ qudits form approximate 2-designs in $O(n)$ depth, as is known. Furthermore, we argue that random circuits form approximate unitary $k$-designs in $O(nk)$ depth and are thus essentially optimal in both $n$ and $k$. We can show this in the limit of large local dimension, but more generally rely on a conjecture about the dominance of certain domain wall configurations.

categories: quant-ph cond-mat.stat-mech cond-mat.str-el hep-th


title: Spatial rogue waves in photorefractive SBN crystals

abstract: We report on the excitation of large-amplitude waves, with a probability of around 1% of total peaks, on a photorefractive SBN crystal by using a simple experimental setup at room temperature. We excite the system using a narrow Gaussian beam and observe different dynamical regimes tailored by the value and time rate of an applied voltage. We identify two main dynamical regimes: a caustic one for energy spreading and a speckling one for peak emergence. Our observations are well described by a two-dimensional Schr\"odinger model with saturable local nonlinearity.


title: Superconvergent interpolatory HDG methods for reaction diffusion equations I: An HDG$_{k}$ method

abstract: In our earlier work [8], we approximated solutions of a general class of scalar parabolic semilinear PDEs by an interpolatory hybridizable discontinuous Galerkin (Interpolatory HDG) method. This method reduces the computational cost compared to standard HDG since the HDG matrices are assembled once before the time integration. Interpolatory HDG also achieves optimal convergence rates; however, we did not observe superconvergence after an element-by-element postprocessing. In this work, we revisit the Interpolatory HDG method for reaction diffusion problems, and use the postprocessed approximate solution to evaluate the nonlinear term. We prove this simple change restores the superconvergence and keeps the computational advantages of the Interpolatory HDG method. We present numerical results to illustrate the convergence theory and the performance of the method.


title: Information-Theoretic Registration with Explicit Reorientation of Diffusion-Weighted Images

abstract: We present an information-theoretic approach to the registration of images with directional information, and especially for diffusion-Weighted Images (DWI), with explicit optimization over the directional scale. We call it Locally Orderless Registration with Directions (LORD). We focus on normalized mutual information as a robust information-theoretic similarity measure for DWI. The framework is an extension of the LOR-DWI density-based hierarchical scale-space model that varies and optimizes the integration, spatial, directional, and intensity scales. As affine transformations are insufficient for inter-subject registration, we extend the model to non-rigid deformations. We illustrate that the proposed model deforms orientation distribution functions (ODFs) correctly and is capable of handling the classic complex challenges in DWI-registrations, such as the registration of fiber-crossings along with kissing, fanning, and interleaving fibers. Our experimental results clearly illustrate a novel promising regularizing effect, that comes from the nonlinear orientation-based cost function. We show the properties of the different image scales and, we show that including orientational information in our model makes the model better at retrieving deformations in contrast to standard scalar-based registration.


title: Hypercyclic bilinear operators on Banach spaces

abstract: We study the dynamics induced by an $m$-linear operator. We answer a question of B\`es and Conejero showing an example of an $m$-linear hypercyclic operator acting on a Banach space. Moreover we prove the existence of $m$-linear hypercyclic operators on arbitrary infinite dimensional separable Banach spaces. We also prove an existence result about symmetric bihypercyclic bilinear operators, answering a question by Grosse-Erdman and Kim.


title: A common origin for dynamically associated near-Earth asteroid pairs

abstract: Though pairs of dynamically associated asteroids in the Main Belt have been identified and studied for over a decade, very few pair systems have been identified in the near-Earth asteroid population. We present data and analysis that supports the existence of two genetically related pairs in near-Earth space. The members of the individual systems, 2015 EE7 -- 2015 FP124 and 2017 SN16 -- 2018 RY7, are found to be of the same spectral taxonomic class, and both pairs are interpreted to have volatile-poor compositions. In conjunction with dynamical arguments, this suggests that these two systems formed via YORP spin-up and/or dissociation of a binary precursor. Backwards orbital integrations suggest a separation age of <10 kyr for the pair 2017 SN16 -- 2018 RY7, making these objects amongst the youngest multiple asteroid systems known to date. A unique separation age was not realized for 2015 EE7 -- 2015 FP124 due to large uncertainties associated with these objects' orbits. Determining the ages of such young pairs is of great value for testing models of space weathering and asteroid spin-state evolution. As the NEO catalog continues to grow with current and future discovery surveys, it is expected that more NEO pairs will be found, thus providing an ideal laboratory for studying time dependent evolutionary processes that are relevant to asteroids throughout the Solar System.


title: The first eigenvalue and eigenfunction of a nonlinear elliptic system

abstract: In this paper, we study the first eigenvalue of a nonlinear elliptic system involving $p$-Laplacian as the differential operator. The principal eigenvalue of the system and the corresponding eigenfunction are investigated both analytically and numerically. An alternative proof to show the simplicity of the first eigenvalue is given. In addition, the upper and lower bounds of the first eigenvalue are provided. Then, a numerical algorithm is developed to approximate the principal eigenvalue. This algorithm generates a decreasing sequence of positive numbers and various examples numerically indicate its convergence. Further, the algorithm is generalized to a class of gradient quasilinear elliptic systems.


title: Discovery of Soft X-ray Pulsations from PSR J1231-1411 using NICER

abstract: We report the discovery of soft X-ray pulsations from the nearby millisecond pulsar PSR J1231$-$1411 using NICER. The pulsed emission is characterized by a broad and asymmetric main pulse and a much fainter secondary interpulse, with a total pulsed count rate of 0.055 c s$^{-1}$ in the 0.35-1.5 keV band. We analyzed Fermi LAT data to update the pulse timing model covering 10 years of data and used that model to coherently combine NICER data over a year of observations. Spectral modeling suggests that the flux is dominated by thermal emission from a hot spot (or spots) on the neutron star surface. The phase relationship between the X-ray pulse and the radio and $\gamma$ rays provides insight into the geometry of the system.


title: On Evaluating the Effectiveness of the HoneyBot: A Case Study

abstract: In recent years, cyber-physical system (CPS) security as applied to robotic systems has become a popular research area. Mainly because robotics systems have traditionally emphasized the completion of a specific objective and lack security oriented design. Our previous work, HoneyBot \cite{celine}, presented the concept and prototype of the first software hybrid interaction honeypot specifically designed for networked robotic systems. The intuition behind HoneyBot was that it would be a remotely accessible robotic system that could simulate unsafe actions and physically perform safe actions to fool attackers. Unassuming attackers would think they were connected to an ordinary robotic system, believing their exploits were being successfully executed. All the while, the HoneyBot is logging all communications and exploits sent to be used for attacker attribution and threat model creation. In this paper, we present findings from the result of a user study performed to evaluate the effectiveness of the HoneyBot framework and architecture as it applies to real robotic systems. The user study consisted of 40 participants, was conducted over the course of several weeks, and drew from a wide range of participants aged between 18-60 with varying level of technical expertise. From the study we found that research subjects could not tell the difference between the simulated sensor values and the real sensor values coming from the HoneyBot, meaning the HoneyBot convincingly spoofed communications.


title: Array-RQMC for option pricing under stochastic volatility models

abstract: Array-RQMC has been proposed as a way to effectively apply randomized quasi-Monte Carlo (RQMC) when simulating a Markov chain over a large number of steps to estimate an expected cost or reward. The method can be very effective when the state of the chain has low dimension. For pricing an Asian option under an ordinary geometric Brownian motion model, for example, Array-RQMC reduces the variance by huge factors. In this paper, we show how to apply this method and we study its effectiveness in case the underlying process has stochastic volatility. We show that Array-RQMC can also work very well for these models, even if it requires RQMC points in larger dimension. We examine in particular the variance-gamma, Heston, and Ornstein-Uhlenbeck stochastic volatility models, and we provide numerical results.


title: Putting Strong Linearizability in Context: Preserving Hyperproperties in Programs that Use Concurrent Objects

abstract: It has been observed that linearizability, the prevalent consistency condition for implementing concurrent objects, does not preserve some probability distributions. A stronger condition, called strong linearizability has been proposed, but its study has been somewhat ad-hoc. This paper investigates strong linearizability by casting it in the context of observational refinement of objects. We present a strengthening of observational refinement, which generalizes strong linearizability, obtaining several important implications. When a concrete concurrent object refining another, more abstract object - often sequential - the correctness of a program employing the concrete object can be verified by considering its behaviors when using the more abstract object. This means that trace properties of a program using the concrete object can be proved by considering the program with the abstract object. This, however, does not hold for hyperproperties, including many security properties and probability distributions of events. We define strong observational refinement, a strengthening of refinement that preserves hyperproperties, and prove that it is equivalent to the existence of forward simulations. We show that strong observational refinement generalizes strong linearizability. This implies that strong linearizability is also equivalent to forward simulation, and shows that strongly linearizable implementations can be composed both horizontally (i.e., locality) and vertically (i.e., with instantiation). For situations where strongly linearizable implementations do not exist (or are less efficient), we argue that reasoning about hyperproperties of programs can be simplified by strong observational refinement of abstract objects that are not necessarily sequential.


title: Zero temperature ordering dynamics in two dimensional BNNNI model

abstract: We investigate the dynamics of a two dimensional bi-axial next nearest neighbour Ising (BNNNI) model following a quench to zero temperature. The Hamiltonian is given by $H = -J_0\sum_{i,j=1}^L [(S_{i,j}S_{i+1,j}+S_{i,j}S_{i,j+1}) -\kappa (S_{i,j}S_{i+2,j} + S_{i,j}S_{i,j+2})]$ . For $\kappa <1$, the system does not reach the equilibrium ground state and keep evolving in active states for ever. For $\kappa \geq 1$, though the system reaches a final state, but it do not reach the ground state always and freezes to a striped state with a finite probability like two dimensional ferromagnetic Ising model and ANNNI model. The overall dynamical behaviour for $\kappa > 1$ and $\kappa =1$ is quite different. The residual energy decays in a power law for both $ \kappa >1$ and $\kappa =1$ from which the dynamical exponent $z$ have been estimated. The persistence probability shows algebraic decay for $\kappa > 1$ with an exponent $\theta = 0.22 \pm 0.002$ while the dynamical exponent for ordering $z=2.33\pm 0.01$. For $\kappa =1$, the system belongs to a completely different dynamical class with $\theta = 0.332 \pm 0.002$ and $z=2.47\pm 0.04$. We have computed the freezing probability for different values of $\kappa$. We have also studied the decay of autocorrelation function with time for different regime of $\kappa$ values. The results have been compared with that of the two dimensional ANNNI model.


title: Automatic Ambiguity Detection

abstract: Most work on sense disambiguation presumes that one knows beforehand -- e.g. from a thesaurus -- a set of polysemous terms. But published lists invariably give only partial coverage. For example, the English word tan has several obvious senses, but one may overlook the abbreviation for tangent. In this paper, we present an algorithm for identifying interesting polysemous terms and measuring their degree of polysemy, given an unlabeled corpus. The algorithm involves: (i) collecting all terms within a k-term window of the target term; (ii) computing the inter-term distances of the contextual terms, and reducing the multi-dimensional distance space to two dimensions using standard methods; (iii) converting the two-dimensional representation into radial coordinates and using isotonic/antitonic regression to compute the degree to which the distribution deviates from a single-peak model. The amount of deviation is the proposed polysemy index


title: Emerging trends in metallicity and lithium properties of debris disc stars

abstract: Dwarf stars with debris discs and planets appear to be excellent laboratories to study the core accretion theory of planets formation. These systems are however, insufficiently studied. In this paper we present the main metallicity and lithium abundance properties of these stars together with stars with only debris discs and stars with only planets. Stars without detected planets nor discs are also considered. The analysed sample is formed by main-sequence FGK field single stars. Apart from the basic stellar parameters, we include the use of dusty discs masses. The main results show for the first time that the dust mass of debris disc stars with planets correlate with metallicity. We confirm that these disc dust masses are related to their central stellar masses. Separately, the masses of stars and those of planets also correlate with metallicity. We conclude that two conditions are necessary to form giant planets: to have a sufficient metallicity and also a sufficient protoplanetary mass of gas and dust. The debris discs masses of stars without giant planets do not correlate with metallicity, because they do not fulfil these two conditions. Concerning lithium, by adopting a stellar model for lithium depletion based on a strong interaction between the star and a protoplanetary disc, we found that in agreement with the model predictions, observations indicate that the main lithium depletion occurs during this initial protoplanetary evolution stage. We show that the ultimately lithium depletion is independent of the presence or absence of planets and appears to be only age dependent.


title: On the identification of a nonlinear term in a reaction-diffusion equation

abstract: Reaction-diffusion equations are one of the most common partial differential equations used to model physical phenomenon. They arise as the combination of two physical processes: a driving force $f(u)$ that depends on the state variable $u$ and a diffusive mechanism that spreads this effect over a spatial domain. The canonical form is $u_t - \triangle u = f(u)$. Application areas include chemical processes, heat flow models and population dynamics. The direct or forwards problem for such equations is now very well-developed and understood. However, our interest lies in the inverse problem of recovering the reaction term $f(u)$ not just at the level of determining a few parameters in a known functional form, but recovering the complete functional form itself. To achieve this we set up the usual paradigm for the parabolic equation where $u$ is subject to both given initial and boundary data, then prescribe overposed data consisting of the solution at a later time $T$. For example, in the case of a population model this amounts to census data at a fixed time. Our approach will be two-fold.First we will transform the inverse problem into an equivalent nonlinear mapping from which we seek a fixed point. We will be able to prove important features of this map such as a self-mapping property and give conditions under which it is contractive. Second, we consider the direct map from $f$ through the partial differential operator to the overposed data. We will investigate Newton schemes for this case. In recent decades various anomalous processes have been used to generalize classical Brownian diffusion. Amongst the most popular is one that replaces the usual time derivative by a fractional one of order $\alpha \leq 1$. We will also include this model in our analysis. The final section of the paper shows numerical reconstructions that demonstrate the viability of the suggested approaches.


title: Quantum Circuits for Dynamic Runtime Assertions in Quantum Computation

abstract: In this paper, we propose quantum circuits for runtime assertions, which can be used for both software debugging and error detection. Runtime assertion is challenging in quantum computing for two key reasons. First, a quantum bit (qubit) cannot be copied, which is known as the non-cloning theorem. Second, when a qubit is measured, its superposition state collapses into a classical state, losing the inherent parallel information. In this paper, we overcome these challenges with runtime computation through ancilla qubits, which are used to indirectly collect the information of the qubits of interest. We design quantum circuits to assert classical states, entanglement, and superposition states.


title: SEMA: an Extended Semantic Evaluation Metric for AMR

abstract: Abstract Meaning Representation (AMR) is a recently designed semantic representation language intended to capture the meaning of a sentence, which may be represented as a single-rooted directed acyclic graph with labeled nodes and edges. The automatic evaluation of this structure plays an important role in the development of better systems, as well as for semantic annotation. Despite there is one available metric, smatch, it has some drawbacks. For instance, smatch creates a self-relation on the root of the graph, has weights for different error types, and does not take into account the dependence of the elements in the AMR structure. With these drawbacks, smatch masks several problems of the AMR parsers and distorts the evaluation of the AMRs. In view of this, in this paper, we introduce an extended metric to evaluate AMR parsers, which deals with the drawbacks of the smatch metric. Finally, we compare both metrics, using four well-known AMR parsers, and we argue that our metric is more refined, robust, fairer, and faster than smatch.


title: Extreme sensitivity and climate tipping points

abstract: A climate state close to a tipping point will have a degenerate linear response to perturbations, which can be associated with extreme values of the equilibrium climate sensitivity (ECS). In this paper we contrast linearized (`instantaneous') with fully nonlinear geometric (`two-point') notions of ECS, in both presence and absence of tipping points. For a stochastic energy balance model of the global mean surface temperature with two stable regimes, we confirm that tipping events cause the appearance of extremes in both notions of ECS. Moreover, multiple regimes with different mean sensitivities are visible in the two-point ECS. We confirm some of our findings in a physics-based multi-box model of the climate system.


title: Guarantees for Sound Abstractions for Generalized Planning (Extended Paper)

abstract: Generalized planning is about finding plans that solve collections of planning instances, often infinite collections, rather than single instances. Recently it has been shown how to reduce the planning problem for generalized planning to the planning problem for a qualitative numerical problem; the latter being a reformulation that simultaneously captures all the instances in the collection. An important thread of research thus consists in finding such reformulations, or abstractions, automatically. A recent proposal learns the abstractions inductively from a finite and small sample of transitions from instances in the collection. However, as in all inductive processes, the learned abstraction is not guaranteed to be correct for the whole collection. In this work we address this limitation by performing an analysis of the abstraction with respect to the collection, and show how to obtain formal guarantees for generalization. These guarantees, in the form of first-order formulas, may be used to 1) define subcollections of instances on which the abstraction is guaranteed to be sound, 2) obtain necessary conditions for generalization under certain assumptions, and 3) do automated synthesis of complex invariants for planning problems. Our framework is general, it can be extended or combined with other approaches, and it has applications that go beyond generalized planning.


title: Second law of thermodynamics for batteries with vacuum state

abstract: In stochastic thermodynamics work is a random variable whose average is bounded by the change in the free energy of the system. In most treatments, however, the work reservoir that absorbs this change is either tacitly assumed or modelled using unphysical systems with unbounded Hamiltonians (i.e. the ideal weight). In this work we describe the consequences of introducing the ground state of the battery and hence -- of breaking its translational symmetry. The most striking consequence of this shift is the fact that the Jarzynski identity is replaced by a family of inequalities. Using these inequalities we obtain corrections to the second law of thermodynamics which vanish exponentially with the distance of the initial state of the battery to the bottom of its spectrum. Finally, we study an exemplary thermal operation which realizes the approximate Landauer erasure and demonstrate the consequences which arise when the ground state of the battery is explicitly introduced. In particular, we show that occupation of the vacuum state of any physical battery sets a lower bound on fluctuations of work, while batteries without vacuum state allow for fluctuation-free erasure.


title: Quantum Cluster Quasicrystals

abstract: Quasicrystals remain among the most intriguing materials in physics and chemistry. Their structure results in many unusual properties including anomalously low friction as well as poor electrical and thermal conductivity but it also supports superconductivity, which shows that quantum effects in quasicrystals can be quite unique. Here we theoretically study superfluidity in a model quantum cluster quasicrystal. Using path-integral Monte Carlo simulations, we explore a 2D ensemble of bosons with the Lifshitz-Petrich-Gaussian pair potential, finding that moderate quantum fluctuations do not destroy the dodecagonal quasicrystalline order. This quasicrystal is characterized by a small yet finite superfluidity, demonstrating that particle clustering combined with the local cogwheel structure can underpin superfluidity even in the almost classical regime. This type of distributed superfluidity may also be expected in certain open crystalline lattices. Large quantum fluctuations are shown to induce transitions to cluster solids, supersolids and superfluids, which we characterize fully quantum-mechanically.


title: Inverse approximation and GBS of bivariate Kantorovich type sampling series

abstract: In this paper, we derive an inverse result for bivariate Kantorovich type sampling series for the space of all continuous functions with upto second order partial derivatives are continuous and bounded on $R^2.$ Further, we prove the rate of approximation in the Bogel space of continuous functions for the GBS (Generalized Boolean Sum) of these operators. Finally, we give some examples for the kernel to which the theory can be applied


title: Regimes of classical simulability for noisy Gaussian boson sampling

abstract: As a promising candidate for exhibiting quantum computational supremacy, Gaussian Boson Sampling (GBS) is designed to exploit the ease of experimental preparation of Gaussian states. However, sufficiently large and inevitable experimental noise might render GBS classically simulable. In this work, we formalize this intuition by establishing a sufficient condition for approximate polynomial-time classical simulation of noisy GBS --- in the form of an inequality between the input squeezing parameter, the overall transmission rate and the quality of photon detectors. Our result serves as a non-classicality test that must be passed by any quantum computationalsupremacy demonstration based on GBS. We show that, for most linear-optical architectures, where photon loss increases exponentially with the circuit depth, noisy GBS loses its quantum advantage in the asymptotic limit. Our results thus delineate intermediate-sized regimes where GBS devices might considerably outperform classical computers for modest noise levels. Finally, we find that increasing the amount of input squeezing is helpful to evade our classical simulation algorithm, which suggests a potential route to mitigate photon loss.


title: Ionic activity in concentrated electrolytes: solvent structure effect revisited

abstract: We revisit the role of the local solvent structure on the activity coefficient of electrolytes with a general non-local dielectric function approach. We treat the concentrated electrolyte as a dielectric medium and suggest an interpolated formula for the dielectric response. The pure water limit is calibrated based on MD simulations and experimental data. Solving our model around a central ion, we find strong over-screening and oscillations in the potential, which are absent in the standard "primitive model" predictions. We obtain mathematically tractable closed-form expressions for the activity coefficients and show reasonable agreement with experimental data.


title: Verification and Control for Finite-Time Safety of Stochastic Systems via Barrier Functions

abstract: This paper studies the problem of enforcing safety of a stochastic dynamical system over a finite time horizon. We use stochastic barrier functions as a means to quantify the probability that a system exits a given safe region of the state space in finite time. A barrier certificate condition that bounds the infinitesimal generator of the system, and hence bounds the expected value of the barrier function over the time horizon, is recast as a sum-of-squares optimization problem for efficient numerical computation. Unlike prior works, the proposed certificate condition includes a state-dependent bound on the infinitesimal generator, allowing for tighter probability bounds. Moreover, for stochastic systems for which the drift dynamics are affine-in-control, we propose a method for synthesizing polynomial state feedback controllers that achieve a specified probability of safety. Two case studies are presented that benchmark and illustrate the performance of our method.


title: Stationary states for underdamped anharmonic oscillators driven by Cauchy noise

abstract: Using methods of stochastic dynamics, we have studied stationary states in the underdamped anharmonic stochastic oscillators driven by Cauchy noise. Shape of stationary states depend both on the potential type and the damping. If the damping is strong enough, for potential wells which in the overdamped regime produce multimodal stationary states, stationary states in the underdamped regime can be multimodal with the same number of modes like in the overdamped regime. For the parabolic potential, the stationary density is always unimodal and it is given by the two dimensional $\alpha$-stable density. For the mixture of quartic and parabolic single-well potentials the stationary density can be bimodal. Nevertheless, the parabolic addition, which is strong enough, can destroy bimodlity of the stationary state.


title: Probabilistic Category-Level Pose Estimation via Segmentation and Predicted-Shape Priors

abstract: We introduce a new method for category-level pose estimation which produces a distribution over predicted poses by integrating 3D shape estimates from a generative object model with segmentation information. Given an input depth-image of an object, our variable-time method uses a mixture density network architecture to produce a multi-modal distribution over 3DOF poses; this distribution is then combined with a prior probability encouraging silhouette agreement between the observed input and predicted object pose. Our approach significantly outperforms the current state-of-the-art in category-level 3DOF pose estimation---which outputs a point estimate and does not explicitly incorporate shape and segmentation information---as measured on the Pix3D and ShapeNet datasets.


title: Non-normal Recurrent Neural Network (nnRNN): learning long time dependencies while improving expressivity with transient dynamics

abstract: A recent strategy to circumvent the exploding and vanishing gradient problem in RNNs, and to allow the stable propagation of signals over long time scales, is to constrain recurrent connectivity matrices to be orthogonal or unitary. This ensures eigenvalues with unit norm and thus stable dynamics and training. However this comes at the cost of reduced expressivity due to the limited variety of orthogonal transformations. We propose a novel connectivity structure based on the Schur decomposition and a splitting of the Schur form into normal and non-normal parts. This allows to parametrize matrices with unit-norm eigenspectra without orthogonality constraints on eigenbases. The resulting architecture ensures access to a larger space of spectrally constrained matrices, of which orthogonal matrices are a subset. This crucial difference retains the stability advantages and training speed of orthogonal RNNs while enhancing expressivity, especially on tasks that require computations over ongoing input sequences.


title: Semi-Supervised Learning, Causality and the Conditional Cluster Assumption

abstract: While the success of semi-supervised learning (SSL) is still not fully understood, Sch\"olkopf et al. (2012) have established a link to the principle of independent causal mechanisms. They conclude that SSL should be impossible when predicting a target variable from its causes, but possible when predicting it from its effects. Since both these cases are somewhat restrictive, we extend their work by considering classification using cause and effect features at the same time, such as predicting disease from both risk factors and symptoms. While standard SSL exploits information contained in the marginal distribution of all inputs (to improve the estimate of the conditional distribution of the target given inputs), we argue that in our more general setting we should use information in the conditional distribution of effect features given causal features. We explore how this insight generalises the previous understanding, and how it relates to and can be exploited algorithmically for SSL.


title: Memory Integrity of CNNs for Cross-Dataset Facial Expression Recognition

abstract: Facial expression recognition is a major problem in the domain of artificial intelligence. One of the best ways to solve this problem is the use of convolutional neural networks (CNNs). However, a large amount of data is required to train properly these networks but most of the datasets available for facial expression recognition are relatively small. A common way to circumvent the lack of data is to use CNNs trained on large datasets of different domains and fine-tuning the layers of such networks to the target domain. However, the fine-tuning process does not preserve the memory integrity as CNNs have the tendency to forget patterns they have learned. In this paper, we evaluate different strategies of fine-tuning a CNN with the aim of assessing the memory integrity of such strategies in a cross-dataset scenario. A CNN pre-trained on a source dataset is used as the baseline and four adaptation strategies have been evaluated: fine-tuning its fully connected layers; fine-tuning its last convolutional layer and its fully connected layers; retraining the CNN on a target dataset; and the fusion of the source and target datasets and retraining the CNN. Experimental results on four datasets have shown that the fusion of the source and the target datasets provides the best trade-off between accuracy and memory integrity.


title: EasySched: a multi-agent architecture for the predictive and reactive scheduling of Industry 4.0 production systems based on the available renewable energy

abstract: Industry 4.0 is concerned with sustainable development constraints. In this context, we propose a multi-agent architecture, named EasySched, aiming at elaborating predictive and reactive scheduling as the result of a coordination between systems producing goods and systems producing renewable energy. The validation of this architecture is original, and was conducted in a completely and physically distributed way, using networked embedded systems. This validation was done on a series of instances inspired by the literature. The results showed that EasySched succeeds in adapting the production of goods according to the available renewable energy


title: Attractive Kronig-Penney Band Structures and Wave Functions

abstract: The repulsive-potential Kronig-Penney (KP) model for a one-dimensional band structure is well known. However, real metals contain positively-charged ions resulting in attractive potential wells seen by the metallic electrons. Here we consider the latter case in detail. The square-well version of the KP model is considered first, for which the band structure and wave functions for different potential-well depths are derived. Then an extended treatment of the attractive Dirac-comb version of the KP model is presented. For the nearly-free-electron case, the band structure exhibits a negative-energy band in addition to positive-energy bands. The wave functions, electron densities of states, effective masses, and group velocities are derived for the positive-energy band states. The wave functions of the negative-energy band states are also calculated and found to be quite different from the sinusoidal wave functions for the positive-energy band and band-gap states. High-degeneracy bound states are found at negative energies and their wave functions are derived.


title: Zero mode protection at particle-hole symmetry: a geometric interpretation

abstract: The properties of zero modes in particle-hole symmetric systems are analyzed in the presence of strong random scattering by a disordered environment. The study is based on the calculation of the time-averaged density distribution on a lattice. In particular, a flat distribution is found for strong random scattering. This result is compared with a decaying distribution for weak random scattering by an analysis of the scattering paths. In the calculation we consider the invariant measure of the average two-particle Green's function, which is related to lattice-covering self-avoiding (LCSA) strings. In particular, strong scattering is associated with LCSA loops, whereas weaker scattering is associated with open LCSA strings. Our results are a generalization of the delocalized state observed at the band center of a one-dimensional tight-binding model with random hopping by Dyson in 1953.


title: Riemann solver with internal reconstruction (RSIR) for compressible single-phase and non-equilibrium two-phase flows

abstract: A new Riemann solver is built to address numerical resolution of complex flow models. The research direction is closely linked to a variant of the Baer and Nunziato (1986) model developed in Saurel et al. (2017a). This recent model provides a link between the Marble (1963) model for two-phase dilute suspensions and dense mixtures. As in the Marble model, Saurel et al. system is weakly hyperbolic with the same 4 characteristic waves, while the system involves 7 partial differential equations. It poses serious theoretical and practical issues to built simple and accurate flow solver. To overcome related difficulties the Riemann solver of Linde (2002) is revisited. The method is first examined in the simplified context of compressible Euler equations. Physical considerations are introduced in the solver improving robustness and accuracy of the Linde method. With these modifications the flow solver appears as accurate as the HLLC solver of Toro et al. (1994). Second the two-phase flow model is considered. A locally conservative formulation is built and validated removing issues related to non-conservative terms. However, two extra major issues appear from numerical experiments: The solution appears not self-similar and multiple contact waves appear in the dispersed phase. Building HLLC-type or any other solver appears consequently challenging. The modified Linde (2002) method is thus examined for the considered flow model. Some basic properties of the equations are used, such as shock relations of the dispersed phase and jump conditions across the contact wave. Thanks to these ingredients the new Riemann solver with internal reconstruction (RSIR), modification of the Linde method, handles stationary volume fraction discontinuities, presents low dissipation for transport waves and handles shocks and expansion waves accurately. It is validated on various test problems showing method's accuracy and versatility for complex flow models. Its capabilities are illustrated on a difficult two-phase flow instability problem, unresolved before.


title: A novel application of solid state detectors for high precision, low systematics measurement of beta decay energy spectra of interest for neutrino and nuclear physics

abstract: This project is focused on the development of a novel strategy for the precise measurement of beta decay energy spectra. The exact determination of beta spectra has wide implications in the particle and nuclear physics fields: from sterile neutrino searches \cite{Adhikari 2016} to validation of nuclear models of interest for double beta decay searches \cite{Suhonen 2017}, to reactor neutrino experiments. The experimental strategy is focused on the mitigation of the systematic uncertainty in the determination of the spectral shape related to the energy response of the detector. The plan is to use Silicon Drift Detectors (SDDs), exploiting their excellent energy resolution and response uniformity, the thin dead layer and the fast signal that allows high rate operations. A complete modeling of SDDs response to the interaction of beta particles via numerical simulations, validated with dedicated measurements, aims at identifying and canceling the detector-related systematics. Similar attention will be devoted to the integration of beta radioactive sources that preserve the energy information carried by the emitted electrons, and to a veto system with a full solid angle coverage to intercept any escaping particle potentially carrying a fraction of the original electron energy.


title: Reversed electromagnetic Vavilov-\v{C}erenkov radiation in naturally existing magnetoelectric media

abstract: We consider two semi-infinite magnetoelecteric media separated by a planar interface whose electromagnetic response is described by axion-electrodynamics. The time dependent Green's function characterizing this geometry is obtained by a method which can be directly generalized to cylindrical and spherical configurations of two magnetoelectrics separated by an interface. We establish the far-field approximation of the Green's function and apply these results to the case of a charged particle moving from one medium to the other at a high constant velocity perpendicular to the interface. From the resulting angular distribution of the radiated energy per unit frequency we provide theoretical evidence for the emergence of reversed Vavilov-\v{C}erenkov radiation in naturally existing magnetoelectric media. In the case where one of the magnetolectrics is a 3D topological insulator, ${\rm TlBiSe}_2$ for example, located in front of a regular insulator, we estimate that an average forward Vavilov-\v{C}erenkov radiation with frequency $\sim 2.5 \,\, {\rm eV}$ ($\sim 500\,\, {\rm nm}$) will produce a highly suppressed reversed Vavilov-\v{C}erenkov radiation which can be characterized by an effective frequency in the range of $\sim (4\times 10^{-3}-0.5) \,\, {\rm meV}$. However, this value compares favorably with recent measurements in left-handed metamaterials yielding reversed Vavilov-\v{C}erenkov radiation with frequencies of the order of $(1.2-3.9)\times 10^{-2}\,\, {\rm meV}$.


title: Weak and Compact Radio Emission in Early High-Mass Star Forming Regions: II. The Nature of the Radio Sources

abstract: In this study we analyze 70 radio continuum sources associated with dust clumps and considered to be candidates for the earliest stages of high-mass star formation. The detection of these sources was reported by Rosero et al. (2016), who found most of them to show weak (${\scriptstyle <}$1 mJy) and compact (${\scriptstyle <}\,$0.6$^{\prime \prime}$) radio emission. Herein, we used the observed parameters of these sources to investigate the origin of the radio continuum emission. We found that at least $\sim 30\%$ of these radio detections are most likely ionized jets associated with high-mass protostars, but for the most compact sources we cannot discard the scenario that they represent pressure-confined HII regions. This result is highly relevant for recent theoretical models based on core accretion that predict the first stages of ionization from high-mass stars to be in the form of jets. Additionally, we found that properties such as the radio luminosity as a function of the bolometric luminosity of ionized jets from low and high-mass stars are extremely well-correlated. Our data improve upon previous studies by providing further evidence of a common origin for jets independently of luminosity.


title: Efficient Amortised Bayesian Inference for Hierarchical and Nonlinear Dynamical Systems

abstract: We introduce a flexible, scalable Bayesian inference framework for nonlinear dynamical systems characterised by distinct and hierarchical variability at the individual, group, and population levels. Our model class is a generalisation of nonlinear mixed-effects (NLME) dynamical systems, the statistical workhorse for many experimental sciences. We cast parameter inference as stochastic optimisation of an end-to-end differentiable, block-conditional variational autoencoder. We specify the dynamics of the data-generating process as an ordinary differential equation (ODE) such that both the ODE and its solver are fully differentiable. This model class is highly flexible: the ODE right-hand sides can be a mixture of user-prescribed or "white-box" sub-components and neural network or "black-box" sub-components. Using stochastic optimisation, our amortised inference algorithm could seamlessly scale up to massive data collection pipelines (common in labs with robotic automation). Finally, our framework supports interpretability with respect to the underlying dynamics, as well as predictive generalization to unseen combinations of group components (also called "zero-shot" learning). We empirically validate our method by predicting the dynamic behaviour of bacteria that were genetically engineered to function as biosensors. Our implementation of the framework, the dataset, and all code to reproduce the experimental results is available at https://www.github.com/Microsoft/vi-hds .


title: Approximate Guarantees for Dictionary Learning

abstract: In the dictionary learning (or sparse coding) problem, we are given a collection of signals (vectors in $\mathbb{R}^d$), and the goal is to find a "basis" in which the signals have a sparse (approximate) representation. The problem has received a lot of attention in signal processing, learning, and theoretical computer science. The problem is formalized as factorizing a matrix $X (d \times n)$ (whose columns are the signals) as $X = AY$, where $A$ has a prescribed number $m$ of columns (typically $m \ll n$), and $Y$ has columns that are $k$-sparse (typically $k \ll d$). Most of the known theoretical results involve assuming that the columns of the unknown $A$ have certain incoherence properties, and that the coefficient matrix $Y$ has random (or partly random) structure. The goal of our work is to understand what can be said in the absence of such assumptions. Can we still find $A$ and $Y$ such that $X \approx AY$? We show that this is possible, if we allow violating the bounds on $m$ and $k$ by appropriate factors that depend on $k$ and the desired approximation. Our results rely on an algorithm for what we call the threshold correlation problem, which turns out to be related to hypercontractive norms of matrices. We also show that our algorithmic ideas apply to a setting in which some of the columns of $X$ are outliers, thus giving similar guarantees even in this challenging setting.


title: Instanton sheaves and representations of quivers

abstract: We study the moduli space of rank 2 instanton sheaves on $\p3$ in terms of representations of a quiver consisting of 3 vertices and 4 arrows between two pairs of vertices. Aiming at an alternative compactification for the moduli space of instanton sheaves, we show that for each rank 2 instanton sheaf, there is a stability parameter $\theta$ for which the corresponding quiver representation is $\theta$-stable (in the sense of King), and that the space of stability parameters has a non trivial wall-and-chamber decomposition. Looking more closely at instantons of low charge, we prove that there are stability parameters with respect to which every representation corresponding to a rank 2 instanton sheaf of charge 2 is stable, and provide a complete description of the wall-and-chamber decomposition for representation corresponding to a rank 2 instanton sheaf of charge 1.


title: Morphology and kinetics of asphalt binder microstructure at gas, liquid, and solid interfaces

abstract: We combined optical and atomic force microscopy to observe morphology and kinetics of microstructures that formed at free surfaces of unmodified pavement-grade 64-22 asphalt binders upon cooling from 150 $^{\circ}$C to room temperature (RT) at 5 $^{\circ}$C/min, and changes in these microstructures when the surface was terminated with a transparent solid (glass) or liquid (glycerol) over-layer. The main findings are: (1) At free binder surfaces, wrinkled microstructures started to form near the wax crystallization temperature ($\sim$45 $^{\circ}$C), then grew to $\sim$5 $\mu$m diameter, $\sim$25 nm wrinkle amplitude and 10-30$\%$ surface area coverage upon cooling to RT, where they persisted indefinitely without observable change in shape or density. (2) Glycerol coverage of the binder surface during cooling reduced wrinkled area and wrinkle amplitude three-fold compared to free binder surfaces upon initial cooling to RT; continued glycerol coverage at RT eliminated most surface microstructures within $\sim$4 hours. (3) No surface microstructures were observed to form at binder surfaces covered with glass. (4) Sub-micron bulk microstructures were observed by near-infrared microscopy beneath the surfaces of all binder samples, with size, shape and density independent of surface coverage. No tendency of such structures to float to the top or sink to the bottom of mm-thick samples was observed. (5) We attribute the dependence of surface wrinkling on surface coverage to variation in interface tension, based on a thin-film continuum mechanics model.


title: Quantum entropic self-localization with ultracold fermions

abstract: We study a driven, spin-orbit coupled fermionic system in a lattice at the resonant regime where the drive frequency equals the Hubbard repulsion, for which non-trivial constrained dynamics emerge at fast timescales. An effective density-dependent tunneling model is derived, and examined in the sparse filling regime in 1D. The system exhibits entropic self-localization, where while even numbers of atoms propagate ballistically, odd numbers form localized bound states induced by an effective attraction from a higher configurational entropy. These phenomena occur in the strong coupling limit where interactions only impose a constraint with no explicit Hamiltonian term. We show how the constrained dynamics lead to quantum few-body scars and map to an Anderson impurity model with an additional intriguing feature of non-reciprocal scattering. Connections to many-body scars and localization are also discussed.

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


title: On Linear Programming for Constrained and Unconstrained Average-Cost Markov Decision Processes with Countable Action Spaces and Strictly Unbounded Costs

abstract: We consider the linear programming approach for constrained and unconstrained Markov decision processes (MDPs) under the long-run average cost criterion, where the class of MDPs in our study have Borel state spaces and discrete countable action spaces. Under a strict unboundedness condition on the one-stage costs and a recently introduced majorization condition on the state transition stochastic kernel, we study infinite-dimensional linear programs for the average-cost MDPs and prove the absence of a duality gap and other optimality results. Our results do not require a lower-semicontinuous MDP model. Thus, they can be applied to countable action space MDPs where the dynamics and one-stage costs are discontinuous in the state variable. Our proofs make use of the continuity property of Borel measurable functions asserted by Lusin's theorem.


title: Planning with State Abstractions for Non-Markovian Task Specifications

abstract: Often times, we specify tasks for a robot using temporal language that can also span different levels of abstraction. The example command ``go to the kitchen before going to the second floor'' contains spatial abstraction, given that ``floor'' consists of individual rooms that can also be referred to in isolation ("kitchen", for example). There is also a temporal ordering of events, defined by the word "before". Previous works have used Linear Temporal Logic (LTL) to interpret temporal language (such as "before"), and Abstract Markov Decision Processes (AMDPs) to interpret hierarchical abstractions (such as "kitchen" and "second floor"), separately. To handle both types of commands at once, we introduce the Abstract Product Markov Decision Process (AP-MDP), a novel approach capable of representing non-Markovian reward functions at different levels of abstractions. The AP-MDP framework translates LTL into its corresponding automata, creates a product Markov Decision Process (MDP) of the LTL specification and the environment MDP, and decomposes the problem into subproblems to enable efficient planning with abstractions. AP-MDP performs faster than a non-hierarchical method of solving LTL problems in over 95% of tasks, and this number only increases as the size of the environment domain increases. We also present a neural sequence-to-sequence model trained to translate language commands into LTL expression, and a new corpus of non-Markovian language commands spanning different levels of abstraction. We test our framework with the collected language commands on a drone, demonstrating that our approach enables a robot to efficiently solve temporal commands at different levels of abstraction.


title: Interpolation over ZZ and torsion in class groups

abstract: We prove an interpolation result for homogeneous polynomials over the integers, or more generally for PIDs with finite residue fields. Previous proofs of this result use the well-known but nontrivial fact that class groups of rings of integers are torsion. In this short note, we provide an independent proof using elementary techniques.


title: Transferable and extensible machine learning derived atomic charges for modeling hybrid nanoporous materials

abstract: Nanoporous materials have attracted significant interest as an emerging platform for adsorption-related applications. The high-throughput computational screening became a standard technique to access the performance of thousands of candidates, but its accuracy is highly dependent on a partial charge assignment method. In this study, we propose a machine learning model that can reconcile the benefits of two main approaches-the high accuracy of density-derived electrostatic and chemical charge (DDEC) method and the scalability of charge equilibration (Qeq) method. The mean absolute deviation of predicted partial charges from the original DDEC counterparts archive an excellent level of 0.01e. The model, initially designed for metal-organic frameworks, is also capable of assigning charges to another class of nanoporous materials, covalent organic frameworks, with acceptable accuracy. Adsorption properties of carbon dioxide, calculated by means of machine learning derived charges, are consistent with the reference data obtained with DDEC charges.


title: Parallax: Visualizing and Understanding the Semantics of Embedding Spaces via Algebraic Formulae

abstract: Embeddings are a fundamental component of many modern machine learning and natural language processing models. Understanding them and visualizing them is essential for gathering insights about the information they capture and the behavior of the models. State of the art in analyzing embeddings consists in projecting them in two-dimensional planes without any interpretable semantics associated to the axes of the projection, which makes detailed analyses and comparison among multiple sets of embeddings challenging. In this work, we propose to use explicit axes defined as algebraic formulae over embeddings to project them into a lower dimensional, but semantically meaningful subspace, as a simple yet effective analysis and visualization methodology. This methodology assigns an interpretable semantics to the measures of variability and the axes of visualizations, allowing for both comparisons among different sets of embeddings and fine-grained inspection of the embedding spaces. We demonstrate the power of the proposed methodology through a series of case studies that make use of visualizations constructed around the underlying methodology and through a user study. The results show how the methodology is effective at providing more profound insights than classical projection methods and how it is widely applicable to many other use cases.


title: Using local plasticity rules to train recurrent neural networks

abstract: To learn useful dynamics on long time scales, neurons must use plasticity rules that account for long-term, circuit-wide effects of synaptic changes. In other words, neural circuits must solve a credit assignment problem to appropriately assign responsibility for global network behavior to individual circuit components. Furthermore, biological constraints demand that plasticity rules are spatially and temporally local; that is, synaptic changes can depend only on variables accessible to the pre- and postsynaptic neurons. While artificial intelligence offers a computational solution for credit assignment, namely backpropagation through time (BPTT), this solution is wildly biologically implausible. It requires both nonlocal computations and unlimited memory capacity, as any synaptic change is a complicated function of the entire history of network activity. Similar nonlocality issues plague other approaches such as FORCE (Sussillo et al. 2009). Overall, we are still missing a model for learning in recurrent circuits that both works computationally and uses only local updates. Leveraging recent advances in machine learning on approximating gradients for BPTT, we derive biologically plausible plasticity rules that enable recurrent networks to accurately learn long-term dependencies in sequential data. The solution takes the form of neurons with segregated voltage compartments, with several synaptic sub-populations that have different functional properties. The network operates in distinct phases during which each synaptic sub-population is updated by its own local plasticity rule. Our results provide new insights into the potential roles of segregated dendritic compartments, branch-specific inhibition, and global circuit phases in learning.


title: Differential Privacy Has Disparate Impact on Model Accuracy

abstract: Differential privacy (DP) is a popular mechanism for training machine learning models with bounded leakage about the presence of specific points in the training data. The cost of differential privacy is a reduction in the model's accuracy. We demonstrate that in the neural networks trained using differentially private stochastic gradient descent (DP-SGD), this cost is not borne equally: accuracy of DP models drops much more for the underrepresented classes and subgroups. For example, a gender classification model trained using DP-SGD exhibits much lower accuracy for black faces than for white faces. Critically, this gap is bigger in the DP model than in the non-DP model, i.e., if the original model is unfair, the unfairness becomes worse once DP is applied. We demonstrate this effect for a variety of tasks and models, including sentiment analysis of text and image classification. We then explain why DP training mechanisms such as gradient clipping and noise addition have disproportionate effect on the underrepresented and more complex subgroups, resulting in a disparate reduction of model accuracy.


title: Conservative numerical schemes with optimal dispersive wave relations -- Part I. Derivations and analyses

abstract: An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived, by approximating the Hamiltonian formulation, based on the Poisson brackets and the vorticity-divergence variables, of the inviscid shallow water flows. The conservation of the energy and/or enstrophy stems from skew-symmetry of the Poisson brackets, which is retained in the discrete approximations. These schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme.


title: Competitive Gradient Descent

abstract: We introduce a new algorithm for the numerical computation of Nash equilibria of competitive two-player games. Our method is a natural generalization of gradient descent to the two-player setting where the update is given by the Nash equilibrium of a regularized bilinear local approximation of the underlying game. It avoids oscillatory and divergent behaviors seen in alternating gradient descent. Using numerical experiments and rigorous analysis, we provide a detailed comparison to methods based on \emph{optimism} and \emph{consensus} and show that our method avoids making any unnecessary changes to the gradient dynamics while achieving exponential (local) convergence for (locally) convex-concave zero sum games. Convergence and stability properties of our method are robust to strong interactions between the players, without adapting the stepsize, which is not the case with previous methods. In our numerical experiments on non-convex-concave problems, existing methods are prone to divergence and instability due to their sensitivity to interactions among the players, whereas we never observe divergence of our algorithm. The ability to choose larger stepsizes furthermore allows our algorithm to achieve faster convergence, as measured by the number of model evaluations.


title: Heuristics in Multi-Winner Approval Voting

abstract: In many real world situations, collective decisions are made using voting. Moreover, scenarios such as committee or board elections require voting rules that return multiple winners. In multi-winner approval voting (AV), an agent may vote for as many candidates as they wish. Winners are chosen by tallying up the votes and choosing the top-$k$ candidates receiving the most votes. An agent may manipulate the vote to achieve a better outcome by voting in a way that does not reflect their true preferences. In complex and uncertain situations, agents may use heuristics to strategize, instead of incurring the additional effort required to compute the manipulation which most favors them. In this paper, we examine voting behavior in multi-winner approval voting scenarios with complete information. We show that people generally manipulate their vote to obtain a better outcome, but often do not identify the optimal manipulation. Instead, voters tend to prioritize the candidates with the highest utilities. Using simulations, we demonstrate the effectiveness of these heuristics in situations where agents only have access to partial information.


title: Certifiably Robust Interpretation in Deep Learning

abstract: Deep learning interpretation is essential to explain the reasoning behind model predictions. Understanding the robustness of interpretation methods is important especially in sensitive domains such as medical applications since interpretation results are often used in downstream tasks. Although gradient-based saliency maps are popular methods for deep learning interpretation, recent works show that they can be vulnerable to adversarial attacks. In this paper, we address this problem and provide a certifiable defense method for deep learning interpretation. We show that a sparsified version of the popular SmoothGrad method, which computes the average saliency maps over random perturbations of the input, is certifiably robust against adversarial perturbations. We obtain this result by extending recent bounds for certifiably robust smooth classifiers to the interpretation setting. Experiments on ImageNet samples validate our theory.


title: EM Converges for a Mixture of Many Linear Regressions

abstract: We study the convergence of the Expectation-Maximization (EM) algorithm for mixtures of linear regressions with an arbitrary number $k$ of components. We show that as long as signal-to-noise ratio (SNR) is $\tilde{\Omega}(k)$, well-initialized EM converges to the true regression parameters. Previous results for $k \geq 3$ have only established local convergence for the noiseless setting, i.e., where SNR is infinitely large. Our results enlarge the scope to the environment with noises, and notably, we establish a statistical error rate that is independent of the norm (or pairwise distance) of the regression parameters. In particular, our results imply exact recovery as $\sigma \rightarrow 0$, in contrast to most previous local convergence results for EM, where the statistical error scaled with the norm of parameters. Standard moment-method approaches may be applied to guarantee we are in the region where our local convergence guarantees apply.


title: SpArSe: Sparse Architecture Search for CNNs on Resource-Constrained Microcontrollers

abstract: The vast majority of processors in the world are actually microcontroller units (MCUs), which find widespread use performing simple control tasks in applications ranging from automobiles to medical devices and office equipment. The Internet of Things (IoT) promises to inject machine learning into many of these every-day objects via tiny, cheap MCUs. However, these resource-impoverished hardware platforms severely limit the complexity of machine learning models that can be deployed. For example, although convolutional neural networks (CNNs) achieve state-of-the-art results on many visual recognition tasks, CNN inference on MCUs is challenging due to severe finite memory limitations. To circumvent the memory challenge associated with CNNs, various alternatives have been proposed that do fit within the memory budget of an MCU, albeit at the cost of prediction accuracy. This paper challenges the idea that CNNs are not suitable for deployment on MCUs. We demonstrate that it is possible to automatically design CNNs which generalize well, while also being small enough to fit onto memory-limited MCUs. Our Sparse Architecture Search method combines neural architecture search with pruning in a single, unified approach, which learns superior models on four popular IoT datasets. The CNNs we find are more accurate and up to $4.35\times$ smaller than previous approaches, while meeting the strict MCU working memory constraint.


title: Fundamental differences in the radio properties of red and blue quasars: evolution strongly favoured over orientation

abstract: A minority of the optically selected quasar population are red at optical wavelengths due to the presence of dust along the line-of-sight. A key focus of many red quasar studies is to understand their relationship with the overall quasar population: are they blue quasars observed at a (slight) inclination angle or do they represent a transitional phase in the evolution of quasars? Identifying fundamental differences between red and blue quasars is key to discriminate between these two paradigms. To robustly explore this, we have uniformly selected quasars from the Sloan Digital Sky Survey with mid-infrared counterparts, carefully controlling for luminosity and redshift effects. We take a novel approach to distinguish between colour-selected quasars in the redshift range of $0.2 < z < 2.4$ by constructing redshift-sensitive $g^* - i^*$ colour cuts. From cross-matching this sample to the Faint Images of the Radio Sky at Twenty-centimeters (FIRST) survey, we have found a factor $\approx$ 3 larger fraction of radio-detected red quasars with respect to that of blue quasars. Through a visual inspection of the FIRST images and an assessment of the radio luminosities (rest-frame ${L_{\rm 1.4~GHz}}$ and ${L_{\rm 1.4~GHz}}/{L_{\rm 6\mu m}}$) we find that the radio-detection excess for red quasars is primarily due to compact and radio-faint systems (around the radio quiet--radio loud threshold). We show that our results rule out orientation as the origin for the differences between red and blue quasars and argue that they provide broad agreement with an evolutionary model.


title: What They Don't Want: An Analysis of Brexit's First Round of Indicative Votes

abstract: Since the result of the 2016 referendum, Brexit has been an unpredictable democratic adventure, the finale of which remains unclear. This year, in the final days of March, parliamentarians seized control of the order paper from the Government and held their own indicative votes in an attempt to break the deadlock. In this paper we analyse the results of this unusual cardinal ballot. We express the various motions in terms of `how much Brexit' they deliver, and employ Monte Carlo in an attempt to determine this from the data. We find solutions which reproduce our intuitive understanding of the debate. Finally, we construct hypothetical ordinal ballots for the various Brexit scenarios, using three different processes. The results suggest that the Government would be more successful taking a softer position, and we quantify this. Additionally, there is some discussion of how tactical voting might manifest itself in the event of such an exercise.


title: Inducing Uniform Asymptotic Stability in Non-Autonomous Accelerated Optimization Dynamics via Hybrid Regularization

abstract: There have been many recent efforts to study accelerated optimization algorithms from the perspective of dynamical systems. In this paper, we focus on the robustness properties of the time-varying continuous-time version of these dynamics. These properties are critical for the implementation of accelerated algorithms in feedback-based control and optimization architectures. We show that a family of dynamics related to the continuous-time limit of Nesterov's accelerated gradient method can be rendered unstable under arbitrarily small bounded disturbances. Indeed, while solutions of these dynamics may converge to the set of optimizers, in general, this set may not be uniformly asymptotically stable. To induce uniformity, and robustness as a byproduct, we propose a framework where we regularize the dynamics by using resetting mechanisms that are modeled by well-posed hybrid dynamical systems. For these hybrid dynamics, we establish uniform asymptotic stability and robustness properties, as well as convergence rates that are similar to those of the non-hybrid dynamics. We finish by characterizing a family of discretization mechanisms that retain the main stability and robustness properties of the hybrid algorithms.


title: Analyzing and Supporting Adaptation of Online Code Examples

abstract: Developers often resort to online Q&A forums such as Stack Overflow (SO) for filling their programming needs. Although code examples on those forums are good starting points, they are often incomplete and inadequate for developers' local program contexts; adaptation of those examples is necessary to integrate them to production code. As a consequence, the process of adapting online code examples is done over and over again, by multiple developers independently. Our work extensively studies these adaptations and variations, serving as the basis for a tool that helps integrate these online code examples in a target context in an interactive manner. We perform a large-scale empirical study about the nature and extent of adaptations and variations of SO snippets. We construct a comprehensive dataset linking SO posts to GitHub counterparts based on clone detection, time stamp analysis, and explicit URL references. We then qualitatively inspect 400 SO examples and their GitHub counterparts and develop a taxonomy of 24 adaptation types. Using this taxonomy, we build an automated adaptation analysis technique on top of GumTree to classify the entire dataset into these types. We build a Chrome extension called ExampleStack that automatically lifts an adaptation-aware template from each SO example and its GitHub counterparts to identify hot spots where most changes happen. A user study with sixteen programmers shows that seeing the commonalities and variations in similar GitHub counterparts increases their confidence about the given SO example, and helps them grasp a more comprehensive view about how to reuse the example differently and avoid common pitfalls.


title: Multiple Boris integrators for particle-in-cell simulation

abstract: We construct Boris-type schemes for integrating the motion of charged particles in particle-in-cell (PIC) simulation. The new solvers virtually combine the 2-step Boris procedure arbitrary n times in the Lorentz-force part, and therefore we call them the multiple Boris solvers. Using Chebyshev polynomials, a one-step form of the new solvers is provided. The new solvers give n^2 times smaller errors, allow larger timesteps, and have a long-term stability. We present numerical tests of the new solvers, in comparison with other particle integrators.

categories: physics.comp-ph astro-ph.HE physics.plasm-ph physics.space-ph


title: Hyperelliptic limits of quadrics through canonical curves and ribbons

abstract: We describe explicitly all hyperelliptic limits of quadrics through smooth canonical curves of genus $g$ in ${\mathbb P}^{g-1}$. Also, we construct an open embedding of the blow up of a ${\rm PGL}_g$-bundle over the moduli space of curves of genus $g$ along the hyperelliptic locus into the blow up of the canonical Hilbert scheme of ${\mathbb P}^{g-1}$ along the closure of the locus of canonical ribbons, which are certain double thickenings of rational normal curves introduced and studied by Bayer and Eisenbud.


title: Theory of coherent active convolved illumination for superresolution enhancement

abstract: Recently an optical amplification process called the plasmon injection scheme was introduced as an effective solution to overcoming losses in metamaterials. Implementations with near-field imaging applications have indicated substantial performance enhancements even in the presence of noise. This powerful and versatile compensation technique, which has since been renamed to a more generalized active convolved illumination, offers new possibilities of improving the performance of many previously conceived metamaterial-based devices and conventional imaging systems. In this work, we present the first comprehensive mathematical breakdown of active convolved illumination for coherent imaging. Our analysis highlights the distinctive features of active convolved illumination, such as selective spectral amplification and correlations, and provides a rigorous understanding of the loss compensation process. These features are achieved by an auxiliary source coherently superimposed with the object field. The auxiliary source is designed to have three important properties. First, it is correlated with the object field. Second, it is defined over a finite spectral bandwidth. Third, it is amplified over that selected bandwidth. We derive the variance for the image spectrum and show that utilizing the auxiliary source with the above properties can significantly improve the spectral signal-to-noise ratio and resolution limit. Besides enhanced superresolution imaging, the theory can be potentially generalized to the compensation of information or photon loss in a wide variety of coherent and incoherent linear systems including those, for example, in atmospheric imaging, time-domain spectroscopy, ${\cal PT}$ symmetric non-Hermitian photonics, and even quantum computing.


title: AdaOja: Adaptive Learning Rates for Streaming PCA

abstract: Oja's algorithm has been the cornerstone of streaming methods in Principal Component Analysis (PCA) since it was first proposed in 1982. However, Oja's algorithm does not have a standardized choice of learning rate (step size) that both performs well in practice and truly conforms to the online streaming setting. In this paper, we propose a new learning rate scheme for Oja's method called AdaOja. This new algorithm requires only a single pass over the data and does not depend on knowing properties of the data set a priori. AdaOja is a novel variation of the Adagrad algorithm to Oja's algorithm in the single eigenvector case and extended to the multiple eigenvector case. We demonstrate for dense synthetic data, sparse real-world data and dense real-world data that AdaOja outperforms common learning rate choices for Oja's method. We also show that AdaOja performs comparably to state-of-the-art algorithms (History PCA and Streaming Power Method) in the same streaming PCA setting.


title: Harnessing Slow Dynamics in Neuromorphic Computation

abstract: Neuromorphic Computing is a nascent research field in which models and devices are designed to process information by emulating biological neural systems. Thanks to their superior energy efficiency, analog neuromorphic systems are highly promising for embedded, wearable, and implantable systems. However, optimizing neural networks deployed on these systems is challenging. One main challenge is the so-called timescale mismatch: Dynamics of analog circuits tend to be too fast to process real-time sensory inputs. In this thesis, we propose a few working solutions to slow down dynamics of on-chip spiking neural networks. We empirically show that, by harnessing slow dynamics, spiking neural networks on analog neuromorphic systems can gain non-trivial performance boosts on a battery of real-time signal processing tasks.


title: The Evolution of Brightest Cluster Galaxies in the Nearby Universe I: Colours and Stellar Masses from the Sloan Digital Sky Survey and Wide Infrared Survey Explorer

abstract: We present a study of the evolution of brightest cluster galaxies (BCGs) in a sample of clusters at $0.05 \leq z<0.35$ from the SDSS and WISE with halo masses in the range $6 \times 10^{13}M_\odot$ (massive groups) - $10^{15.5}M_\odot$ (Coma-like clusters). We analyse optical and infrared colours and stellar masses of BCGs as a function of the mass of their host haloes. We find that BCGs are mostly red and quiescent galaxies and that a minority ($\sim 9$\%) of them are star-forming. We find that the optical $g-r$ colours are consistent with those of red sequence galaxies at the same redshifts; however, we detect the presence of a tail of blue and mostly star-forming BCGs preferentially located in low-mass clusters and groups. Although the blue tail is dominated by star-forming galaxies, we find that star-forming BCGs may also have red $g-r$ colours, indicating dust-enshrouded star formation. The fraction of star-forming BCGs increases with redshift and decreases with cluster mass and BCG stellar mass. We find that cool-core clusters host both star-forming and quiescent BCGs; however, non cool-core clusters are dominated by quiescent BCGs. Star formation appears thus as the result of processes that depend on stellar mass, cluster mass and cooling state of the intra-cluster medium. Our results suggest no significant stellar mass growth at $z<0.35$, supporting the notion that BCGs had accreted most of their mass by $z = 0.35$. Overall we find a low (1\%) AGN fraction detected at IR wavelengths.


title: Cyclicality, Periodicity and the Topology of Time Series

abstract: Periodic and semi periodic patterns are very common in nature. In this paper we introduce a topological toolbox aiming in detecting and quantifying periodicity. The presented technique is of a general nature and may be employed wherever there is suspected cyclic behaviour in a time series with no trend. The approach is tested on a number of real-world examples enabling us to consistently demonstrate an ability to recognise periodic behaviour where conventional techniques fail to do so. Quicker to react to changes in time series behaviour, and with a high robustness to noise, the toolbox offers a powerful way to deeper understanding of time series dynamics.


title: Order reduction methods for solving large-scale differential matrix Riccati equations

abstract: We consider the numerical solution of large-scale symmetric differential matrix Riccati equations. Under certain hypotheses on the data, reduced order methods have recently arisen as a promising class of solution strategies, by forming low-rank approximations to the sought after solution at selected timesteps. We show that great computational and memory savings are obtained by a reduction process onto rational Krylov subspaces, as opposed to current approaches. By specifically addressing the solution of the reduced differential equation and reliable stopping criteria, we are able to obtain accurate final approximations at low memory and computational requirements. This is obtained by employing a two-phase strategy that separately enhances the accuracy of the algebraic approximation and the time integration. The new method allows us to numerically solve much larger problems than in the current literature. Numerical experiments on benchmark problems illustrate the effectiveness of the procedure with respect to existing solvers.


title: Deep Dilated Convolutional Nets for the Automatic Segmentation of Retinal Vessels

abstract: The reliable segmentation of retinal vasculature can provide the means to diagnose and monitor the progression of a variety of diseases affecting the blood vessel network, including diabetes and hypertension. We leverage the power of convolutional neural networks to devise a reliable and fully automated method that can accurately detect, segment, and analyze retinal vessels. In particular, we propose a novel, fully convolutional deep neural network with an encoder-decoder architecture that employs dilated spatial pyramid pooling with multiple dilation rates to recover the lost content in the encoder and add multiscale contextual information to the decoder. We also propose a simple yet effective way of quantifying and tracking the widths of retinal vessels through direct use of the segmentation predictions. Unlike previous deep-learning-based approaches to retinal vessel segmentation that mainly rely on patch-wise analysis, our proposed method leverages a whole-image approach during training and inference, resulting in more efficient training and faster inference through the access of global content in the image. We have tested our method on two publicly available datasets, and our state-of-the-art results on both the DRIVE and CHASE-DB1 datasets attest to the effectiveness of our approach.


title: An Investigation of Data Poisoning Defenses for Online Learning

abstract: Data poisoning attacks -- where an adversary can modify a small fraction of training data, with the goal of forcing the trained classifier to high loss -- are an important threat for machine learning in many applications. While a body of prior work has developed attacks and defenses, there is not much general understanding on when various attacks and defenses are effective. In this work, we undertake a rigorous study of defenses against data poisoning for online learning. First, we study four standard defenses in a powerful threat model, and provide conditions under which they can allow or resist rapid poisoning. We then consider a weaker and more realistic threat model, and show that the success of the adversary in the presence of data poisoning defenses there depends on the "ease" of the learning problem.


title: Deep Learning Moment Closure Approximations using Dynamic Boltzmann Distributions

abstract: The moments of spatial probabilistic systems are often given by an infinite hierarchy of coupled differential equations. Moment closure methods are used to approximate a subset of low order moments by terminating the hierarchy at some order and replacing higher order terms with functions of lower order ones. For a given system, it is not known beforehand which closure approximation is optimal, i.e. which higher order terms are relevant in the current regime. Further, the generalization of such approximations is typically poor, as higher order corrections may become relevant over long timescales. We have developed a method to learn moment closure approximations directly from data using dynamic Boltzmann distributions (DBDs). The dynamics of the distribution are parameterized using basis functions from finite element methods, such that the approach can be applied without knowing the true dynamics of the system under consideration. We use the hierarchical architecture of deep Boltzmann machines (DBMs) with multinomial latent variables to learn closure approximations for progressively higher order spatial correlations. The learning algorithm uses a centering transformation, allowing the dynamic DBM to be trained without the need for pre-training. We demonstrate the method for a Lotka-Volterra system on a lattice, a typical example in spatial chemical reaction networks. The approach can be applied broadly to learn deep generative models in applications where infinite systems of differential equations arise.


title: Higher Correlations and the Alternative Hypothesis

abstract: The Alternative Hypothesis concerns a hypothetical and unlikely picture of how zeros of the Riemann zeta function are spaced which one would like to rule out. In the Alternative Hypothesis, the renormalized distance between nontrivial zeros is supposed to always lie at a half integer. It is known that the Alternative Hypothesis is compatible with what is known about the pair correlation function of zeta zeros. We ask whether what is currently known about higher correlation functions of the zeros is sufficient to rule out the Alternative Hypothesis and show by construction of an explicit counterexample point process that it is not. A similar result was recently independently obtained by T. Tao, using slightly different methods. We also apply the ergodic theorem to this point process to show there exists a deterministic collection of points lying in $\tfrac{1}{2}\mathbb{Z}$ which satisfy the Alternative Hypothesis spacing but mimic all statistics which are currently known about zeros of the zeta function.


title: Moduli of flat connections on smooth varieties

abstract: We study the moduli functor of flat bundles on smooth, possibly non-proper, algebraic variety $X$ (over a field of characteristic zero). For this we introduce the notion of \emph{formal boundary} of $X$, denoted by $\partial X$, which is a formal analogue of the boundary at infinity of the Betti topological space associated to $X$. We explain how to construct two derived moduli functors $Vect^{\nabla}(X)$ and $Vect^{\nabla}(\partial X)$, of flat bundles on $X$ and on $\partial X$, as well as a restriction map $R : Vect^{\nabla}(X) \rightarrow Vect^\nabla(\partial X)$ from the former to the later. This work contains two main results. First we prove that the morphism R comes equipped with a canonical shifted Lagrangian structure in the sense of [PTVV]. This first result can be understood as the de Rham analogue of the existence of Poisson structures on moduli of local systems previously studied by the authors. As a second statement, we prove that the geometric fibers of $R$ are representable by "quasi-algebraic spaces", a slight weakening of the notion of algebraic spaces.