title: An Extreme Precision Radial Velocity Pipeline: First Radial Velocities from EXPRES

abstract: The EXtreme PREcision Spectrograph (EXPRES) is an environmentally stabilized, fiber-fed, $R=137,500$, optical spectrograph. It was recently commissioned at the 4.3-m Lowell Discovery Telescope (LDT) near Flagstaff, Arizona. The spectrograph was designed with a target radial-velocity (RV) precision of 30$\mathrm{~cm~s^{-1}}$. In addition to instrumental innovations, the EXPRES pipeline, presented here, is the first for an on-sky, optical, fiber-fed spectrograph to employ many novel techniques---including an "extended flat" fiber used for wavelength-dependent quantum efficiency characterization of the CCD, a flat-relative optimal extraction algorithm, chromatic barycentric corrections, chromatic calibration offsets, and an ultra-precise laser frequency comb for wavelength calibration. We describe the reduction, calibration, and radial-velocity analysis pipeline used for EXPRES and present an example of our current sub-meter-per-second RV measurement precision, which reaches a formal, single-measurement error of 0.3$\mathrm{~m~s^{-1}}$ for an observation with a per-pixel signal-to-noise ratio of 250. These velocities yield an orbital solution on the known exoplanet host 51 Peg that matches literature values with a residual RMS of 0.895$\mathrm{~m~s^{-1}}$.


title: Performance Verification of the EXtreme PREcision Spectrograph

abstract: The EXtreme PREcision Spectrograph (EXPRES) is a new Doppler spectrograph designed to reach a radial velocity measurement precision sufficient to detect Earth-like exoplanets orbiting nearby, bright stars. We report on extensive laboratory testing and on-sky observations to quantitatively assess the instrumental radial velocity measurement precision of EXPRES, with a focused discussion of individual terms in the instrument error budget. We find that EXPRES can reach a single-measurement instrument calibration precision better than 10 cm/s, not including photon noise from stellar observations. We also report on the performance of the various environmental, mechanical, and optical subsystems of EXPRES, assessing any contributions to radial velocity error. For atmospheric and telescope related effects, this includes the fast tip-tilt guiding system, atmospheric dispersion compensation, and the chromatic exposure meter. For instrument calibration, this includes the laser frequency comb (LFC), flat-field light source, CCD detector, and effects in the optical fibers. Modal noise is mitigated to a negligible level via a chaotic fiber agitator, which is especially important for wavelength calibration with the LFC. Regarding detector effects, we empirically assess the impact on radial velocity precision due to pixel-position non-uniformities (PPNU) and charge transfer inefficiency (CTI). EXPRES has begun its science survey to discover exoplanets orbiting G-dwarf and K-dwarf stars, in addition to transit spectroscopy and measurements of the Rossiter-McLaughlin effect.


title: Semi-closed form solutions for barrier and American options written on a time-dependent Ornstein Uhlenbeck process

abstract: In this paper we develop a semi-closed form solutions for the barrier (perhaps, time-dependent) and American options written on the underlying stock which follows a time-dependent OU process with a log-normal drift. This model is equivalent to the familiar Hull-White model in FI, or a time dependent OU model in FX. Semi-closed form means that given the time-dependent interest rate, continuous dividend and volatility functions, one need to solve numerically a linear (for the barrier option) or nonlinear (for the American option) Fredholm equation of the first kind. After that the option prices in all cases are presented as one-dimensional integrals of combination of the above solutions and Jacobi theta functions. We also demonstrate that computationally our method is more efficient than the backward finite difference method used for solving these problems, and can also be as efficient as the forward finite difference solver while providing better accuracy and stability.


title: Goal-Conditioned End-to-End Visuomotor Control for Versatile Skill Primitives

abstract: Visuomotor control (VMC) is an effective means of achieving basic manipulation tasks such as pushing or pick-and-place from raw images. Conditioning VMC on desired goal states is a promising way of achieving versatile skill primitives. However, common conditioning schemes either rely on task-specific fine tuning - e.g. using one-shot imitation learning (IL) - or on sampling approaches using a forward model of scene dynamics i.e. model-predictive control (MPC), leaving deployability and planning horizon severely limited. In this paper we propose a conditioning scheme which avoids these pitfalls by learning the controller and its conditioning in an end-to-end manner. Our model predicts complex action sequences based directly on a dynamic image representation of the robot motion and the distance to a given target observation. In contrast to related works, this enables our approach to efficiently perform complex manipulation tasks from raw image observations without predefined control primitives or test time demonstrations. We report significant improvements in task success over representative MPC and IL baselines. We also demonstrate our model's generalisation capabilities in challenging, unseen tasks featuring visual noise, cluttered scenes and unseen object geometries.


title: Integrated Power and Thermal Management of Connected HEVs via Multi-Horizon MPC

abstract: In this paper, a multi-horizon model predictive controller (MH-MPC) is developed for integrated power and thermal management (iPTM) of a power-split hybrid electric vehicle (HEV). The proposed MH-MPC leverages an accurate short-horizon vehicle speed preview and an approximate forecast over a longer shrinking horizon till the end of the driving cycle. This multiple-horizon scheme is developed to cope with fast and slow dynamics associated with power and thermal responses. The main objective of the proposed MH-MPC is to minimize fuel consumption and enforce the power and thermal constraints on the battery state-of-charge and engine coolant temperature, while meeting the driving (traction) and cabin air conditioning (heating) demands. The proposed MH-MPC allows for exploiting the engine coolant as thermal energy storage, providing more flexibility for the HEV energy flow optimization. The simulation results show that the proposed MH-MPC provides near-optimal results in reference to the Dynamic Programming (DP) solution with an affordable computational cost. Moreover, compared with a more conventional MPC strategy, the MH-MPC can leverage the speed previews with different resolutions effectively to achieve the desired performance with satisfactory robustness.


title: Distance evolutions in growing preferential attachment graphs

abstract: We study the evolution of the graph distance and weighted distance between two fixed vertices in dynamically growing random graph models. More precisely, we consider preferential attachment models with power-law exponent $\tau\in(2,3)$, sample two vertices $u_t,v_t$ uniformly at random when the graph has $t$ vertices, and study the evolution of the graph distance between these two fixed vertices as the surrounding graph grows. This yields a discrete-time stochastic process in $t'\geq t$, called the distance evolution. We show that there is a tight strip around the function $4\frac{\log\log(t)-\log(\log(t'/t)\vee1)}{|\log(\tau-2)|}\vee 2$ that the distance evolution never leaves with high probability as $t$ tends to infinity. We extend our results to weighted distances, where every edge is equipped with an i.i.d. copy of a non-negative random variable $L$.


title: Examples of property (T) II$_1$ factors with trivial fundamental group

abstract: In this article we provide the first examples of property (T) $\rm II_1$ factors $\mathcal N$ with trivial fundamental group, $\mathcal F (\mathcal N)=1$. Our examples arise as group factors $\mathcal N=\mathcal L(G)$ where $G$ belong to two distinct families of property (T) groups previously studied in the literature: the groups introduced by Valette in \cite{Va04} and the ones introduced recently in \cite{CDK19} using the Belegradek-Osin Rips construction from \cite{BO06}. In particular, our results provide a continuum of explicit pairwise non-isomorphic property (T) factors.


title: Nonparametric estimation of variable productivity Hawkes processes

abstract: An extension of the Hawkes model where the productivity is variable is considered. In particular, the case is considered where each point may have its own productivity and a simple analytic formula is derived for the maximum likelihood estimators of these productivities. This estimator is compared with an empirical estimator and ways are explored of stabilizing both estimators by lower truncating, smoothing, and rescaling the estimates. Properties of the estimators are explored in simulations, and the methods are applied to seismological and epidemic datasets to show and quantify substantial variation in productivity.


title: New, multipole solutions of relativistic, viscous hydrodynamics

abstract: We present a new class of exact fireball solutions of relativistic dissipative hydrodynamics. We describe new exact solutions both for the relativistic Navier-Stokes and for the Israel-Stewart theory, for arbitrary shear and bulk viscosities, as well as for other dissipative coefficients. The common property of these solutions is the presence of the relativistic Hubble flow. Our results generalize the recently found first solution in these classes, for an arbitrary temperature dependent speed of sound, shear and bulk viscosity, heat conduction and fluctuating initial temperature profiles. These solutions are causal and not only stable but also asymptotically perfect. A strong and narrow peak in the kinematic bulk viscosity is shown to imitate the effects of a first order phase transition. The new class of asymptotically perfect solutions is thus found to be very rich, but at the same time mostly academic as the solutions are limited by the spherical symmetry of the Hubble flow field.


title: Adaptive motion control of parallel robots with kinematic and dynamic uncertainties

abstract: One of the most challenging issues in adaptive control of robot manipulators with kinematic uncertainties is requirement of the inverse of Jacobian matrix in regressor form. This requirement is inevitable in the case of the control of parallel robots, whose dynamic equations are written directly in the task space. In this paper, an adaptive controller is designed for parallel robots based on representation of Jacobian matrix in regressor form, such that asymptotic trajectory tracking is ensured. The main idea is separation of determinant and adjugate of Jacobian matrix and then organize new regressor forms. Simulation and experimental results on a 2--DOF R\underline{P}R and 3--DOF redundant cable driven robot, verify promising performance of the proposed methods.


title: Face-Off: Adversarial Face Obfuscation

abstract: Advances in deep learning have made face recognition technologies pervasive. While useful to social media platforms and users, this technology carries significant privacy threats. Coupled with the abundant information they have about users, service providers can associate users with social interactions, visited places, activities, and preferences--some of which the user may not want to share. Additionally, facial recognition models used by various agencies are trained by data scraped from social media platforms. Existing approaches to mitigate these privacy risks from unwanted face recognition result in an imbalanced privacy-utility trade-off to users. In this paper, we address this trade-off by proposing Face-Off, a privacy-preserving framework that introduces strategic perturbations to the user's face to prevent it from being correctly recognized. To realize Face-Off, we overcome a set of challenges related to the black-box nature of commercial face recognition services, and the scarcity of literature for adversarial attacks on metric networks. We implement and evaluate Face-Off to find that it deceives three commercial face recognition services from Microsoft, Amazon, and Face++. Our user study with 423 participants further shows that the perturbations come at an acceptable cost for the users.


title: The focusing NLS equation with step-like oscillating background: scenarios of long-time asymptotics

abstract: We consider the Cauchy problem for the focusing nonlinear Schr\"odinger equation with initial data approaching two different plane waves $A_j\mathrm{e}^{\mathrm{i}\phi_j}\mathrm{e}^{-2\mathrm{i}B_jx}$, $j=1,2$ as $x\to\pm\infty$. Using Riemann-Hilbert techniques and Deift-Zhou steepest descent arguments, we study the long-time asymptotics of the solution. We detect that each of the cases $B_1<B_2$, $B_1>B_2$, and $B_1=B_2$ deserves a separate analysis. Focusing mainly on the first case, the so-called shock case, we show that there is a wide range of possible asymptotic scenarios. We also propose a method for rigorously establishing the existence of certain higher-genus asymptotic sectors.


title: Towards a Computer Vision Particle Flow

abstract: In High Energy Physics experiments Particle Flow (PFlow) algorithms are designed to provide an optimal reconstruction of the nature and kinematic properties of the particles produced within the detector acceptance during collisions. At the heart of PFlow algorithms is the ability to distinguish the calorimeter energy deposits of neutral particles from those of charged particles, using the complementary measurements of charged particle tracking devices, to provide a superior measurement of the particle content and kinematics. In this paper, a computer vision approach to this fundamental aspect of PFlow algorithms, based on calorimeter images, is proposed. A comparative study of the state of the art deep learning techniques is performed. A significantly improved reconstruction of the neutral particle calorimeter energy deposits is obtained in a context of large overlaps with the deposits from charged particles. Calorimeter images with augmented finer granularity are also obtained using super-resolution techniques.

categories: physics.data-an hep-ex hep-ph physics.ins-det stat.ML


title: Interaction induced dynamical $\mathcal{PT}$ symmetry breaking in dissipative Fermi-Hubbard models

abstract: We investigate the dynamical properties of one-dimensional dissipative Fermi-Hubbard models, which are described by the Lindblad master equations with site-dependent jump operators. The corresponding non-Hermitian effective Hamiltonians with pure loss terms possess parity-time ($\mathcal{PT}$) symmetry if we compensate the system additionally an overall gain term. By solving the two-site Lindblad equation with fixed dissipation exactly, we find that the dynamics of rescaled density matrix shows an instability as the interaction increases over a threshold, which can be equivalently described in the scheme of non-Hermitian effective Hamiltonians. This instability is also observed in multi-site systems and closely related to the $\mathcal{PT}$ symmetry breaking accompanied by appearance of complex eigenvalues of the effective Hamiltonian. Moreover, we unveil that the dynamical instability of the anti-ferromagnetic Mott phase comes from the $\mathcal{PT}$ symmetry breaking in highly excited bands, although the low-energy effective model of the non-Hermitian Hubbard model in the strongly interacting regime is always Hermitian. We also provide a quantitative estimation of the time for the observation of dynamical $\mathcal{PT}$ symmetry breaking which could be probed in experiments.


title: DRST: Deep Residual Shearlet Transform for Densely Sampled Light Field Reconstruction

abstract: The Image-Based Rendering (IBR) approach using Shearlet Transform (ST) is one of the most effective methods for Densely-Sampled Light Field (DSLF) reconstruction. The ST-based DSLF reconstruction typically relies on an iterative thresholding algorithm for Epipolar-Plane Image (EPI) sparse regularization in shearlet domain, involving dozens of transformations between image domain and shearlet domain, which are in general time-consuming. To overcome this limitation, a novel learning-based ST approach, referred to as Deep Residual Shearlet Transform (DRST), is proposed in this paper. Specifically, for an input sparsely-sampled EPI, DRST employs a deep fully Convolutional Neural Network (CNN) to predict the residuals of the shearlet coefficients in shearlet domain in order to reconstruct a densely-sampled EPI in image domain. The DRST network is trained on synthetic Sparsely-Sampled Light Field (SSLF) data only by leveraging elaborately-designed masks. Experimental results on three challenging real-world light field evaluation datasets with varying moderate disparity ranges (8 - 16 pixels) demonstrate the superiority of the proposed learning-based DRST approach over the non-learning-based ST method for DSLF reconstruction. Moreover, DRST provides a 2.4x speedup over ST, at least.


title: Spatially Adaptive Inference with Stochastic Feature Sampling and Interpolation

abstract: In the feature maps of CNNs, there commonly exists considerable spatial redundancy that leads to much repetitive processing. Towards reducing this superfluous computation, we propose to compute features only at sparsely sampled locations, which are probabilistically chosen according to activation responses, and then densely reconstruct the feature map with an efficient interpolation procedure. With this sampling-interpolation scheme, our network avoids expending computation on spatial locations that can be effectively interpolated, while being robust to activation prediction errors through broadly distributed sampling. A technical challenge of this sampling-based approach is that the binary decision variables for representing discrete sampling locations are non-differentiable, making them incompatible with backpropagation. To circumvent this issue, we make use of a reparameterization trick based on the Gumbel-Softmax distribution, with which backpropagation can iterate these variables towards binary values. The presented network is experimentally shown to save substantial computation while maintaining accuracy over a variety of computer vision tasks.


title: Analysis of a fully discrete approximation for the classical Keller--Segel model: lower and a priori bounds

abstract: This paper is devoted to constructing approximate solutions for the classical Keller--Segel model governing \emph{chemotaxis}. It consists of a system of nonlinear parabolic equations, where the unknowns are the average density of cells (or organisms), which is a conserved variable, and the average density of chemoattractant. The numerical proposal is made up of a crude finite element method together with a mass lumping technique and a semi-implicit Euler time integration. The resulting scheme turns out to be linear and decouples the computation of variables. The approximate solutions keep lower bounds -- positivity for the cell density and nonnegativity for the chemoattractant density --, are bounded in the $L^1(\Omega)$-norm, satisfy a discrete energy law, and have \emph{ a priori} energy estimates. The latter is achieved by means of a discrete Moser--Trudinger inequality. As far as we know, our numerical method is the first one that can be encountered in the literature dealing with all of the previously mentioned properties at the same time. Furthermore, some numerical examples are carried out to support and complement the theoretical results.


title: Ostwald growth rate in controlled Covid-19 epidemic spreading as in arrested growth in quantum complex matter

abstract: Here, we focus on the data analysis of the growth of epidemic spread of Covid-19 in countries where different policies of containment were activated. It is known that the growth of pandemic spread at its threshold is exponential, but it is not known how to quantify the success of different containment policies. We identify that a successful approach gives an arrested phase regime following the Ostwald growth, where, over the course of time, one phase transforms into another metastable phase with a similar free energy as observed in oxygen interstitial diffusion in quantum complex matter and in crystallization of proteins. We introduce the s factor which provides a quantitative measure of the efficiency and speed of the adopted containment policy, which is very helpful not only to monitor the Covid-19 pandemic spread but also for other countries to choose the best containment policy. The results show that a policy based on joint confinement, targeted tests, and tracking positive cases is the most rapid pandemic containment policy; in fact, we found values of 9, 5, and 31 for the success s factor for China, South Korea, and Italy, respectively, where the lowest s factor indicates the best containment policy


title: Sensitive capacitive pressure sensors based on graphene membrane arrays

abstract: The high flexibility, impermeability and strength of graphene membranes are key properties that can enable the next generation of nanomechanical sensors. However, for capacitive pressure sensors the sensitivity offered by a single suspended graphene membrane is too small to compete with commercial sensors. Here, we realize highly sensitive capacitive pressure sensors consisting of arrays of nearly ten thousand small, freestanding double-layer graphene membranes. We fabricate large arrays of small diameter membranes using a procedure that maintains the superior material and mechanical properties of graphene, even after high-temperature anneals. These sensors are readout using a low cost battery-powered circuit board, with a responsivity of up to 47.8 aF Pa$^{-1}$ mm$^{-2}$, thereby outperforming commercial sensors.


title: Brain tumor segmentation with missing modalities via latent multi-source correlation representation

abstract: Multimodal MR images can provide complementary information for accurate brain tumor segmentation. However, it's common to have missing imaging modalities in clinical practice. Since there exists a strong correlation between multi modalities, a novel correlation representation block is proposed to specially discover the latent multi-source correlation. Thanks to the obtained correlation representation, the segmentation becomes more robust in the case of missing modalities. The model parameter estimation module first maps the individual representation produced by each encoder to obtain independent parameters, then, under these parameters, the correlation expression module transforms all the individual representations to form a latent multi-source correlation representation. Finally, the correlation representations across modalities are fused via the attention mechanism into a shared representation to emphasize the most important features for segmentation. We evaluate our model on BraTS 2018 datasets, it outperforms the current state-of-the-art method and produces robust results when one or more modalities are missing.


title: Impact of spectral effects on photovoltaic energy production: A case study in the United States

abstract: The time averaged efficiency of photovoltaic modules in the field is generally lower than the efficiency measured in the laboratory under standard testing conditions due to the combined effects of temperature and spectral variability, affecting the bankability of power plant projects. We report correction factors to account for spectral effects ranging from -2% to 1.3% of the produced energy for silicon modules depending on location and collector geometry. In high irradiance locations, the energy yield advantage of trackers is underestimated by 0.4% if spectral sensitivity effects are neglected. We find a correlation between the locations most favourable for tracking, and those most favourable for multijunctions. As the photovoltaic market grows to a multi-terawatt size, these seemingly small effects are expected to have an economic impact equivalent to tens of billions of dollars in the next few decades, far out-weighting the cost of the required research effort.


title: Across Scales & Across Dimensions: Temporal Super-Resolution using Deep Internal Learning

abstract: When a very fast dynamic event is recorded with a low-framerate camera, the resulting video suffers from severe motion blur (due to exposure time) and motion aliasing (due to low sampling rate in time). True Temporal Super-Resolution (TSR) is more than just Temporal-Interpolation (increasing framerate). It can also recover new high temporal frequencies beyond the temporal Nyquist limit of the input video, thus resolving both motion-blur and motion-aliasing effects that temporal frame interpolation (as sophisticated as it maybe) cannot undo. In this paper we propose a "Deep Internal Learning" approach for true TSR. We train a video-specific CNN on examples extracted directly from the low-framerate input video. Our method exploits the strong recurrence of small space-time patches inside a single video sequence, both within and across different spatio-temporal scales of the video. We further observe (for the first time) that small space-time patches recur also across-dimensions of the video sequence - i.e., by swapping the spatial and temporal dimensions. In particular, the higher spatial resolution of video frames provides strong examples as to how to increase the temporal resolution of that video. Such internal video-specific examples give rise to strong self-supervision, requiring no data but the input video itself. This results in Zero-Shot Temporal-SR of complex videos, which removes both motion blur and motion aliasing, outperforming previous supervised methods trained on external video datasets.


title: Adsorption of H$_2$ on Amorphous Solid Water Studied with Molecular Dynamics Simulations

abstract: We investigated the behavior of H$_2$, main constituent of the gas phase in dense clouds, after collision with amorphous solid water (ASW) surfaces, one of the most abundant chemical species of interstellar ices. We developed a general framework to study the adsorption dynamics of light species on interstellar ices. We provide binding energies and their distribution, sticking probabilities for incident energies between 1 meV and 60 meV, and thermal sticking coefficients between 10 and 300 K for surface temperatures from 10 to 110 K. We found that the sticking probability depends strongly on the adsorbate kinetic energy and the surface temperature, but hardly on the angle of incidence. We observed finite sticking probabilities above the thermal desorption temperature. Adsorption and thermal desorption should be considered as separate events with separate time scales. Laboratory results for these species have shown a gap in the trends attributed to the differently employed experimental techniques. Our results complement observations and extend them, increasing the range of gas temperatures under consideration. We plan to employ our method to study a variety of adsorbates, including radical and charged species.


title: On the Detectability of Conflict: a Remote Sensing Study of the Rohingya Conflict

abstract: The detection and quantification of conflict through remote sensing modalities represents a challenging but crucial aspect of human rights monitoring. In this work we demonstrate how utilizing multi-modal data sources can help build a comprehensive picture of conflict and human displacement, using the Rohingya conflict in the state of Rakhine, Myanmar as a case study. We show that time series analysis of fire detections from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Visible Infrared Imaging Radiometer Suite (VIIRS) can reveal anomalous spatial and temporal distributions of fires related to conflict. This work also shows that Synthetic Aperture Radar (SAR) backscatter and coherence data can detect the razing and burning of buildings and villages, even in cloudy conditions. These techniques may be further developed in the future to enable the monitoring and detection of signals originating from these types of conflict.


title: Beheshti-NER: Persian Named Entity Recognition Using BERT

abstract: Named entity recognition is a natural language processing task to recognize and extract spans of text associated with named entities and classify them in semantic Categories. Google BERT is a deep bidirectional language model, pre-trained on large corpora that can be fine-tuned to solve many NLP tasks such as question answering, named entity recognition, part of speech tagging and etc. In this paper, we use the pre-trained deep bidirectional network, BERT, to make a model for named entity recognition in Persian. We also compare the results of our model with the previous state of the art results achieved on Persian NER. Our evaluation metric is CONLL 2003 score in two levels of word and phrase. This model achieved second place in NSURL-2019 task 7 competition which associated with NER for the Persian language. our results in this competition are 83.5 and 88.4 f1 CONLL score respectively in phrase and word level evaluation.


title: Learning to Fly via Deep Model-Based Reinforcement Learning

abstract: Learning to control robots without requiring engineered models has been a long-term goal, promising diverse and novel applications. Yet, reinforcement learning has only achieved limited impact on real-time robot control due to its high demand of real-world interactions. In this work, by leveraging a learnt probabilistic model of drone dynamics, we learn a thrust-attitude controller for a quadrotor through model-based reinforcement learning. No prior knowledge of the flight dynamics is assumed; instead, a sequential latent variable model, used generatively and as an online filter, is learnt from raw sensory input. The controller and value function are optimised entirely by propagating stochastic analytic gradients through generated latent trajectories. We show that "learning to fly" can be achieved with less than 30 minutes of experience with a single drone, and can be deployed solely using onboard computational resources and sensors, on a self-built drone.


title: Abstraction, Up-to Techniques and Games for Systems of Fixpoint Equations

abstract: Systems of fixpoint equations over complete lattices, consisting of (mixed) least and greatest fixpoint equations, allow one to express a number of verification tasks such as model-checking of various kinds of specification logics or the check of coinductive behavioural equivalences. In this paper we develop a theory of approximation for systems of fixpoint equations in the style of abstract interpretation: a system over some concrete domain is abstracted to a system in a suitable abstract domain, with conditions ensuring that the abstract solution represents a sound/complete overapproximation of the concrete solution. Interestingly, up-to techniques, a classical approach used in coinductive settings to obtain easier or feasible proofs, can be interpreted as abstractions in a way that they naturally fit in our framework and extend to systems of equations. Additionally, relying on the approximation theory, we can provide a characterisation of the solution of systems of fixpoint equations over complete lattices in terms of a suitable parity game, generalising some recent work that was restricted to continuous lattices. The game view opens the way to the development of on-the-fly algorithms for characterising the solution of such equation systems.


title: Uniform Annihilators of Systems of parameters

abstract: In this paper we give uniform annihilators for some relations of all systems of parameters in a local ring. Our results provide simple proofs for some classical results.


title: Revisiting the Common Neighbour Analysis and the Centrosymmetry Parameter

abstract: We review two standard methods for structural classification in simulations of crystalline phases, the Common Neighbour Analysis and the Centrosymmetry Parameter. We explore the definitions and implementations of each of their common variants, and investigate their respective failure modes and classification biases. Simple modifications to both methods are proposed, which improve their robustness, interpretability, and applicability. We denote these variants the Interval Common Neighbour Analysis, and the Minimum-Weight Matching Centrosymmetry Parameter.


title: Augmented Lagrangian based first-order methods for convex-constrained programs with weakly-convex objective

abstract: First-order methods (FOMs) have been widely used for solving large-scale problems. A majority of existing works focus on problems without constraint or with simple constraints. Several recent works have studied FOMs for problems with complicated functional constraints. In this paper, we design a novel augmented Lagrangian (AL) based FOM for solving problems with non-convex objective and convex constraint functions. The new method follows the framework of the proximal point (PP) method. On approximately solving PP subproblems, it mixes the usage of the inexact AL method (iALM) and the quadratic penalty method, while the latter is always fed with estimated multipliers by the iALM. We show a complexity result of $O(\varepsilon^{-\frac{5}{2}}|\log\varepsilon|)$ for the proposed method to achieve an $\varepsilon$-KKT point. This is the best known result. Theoretically, the hybrid method has lower iteration-complexity requirement than its counterpart that only uses iALM to solve PP subproblems, and numerically, it can perform significantly better than a pure-penalty-based method. Numerical experiments are conducted on nonconvex linearly constrained quadratic programs and nonconvex QCQP. The numerical results demonstrate the efficiency of the proposed methods over existing ones.


title: Characterizing multipartite entanglement by violation of CHSH inequalities

abstract: Entanglement of high-dimensional and multipartite quantum systems offer promising perspectives in quantum information processing. However, the characterization and measure of such kind of entanglement is of great challenge. Here we consider the overlaps between the maximal quantum mean values and the classical bound of the CHSH inequalities for pairwise-qubit states in two-dimensional subspaces. We show that the concurrence of a pure state in any high-dimensional multipartite system can be equivalently represented by these overlaps. Here we consider the projections of an arbitrary high-dimensional multipartite state to two-qubit states. We investigate the non-localities of these projected two-qubit sub-states by their violations of CHSH inequalities. From these violations, the overlaps between the maximal quantum mean values and the classical bound of the CHSH inequality, we show that the concurrence of a high-dimensional multipartite pure state can be exactly expressed by these overlaps. We further derive a lower bound of the concurrence for any quantum states, which is tight for pure states. The lower bound not only imposes restriction on the non-locality distributions among the pairwise qubit states, but also supplies a sufficient condition for distillation of bipartite entanglement. Effective criteria for detecting genuine tripartite entanglement and the lower bound of concurrence for genuine tripartite entanglement are also presented based on such non-localities.


title: MENA Compared to Europe: The Influence of Land Use, Nuclear Power, and Transmission Expansion on Renewable Electricity System Costs

abstract: Most studies that examine CO2-neutral, or near CO2-neutral, power systems by using energy system models investigate Europe or the United States, while similar studies for other regions are rare. In this paper, we focus on the Middle East and North Africa (MENA), where weather conditions differ substantially from those in Europe. We use a green-field linear capacity expansion model with over-night investment to assess the effect on the system cost of (i) limiting/expanding the amount of land available for wind and solar farms, (ii) allowing for nuclear power and (iii) disallowing for international transmission. The assessment is done under three different cost regimes for solar PV and battery storage. First, we find that the amount of available land for wind and solar farms can have a significant impact on the system cost, with a cost increase of 0-50% as a result of reduced available land. Second, allowing for nuclear power has a minor effect in MENA, while it may decrease the system cost in Europe by up to 20%. In Europe, the effect on system cost from allowing for nuclear power is highly dependent on the PV and battery cost regime. Third, disallowing for international transmission increases the system cost by up to 25% in both Europe and MENA. The impacts on system cost from these three controversial and policy-relevant factors in a decarbonized power system thus play out differently, depending on (i) the region and (ii) uncertain future investment costs for solar PV and storage. We conclude that a renewable power system in MENA is likely to be less costly than one in Europe, irrespective of future uncertainties regarding investment cost for PV and batteries, and policies surrounding nuclear power, transmission, and land available for wind and solar farms.


title: Interaction-round-a-face and consistency-around-a-face-centered-cube

abstract: There is a correspondence between integrable lattice models of statistical mechanics and discrete integrable equations which satisfy multidimensional consistency, where the latter may be found in a quasi-classical expansion of the former. This paper extends this correspondence to interaction-round-a-face (IRF) models, resulting in a new formulation of the consistency-around-a-cube (CAC) integrability condition applicable to five-point equations in the square lattice. Multidimensional consistency for these equations is formulated as consistency-around-a-face-centered-cube (CAFCC), which namely involves satisfying an overdetermined system of fourteen five-point lattice equations for eight unknown variables on the face-centered cubic unit cell. From the quasi-classical limit of IRF models, which are constructed from the continuous spin solutions of the star-triangle relations associated to the Adler-Bobenko-Suris (ABS) list, fifteen sets of equations are obtained which satisfy CAFCC.


title: Geometrically finite transcendental entire functions

abstract: For polynomials, local connectivity of Julia sets is a much-studied and important property. Indeed, when the Julia set of a polynomial of degree $d\geq 2$ is locally connected, the topological dynamics can be completely described as a quotient of a much simpler system: angle $d$-tupling on the circle. For a transcendental entire function, local connectivity is less significant, but we may still ask for a description of the topological dynamics as the quotient of a simpler system. To this end, we introduce the notion of "docile" functions: a transcendental entire function with bounded postsingular set is docile if it is the quotient of a suitable disjoint-type function. Moreover, we prove docility for the large class of geometrically finite transcendental entire functions with bounded criticality on the Julia set. This can be seen as an analogue of the local connectivity of Julia sets for geometrically finite polynomials, first proved by Douady and Hubbard, and extends previous work of the second author and of Mihaljevi\'c for more restrictive classes of entire functions.


title: Unique Geometry and Texture from Corresponding Image Patches

abstract: We present a sufficient condition for recovering unique texture and viewpoints from unknown orthographic projections of a flat texture process. We show that four observations are sufficient in general, and we characterize the ambiguous cases. The results are applicable to shape from texture and texture-based structure from motion.


title: Native point defects and low $p$-doping efficiency in $Mg_2 (Si,Sn)$ solid solutions: A hybrid-density functional study

abstract: We perform hybrid-density functional calculations to investigate the charged defect formation energy of native point defects in $Mg_2 Si$, $Mg_2 Sn$, and their solid solutions. The band gap correction by hybrid-density functional is found to be critical to determine the charged defect density in these materials. For $Mg_2 Si$, $Mg$ interstitials are dominant and provide unintentional $n$-type conductivity. Additionally, as the $Mg$ vacancies can dominate in $Mg$-poor $Mg_2 Sn$, $p$-type conductivity is possible for $Mg_2 Sn$. However, the existence of low formation energy defects such as $Mg_{Sn}^{1+}$ and $I_{Mg}^{2+}$ in $Mg_2 Sn$ and their diffusion can cause severe charge compensation of hole carriers resulting in low $p$-type doping efficiency and thermal degradation. Our results indicate that, in addition to the extrinsic doping strategy, alloying of $Mg_2 Si$ with $Mg_2 Sn$ under $Mg$-poor conditions would be necessary to enhance the $p$-type conductivity with less charge compensation.


title: Dependence of nonthermal metallization kinetics on bond ionicity of compounds

abstract: It is known that covalently bonded materials undergo nonthermal structure transformations upon ultrafast excitation of an electronic system, whereas metals exhibit phonon hardening. Here we study how ionic bonds react to electronic excitation. Density-functional molecular dynamics predicts that ionic crystals may melt nonthermally, however, into an electronically insulating state, in contrast to covalent materials. We demonstrate that the band gap behavior during nonthermal transitions depends on a bonding type: it is harder to collapse the band gap in more ionic compounds, which is illustrated by transformations in Y2O3 vs. NaCl, LiF and KBr.


title: A numerical approach toward the $p$-adic Beilinson conjecture for elliptic curves over $\mathbb{Q}$

abstract: Restricting ourselves to elliptic curves over $\mathbb{Q}$, we reformulate the $p$-adic Beilinson conjecture due to Perrin-Riou, which is customized to our computational approach. We then develop a new algorithm for numerical verifications of the $p$-adic Beilinson conjecture, which is based on the theory of rigid cohomology and $F$-isocrystals.


title: Infinite-pressure phase diagram of binary mixtures of (non)additive hard disks

abstract: One versatile route to the creation of two-dimensional crystal structures on the nanometer to micrometer scale is the self-assembly of colloidal particles at an interface. Here, we explore the crystal phases that can be expected from the self-assembly of mixtures of spherical particles of two different sizes, which we map to (additive or non-additive) hard-disk mixtures. We map out the infinite-pressure phase diagram for these mixtures, using Floppy Box Monte Carlo simulations to systematically sample candidate crystal structures with up to 12 disks in the unit cell. As a function of the size ratio and number ratio of the two species of particles, we find a rich variety of periodic crystal structures. Additionally, we identify random tiling regions to predict random tiling quasicrystal stability ranges. Increasing non-additivity both gives rise to additional crystal phases and broadens the stability regime for crystal structures involving a large number of large-small contacts, including random tilings. Our results provide useful guidelines for controlling the self-assembly of colloidal particles at interfaces.


title: Local Rotation Invariance in 3D CNNs

abstract: Locally Rotation Invariant (LRI) image analysis was shown to be fundamental in many applications and in particular in medical imaging where local structures of tissues occur at arbitrary rotations. LRI constituted the cornerstone of several breakthroughs in texture analysis, including Local Binary Patterns (LBP), Maximum Response 8 (MR8) and steerable filterbanks. Whereas globally rotation invariant Convolutional Neural Networks (CNN) were recently proposed, LRI was very little investigated in the context of deep learning. LRI designs allow learning filters accounting for all orientations, which enables a drastic reduction of trainable parameters and training data when compared to standard 3D CNNs. In this paper, we propose and compare several methods to obtain LRI CNNs with directional sensitivity. Two methods use orientation channels (responses to rotated kernels), either by explicitly rotating the kernels or using steerable filters. These orientation channels constitute a locally rotation equivariant representation of the data. Local pooling across orientations yields LRI image analysis. Steerable filters are used to achieve a fine and efficient sampling of 3D rotations as well as a reduction of trainable parameters and operations, thanks to a parametric representations involving solid Spherical Harmonics (SH), which are products of SH with associated learned radial profiles.Finally, we investigate a third strategy to obtain LRI based on rotational invariants calculated from responses to a learned set of solid SHs. The proposed methods are evaluated and compared to standard CNNs on 3D datasets including synthetic textured volumes composed of rotated patterns, and pulmonary nodule classification in CT. The results show the importance of LRI image analysis while resulting in a drastic reduction of trainable parameters, outperforming standard 3D CNNs trained with data augmentation.


title: Trapped Rydberg ions: a new platform for quantum information processing

abstract: In this chapter, we present an overview of experiments with trapped Rydberg ions and outline the advantages and challenges of developing applications of this new platform for quantum computing, sensing and simulation. Trapped Rydberg ions feature several important properties, unique in their combination: they are tightly bound in a harmonic potential of a Paul trap, in which their internal and external degrees of freedom can be controlled in a precise fashion. High fidelity state preparation of both internal and motional states of the ions has been demonstrated, and the internal states have been employed to store and manipulate qubit information. Furthermore, strong dipolar interactions can be realised between ions in Rydberg states and be explored for investigating correlated many-body systems. By laser coupling to Rydberg states, the polarisability of the ions can be both enhanced and tuned. This can be used to control the interactions with the trapping fields in a Paul trap as well as dipolar interactions between the ions. Thus, trapped Rydberg ions present an attractive alternative for fast entangling operations as compared to those mediated by normal modes of trapped ions, which are advantageous for a future quantum computer or a quantum simulator.


title: Phase transition dynamics and gravitational wave spectra of strong first-order phase transition in supercooled universe

abstract: Phase transition dynamics may play important roles in the evolution history of the early universe, such as its possible roles in electroweak baryogenesis and dark matter.We systematically discuss and clarify the important details of the phase transition dynamics during a strong first-order phase transition (SFOPT). We classify the SFOPT into four types: slight supercooling, mild supercooling, strong supercooling, and ultra supercooling. Using different characteristic temperatures, length scales and bubble wall velocities, the corresponding gravitational wave (GW) spectra are investigated in details. We emphasize the essential importance of using the correct characteristic temperature and length scale when the phase transition dynamics and GW spectra are calculated. Especially, for strong supercooling and ultra supercooling cases, there are obvious differences of the phase transition strength and GW spectra between the results calculated at the nucleation temperature and those derived at the percolation temperature. For ultra supercooling case, we propose a criterion to quantify whether the phase transition can terminate. Besides the model-independent discussions, we also study three representative models as concrete examples to clearly show the subtle points therein.


title: Phonon-Phonon Interactions in Strongly Bonded Solids: Selection Rules and Higher-Order Processes

abstract: We show that the commonly used lowest-order theory of phonon-phonon interactions frequently fails to accurately describe the anharmonic phonon decay rates and thermal conductivity ($\kappa$), even among strongly bonded crystals. Applying a first principles theory that includes both the lowest-order three-phonon and the higher-order four-phonon processes to seventeen zinc blende semiconductors, we find that the lowest-order theory drastically overestimates the measured $\kappa$ for many of these materials, while inclusion of four-phonon scattering gives significantly improved agreement with measurements. We have identified new selection rules on three-phonon processes that help explain many of these failures in terms of anomalously weak anharmonic phonon decay rates predicted by the lowest-order theory competing with four-phonon processes. We also show that zinc blende compounds containing boron (B), carbon (C) or nitrogen (N) atoms have exceptionally weak four-phonon scattering, much weaker than in compounds that do not contain B, C or N atoms. This new understanding helps explain the ultrahigh $\kappa$ in several technologically important materials like cubic boron arsenide, boron phosphide and silicon carbide. At the same time, it not only makes the possibility of achieving high $\kappa$ in materials without B, C or N atoms unlikely, but it also suggests that it may be necessary to include four-phonon processes in many future studies. Our work gives new insights into the nature of anharmonic processes in solids and demonstrates the broad importance of higher-order phonon-phonon interactions in assessing the thermal properties of materials.


title: Degenerations of representations and thin triangles

abstract: There is a compactification of the space of representations of a finitely generated group into the groups of isometries of all spaces with $\Delta$-thin triangles. The ideal points are actions on $\mathbb R$-trees. It is a geometric reformulation and extension of the Culler-Morgan-Shalen theory concerning limits of representations into $\operatorname{SL}(2,{\mathbb C})$ and more generally $\operatorname{O}(n, 1)$. This paper was written and circulated in the early 90's, but never published.


title: Bosonic quantum communication across arbitrarily high loss channels

abstract: A general attenuator $\Phi_{\lambda, \sigma}$ is a bosonic quantum channel that acts by combining the input with a fixed environment state $\sigma$ in a beam splitter of transmissivity $\lambda$. If $\sigma$ is a thermal state the resulting channel is a thermal attenuator, whose quantum capacity vanishes for $\lambda\leq 1/2$. We study the quantum capacity of these objects for generic $\sigma$, proving a number of unexpected results. Most notably, we show that for any arbitrary value of $\lambda>0$ there exists a suitable single-mode state $\sigma(\lambda)$ such that the quantum capacity of $\Phi_{\lambda,\sigma(\lambda)}$ is larger than a universal constant $c>0$. Our result holds even when we fix an energy constraint at the input of the channel, and implies that quantum communication at a constant rate is possible even in the limit of arbitrarily low transmissivity, provided that the environment state is appropriately controlled. We also find examples of states $\sigma$ such that the quantum capacity of $\Phi_{\lambda,\sigma}$ is not monotonic in $\lambda$. These findings may have implications for the study of communication lines running across integrated optical circuits, of which general attenuators provide natural models.


title: Contaminants removal and bacterial activity enhancement along the flow path of constructed wetland microbial fuel cells

abstract: Microbial fuel cells implemented in constructed wetlands (CW-MFCs), albeit a relatively new technology still under study, have shown to improve treatment efficiency of urban wastewater. So far the vast majority of CW-MFC systems investigated were designed as lab-scale systems working under rather unrealistic hydraulic conditions using synthetic wastewater. The main objective of this work was to quantify CW-MFCs performance operated under different conditions in a more realistic setup using meso-scale systems with horizontal flow fed with real urban wastewater. Operational conditions tested were organic loading rate (4.9+-1.6, 6.7+-1.4 and 13.6+-3.2 g COD/m2.day) and hydraulic regime (continuous vs intermittent feeding) as well as different electrical connections: CW control (conventional CW without electrodes), open-circuit CW-MFC (external circuit between anode and cathode not connected) and closed-circuit CW-MFC (external circuit connected). Eight horizontal subsurface flow CWs were operated for about four months. Each wetland consisted of a PVC reservoir of 0.193 m2 filled with 4/8 mm granitic riverine gravel. All wetlands had intermediate sampling points for gravel and interstitial liquid sampling. The CW-MFCs were designed as three MFCs incorporated one after the other along the flow path of the CWs. Results showed no significant differences between tested organic loading rates, hydraulic regimes or electrical connections, however, on average, systems operated in closed-circuit CW-MFC mode under continuous flow outperformed the other experimental conditions. Closed-circuit CW-MFC compared to conventional CW control systems showed around 5% and 22% higher COD and ammonium removal, respectively. Correspondingly, overall bacteria activity, as measured by the fluorescein diacetate technique, was higher (4% to 34%) in closed-circuit systems when compared to CW control systems.


title: Normalized and Geometry-Aware Self-Attention Network for Image Captioning

abstract: Self-attention (SA) network has shown profound value in image captioning. In this paper, we improve SA from two aspects to promote the performance of image captioning. First, we propose Normalized Self-Attention (NSA), a reparameterization of SA that brings the benefits of normalization inside SA. While normalization is previously only applied outside SA, we introduce a novel normalization method and demonstrate that it is both possible and beneficial to perform it on the hidden activations inside SA. Second, to compensate for the major limit of Transformer that it fails to model the geometry structure of the input objects, we propose a class of Geometry-aware Self-Attention (GSA) that extends SA to explicitly and efficiently consider the relative geometry relations between the objects in the image. To construct our image captioning model, we combine the two modules and apply it to the vanilla self-attention network. We extensively evaluate our proposals on MS-COCO image captioning dataset and superior results are achieved when comparing to state-of-the-art approaches. Further experiments on three challenging tasks, i.e. video captioning, machine translation, and visual question answering, show the generality of our methods.


title: Classification of topological phases with finite internal symmetries in all dimensions

abstract: We develop a mathematical theory of symmetry protected trivial (SPT) orders and anomaly-free symmetry enriched topological (SET) orders in all dimensions via two different approaches with an emphasis on the second approach. The first approach is to gauge the symmetry in the same dimension by adding topological excitations as it was done in the 2d case, in which the gauging process is mathematically described by the minimal modular extensions of unitary braided fusion 1-categories. This 2d result immediately generalizes to all dimensions except in 1d, which is treated with special care. The second approach is to use the 1-dimensional higher bulk of the SPT/SET order and the boundary-bulk relation. This approach also leads us to a precise mathematical description and a classification of SPT/SET orders in all dimensions. The equivalence of these two approaches, together with known physical results, provides us with many precise mathematical predictions.

categories: math-ph cond-mat.str-el hep-th math.CT math.MP math.QA


title: Highly-efficient terahertz radiation generated by surface electrons from laser-foil interactions

abstract: A novel scheme for generating powerful terahertz (THz) radiation based on laser-solid interactions is proposed. When a $p$-polarized femtosecond laser impinges obliquely on a plane solid target and the target partially blocks the laser energy, surface electrons are extracted out and accelerated by the laser fields, forming a low-divergence electron beam. A half-cycle THz radiation pulse is emitted simultaneously as the beam passes by the edge of the target, due to coherent diffraction radiation. Our particle-in-cell simulations show that the relativistic THz pulse can have an energy of a few tens of millijoule and the conversion efficiency can be over 1$\%$ with existing $\sim$J level femtosecond laser sources.


title: Discrete Hirota reductions associated with the lattice KdV equation

abstract: We study the integrability of a family of birational maps obtained as reductions of the discrete Hirota equation, which are related to travelling wave solutions of the lattice KdV equation. In particular, for reductions corresponding to waves moving with rational speed N/M on the lattice, where N,M are coprime integers, we prove the Liouville integrability of the maps when N + M is odd, and prove various properties of the general case. There are two main ingredients to our construction: the cluster algebra associated with each of the Hirota bilinear equations, which provides invariant (pre)symplectic and Poisson structures; and the connection of the monodromy matrices of the dressing chain with those of the KdV travelling wave reductions.


title: Testing the isotropy of the Dark Energy Survey's extreme trans-Neptunian objects

abstract: We test whether the population of "extreme" trans-Neptunian objects (eTNOs) detected in the Y4 Dark Energy Survey (DES) data exhibit azimuthal asymmetries which might be evidence of gravitational perturbations from an unseen super-Earth in a distant orbit. By rotating the orbits of the detected eTNOs, we construct a synthetic population which, when subject to the DES selection function, reproduces the detected distribution of eTNOs in the orbital elements $a,e,$ and $i$ as well as absolute magnitude $H$, but has uniform distributions in mean anomaly $M$, longitude of ascending node $\Omega,$ and argument of perihelion $\omega.$ We then compare the detected distributions in each of $\Omega, \omega,$ and $\varpi\equiv\Omega+\omega$ to those expected from the isotropic population, using Kuiper's variant of the Kolmogorov-Smirnov test. The three angles are tested for each of 4 definitions of the eTNO population, choosing among $a>(150,250)$ AU and perihelion $q>(30,37)$ AU. These choices yield 3--7 eTNOs in the DES Y4 sample. Among the twelve total tests, two have the likelihood of drawing the observed angles from the isotropic population at $p<0.05.$ The 3 detections at $a>250, q>37$ AU, and the 4 detections at $a>250, q>30$ AU, have $\Omega$ distribution with $p=0.03$ of coming from the isotropic construction, but this is not strong evidence of anisotropy given the 12 different tests. The DES data taken on their own are thus consistent with azimuthal isotropy and do not require a "Planet 9" hypothesis. The limited sky coverage and object count mean, however, that the DES data by no means falsify this hypothesis.


title: Fast proximal algorithms for nonsmooth convex optimization

abstract: In the lines of our approach in \cite{Ouorou2019}, where we exploit Nesterov fast gradient concept \cite{Nesterov1983} to the Moreau-Yosida regularization of a convex function, we devise new proximal algorithms for nonsmooth convex optimization. These algorithms need no bundling mechanism to update the stability center while preserving the complexity estimates established in \cite{Ouorou2019}. We report some preliminary computational results on some academic test problem to give a first estimate of their performance in relation with the classical proximal bundle algorithm.


title: The $p$-Zassenhaus Filtration of a Free Profinite Group and Shuffle Relations

abstract: For a prime number $p$ and a free profinite group $S$ on the basis $X$, let $S_{(n,p)}$, $n=1,2,\ldots,$ be the $p$-Zassenhaus filtration of $S$. For $p>n$, we give a word-combinatorial description of the cohomology group $H^2(S/S_{(n,p)},\mathbb{Z}/p)$ in terms of the shuffle algebra on $X$. We give a natural linear basis for this cohomology group, which is constructed by means of unitriangular representations arising from Lyndon words.


title: RAB: Provable Robustness Against Backdoor Attacks

abstract: Recent studies have shown that deep neural networks (DNNs) are vulnerable to adversarial attacks, including evasion and backdoor (poisoning) attacks. On the defense side, there have been intensive efforts on improving both empirical and provable robustness against evasion attacks; however, the provable robustness against backdoor attacks still remains largely unexplored. In this paper, we focus on certifying the machine learning model robustness against general threat models, especially backdoor attacks. We first provide a unified framework via randomized smoothing techniques and show how it can be instantiated to certify the robustness against both evasion and backdoor attacks. We then propose the first robust training process, RAB, to smooth the trained model and certify its robustness against backdoor attacks. We prove the robustness bound for machine learning models trained with RAB and prove that our robustness bound is tight. In addition, we theoretically show that it is possible to train the robust smoothed models efficiently for simple models such as K-nearest neighbor classifiers, and we propose an exact smooth-training algorithm that eliminates the need to sample from a noise distribution for such models. Empirically, we conduct comprehensive experiments for different machine learning (ML) models such as DNNs, support vector machines, and K-NN models on MNIST, CIFAR-10, and ImageNette datasets and provide the first benchmark for certified robustness against backdoor attacks. In addition, we evaluate K-NN models on a spambase tabular dataset to demonstrate the advantages of the proposed exact algorithm. Both the theoretic analysis and the comprehensive evaluation on diverse ML models and datasets shed light on further robust learning strategies against general training time attacks.


title: Asteroseismic sensitivity to internal rotation along the red-giant branch

abstract: Transport of angular momentum in stellar interiors is currently not well understood. Asteroseismology can provide us with estimates of internal rotation of stars and thereby advances our understanding of angular momentum transport. We can measure core-rotation rates in red-giant stars and we can place upper bounds on surface-rotation rates using measurements of dipole ($l=1$) modes. Here, we aim to determine the theoretical sensitivity of modes of different spherical degree towards the surface rotation. Additionally, we aim to identify modes that can potentially add sensitivity at intermediate radii. We used asteroseismic rotational inversions to probe the internal stellar rotation profiles in red-giant models from the base of the red-giant branch up to the luminosity bump. We used the inversion method of multiplicative optimally localised averages (MOLA) to assess how well internal and surface rotation rates can be recovered from different mode sets and different synthetic rotation profiles. We confirm that dipole mixed modes are sufficient to set constraints on the average core-rotation rates in red giants. However, surface-rotation rates estimated with only dipole mixed modes are contaminated by the core rotation. We show that the sensitivity to the surface rotation decreases from the base of the red-giant branch until it reaches a minimum at 0.6-0.8$L_\text{bump}$ due to a glitch in the buoyancy frequency. Thereafter a narrow range of increased surface sensitivity just below the bump luminosity exists. Quadrupole and octopole modes have more sensitivity in the outer parts of the star. If observed, quadrupole and octopole modes enable us to distinguish between differential and solid body rotation in the convection zone. To obtain accurate estimates of rotation rates at intermediate radii, acoustic oscillation modes with a spherical degree of $l\approx10$ are needed.


title: A Lefschetz theorem for crystalline representations

abstract: As a corollary of nonabelian Hodge theory, Simpson proved a strong Lefschetz theorem for complex polarized variations of Hodge structure. We show an arithmetic analog. Our primary technique is $p$-adic nonabelian Hodge theory.


title: Overinterpretation reveals image classification model pathologies

abstract: Image classifiers are typically scored on their test set accuracy, but high accuracy can mask a subtle type of model failure. We find that high scoring convolutional neural networks (CNNs) on popular benchmarks exhibit troubling pathologies that allow them to display high accuracy even in the absence of semantically salient features. When a model provides a high-confidence decision without salient supporting input features, we say the classifier has overinterpreted its input, finding too much class-evidence in patterns that appear nonsensical to humans. Here, we demonstrate that neural networks trained on CIFAR-10 and ImageNet suffer from overinterpretation, and we find models on CIFAR-10 make confident predictions even when 95% of input images are masked and humans cannot discern salient features in the remaining pixel-subsets. We introduce Batched Gradient SIS, a new method for discovering sufficient input subsets for complex datasets, and use this method to show the sufficiency of border pixels in ImageNet for training and testing. Although these patterns portend potential model fragility in real-world deployment, they are in fact valid statistical patterns of the benchmark that alone suffice to attain high test accuracy. Unlike adversarial examples, overinterpretation relies upon unmodified image pixels. We find ensembling and input dropout can each help mitigate overinterpretation.


title: Primordial flat frame -- a new view on inflation

abstract: Models of inflationary cosmology admit a choice of the metric for which the geometry of homogeneous isotropic solutions becomes flat Minkowski space in the infinite past. In this primordial flat frame all mass scales vanish in the infinite past and quantum scale symmetry is realized. The cosmological evolution is dominantly described by the slow increase of a scalar field which sets the scale of all masses. We construct the primordial flat frame for standard models of inflation as Starobinsky inflation or chaotic inflation. In particular, we discuss the evolution of inhomogeneous solutions in the neighborhood of the homogeneous isotropic background solution and their relation to the observable primordial fluctuation spectrum. If the propagators for the graviton and scalar field remain regular, our observed inhomogeneous Universe can be extrapolated back to the infinite past in physical time. In this case there is no physical big-bang singularity -- the latter reflects only a singular choice of ``field coordinates". Independently of the issue of singularity the primordial flat frame offers a new view on the physical properties of the inflationary universe, which can be characterized as a very slowly evolving almost empty vacuum state with approximate scale symmetry.


title: Quiver Yangian from Crystal Melting

abstract: We find a new infinite class of infinite-dimensional algebras acting on BPS states for non-compact toric Calabi-Yau threefolds. In Type IIA superstring compactification on a toric Calabi-Yau threefold, the D-branes wrapping holomorphic cycles represent the BPS states, and the fixed points of the moduli spaces of BPS states are described by statistical configurations of crystal melting. Our algebras are "bootstrapped" from the molten crystal configurations, hence they act on the BPS states. We discuss the truncation of the algebra and its relation with D4-branes. We illustrate our results in many examples, with and without compact $4$-cycles.


title: Formal Synthesis of Lyapunov Neural Networks

abstract: We propose an automatic and formally sound method for synthesising Lyapunov functions for the asymptotic stability of autonomous non-linear systems. Traditional methods are either analytical and require manual effort or are numerical but lack of formal soundness. Symbolic computational methods for Lyapunov functions, which are in between, give formal guarantees but are typically semi-automatic because they rely on the user to provide appropriate function templates. We propose a method that finds Lyapunov functions fully automatically$-$using machine learning$-$while also providing formal guarantees$-$using satisfiability modulo theories (SMT). We employ a counterexample-guided approach where a numerical learner and a symbolic verifier interact to construct provably correct Lyapunov neural networks (LNNs). The learner trains a neural network that satisfies the Lyapunov criteria for asymptotic stability over a samples set; the verifier proves via SMT solving that the criteria are satisfied over the whole domain or augments the samples set with counterexamples. Our method supports neural networks with polynomial activation functions and multiple depth and width, which display wide learning capabilities. We demonstrate our method over several non-trivial benchmarks and compare it favourably against a numerical optimisation-based approach, a symbolic template-based approach, and a cognate LNN-based approach. Our method synthesises Lyapunov functions faster and over wider spatial domains than the alternatives, yet providing stronger or equal guarantees.


title: Projection and rescaling algorithm for finding maximum support solutions to polyhedral conic systems

abstract: We propose a simple projection and rescaling algorithm that finds maximum support solutions to the pair of feasibility problems \[ \text{find} \; x\in L\cap\mathbb{R}^n_{+} \;\;\;\; \text{ and } \; \;\;\;\; \text{find} \; \hat x\in L^\perp\cap\mathbb{R}^n_{+}, \] where $L$ is a linear subspace of $\mathbb{R}^n$ and $L^\perp$ is its orthogonal complement. The algorithm complements a basic procedure that involves only projections onto $L$ and $L^\perp$ with a periodic rescaling step. The number of rescaling steps and thus overall computational work performed by the algorithm are bounded above in terms of a condition measure of the above pair of problems. Our algorithm is a natural but significant extension of a previous projection and rescaling algorithm that finds a solution to the problem \[ \text{find} \; x\in L\cap\mathbb{R}^n_{++} \] when this problem is feasible. As a byproduct of our new developments, we obtain a sharper analysis of the projection and rescaling algorithm in the latter special case.


title: Morse Quasiflats II

abstract: This is the second in a two part series of papers concerning Morse quasiflats - higher dimensional analogs of Morse quasigeodesics. Our focus here is on their asymptotic structure. In metric spaces with convex geodesic bicombings, we prove asymptotic conicality, uniqueness of tangent cones at infinity and Euclidean volume growth rigidity for Morse quasiflats. Moreover, we provide some immediate consequences.


title: Step-wise target controllability identifies dysregulated pathways of macrophage networks in multiple sclerosis

abstract: Identifying the nodes that have the potential to influence the state of a network is a relevant question for many complex systems. In many applications it is often essential to test the ability of an individual node to control a specific target subset of the network. In biological networks, this might provide precious information on how single genes regulate the expression of specific groups of molecules in the cell. Taking into account these constraints, we propose an optimized heuristic based on the Kalman rank condition to quantify the centrality of a node as the number of target nodes it can control. By introducing a hierarchy among the nodes in the target set, and performing a step-wise research, we ensure for sparse and directed networks the identification of a controllable driver-target configuration in a significantly reduced space and time complexity. We show how the method works for simple network configurations, then we use it to characterize the inflammatory pathways in molecular gene networks associated with macrophage dysfunction in patients with multiple sclerosis. Results indicate that the targeted secreted molecules can in general be controlled by a large number of driver nodes (51%) involved in different cell functions, i.e. sensing, signaling and transcription. However, during the inflammatory response only a moderate fraction of all the possible driver-target pairs are significantly coactivated, as measured by gene expression data obtained from human blood samples. Notably, they differ between multiple sclerosis patients and healthy controls, and we find that this is related to the presence of dysregulated genes along the controllable walks. Our method, that we name step-wise target controllability, represents a practical solution to identify controllable driver-target configurations in directed complex networks and test their relevance from a functional perspective.


title: Idoneal genera and K3 surfaces covering an Enriques surface

abstract: Idoneal genera are a generalization of Euler's idoneal numbers. We enumerate all idoneal genera by means of the Smith--Minkowski--Siegel mass formula. As an application, we classify transcendental lattices of K3 surfaces covering an Enriques surface.


title: Automatically Proving Microkernels Free from Privilege Escalation from their Executable

abstract: Operating system kernels are the security keystone of most computer systems, as they provide the core protection mechanisms. Kernels are in particular responsible for their own security, i.e. they must prevent untrusted user tasks from reaching their level of privilege. We demonstrate that proving such absence of privilege escalation is a pre-requisite for any definitive security proof of the kernel. While prior OS kernel formal verifications were performed either on source code or crafted kernels, with manual or semi-automated methods requiring significant human efforts in annotations or proofs, we show that it is possible to compute such kernel security proofs using fully-automated methods and starting from the executable code of an existing microkernel with no modification, thus formally verifying absence of privilege escalation with high confidence for a low cost. We applied our method on two embedded microkernels, including the industrial kernel AnonymOS: with only 58 lines of annotation and less than 10 minutes of computation, our method finds a vulnerability in a first (buggy) version of AnonymOS and verifies absence of privilege escalation in a second (secure) version.


title: Poisson quasi-Nijenhuis manifolds and the Toda system

abstract: The notion of Poisson quasi-Nijenhuis manifold generalizes that of Poisson-Nijenhuis manifold. The relevance of the latter in the theory of completely integrable systems is well established since the birth of the bi-Hamiltonian approach to integrability. In this note, we discuss the relevance of the notion of Poisson quasi-Nijenhuis manifold in the context of finite-dimensional integrable systems. Generically (as we show by an example with $3$ degrees of freedom) the Poisson quasi-Nijenhuis structure is largely too general to ensure Liouville integrability of a system. However, we prove that the closed (or periodic) $n$-particle Toda lattice can be framed in such a geometrical structure, and its well-known integrals of the motion can be obtained as spectral invariants of a "quasi-Nijenhuis recursion operator", that is, a tensor field $N$ of type $(1,1)$ defined on the phase space of the lattice. This example and some of its generalizations are used to understand whether one can define in a reasonable sense a notion of {\em involutive\} Poisson quasi-Nijenhuis manifold. A geometrical link between the open (or non periodic) and the closed Toda systems is also framed in the context of a general scheme connecting Poisson quasi-Nijenhuis and Poisson-Nijenhuis manifolds.


title: A Real Polynomial for Bipartite Graph Minimum Weight Perfect Matchings

abstract: In a recent paper, Beniamini and Nisan gave a closed-form formula for the unique multilinear polynomial for the Boolean function determining whether a given bipartite graph $G \subseteq K_{n,n}$ has a perfect matching, together with an efficient algorithm for computing the coefficients of the monomials of this polynomial. We give the following generalization: Given an arbitrary non-negative weight function $w$ on the edges of $K_{n,n}$, consider its set of minimum weight perfect matchings. We give the real multilinear polynomial for the Boolean function which determines if a graph $G \subseteq K_{n,n}$ contains one of these minimum weight perfect matchings.


title: Delay time of waves performing L\'evy walks in 1D random media

abstract: The time that waves spend inside 1D random media with the possibility of performing L\'evy walks is experimentally and theoretically studied. The dynamics of quantum and classical wave diffusion has been investigated in canonical disordered systems via the delay time. We show that a wide class of disorder--L\'evy disorder--leads to strong random fluctuations of the delay time; nevertheless, some statistical properties such as the tail of the distribution and the average of the delay time are insensitive to L\'evy walks. Our results reveal a universal character of wave propagation that goes beyond standard Brownian wave-diffusion.


title: Scalable integrated single-photon source

abstract: Photonic qubits are key enablers for quantum-information processing deployable across a distributed quantum network. An on-demand and truly scalable source of indistinguishable single photons is the essential component enabling high-fidelity photonic quantum operations. A main challenge is to overcome noise and decoherence processes in order to reach the steep benchmarks on generation efficiency and photon indistinguishability required for scaling up the source. We report on the realization of a deterministic single-photon source featuring near-unity indistinguishability using a quantum dot in an 'on-chip' planar nanophotonic waveguide circuit. The device produces long strings of $>100$ single photons without any observable decrease in the mutual indistinguishability between photons. A total generation rate of $122$ million photons per second is achieved corresponding to an 'on-chip' source efficiency of $84 \%$. These specifications of the single-photon source are benchmarked for boson sampling and found to enable scaling into the regime of quantum advantage.


title: Parameter estimation for discretely sampled stochastic heat equation driven by space-only noise

abstract: We derive consistent and asymptotically normal estimators for the drift and volatility parameters of the stochastic heat equation driven by an additive space-only white noise when the solution is sampled discretely in the physical domain. We consider both the full space and the bounded domain. We establish the exact spatial regularity of the solution, which in turn, using power-variation arguments, allows building the desired estimators. We show that naive approximations of the derivatives appearing in the power-variation based estimators may create nontrivial biases, which we compute explicitly. The proofs are rooted in Malliavin-Stein's method.


title: Multilevel Emulation for Stochastic Computer Models with Application to Large Offshore Wind farms

abstract: Renewable energy projects, such as large offshore wind farms, are critical to achieving low-emission targets set by governments. Stochastic computer models allow us to explore future scenarios to aid decision making whilst considering the most relevant uncertainties. Complex stochastic computer models can be prohibitively slow and thus an emulator may be constructed and deployed to allow for efficient computation. We present a novel heteroscedastic Gaussian Process emulator which exploits cheap approximations to a stochastic offshore wind farm simulator. We also conduct a probabilistic sensitivity analysis to understand the influence of key parameters in the wind farm model which will help us to plan a probability elicitation in the future.


title: A Sub-Solar Metallicity Progenitor for Cassiopeia A, the remnant of a Type IIb Supernova

abstract: We report, for the first time, the detection of the Mn-K$\alpha$ line in the Type IIb supernova (SN IIb) remnant, Cassiopeia A. Manganese ($^{55}$Mn after decay of $^{55}$Co), a neutron-rich element, together with chromium ($^{52}$Cr after decay of $^{52}$Fe), is mainly synthesized at the explosive incomplete Si burning regime. Therefore, the Mn/Cr mass ratio with its neutron excess reflects the neutronization at the relevant burning layer during the explosion. Chandra's archival X-ray data of Cassiopeia A indicate a low Mn/Cr mass ratio with values in the range 0.10--0.66, which, when compared to one-dimensional SN explosion models, requires that the electron fraction be 0.4990 $\lesssim Y_{\rm e} \lesssim$ 0.5 at the incomplete Si burning layer. An explosion model assuming a solar-metallicity progenitor with a typical explosion energy ($1 \times 10^{51}$ erg) fails to reproduce such a high electron fraction. In such models, the explosive Si-burning regime extends only to the Si/O layer established during the progenitor's hydrostatic evolution; the $Y_e$ in the Si/O layer is lower than the value required by our observational constraints. We can satisfy the observed Mn/Cr mass ratio if the explosive Si-burning regime were to extend into the O/Ne hydrostatic layer, which has a higher $Y_{\rm e}$. This would require an energetic ($> 2 \times 10^{51}$ erg) and/or asymmetric explosion of a sub-solar metallicity progenitor ($Z \lesssim 0.5Z_{\odot}$) for Cassiopeia A. The low initial metallicity can be used to rule out a single-star progenitor, leaving the possibility of a binary progenitor with a compact companion (white dwarf, neutron star or black hole). We discuss the detectability of X-rays from Bondi accretion onto such a compact companion around the explosion site. We also discuss other possible mass-loss scenarios for the progenitor system of Cassiopeia A.


title: RF-Rhythm: Secure and Usable Two-Factor RFID Authentication

abstract: Passive RFID technology is widely used in user authentication and access control. We propose RF-Rhythm, a secure and usable two-factor RFID authentication system with strong resilience to lost/stolen/cloned RFID cards. In RF-Rhythm, each legitimate user performs a sequence of taps on his/her RFID card according to a self-chosen secret melody. Such rhythmic taps can induce phase changes in the backscattered signals, which the RFID reader can detect to recover the user's tapping rhythm. In addition to verifying the RFID card's identification information as usual, the backend server compares the extracted tapping rhythm with what it acquires in the user enrollment phase. The user passes authentication checks if and only if both verifications succeed. We also propose a novel phase-hopping protocol in which the RFID reader emits Continuous Wave (CW) with random phases for extracting the user's secret tapping rhythm. Our protocol can prevent a capable adversary from extracting and then replaying a legitimate tapping rhythm from sniffed RFID signals. Comprehensive user experiments confirm the high security and usability of RF-Rhythm with false-positive and false-negative rates close to zero.


title: Interpreting solar wind turbulent spectra beyond Taylor's Hypothesis

abstract: In this paper we apply a methodology, recently proposed by Bourouaine and Perez (2019) (BP19), to interpret solar-wind turbulent power spectra beyond Taylor approximation (TA). The turbulent power spectra were measured using \emph{Helios} spacecraft data near 0.6 au. We use the model proposed in BP19 to reproduce the field-perpendicular power spectrum $E(k_\perp)$ of anti-sunward Alfv\'enic fluctuations in the plasma frame (where $k_\perp$ is the field-perpendicular wavenumber) from the corresponding measured frequency power spectrum $P_{\rm sc}(\omega,\theta_b)$ along the sampling angle $\theta_b$, which is the angle between the local magnetic field and the sampling direction. Here $\omega=2\pi f$ and $f$ is the frequency of the time signal. Interestingly enough, we found that for all corresponding measured frequency power spectrum $P_{\rm sc}(\omega,\theta_b)$ the reproduced field-perpendicular power spectrum $E(k_\perp)$ is the same and independent of the considered sampling angle $\theta_b$. This finding is consistent with the fact that the analyzed turbulence is strong and highly anisotropic with $k_\| \ll k_\perp$ (where $k_\|$ is the field-parallel wavenumber). Furthermore, for this specific time signal we found that the commonly used TA is still approximately valid with the important difference that a broadening in $k_\perp$ for each angular frequency $\omega$ is present. This broadening can be described in the context of the methodology proposed in BP19.


title: Utilizing Language Relatedness to improve Machine Translation: A Case Study on Languages of the Indian Subcontinent

abstract: In this work, we present an extensive study of statistical machine translation involving languages of the Indian subcontinent. These languages are related by genetic and contact relationships. We describe the similarities between Indic languages arising from these relationships. We explore how lexical and orthographic similarity among these languages can be utilized to improve translation quality between Indic languages when limited parallel corpora is available. We also explore how the structural correspondence between Indic languages can be utilized to re-use linguistic resources for English to Indic language translation. Our observations span 90 language pairs from 9 Indic languages and English. To the best of our knowledge, this is the first large-scale study specifically devoted to utilizing language relatedness to improve translation between related languages.


title: Hausdorff and packing dimensions and measures for nonlinear transversally non-conformal thin solenoids

abstract: We extend results by B. Hasselblatt, J. Schmeling in \emph{Dimension product structure of hyperbolic sets} (2004), and by the third author and K. Simon in \emph{Hausdorff and packing measures for solenoids} (2003), for $C^{1+\varepsilon}$ hyperbolic, (partially) linear solenoids $\Lambda$ over the circle embedded in $\mathbb{R}^3$ non-conformally attracting in the stable discs $W^s$ direction, to nonlinear ones. Under an assumption of transversality and assumptions on Lyapunov exponents for an appropriate Gibbs measure imposing \emph{thinness}, assuming also there is an invariant $C^{1+\varepsilon}$ strong stable foliation, we prove that Hausdorff dimension ${\rm HD}(\Lambda\cap W^s)$ is the same quantity $t_0$ for all $W^s$ and else ${\rm HD}(\Lambda)=t_0+1$. We prove also that for the packing measure $0<\Pi_{t_0}(\Lambda \cap W^s)<\infty$ but for Hausdorff measure ${\rm HM}_{t_0}(\Lambda\cap W^s)=0$ for all $W^s$. Also $0<\Pi_{1+t_0}(\Lambda) <\infty$ and ${\rm HM}_{1+t_0}(\Lambda)=0$. A technical part says that the holonomy along unstable foliation is locally Lipschitz, except for a set of unstable leaves whose intersection with every $W^s$ has measure ${\rm HM}_{t_0}$ equal to 0 and even Hausdorff dimension less than $t_0$. The latter holds due to a large deviations phenomenon.


title: Rotating black holes in $4D$ Einstein-Gauss-Bonnet gravity and its shadow

abstract: Recently there has been a surge of interest in regularizing, a $ D \to 4 $ limit of, the Einstein-Gauss-Bonnet (EGB) gravity, and the resulting regularized $4D$ EGB gravity has nontrivial dynamics. The theory admits spherically symmetric black holes generalizing the Schwarzschild black holes. We consider the rotating black hole in regularized $4D$ EGB gravity and discuss their horizon properties and shadow cast. The effects of the GB coupling parameter on the shape and size of shadows are investigated in the context of recent M87* observations from the EHT. Interestingly, for a given spin parameter, the apparent size of the shadow decreases and gets more distorted due to the GB coupling parameter. We find that within the finite parameter space, e.g. for $a=0.1M$, $\alpha\leq 0.00394M^2$, and within the current observational uncertainties, the rotating black holes of the $4D$ EGB gravity are consistent with the inferred features of M87* black hole shadow.


title: A Systematic Analysis of the Phase Lags Associated with the Type-C Quasi-periodic Oscillation in GRS 1915+105

abstract: We present a systematic analysis of the phase lags associated with the type-C QPOs in GRS 1915+105 using RXTE data. Our sample comprises of 620 RXTE observations with type-C QPOs ranging from ~0.4 Hz to ~6.3 Hz. Based on our analysis, we confirm that the QPO phase lags decrease with QPO frequency, and change sign from positive to negative at a QPO frequency of ~2 Hz. In addition, we find that the slope of this relation is significantly different between QPOs below and above 2 Hz. The relation between the QPO lags and QPO rms can be well fitted with a broken line: as the QPO lags go from negative to positive, the QPO rms first increases, reaching its maximum at around zero lag, and then decreases. The phase-lag behaviour of the subharmonic of the QPO is similar to that of the QPO fundamental, where the subharmonic lags decrease with subharmonic frequency and change sign from positive to negative at a subharmonic frequency of ~1 Hz; on the contrary, the second harmonic of the QPO shows a quite different phase-lag behaviour, where all the second harmonics show hard lags that remain more or less constant. For both the QPO and its (sub)harmonics, the slope of the lag-energy spectra shows a similar evolution with frequency as the average phase lags. This suggests that the lag-energy spectra drives the average phase lags. We discuss the possibility for the change in lag sign, and the physical origin of the QPO lags.


title: Faster Divergence Maximization for Faster Maximum Flow

abstract: In this paper we provide an algorithm which given any $m$-edge $n$-vertex directed graph with integer capacities at most $U$ computes a maximum $s$-$t$ flow for any vertices $s$ and $t$ in $m^{4/3+o(1)}U^{1/3}$ time. This improves upon the previous best running times of $m^{11/8+o(1)}U^{1/4}$ (Liu Sidford 2019), $\tilde{O}(m \sqrt{n} \log U)$ (Lee Sidford 2014), and $O(mn)$ (Orlin 2013) when the graph is not too dense or has large capacities. To achieve the results this paper we build upon previous algorithmic approaches to maximum flow based on interior point methods (IPMs). In particular, we overcome a key bottleneck of previous advances in IPMs for maxflow (M\k{a}dry 2013, M\k{a}dry 2016, Liu Sidford 2019), which make progress by maximizing the energy of local $\ell_2$ norm minimizing electric flows. We generalize this approach and instead maximize the divergence of flows which minimize the Bregman divergence distance with respect to the weighted logarithmic barrier. This allows our algorithm to avoid dependencies on the $\ell_4$ norm that appear in other IPM frameworks (e.g. Cohen M\k{a}dry Sankowski Vladu 2017, Axiotis M\k{a}dry Vladu 2020). Further, we show that smoothed $\ell_2$-$\ell_p$ flows (Kyng, Peng, Sachdeva, Wang 2019), which we previously used to efficiently maximize energy (Liu Sidford 2019), can also be used to efficiently maximize divergence, thereby yielding our desired runtimes. We believe both this generalized view of energy maximization and generalized flow solvers we develop may be of further interest.


title: Validation of secondary fluorescence excitation in quantitative X-ray fluorescence analysis of thin alloy films

abstract: X-ray fluorescence (XRF) analysis is a widely applied technique for the quantitative analysis of thin films up to the $\mu$m scale because of its non-destructive nature and because it is easily automated. When low uncertainties of the analytical results in the few percent range are required, the non-linear secondary fluorescence effect in multi-elemental samples may complicate an otherwise straightforward quantification, since it can easily exceed a relative contribution of 20%. The conventional solution, to rely on good performing reference samples, is hindered by their low availability, especially for thin film applications. To address this challenge, we demonstrate a flexible production method of multilayered, alloyed thin films with significant secondary fluorescence contributions. We use reference-free XRF analysis to validate the reliability of the physical model for secondary fluorescence, which includes a thorough uncertainty estimation. The investigated specimens are being qualified as calibration samples for XRF or other quantitative analyses.


title: Integrable Lattice Models and Holography

abstract: We study four-dimensional Chern-Simons theory on $D \times \mathbb{C}$ (where $D$ is a disk), which is understood to describe rational solutions of the Yang-Baxter equation from the work of Costello, Witten and Yamazaki. We find that the theory is dual to a boundary theory, that is a three-dimensional analogue of the two-dimensional chiral WZW model. This boundary theory gives rise to a current algebra that turns out to be an "analytically-continued" toroidal Lie algebra. In addition, we show how certain bulk correlation functions of two and three Wilson lines can be captured by boundary correlation functions of local operators in the three-dimensional WZW model. In particular, we reproduce the leading and subleading nontrivial contributions to the rational R-matrix purely from the boundary theory.

categories: hep-th cond-mat.stat-mech math-ph math.MP


title: GIQA: Generated Image Quality Assessment

abstract: Generative adversarial networks (GANs) have achieved impressive results today, but not all generated images are perfect. A number of quantitative criteria have recently emerged for generative model, but none of them are designed for a single generated image. In this paper, we propose a new research topic, Generated Image Quality Assessment (GIQA), which quantitatively evaluates the quality of each generated image. We introduce three GIQA algorithms from two perspectives: learning-based and data-based. We evaluate a number of images generated by various recent GAN models on different datasets and demonstrate that they are consistent with human assessments. Furthermore, GIQA is available to many applications, like separately evaluating the realism and diversity of generative models, and enabling online hard negative mining (OHEM) in the training of GANs to improve the results.


title: DELTAS: Depth Estimation by Learning Triangulation And densification of Sparse points

abstract: Multi-view stereo (MVS) is the golden mean between the accuracy of active depth sensing and the practicality of monocular depth estimation. Cost volume based approaches employing 3D convolutional neural networks (CNNs) have considerably improved the accuracy of MVS systems. However, this accuracy comes at a high computational cost which impedes practical adoption. Distinct from cost volume approaches, we propose an efficient depth estimation approach by first (a) detecting and evaluating descriptors for interest points, then (b) learning to match and triangulate a small set of interest points, and finally (c) densifying this sparse set of 3D points using CNNs. An end-to-end network efficiently performs all three steps within a deep learning framework and trained with intermediate 2D image and 3D geometric supervision, along with depth supervision. Crucially, our first step complements pose estimation using interest point detection and descriptor learning. We demonstrate state-of-the-art results on depth estimation with lower compute for different scene lengths. Furthermore, our method generalizes to newer environments and the descriptors output by our network compare favorably to strong baselines. Code is available at https://github.com/magicleap/DELTAS


title: NeRF: Representing Scenes as Neural Radiance Fields for View Synthesis

abstract: We present a method that achieves state-of-the-art results for synthesizing novel views of complex scenes by optimizing an underlying continuous volumetric scene function using a sparse set of input views. Our algorithm represents a scene using a fully-connected (non-convolutional) deep network, whose input is a single continuous 5D coordinate (spatial location $(x,y,z)$ and viewing direction $(\theta, \phi)$) and whose output is the volume density and view-dependent emitted radiance at that spatial location. We synthesize views by querying 5D coordinates along camera rays and use classic volume rendering techniques to project the output colors and densities into an image. Because volume rendering is naturally differentiable, the only input required to optimize our representation is a set of images with known camera poses. We describe how to effectively optimize neural radiance fields to render photorealistic novel views of scenes with complicated geometry and appearance, and demonstrate results that outperform prior work on neural rendering and view synthesis. View synthesis results are best viewed as videos, so we urge readers to view our supplementary video for convincing comparisons.