title: Max-value Entropy Search for Multi-Objective Bayesian Optimization with Constraints

abstract: We consider the problem of constrained multi-objective blackbox optimization using expensive function evaluations, where the goal is to approximate the true Pareto set of solutions satisfying a set of constraints while minimizing the number of function evaluations. For example, in aviation power system design applications, we need to find the designs that trade-off total energy and the mass while satisfying specific thresholds for motor temperature and voltage of cells. This optimization requires performing expensive computational simulations to evaluate designs. In this paper, we propose a new approach referred as {\em Max-value Entropy Search for Multi-objective Optimization with Constraints (MESMOC)} to solve this problem. MESMOC employs an output-space entropy based acquisition function to efficiently select the sequence of inputs for evaluation to uncover high-quality pareto-set solutions while satisfying constraints. We apply MESMOC to two real-world engineering design applications to demonstrate its effectiveness over state-of-the-art algorithms.


title: What Makes Agile Test Artifacts Useful? An Activity-Based Quality Model from a Practitioners' Perspective

abstract: Background: The artifacts used in Agile software testing and the reasons why these artifacts are used are fairly well-understood. However, empirical research on how Agile test artifacts are eventually designed in practice and which quality factors make them useful for software testing remains sparse. Aims: Our objective is two-fold. First, we identify current challenges in using test artifacts to understand why certain quality factors are considered good or bad. Second, we build an Activity-Based Artifact Quality Model that describes what Agile test artifacts should look like. Method: We conduct an industrial survey with 18 practitioners from 12 companies operating in seven different domains. Results: Our analysis reveals nine challenges and 16 factors describing the quality of six test artifacts from the perspective of Agile testers. Interestingly, we observed mostly challenges regarding language and traceability, which are well-known to occur in non-Agile projects. Conclusions: Although Agile software testing is becoming the norm, we still have little confidence about general do's and don'ts going beyond conventional wisdom. This study is the first to distill a list of quality factors deemed important to what can be considered as useful test artifacts.


title: The Well-Tempered Cosmological Constant: Fugue in B$^\flat$

abstract: Zero point fluctuations of quantum fields should generate a large cosmological constant energy density in any spacetime. How then can we have anything other than de Sitter space without fine tuning? Well tempering -- dynamical cancellation of the cosmological constant using degeneracy within the field equations -- can replace a large cosmological constant with a much lower energy state. Here we give an explicit mechanism to obtain a Minkowski solution, replacing the cosmological constant with zero, and testing its attractor nature and persistence through a vacuum phase transition. We derive the general conditions that Horndeski scalar-tensor gravity must possess, and evolve in a fugue of functions, to deliver nothing and make the universe be flat.


title: Symmetry and scaling limits for matching of implicit surfaces based on thin shell energies

abstract: In a recent paper by Iglesias, Rumpf and Scherzer (Found. Comput. Math. 18(4), 2018) a variational model for deformations matching a pair of shapes given as level set functions was proposed. Its main feature is the presence of anisotropic energies active only in a narrow band around the hypersurfaces that resemble the behavior of elastic shells. In this work we consider some extensions and further analysis of that model. First, we present a symmetric energy functional such that given two particular shapes, it assigns the same energy to any given deformation as to its inverse when the roles of the shapes are interchanged, and introduce the adequate parameter scaling to recover a surface problem when the width of the narrow band vanishes. Then, we obtain existence of minimizing deformations for the symmetric energy in classes of bi-Sobolev homeomorphisms for small enough widths, and prove a $\Gamma$-convergence result for the corresponding non-symmetric energies as the width tends to zero. Finally, numerical results on realistic shape matching applications demonstrating the effect of the symmetric energy are presented.


title: An application of tomographic PIV to investigate the spray-induced turbulence in a direct-injection engine

abstract: Fuel sprays produce high-velocity, jet-like flows that impart turbulence onto the ambient flow field. The spray-induced turbulence augments fuel-air mixing, which has a primary role in controlling pollutant formation and cyclic variability in engines. This paper presents tomographic particle image velocimetry (TPIV) measurements to analyse the 3D spray-induced turbulence during the intake stroke of a direct-injection engine. The spray produces a strong spray-induced jet in the far field, which travels through the cylinder and imparts turbulence onto the surrounding flow. Planar high-speed PIV measurements at 4.8 kHz are combined with TPIV at 3.3 Hz to evaluate spray particle distributions and validate TPIV measurements in the particle-laden flow. An uncertainty analysis is performed to assess the uncertainty associated with vorticity and strain rate components. TPIV analyses quantify the spatial domain of the turbulence in relation to the SIJ and describe how turbulent flow features such as turbulent kinetic energy, strain rate and vorticity evolve into the surrounding flow field. Access to the full tensors facilitate the evaluation of turbulence for individual spray events. TPIV images reveal the presence of strong shear layers (visualized by high S magnitudes) and pockets of elevated vorticity along the immediate boundary of the SIJ. Values are extracted from spatial domains extending in 1mm increments from the SIJ. Turbulence levels are greatest within the 0-1mm region from the SIJ boarder and dissipate with radial distance. Individual strain rate and vorticity components are analyzed in detail to describe the relationship between local strain rates and 3D vortical structures produced within strong shear layers of the SIJ. Analyses are intended to understand the flow features responsible for rapid fuel-air mixing and provide valuable data for the development of numerical models.


title: Survival Estimation for Missing not at Random Censoring Indicators based on Copula Models

abstract: In the presence of right-censored data with covariates, the conditional Kaplan-Meier estimator (also known as the Beran estimator) consistently estimates the conditional survival function of the random follow-up for the event of interest. However, a necessary condition is the unambiguous knowledge of whether each individual is censored or not, which may be incomplete in practice. We therefore propose a study of the Beran estimator when the censoring indicators are generic random variables and discuss necessary conditions for the efficiency of the Beran estimator. From this, we provide a new estimator for the conditional survival function with missing not at random (MNAR) censoring indicators based on a conditional copula model for the missingness mechanism. In addition to the theoretical results, we illustrate how the estimators work for small samples through a simulation study and show their practical applicability by analyzing synthetic and real data.


title: Parallel selective nuclear spin addressing for fast high-fidelity quantum gates

abstract: Due to their long coherence times, nuclear spins have gained considerable attention as physical qubits. Two-qubit gates between nuclear spins of distinct resonance frequencies can be mediated by electron spins, usually employing a sequence of electron-nuclear gates. Here we present a different approach inspired by, but not limited to, NV centers in diamond and discuss possible applications. To this end we generalize external electron spin control sequences for nuclear spin initialization and hyperpolarization to achieve the simultaneous control of distinct nuclear spins via an electron spin. This approach results in efficient entangling gates that, compared to standard techniques, reduce the gate time by more than 50% when the gate time is limited by off-resonant coupling to other spins, and by up to 22% when the gate time is limited by small electron-nuclear coupling.


title: The GW190521 Mass Gap Event and the Primordial Black Hole Scenario

abstract: The LIGO/Virgo Collaboration has recently observed GW190521, the first binary black hole merger with at least the primary component mass in the mass gap predicted by the pair-instability supernova theory. This observation disfavors the standard stellar-origin formation scenario for the heavier black hole, motivating alternative hypotheses. We show that GW190521 cannot be explained within the Primordial Black Hole (PBH) scenario if PBHs do not accrete during their cosmological evolution, since this would require an abundance which is already in tension with current constraints. On the other hand, GW190521 may have a primordial origin if PBHs accrete efficiently before the reionization epoch.


title: MIPGAN -- Generating Strong and High Quality Morphing Attacks Using Identity Prior Driven GAN

abstract: Face morphing attacks target to circumvent Face Recognition Systems (FRS) by employing face images derived from multiple data subjects (e.g., accomplices and malicious actors). Morphed images can be verified against contributing data subjects with a reasonable success rate, given they have a high degree of facial resemblance. The success of morphing attacks is directly dependent on the quality of the generated morph images. We present a new approach for generating strong attacks extending our earlier framework for generating face morphs. We present a new approach using an Identity Prior Driven Generative Adversarial Network, which we refer to as MIPGAN (Morphing through Identity Prior driven GAN). The proposed MIPGAN is derived from the StyleGAN with a newly formulated loss function exploiting perceptual quality and identity factor to generate a high quality morphed facial image with minimal artefacts and with high resolution. We demonstrate the proposed approach's applicability to generate strong morphing attacks by evaluating its vulnerability against both commercial and deep learning based Face Recognition System (FRS) and demonstrate the success rate of attacks. Extensive experiments are carried out to assess the FRS's vulnerability against the proposed morphed face generation technique on three types of data such as digital images, re-digitized (printed and scanned) images, and compressed images after re-digitization from newly generated MIPGAN Face Morph Dataset. The obtained results demonstrate that the proposed approach of morph generation poses a high threat to FRS.


title: Bayesian Perceptron: Towards fully Bayesian Neural Networks

abstract: Artificial neural networks (NNs) have become the de facto standard in machine learning. They allow learning highly nonlinear transformations in a plethora of applications. However, NNs usually only provide point estimates without systematically quantifying corresponding uncertainties. In this paper a novel approach towards fully Bayesian NNs is proposed, where training and predictions of a perceptron are performed within the Bayesian inference framework in closed-form. The weights and the predictions of the perceptron are considered Gaussian random variables. Analytical expressions for predicting the perceptron's output and for learning the weights are provided for commonly used activation functions like sigmoid or ReLU. This approach requires no computationally expensive gradient calculations and further allows sequential learning.


title: Strong Consistency and Thomas Decomposition of Finite Difference Approximations to Systems of Partial Differential Equations

abstract: For a wide class of polynomially nonlinear systems of partial differential equations we suggest an algorithmic approach that combines differential and difference algebra to analyze s(trong)-consistency of finite difference approximations. Our approach is applicable to regular solution grids. For the grids of this type we give a new definition of s-consistency for finite difference approximations which generalizes our definition given earlier for Cartesian grids. The algorithmic verification of s-consistency presented in the paper is based on the use of both differential and difference Thomas decomposition. First, we apply the differential decomposition to the input system, resulting in a partition of its solution space. Then, to the output subsystem that contains a solution of interest we apply a difference analogue of the differential Thomas decomposition which allows to check the s-consistency. For linear and some quasi-linear differential systems one can also apply difference \Gr bases for the s-consistency analysis. We illustrate our methods and algorithms by a number of examples, which include Navier-Stokes equations for viscous incompressible flow.

categories: cs.SC cs.NA math.AP math.NA math.RA physics.flu-dyn


title: The Generic Sudden Singularity in Brans-Dicke Theory

abstract: We construct a formal asymptotic series expansion for a general solution of the Brans-Dicke equations with a fluid source near a sudden singularity. This solution contains eleven independent arbitrary functions of the spatial coordinates as required by the Cauchy problem of the theory. We show that the solution is geodesically complete and has the character of a shock wave in the sudden asymptotic region. This solution is weak in the senses of Tipler and Krolak as in the corresponding case of general relativity.


title: Variation of TiO2/SiO2 mixed layers induced by different Xe+ ion energies

abstract: The broadening and optical parameters of TiO2/SiO2 transition layers depending on the ion energy have been investigated using the Rutherford Backscattering Spectrometry (RBS) and Ellipstrometry Spectroscopy (ES) methods. The TiO2/SiO2 samples were irradiated by Xe+ ions with energies of 100, 150, 200 and 250 keV. The depth profiles of the elements determined by the RBS spectra show the structure and thickness of the TiO2/SiO2 transition layers before and after implantation. We have found that the thickness of the transition region between the TiO2 and SiO2 layers increases with the increasing incident ion energy. This phenomenon indicates an increasing amount of atomic mixing at the TiO2/SiO2 interface. In addition, the variation of transition layers could be explained by defect depth profiles and ions energy transferred in the mixed layers by means of SRIM calculations. The thickness obtained from the RBS is in good agreement with that measured using the ES method. Based on these obtained results, we have also investigated the optical constants of implanted and non-implanted TiO2/SiO2 structures. The wave forms measured with varying incident angles suggest that the measurements were made close to near the main principle point. The yields of \psi and \delta bands vary at different incident angles, is associated with interference processes of the light reflected from the structures examined. The refractive index and the extinction coefficient were found to increase after implantation taking place up to 200-keV Xe and then decrease at 250 keV.


title: Computation Of Microbial Ecosystems in Time and Space (COMETS): An open source collaborative platform for modeling ecosystems metabolism

abstract: Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are also emerging as a valuable avenue for predicting, understanding and designing microbial communities. COMETS (Computation Of Microbial Ecosystems in Time and Space) was initially developed as an extension of dynamic flux balance analysis, which incorporates cellular and molecular diffusion, enabling simulations of multiple microbial species in spatially structured environments. Here we describe how to best use and apply the most recent version of this platform, COMETS 2, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, as well as several new biological simulation modules, including evolutionary dynamics and extracellular enzyme activity. COMETS 2 provides user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, and comprehensive documentation and tutorials, facilitating the use of COMETS for researchers at all levels of expertise with metabolic simulations. This protocol provides a detailed guideline for installing, testing and applying COMETS 2 to different scenarios, with broad applicability to microbial communities across biomes and scales.


title: Generation of zonal flows in convective systems by travelling thermal waves

abstract: This work addresses the effect of travelling thermal waves applied at the fluid layer surface, on the formation of global flow structures in 2D and 3D convective systems. For a broad range of Rayleigh numbers ($10^3\leq Ra \leq 10^7$) and thermal wave frequencies ($10^{-4}\leq \Omega \leq 10^{0}$), we investigate flows with and without imposed mean temperature gradients. Our results confirm that the travelling thermal waves can cause zonal flows, i.e. strong mean horizontal flows. We show that the zonal flows in diffusion dominated regimes are driven purely by the Reynolds stresses, always travelling retrograde, while in convection dominated regimes, mean flow advection, caused by tilted convection cells, becomes dominant, which generally leads to prograde mean zonal flows. By means of direct numerical simulations we validate theoretical predictions made for the diffusion dominated regime. Furthermore, we make use of the linear stability analysis and explain the existence of the tilted convection cell mode. Our extensive 3D simulations support the results for 2D flows and thus confirm the relevance of the findings for geopyhsical and astrophysical systems.


title: Jet quenching in the hadron gas: An exploratory study

abstract: In most calculations of hard particle suppression in heavy-ion reactions the hadronic stage has been neglected due to formation time arguments. Most of the hard particle shower exits the hot and dense medium before the system enters the hadronic evolution. In this contribution, a first assessment within the hadronic transport approach SMASH (Simulating Many Accelerated Strongly-interacting Hadrons) of rescattering effects on hard particles is presented. In particular, it is shown that the hadronic energy loss depends on the particle species as well as the energy of the probe. A parametrization for the $\langle \tilde{q} \rangle$ parameter as a function of temperature and particle energy is given for pions. Overall, major effects of the hadronic stage are expected in the transverse momentum range from 2-10 GeV and therefore jet sub-structure analysis and (hard-soft) correlation observables might be affected.


title: Tuning the pull-in instability of soft dielectric elastomers through loading protocols

abstract: Pull-in (or electro-mechanical) instability occurs when a drastic decrease in the thickness of a dielectric elastomer results in electrical breakdown, which limits the applications of dielectric devices. Here we derive the criteria for determining the pull-in instability of dielectrics actuated by different loading methods: voltage-control, charge-control, fixed pre-stress and fixed pre-stretch, by analyzing the free energy of the actuated systems. The Hessian criterion identifies a maximum in the loading curve beyond which the elastomer will stretch rapidly and lose stability, and can be seen as a path to failure. We present numerical calculations for neo-Hookean ideal dielectrics, and obtain the maximum allowable actuation stretch of a dielectric before failure by electrical breakdown. We find that applying a fixed pre-stress or a fixed pre-stretch to a charge-driven dielectric may decrease the stretchability of the elastomer, a scenario which is the opposite of what happens in the case of a voltage-driven dielectric. Results show that a reversible large actuation of a dielectric elastomer, free of the pull-in instability, can be achieved by tuning the actuation method.


title: Modified multiplicative decomposition model for tissue growth: Beyond the initial stress-free state

abstract: The multiplicative decomposition model is widely employed for predicting residual stresses and morphologies of biological tissues due to growth. However, it relies on the assumption that the tissue is initially in a stress-free state, which conflicts with the observations that any growth state of tissue is under a significant level of residual stresses that helps to maintain its ideal mechanical conditions. Here, we propose a modified multiplicative decomposition model where the initial state (or reference configuration) of biological tissue is endowed with residual stress instead of being stress-free. Releasing theoretically the initial residual stress, the initially stressed state is first transmitted into a virtual stress-free state, resulting in an initial elastic deformation. The initial virtual stress-free state subsequently grows to another counterpart with a growth deformation and the latter is further integrated into its natural configuration with an excessive elastic deformation that ensures tissue compatibility. With this decomposition, the total deformation may be expressed as the product of elastic deformation, growth deformation and initial elastic deformation, while the corresponding free energy density depends on the initial residual stress and the total deformation. We address three key issues: explicit expression of the free energy density,predetermination of the initial elastic deformation, and initial residual stress. Finally, we consider a tubular organ to demonstrate the effects of the proposed initial residual stress on stress distribution and on shape formation through an incremental stability analysis. Our results suggest that the initial residual stress exerts a major influence on the growth stress and the morphology of tissues.


title: Flame/flow dynamics at the piston surface of an IC engine measured by high-speed PLIF and PTV

abstract: Resolving fluid transport at engine surfaces is required to predict transient heat loss, which is becoming increasingly important for the development of high-efficiency internal combustion engines (ICE). The limited number of available investigations have focused on non-reacting flows near engine surfaces, while this work focuses on the near-wall flow field dynamics in response to a propagating flame front. Flow-field and flame distributions were measured simultaneously at kHz repetition rates using particle tracking velocimetry (PTV) and planar laser induced fluorescence (PLIF) of sulfur dioxide (SO2). Measurements were performed near the piston surface of an optically accessible engine operating at 800 rpm with homogeneous, stoichiometric isooctane-air mixtures. High-speed measurements reveal a strong interdependency between near-wall flow and flame development which also influences subsequent combustion. A conditional analysis is performed to analyze flame/flow dynamics at the piston surface for cycles with weak and strong flow velocities parallel to the surface. Faster flame propagation associated with higher velocities before ignition demonstrates a stronger flow acceleration ahead of the flame. Flow acceleration associated with an advancing flame front is a transient feature that strongly influences boundary layer development. The distance from the wall to 75% maximum velocity ({\delta}75) is analyzed to compare boundary layer development between fired and motored datasets. Decreases in {\delta}75 are strongly related to flow acceleration produced by an approaching flame front. Measurements reveal strong deviations of the boundary layer flow between fired and motored datasets, emphasizing the need to consider transient flow behavior when modelling boundary layer physics for reacting flows.


title: Effects of nerve bundle geometry on neurotrauma evaluation

abstract: Objective: We confirm that alteration of a neuron structure can induce abnormalities in signal propagation for nervous systems, as observed in brain damage. Here, we investigate the effects of geometrical changes and damage of a neuron structure in 2 scaled nerve bundle models, made of myelinated nerve fibers or unmyelinated nerve fibers. Methods: We propose a 3D finite element model of nerve bundles, combining a real-time full electromechanical coupling, a modulated threshold for spiking activation, and independent alteration of the electrical properties for each fiber. We then simulate mechanical compression and tension to induce damage at the membrane of a nerve bundle made of 4 fibers. We examine the resulting changes in strain and neural activity by considering in turn the cases of intact and traumatized nerve membranes. Results: Our results show lower strain and lower electrophysiological impairments in unmyelinated fibers than in myelinated fibers, higher deformation levels in larger bundles, and higher electrophysiological impairments in smaller bundles. Conclusion: We conclude that the insulation sheath of myelin constricts the membrane deformation and scatters plastic strains within the bundle, that larger bundles deform more than small bundles, and that small fibers tolerate a higher level of elongation before mechanical failure.


title: Characterization of Curvature positivity of Riemannian metrics on flat vector bundles

abstract: We give a characterization of Nakano positivity of Riemannian flat vector bundles over bounded domains $D\subset\mathbb{R}^n$ in terms of solvability of the $d$ equation with certain good $L^2$ estimate condition. As an application, we give an alternative proof of the matrix-valued Prekopa's theorem that is originally proved by Raufi. Our methods are inspired by the recent works of Deng-Ning-Wang-Zhou on characterization of Nakano positivity of Hermitian holomorphic vector bundles and positivity of direct image sheaves associated to holomorphic fibrations.


title: Online Estimation and Community Detection of Network Point Processes for Event Streams

abstract: A common goal in network modeling is to uncover the latent community structure present among nodes. For many real-world networks, the true connections consist of events arriving as streams, which are then aggregated to form edges, ignoring the dynamic temporal component. A natural way to take account of these temporal dynamics of interactions is to use point processes as the foundation of network models for community detection. Computational complexity hampers the scalability of such approaches to large sparse networks. To circumvent this challenge, we propose a fast online variational inference algorithm for estimating the latent structure underlying dynamic event arrivals on a network, using continuous-time point process latent network models. We describe this procedure for networks models capturing community structure. This structure can be learned as new events are observed on the network, updating the inferred community assignments. We investigate the theoretical properties of such an inference scheme, and provide regret bounds on the loss function of this procedure. The proposed inference procedure is then thoroughly compared, using both simulation studies and real data, to non-online variants. We demonstrate that online inference can obtain comparable performance, in terms of community recovery, to non-online variants, while realising computational gains. Our proposed inference framework can also be readily modified to incorporate other popular network structures.


title: Thermal effects in spin torque switching of perpendicular magnetic tunnel junctions at cryogenic temperatures

abstract: Temperature plays an important role in spin torque switching of magnetic tunnel junctions causing magnetization fluctuations that decrease the switching voltage but also introduce switching errors. Here we present a systematic study of the temperature dependence of the spin torque switching probability of state-of-the-art perpendicular magnetic tunnel junction nanopillars (40 to 60 nm in diameter) from room temperature down to 4 K, sampling up to a million switching events. The junction temperature at the switching voltage---obtained from the thermally assisted spin torque switching model---saturates at temperatures below about 75 K, showing that junction heating is significant below this temperature and that spin torque switching remains highly stochastic down to 4 K. A model of heat flow in a nanopillar junction shows this effect is associated with the reduced thermal conductivity and heat capacity of the metals in the junction.


title: Dusuk Reynolds Sayili akislarda kanat uzerinde kararsiz aerodinamigin incelenmesi

abstract: The number of studies related to unmanned air vehicles (UAV), micro air vehicles (MAV), wind turbines, birds and insects which operate at low Reynolds (Re) numbers have rapidly increased, recently. In this thesis, unsteady aerodynamics of low Reynolds number flows (Re=25000, 50000 and 75000) on NACA4412 airfoil shaped wing which is mainly used in wind turbines and micro air vehicles and operates at low Reynolds numbers and has aspect ratio of 1 (AR=1) was investigated by means of numerous experimental analyses at various angles of attack changing from 0{\deg} to 45{\deg}. Due to the obtained force measurement results, maximum lift coefficients and stall angles were determined as 1.19, 1.21, 1.24 and 39{\deg}, 38{\deg}, 37{\deg} for Re=25000, 50000 and 75000, respectively. Flow visualization with smoke-wire and oil, velocity distribution experiments were conducted at different planes (z/s= +0.4, +0.2, +0.1, 0, -0.1, -0.2) on the wing. Furthermore, the effects of tip vortices on the wing were determined. The tip vortices were observed at the tip of the wing at low angles of attack and developed through the wake region moved on the wing horizontally with increasing angle of attack and effect of these vortices on the wing increased. It was seen that the tip vortices caused to pressure differences in the wake region and moved to the spanwise direction and this situation led to decrease separation bubble, which is occurred at low Re numbers, and increase stall angle. Additionally, flow condition became steadier due to the tip vortices.


title: GTAdam: Gradient Tracking with Adaptive Momentum for Distributed Online Optimization

abstract: This paper deals with a network of computing agents aiming to solve an online optimization problem in a distributed fashion, i.e., by means of local computation and communication, without any central coordinator. We propose the gradient tracking with adaptive momentum estimation (GTAdam) distributed algorithm, which combines a gradient tracking mechanism with first and second order momentum estimates of the gradient. The algorithm is analyzed in the online setting for strongly convex cost functions with Lipschitz continuous gradients. We provide an upper bound for the dynamic regret given by a term related to the initial conditions and another term related to the temporal variations of the objective functions. Moreover, a linear convergence rate is guaranteed in the static setup. The algorithm is tested on a time-varying classification problem, on a (moving) target localization problem, and in a stochastic optimization setup from image classification. In these numerical experiments from multi-agent learning, GTAdam outperforms state-of-the-art distributed optimization methods.


title: Dynamical Phase Transition of ASEP in the KPZ regime

abstract: We consider the asymmetric simple exclusion process (ASEP) on $\mathbb{Z}$. For continuous densities, ASEP is in local equilibrium for large times, at discontinuities however, one expects to see a dynamical phase transition, i.e. a mixture of different equilibriums. We consider ASEP with deterministic initial data such that at large times, two rarefactions come together at the origin, and the density jumps from $0$ to $1$. Shifting the measure on the KPZ $1/3$ scale, we show that the law of ASEP converges to a mixture of the Dirac measures with only holes resp. only particles. The parameter of that mixture is the probability that the second class particle, which is distributed as the difference of two independent GUEs, stays to the left of the shift. This should be compared with the results of Ferrari and Fontes from 1994 \cite{FF94b}, who obtained a mixture of Bernoulli product measures at discontinuities created by random initial data, with the GUEs replaced by Gaussians.


title: How to evaluate data visualizations across different levels of understanding

abstract: Understanding a visualization is a multi-level process. A reader must extract and extrapolate from numeric facts, understand how those facts apply to both the context of the data and other potential contexts, and draw or evaluate conclusions from the data. A well-designed visualization should support each of these levels of understanding. We diagnose levels of understanding of visualized data by adapting Bloom's taxonomy, a common framework from the education literature. We describe each level of the framework and provide examples for how it can be applied to evaluate the efficacy of data visualizations along six levels of knowledge acquisition - knowledge, comprehension, application, analysis, synthesis, and evaluation. We present three case studies showing that this framework expands on existing methods to comprehensively measure how a visualization design facilitates a viewer's understanding of visualizations. Although Bloom's original taxonomy suggests a strong hierarchical structure for some domains, we found few examples of dependent relationships between performance at different levels for our three case studies. If this level-independence holds across new tested visualizations, the taxonomy could serve to inspire more targeted evaluations of levels of understanding that are relevant to a communication goal.


title: Central points of the double heptagon translation surface are not connection points

abstract: We consider flow directions on the translation surfaces formed from double $(2n+1)$-gons, and give a sufficient condition in terms of a natural gcd algorithm for a direction to be hyperbolic in the sense that it is the fixed direction for some hyperbolic element of the Veech group of the surface. In particular, we give explicit points in the holonomy field of the double heptagon translation surface which are not so-called connection points. Among these are the central points of the heptagons, giving a negative answer to a question by P.Hubert and T.Schmidt.


title: Using Household Grants to Benchmark the Cost Effectiveness of a USAID Workforce Readiness Program

abstract: We use a randomized experiment to compare a workforce training program to cash transfers in Rwanda. Conducted in a sample of poor and underemployed youth, this study measures the impact of the training program not only relative to a control group but relative to the counterfactual of simply disbursing the cost of the program directly to beneficiaries. While the training program was successful in improving a number of core outcomes (productive hours, assets, savings, and subjective well-being), cost-equivalent cash transfers move all these outcomes as well as consumption, income, and wealth. In the head-to-head costing comparison cash proves superior across a number of economic outcomes, while training outperforms cash only in the production of business knowledge. We find little evidence of complementarity between human and physical capital interventions, and no signs of heterogeneity or spillover effects.


title: Electron capture in stars

abstract: Electron captures on nuclei play an essential role for the dynamics of several astrophysical objects. The capture rate can be derived in perturbation theory where allowed nuclear transitions (Gamow-Teller transitions) dominate, except at the higher temperatures achieved in core-collapse supernovae where also forbidden transitions contribute significantly to the rates. There has been decisive progress in recent years in measuring Gamow-Teller (GT) strength distributions using novel experimental techniques based on charge-exchange reactions. These measurements provide not only data for the GT distributions of ground states for many relevant nuclei, but also serve as valuable constraints for nuclear models which are needed to derive the capture rates for the many nuclei, for which no data exist yet. In particular models are needed to evaluate the stellar capture rates at finite temperatures, where the capture can also occur on excited nuclear states. There has also been significant progress in recent years in the modelling of stellar capture rates. This has been made possible by advances in nuclear many-body models as well as in computer soft- and hardware. Specifically to derive reliable capture rates for core-collapse supernovae a dedicated strategy has been developed based on a hierarchy of nuclear models specifically adapted to the abundant nuclei and astrophysically conditions present at the various collapse conditions. This manuscript reviews the experimental and theoretical progress achieved recently in deriving stellar electron capture rates. It also discusses the impact these improved rates have on the various astrophysical objects. (Abridged)


title: Model-Free Design of Control Systems over Wireless Fading Channels

abstract: Wireless control systems replace traditional wired communication with wireless networks to exchange information between actuators, plants and sensors in a control system. The noise in wireless channels renders ideal control policies suboptimal, and their performance is moreover directly dependent on the way in which wireless resources are allocated between control loops. Proper design of the control policy and the resource allocation policy based on both plant states and wireless fading states is then critical to achieve good performance. The resulting problem of co-designing control-aware resource allocation policies and communication-aware controllers, however, is challenging due to its infinite dimensionality, existence of system constraints and need for explicit knowledge of the plants and wireless network models. To overcome those challenges, we rely on constrained reinforcement learning algorithms to propose a model-free approach to the design of wireless control systems. We demonstrate the near optimality of control system performance and stability using near-universal policy parametrizations and present a practical model-free algorithm to learn the co-design policy. Numerical experiments show the strong performance of learned policies over baseline solutions.


title: Dual-probe 1-d hybrid fs/ps rotational CARS for simultaneous single-shot temperature, pressure, and O2/N2 measurements

abstract: We employ dual-probe 1-d fs/ps hybrid rotational coherent anti-Stokes Raman spectroscopy to investigate simultaneous temperature, pressure, and O2/N2 measurements for gas-phase diagnostics. The dual-probe HRCARS technique allows for simultaneous measurements from the time and frequency-domain. A novel approach for measuring pressure, which offers high accuracy (<1%) and precision (0.42%) is presented. The technique is first demonstrated in a chamber for a range of pressures (1-1.5 bar). This technique shows an impressive capability of resolving 1-d pressure gradients arising from a N2 jet impinging on a surface, both in laminar and turbulent conditions. The technique is shown to be capable of resolving single-shot pressure gradients (0.04 bar/mm) originating from kinetic energy conversion to pressure, and resolves characteristic O2/N2 structures from laminar and turbulent mixing.


title: Optimal Streaming of 360 VR Videos with Perfect, Imperfect and Unknown FoV Viewing Probabilities

abstract: In this paper, we investigate wireless streaming of multi-quality tiled 360 virtual reality (VR) videos from a multi-antenna server to multiple single-antenna users in a multi-carrier system. To capture the impact of field-of-view (FoV) prediction, we consider three cases of FoV viewing probability distributions, i.e., perfect, imperfect and unknown FoV viewing probability distributions, and use the average total utility, worst average total utility and worst total utility as the respective performance metrics. We adopt rate splitting with successive decoding for efficient transmission of multiple sets of tiles of different 360 VR videos to their requesting users. In each case, we optimize the encoding rates of the tiles, minimum encoding rates of the FoVs, rates of the common and private messages and transmission beamforming vectors to maximize the total utility. The problems in the three cases are all challenging nonconvex optimization problems. We successfully transform the problem in each case into a difference of convex (DC) programming problem with a differentiable objective function, and obtain a suboptimal solution using concave-convex procedure (CCCP). Finally, numerical results demonstrate the proposed solutions achieve notable gains over existing schemes in all three cases. To the best of our knowledge, this is the first work revealing the impact of FoV prediction and its accuracy on the performance of streaming of multi-quality tiled 360 VR videos.


title: A tale of two Nekrasov's integral equations

abstract: Just 100 years ago, Nekrasov published the widely cited paper \cite{N1}, in which he derived the first of his two integral equations describing steady periodic waves on the free surface of water. We examine how Nekrasov arrived at these equations and his approach to investigating their solutions. In this connection, Nekrasov's life after 1917 is briefly outlined, in particular, how he became a victim of Stalin's terror. Further results concerning Nekrasov's equations and related topicz are surveyed.


title: Group actions of $A_5$ on contractible $2$-complexes

abstract: We prove that every action of $A_5$ on a finite $2$-dimensional contractible complex has a fixed point.


title: A Family of Horizon-penetrating Coordinate Systems for the Schwarzschild Black Hole Geometry with Cauchy Temporal Functions

abstract: We introduce a new family of horizon-penetrating coordinate systems for the Schwarzschild black hole geometry that feature time coordinates, which are specific Cauchy temporal functions, i.e., the level sets of these time coordinates are smooth, asymptotically flat, spacelike Cauchy hypersurfaces. Coordinate systems of this kind are well suited for the study of the temporal evolution of matter and radiation fields in the joined exterior and interior regions of the Schwarzschild black hole geometry, whereas the associated foliations can be employed as initial data sets for the globally hyperbolic development under the Einstein flow. For their construction, we formulate an explicit method that utilizes the geometry of - and structures inherent in - the Penrose diagram of the Schwarzschild black hole geometry, thus relying on the corresponding metrical product structure. As an example, we consider an integrated algebraic sigmoid function as the basis for the determination of such a coordinate system. Finally, we generalize our results to the Reissner-Nordstr\"om black hole geometry up to the Cauchy horizon. The geometric construction procedure presented here can be adapted to yield similar coordinate systems for various other spacetimes with the same metrical product structure.


title: Surrounding the solution of a Linear System of Equations from all sides

abstract: Suppose $A \in \mathbb{R}^{n \times n}$ is invertible and we are looking for the solution of $Ax = b$. Given an initial guess $x_1 \in \mathbb{R}$, we show that by reflecting through hyperplanes generated by the rows of $A$, we can generate an infinite sequence $(x_k)_{k=1}^{\infty}$ such that all elements have the same distance to the solution, i.e. $\|x_k - x\| = \|x_1 - x\|$. If the hyperplanes are chosen at random, averages over the sequence converge and $$ \mathbb{E} \left\| x - \frac{1}{m} \sum_{k=1}^{m}{ x_k} \right\| \leq \frac{1 + \|A\|_F \|A^{-1}\|}{\sqrt{m}} \cdot\|x-x_1\|.$$ The bound does not depend on the dimension of the matrix. This introduces a purely geometric way of attacking the problem: are there fast ways of estimating the location of the center of a sphere from knowing many points on the sphere? Our convergence rate (coinciding with that of the Random Kaczmarz method) comes from averaging, can one do better?


title: Time-resolved study of resonant interatomic Coulombic decay in helium nanodroplets

abstract: When weakly-bound complexes are multiply excited by intense electromagnetic radiation, energy can be exchanged between neighboring atoms through a type of resonant interatomic Coulombic decay (ICD). This decay mechanism due to multiple excitations has been predicted to be relatively slow, typically lasting tens to hundreds of picoseconds. Here, we directly measure the ICD timescale in resonantly excited helium droplets using a high resolution, tunable, extreme ultraviolet free electron laser. Over an extensive range of droplet sizes and laser intensities, we discover the decay to be surprisingly fast, with decay times as fast as 400 femtoseconds, and to only present a weak dependence on the density of the excited states. Using a combination of time dependent density functional theory and ab initio quantum chemistry calculations, we elucidate the mechanisms of this ultrafast decay process where pairs of excited helium atoms in one droplet strongly attract each other and form merging void bubbles which drastically accelerates ICD.


title: Intra-Utterance Similarity Preserving Knowledge Distillation for Audio Tagging

abstract: Knowledge Distillation (KD) is a popular area of research for reducing the size of large models while still maintaining good performance. The outputs of larger teacher models are used to guide the training of smaller student models. Given the repetitive nature of acoustic events, we propose to leverage this information to regulate the KD training for Audio Tagging. This novel KD method, "Intra-Utterance Similarity Preserving KD" (IUSP), shows promising results for the audio tagging task. It is motivated by the previously published KD method: "Similarity Preserving KD" (SP). However, instead of preserving the pairwise similarities between inputs within a mini-batch, our method preserves the pairwise similarities between the frames of a single input utterance. Our proposed KD method, IUSP, shows consistent improvements over SP across student models of different sizes on the DCASE 2019 Task 5 dataset for audio tagging. There is a 27.1% to 122.4% percent increase in improvement of micro AUPRC over the baseline relative to SP's improvement of over the baseline.


title: Classical variational simulation of the Quantum Approximate Optimization Algorithm

abstract: A key open question in quantum computing is whether quantum algorithms can potentially offer a significant advantage over classical algorithms for tasks of practical interest. Understanding the limits of classical computing in simulating quantum systems is an important component of addressing this question. We introduce a method to simulate layered quantum circuits consisting of parametrized gates, an architecture behind many variational quantum algorithms suitable for near-term quantum computers. A neural-network parametrization of the many-qubit wave function is used, focusing on states relevant for the Quantum Approximate Optimization Algorithm (QAOA). For the largest circuits simulated, we reach 54 qubits at 4 QAOA layers, approximately implementing 324 RZZ gates and 216 RX gates without requiring large-scale computational resources. For larger systems, our approach can be used to provide accurate QAOA simulations at previously unexplored parameter values and to benchmark the next generation of experiments in the Noisy Intermediate-Scale Quantum (NISQ) era.


title: Coulomb sum rule for $^4$He and $^{16}$O from coupled-cluster theory

abstract: We demonstrate the capability of coupled-cluster theory to compute the Coulomb sum rule for the $^4$He and $^{16}$O nuclei using interactions from chiral effective field theory. We perform several checks, including a few-body benchmark for $^4$He. We provide an analysis of the center-of-mass contaminations, which we are able to safely remove. We then compare with other theoretical results and experimental data available in the literature, obtaining a fair agreement. This is a first and necessary step towards initiating a program for computing neutrino-nucleus interactions from first principles and supporting the experimental long-baseline neutrino program with a state-of-the-art theory that can reach medium-mass nuclei.


title: A rational Even-IRA algorithm for the solution of T-even polynomial eigenvalue problems

abstract: In this work we present a rational Krylov subspace method for solving real large-scale polynomial eigenvalue problems with T-even (that is, symmetric/skew-symmetric) structure. Our method is based on the Even-IRA algorithm. To preserve the structure, a sparse T-even linearization from the class of block minimal bases pencils is applied. Due to this linearization, the Krylov basis vectors can be computed in a cheap way. A rational decomposition is derived so that our method explicitly allows for changes of the shift during the iteration. This leads to a method that is able to compute parts of the spectrum of a T-even matrix polynomial in a fast and reliable way.


title: Topological Nature of High Temperature Superconductivity

abstract: The key to unraveling the nature of high-temperature superconductivity (HTS) lies in resolving the enigma of the pseudogap state. The pseudogap state in the underdoped region is a distinct thermodynamic phase characterized by nematicity, temperature-quadratic resistive behavior, and magnetoelectric effects. Till present, a general description of the observed universal features of the pseudogap phase and their connection with HTS was lacking. The proposed work constructs a unifying effective field theory capturing all universal characteristics of HTS materials and explaining the observed phase diagram. The pseudogap state is established to be a phase where a charged magnetic monopole condensate confines Cooper pairs to form an oblique version of a superinsulator. The HTS phase diagram is dominated by a tricritical point (TCP) at which the first order transition between a fundamental Cooper pair condensate and a charged magnetic monopole condensate merges with the continuous superconductor-normal metal and superconductor-pseudogap state phase transitions. The universality of the HTS phase diagram reflects a unique topological mechanism of competition between the magnetic monopole condensate, inherent to antiferromagnetic-order-induced Mott insulators and the Cooper pair condensate. The obtained results establish the topological nature of the HTS and provide a platform for devising materials with the enhanced superconducting transition temperature.


title: IrRep: symmetry eigenvalues and irreducible representations of ab initio band structures

abstract: We present IrRep - a Python code that calculates the symmetry eigenvalues of electronic Bloch states in crystalline solids and the irreducible representations under which they transform. As input it receives bandstructures computed with state-of-the-art Density Functional Theory codes such as VASP, Quantum Espresso, or Abinit, as well as any other code that has an interface to Wannier90. Our code is applicable to materials in any of the 230 space groups and double groups preserving time-reversal symmetry with or without spin-orbit coupling included, for primitive or conventional unit cells. This makes IrRep a powerful tool to systematically analyze the connectivity and topological classification of bands, as well as to detect insulators with non-trivial topology, following the Topological Quantum Chemistry formalism: IrRep can generate the input files needed to calculate the (physical) elementary band representations and the symmetry-based indicators using the CheckTopologicalMat routine of the Bilbao Crystallographic Server. It is also particularly suitable for interfaces with other plane-waves based codes, due to its flexible structure.


title: Optimal Load Balanced Demand Distribution under Overload Penalties

abstract: Input to the Load Balanced Demand Distribution (LBDD) consists of the following: (a) a set of service centers; (b) a set of demand nodes and; (c) a cost matrix containing the cost of assignment for each (demand node, service center) pair. In addition, each service center is also associated with a notion of capacity and a penalty which is incurred if it gets overloaded. Given the input, the LBDD problem determines a mapping from the set of n demand vertices to the set of k service centers, n being much larger than k. The objective is to determine a mapping that minimizes the sum of the following two terms: (i) the total cost between demand units and their allotted service centers and, (ii) total penalties incurred. The problem of LBDD has a variety of applications. An instance of the LBDD problem can be reduced to an instance of the min-cost bi-partite matching problem. The best known algorithm for min-cost matching in an unbalanced bipartite graph yields a complexity of O($n^3k$). This paper proposes novel allotment subspace re-adjustment based approach which allows us to characterize the optimality of the mapping without invoking matching or mincost flow. This approach yields an optimal solution with time complexity $O(nk^3 +nk^2 log n)$, and also allows us to efficiently maintain an optimal allotment under insertions and deletions.


title: Synthetic-to-Real Unsupervised Domain Adaptation for Scene Text Detection in the Wild

abstract: Deep learning-based scene text detection can achieve preferable performance, powered with sufficient labeled training data. However, manual labeling is time consuming and laborious. At the extreme, the corresponding annotated data are unavailable. Exploiting synthetic data is a very promising solution except for domain distribution mismatches between synthetic datasets and real datasets. To address the severe domain distribution mismatch, we propose a synthetic-to-real domain adaptation method for scene text detection, which transfers knowledge from synthetic data (source domain) to real data (target domain). In this paper, a text self-training (TST) method and adversarial text instance alignment (ATA) for domain adaptive scene text detection are introduced. ATA helps the network learn domain-invariant features by training a domain classifier in an adversarial manner. TST diminishes the adverse effects of false positives~(FPs) and false negatives~(FNs) from inaccurate pseudo-labels. Two components have positive effects on improving the performance of scene text detectors when adapting from synthetic-to-real scenes. We evaluate the proposed method by transferring from SynthText, VISD to ICDAR2015, ICDAR2013. The results demonstrate the effectiveness of the proposed method with up to 10% improvement, which has important exploration significance for domain adaptive scene text detection. Code is available at https://github.com/weijiawu/SyntoReal_STD


title: Nonminimal Lorentz-Violating Extensions of Gauge Field Theories

abstract: A general method is presented to build all gauge-invariant terms in gauge field theories, including quantum electrodynamics and quantum chromodynamics. It is applied to two experiments, light-by-light scattering and deep inelastic scattering, to extract first bounds on certain nonminimal coefficients for Lorentz violation.


title: Design of Gallium Nitride Resonant Cavity Light Emitting Diodes on Si Substrates

abstract: A GaN resonant cavity light emitting diode was built on a GaN-AlN distributed Bragg reflector grown on a silicon substrate. The electroluminescence output increased by 2.5 times for a GaN diode coupled to a properly designed resonant cavity. Theoretical calculations showed that this enhancement could increase up to four times for transmission through a sem-transparent metal contact design, up to eight times for a flip-chip design


title: HyperBench: A Benchmark and Tool for Hypergraphs and Empirical Findings

abstract: To cope with the intractability of answering Conjunctive Queries (CQs) and solving Constraint Satisfaction Problems (CSPs), several notions of hypergraph decompositions have been proposed -- giving rise to different notions of width, noticeably, plain, generalized, and fractional hypertree width (hw, ghw, and fhw). Given the increasing interest in using such decomposition methods in practice, a publicly accessible repository of decomposition software, as well as a large set of benchmarks, and a web-accessible workbench for inserting, analyzing, and retrieving hypergraphs are called for. We address this need by providing (i) concrete implementations of hypergraph decompositions (including new practical algorithms), (ii) a new, comprehensive benchmark of hypergraphs stemming from disparate CQ and CSP collections, and (iii) HyperBench, our new web-inter\-face for accessing the benchmark and the results of our analyses. In addition, we describe a number of actual experiments we carried out with this new infrastructure.


title: Exterior bundles in diffeology

abstract: We explore several notions of $k$-form at a point in a diffeological space, construct bundles of such $k$-forms, and compare sections of these bundles to differential forms. As they are defined locally, our $k$-forms can contain more information than the values of differential forms contain, and we illustrate this with many examples. To organize our work, we develop the basic theory of diffeological vector pseudo-bundles, including a detailed understanding of their limits and colimits, as well as a variety of fibrewise operations such as products, direct sums, tensor products, exterior powers and dual bundles.


title: Unified Framework for $B$-Anomalies, Muon $g-2$, and Neutrino Masses

abstract: We present a model of radiative neutrino masses which also resolves anomalies reported in $B$-meson decays, $R_{D^{(\star)}}$ and $R_{K^{(\star)}}$, as well as in muon $g-2$ measurement, $\Delta a_\mu$. Neutrino masses arise in the model through loop diagrams involving TeV-scale leptoquark (LQ) scalars $R_2$ and $S_3$. Fits to neutrino oscillation parameters are obtained satisfying all flavor constraints which also explain the anomalies in $R_{D^{(\star)}}$, $R_{K^{(\star)}}$ and $\Delta a_\mu$ within $1\, \sigma$. An isospin-3/2 Higgs quadruplet plays a crucial role in generating neutrino masses; we point out that the doubly-charged scalar contained therein can be produced in the decays of the $S_3$ LQ, which enhances its reach to 1.1 (6.2) TeV at $\sqrt s=14$ TeV high-luminosity LHC ($\sqrt s=100$ TeV FCC-hh). We also present flavor-dependent upper limits on the Yukawa couplings of the LQs to the first two family fermions, arising from non-resonant dilepton ($pp \rightarrow \ell^+ \ell^-$) processes mediated by $t$-channel LQ exchange, which for 1 TeV LQ mass, are found to be in the range $(0.15 - 0.36)$. These limits preclude any explanation of $R_{D^{(\star)}}$ through LQ-mediated $B$-meson decays involving $\nu_e$ or $\nu_\mu$ in the final state. We also find that the same Yukawa couplings responsible for the chirally-enhanced contribution to $\Delta a_\mu$ give rise to new contributions to the SM Higgs decays to muon and tau pairs, with the modifications to the corresponding branching ratios being at (2-6)% level, which could be tested at future hadron colliders, such as HL-LHC and FCC-hh.


title: Reading In-Between the Lines: An Analysis of Dissenter

abstract: Efforts by content creators and social networks to enforce legal and policy-based norms, e.g. blocking hate speech and users, has driven the rise of unrestricted communication platforms. One such recent effort is Dissenter, a browser and web application that provides a conversational overlay for any web page. These conversations hide in plain sight - users of Dissenter can see and participate in this conversation, whereas visitors using other browsers are oblivious to their existence. Further, the website and content owners have no power over the conversation as it resides in an overlay outside their control. In this work, we obtain a history of Dissenter comments, users, and the websites being discussed, from the initial release of Dissenter in Feb. 2019 through Apr. 2020 (14 months). Our corpus consists of approximately 1.68M comments made by 101k users commenting on 588k distinct URLs. We first analyze macro characteristics of the network, including the user-base, comment distribution, and growth. We then use toxicity dictionaries, Perspective API, and a Natural Language Processing model to understand the nature of the comments and measure the propensity of particular websites and content to elicit hateful and offensive Dissenter comments. Using curated rankings of media bias, we examine the conditional probability of hateful comments given left and right-leaning content. Finally, we study Dissenter as a social network, and identify a core group of users with high comment toxicity.


title: Short duration gamma-ray bursts and their outflows in light of GW170817

abstract: The detection of GW170817, it's extensive multi-wavelength follow-up campaign, and the large amount of theoretical development and interpretation that followed, have resulted in a significant step forward in the understanding of the binary neutron star merger phenomenon as a whole. One of its aspects is seeing the merger as a progenitor of short gamma-ray bursts (SGRB), which will be the subject of this review. On the one hand, GW170817 observations have confirmed some theoretical expectations, exemplified by the confirmation that binary neutron star mergers are the progenitors of SGRBs. In addition, the multimessenger nature of GW170817 has allowed for gathering of unprecedented data, such as the trigger time of the merger, the delay with which the gamma-ray photons were detected, and the brightening afterglow of an off-axis event. All together, the incomparable richness of the data from GW170817 has allowed us to paint a fairly detailed picture of at least one SGRB. I will detail what we learned, what new questions have arisen, and the perspectives for answering them when a sample of GW170817-comparable events have been studied.


title: Quantitative inverse theory of Gowers uniformity norms

abstract: (This text is a survey written for the Bourbaki seminar on the work of F. Manners.) Gowers uniformity norms are the central objects of higher order Fourier analysis, one of the cornerstones of additive combinatorics, and play an important role in both Gowers' proof of Szemer\'{e}di's theorem and the Green-Tao theorem. The inverse theorem states that if a function has a large uniformity norm, which is a robust combinatorial measure of structure, then it must correlate with a nilsequence, which is a highly structured algebraic object. This was proved in a qualitative sense by Green, Tao, and Ziegler, but with a proof that was incapable of providing reasonable bounds. In 2018 Manners achieved a breakthrough by giving a new proof of the inverse theorem. Not only does this new proof contain a wealth of new insights but it also, for the first time, provides quantitative bounds, that are at worst only doubly exponential. This talk will give a high-level overview of what the inverse theorem says, why it is important, and the new proof of Manners.


title: On the impact of baryons on the halo mass function, bias, and cluster cosmology

abstract: Luminous matter produces very energetic events, such as active galactic nuclei and supernova explosions, that significantly affect the internal regions of galaxy clusters. Although the current uncertainty in the effect of baryonic physics on cluster statistics is subdominant as compared to other systematics, the picture is likely to change soon as the amount of high-quality data is growing fast, urging the community to keep theoretical systematic uncertainties below the ever-growing statistical precision. In this paper, we study the effect of baryons on galaxy clusters, and their impact on the cosmological applications of clusters, using the Magneticum suite of cosmological hydrodynamical simulations. We show that the impact of baryons on the halo mass function can be recast in terms on a variation of the mass of the halos simulated with pure N-body, when baryonic effects are included. The halo mass function and halo bias are only indirectly affected. Finally, we demonstrate that neglecting baryonic effects on halos mass function and bias would significantly alter the inference of cosmological parameters from high-sensitivity next-generations surveys of galaxy clusters.


title: HiFiSinger: Towards High-Fidelity Neural Singing Voice Synthesis

abstract: High-fidelity singing voices usually require higher sampling rate (e.g., 48kHz) to convey expression and emotion. However, higher sampling rate causes the wider frequency band and longer waveform sequences and throws challenges for singing voice synthesis (SVS) in both frequency and time domains. Conventional SVS systems that adopt small sampling rate cannot well address the above challenges. In this paper, we develop HiFiSinger, an SVS system towards high-fidelity singing voice. HiFiSinger consists of a FastSpeech based acoustic model and a Parallel WaveGAN based vocoder to ensure fast training and inference and also high voice quality. To tackle the difficulty of singing modeling caused by high sampling rate (wider frequency band and longer waveform), we introduce multi-scale adversarial training in both the acoustic model and vocoder to improve singing modeling. Specifically, 1) To handle the larger range of frequencies caused by higher sampling rate, we propose a novel sub-frequency GAN (SF-GAN) on mel-spectrogram generation, which splits the full 80-dimensional mel-frequency into multiple sub-bands and models each sub-band with a separate discriminator. 2) To model longer waveform sequences caused by higher sampling rate, we propose a multi-length GAN (ML-GAN) for waveform generation to model different lengths of waveform sequences with separate discriminators. 3) We also introduce several additional designs and findings in HiFiSinger that are crucial for high-fidelity voices, such as adding F0 (pitch) and V/UV (voiced/unvoiced flag) as acoustic features, choosing an appropriate window/hop size for mel-spectrogram, and increasing the receptive field in vocoder for long vowel modeling. Experiment results show that HiFiSinger synthesizes high-fidelity singing voices with much higher quality: 0.32/0.44 MOS gain over 48kHz/24kHz baseline and 0.83 MOS gain over previous SVS systems.


title: A small-step approach to multi-trace checking against interactions

abstract: Interaction models describe the exchange of messages between the different components of distributed systems. We have previously defined a small-step operational semantics for interaction models. The paper extends this work by presenting an approach for checking the validity of multi-traces against interaction models. A multi-trace is a collection of traces (sequences of emissions and receptions), each representing a local view of the same global execution of the distributed system. We have formally proven our approach, studied its complexity, and implemented it in a prototype tool. Finally, we discuss some observability issues when testing distributed systems via the analysis of multi-traces.


title: Reconstructing 2D spatial modes for classical and quantum light

abstract: We propose a method for finding 2D spatial modes of thermal field through a direct measurement of the field intensity and an offline analysis of its spatial fluctuations. Using this method, in a simple and efficient way we reconstruct the modes of a multimode fiber and the spatial Schmidt modes of squeezed vacuum generated via high-gain parametric down conversion. The reconstructed shapes agree with the theoretical results.


title: Two-phase interfacial structure of bubbly-to-slug transition flows in a 12.7 mm ID vertical tube

abstract: This experimental study focuses on the characteristics of air-water two-phase interfacial structure. Interfacial parameters including void fraction, interfacial area concentration, and bubble interfacial velocity are measured using four-sensor electrical conductivity probe on a 12.7 mm ID vertical tube. The tube size is approximately equal to the maximum distorted bubble size. Therefore, the bubbly-to-slug transition characteristics can be different from in other sizes of tubes. Comparing with previous studies, this study provides an experimental database with a wide range on the bubbly-to-slug transition flows, with 4 different superficial liquid velocities (0.3, 0.5, 1.0, and 2.0 m/s) and void fraction ranging from 0.07 to 0.66. Experimental results show that the wall-peak void distribution does not appear in a small diameter tube under the bubbly-to-slug flow transition flow. The distribution is related with both void fraction and the relative bubble size to the tube size. In this sense, a new correlation of distribution parameter in the Drift Flux model is proposed based on the previous studies by Ishii and Hibiki et al. This experimental study can be a good reference for the model development of flow regime transition and the Interfacial Area Transport Equation.


title: Topological phase transition in coupled rock-paper-scissor cycles

abstract: A hallmark of topological phases is the occurrence of topologically protected modes at the system`s boundary. Here we find topological phases in the antisymmetric Lotka-Volterra equation (ALVE). The ALVE is a nonlinear dynamical system and describes, e.g., the evolutionary dynamics of a rock-paper-scissors cycle. On a one-dimensional chain of rock-paper-scissor cycles, topological phases become manifest as robust polarization states. At the transition point between left and right polarization, solitonic waves are observed. This topological phase transition lies in symmetry class $D$ within the "ten-fold way" classification as also realized by 1D topological superconductors.


title: Role of the conserved charges in the chiral symmetry restoration phase transition

abstract: The effect of conserved baryon, isospin and strangeness charges on the behavior of phase transitions in dense matter is studied. Baryonic matter is described within the three-flavor Polyakov$-$Nambu$-$Jona-Lasinio model and several charge fractions $Y_Q$ are considered. The role of the vector interaction, which can be important to describe dense systems, is discussed. Special attention is given to the case with charge fraction $Y_Q=0.4$, due to its importance in heavy-ion collisions and core-collapse supernova matter. It is shown that the possible formation of chiral-symmetric quark matter in the laboratory will be favored in asymmetric matter. Besides, the inclusion of the vector interaction reinforces the formation of quark matter at lower densities.


title: Limited View Tomographic Reconstruction Using a Deep Recurrent Framework with Residual Dense Spatial-Channel Attention Network and Sinogram Consistency

abstract: Limited view tomographic reconstruction aims to reconstruct a tomographic image from a limited number of sinogram or projection views arising from sparse view or limited angle acquisitions that reduce radiation dose or shorten scanning time. However, such a reconstruction suffers from high noise and severe artifacts due to the incompleteness of sinogram. To derive quality reconstruction, previous state-of-the-art methods use UNet-like neural architectures to directly predict the full view reconstruction from limited view data; but these methods leave the deep network architecture issue largely intact and cannot guarantee the consistency between the sinogram of the reconstructed image and the acquired sinogram, leading to a non-ideal reconstruction. In this work, we propose a novel recurrent reconstruction framework that stacks the same block multiple times. The recurrent block consists of a custom-designed residual dense spatial-channel attention network. Further, we develop a sinogram consistency layer interleaved in our recurrent framework in order to ensure that the sampled sinogram is consistent with the sinogram of the intermediate outputs of the recurrent blocks. We evaluate our methods on two datasets. Our experimental results on AAPM Low Dose CT Grand Challenge datasets demonstrate that our algorithm achieves a consistent and significant improvement over the existing state-of-the-art neural methods on both limited angle reconstruction (over 5dB better in terms of PSNR) and sparse view reconstruction (about 4dB better in term of PSNR). In addition, our experimental results on Deep Lesion datasets demonstrate that our method is able to generate high-quality reconstruction for 8 major lesion types.


title: Quantum Long Short-Term Memory

abstract: Long short-term memory (LSTM) is a kind of recurrent neural networks (RNN) for sequence and temporal dependency data modeling and its effectiveness has been extensively established. In this work, we propose a hybrid quantum-classical model of LSTM, which we dub QLSTM. We demonstrate that the proposed model successfully learns several kinds of temporal data. In particular, we show that for certain testing cases, this quantum version of LSTM converges faster, or equivalently, reaches a better accuracy, than its classical counterpart. Due to the variational nature of our approach, the requirements on qubit counts and circuit depth are eased, and our work thus paves the way toward implementing machine learning algorithms for sequence modeling on noisy intermediate-scale quantum (NISQ) devices.


title: Device-independent quantum key distribution from generalized CHSH inequalities

abstract: Device-independent quantum key distribution aims at providing security guarantees even when using largely uncharacterised devices. In the simplest scenario, these guarantees are derived from the CHSH score, which is a simple linear combination of four correlation functions. We here derive a security proof from a generalisation of the CHSH score, which effectively takes into account the individual values of two correlation functions. We show that this additional information, which is anyway available in practice, allows one to get higher key rates than with the CHSH score. We discuss the potential advantage of this technique for realistic photonic implementations of device-independent quantum key distribution.


title: Data-First Visualization Design Studies

abstract: We introduce the notion of a data-first design study which is triggered by the acquisition of real-world data instead of specific stakeholder analysis questions. We propose an adaptation of the design study methodology framework to provide practical guidance and to aid transferability to other data-first design processes. We discuss opportunities and risks by reflecting on two of our own data-first design studies. We review 64 previous design studies and identify 16 of them as edge cases with characteristics that may indicate a data-first design process in action.


title: Computational Analysis of Deformable Manifolds: from Geometric Modelling to Deep Learning

abstract: Leo Tolstoy opened his monumental novel Anna Karenina with the now famous words: Happy families are all alike; every unhappy family is unhappy in its own way A similar notion also applies to mathematical spaces: Every flat space is alike; every unflat space is unflat in its own way. However, rather than being a source of unhappiness, we will show that the diversity of non-flat spaces provides a rich area of study. The genesis of the so-called big data era and the proliferation of social and scientific databases of increasing size has led to a need for algorithms that can efficiently process, analyze and, even generate high dimensional data. However, the curse of dimensionality leads to the fact that many classical approaches do not scale well with respect to the size of these problems. One technique to avoid some of these ill-effects is to exploit the geometric structure of coherent data. In this thesis, we will explore geometric methods for shape processing and data analysis. More specifically, we will study techniques for representing manifolds and signals supported on them through a variety of mathematical tools including, but not limited to, computational differential geometry, variational PDE modeling, and deep learning. First, we will explore non-isometric shape matching through variational modeling. Next, we will use ideas from parallel transport on manifolds to generalize convolution and convolutional neural networks to deformable manifolds. Finally, we conclude by proposing a novel auto-regressive model for capturing the intrinsic geometry and topology of data. Throughout this work, we will use the idea of computing correspondences as a though-line to both motivate our work and analyze our results.


title: Singular solutions to Yamabe-type systems with prescribed asymptotics

abstract: Our primary purpose is to study a class of strongly coupled nonlinear elliptic systems with critical growth in a compact Riemannian manifold with constant scalar curvature. Using a gluing technique and perturbation arguments, we show the existence of singular solutions asymptotic to a Fowler-type solution near the isolated singularity.


title: Black-to-White Hole Scenario: Foundation and Evaporation

abstract: A theory of quantum gravity is expected to change profoundly our understanding of black holes. Quantum theory has already shown, as a first approximation, that a black hole slowly evaporates. Non-perturbative quantum gravity also predicts that such an object could metamorphose into a white hole, its time-reverse. In this thesis, we study the foundations of such a scenario and propose a mathematical model that includes the phenomenon of evaporation.


title: The fully developed remnant of a neutrino-driven supernova: Evolution of ejecta structure and asymmetries in SNR Cassiopeia A

abstract: Abridged. We aim at exploring to which extent the remnant keeps memory of the asymmetries that develop stochastically in the neutrino-heating layer due to hydrodynamic instabilities (e.g., convective overturn and the standing accretion shock instability) during the first second after core bounce. We coupled a 3D HD model of a neutrino-driven SN explosion with 3D MHD/HD simulations of the remnant formation. The simulations cover 2000 years of expansion and include all physical processes relevant to describe the complexities in the SN evolution and the subsequent interaction of the stellar debris with the wind of the progenitor star. The interaction of large-scale asymmetries left from the earliest phases of the explosion with the reverse shock produces, at the age of $\approx 350$~years, an ejecta structure and a remnant morphology which are remarkably similar to those observed in Cas A. Small-scale structures in the large-scale Fe-rich plumes created during the initial stages of the SN, combined with HD instabilities that develop after the passage of the reverse shock, naturally produce a pattern of ring- and crown-like structures of shocked ejecta. The consequence is a spatial inversion of the ejecta layers with Si-rich ejecta being physically interior to Fe-rich ejecta. The full-fledged remnant shows voids and cavities in the innermost unshocked ejecta resulting from the expansion of Fe-rich plumes and their inflation due to the decay of radioactive species. The asymmetric distributions of $^{44}$Ti and $^{56}$Fe and their abundance ratio are both compatible with those inferred from high-energy observations of Chandra and NuSTAR. The main asymmetries observed in the ejecta distribution of Cas A can be explained by the interaction of the reverse shock with the large-scale asymmetries that developed from stochastic processes that originate during the first seconds of the SN blast.


title: Private Weighted Random Walk Stochastic Gradient Descent

abstract: We consider a decentralized learning setting in which data is distributed over nodes in a graph. The goal is to learn a global model on the distributed data without involving any central entity that needs to be trusted. While gossip-based stochastic gradient descent (SGD) can be used to achieve this learning objective, it incurs high communication and computation costs, since it has to wait for all the local models at all the nodes to converge. To speed up the convergence, we propose instead to study random walk based SGD in which a global model is updated based on a random walk on the graph. We propose two algorithms based on two types of random walks that achieve, in a decentralized way, uniform sampling and importance sampling of the data. We provide a non-asymptotic analysis on the rate of convergence, taking into account the constants related to the data and the graph. Our numerical results show that the weighted random walk based algorithm has a better performance for high-variance data. Moreover, we propose a privacy-preserving random walk algorithm that achieves local differential privacy based on a Gamma noise mechanism that we propose. We also give numerical results on the convergence of this algorithm and show that it outperforms additive Laplace-based privacy mechanisms.


title: Action and Perception as Divergence Minimization

abstract: To learn directed behaviors in complex environments, intelligent agents need to optimize objective functions. Various objectives are known for designing artificial agents, including task rewards and intrinsic motivation. However, it is unclear how the known objectives relate to each other, which objectives remain yet to be discovered, and which objectives better describe the behavior of humans. We introduce the Action Perception Divergence (APD), an approach for categorizing the space of possible objective functions for embodied agents. We show a spectrum that reaches from narrow to general objectives. While the narrow objectives correspond to domain-specific rewards as typical in reinforcement learning, the general objectives maximize information with the environment through latent variable models of input sequences. Intuitively, these agents use perception to align their beliefs with the world and use actions to align the world with their beliefs. They infer representations that are informative of past inputs, explore future inputs that are informative of their representations, and select actions or skills that maximally influence future inputs. This explains a wide range of unsupervised objectives from a single principle, including representation learning, information gain, empowerment, and skill discovery. Our findings suggest leveraging powerful world models for unsupervised exploration as a path toward highly adaptive agents that seek out large niches in their environments, rendering task rewards optional.


title: Accessing the high-$\ell$ frontier under the Reduced Shear Approximation with $k$-cut Cosmic Shear

abstract: The precision of Stage IV cosmic shear surveys will enable us to probe smaller physical scales than ever before, however, model uncertainties from baryonic physics and non-linear structure formation will become a significant concern. The $k$-cut method -- applying a redshift-dependent $\ell$-cut after making the Bernardeau-Nishimichi-Taruya transform -- can reduce sensitivity to baryonic physics; allowing Stage IV surveys to include information from increasingly higher $\ell$-modes. Here we address the question of whether it can also mitigate the impact of making the reduced shear approximation; which is also important in the high-$\kappa$, small-scale regime. The standard procedure for relaxing this approximation requires the repeated evaluation of the convergence bispectrum, and consequently can be prohibitively computationally expensive when included in Monte Carlo analyses. We find that the $k$-cut cosmic shear procedure suppresses the $w_0w_a$CDM cosmological parameter biases expected from the reduced shear approximation for Stage IV experiments, when $\ell$-modes up to $5000$ are probed. The maximum cut required for biases from the reduced shear approximation to be below the threshold of significance is at $k = 5.37 \, h{\rm Mpc}^{-1}$. With this cut, the predicted $1\sigma$ constraints increase, relative to the case where the correction is directly computed, by less than $10\%$ for all parameters. This represents a significant improvement in constraints compared to the more conservative case where only $\ell$-modes up to 1500 are probed, and no $k$-cut is used. We also repeat this analysis for a hypothetical, comparable kinematic weak lensing survey. The key parts of code used for this analysis are made publicly available.


title: Optimal approximants and orthogonal polynomials in several variables II: families of polynomials in the unit ball

abstract: We obtain closed expressions for weighted orthogonal polynomials and optimal approximants associated with the function $f(z)=1-\frac{1}{\sqrt{2}}(z_1+z_2)$ and a scale of Hilbert function spaces in the unit $2$-ball having reproducing kernel $(1-\langle z,w\rangle)^{-\gamma}$, $\gamma>0$. Our arguments are elementary but do not rely on reduction to the one-dimensional case.


title: Note on stability and microcausality in Lorentz violating antisymmetric tensor field

abstract: The fundamental issues of microcausality and energy positivity conditions are important to be investigated in the context of spontaneous Lorentz violating theories. We check the microcausality and energy positivity conditions for a free rank-2 antisymmetric tensor field with spontaneous Lorentz violating term and a classically equivalent vector theory from dispersion relations and propagators, and find that the two theories lead to different conditions. While the antisymmetric tensor theory satisfies energy positivity condition and can satisfy microcausality condition for a particular choice of vacuum value, the vector theory violates both energy positivity and microcausality conditions.


title: Totally p-adic Numbers of Degree 3

abstract: The height of an algebraic number $\alpha$ is a measure of how arithmetically complicated $\alpha$ is. We say $\alpha$ is totally $p$-adic if the minimal polynomial of $\alpha$ splits completely over the field $\mathbb{Q}_p$ of $p$-adic numbers. In this paper, we investigate what can be said about the smallest nonzero height of a degree $3$ totally $p$-adic number. In particular, we provide an algorithm to determine, given a prime $p$, the smallest height of a degree $3$ totally $p$-adic number.


title: On-demand quantum storage of photonic qubits in an on-chip waveguide

abstract: Photonic quantum memory is the core element in quantum information processing (QIP). For the scalable and convenient practical applications, great efforts have been devoted to the integrated quantum memory based on various waveguides fabricated in solids. However, on-demand storage of qubits, which is an essential requirement for QIP, is still challenging to be implemented using such integrated quantum memory. Here we report the on-demand storage of time-bin qubits in an on-chip waveguide memory on the surface of a $^{151}$Eu$^{3+}$:Y$_2$SiO$_5$ crystal, utilizing the Stark modulated atomic frequency comb protocol. A qubit storage fidelity of $99.3\%\pm0.2\%$ is obtained with a input of 0.5 photons per pulse, far beyond the highest fidelity achievable using the classical measure-and-prepare strategy. The developed integrated quantum memory with the on-demand retrieval capability, represents an important step towards practical applications of integrated quantum nodes in quantum networks.


title: A Wholistic View of Continual Learning with Deep Neural Networks: Forgotten Lessons and the Bridge to Active and Open World Learning

abstract: Current deep learning methods are regarded as favorable if they empirically perform well on dedicated test sets. This mentality is seamlessly reflected in the resurfacing area of continual learning, where consecutively arriving data is investigated. The core challenge is framed as protecting previously acquired representations from being catastrophically forgotten. However, comparison of individual methods is nevertheless performed in isolation from the real world by monitoring accumulated benchmark test set performance. The closed world assumption remains predominant, i.e. models are evaluated on data that is guaranteed to originate from the same distribution as used for training. This poses a massive challenge as neural networks are well known to provide overconfident false predictions on unknown and corrupted instances. In this work we critically survey the literature and argue that notable lessons from open set recognition, identifying unknown examples outside of the observed set, and the adjacent field of active learning, querying data to maximize the expected performance gain, are frequently overlooked in the deep learning era. Hence, we propose a consolidated view to bridge continual learning, active learning and open set recognition in deep neural networks. Finally, the established synergies are supported empirically, showing joint improvement in alleviating catastrophic forgetting, querying data, selecting task orders, while exhibiting robust open world application.


title: Ramifications of Approximate Posterior Inference for Bayesian Deep Learning in Adversarial and Out-of-Distribution Settings

abstract: Deep neural networks have been successful in diverse discriminative classification tasks, although, they are poorly calibrated often assigning high probability to misclassified predictions. Potential consequences could lead to trustworthiness and accountability of the models when deployed in real applications, where predictions are evaluated based on their confidence scores. Existing solutions suggest the benefits attained by combining deep neural networks and Bayesian inference to quantify uncertainty over the models' predictions for ambiguous datapoints. In this work we propose to validate and test the efficacy of likelihood based models in the task of out of distribution detection (OoD). Across different datasets and metrics we show that Bayesian deep learning models on certain occasions marginally outperform conventional neural networks and in the event of minimal overlap between in/out distribution classes, even the best models exhibit a reduction in AUC scores in detecting OoD data. Preliminary investigations indicate the potential inherent role of bias due to choices of initialisation, architecture or activation functions. We hypothesise that the sensitivity of neural networks to unseen inputs could be a multi-factor phenomenon arising from the different architectural design choices often amplified by the curse of dimensionality. Furthermore, we perform a study to find the effect of the adversarial noise resistance methods on in and out-of-distribution performance, as well as, also investigate adversarial noise robustness of Bayesian deep learners.


title: Globally-centered autocovariances in MCMC

abstract: Autocovariances are a fundamental quantity of interest in Markov chain Monte Carlo (MCMC) simulations with autocorrelation function (ACF) plots being an integral visualization tool for performance assessment. Unfortunately, for slow-mixing Markov chains, the empirical autocovariance can highly underestimate the truth. For multiple-chain MCMC sampling, we propose a globally-centered estimator of the autocovariance function (G-ACvF) that exhibits significant theoretical and empirical improvements. We show that the bias of the G-ACvF estimator is smaller than the bias of the current state-of-the-art. The impact of this improved estimator is evident in three critical output analysis applications: (1) ACF plots, (2) estimates of the Monte Carlo asymptotic covariance matrix, and (3) estimates of the effective sample size. Under weak conditions, we establish strong consistency of our improved asymptotic covariance estimator, and obtain its large-sample bias and variance. The performance of the new estimators is demonstrated through various examples.


title: Inaccuracy measures for concomitants of GOS in Morgenstern family

abstract: In this paper, we obtain a measure of inaccuracy between rth concomitant of generalized order statistic and the parent random variable in Morgenstern family. Applications of this result are given for concomitants of order statistics and record values. We also study some results of cumulative past inaccuracy (CPI) between the distribution function of rth concomitant of order statistic (record value) and the distribution function of parent random variable. Finally, we discuss on a problem of estimating the CPI by means of the empirical CPI in concomitants of generalized order statistics.


title: Local marker for interacting topological insulators

abstract: Topological states of matter were first introduced for non-interacting fermions on an infinite uniform lattice. Since then, substantial effort has been made to generalize these concepts to more complex settings. Recently, local markers have been developed that can describe the topological state of systems without translational symmetry and well-defined gap. However, no local marker for interacting matter has been proposed yet that is capable of directly addressing an interacting system. Here we suggest such a many-body local marker based on the single-particle Green's function. Using this marker we identify topological transitions in finite lattices of a Chern insulator with Anderson disorder and Hubbard interactions. Importantly, our proposal can be straightforwardly generalised to non-equilibrium systems.


title: Bipartite Matching in Nearly-linear Time on Moderately Dense Graphs

abstract: We present an $\tilde O(m+n^{1.5})$-time randomized algorithm for maximum cardinality bipartite matching and related problems (e.g. transshipment, negative-weight shortest paths, and optimal transport) on $m$-edge, $n$-node graphs. For maximum cardinality bipartite matching on moderately dense graphs, i.e. $m = \Omega(n^{1.5})$, our algorithm runs in time nearly linear in the input size and constitutes the first improvement over the classic $O(m\sqrt{n})$-time [Dinic 1970; Hopcroft-Karp 1971; Karzanov 1973] and $\tilde O(n^\omega)$-time algorithms [Ibarra-Moran 1981; Mucha-Sankowski 2004] (where currently $\omega\approx 2.373$). On sparser graphs, i.e. when $m = n^{9/8 + \delta}$ for any constant $\delta>0$, our result improves upon the recent advances of [Madry 2013] and [Liu-Sidford 2020b, 2020a] which achieve an $\tilde O(m^{4/3+o(1)})$ runtime. We obtain these results by combining and advancing recent lines of research in interior point methods (IPMs) and dynamic graph algorithms. First, we simplify and improve the IPM of [v.d.Brand-Lee-Sidford-Song 2020], providing a general primal-dual IPM framework and new sampling-based techniques for handling infeasibility induced by approximate linear system solvers. Second, we provide a simple sublinear-time algorithm for detecting and sampling high-energy edges in electric flows on expanders and show that when combined with recent advances in dynamic expander decompositions, this yields efficient data structures for maintaining the iterates of both [v.d.Brand et al.] and our new IPMs. Combining this general machinery yields a simpler $\tilde O(n \sqrt{m})$ time algorithm for matching based on the logarithmic barrier function, and our state-of-the-art $\tilde O(m+n^{1.5})$ time algorithm for matching based on the [Lee-Sidford 2014] barrier (as regularized in [v.d.Brand et al.]).


title: Sparse Meta Networks for Sequential Adaptation and its Application to Adaptive Language Modelling

abstract: Training a deep neural network requires a large amount of single-task data and involves a long time-consuming optimization phase. This is not scalable to complex, realistic environments with new unexpected changes. Humans can perform fast incremental learning on the fly and memory systems in the brain play a critical role. We introduce Sparse Meta Networks -- a meta-learning approach to learn online sequential adaptation algorithms for deep neural networks, by using deep neural networks. We augment a deep neural network with a layer-specific fast-weight memory. The fast-weights are generated sparsely at each time step and accumulated incrementally through time providing a useful inductive bias for online continual adaptation. We demonstrate strong performance on a variety of sequential adaptation scenarios, from a simple online reinforcement learning to a large scale adaptive language modelling.


title: Intriguing relaxor ferroelectricity in a polar CT crystal at the neutral-ionic interface

abstract: We investigated the mixed-stack charge-transfer crystal, N,N'-dimethylphenazine-TCNQ (M$_2$P-TCNQ), which is polar at room temperature and just at the neutral-to-ionic interface (ionicity $\rho \approx 0.5$). We detect the typical dielectric signature of a relaxor ferroelectric and an asymmetric positive-up-negative-down behavior. While relaxor ferroelectricity is usually ascribed to disorder in the crystal, we find no evidence for structural disorder in the investigated crystals. To elucidate the origin of M$_2$P-TCNQ's dielectric properties we perform parallel structural and spectroscopic measurements, associated with theoretical modeling and quantum-mechanical calculations. Our combined effort points to a highly polarizable electronic system that is strongly coupled to lattice vibrations. The found indications for polarization reversal imply flipping of the bent conformation of the M$_2$P molecule with an associated energy barrier of a few tens of an eV, broadly consistent with an Arrhenius fit of the dielectric relaxation times. While the polarization is mostly of electronic origin, its possible reversal implies slow collective motions that are affected by solid-state intermolecular interactions.