title: Functional SAC model: With application to spatial econometrics

abstract: Spatial autoregressive combined (SAC) model has been widely studied in the literature for the analysis of spatial data in various areas such as geography, economics, demography, regional sciences. This is a linear model with scalar response, scalar explanatory variables and which allows for spatial interactions in the dependent variable and the disturbances. In this work we extend this modeling approach from scalar to functional covariate. The parameters of the model are estimated via the maximum likelihood estimation method. A simulation study is conducted to evaluate the performance of the proposed methodology. As an illustration, the model is used to establish the relationship between unemployment and illiteracy in Senegal.


title: Magnetic field generation from composition gradients in inertial confinement fusion fuel

abstract: Experimental asymmetries in fusion implosions can lead to magnetic field generation in the hot plasma core. For typical parameters, previous studies found that the magnetisation Hall parameter, given by the product of the electron gyro-frequency and Coulomb collision time, can exceed one. This will affect the hydrodynamics through inhibition and deflection of the electron heat flux. The magnetic field source is the collisionless Biermann term, which arises from the Debye shielding potential in electron pressure gradients. We show that there is an additional source term due to the Z dependence of the Coulomb collision operator. If there are ion composition gradients, such as jets of carbon ablator mix entering the hot-spot, this source term can rapidly exceed the Biermann fields. In addition, the Biermann fields are enhanced due to the increased temperature gradients from carbon radiative cooling. With even stronger self-generated fields, heat loss to the carbon regions will be reduced, potentially reducing the negative effect of carbon mix.


title: Accretion Disc Winds in Tidal Disruption Events: Ultraviolet Spectral Lines as Orientation Indicators

abstract: Some tidal disruption events (TDEs) exhibit blueshifted broad absorption lines (BALs) in their rest-frame ultraviolet (UV) spectra, while others display broad emission lines (BELs). Similar phenomenology is observed in quasars and accreting white dwarfs, where it can be interpreted as an orientation effect associated with line formation in an accretion disc wind.We propose and explore a similar unification scheme for TDEs. We present synthetic UV spectra for disc and wind-hosting TDEs, produced by a state-of-the-art Monte Carlo ionization and radiative transfer code. Our models cover a wide range of disc wind geometries and kinematics. Such winds naturally reproduce both BALs and BELs. In general, sight lines looking into the wind cone preferentially produce BALs, while other orientations preferentially produce BELs. We also study the effect of wind clumping and CNO-processed abundances on the observed spectra. Clumpy winds tend to produce stronger UV emission and absorption lines, because clumping increases both the emission measure and the abundances of the relevant ionic species, the latter by reducing the ionization state of the outflow. The main effect of adopting CNO-processed abundances is a weakening of C~{\sc iv}~1550~\AA~ and an enhancement of N \textsc{v}~1240~\AA~ in the spectra. We conclude that line formation in an accretion disc wind is a promising mechanism for explaining the diverse UV spectra of TDEs. If this is correct, the relative number of BAL and BEL TDEs can be used to estimate the covering factor of the outflow. The models in this work are publicly available online and upon request.


title: Image Quality Assessment: Unifying Structure and Texture Similarity

abstract: Objective measures of image quality generally operate by comparing pixels of a "degraded" image to those of the original. Relative to human observers, these measures are overly sensitive to resampling of texture regions (e.g., replacing one patch of grass with another). Here, we develop the first full-reference image quality model with explicit tolerance to texture resampling. Using a convolutional neural network, we construct an injective and differentiable function that transforms images to multi-scale overcomplete representations. We demonstrate empirically that the spatial averages of the feature maps in this representation capture texture appearance, in that they provide a set of sufficient statistical constraints to synthesize a wide variety of texture patterns. We then describe an image quality method that combines correlations of these spatial averages ("texture similarity") with correlations of the feature maps ("structure similarity"). The parameters of the proposed measure are jointly optimized to match human ratings of image quality, while minimizing the reported distances between subimages cropped from the same texture images. Experiments show that the optimized method explains human perceptual scores, both on conventional image quality databases, as well as on texture databases. The measure also offers competitive performance on related tasks such as texture classification and retrieval. Finally, we show that our method is relatively insensitive to geometric transformations (e.g., translation and dilation), without use of any specialized training or data augmentation. Code is available at https://github.com/dingkeyan93/DISTS.


title: Estimating quantum coherence by noncommutativity of any observable and its incoherent part

abstract: We establish a lower bound on the quantum coherence of an arbitrary quantum state in arbitrary dimension, using a noncommutativity estimator of an arbitrary observable of sub-unit norm, where the estimator is the commutator of the observable and its incoherent or classical part. The relation provides a direct method of obtaining an estimate of the quantum coherence of an arbitrary quantum state, without resorting to quantum state tomography or the existing witness operators.


title: Typical knots: size, link component count, and writhe

abstract: We model the typical behavior of knots and links using grid diagrams. Links are ubiquitous in the sciences, and their "normal" or "typical" behavior is of significant importance in understanding situations such as the topological state of DNA or the statistical mechanics of ring polymers. We examine three invariants: the expected size of a random knot; the expected number of components of a random link; and the expected writhe of a random knot. We investigate the first two numerically and produce generating functions which codify the observed patterns: knot size is uniformly distributed and linearly dependent upon grid size, and the number of components follows a distribution whose mean and variance grow with log_2 of grid size; in particular, for any fixed k, the k-component links grow vanishingly rare as grid size increases. Finally, we observe that the odd moments of writhe vanish, and we perform an exploratory data analysis to discover that variance grows with the square of grid size and kurtosis is constant at approximately 3.5. We continue this project in a future work, where we investigate genus and the effects of crossing change on it.

categories: math.GT cond-mat.stat-mech math.CO q-bio.QM


title: Experimental Data from a Quantum Computer Verifies the Generalized Pauli Exclusion Principle

abstract: "What are the consequences ... that Fermi particles cannot get into the same state ... " R. P. Feynman wrote of the Pauli exclusion principle, "In fact, almost all the peculiarities of the material world hinge on this wonderful fact." In 1972 Borland and Dennis showed that there exist powerful constraints beyond the Pauli exclusion principle on the orbital occupations of Fermi particles, providing important restrictions on quantum correlation and entanglement. Here we use computations on quantum computers to experimentally verify the existence of these additional constraints. Quantum many-fermion states are randomly prepared on the quantum computer and tested for constraint violations. Measurements show no violation and confirm the generalized Pauli exclusion principle with an error of one part in one quintillion.


title: Multigrid for Chiral Lattice Fermions: Domain Wall

abstract: Critical slowing down for the Krylov Dirac solver presents a major obstacle to further advances in lattice field theory as it approaches the continuum solution. We propose a new multi-grid approach for chiral fermions, applicable to both the 5-d domain wall or 4-d Overlap operator. The central idea is to directly coarsen the 4-d Wilson kernel, giving an effective domain wall or overlap operator on each level. We provide here an explicit construction for the Shamir domain wall formulation with numerical tests for the 2-d Schwinger prototype, demonstrating near ideal multi-grid scaling. The framework is designed for a natural extension to 4-d lattice QCD chiral fermions, such as the M\"obius, Zolotarev or Borici domain wall discretizations or directly to a rational expansion of the 4-d Overlap operator. For the Shamir operator, the effective overlap operator is isolated by the use of a Pauli-Villars preconditioner in the spirit of the K\"ahler-Dirac spectral map used in a recent staggered MG algorithm [1].


title: Bridging the Gap: FPGAs as Programmable Switches

abstract: The emergence of P4, a domain specific language, coupled to PISA, a domain specific architecture, is revolutionizing the networking field. P4 allows to describe how packets are processed by a programmable data plane, spanning ASICs and CPUs, implementing PISA. Because the processing flexibility can be limited on ASICs, while the CPUs performance for networking tasks lag behind, recent works have proposed to implement PISA on FPGAs. However, little effort has been dedicated to analyze whether FPGAs are good candidates to implement PISA. In this work, we take a step back and evaluate the micro-architecture efficiency of various PISA blocks. We demonstrate, supported by a theoretical and experimental analysis, that the performance of a few PISA blocks is severely limited by the current FPGA architectures. Specifically, we show that match tables and programmable packet schedulers represent the main performance bottlenecks for FPGA-based programmable switches. Thus, we explore two avenues to alleviate these shortcomings. First, we identify network applications well tailored to current FPGAs. Second, to support a wider range of networking applications, we propose modifications to the FPGA architectures which can also be of interest out of the networking field.


title: Stationary sets and asymptotic behavior of the mean curvature flow with forcing in the plane

abstract: We consider the flat flow solutions of the mean curvature equation with a forcing term in the plane. We prove that for every constant forcing term the stationary sets are given by a finite union of disks with equal radii and disjoint closures. On the other hand for every bounded forcing term tangent disks are never stationary. Finally in the case of an asymptotically constant forcing term we show that the only possible long time limit sets are given by disjoint unions of disks with equal radii and possibly tangent.


title: Maximum likelihood degree, complete quadrics and ${\mathbb C}^*$-action

abstract: We study the maximum likelihood (ML) degree of linear concentration models in algebraic statistics. We relate it to an intersection problem on the variety of complete quadrics. This allows us to provide an explicit, basic, albeit of high computational complexity, formula for the ML-degree. The variety of complete quadrics is an exact analog for symmetric matrices of the permutohedron variety for the diagonal matrices.


title: Machine learning for multiple yield curve markets: fast calibration in the Gaussian affine framework

abstract: Calibration is a highly challenging task, in particular in multiple yield curve markets. This paper is a first attempt to study the chances and challenges of the application of machine learning techniques for this. We employ Gaussian process regression, a machine learning methodology having many similarities with extended Kalman filtering - a technique which has been applied many times to interest rate markets and term structure models. We find very good results for the single curve markets and many challenges for the multi curve markets in a Vasicek framework. The Gaussian process regression is implemented with the Adam optimizer and the non-linear conjugate gradient method, where the latter performs best. We also point towards future research.


title: Cross-lingual Contextualized Topic Models with Zero-shot Learning

abstract: Many data sets (e.g., reviews, forums, news, etc.) exist parallelly in multiple languages. They all cover the same content, but the linguistic differences make it impossible to use traditional, bag-of-word-based topic models. Models have to be either single-language or suffer from a huge, but extremely sparse vocabulary. Both issues can be addressed by transfer learning. In this paper, we introduce a zero-shot cross-lingual topic model. Our model learns topics on one language (here, English), and predicts them for unseen documents in different languages (here, Italian, French, German, and Portuguese). We evaluate the quality of the topic predictions for the same document in different languages. Our results show that the transferred topics are coherent and stable across languages, which suggests exciting future research directions.


title: Inversion of a SIR-based model: a critical analysis about the application to COVID-19 epidemic

abstract: Calibration of a SIR (Susceptibles-Infected-Recovered) model with official international data for the COVID-19 pandemics provides a good example of the difficulties inherent the solution of inverse problems. Inverse modeling is set up in a framework of discrete inverse problems, which explicitly considers the role and the relevance of data. Together with a physical vision of the model, the present work addresses numerically the issue of parameters calibration in SIR models, it discusses the uncertainties in the data provided by international authorities, how they influence the reliability of calibrated model parameters and, ultimately, of model predictions.


title: Quantum-classical hybrid algorithm using an error-mitigating $N$-representability condition to compute the Mott metal-insulator transition

abstract: Quantum algorithms for molecular electronic structure have been developed with lower computational scaling than their classical counterparts, but emerging quantum hardware is far from being capable of the coherence,connectivity and gate errors required for their experimental realization. Here we propose a class of quantum-classical hybrid algorithms that compute the energy from a two-electron reduced density matrix (2-RDM). The 2-RDM is constrained by $N$-representability conditions, conditions for representing an $N$-electron wavefunction, that mitigates noise from the quantum circuit. We compute the strongly correlated dissociation of doublet H$_{3}$ into three hydrogen atoms. The hybrid quantum-classical computer matches the energies from full configuration interaction to 0.1 kcal/mol, one-tenth of "chemical accuracy," even in the strongly correlated limit of dissociation. Furthermore, the spatial locality of the computed one-electron RDM reveals that the quantum computer accurately predicts the Mott metal-insulator transition.


title: Really Useful Synthetic Data -- A Framework to Evaluate the Quality of Differentially Private Synthetic Data

abstract: Recent advances in generating synthetic data that allow to add principled ways of protecting privacy -- such as Differential Privacy -- are a crucial step in sharing statistical information in a privacy preserving way. But while the focus has been on privacy guarantees, the resulting private synthetic data is only useful if it still carries statistical information from the original data. To further optimise the inherent trade-off between data privacy and data quality, it is necessary to think closely about the latter. What is it that data analysts want? Acknowledging that data quality is a subjective concept, we develop a framework to evaluate the quality of differentially private synthetic data from an applied researcher's perspective. Data quality can be measured along two dimensions. First, quality of synthetic data can be evaluated against training data or against an underlying population. Second, the quality of synthetic data depends on general similarity of distributions or specific tasks such as inference or prediction. It is clear that accommodating all goals at once is a formidable challenge. We invite the academic community to jointly advance the privacy-quality frontier.


title: On the bounds of coefficients of Daubechie's orthonormal wavelets

abstract: In the present work we provide the bounds for Daubechies orthonormal wavelet coefficients for function spaces $\mathcal{A}_k^p:=\{f: \|(i \omega)^k\hat{f}(\omega)\|_p< \infty\}$, $k\in\mathbf{N}\cup\{0\}$, $p\in(1,\infty)$.


title: CO.ME.T.A. -- covid-19 media textual analysis. A dashboard for media monitoring

abstract: The focus of this paper is to trace how mass media, particularly newspapers, have addressed the issues about the containment of contagion or the explanation of epidemiological evolution. We propose an interactive dashboard: CO.ME.T.A.. During crises it is important to shape the best communication strategies in order to respond to critical situations. In this regard, it is important to monitor the information that mass media and social platforms convey. The dashboard allows to explore the mining of contents extracted and study the lexical structure that links the main discussion topics. The dashboard merges together four methods: text mining, sentiment analysis, textual network analysis and latent topic models. Results obtained on a subset of documents show not only a health-related semantic dimension, but it also extends to social-economic dimensions.


title: BETS: The dangers of selection bias in early analyses of the coronavirus disease (COVID-19) pandemic

abstract: The coronavirus disease 2019 (COVID-19) has quickly grown from a regional outbreak in Wuhan, China to a global pandemic. Early estimates of the epidemic growth and incubation period of COVID-19 may have been biased due to sample selection. Using detailed case reports from 14 locations in and outside mainland China, we obtained 378 Wuhan-exported cases who left Wuhan before an abrupt travel quarantine. We developed a generative model we call BETS for four key epidemiological events---Beginning of exposure, End of exposure, time of Transmission, and time of Symptom onset (BETS)---and derived explicit formulas to correct for the sample selection. We gave a detailed illustration of why some early and highly influential analyses of the COVID-19 pandemic were severely biased. All our analyses, regardless of which subsample and model were being used, point to an epidemic doubling time of 2 to 2.5 days during the early outbreak in Wuhan. A Bayesian nonparametric analysis further suggests that about 5% of the symptomatic cases may not develop symptoms within 14 days of infection and that men may be much more likely than women to develop symptoms within 2 days of infection.


title: Direct Astrophysical Tests of Chiral Effective Field Theory at Supranuclear Densities

abstract: Recent observations of neutron stars with gravitational waves and X-ray timing provide unprecedented access to the equation of state (EoS) of cold dense matter at densities difficult to realize in terrestrial experiments. At the same time, predictions for the EoS with reliable uncertainty estimates from chiral effective field theory ($\chi$EFT) bound our theoretical ignorance. In this work, we analyze astrophysical data using a nonparametric representation of the neutron-star EoS conditioned on $\chi$EFT to directly constrain the underlying physical properties of the compact objects. We discuss how the data alone constrain the EoS at high densities when we condition on $\chi$EFT at low densities. We also demonstrate how to exploit astrophysical data to directly test the predictions of $\chi$EFT for the EoS up to twice nuclear saturation density, and estimate the density at which these predictions might break down. We find that the existence of massive pulsars, gravitational waves from GW170817, and NICER observations of PSR J0030+0451 favor $\chi$EFT predictions for the EoS up to nuclear saturation density over a more agnostic analysis by as much as a factor of 7 for the quantum Monte Carlo (QMC) calculations used in this work. While $\chi$EFT predictions using QMC are fully consistent with gravitational-wave data up to twice nuclear saturation density, NICER observations suggest that the EoS stiffens relative to these predictions at nuclear saturation density. Additionally, we marginalize over the uncertainty in the density at which $\chi$EFT begins to break down, constraining the radius of a $1.4\,M_\odot$ neutron star to $R_{1.4}=11.40^{+1.38}_{-1.04}$ ($12.54^{+0.71}_{-0.63}$) km and the pressure at twice nuclear saturation density to $p(2n_\mathrm{sat})=14.2^{+18.1}_{-8.4}$ ($28.7^{+15.3}_{-15.0}$) MeV/fm$^3$ with massive pulsar and gravitational-wave (and NICER) data.


title: Blur Aware Calibration of Multi-Focus Plenoptic Camera

abstract: This paper presents a novel calibration algorithm for Multi-Focus Plenoptic Cameras (MFPCs) using raw images only. The design of such cameras is usually complex and relies on precise placement of optic elements. Several calibration procedures have been proposed to retrieve the camera parameters but relying on simplified models, reconstructed images to extract features, or multiple calibrations when several types of micro-lens are used. Considering blur information, we propose a new Blur Aware Plenoptic (BAP) feature. It is first exploited in a pre-calibration step that retrieves initial camera parameters, and secondly to express a new cost function for our single optimization process. The effectiveness of our calibration method is validated by quantitative and qualitative experiments.


title: Magnetic resonance imaging of spin-wave transport and interference in a magnetic insulator

abstract: Spin waves - the elementary excitations of magnetic materials - are prime candidate signal carriers for low dissipation information processing. Being able to image coherent spin-wave transport is crucial for developing interference-based spin-wave devices. We introduce a platform for probing coherent spin waves based on magnetic resonance imaging with electron spins in diamond. Focusing on a thin-film magnetic insulator, we quantify spin-wave amplitudes, visualize the dispersion, and demonstrate time-domain measurements of spin-wave packets. We use our platform to study spin-wave interference, revealing uni-directional, autofocused spin-wave patterns with frequency-controlled numerical apertures. A theoretical analysis explains the patterns in terms of chiral spin-wave excitation and stray-field coupling to the sensor spins. These results pave the way for probing spin waves in atomically thin magnets, even when embedded between opaque materials.


title: Helical velocity fluctuations and large-scale circulation in turbulent thermal convection

abstract: It is shown that helicity dynamics dominates spectral distribution of the velocity fluctuations in turbulent thermal convection at moderate and large values of the Rayleigh number (distributed chaos and scaling respectively). The large-scale circulation (wind) in the turbulent Rayleigh-B\'{e}nard convection is especially sensitive to the helicity dynamics and its chaotic reversals can be associated with the reversals of the mean helicity sign (the ergodic restoration of the reflectional symmetry in a long run). Results of the laboratory experiments and measurements in the free atmospheric convection have been used in order to support the theoretical considerations based on the adiabatic invariance of the second order moment of the helicity distribution. The geomagnetic field dipole variability and reversals have been also discussed in this context and good agreement with results of numerical simulations and observational data has been established.


title: Light-matter interactions with photonic quasiparticles

abstract: Interactions between light and matter play an instrumental role in many fields of science, giving rise to important applications in spectroscopy, sensing, quantum information processing, and lasers. In most of these applications, light is considered in terms of electromagnetic plane waves that propagate at the speed of light in vacuum. As a result, light-matter interactions can usually be treated as very weak, and captured at the lowest order in quantum electrodynamics (QED). However, recent progress in coupling photons to material quasiparticles (e.g., plasmons, phonons, and excitons) forces us to generalize the way we picture the photon at the core of every light-matter interaction. In this new picture, the photon, now of partly matter-character, can have greatly different polarization and dispersion, and be confined to the scale of a few nanometers. Such photonic quasiparticles enable a wealth of light-matter interaction phenomena that could not have been observed before, both in interactions with bound electrons and with free electrons. This Review focuses on exciting theoretical and experimental developments in realizing new light-matter interactions with photonic quasiparticles. As just a few examples, we discuss how photonic quasiparticles enable room-temperature strong coupling, ultrafast "forbidden transitions" in atoms, and new applications of the Cherenkov effect, as well as breakthroughs in ultrafast electron microscopy and new concepts for compact X-ray sources.


title: Removing Algebraic Data Types from Constrained Horn Clauses Using Difference Predicates

abstract: We address the problem of proving the satisfiability of Constrained Horn Clauses (CHCs) with Algebraic Data Types (ADTs), such as lists and trees. We propose a new technique for transforming CHCs with ADTs into CHCs where predicates are defined over basic types, such as integers and booleans, only. Thus, our technique avoids the explicit use of inductive proof rules during satisfiability proofs. The main extension over previous techniques for ADT removal is a new transformation rule, called differential replacement, which allows us to introduce auxiliary predicates corresponding to the lemmas that are often needed when making inductive proofs. We present an algorithm that uses the new rule, together with the traditional folding/unfolding transformation rules, for the automatic removal of ADTs. We prove that if the set of the transformed clauses is satisfiable, then so is the set of the original clauses. By an experimental evaluation, we show that the use of the differential replacement rule significantly improves the effectiveness of ADT removal, and we show that our transformation-based approach is competitive with respect to a well-established technique that extends the CVC4 solver with induction.


title: Extracting the effective contact rate of COVID-19 pandemic

abstract: In the absence of any available vaccines or drugs, prevention of the spread of Coronavirus Disease 2019 (COVID-19) pandemic is being achieved by putting many mitigation measures in place. It is indispensable to have robust and reliable ways of evaluating the effectiveness of these measures. In this work, we assume that, at a very coarse-grained level of description, the overall effect of all the mitigation measures is that we can still describe the spread of the pandemic using the most basic Susceptible-Exposed-Infectious-Removed ($SEIR$) model but with an "effective contact rate" ($\beta$) which is time-dependent. We then use the time series data of the number of infected individuals in the population to extract the instantaneous effective contact rate which is the result of various social interventions put in place. This approach has the potential to be significantly useful while evaluating the impact of mitigation measures on the spread of COVID-19 in near future.


title: Strichartz estimates for mixed homogeneous surfaces in three dimensions

abstract: We obtain sharp mixed norm Strichartz estimates associated to mixed homogeneous surfaces in $\mathbb{R}^3$. Both cases with and without a damping factor are considered. In the case when a damping factor is considered our results yield a wide generalization of a result of Carbery, Kenig, and Ziesler [CKZ13]. The approach we use is to first classify all possible singularities locally, after which one can tackle the problem by appropriately modifying the methods from the paper of Ginibre and Velo [GV92], and by using the recently developed methods by Ikromov and M\"uller [IM16].


title: Higher derivative scalar-tensor theory and spatially covariant gravity: the correspondence

abstract: We investigate the correspondence between generally covariant higher derivative scalar-tensor theory and spatially covariant gravity theory. The building blocks are the scalar field and spacetime curvature tensor together with their generally covariant derivatives for the former, and the spatially covariant geometric quantities together with their spatially covariant derivatives for the later. In the case of a single scalar degree of freedom, they are transformed to each other by gauge fixing and recovering procedures, of which we give the explicit expressions. We make a systematic classification of all the scalar monomials in the spatially covariant gravity according to the total number of derivatives up to $d=4$, and their correspondence to the scalar-tensor monomials. We discusse the possibility of using spatially covariant monomials to generate ghostfree higher derivative scalar-tensor theories. We also derive the covariant 3+1 decomposition without fixing any specific coordinate, which will be useful when performing a covariant Hamiltonian analysis.


title: The Impact of 5G Channel Models on the Performance of Intelligent Reflecting Surfaces and Decode-and-Forward Relaying

abstract: An intelligent reflecting surface (IRS) is an array of discrete elements with configurable scattering properties. It has the capability to beamform arriving radio waves to an intended receiver, making it an attractive candidate technology for fifth-generation (5G) communications. A recent study debunked the notion that IRSs can replace relays because a large number of IRS elements is required even to approach the performance of simple single-antenna decode-and-forward (DF) relays. The study introduced 4G channel models into a theoretical framework to obtain simulation results, based on which comparisons between the two schemes were carried out. In this paper, we consider 5G channel models, reflect on the revised results, and argue that IRSs and DF relays can complement each other's strengths and can both have a place in 5G and beyond 5G architectures.


title: Investigating Efficient Learning and Compositionality in Generative LSTM Networks

abstract: When comparing human with artificial intelligence, one major difference is apparent: Humans can generalize very broadly from sparse data sets because they are able to recombine and reintegrate data components in compositional manners. To investigate differences in efficient learning, Joshua B. Tenenbaum and colleagues developed the character challenge: First an algorithm is trained in generating handwritten characters. In a next step, one version of a new type of character is presented. An efficient learning algorithm is expected to be able to re-generate this new character, to identify similar versions of this character, to generate new variants of it, and to create completely new character types. In the past, the character challenge was only met by complex algorithms that were provided with stochastic primitives. Here, we tackle the challenge without providing primitives. We apply a minimal recurrent neural network (RNN) model with one feedforward layer and one LSTM layer and train it to generate sequential handwritten character trajectories from one-hot encoded inputs. To manage the re-generation of untrained characters, when presented with only one example of them, we introduce a one-shot inference mechanism: the gradient signal is backpropagated to the feedforward layer weights only, leaving the LSTM layer untouched. We show that our model is able to meet the character challenge by recombining previously learned dynamic substructures, which are visible in the hidden LSTM states. Making use of the compositional abilities of RNNs in this way might be an important step towards bridging the gap between human and artificial intelligence.


title: Accelerating complex control schemes on a heterogeneous MPSoC platform for quantum computing

abstract: Control and readout of superconducting quantum bits (qubits) require microwave pulses with gigahertz frequencies and nanosecond precision. To generate and analyze these microwave pulses, we developed a versatile FPGA-based electronics platform. While basic functionality is directly handled within the FPGA, guaranteeing highest accuracy on the nanosecond timescale, more complex control schemes render impractical to implement in hardware. To provide deterministic timing and low latency with high flexibility, we developed the Taskrunner framework. It enables the execution of complex control schemes, so-called user tasks, on the real-time processing unit (RPU) of a heterogeneous Multiprocessor System-on-Chip (MPSoC). These user tasks are specified conveniently using standard C language and are compiled automatically by the MPSoC platform when loaded onto the RPU. We present the architecture of the Taskrunner framework as well as timing benchmarks and discuss applications in the field of quantum computing.


title: Light-Front Holography and Supersymmetric Conformal Algebra: A Novel Approach to Hadron Spectroscopy, Structure, and Dynamics

abstract: We give an overview of recent progress into the infrared structure of QCD based on the gauge/gravity correspondence and light-front quantization, where the color confining interaction for mesons and baryons is determined by an underlying superconformal algebraic structure. This new approach to hadron physics gives remarkable connections and predictions across the entire mass spectrum of hadrons and also describes the infrared behavior of the strong coupling. More recently, an extensive study of form factors, polarized and unpolarized parton distributions and the sea quark contribution to the nucleon has been carried out by extending the holographic formalism to incorporate the nonperturbative structure of Veneziano amplitudes. Contribution to the report "Strong QCD from Hadron Structure Experiments" Jefferson Lab Workshop, 5-9 November, 2019


abstract: In this letter, we present a novel Gaussian Process Learning-based Probabilistic Optimal Power Flow (GP-POPF) for solving POPF under renewable and load uncertainties of arbitrary distribution. The proposed method relies on a non-parametric Bayesian inference-based uncertainty propagation approach, called Gaussian Process (GP). We also suggest a new type of sensitivity called Subspace-wise Sensitivity, using observations on the interpretability of GP-POPF hyperparameters. The simulation results on 14-bus and 30-bus systems show that the proposed method provides reasonably accurate solutions when compared with Monte-Carlo Simulations (MCS) solutions at different levels of uncertain renewable penetration as well as load uncertainties, while requiring much less number of samples and elapsed time.


title: Darboux transformation and multi-soliton solutions of discrete sine-Gordon equation

abstract: We study a discrete Darboux transformation and construct the multi-soliton solutions in terms of ratio of determinants for integrable discrete sine-Gordon equation. We also calculate explicit expressions of single, double, triple, quad soliton solutions as well as single and double breather solutions of discrete sine-Gordon equation. Dynamical features of discrete kinks and breathers have also been illustrated.


title: Spectral Index of OAM-Carrying Ultrafast Localized Pulses

abstract: We investigate the spectral degree of freedom of OAM-carrying localized waves, and its influence on their transverse intensity distribution. In particular, we focus our attention on two different families of spectra, namely exponentially decaying spectra, which are very tightly connected to fundamental X-waves, and Bessel-modulated spectra. For each class we show how it is possible, by suitable manipulating their spectrum, to structure their transverse intensity distribution, thus creating a radial structure similar to that of Laguerre-Gaussian beams. To complete our analysis, we investigate the impact of the two spectral functions on linear and angular momentum of such localized waves.


title: Decoding Probability Analysis of Network-Coded Data Collection and Delivery by Relay Drones

abstract: Relay drones in delay-tolerant applications are dispatched to remote locations in order to gather data transmitted by a source node. Collected data are stored on the drones and delivered to one or multiple bases. This paper considers two schemes for broadcasting data to drones when feedback channels are not available: a data carousel and systematic random linear network coding (RLNC). We propose a theoretical framework for the calculation of the probability that a base will fully or partially recover the transmitted data and the probability that all involved bases will successfully obtain the data, when the bases are either isolated or interconnected. Theoretical results are validated through simulations. Design considerations are also discussed, including the relationship among the field size used by RLNC, the number of relay drones and the requirement for full data recovery or the retrieval of at least part of the data.


title: Deep Generation of Coq Lemma Names Using Elaborated Terms

abstract: Coding conventions for naming, spacing, and other essentially stylistic properties are necessary for developers to effectively understand, review, and modify source code in large software projects. Consistent conventions in verification projects based on proof assistants, such as Coq, increase in importance as projects grow in size and scope. While conventions can be documented and enforced manually at high cost, emerging approaches automatically learn and suggest idiomatic names in Java-like languages by applying statistical language models on large code corpora. However, due to its powerful language extension facilities and fusion of type checking and computation, Coq is a challenging target for automated learning techniques. We present novel generation models for learning and suggesting lemma names for Coq projects. Our models, based on multi-input neural networks, are the first to leverage syntactic and semantic information from Coq's lexer (tokens in lemma statements), parser (syntax trees), and kernel (elaborated terms) for naming; the key insight is that learning from elaborated terms can substantially boost model performance. We implemented our models in a toolchain, dubbed Roosterize, and applied it on a large corpus of code derived from the Mathematical Components family of projects, known for its stringent coding conventions. Our results show that Roosterize substantially outperforms baselines for suggesting lemma names, highlighting the importance of using multi-input models and elaborated terms.


title: Adaptive BEM for elliptic PDE systems, Part I: Abstract framework for weakly-singular integral equations

abstract: In the present work, we consider weakly-singular integral equations arising from linear second-order strongly-elliptic PDE systems with constant coefficients, including, e.g., linear elasticity. We introduce a general framework for optimal convergence of adaptive Galerkin BEM. We identify certain abstract properties for the underlying meshes, the corresponding mesh-refinement strategy, and the ansatz spaces that guarantee convergence at optimal algebraic rate of an adaptive algorithm driven by the weighted-residual error. These properties are satisfied, e.g., for discontinuous piecewise polynomials on simplicial meshes as well as certain ansatz spaces used for isogeometric analysis. Technical contributions include local inverse estimates for the (non-local) boundary integral operators associated to the PDE system.


title: A Magnon Scattering Platform

abstract: Scattering experiments have revolutionized our understanding of nature. Examples include the discovery of the nucleus, crystallography, and the discovery of the double helix structure of DNA. Scattering techniques differ by the type of the particles used, the interaction these particles have with target materials and the range of wavelengths used. Here, we demonstrate a new 2-dimensional table-top scattering platform for exploring magnetic properties of materials on mesoscopic length scales. Long lived, coherent magnonic excitations are generated in a thin film of YIG and scattered off a magnetic target deposited on its surface. The scattered waves are then recorded using a scanning NV center magnetometer that allows sub-wavelength imaging and operation under conditions ranging from cryogenic to ambient environment. While most scattering platforms measure only the intensity of the scattered waves, our imaging method allows for spatial determination of both amplitude and phase of the scattered waves thereby allowing for a systematic reconstruction of the target scattering potential. Our experimental results are consistent with theoretical predictions for such a geometry and reveal several unusual features of the magnetic response of the target, including suppression near the target edges and gradient in the direction perpendicular to the direction of surface wave propagation. Our results establish magnon scattering experiments as a new platform for studying correlated many-body systems.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el cond-mat.supr-con physics.app-ph


title: Quantum $\mathcal{R}$-matrices as universal qubit gates

abstract: We study the Chern-Simons approach to the topological quantum computing. We use quantum $\mathcal{R}$-matrices as universal quantum gates and study the approximations of some one-qubit operations. We make some modifications to the known Solovay-Kitaev algorithm suitable for our particular problem.


title: Cost-effectiveness Analysis of Antiepidemic Policies and Global Situation Assessment of COVID-19

abstract: With a two-layer contact-dispersion model and data in China, we analyze the cost-effectiveness of three types of antiepidemic measures for COVID-19: regular epidemiological control, local social interaction control, and inter-city travel restriction. We find that: 1) intercity travel restriction has minimal or even negative effect compared to the other two at the national level; 2) the time of reaching turning point is independent of the current number of cases, and only related to the enforcement stringency of epidemiological control and social interaction control measures; 3) strong enforcement at the early stage is the only opportunity to maximize both antiepidemic effectiveness and cost-effectiveness; 4) mediocre stringency of social interaction measures is the worst choice. Subsequently, we cluster countries/regions into four groups based on their control measures and provide situation assessment and policy suggestions for each group.


title: A Canonical Polynomial Van der Waerden's Theorem

abstract: We prove a canonical polynomial Van der Waerden's Theorem. More precisely, we show the following. Let $\{p_1(x),\ldots,p_k(x)\}$ be a set of polynomials such that $p_i(x)\in \mathbb{Z}[x]$ and $p_i(0)=0$, for every $i\in \{1,\ldots,k\}$. Then, in any colouring of $\mathbb{Z}$, there exist $a,d\in \mathbb{Z}$ such that $\{a+p_1(d),\ldots,a+p_{k}(d)\}$ forms either a monochromatic or a rainbow set.


title: Sliding Window Polar Codes

abstract: We propose a novel coupling technique for the design of polar codes of length N, making them decodable through a sliding window of size M < N. This feature allows to reduce the computational complexity of the decoder, an important possibility in wireless communication downlink scenarios. Our approach is based on the design of an ad-hoc kernel to be inserted in a multi-kernel polar code framework; this structure enables the sliding window decoding of the code. Simulation results show that the proposed sliding window polar codes outperform the independent blocks transmission in the proposed scenario, at the cost of a negligible decoding overhead.


title: HST/FGS Trigonometric Parallaxes of M-dwarf Eclipsing Binaries

abstract: Hubble Space Telescope (HST) Fine Guidance Sensor (FGS) trigonometric parallax observations were obtained to directly determine distances to five nearby M-dwarf / M-dwarf eclipsing binary systems. These systems are intrinsically interesting as benchmark systems for establishing basic physical parameters for low-mass stars, such as luminosity L, and radius R. HST/FGS distances are also one of the few direct checks on Gaia trigonometric parallaxes, given the comparable sensitivity in both magnitude limit and determination of parallactic angles. A spectral energy distribution (SED) fit of each system's blended flux output was carried out, allowing for estimation of the bolometric flux from the primary and secondary components of each system. From the stellar M, L, and R values, the low-mass star relationships between L and M, and R and M, are compared against idealized expectations for such stars. An examination on the inclusion of these close M-dwarf/M-dwarf pairs in higher-order common proper motion (CPM) pairs is analysed; each of the 5 systems has indications of being part of a CPM system. Unexpected distances on interesting objects found within the grid of parallactic reference stars are also presented, including a nearby M dwarf and a white dwarf.


title: SCOUT: Self-aware Discriminant Counterfactual Explanations

abstract: The problem of counterfactual visual explanations is considered. A new family of discriminant explanations is introduced. These produce heatmaps that attribute high scores to image regions informative of a classifier prediction but not of a counter class. They connect attributive explanations, which are based on a single heat map, to counterfactual explanations, which account for both predicted class and counter class. The latter are shown to be computable by combination of two discriminant explanations, with reversed class pairs. It is argued that self-awareness, namely the ability to produce classification confidence scores, is important for the computation of discriminant explanations, which seek to identify regions where it is easy to discriminate between prediction and counter class. This suggests the computation of discriminant explanations by the combination of three attribution maps. The resulting counterfactual explanations are optimization free and thus much faster than previous methods. To address the difficulty of their evaluation, a proxy task and set of quantitative metrics are also proposed. Experiments under this protocol show that the proposed counterfactual explanations outperform the state of the art while achieving much higher speeds, for popular networks. In a human-learning machine teaching experiment, they are also shown to improve mean student accuracy from chance level to 95\%.


title: Reinforcement learning approach to non-equilibrium quantum thermodynamics

abstract: We use a reinforcement learning approach to reduce entropy production in a closed quantum system brought out of equilibrium. Our strategy makes use of an external control Hamiltonian and a policy gradient technique. Our approach bears no dependence on the quantitative tool chosen to characterize the degree of thermodynamic irreversibility induced by the dynamical process being considered, require little knowledge of the dynamics itself and does not need the tracking of the quantum state of the system during the evolution, thus embodying an experimentally non-demanding approach to the control of non-equilibrium quantum thermodynamics. We successfully apply our methods to the case of single- and two-particle systems subjected to time-dependent driving potentials.


title: Quantum Gravity, Constant Negative Curvatures, and Black Holes

abstract: For purposes of quantization, classical gravity is normally expressed by canonical variables, namely the metric $g_{ab}(x)$ and the momentum $\pi^{cd}(x)$. Canonical quantization requires a proper promotion of these classical variables to quantum operators, which, according to Dirac, the favored operators should be those arising from classical variables that formed Cartesian coordinates; sadly, in this case, that is not possible. However, an affine quantization features promoting the metric $g_{ab}(x)$ and the momentric $\pi^c_d(x)\;[\equiv \pi^{ce}(x) \,g_{de}(x)]$ to operators. Instead of these classical variables belonging to a constant zero curvature space (i.e., instead of a flat space), they belong to a space of constant negative curvatures. This feature may even have its appearance in black holes, which could strongly point toward an affine quantization approach to quantize gravity.


title: QZE and IZE in a simple approach and the neutron decay

abstract: We discuss a simple and analytically solvable measurement model which describes the famous Quantum Zeno Effect (QZE) and Inverse Zeno Effect (IZE), that correspond to the slow down and to the increase of the decay rate caused by measurements (or, more in general, by the interaction of an unstable state with the detector and the environment). Within this model one can understand quite general features of the QZE and IZE: by considering an unstable quantum state, such as an unstable particle, whose decay width as function of energy is $\Gamma(\omega)=g^{2}\omega^{\alpha},$ then -- under quite general assumptions -- the QZE occurs for $\alpha\in(0,1)$, while the IZE for $\alpha \in(-\infty,0)\cup(1,\infty).$ This result is valid also for more realistic measurement models than the one described in this work. We then apply these considerations to the decay of the neutron, for which $\alpha=5.$ Hence, the realization of the IZE for the neutron decay (and for the majority of weak decays) is in principle possible. Indeed, trap experiments find a lifetime that is $8.7\pm2.1$ s shorter than beam experiments, suggesting that the IZE could have taken place.


title: Accelerated multimodel Newton-type algorithms for faster convergence of ground and excited state coupled cluster equations

abstract: We introduce a multimodel approach to solve coupled cluster equations, employing a quasi Newton algorithm for the ground state and an Olsen algorithm for the excited states. In these algorithms, both of which can be viewed as Newton algorithms, the Jacobian matrix of a lower level coupled cluster model is used in Newton equations associated with the target model. Improvements in convergence then implies savings for sufficiently large molecular systems, since the computational cost of macroiterations scales more steeply with system size than the cost of microiterations. The multimodel approach is suitable when there is a lower level Jacobian matrix that is much more accurate than the zeroth order approximation. Applying the approach to the CC3 equations, using the CCSD approximation of the Jacobian, we show that the time spent to determine the ground and valence excited states can be significantly reduced. We also find improved convergence for core excited states, indicating that similar savings will be obtained with an explicit implementation of the core-valence separated CCSD Jacobian transformation.


title: Computing all identifiable functions of parameters for ODE models

abstract: Parameter identifiability is a structural property of an ODE model for recovering the values of parameters from the data (i.e., from the input and output variables). This property is a prerequisite for meaningful parameter identification in practice. In the presence of nonidentifiability, it is important to find all functions of the parameters that are identifiable. The existing algorithms check whether a given function of parameters is identifiable or, under the solvability condition, find all identifiable functions. However, this solvability condition is not always satisfied, which presents a challenge. Our first main result is an algorithm that computes all identifiable functions without any additional assumptions, which is the first such algorithm as far as we know. Our second main result concerns the identifiability from multiple experiments (with generically different inputs and initial conditions among the experiments). For this problem, we prove that the set of functions identifiable from multiple experiments is what would actually be computed by input-output equation-based algorithms (whether or not the solvability condition is fulfilled), which was not known before. We give an algorithm that not only finds these functions but also provides an upper bound for the number of experiments to be performed to identify these functions. We provide an implementation of the presented algorithms.


title: Quantification of plasticity via particle dynamics above and below yield in a 2D jammed suspension

abstract: Failure of amorphous materials is characterized by the emergence of dissipation. The connection between particle dynamics, dissipation, and overall material rheology, however, has still not been elucidated. Here, we take a new approach relating trajectories to yielding, using a custom built interfacial stress rheometer, which allows for measurement of shear moduli (G',G'') of a dense athermal suspension's microstructure while simultaneously tracking particle trajectories undergoing cyclic shear. We find an increase in total area traced by particle trajectories as the system is stressed well below to well above yield. Trajectories may be placed into three categories: reversibly elastic paths; reversibly plastic paths, associated with smooth limit cycles; and irreversibly plastic paths, in which particles do not return to their original position. We find that above yield, reversibly plastic trajectories are predominantly found near to the shearing surface, whereas reversibly elastic paths are more prominent near the stationary wall. This spatial transition between particles acting as solids to those acting as liquids is characteristic of a 'melting front', which is observed to shift closer to the wall with increasing strain. We introduce a non-dimensional measure of plastic dissipation based on particle trajectories that scales linearly with strain amplitude both above and below yield, and that is unity at the rheological yield point. Surprisingly, this relation collapses for three systems of varying degrees of disorder.


title: Kvistur 2.0: a BiLSTM Compound Splitter for Icelandic

abstract: In this paper, we present a character-based BiLSTM model for splitting Icelandic compound words, and show how varying amounts of training data affects the performance of the model. Compounding is highly productive in Icelandic, and new compounds are constantly being created. This results in a large number of out-of-vocabulary (OOV) words, negatively impacting the performance of many NLP tools. Our model is trained on a dataset of 2.9 million unique word forms and their constituent structures from the Database of Icelandic Morphology. The model learns how to split compound words into two parts and can be used to derive the constituent structure of any word form. Knowing the constituent structure of a word form makes it possible to generate the optimal split for a given task, e.g., a full split for subword tokenization, or, in the case of part-of-speech tagging, splitting an OOV word until the largest known morphological head is found. The model outperforms other previously published methods when evaluated on a corpus of manually split word forms. This method has been integrated into Kvistur, an Icelandic compound word analyzer.


title: Gaze-Net: Appearance-Based Gaze Estimation using Capsule Networks

abstract: Recent studies on appearance based gaze estimation indicate the ability of Neural Networks to decode gaze information from facial images encompassing pose information. In this paper, we propose Gaze-Net: A capsule network capable of decoding, representing, and estimating gaze information from ocular region images. We evaluate our proposed system using two publicly available datasets, MPIIGaze (200,000+ images in the wild) and Columbia Gaze (5000+ images of users with 21 gaze directions observed at 5 camera angles/positions). Our model achieves a Mean Absolute Error (MAE) of 2.84$^\circ$ for Combined angle error estimate within dataset for MPI-IGaze dataset. Further, model achieves a MAE of 10.04$^\circ$ for across dataset gaze estimation error for Columbia gaze dataset. Through transfer learning, the error is reduced to 5.9$^\circ$. The results show this approach is promising with implications towards using commodity webcams to develop low-cost multi-user gaze tracking systems.


title: Privacy-Preserving Policy Synthesis in Markov Decision Processes

abstract: In decision-making problems, the actions of an agent may reveal sensitive information that drives its decisions. For instance, a corporation's investment decisions may reveal its sensitive knowledge about market dynamics. To prevent this type of information leakage, we introduce a policy synthesis algorithm that protects the privacy of the transition probabilities in a Markov decision process. We use differential privacy as the mathematical definition of privacy. The algorithm first perturbs the transition probabilities using a mechanism that provides differential privacy. Then, based on the privatized transition probabilities, we synthesize a policy using dynamic programming. Our main contribution is to bound the "cost of privacy," i.e., the difference between the expected total rewards with privacy and the expected total rewards without privacy. We also show that computing the cost of privacy has time complexity that is polynomial in the parameters of the problem. Moreover, we establish that the cost of privacy increases with the strength of differential privacy protections, and we quantify this increase. Finally, numerical experiments on two example environments validate the established relationship between the cost of privacy and the strength of data privacy protections.


title: Grothendieck-Witt groups of some singular schemes

abstract: We establish some structural results for the Witt and Grothendieck-Witt groups of schemes over $\mathbb{Z}[1/2]$, including homotopy invariance for Witt groups and a formula for the Witt and Grothendieck-Witt groups of punctured affine spaces over a scheme. All these results hold for singular schemes and at the level of spectra.


title: Shortcut Learning in Deep Neural Networks

abstract: Deep learning has triggered the current rise of artificial intelligence and is the workhorse of today's machine intelligence. Numerous success stories have rapidly spread all over science, industry and society, but its limitations have only recently come into focus. In this perspective we seek to distill how many of deep learning's problems can be seen as different symptoms of the same underlying problem: shortcut learning. Shortcuts are decision rules that perform well on standard benchmarks but fail to transfer to more challenging testing conditions, such as real-world scenarios. Related issues are known in Comparative Psychology, Education and Linguistics, suggesting that shortcut learning may be a common characteristic of learning systems, biological and artificial alike. Based on these observations, we develop a set of recommendations for model interpretation and benchmarking, highlighting recent advances in machine learning to improve robustness and transferability from the lab to real-world applications.


title: Finite Sample Theory for High-Dimensional Functional/Scalar Time Series with Applications

abstract: Statistical analysis of high-dimensional functional times series arises in various applications. Under this scenario, in addition to the intrinsic infinite-dimensionality of functional data, the number of functional variables can grow with the number of serially dependent observations. In this paper, we focus on the theoretical analysis of relevant estimated cross-(auto)covariance terms between two multivariate functional time series or a mixture of multivariate functional and scalar time series beyond the Gaussianity assumption. We introduce a new perspective on dependence by proposing functional cross-spectral stability measure to characterize the effect of dependence on these estimated cross terms, which are essential in the estimates for additive functional linear regressions. With the proposed functional cross-spectral stability measure, we develop useful concentration inequalities for estimated cross-(auto)covariance matrix functions to accommodate more general sub-Gaussian functional linear processes and, furthermore, establish finite sample theory for relevant estimated terms under a commonly adopted functional principal component analysis framework. Using our derived non-asymptotic results, we investigate the convergence properties of the regularized estimates for two additive functional linear regression applications under sparsity assumptions including functional linear lagged regression and partially functional linear regression in the context of high-dimensional functional/scalar time series.


title: Network-principled deep generative models for designing drug combinations as graph sets

abstract: Combination therapy has shown to improve therapeutic efficacy while reducing side effects. Importantly, it has become an indispensable strategy to overcome resistance in antibiotics, anti-microbials, and anti-cancer drugs. Facing enormous chemical space and unclear design principles for small-molecule combinations, the computational drug-combination design has not seen generative models to meet its potential to accelerate resistance-overcoming drug combination discovery. We have developed the first deep generative model for drug combination design, by jointly embedding graph-structured domain knowledge and iteratively training a reinforcement learning-based chemical graph-set designer. First, we have developed Hierarchical Variational Graph Auto-Encoders (HVGAE) trained end-to-end to jointly embed gene-gene, gene-disease, and disease-disease networks. Novel attentional pooling is introduced here for learning disease-representations from associated genes' representations. Second, targeting diseases in learned representations, we have recast the drug-combination design problem as graph-set generation and developed a deep learning-based model with novel rewards. Specifically, besides chemical validity rewards, we have introduced a novel generative adversarial award, being generalized sliced Wasserstein, for chemically diverse molecules with distributions similar to known drugs. We have also designed a network principle-based reward for drug combinations. Numerical results indicate that, compared to graph embedding methods, HVGAE learns more informative and generalizable disease representations. Case studies on four diseases show that network-principled drug combinations tend to have low toxicity. The generated drug combinations collectively cover the disease module similar to FDA-approved drug combinations and could potentially suggest novel systems-pharmacology strategies.


title: TOI-1338: TESS' First Transiting Circumbinary Planet

abstract: We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars.


title: Stability and instability issues of the Weinstock inequality

abstract: Given two planar, conformal, smooth open sets $\Omega$ and $\omega$, we prove the existence of a sequence of smooth sets $\Omega_n$ which geometrically converges to $\Omega$ and such that the (perimeter normalized) Steklov eigenvalues of $\Omega_n$ converge to the ones of $\omega$. As a consequence, we answer a question raised by Girouard and Polterovich on the stability of the Weinstock inequality and prove that the inequality is genuinely unstable. However, under some a priori knowledge of the geometry related to the oscillations of the boundaries, stability may occur.


title: Quantum erasing of laser emission in N$_2^+$

abstract: Cavity-free lasing of N$_2^+$ induced by a femtosecond laser pulse at 800~nm is nearly totally suppressed by a delayed twin pump pulse. We explain this surprising effect within the V-scheme of lasing without population inversion. A fast transfer of population between nitrogen ionic states X$^2\Sigma_g^+$ and A$^2\Pi_u$, induced by the second pulse, terminates the conditions for amplification in the system. Appearance of short lasing bursts at delays corresponding to revivals of rotational wave packets are explained along the same lines.


title: Multi-Object Tracking with Siamese Track-RCNN

abstract: Multi-object tracking systems often consist of a combination of a detector, a short term linker, a re-identification feature extractor and a solver that takes the output from these separate components and makes a final prediction. Differently, this work aims to unify all these in a single tracking system. Towards this, we propose Siamese Track-RCNN, a two stage detect-and-track framework which consists of three functional branches: (1) the detection branch localizes object instances; (2) the Siamese-based track branch estimates the object motion and (3) the object re-identification branch re-activates the previously terminated tracks when they re-emerge. We test our tracking system on two popular datasets of the MOTChallenge. Siamese Track-RCNN achieves significantly higher results than the state-of-the-art, while also being much more efficient, thanks to its unified design.


title: Bifurcation of perturbations of non-generic closed self-shrinkers

abstract: We discover a bifurcation of the perturbations of non-generic closed self-shrinkers. If the generic perturbation is outward, then the next mean curvature flow singularity is cylindrical and collapsing from outside; if the generic perturbation is inward, then the next mean curvature flow singularity is cylindrical and collapsing from inside.


title: RGBD-Dog: Predicting Canine Pose from RGBD Sensors

abstract: The automatic extraction of animal \reb{3D} pose from images without markers is of interest in a range of scientific fields. Most work to date predicts animal pose from RGB images, based on 2D labelling of joint positions. However, due to the difficult nature of obtaining training data, no ground truth dataset of 3D animal motion is available to quantitatively evaluate these approaches. In addition, a lack of 3D animal pose data also makes it difficult to train 3D pose-prediction methods in a similar manner to the popular field of body-pose prediction. In our work, we focus on the problem of 3D canine pose estimation from RGBD images, recording a diverse range of dog breeds with several Microsoft Kinect v2s, simultaneously obtaining the 3D ground truth skeleton via a motion capture system. We generate a dataset of synthetic RGBD images from this data. A stacked hourglass network is trained to predict 3D joint locations, which is then constrained using prior models of shape and pose. We evaluate our model on both synthetic and real RGBD images and compare our results to previously published work fitting canine models to images. Finally, despite our training set consisting only of dog data, visual inspection implies that our network can produce good predictions for images of other quadrupeds -- e.g. horses or cats -- when their pose is similar to that contained in our training set.


title: The onset of instability in resonant chains

abstract: There is evidence that most chains of mean motion resonances of type $k$:$k-1$ among exoplanets become unstable once the dissipative action from the gas is removed from the system, particularly for large $N$ (the number of planets) and $k$ (indicating how compact the chain is). We present a novel dynamical mechanism that can explain the origin of these instabilities and thus the dearth of resonant systems in the exoplanet sample. It relies on the emergence of secondary resonances between a fraction of the synodic frequency $2 \pi (1/P_1-1/P_2)$ and the libration frequencies in the mean motion resonance. These secondary resonances excite the amplitudes of libration of the mean motion resonances thus leading to an instability. We detail the emergence of these secondary resonances by carrying out an explicit perturbative scheme to second order in the planetary masses and isolating the harmonic terms that are associated with them. Focusing on the case of three planets in the 3:2 -- 3:2 mean motion resonance as an example, a simple but general analytical model of one of these resonances is obtained which describes the initial phase of the activation of one such secondary resonance. The dynamics of the excited system is also briefly described. This scheme shows how one can obtain analytical insight into the emergence of these resonances, and into the dynamics that they trigger. Finally, a generalisation of this dynamical mechanism is obtained for arbitrary $N$ and $k$. This leads to an explanation of previous numerical experiments on the stability of resonant chains, showing why the critical planetary mass allowed for stability decreases with increasing $N$ and $k$.


title: Avoiding the Hypothesis-Only Bias in Natural Language Inference via Ensemble Adversarial Training

abstract: Natural Language Inference (NLI) datasets contain annotation artefacts resulting in spurious correlations between the natural language utterances and their respective entailment classes. These artefacts are exploited by neural networks even when only considering the hypothesis and ignoring the premise, leading to unwanted biases. Belinkov et al. (2019b) proposed tackling this problem via adversarial training, but this can lead to learned sentence representations that still suffer from the same biases. We show that the bias can be reduced in the sentence representations by using an ensemble of adversaries, encouraging the model to jointly decrease the accuracy of these different adversaries while fitting the data. This approach produces more robust NLI models, outperforming previous de-biasing efforts when generalised to 12 other datasets (Belinkov et al., 2019a; Mahabadi et al., 2020). In addition, we find that the optimal number of adversarial classifiers depends on the dimensionality of the sentence representations, with larger sentence representations being more difficult to de-bias while benefiting from using a greater number of adversaries.


title: Investigation of background electron emission in the LUX detector

abstract: Dual-phase xenon detectors, as currently used in direct detection dark matter experiments, have observed elevated rates of background electron events in the low energy region. While this background negatively impacts detector performance in various ways, its origins have only been partially studied. In this paper we report a systematic investigation of the electron pathologies observed in the LUX dark matter experiment. We characterize different electron populations based on their emission intensities and their correlations with preceding energy depositions in the detector. By studying the background under different experimental conditions, we identified the leading emission mechanisms, including photoionization and the photoelectric effect induced by the xenon luminescence, delayed emission of electrons trapped under the liquid surface, capture and release of drifting electrons by impurities, and grid electron emission. We discuss how these backgrounds can be mitigated in LUX and future xenon-based dark matter experiments.


title: Scalar Field Wave (Fuzzy) Dark Matter and the Formation of Galaxies

abstract: We investigate Scalar Field Wave Dark Matter in the context of galactic Dark Matter halos. In particular, we offer an analysis of the Baryonic Tully-Fisher Relation (BTFR). We detail a particular family of excited state solutions to the Einstein-Klein-Gordon equations, and use it to provide a novel theoretical model for producing the BTFR. We then solve this model computationally to simulate the BTFR. Interpreting the Dark Matter scalar field as an ultralight boson, this returns a conservative mass constraint of $m\geq 10^{-23}eV$. Assuming slightly stronger conditions suggests $m\geq 10^{-22}eV$ to be more compatible with the BTFR. We provide a discussion of Scalar Field Dark Matter rotation curves and the structure of Scalar Field Dark Matter halos. Compatibility with the BTFR requires the excited state solutions to obey particular boundary conditions; this may have implications for the behavior of Dark Matter halos and the formation of galaxies.


title: Trajectory Planning and Control for Automatic Docking of ASVs with Full-Scale Experiments

abstract: We propose a method for performing automatic docking of a small autonomous surface vehicle (ASV) by interconnecting an optimization-based trajectory planner with a dynamic positioning (DP) controller for trajectory tracking. The trajectory planner provides collision-free trajectories by considering a map with static obstacles, and produces feasible trajectories through inclusion of a mathematical model of the ASV and its actuators. The DP controller tracks the time-parametrized position, velocity and acceleration produced by the trajectory planner using proportional-integral-derivative feedback with velocity and acceleration feed forward. The method's performance is tested on a small ASV in confined waters in Trondheim, Norway. The ASV performs collision-free docking maneuvers with respect to static obstacles when tracking the generated reference trajectories and achieves successful docking.


title: Global Existence of Non-cutoff Boltzmann Equation in Weighted Sobolev Space

abstract: This article presents a new approach of semigroup analysis and pseudo-differential calculus for deriving the regularizing estimate on non-cutoff linearized Boltzmann equation. We are able to obtain regularizing estimate of semigroup $e^{tB}$ that is continuous from weighted Sobolev space $H(a^{-1/2})H^m_x$ to $H(a^{1/2})H^m_x$ with a sharp large time decay. With these properties, we prove the existence of global-in-time unique solution to the non-cutoff Boltzmann equation for hard potential on the whole space with weak regularity assumption on initial data. We consider the hard potential case since $H(a^{1/2})$ can be embedded in $L^2$. This work develops the application of pseudo-differential calculus, spectrum analysis and semigroup theory to non-cutoff Boltzmann equation.


title: Comment on "The generalized Boltzmann distribution is the only distribution in which the Gibbs-Shannon entropy equals the thermodynamic entropy" by X. Gao, E. Gallicchio and A.E. Roitberg [J. Chem. Phys. 151, 034113 (2019)]

abstract: The title of the paper leads to an incorrect conclusion as we show that the equilibrium result of the paper is a special limit of a general result for nonequilibrium systems in internal equilibrium already available in the literature. We also point out some of the limitations of the approach taken by the authors.


title: Detecting many-body Bell non-locality by solving Ising models

abstract: Bell non-locality represents the ultimate consequence of quantum entanglement, fundamentally undermining the classical tenet that spatially-separated degrees of freedom possess objective attributes independently of the act of their measurement. Despite its importance, probing Bell non-locality in many-body systems is considered to be a formidable challenge, with a computational cost scaling exponentially with system size. Here we propose and validate an efficient variational scheme, based on the solution of inverse classical Ising problems, which in polynomial time can probe whether an arbitrary set of quantum data is compatible with a local theory; and, if not, it delivers a many-body Bell inequality violated by the quantum data. We use our approach to unveil new many-body Bell inequalities, violated by suitable measurements on paradigmatic quantum states (the low-energy states of Heisenberg antiferromagnets), paving the way to systematic Bell tests in the many-body realm.


title: Non-Gaussian correlations imprinted by local dephasing in fermionic wires

abstract: We study the behavior of an extended fermionic wire coupled to a local stochastic field. Since the quantum jump operator is Hermitian and quadratic in fermionic operators, it renders the model soluble, allowing investigation of the properties of the non-equilibrium steady-state and the role of dissipation-induced fluctuations. We derive a closed set of equations of motion solely for the two-point correlator; on the other hand, we find, surprisingly, that the many-body state exhibits non-Gaussian correlations. Density-density correlation function demonstrates a crossover from a regime of weak dissipation characterized by moderate heating and stimulated fluctuations to a quantum Zeno regime ruled by strong dissipation, which tames quantum fluctuations. Instances of soluble dissipative impurities represent an experimentally viable platform to understand the interplay between dissipation and Hamiltonian dynamics in many-body quantum systems.


title: Extending the Reach of the Point-to-Set Principle

abstract: The point-to-set principle of J. Lutz and N. Lutz (2018) has recently enabled the theory of computing to be used to answer open questions about fractal geometry in Euclidean spaces $\mathbb{R}^n$. These are classical questions, meaning that their statements do not involve computation or related aspects of logic. In this paper we extend the reach of the point-to-set principle from Euclidean spaces to arbitrary separable metric spaces $X$. We first extend two fractal dimensions--computability-theoretic versions of classical Hausdorff and packing dimensions that assign dimensions $\dim(x)$ and $\textrm{Dim}(x)$ to individual points $x\in X$--to arbitrary separable metric spaces and to arbitrary gauge families. Our first two main results then extend the point-to-set principle to arbitrary separable metric spaces and to a large class of gauge families. We demonstrate the power of our extended point-to-set principle by using it to prove new theorems about classical fractal dimensions in hyperspaces. (For a concrete computational example, the stages $E_0, E_1, E_2, \ldots$ used to construct a self-similar fractal $E$ in the plane are elements of the hyperspace of the plane, and they converge to $E$ in the hyperspace.) Our third main result, proven via our extended point-to-set principle, states that, under a wide variety of gauge families, the classical packing dimension agrees with the classical upper Minkowski dimension on all hyperspaces of compact sets. We use this theorem to give, for all sets $E$ that are analytic, i.e., $\mathbf{\Sigma}^1_1$, a tight bound on the packing dimension of the hyperspace of $E$ in terms of the packing dimension of $E$ itself.


title: Lorentz-violating scenarios in a thermal reservoir

abstract: In this work, we analyze the thermodynamic properties of the graviton and the generalized model involving anisotropic Podolsky and Lee$-$Wick terms with Lorentz violation. We build up the so-called partition function from the accessible states of the system seeking the following thermodynamic functions: spectral radiance, mean energy, Helmholtz free energy, entropy and heat capacity. Besides, we verify that, when the temperature rises, the spectral radiance $\chi(\nu)$ tends to attenuate for fixed values of $\xi$. Notably, when parameter $\eta_{2}$ increases, the spectral radiance $\bar{\chi}(\eta_{2},\nu)$ weakens until reaching a flat characteristic. Finally, for both theories, we perform the calculation of the modified black body radiation and the correction to the $\textit{Stefan-Boltzmann}$ law in the inflationary era of the universe.


title: Leveraging GANs to Improve Continuous Path Keyboard Input Models

abstract: Continuous path keyboard input has higher inherent ambiguity than standard tapping, because the path trace may exhibit not only local overshoots/undershoots (as in tapping) but also, depending on the user, substantial mid-path excursions. Deploying a robust solution thus requires a large amount of high-quality training data, which is difficult to collect/annotate. In this work, we address this challenge by using GANs to augment our training corpus with user-realistic synthetic data. Experiments show that, even though GAN-generated data does not capture all the characteristics of real user data, it still provides a substantial boost in accuracy at a 5:1 GAN-to-real ratio. GANs therefore inject more robustness in the model through greatly increased word coverage and path diversity.


title: $L^p-L^q$ boundedness of Fourier multipliers associated with the anharmonic Oscillator

abstract: In this paper we study the $L^p$-$L^q$ boundedness of the Fourier multipliers in the setting where the underlying Fourier analysis is introduced with respect to the eigenfunctions of an anharmonic oscillator $A$. Using the notion of a global symbol that arises from this analysis, we extend a version of the Hausdorff-Young-Paley inequality that guarantees the $L^p$-$L^q$ boundedness of these operators for the range $1<p \leq 2 \leq q <\infty$. The boundedness results for spectral multipliers acquired, yield as particular cases Sobolev embedding theorems and time asymptotics for the $L^p$-$L^q$ norms of the heat kernel associated with the anharmonic oscillator. Additionally, we consider functions $f(A)$ of the anharmonic oscillator on modulation spaces and prove that Linsk\u ii's trace formula holds true even when $f(A)$ is simply a nuclear operator.


title: Geometry-Aware Gradient Algorithms for Neural Architecture Search

abstract: Recent state-of-the-art methods for neural architecture search (NAS) exploit gradient-based optimization by relaxing the problem into continuous optimization over architectures and shared-weights, a noisy process that remains poorly understood. We argue for the study of single-level empirical risk minimization to understand NAS with weight-sharing, reducing the design of NAS methods to devising optimizers and regularizers that can quickly obtain high-quality solutions to this problem. Invoking the theory of mirror descent, we present a geometry-aware framework that exploits the underlying structure of this optimization to return sparse architectural parameters, leading to simple yet novel algorithms that enjoy fast convergence guarantees and achieve state-of-the-art accuracy on the latest NAS benchmarks in computer vision. Notably, we exceed the best published results for both CIFAR and ImageNet on both the DARTS search space and NAS-Bench201; on the latter we achieve near-oracle-optimal performance on CIFAR-10 and CIFAR-100. Together, our theory and experiments demonstrate a principled way to co-design optimizers and continuous relaxations of discrete NAS search spaces.


title: Refining the Uluzzian through a new lithic assemblage from Roccia San Sebastiano (Mondragone, southern Italy)

abstract: Roccia San Sebastiano is a tectonic-karstic cave located at the foot of the southern slope of Mt. Massico, in the territory of Mondragone (Caserta) in Campania (southern Italy). Systematic excavation has been carried out since 2001, leading to the partial exploration of an important Pleistocene deposit, extraordinarily rich in lithic and faunal remains. The aim of this paper is to (1) present the stratigraphic sequence of Roccia San Sebastiano, and (2) technologically describe the lithic materials of squares F14 t18, t19, t20; E16 t16, t17, t18 recently recognised as Uluzzian. The stratigraphic sequence is more than 3 metres thick and dates from the Middle to the Upper Palaeolithic. It contains different techno-complexes: Gravettian, Aurignacian, Uluzzian and Mousterian. In the Uluzzian lithic assemblage mostly local pebbles of chert were used in order to produce small-sized objects. The concept of debitage mainly deals with unidirectional debitage with absent or fairly accurate management of the convexities and angles; the striking platforms are usually natural or made by one stroke. It is attested the use of both direct freehand percussion and bipolar technique on anvil in the same reduction sequence. Amongst the retouched tools the presence of two lunates is of note. This study of the Roccia San Sebastiano Uluzzian lithic complexes is significant for understanding the dynamics of the transition from Middle to Upper Palaeolithic in the Tyrrhenian margin of southern Italy.


title: A Game Theoretic Framework for Model Based Reinforcement Learning

abstract: Model-based reinforcement learning (MBRL) has recently gained immense interest due to its potential for sample efficiency and ability to incorporate off-policy data. However, designing stable and efficient MBRL algorithms using rich function approximators have remained challenging. To help expose the practical challenges in MBRL and simplify algorithm design from the lens of abstraction, we develop a new framework that casts MBRL as a game between: (1) a policy player, which attempts to maximize rewards under the learned model; (2) a model player, which attempts to fit the real-world data collected by the policy player. For algorithm development, we construct a Stackelberg game between the two players, and show that it can be solved with approximate bi-level optimization. This gives rise to two natural families of algorithms for MBRL based on which player is chosen as the leader in the Stackelberg game. Together, they encapsulate, unify, and generalize many previous MBRL algorithms. Furthermore, our framework is consistent with and provides a clear basis for heuristics known to be important in practice from prior works. Finally, through experiments we validate that our proposed algorithms are highly sample efficient, match the asymptotic performance of model-free policy gradient, and scale gracefully to high-dimensional tasks like dexterous hand manipulation. Additional details and code can be obtained from the project page at https://sites.google.com/view/mbrl-game


title: Quantifying uncertainties and correlations in the nuclear-matter equation of state

abstract: We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory ($\chi$EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the $\chi$EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. We then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium $\chi$EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.


title: MobiAxis: An Embodied Learning Task for Teaching Multiplication with a Social Robot

abstract: The use of robots in educational settings is growing increasingly popular. Yet, many of the learning tasks involving social robots do not take full advantage of their physical embodiment. MobiAxis is a proposed learning task which uses the physical capabilities of a Pepper robot to teach the concepts of positive and negative multiplication along a number line. The robot is embodied with a number of multi-modal socially intelligent features and behaviours which are designed to enhance learning. This paper is a position paper describing the technical and theoretical implementation of the task, as well as proposed directions for future studies.


title: Classification Benchmarks for Under-resourced Bengali Language based on Multichannel Convolutional-LSTM Network

abstract: Exponential growths of social media and micro-blogging sites not only provide platforms for empowering freedom of expressions and individual voices but also enables people to express anti-social behaviour like online harassment, cyberbullying, and hate speech. Numerous works have been proposed to utilize these data for social and anti-social behaviours analysis, document characterization, and sentiment analysis by predicting the contexts mostly for highly resourced languages such as English. However, there are languages that are under-resources, e.g., South Asian languages like Bengali, Tamil, Assamese, Telugu that lack of computational resources for the NLP tasks. In this paper, we provide several classification benchmarks for Bengali, an under-resourced language. We prepared three datasets of expressing hate, commonly used topics, and opinions for hate speech detection, document classification, and sentiment analysis, respectively. We built the largest Bengali word embedding models to date based on 250 million articles, which we call BengFastText. We perform three different experiments, covering document classification, sentiment analysis, and hate speech detection. We incorporate word embeddings into a Multichannel Convolutional-LSTM (MConv-LSTM) network for predicting different types of hate speech, document classification, and sentiment analysis. Experiments demonstrate that BengFastText can capture the semantics of words from respective contexts correctly. Evaluations against several baseline embedding models, e.g., Word2Vec and GloVe yield up to 92.30%, 82.25%, and 90.45% F1-scores in case of document classification, sentiment analysis, and hate speech detection, respectively during 5-fold cross-validation tests.


title: Mathematical analysis of the acoustic imaging modality using bubbles as contrast agents at nearly resonating frequencies

abstract: We analyze mathematically the acoustic imaging modality using bubbles as contrast agents. These bubbles are modeled by mass densities and bulk moduli enjoying contrasting scales. These contrasting scales allow them to resonate at certain incident frequencies. We consider two types of such contrasts. In the first one, the bubbles are light with small bulk modulus, as compared to the ones of the background, so that they generate the Minnaert resonance (corresponding to a local surface wave). In the second one, the bubbles have moderate mass density but still with small bulk modulus so that they generate a sequence of resonances (corresponding to local body waves). We propose to use as measurements the far-fields collected before and after injecting a bubble, set at a given location point in the target domain, generated at a band of incident frequencies and at a fixed single backscattering direction. Then, we scan the target domain with such bubbles and collect the corresponding far-fields. The goal is to reconstruct both the, variable, mass density and bulk modulus of the background in the target region. 1. We show that, for each fixed used bubble, the contrasted far-fields reach their maximum value at, incident, frequencies close to the Minneart resonance (or the body-wave resonances depending on the types of bubbles we use). Hence, we can reconstruct this resonance from our data. The explicit dependence of the Minnaert resonance in terms of the background mass density of the background allows us to recover it, i.e. the mass density, in a straightforward way. 2. In addition, this measured contrasted far-fields allow us to recover the total field at the location points of the bubbles (i.e. the total field in the absence of the bubbles). A numerical differentiation argument, for instance, allows us to recover the bulk modulus of the targeted region as well.


title: Security of quantum key distribution with detection-efficiency mismatch in the multiphoton case

abstract: Detection-efficiency mismatch is a common problem in practical quantum key distribution (QKD) systems. Current security proofs of QKD with detection-efficiency mismatch rely either on the assumption of the single-photon light source on the sender side or on the assumption of the single-photon input of the receiver side. These assumptions impose restrictions on the class of possible eavesdropping strategies. Here we present a rigorous security proof without these assumptions and, thus, solve this important problem and prove the security of QKD with detection-efficiency mismatch against general attacks (in the asymptotic regime). In particular, we adapt the decoy state method to the case of detection-efficiency mismatch.


title: Harmonic based model predictive control for set-point tracking

abstract: This paper presents a novel model predictive control (MPC) formulation for set-point tracking. Stabilizing predictive controllers based on terminal ingredients may exhibit stability and feasibility issues in the event of a reference change for small to moderate prediction horizons. In the MPC for tracking formulation, these issues are solved by the addition of an artificial equilibrium point as a new decision variable, providing a significantly enlarged domain of attraction and guaranteeing recursive feasibility for any reference change. However, it may suffer from performance issues if the prediction horizon is not large enough. This paper presents an extension of this formulation where a harmonic artificial reference is used in place of the equilibrium point. The proposed formulation achieves even greater domains of attraction and can significantly outperform other MPC formulations when the prediction horizon is small. We prove the asymptotic stability and recursive feasibility of the proposed controller, as well as provide guidelines for the design of its main ingredients. Finally, we highlight its advantages with a case study of a ball and plate system.


title: Minimising the impact of scale-dependent galaxy bias on the joint cosmological analysis of large scale structures

abstract: We present a mitigation strategy to reduce the impact of non-linear galaxy bias on the joint `$3 \times 2 $pt' cosmological analysis of weak lensing and galaxy surveys. The $\Psi$-statistics that we adopt are based on Complete Orthogonal Sets of E/B Integrals (COSEBIs). As such they are designed to minimise the contributions to the observable from the smallest physical scales where models are highly uncertain. We demonstrate that $\Psi$-statistics carry the same constraining power as the standard two-point galaxy clustering and galaxy-galaxy lensing statistics, but are significantly less sensitive to scale-dependent galaxy bias. Using two galaxy bias models, motivated by halo-model fits to data and simulations, we quantify the error in a standard $3 \times 2$pt analysis where constant galaxy bias is assumed. Even when adopting conservative angular scale cuts, that degrade the overall cosmological parameter constraints, we find of order $1 \sigma$ biases for Stage III surveys on the cosmological parameter $S_8 = \sigma_8(\Omega_{\rm m}/0.3)^{\alpha}$. This arises from a leakage of the smallest physical scales to all angular scales in the standard two-point correlation functions. In contrast, when analysing $\Psi$-statistics under the same approximation of constant galaxy bias, we show that the bias on the recovered value for $S_8$ can be decreased by a factor of $\sim 2$, with less conservative scale cuts. Given the challenges in determining accurate galaxy bias models in the highly non-linear regime, we argue that $3 \times 2$pt analyses should move towards new statistics that are less sensitive to the smallest physical scales.