title: Session: A Model for End-To-End Encrypted Conversations With Minimal Metadata Leakage

abstract: Session is an open-source, public-key-based secure messaging application which uses a set of decentralised storage servers and an onion routing protocol to send end-to-end encrypted messages with minimal exposure of user metadata. It does this while also providing common features of mainstream messaging applications.


title: A method for mapping the aliphatic hydrocarbon content of interstellar dust towards the Galactic Centre

abstract: In the interstellar medium, the cosmic elemental carbon abundance includes the total carbon in both gas and solid phases. The aim of the study was to trial a new method for measuring the amount and distribution of aliphatic carbon within interstellar dust over wide fields of view of our Galaxy. This method is based on measurement of the 3.4 $\mu$m absorption feature from aliphatic carbonaceous matter. This can readily be achieved for single sources using IR spectrometers. However, making such measurements over wide fields requires an imaging IR camera, equipped with narrow-band filters that are able to sample the spectrum. While this cannot produce as good a determination of the spectra, the technique can be applied to potentially tens to hundreds of sources simultaneously, over the field of view of the camera. We examined this method for a field in the centre of the Galaxy, and produced a map showing the variation of 3.4 $\mu$m optical depth across it.


title: The role of delocalized chemical bonding in square-net-based topological semimetals

abstract: Principles that predict reactions or properties of materials define the discipline of chemistry. In this work we derive chemical rules, based on atomic distances and chemical bond character, which predict topological materials in compounds that feature the structural motif of a square-net. Using these rules we identify over 300 potential new topological materials. We show that simple chemical heuristics can be a powerful tool to characterize topological matter. In contrast to previous database-driven materials categorization our approach allows us to identify candidates that are alloys, solid-solutions, or compounds with statistical vacancies. While previous material searches relied on density functional theory, our approach is not limited by this method and could also be used to discover magnetic and statistically-disordered topological semimetals.


title: Minimum Hardware Requirements for Hybrid Quantum-Classical DMFT

abstract: We numerically emulate noisy intermediate-scale quantum (NISQ) devices and determine the minimal hardware requirements for two-site hybrid quantum-classical dynamical mean-field theory (DMFT). We develop a circuit recompilation algorithm which significantly reduces the number of quantum gates of the DMFT algorithm and find that the quantum-classical algorithm converges if the two-qubit gate fidelities are larger than 99%. The converged results agree with the exact solution within 10%, and perfect agreement within noise-induced error margins can be obtained for two-qubit gate fidelities exceeding 99.9%. By comparison, the quantum-classical algorithm without circuit recompilation requires a two-qubit gate fidelity of at least 99.999% to achieve perfect agreement with the exact solution. We thus find quantum-classical DMFT calculations can be run on the next generation of NISQ devices if combined with the recompilation techniques developed in this work.


title: Neural network wave functions and the sign problem

abstract: Neural quantum states (NQS) are a promising approach to study many-body quantum physics. However, they face a major challenge when applied to lattice models: Convolutional networks struggle to converge to ground states with a nontrivial sign structure. We tackle this problem by proposing a neural network architecture with a simple, explicit, and interpretable phase ansatz, which can robustly represent such states and achieve state-of-the-art variational energies for both conventional and frustrated antiferromagnets. In the latter case, our approach uncovers low-energy states that exhibit the Marshall sign rule and are therefore inconsistent with the expected ground state. Such states are the likely cause of the obstruction for NQS-based variational Monte Carlo to access the true ground states of these systems. We discuss the implications of this observation and suggest potential strategies to overcome the problem.

categories: cond-mat.str-el cond-mat.dis-nn cs.LG physics.comp-ph quant-ph


title: Polarization constraints on the geometry of the magnetic field in the external shock of gamma-ray bursts

abstract: We study the ensemble of linear polarization measurement in the optical afterglows of long-duration gamma-ray bursts. We assume a non sideways-expanding top-hat jet geometry and use the relatively large number of measurements under the assumption that they represent a statistically unbiased sample. This allows us to constrain the ratio between the maximum predicted polarization and the measured one, which is an indicator of the geometry of the magnetic field in the downstream region of the external shock. We find that the measured polarization is substantially suppressed with respect to the maximum possible for either a completely ordered magnetic field parallel to the shock normal or to a field that is entirely contained in the shock plane. The measured polarization is limited, on average, to between 25 and 30% of the maximum theoretically possible value. This reduction requires the perpendicular component of the magnetic field to be dominant in energy with respect to the component parallel to the shock front, as expected for a shock generated and/or shock compressed field. We find, however, that the data only marginally support the assumption of a simple top-hat jet, pointing towards a more complex geometry for the outflow.


title: New Sensitivity Curves for Gravitational-Wave Signals from Cosmological Phase Transitions

abstract: Gravitational waves (GWs) from strong first-order phase transitions (SFOPTs) in the early Universe are a prime target for upcoming GW experiments. In this paper, I construct novel peak-integrated sensitivity curves (PISCs) for these experiments, which faithfully represent their projected sensitivities to the GW signal from a cosmological SFOPT by explicitly taking into account the expected shape of the signal. Designed to be a handy tool for phenomenologists and model builders, PISCs allow for a quick and systematic comparison of theoretical predictions with experimental sensitivities, as I illustrate by a large range of examples. PISCs also offer several advantages over the conventional power-law-integrated sensitivity curves (PLISCs); in particular, they directly encode information on the expected signal-to-noise ratio for the GW signal from a SFOPT. I provide semianalytical fit functions for the exact numerical PISCs of LISA, DECIGO, and BBO. In an appendix, I moreover present a detailed review of the strain noise power spectra of a large number of GW experiments. The numerical results for all PISCs, PLISCs, and strain noise power spectra presented in this paper can be downloaded from the Zenodo online repository [https://doi.org/10.5281/zenodo.3689582]. In a companion paper [1909.11356], the concept of PISCs is used to perform an in-depth study of the GW signal from the cosmological phase transition in the real-scalar-singlet extension of the standard model. The PISCs presented in this paper will need to be updated whenever new theoretical results on the expected shape of the signal become available. The PISC approach is therefore suited to be used as a bookkeeping tool to keep track of the theoretical progress in the field.


title: Collider signatures of dark CP-violation

abstract: We study an extension of the Standard Model (SM) in which two copies of the SM-Higgs doublet which do not acquire a vacuum expectation value, and hence are {\it inert}, are added to the scalar sector. The lightest particle from the inert sector, which is protected from decaying to SM particles through the conservation of a $Z_2$ symmetry, is a viable dark matter candidate. We allow for CP-violation in this extended dark sector and evaluate the $ZZZ$ vertex and its CP-violating form factor in several benchmark scenarios. We provide collider signatures of this dark CP-violation in the form of potentially observable asymmetries and cross sections for the $f\bar f\to Z^*\to ZZ$ process at both leptonic and hadronic machines.


title: Cusp and collinear anomalous dimensions in four-loop QCD from form factors

abstract: We calculate the complete quark and gluon cusp anomalous dimensions in four-loop massless QCD analytically from first principles. In addition, we determine the complete matter dependence of the quark and gluon collinear anomalous dimensions. Our approach is to Laurent expand four-loop quark and gluon form factors in the parameter of dimensional regularization. We employ finite field and syzygy techniques to reduce the relevant Feynman integrals to a basis of finite integrals, and subsequently evaluate the basis integrals directly from their standard parametric representations.


title: Multiwavelength Follow-up of the Hyperluminous Intermediate-mass Black Hole Candidate 3XMM J215022.4-055108

abstract: We recently discovered the X-ray/optical outbursting source 3XMM J215022.4-055108. It was best explained as the tidal disruption of a star by an intermediate-mass black hole of mass of a few tens of thousand solar masses in a massive star cluster at the outskirts of a large barred lenticular galaxy at D_L=247 Mpc. However, we could not completely rule out a Galactic cooling neutron star as an alternative explanation for the source. In order to further pin down the nature of the source, we have obtained new multiwavelength observations by XMM-Newton and Hubble Space Telescope (HST). The optical counterpart to the source in the new HST image is marginally resolved, which rules out the Galactic cooling neutron star explanation for the source and suggests a star cluster of half-light radius ~27 pc. The new XMM-Newton observation indicates that the luminosity was decaying as expected for a tidal disruption event and that the disk was still in the thermal state with a super-soft X-ray spectrum. Therefore, the new observations confirm the source as one of the best intermediate-mass black hole candidates.


title: On the consistency of a class of R-symmetry gauged 6D N=(1,0) supergravities

abstract: R-symmetry gauged 6D (1,0) supergravities free from all local anomalies, with gauge groups $G\times G_R$ where $G_R$ is the R-symmetry group and $G$ is semisimple with rank greater than one, and which have no hypermultiplet singlets, are extremely rare. There are three such models known in which the gauge symmetry group is $G_1\times G_2 \times U(1)_R$ where the first two factors are $ \left(E_6/{\mathbb{Z}_3}\right) \times E_7$, $ G_2 \times E_7 $ and $F_4 \times Sp(9)$. These are models with single tensor multiplet, and hyperfermions in the $(1,912)$, $(14,56)$ and $(52,18)$ dimensional representations of $G_1\times G_2$, respectively. So far it is not known if these models follow from string theory. We highlight key properties of these theories, and examine constraints which may arise from the consistency of the quantization of anomaly coefficients formulated in their strongest form by Monnier and Moore. Assuming that the gauged models accommodate dyonic string excitations, we find that these constraints are satisfied only by the model with the $F_4 \times Sp(9)\times U(1)_R$ symmetry. We also discuss aspects of dyonic strings and potential caveats they may pose in applying the stated consistency conditions to the $R$-symmetry gauged models.


title: Identification of symmetry-protected topological states on noisy quantum computers

abstract: Identifying topological properties is a major challenge because, by definition, topological states do not have a local order parameter. While a generic solution to this challenge is not available yet, a broad class of topological states, namely symmetry-protected topological (SPT) states, can be identified by distinctive degeneracies in their entanglement spectrum. Here, we propose and realize two complementary protocols to probe these degeneracies based on, respectively, symmetry-resolved entanglement entropies and measurement-based computational algorithms. The two protocols link quantum information processing to the classification of SPT phases of matter. They invoke the creation of a cluster state, and are implemented on an IBM quantum computer. The experimental findings are compared to noisy simulations, allowing us to study the stability of topological states to perturbations and noise.


title: Dynamical friction-driven orbital circularisation in rotating discs: a semi-analytical description

abstract: We present and validate a novel semi-analytical approach to study the effect of dynamical friction on the orbits of massive perturbers in rotating stellar discs. We find that dynamical friction efficiently circularises the orbit of co-rotating perturbers, while it constantly increases the eccentricity of counter-rotating ones until their angular momenta reverse, then once again promoting circularisation. Such "drag toward circular corotation" could shape the distribution of orientations of kinematically decoupled cores in disc galaxies, naturally leading to the observed larger fraction of co-rotating cores.


title: Keeping it Cool: Much Orbit Migration, yet Little Heating, in the Galactic Disk

abstract: A star in the Milky Way's disk can now be at a Galactocentric radius quite distant from its birth radius for two reasons: either its orbit has become eccentric through radial heating, which increases its radial action $J_R$ (`blurring'); or merely its angular momentum $L_z$ has changed and thereby its guiding radius (`churning'). We know that radial orbit migration is strong in the Galactic low-$\alpha$ disk and set out to quantify the relative importance of these two effects, by devising and applying a parameterized model for the distribution $p(L_z, J_R, \tau, \mathrm[Fe/H])$ in the stellar disk. This model describes the orbit evolution for stars of age $\tau$ and metallicity [Fe/H], presuming coeval stars were initially born on (near-)circular orbits, and with a unique [Fe/H] at a given birth angular momentum and age. We fit this model to APOGEE red clump stars, accounting for the complex selection function of the survey. The best fit model implies changes of angular momentum of $\sqrt{\langle \Delta L_z \rangle^2} \approx 619\, \mathrm{kpc~km/s~}(\tau/\mathrm{6~Gyr})^{0.5}$, and changes of radial action as $\sqrt{\langle \Delta J_R \rangle^2} \approx 63\, \mathrm{kpc~km/s~} (\tau/\mathrm{6~Gyr})^{0.6}$ at 8 kpc. This suggests that the secular orbit evolution of the disk is dominated by diffusion in angular momentum, with radial heating being an order of magnitude lower.


title: Quark Flavor Phenomenology of the QCD Axion

abstract: Axion models with generation-dependent Peccei-Quinn charges can lead to flavor-changing neutral currents, thus motivating QCD axion searches at precision flavor experiments. We rigorously derive limits on the most general effective flavor-violating couplings from current measurements and assess their discovery potential. For two-body decays we use available experimental data to derive limits on $q\to q' a$ decay rates for all flavor transitions. Axion contributions to neutral-meson mixing are calculated in a systematic way using chiral perturbation theory and operator product expansion. We also discuss in detail baryonic decays and three-body meson decays, which can lead to the best search strategies for some of the couplings. For instance, a strong limit on the $\Lambda\to n a$ transition can be derived from the supernova SN 1987A. In the near future, dedicated searches for $q\to q' a$ decays at ongoing experiments could potentially test Peccei-Quinn breaking scales up to $10^{12}$ GeV at NA62 or KOTO, and up to $10^{9}$ GeV at Belle II or BES III.


title: Tropical fans, scattering equations and amplitudes

abstract: We describe a family of tropical fans related to Grassmannian cluster algebras. These fans are related to the kinematic space of massless scattering processes in a number of ways. For each fan associated to the Grassmannian ${\rm Gr}(k,n)$ there is a notion of a generalised $\phi^3$ amplitude and an associated set of scattering equations which further generalise the ${\rm Gr}(k,n)$ scattering equations that have been recently introduced. Here we focus mostly on the cases related to finite Grassmannian cluster algebras and we explain how face variables for the cluster polytopes are simply related to the scattering equations. For the Grassmannians ${\rm Gr}(4,n)$ the tropical fans we describe are related to the singularities (or symbol letters) of loop amplitudes in planar $\mathcal{N}=4$ super Yang-Mills theory. We show how each choice of tropical fan leads to a natural class of polylogarithms, generalising the notion of cluster adjacency and we describe how the currently known loop data fit into this classification.


title: Dynamics of millicharged dark matter in supernova remnants

abstract: Dark matter could have a small electromagnetic charge, provided the charge-to-mass ratio is much less than that of electrons or protons. This candidate, commonly known as millicharged dark matter (mDM), would form a plasma and interact with the interstellar medium and electromagnetic fields within galaxies. In general, understanding the dynamics of mDM requires consideration of collective plasma effects. It has been proposed that mDM can be accelerated in supernova remnants, forming a dark cosmic ray population that would leave distinct experimental signatures. In this work, we study a microphysical model where mDM is shocked by a supernova remnant and isotropized in the frame of the expanding fluid. We find that for $|q_\chi/m_\chi| \gtrsim 10^{-13} e / {\rm MeV}$, the isotropization length for electromagnetic plasma instabilities is much shorter than the size of the supernova remnant. This is a necessary, though not sufficient, first step for formation of a Fermi-accelerated mDM component, and determining the size of this component requires further study. We discuss additional implications of mDM interactions in supernova remnants.


title: Finding novelty with uncertainty

abstract: Medical images are often used to detect and characterize pathology and disease; however, automatically identifying and segmenting pathology in medical images is challenging because the appearance of pathology across diseases varies widely. To address this challenge, we propose a Bayesian deep learning method that learns to translate healthy computed tomography images to magnetic resonance images and simultaneously calculates voxel-wise uncertainty. Since high uncertainty occurs in pathological regions of the image, this uncertainty can be used for unsupervised anomaly segmentation. We show encouraging experimental results on an unsupervised anomaly segmentation task by combining two types of uncertainty into a novel quantity we call scibilic uncertainty.


title: Robust dynamical exchange cooling with trapped ions

abstract: We investigate theoretically the possibility for robust and fast cooling of a trapped atomic ion by transient interaction with a pre-cooled ion. The transient coupling is achieved through dynamical control of the ions' equilibrium positions. To achieve short cooling times we make use of shortcuts to adiabaticity by applying invariant-based engineering. We design these to take account of imperfections such as stray fields, and trap frequency offsets. For settings appropriate to a currently operational trap in our laboratory, we find that robust performance could be achieved down to $6.3$ motional cycles, comprising $14.2\ \mathrm{\mu s}$ for ions with a $0.44\ \mathrm{MHz}$ trap frequency. This is considerably faster than can be achieved using laser cooling in the weak coupling regime, which makes this an attractive scheme in the context of quantum computing.


title: Topology Classification from Chiral Symmetry: Chiral Phase Index and Spin Correlations in Graphene Nanoribbons

abstract: Topology concepts have significantly deepened of our understanding in recent years of the electronic properties of one-dimensional (1D) nano structures such as the graphene nanoribbons. Controlling topological electronic properties of GNRs has been demonstrated in both theoretical studies and experimental realization. Most previous works rely on classification theory requiring both time reversal and spatial symmetry of a unit cell in the 1D bulk material that is commensurate to its boundary. To access boundary structures that lead to unit cell with no spatial symmetry and to generalize the theory, we propose here another classification scheme, using chiral symmetry, to arrive at a Z classification that is not only applicable to GNRs with arbitrary terminations, but also to any general 1D chiral structures. This theory, combining with Lieb's theorem, moreover enables access to the electron's spin degree of freedom, allowing for investigation of spin physics.


title: Toward the observation of interference effects in nonlinear Compton scattering

abstract: The photon spectrum from electrons scattering on multiple laser pulses exhibits interference effects not present for scattering on a single pulse. We investigate the conditions required for the experimental observation of these interference effects in electron-laser collisions, in particular analysing the roles of the detector resolution and the transverse divergence of the pump electron beam.


title: Single supply of electrical energy in the supply distribution model: using GAMS to model the effects of network parameters on the level of consumer demand

abstract: The proposed model of this study is a single supply of electrical energy and is used in distribution systems. The objective of this study is to optimize the distribution of active and reactive energy in the supply of sub-distribution or the market. The model proposed in this paper accounts for the effects associated with switching capacitors on distribution and dispatching centers. Using capacitors in the power grid has an effect on bus voltage and consequently the amount of active and reactive energy demand of the consumer. The optimized supply is determined through the solution of a mathematical operation model. In this study, the effect of voltage changes on consumer demand for residential,commercial and industrial power are quantified and embodied in an analytical model, which is then solved using General Algebraic Modeling System (GAMS) software.ware.


title: Ask the Experts: What Should Be on an IoT Privacy and Security Label?

abstract: Information about the privacy and security of Internet of Things (IoT) devices is not readily available to consumers who want to consider it before making purchase decisions. While legislators have proposed adding succinct, consumer accessible, labels, they do not provide guidance on the content of these labels. In this paper, we report on the results of a series of interviews and surveys with privacy and security experts, as well as consumers, where we explore and test the design space of the content to include on an IoT privacy and security label. We conduct an expert elicitation study by following a three-round Delphi process with 22 privacy and security experts to identify the factors that experts believed are important for consumers when comparing the privacy and security of IoT devices to inform their purchase decisions. Based on how critical experts believed each factor is in conveying risk to consumers, we distributed these factors across two layers---a primary layer to display on the product package itself or prominently on a website, and a secondary layer available online through a web link or a QR code. We report on the experts' rationale and arguments used to support their choice of factors. Moreover, to study how consumers would perceive the privacy and security information specified by experts, we conducted a series of semi-structured interviews with 15 participants, who had purchased at least one IoT device (smart home device or wearable). Based on the results of our expert elicitation and consumer studies, we propose a prototype privacy and security label to help consumers make more informed IoT-related purchase decisions.


title: Black-Box Optimization with Local Generative Surrogates

abstract: We propose a novel method for gradient-based optimization of black-box simulators using differentiable local surrogate models. In fields such as physics and engineering, many processes are modeled with non-differentiable simulators with intractable likelihoods. Optimization of these forward models is particularly challenging, especially when the simulator is stochastic. To address such cases, we introduce the use of deep generative models to iteratively approximate the simulator in local neighborhoods of the parameter space. We demonstrate that these local surrogates can be used to approximate the gradient of the simulator, and thus enable gradient-based optimization of simulator parameters. In cases where the dependence of the simulator on the parameter space is constrained to a low dimensional submanifold, we observe that our method attains minima faster than baseline methods, including Bayesian optimization, numerical optimization, and approaches using score function gradient estimators.


title: Sun-like Stars Shed Light on Solar Climate Forcing

abstract: Recently published, precise stellar photometry of 72 Sun-like stars obtained at the Fairborn Observatory between 1993 and 2017 is used to set limits on the solar forcing of Earth's atmosphere of $\pm$ 4.5 W m$^{-2}$ since 1750. This compares with the +2.2 $\pm$ 1.1 W m$^{-2}$ IPCC estimate for anthropogenic forcing. Three critical assumptions are made. In decreasing order of importance they are: (a) most of the brightness variations occur within the average time-series length of $\approx$17 years; (b) the Sun seen from the ecliptic behaves as an ensemble of middle-aged solar-like stars; and (c) narrow-band photometry in the Str\"omgren $b$ and $y$ bands are linearly proportional to the total solar irradiance. Assumption (a) can best be relaxed and tested by obtaining more photometric data of Sun-like stars, especially those already observed. Eight stars with near-solar parameters have been observed from 1999, and two since 1993. Our work reveals the importance of continuing and expanding ground-based photometry, to complement expensive solar irradiance measurements from space.


title: Neuroevolution of Neural Network Architectures Using CoDeepNEAT and Keras

abstract: Machine learning is a huge field of study in computer science and statistics dedicated to the execution of computational tasks through algorithms that do not require explicit instructions but instead rely on learning patterns from data samples to automate inferences. A large portion of the work involved in a machine learning project is to define the best type of algorithm to solve a given problem. Neural networks - especially deep neural networks - are the predominant type of solution in the field. However, the networks themselves can produce very different results according to the architectural choices made for them. Finding the optimal network topology and configurations for a given problem is a challenge that requires domain knowledge and testing efforts due to a large number of parameters that need to be considered. The purpose of this work is to propose an adapted implementation of a well-established evolutionary technique from the neuroevolution field that manages to automate the tasks of topology and hyperparameter selection. It uses a popular and accessible machine learning framework - Keras - as the back-end, presenting results and proposed changes concerning the original algorithm. The implementation is available at GitHub (https://github.com/sbcblab/Keras-CoDeepNEAT) with documentation and examples to reproduce the experiments performed for this work.


title: hybrid-Lenstool: A self-consistent algorithm to model galaxy clusters with strong- and weak-lensing simultaneously

abstract: We present a new galaxy cluster lens modeling approach, hybrid-Lenstool, that is implemented in the publicly available modeling software Lenstool. hybrid-Lenstool combines a parametric approach to model the core of the cluster, and a non-parametric (free-form) approach to model the outskirts. hybrid-Lenstool optimizes both strong- and weak-lensing constraints simultaneously (Joint-Fit), providing a self-consistent reconstruction of the cluster mass distribution on all scales. In order to demonstrate the capabilities of the new algorithm, we tested it on a simulated cluster. hybrid-Lenstool yields more accurate reconstructed mass distributions than the former Sequential-Fit approach where the parametric and the non-parametric models are optimized successively. Indeed, we show with the simulated cluster that the mass density profile reconstructed with a Sequential-Fit deviates form the input by $2-3\sigma$ at all scales while the Joint-Fit gives a profile that is within $1-1.5\sigma$ of the true value. This gain in accuracy is consequential for recovering mass distributions exploiting cluster lensing and therefore for all applications of clusters as cosmological probes. Finally we found that the Joint-Fit approach yields shallower slope of the inner density profile than the Sequential-Fit approach, thus revealing possible biases in previous lensing studies.


title: On numerical approximations to fluid-structure interactions involving compressible fluids

abstract: In this paper we introduce a numerical scheme for fluid-structure interaction problems in two or three space dimensions: A flexible elastic plate is interacting with a viscous, compressible barotropic fluid. Hence the physical domain of definition (the domain of Eulerian coordinates) is changing in time. We introduce a fully discrete scheme that is stable, satisfies geometric conservation, mass conservation and the positivity of the density. We also prove that the scheme is consistent with the definition of continuous weak solutions.


title: Modular Hamiltonian for (holographic) excited states

abstract: In this work we study the Tomita-Takesaki construction for a family of excited states that, in a strongly coupled CFT - at large $N$-, correspond to coherent states in an asymptotically AdS spacetime geometry. We compute the modular flow and modular Hamiltonian associated to these excited states in the Rindler wedge and for a ball shaped entangling surface. Using holography, one can compute the bulk modular flow and construct the Tomita-Takesaki theory for these cases. We also discuss generalizations of the entanglement regions in the bulk and how to estimate the modular Hamiltonian in a large N approximation. Finally we present a holographic formula, based on the BDHM prescription, to compute the modular evolution of operators in the corresponding CFT algebra.


title: A polynomial time parallel algorithm for graph isomorphism using a quasipolynomial number of processors

abstract: The Graph Isomorphism (GI) problem is a theoretically interesting problem because it has not been proven to be in P nor to be NP-complete. Babai made a breakthrough in 2015 when announcing a quasipolynomial time algorithm for GI problem. Babai's work gives the most theoretically efficient algorithm for GI, as well as a strong evidence favoring the idea that class GI $\ne$ NP and thus P $\ne$ NP. Based on Babai's algorithm, we prove that GI can further be solved by a parallel algorithm that runs in polynomial time using a quasipolynomial number of processors. We achieve that result by identifying the bottlenecks in Babai's algorithms and parallelizing them. In particular, we prove that color refinement can be computed in parallel logarithmic time using a polynomial number of processors, and the $k$-dimensional WL refinement can be computed in parallel polynomial time using a quasipolynomial number of processors. Our work suggests that Graph Isomorphism and GI-complete problems can be computed efficiently in a parallel computer, and provides insights on speeding up parallel GI programs in practice.


title: Validating uncertainty in medical image translation

abstract: Medical images are increasingly used as input to deep neural networks to produce quantitative values that aid researchers and clinicians. However, standard deep neural networks do not provide a reliable measure of uncertainty in those quantitative values. Recent work has shown that using dropout during training and testing can provide estimates of uncertainty. In this work, we investigate using dropout to estimate epistemic and aleatoric uncertainty in a CT-to-MR image translation task. We show that both types of uncertainty are captured, as defined, providing confidence in the output uncertainty estimates.


title: Debugging Machine Learning Pipelines

abstract: Machine learning tasks entail the use of complex computational pipelines to reach quantitative and qualitative conclusions. If some of the activities in a pipeline produce erroneous or uninformative outputs, the pipeline may fail or produce incorrect results. Inferring the root cause of failures and unexpected behavior is challenging, usually requiring much human thought, and is both time-consuming and error-prone. We propose a new approach that makes use of iteration and provenance to automatically infer the root causes and derive succinct explanations of failures. Through a detailed experimental evaluation, we assess the cost, precision, and recall of our approach compared to the state of the art. Our source code and experimental data will be available for reproducibility and enhancement.


title: Root lattices in number fields

abstract: We explore whether a root lattice may be similar to the lattice $\mathscr O$ of integers of a number field $K$ endowed with the inner product $(x, y):={\rm Trace}_{K/\mathbb Q}(x\cdot\theta(y))$, where $\theta$ is an involution of $K$. We classify all pairs $K$, $\theta$ such that $\mathscr O$ is similar to either an even root lattice or the root lattice $\mathbb Z^{[K:\mathbb Q]}$. We also classify all pairs $K$, $\theta$ such that $\mathscr O$ is a root lattice. In addition to this, we show that $\mathscr O$ is never similar to a positive-definite even unimodular lattice of rank $\leqslant 48$, in particular, $\mathscr O$ is not similar to the Leech lattice. In appendix, we give a general cyclicity criterion for the primary components of the discriminant group of $\mathscr O$.


title: Some Inequalities Related to Strong Convergence of Riesz Logarithmic Means

abstract: In this paper we derive a new strong convergence theorem of Riesz logarithmic means of the one-dimensional Vilenkin-Fourier (Walsh-Fourier) series. The corresponding inequality is pointed out and it is also proved that the inequality is in a sense sharp, at least for the case with Walsh-Fourier series.


title: Did Simon Marius observe Jupiter's satellites on January 8, 1610? An exercise in computation

abstract: The question, Did Simon Marius (1573-1625) observe Jupiter's satellites on January 8, 1610 (December 29, 1609 in the Julian calendar) is moot, for he did not disclose his research method and the instrument he used. To resolve this issue we apply astronomical codes and evaluate the visual performance of a replica of the telescope that Galileo Galilei (1564-1642) had used.


title: Optically tunable microresonator using an azobenzene monolayer

abstract: Photoswitchable organic molecules can undergo reversible structural changes with an external light stimulus. These optically controlled molecules have been used in the development of smart polymers, optical writing of grating films, and even controllable in-vivo drug release. Being the simplest class of photoswitches in terms of structure, azobenzenes have become the most ubiquitous, well-characterized, and implemented organic molecular switch. Given their predictable response, they are ideally suited to create an all-optically controlled switch. However, fabricating a monolithic optical device comprised solely from azobenzene while maintaining the photoswitching functionality is challenging. In this work, we combine integrated photonics with optically switchable organic molecules to create an optically controlled integrated device. A silica toroidal resonant cavity is functionalized with a monolayer of an azobenzene derivative. After functionalization, the loaded cavity Q is above 100,000. When 450 nm light is coupled into cavity resonance, the azobenzene isomerizes from trans-isomer to cis-isomer, inducing a refractive index change. Because the resonant wavelength of the cavity is governed by the index, the resonant wavelength changes in parallel. At the probe wavelength of 1300 nm, the wavelength shift is determined by the duration and intensity of the 450 nm light and the density of azobenzene functional groups on the device surface, providing multiple control mechanisms. Using this photoswitchable device, resonance frequency tuning as far as sixty percent of the cavity free spectral range in the near-IR is demonstrated. The kinetics of the tuning are in agreement with spectroscopic and ellipsometry measurements coupled with finite element method calculations.


title: YZiCS: Unveiling Quenching History of Cluster Galaxies Using Phase-space Analysis

abstract: We used the time since infall (TSI) of galaxies, obtained from the Yonsei Zoom-in Cluster Simulation, and the star formation rate (SFR) from the Sloan Digital Sky Survey (SDSS) Data Release 10 to study how quickly star formation of disk galaxies is quenched in cluster environments. We first confirm that both simulated and observed galaxies are consistently distributed in phase space. We then hypothesize that the TSI and SFR are causally connected; thus, both the TSI and SFR of galaxies at each position of phase space can be associated through abundance matching. Using a flexible model, we derive the star formation history (SFH) of cluster galaxies that best reproduces the relationship between the TSI and SFR at $z\sim 0.08$. According to this SFH, we find that the galaxies with $M_{*} > 10^{9.5} M_{\odot}$ generally follow the so-called "delayed-then-rapid" quenching pattern. Our main results are as following: (i) Part of the quenching takes place outside clusters through mass quenching and pre-processing. The e-folding timescale of this "$ex\text{-}situ$ quenching phase" is roughly 3 Gyr with a strong inverse mass dependence. (ii) The pace of quenching is maintained roughly for 2 Gyr ("delay time") during the first crossing time into the cluster. During the delay time, quenching remains gentle probably because gas loss happens primarily on hot and neutral gases. (iii) Quenching becomes more dramatic (e-folding timescale of roughly 1 Gyr) after delay time, probably because ram pressure stripping is strongest near the cluster center. Counter-intuitively, more massive galaxies show shorter quenching timescales mainly because they enter their clusters with lower gas fractions due to $ex\text{-}situ$ quenching.


title: Short-term variability and mass loss in Be stars V. Space photometry and ground-based spectroscopy of $\gamma$ Cas

abstract: Context. Be stars are physically complex systems that continue to challenge theory to understand their rapid rotation, complex variability and decretion disks. $\gamma$ Cassiopeiae ($\gamma$ Cas) is one such star but is even more curious because of its unexplained hard thermal X-ray emission. Aims. We aim to examine the optical variability of $\gamma$ Cas and thereby to shed more light on its puzzling behaviour. Methods. Three hundred twenty-one archival H$\alpha$ spectra from 2006 to 2017 are analysed to search for frequencies corresponding to the 203.5 day orbit of the companion. Space photometry from the SMEI satellite from 2003 to 2011 and the BRITE-Constellation of nano-satellites between 2015 and 2019 is investigated in the period range from a couple of hours to a few days. Results. The orbital period of the companion of 203.5 days is confirmed with independent measurements from the structure of the H$\alpha$ line emission. A strong blue/red asymmetry in the amplitude distribution across the H$\alpha$ emission line could hint at a spiral structure in the decretion disk. With the space photometry, the known frequency of 0.82 d$^{-1}$ is confirmed in data from the early 2000s. A higher frequency of 2.48 d$^{-1}$ is present in the data from 2015 to 2019 and possibly also in the early 2000s. A third frequency at 1.25 d$^{-1}$ is proposed to exist in both SMEI and BRITE data. The only explanation covering all three rapid variations seems to be nonradial pulsation. The two higher frequencies are incompatible with rotation.


title: A homological approach to pseudoisotopy theory. I

abstract: We construct a zig-zag from the once delooped space of pseudoisotopies of a closed $2n$-disc to the once looped algebraic $K$-theory space of the integers and show that the maps involved are $p$-locally $(2n-4)$-connected for $n>3$ and large primes $p$. The proof uses the computation of the stable homology of the moduli space of high-dimensional handlebodies due to Botvinnik--Perlmutter and is independent of the classical approach to pseudoisotopy theory based on Igusa's stability theorem and work of Waldhausen. Combined with a result of Randal-Williams, one consequence of this identification is a calculation of the rational homotopy groups of $\mathrm{BDiff}_\partial(D^{2n+1})$ in degrees up to $2n-5$.


title: Information Freshness for Timely Detection of Status Changes

abstract: In this paper, we aim to establish the connection between Age of Information (AoI) in network theory, information uncertainty in information theory, and detection delay in time series analysis. We consider a dynamic system whose state changes at discrete time points, and a state change won't be detected until an update generated after the change point is delivered to the destination for the first time. We introduce an information theoretic metric to measure the information freshness at the destination, and name it as generalized Age of Information (GAoI). We show that under any state-independent online updating policy, if the underlying state of the system evolves according to a stationary Markov chain, the GAoI is proportional to the AoI. Besides, the accumulative GAoI and AoI are proportional to the expected accumulative detection delay of all changes points over a period of time. Thus, any (G)AoI-optimal state-independent updating policy equivalently minimizes the corresponding expected change point detection delay, which validates the fundamental role of (G)AoI in real-time status monitoring. Besides, we also investigate a Bayesian change point detection scenario where the underlying state evolution is not stationary. Although AoI is no longer related to detection delay explicitly, we show that the accumulative GAoI is still an affine function of the expected detection delay, which indicates the versatility of GAoI in capturing information freshness in dynamic systems.


title: Effects of Mirror Seeing on High-Contrast Adaptive Optics Instruments

abstract: Ground-based direct imaging surveys like the Gemini Planet Imager Exoplanet Survey (GPIES) rely on Adaptive Optics (AO) systems to image and characterize exoplanets that are up to a million times fainter than their host stars. One factor that can reduce AO performance is turbulence induced by temperature differences in the instrument's immediate surroundings (e.g.: "dome seeing" or "mirror seeing"). In this analysis we use science observations, AO telemetry, and environmental data from September 2014 to February 2017 of the GPIES campaign to quantify the effects of "mirror seeing" on the performance of the GPI instrument. We show that GPI performance is optimal when the primary mirror (M1) is in equilibrium with the outside air temperature. We then examine the characteristics of mirror seeing by calculating the power spectral densities (PSD) of spatial and temporal Fourier modes. Inside the inertial range of the PSDs, we find that the spatial PSD amplitude increases when M1 is out of equilibrium and that the integrated turbulence may exhibit deviations from Kolmogorov atmospheric turbulence models and from the 1-layer frozen flow model. We conclude with an assessment of the current temperature control and ventilation strategy at Gemini South.


title: Optical levitation using broadband light

abstract: The ability to create dynamic, tailored optical potentials has become important across fields ranging from biology to quantum science. We demonstrate a method for the creation of arbitrary optical tweezer potentials using the broadband spectral profile of a superluminescent diode combined with the chromatic aberration of a lens. A tunable filter, typically used for ultra-fast laser pulse shaping, allows us to manipulate the broad spectral profile and therefore the optical tweezer potentials formed by focusing of this light. We characterize these potentials by measuring the Brownian motion of levitated nanoparticles in vacuum and, also demonstrate interferometric detection and feedback cooling of the particle,s motion. This simple and cost-effective technique will enable a wide range of applications and allow rapid modulation of the optical potential landscape in excess of MHz frequencies.

categories: physics.optics nlin.CD physics.atm-clus quant-ph


title: Chemistry Along Accretion Streams in a Viscously-Evolving Protoplanetary Disk

abstract: The composition of a protoplanetary disk is set by a combination of interstellar inheritance and gas and grain surface chemical reactions within the disk. The survival of inherited molecules, as well as the disk in situ chemistry depends on the local temperature, density and irradiation environment, which can change over time due to stellar and disk evolution, as well as transport in the disk. We address one aspect of this coupling between the physical and chemical evolution in disks by following accretion streamlines of gas and small grains in the disk midplane, while simultaneously taking the evolving star into account. This approach is computationally efficient and enables us to take into account changing physical conditions without reducing the chemical network. We find that many species are enhanced in the inner disk midplane in the dynamic model due to inward transport of cosmic-ray driven chemical products, resulting in, e.g., orders-of magnitude hydrocarbon enhancements at 1 au, compared to a static disk. For several other chemical families, there is no difference between the static and dynamic models, indicative of a robust chemical reset, while yet others show differences between static and dynamic models that depend on complex interactions between physics and chemistry during the inward track. The importance of coupling dynamics and chemistry when modeling the chemical evolution of protoplanetary disks is thus depends on what chemistry is of interest.


title: Two major constraints on the inner radii of accretion disks

abstract: The Stefan-Boltzmann law yields a fundamental constraint on the geometry of inner accretion disks in black-hole X-ray binaries. It follows from considering the irradiating flux and the effective temperature of the inner parts of the disk, which implies that a strong quasi-thermal component with the average energy higher than that of a blackbody at the effective temperature has to be present whenever relativistic Fe K fluorescence and reflection features are observed. The apparent absence of such quasi-thermal component with the color temperature of $\sim$1 keV in high-luminosity hard states is not compatible with a strongly irradiated disk extending close to the innermost stable circular orbit. Instead, the disk should be either truncated at a relatively large radius or irradiated by a corona at a large height, which would reduce the effective temperature and bring it to an agreement with the data. We also study constraints on disk/corona models following from comparing the disk densities fitted in literature using variable-density reflection codes with those calculated by us from the ionization parameter, the luminosity and the disk inner radius. We find that the fitted densities are much higher/lower in the hard/soft state of binaries, implying significant problems with the used assumptions and methods.


title: Entropy-stable, high-order summation-by-parts discretizations without interface penalties

abstract: The paper presents high-order accurate, energy-, and entropy-stable discretizations constructed from summation-by-parts (SBP) operators. Notably, the discretizations assemble global SBP operators and use continuous solutions, unlike previous efforts that use discontinuous SBP discretizations. Derivative-based dissipation and local-projection stabilization (LPS) are investigated as options for stabilizing the baseline discretization. These stabilizations are equal up to a multiplicative constant in one dimension, but only LPS remains well conditioned for general, multidimensional SBP operators. Furthermore, LPS is able to take advantage of the additional nodes required by degree $2p$ diagonal-norms, resulting in an element-local stabilization with a bounded spectral radius. An entropy-stable version of LPS is easily obtained by applying the projection on the entropy variables. Numerical experiments with the linear-advection and Euler equations demonstrate the accuracy, efficiency, and robustness of the stabilized discretizations, and the continuous approach compares favorably with the more common discontinuous SBP methods.


title: Analyzing, Exploring, and Visualizing Complex Networks via Hypergraphs using SimpleHypergraphs.jl

abstract: Real-world complex networks are usually being modeled as graphs. The concept of graphs assumes that the relations within the network are binary (for instance, between pairs of nodes); however, this is not always true for many real-life scenarios, such as peer-to-peer communication schemes, paper co-authorship, or social network interactions. For such scenarios, it is often the case that the underlying network is better and more naturally modeled by hypergraphs. A hypergraph is a generalization of a graph in which a single (hyper)edge can connect any number of vertices. Hypergraphs allow modelers to have a complete representation of multi-relational (many-to-many) networks; hence, they are extremely suitable for analyzing and discovering more subtle dependencies in such data structures. Working with hypergraphs requires new software libraries that make it possible to perform operations on them, from basic algorithms (such as searching or traversing the network) to computing significant hypergraph measures, to including more challenging algorithms (such as community detection). In this paper, we present a new software library, SimpleHypergraphs.jl, written in the Julia language and designed for high-performance computing on hypergraphs and propose two new algorithms for analyzing their properties: s-betweenness and modified label propagation. We also present various approaches for hypergraph visualization integrated into our tool. In order to demonstrate how to exploit the library in practice, we discuss two case studies based on the 2019 Yelp Challenge dataset and the collaboration network built upon the Game of Thrones TV series. The results are promising and they confirm the ability of hypergraphs to provide more insight than standard graph-based approaches.


title: Spintronics meets density matrix renormalization group: Quantum spin torque driven nonclassical magnetization reversal and dynamical buildup of long-range entanglement

abstract: We introduce time-dependent density matrix renormalization group (tDMRG) as a solution to long standing problem in spintronics -- how to describe spin-transfer torque (STT) between flowing spins of conduction electrons and localized spins within a magnetic material by treating the dynamics of both spin species fully quantum-mechanically. In contrast to conventional Slonczewski-Berger STT, where the localized spins are viewed as classical vectors obeying the Landau-Lifshitz-Gilbert equation and where their STT-driven dynamics is initiated only when the spin-polarization of flowing electrons and localized spins are noncollinear, quantum STT can occur when these vectors are collinear but antiparallel. Using tDMRG, we simulate the time evolution of a many-body quantum state of electrons and localized spins, where the former are injected as a spin-polarized current pulse while the latter comprise a quantum Heisenberg ferromagnetic metallic (FM) spin-$\frac{1}{2}$ XXZ chain initially in the ground state with spin-polarization antiparallel to that of injected electrons. The quantum STT reverses the direction of localized spins, but without rotation from the initial orientation, when the number of injected electrons exceeds the number of localized spins. Such nonclassical reversal, which is absent from LLG dynamics, is strikingly inhomogeneous across the FM chain and it can be accompanied by reduction of the magnetization associated with localized spins, even to zero at specific locations. This is because quantum STT generates a highly entangled nonequilibrium many-body state of all flowing and localized spins, despite starting from the initially unentangled ground state of a mundane FM. Furthermore, the mutual information between localized spins at the FM edges remains nonzero even at infinite separation as the signature of dynamical buildup of long-range entanglement.


title: Pseudogap Anderson impurity model out of equilibrium: A master equation tensor network approach

abstract: We study equilibrium and nonequilibrium properties of the single-impurity Anderson model with a power-law pseudogap in the density of states. In equilibrium, the model is known to display a quantum phase transition from a generalized Kondo to a local moment phase. In the present work, we focus on the extension of these phases beyond equilibrium, i.e. under the influence of a bias voltage. Within the auxiliary master equation approach combined with a scheme based on matrix product states (MPS) we are able to directly address the current-carrying steady state. Starting with the equilibrium situation, we first corroborate our results by comparing with a direct numerical evaluation of ground state spectral properties of the system by MPS. Here, a scheme to locate the phase boundary by extrapolating the power-law exponent of the self energy produces a very good agreement with previous results obtained by the numerical renormalization group. Our nonequilibrium study as a function of the applied bias voltage is then carried out for two points on either side of the phase boundary. In the Kondo regime the resonance in the spectral function is splitted as a function of the increasing bias voltage. The local moment regime, instead, displays a dip in the spectrum near the position of the chemical potentials. Similar features are observed in the corresponding self energies. The Kondo split peaks approximately obey a power-law behavior as a function of frequency, whose exponents depend only slightly on voltage. Finally, the differential conductance in the Kondo regime shows a peculiar maximum at finite voltages, whose height, however, is below the accuracy level.

categories: cond-mat.str-el cond-mat.mes-hall cond-mat.other quant-ph


title: Geometry of generalized Pauli channels

abstract: We analyze the geometry of the generalized Pauli channels constructed from the mutually unbiased bases. The Choi-Jamio{\l}kowski isomorphism allows us to express the Hilbert-Schmidt line and volume elements in terms of the eigenvalues of the generalized Pauli maps. After determining appropriate regions of integration, we analytically compute the volume of generalized Pauli channels and their important subclasses. In particular, we obtain the volumes of the generalized Pauli channels that can be generated by a legitimate generator and are entanglement breaking. We also provide the upper bound for the volume of positive, trace-preserving generalized Pauli maps.


title: A Non-Intrusive Correction Algorithm for Classification Problems with Corrupted Data

abstract: A novel correction algorithm is proposed for multi-class classification problems with corrupted training data. The algorithm is non-intrusive, in the sense that it post-processes a trained classification model by adding a correction procedure to the model prediction. The correction procedure can be coupled with any approximators, such as logistic regression, neural networks of various architectures, etc. When training dataset is sufficiently large, we prove that the corrected models deliver correct classification results as if there is no corruption in the training data. For datasets of finite size, the corrected models produce significantly better recovery results, compared to the models without the correction algorithm. All of the theoretical findings in the paper are verified by our numerical examples.


title: Resonant scattering of Dice quasiparticles on oscillating quantum dots

abstract: We consider a Dice model with Dirac cones intersected by a topologically flat band at the charge neutrality point and analyze the inelastic scattering of massless pseudospin-1 particles on a circular, gate-defined, oscillating barrier. Focusing on the resonant scattering regime at small energy of the incident wave, we calculate the reflection and transmission coefficients and derive explicit expressions for the time-dependent particle probability, current density and scattering efficiency within (Floquet) Dirac-Weyl theory, both in the near-field and the far-field. We discuss the importance of sideband scattering and Fano resonances in the quantum limit. When resonance conditions are fulfilled, the particle is temporarily trapped in vortices located close to edge of the quantum dot before it gets resubmitted with strong angular dependence. Interestingly even periodically alternating forward and backward radiation may occur. We also demonstrate the revival of resonant scattering related to specific fusiform boundary trapping profiles.


title: Some remarks on equilateral triangulations of surfaces and Belyi functions

abstract: In this paper, following Grothendieck {\it Esquisse d'un programme}, which was motivated by Belyi's work, we study some properties of surfaces $X$ which are triangulated by (possibly ideal) isometric equilateral triangles of one of the spherical, euclidean or hyperbolic geometries. These surfaces have a natural Riemannian metric with conic singularities. In the euclidean case we analyze the closed geodesics and their lengths. Such surfaces can be given the structure of a Riemann surface which, considered as algebraic curves, are defined over $\bar{\mathbb{Q}}$ by a theorem of Belyi. They have been studied by many authors of course. Here we define the notion of connected sum of two Belyi functions and give some concrete examples. In the particular case when $X$ is a torus, the triangulation leads to an elliptic curve and we define the notion of a "peel" obtained from the triangulation (which is a metaphor of an orange peel) and relate this peel with the modulus $\tau$ of the elliptic curve. Many fascinating questions arise regarding the modularity of the elliptic curve and the geometric aspects of the Taniyama-Shimura-Weil theory.


title: A theorem about maximal Cohen-Macaulay modules

abstract: It is shown in a local strongly $F$-regular ring there exits natural number $e_0$ so that if $M$ is any finitely generated maximal Cohen-Macaulay module then the pushforward of $M$ under the $e_0$th iterate of the Frobenius endomorphism contains a free summand. Consequently, the torsion subgroup of the divisor class group of a local strongly $F$-regular ring is finite.


title: Superconductivity with Tc=7 K under pressure for Cu-, and Au-doped BaFe2As2

abstract: It is noteworthy that chemical substitution of BaFe2As2 (122) with the noble elements Cu and Au gives superconductivity with a maximum Tc=3 K, while Ag substitution (Ag-122) stays antiferromagnetic. For Ba(Fe1-xTMx)2As2, TM= Cu, Au, or Ag, and by doping an amount of x=0.04, a-lattice parameter slightly increases (0.4%) for all TM dopants, while c-lattice decreases (-0.2%) for TM=Cu, barely moves (0.05%) for Au, and increases (0.2%) for Ag. Despite the naive expectation that the noble elements of group 11 should affect the quantum properties of 122 similarly, they produce significant differences extending to the character of the ground state. For the Ag-122 crystal, evidence of only a filamentary superconductivity is noted with pressure. However, for Au and Cu doping (x=0.03) we find a substantial improvement in the superconductivity, with Tc increasing to 7 K and 7.5 K, respectively, under 20 kbar of pressure. As with the ambient pressure results, the identity of the dopant therefore has a substantial impact on the ground state properties. Density functional theory calculations corroborate these results and find evidence of strong electronic scattering for Au and Ag dopants, while Cu is comparatively less disruptive to the 122 electronic structure.


title: TDEFSI: Theory Guided Deep Learning Based Epidemic Forecasting with Synthetic Information

abstract: Influenza-like illness (ILI) places a heavy social and economic burden on our society. Traditionally, ILI surveillance data is updated weekly and provided at a spatially coarse resolution. Producing timely and reliable high-resolution spatiotemporal forecasts for ILI is crucial for local preparedness and optimal interventions. We present TDEFSI (Theory Guided Deep Learning Based Epidemic Forecasting with Synthetic Information), an epidemic forecasting framework that integrates the strengths of deep neural networks and high-resolution simulations of epidemic processes over networks. TDEFSI yields accurate high-resolution spatiotemporal forecasts using low-resolution time series data. During the training phase, TDEFSI uses high-resolution simulations of epidemics that explicitly model spatial and social heterogeneity inherent in urban regions as one component of training data. We train a two-branch recurrent neural network model to take both within-season and between-season low-resolution observations as features, and output high-resolution detailed forecasts. The resulting forecasts are not just driven by observed data but also capture the intricate social, demographic and geographic attributes of specific urban regions and mathematical theories of disease propagation over networks. We focus on forecasting the incidence of ILI and evaluate TDEFSI's performance using synthetic and real-world testing datasets at the state and county levels in the USA. The results show that, at the state level, our method achieves comparable/better performance than several state-of-the-art methods. At the county level, TDEFSI outperforms the other methods. The proposed method can be applied to other infectious diseases as well.


title: Universal Average-Case Optimality of Polyak Momentum

abstract: Polyak momentum (PM), also known as the heavy-ball method, is a widely used optimization method that enjoys an asymptotic optimal worst-case complexity on quadratic objectives. However, its remarkable empirical success is not fully explained by this optimality, as the worst-case analysis -- contrary to the average-case -- is not representative of the expected complexity of an algorithm. In this work we establish a novel link between PM and the average-case analysis. Our main contribution is to prove that any optimal average-case method converges in the number of iterations to PM, under mild assumptions. This brings a new perspective on this classical method, showing that PM is asymptotically both worst-case and average-case optimal.


title: FastPET: Near Real-Time PET Reconstruction from Histo-Images Using a Neural Network

abstract: Direct reconstruction of positron emission tomography (PET) data using deep neural networks is a growing field of research. Initial results are promising, but often the networks are complex, memory utilization inefficient, produce relatively small 2D image slices (e.g., 128x128), and low count rate reconstructions are of varying quality. This paper proposes FastPET, a novel direct reconstruction convolutional neural network that is architecturally simple, memory space efficient, works for non-trivial 3D image volumes and is capable of processing a wide spectrum of PET data including low-dose and multi-tracer applications. FastPET uniquely operates on a histo-image (i.e., image-space) representation of the raw data enabling it to reconstruct 3D image volumes 67x faster than Ordered subsets Expectation Maximization (OSEM). We detail the FastPET method trained on whole-body and low-dose whole-body data sets and explore qualitative and quantitative aspects of reconstructed images from clinical and phantom studies. Additionally, we explore the application of FastPET on a neurology data set containing multiple different tracers. The results show that not only are the reconstructions very fast, but the images are high quality and lower noise than iterative reconstructions.


title: Perfect State Transfer on Oriented Graphs

abstract: Quantum walks on undirected graphs have been studied using symmetric matrices, such as the adjacency or Laplacian matrix, and many results about perfect state transfer are known. We extend some of those results to oriented graphs. We also study the phenomena, unique to oriented graphs, of multiple state transfer, where there is a set of vertices such that perfect state transfer occurs between every pair in that set. We give a characterization of multiple state transfer, and a new example of a graph where it occurs.


title: Efficient computation of BSD invariants in genus 2

abstract: Recently, all Birch and Swinnerton-Dyer invariants, except for the order of Sha, have been computed for all curves of genus 2 contained in the L-functions and Modular Forms Database. This report explains the improvements made to the implementation of the algorithm described in arXiv:1711.10409 that were needed to do the computation of the Tamagawa numbers and the real period in reasonable time. We also explain some of the more technical details of the algorithm, and give a brief overview of the methods used to compute the special value of the $L$-function and the regulator.


title: Two-Stage Stochastic Choice Modeling Approach for Electric Vehicle Charging Station Network Design in Urban Communities

abstract: Electric vehicles (EVs) provide a cleaner alternative that not only reduces greenhouse gas emissions but also improves air quality and reduces noise pollution. The consumer market for electrical vehicles is growing very rapidly. Designing a network with adequate capacity and types of public charging stations is a challenge that needs to be addressed to support the current trend in the EV market. In this research, we propose a choice modeling approach embedded in a two-stage stochastic programming model to determine the optimal layout and types of EV supply equipment for a community while considering randomness in demand and drivers' behaviors. Some of the key random data parameters considered in this study are: EV's dwell time at parking {\sv location}, battery's state of charge, distance from home, willingness to walk, drivers' arrival patterns, and traffic on weekdays and weekends. The two-stage model uses the sample average approximation method, which asymptotically converges to an optimal solution. To address the computational challenges for large-scale instances, we propose an outer approximation decomposition algorithm. We conduct extensive computational experiments to quantify the efficacy of the proposed approach. In addition, we present the results and a sensitivity analysis for a case study based on publicly available data sources.


title: Spectroscopy of hot Gamma Doradus and A-F hybrid Kepler candidates close to the hot border of the Delta Scuti instability strip

abstract: If Gamma Dor-type pulsations are driven by the convective blocking mechanism, a convective envelope at a sufficient depth is essential. There are several hot Gamma Dor and hybrid star candidates in which there should not be an adequate convective envelope to excite the{\gamma}Dor-type oscillations. The existence of these hot objects needs an explanation. Therefore, we selected, observed and studied 24 hot{\gamma}Dor and hybrid candidates to investigate their properties. The atmospheric parameters, chemical abundances and vsini values of the candidates were obtained using medium-resolution (R= 46 000)spectra taken with the FIES instrument mounted at the NordicOptical Telescope. We also carried out frequency analyses of theKeplerlong- and short-cadence data to determine the exact pulsation contents. We found only five bona-fide hot{\gamma}Dor and three bona-fide hot hybrid stars in our sample. The other 16 stars were found to benormal{\gamma}Dor,{\delta}Sct, or hybrid variables. No chemical peculiarity was detected in the spectra of the bona-fide hot{\gamma}Dor and hybrid stars. We investigated the interplay between rotation and pulsational modes. We also found that the hot{\gamma}Dor stars havehigherGaialuminosities and larger radii compared to main-sequence A-F stars.


title: Variance Reduced Coordinate Descent with Acceleration: New Method With a Surprising Application to Finite-Sum Problems

abstract: We propose an accelerated version of stochastic variance reduced coordinate descent -- ASVRCD. As other variance reduced coordinate descent methods such as SEGA or SVRCD, our method can deal with problems that include a non-separable and non-smooth regularizer, while accessing a random block of partial derivatives in each iteration only. However, ASVRCD incorporates Nesterov's momentum, which offers favorable iteration complexity guarantees over both SEGA and SVRCD. As a by-product of our theory, we show that a variant of Allen-Zhu (2017) is a specific case of ASVRCD, recovering the optimal oracle complexity for the finite sum objective.


title: Can I Trust You? A User Study of Robot Mediation of a Support Group

abstract: Socially assistive robots have the potential to improve group dynamics when interacting with groups of people in social settings. This work contributes to the understanding of those dynamics through a user study of trust dynamics in the novel context of a robot mediated support group. For this study, a novel framework for robot mediation of a support group was developed and validated. To evaluate interpersonal trust in the multi-party setting, a dyadic trust scale was implemented and found to be uni-factorial, validating it as an appropriate measure of general trust. The results of this study demonstrate a significant increase in average interpersonal trust after the group interaction session, and qualitative post-session interview data report that participants found the interaction helpful and successfully supported and learned from one other. The results of the study validate that a robot-mediated support group can improve trust among strangers and allow them to share and receive support for their academic stress.


title: Object Detection as a Positive-Unlabeled Problem

abstract: As with other deep learning methods, label quality is important for learning modern convolutional object detectors. However, the potentially large number and wide diversity of object instances that can be found in complex image scenes makes constituting complete annotations a challenging task; objects missing annotations can be observed in a variety of popular object detection datasets. These missing annotations can be problematic, as the standard cross-entropy loss employed to train object detection models treats classification as a positive-negative (PN) problem: unlabeled regions are implicitly assumed to be background. As such, any object missing a bounding box results in a confusing learning signal, the effects of which we observe empirically. To remedy this, we propose treating object detection as a positive-unlabeled (PU) problem, which removes the assumption that unlabeled regions must be negative. We demonstrate that our proposed PU classification loss outperforms the standard PN loss on PASCAL VOC and MS COCO across a range of label missingness, as well as on Visual Genome and DeepLesion with full labels.


title: The Einstein condition on nearly K\"ahler six-manifolds

abstract: We review basic facts on the structure of nearly K\"ahler manifolds, focussing in particular on the six-dimensional case. A self-contained proof that nearly K\"ahler six-manifolds are Einstein is given by combining different known results. We finally rephrase the definition of nearly K\"ahler six-manifold in terms of a pair of partial differential equations.


title: Quantum-Enhanced continuous-wave stimulated Raman spectroscopy

abstract: Stimulated Raman spectroscopy has become a powerful tool to study the spatiodynamics of molecular bonds with high sensitivity, resolution and speed. However, sensitivity and speed of state-of-the-art stimulated Raman spectroscopy are currently limited by the shot-noise of the light beam probing the Raman process. Here, we demonstrate an enhancement of the sensitivity of continuous-wave stimulated Raman spectroscopy by reducing the quantum noise of the probing light below the shot-noise limit by means of amplitude squeezed states of light. Probing polymer samples with Raman shifts around 2950 $cm^{-1}$ with squeezed states, we demonstrate a quantum-enhancement of the stimulated Raman signal-to-noise ratio (SNR) of 3.60 dB relative to the shot-noise limited SNR. Our proof-of-concept demonstration of quantum-enhanced Raman spectroscopy paves the way for a new generation of Raman microscopes, where weak Raman transitions can be imaged without the use of markers or an increase in the total optical power.


title: Intra-Horizon Expected Shortfall and Risk Structure in Models with Jumps

abstract: The present article deals with intra-horizon risk in models with jumps. Our general understanding of intra-horizon risk is along the lines of the approach taken in Boudoukh, Richardson, Stanton and Whitelaw (2004), Rossello (2008), Bhattacharyya, Misra and Kodase (2009), Bakshi and Panayotov (2010), and Leippold and Vasiljevi\'c (2019). In particular, we believe that quantifying market risk by strictly relying on point-in-time measures cannot be deemed a satisfactory approach in general. Instead, we argue that complementing this approach by studying measures of risk that capture the magnitude of losses potentially incurred at any time of a trading horizon is necessary when dealing with (m)any financial position(s). To address this issue, we propose an intra-horizon analogue of the expected shortfall for general profit and loss processes and discuss its key properties. Our intra-horizon expected shortfall is well-defined for (m)any popular class(es) of L\'evy processes encountered when modeling market dynamics and constitutes a coherent measure of risk, as introduced in Cheridito, Delbaen and Kupper (2004). On the computational side, we provide a simple method to derive the intra-horizon risk inherent to popular L\'evy dynamics. Our general technique relies on results for maturity-randomized first-passage probabilities and allows for a derivation of diffusion and single jump risk contributions. These theoretical results are complemented with an empirical analysis, where popular L\'evy dynamics are calibrated to S&P 500 index data and an analysis of the resulting intra-horizon risk is presented.


title: Reinforcement Learning Enhanced Quantum-inspired Algorithm for Combinatorial Optimization

abstract: Quantum hardware and quantum-inspired algorithms are becoming increasingly popular for combinatorial optimization. However, these algorithms may require careful hyperparameter tuning for each problem instance. We use a reinforcement learning agent in conjunction with a quantum-inspired algorithm to solve the Ising energy minimization problem, which is equivalent to the Maximum Cut problem. The agent controls the algorithm by tuning one of its parameters with the goal of improving recently seen solutions. We propose a new Rescaled Ranked Reward (R3) method that enables stable single-player version of self-play training that helps the agent to escape local optima. The training on any problem instance can be accelerated by applying transfer learning from an agent trained on randomly generated problems. Our approach allows sampling high-quality solutions to the Ising problem with high probability and outperforms both baseline heuristics and a black-box hyperparameter optimization approach.


title: On a Theorem of Wolff Revisited

abstract: We study $p$-harmonic functions, $ 1 < p\neq 2 < \infty$, in $ \mathbb{R}^{2}_+ = \{ z = x + i y : y > 0, - \infty < x < \infty \} $ and $B( 0, 1 ) = \{ z : |z| < 1 \}$. We first show for fixed $ p$, $1 < p\neq 2 < \infty$, and for all large integers $N\geq N_0$ that there exists $p$-harmonic function, $ V = V ( r e^{i\theta} )$, which is $ 2\pi/N $ periodic in the $ \theta $ variable, and Lipschitz continuous on $ \partial B (0, 1)$ with Lipschitz norm $\leq c N$ on $ \partial B ( 0, 1 )$ satisfying $V(0)=0$ and $ c^{-1} \leq \int_{-\pi}^{\pi} V ( e^{i\theta} ) d \theta \leq c$. In case $2<p<\infty $ we give a more or less explicit example of $V$ and our work is an extension of a result of Wolff on $ \mathbb{R}^{2}_+ $ to $ B (0, 1)$. Using our first result, we extend the work of Wolff on failure of Fatou type theorems for $ \mathbb{R}^{2}_+ $ to $ B (0, 1)$ for $p$-harmonic functions, $1< p\neq 2<\infty$. Finally, we also outline the modifications needed for extending the work of Llorente, Manfredi, and Wu regarding failure of subadditivity of $p$-harmonic measure on $ \partial \mathbb{R}^{2}_+ $ to $\partial B (0, 1)$.


title: Adjusting Image Attributes of Localized Regions with Low-level Dialogue

abstract: Natural Language Image Editing (NLIE) aims to use natural language instructions to edit images. Since novices are inexperienced with image editing techniques, their instructions are often ambiguous and contain high-level abstractions that tend to correspond to complex editing steps to accomplish. Motivated by this inexperience aspect, we aim to smooth the learning curve by teaching the novices to edit images using low-level commanding terminologies. Towards this end, we develop a task-oriented dialogue system to investigate low-level instructions for NLIE. Our system grounds language on the level of edit operations, and suggests options for a user to choose from. Though compelled to express in low-level terms, a user evaluation shows that 25% of users found our system easy-to-use, resonating with our motivation. An analysis shows that users generally adapt to utilizing the proposed low-level language interface. In this study, we identify that object segmentation as the key factor to the user satisfaction. Our work demonstrates the advantages of the low-level, direct language-action mapping approach that can be applied to other problem domains beyond image editing such as audio editing or industrial design.


title: IPBoost -- Non-Convex Boosting via Integer Programming

abstract: Recently non-convex optimization approaches for solving machine learning problems have gained significant attention. In this paper we explore non-convex boosting in classification by means of integer programming and demonstrate real-world practicability of the approach while circumventing shortcomings of convex boosting approaches. We report results that are comparable to or better than the current state-of-the-art.


title: Discretization of the Koch Snowflake Domain with Boundary and Interior Energies

abstract: We study the discretization of a Dirichlet form on the Koch snowflake domain and its boundary with the property that both the interior and the boundary can support positive energy. We compute eigenvalues and eigenfunctions, and demonstrate the localization of high energy eigenfunctions on the boundary via a modification of an argument of Filoche and Mayboroda. H\"older continuity and uniform approximation of eigenfunctions are also discussed.

categories: math.NA cs.NA math-ph math.AP math.MP math.SP


title: Quadratic Optimization with Switching Variables: The Convex Hull for $n = 2$

abstract: We consider quadratic optimization in variables $(x,y)$ where $0\le x\le y$, and $y\in\{0,1\}^n$. Such binary $y$ are commonly refered to as "indicator" or "switching" variables and occur commonly in applications. One approach to such problems is based on representing or approximating the convex hull of the set $\{ (x,xx^T, yy^T) : 0\le x\le y\in\{0,1\}^n\}$. A representation for the case $n=1$ is known and has been widely used. We give an exact representation for the case $n=2$ by starting with a disjunctive representation for the convex hull and then eliminating auxilliary variables and constraints that do not change the projection onto the original variables. An alternative derivation for this representation leads to an appealing conjecture for a simplified representation of the convex hull for $n=2$ when the product term $y_1y_2$ is ignored.


title: On the Mechanical and Thermal Stability of Free-standing Monolayer Amorphous Carbon

abstract: Recently (C.-T. Toh et al., Nature 577, 199 (2020)), the first synthesis of free-standing monolayer amorphous carbon (MAC) was achieved. MAC is a pure carbon structure composed of five, six, seven and eight atom rings randomly distributed. MAC proved to be surprisingly stable and highly fracture resistant. Its electronic properties are similar to boron nitride. In this work, we have investigated the mechanical properties and thermal stability of MAC models using fully-atomistic reactive molecular dynamics simulations. For comparison purposes, the results are contrasted against pristine graphene (PG) models of similar dimensions. Our results show that MAC and PG exhibit distinct mechanical behavior and fracture dynamics patterns. While PG after a critical strain threshold goes directly from elastic to brittle regimes, MAC shows different elastic stages between these two regimes. Remarkably, MAC is thermally stable up to 3600 K, which is close to the PG melting point. These exceptional physical properties make MAC-based materials promising candidates for new technologies, such as flexible electronics.


title: Learning with Out-of-Distribution Data for Audio Classification

abstract: In supervised machine learning, the assumption that training data is labelled correctly is not always satisfied. In this paper, we investigate an instance of labelling error for classification tasks in which the dataset is corrupted with out-of-distribution (OOD) instances: data that does not belong to any of the target classes, but is labelled as such. We show that detecting and relabelling certain OOD instances, rather than discarding them, can have a positive effect on learning. The proposed method uses an auxiliary classifier, trained on data that is known to be in-distribution, for detection and relabelling. The amount of data required for this is shown to be small. Experiments are carried out on the FSDnoisy18k audio dataset, where OOD instances are very prevalent. The proposed method is shown to improve the performance of convolutional neural networks by a significant margin. Comparisons with other noise-robust techniques are similarly encouraging.


title: Controlador LQR y SMC Aplicado a Plataformas Pendulares

abstract: A pendular platform is a robotic structure commonly used in the design of controllers given its nonlinear dynamics; This work presents the modeling, design and implementation of an optimal LQR controller and a Sliding Mode SMC controller applied to two commercial platforms, the Quanser rotary inverted pendulum (RotPen) and the Lego mobile inverted pendulum (NxtWay). The contribution of this work is to present a methodology of implementation of LQR and SMC controllers on pendular platforms, attending the respective restrictions of hardware and software in commercial prototypes. The article presents the behavior of the controller designed on the analytical model compared to its implementation.


title: Learning spatio-temporal representations with temporal squeeze pooling

abstract: In this paper, we propose a new video representation learning method, named Temporal Squeeze (TS) pooling, which can extract the essential movement information from a long sequence of video frames and map it into a set of few images, named Squeezed Images. By embedding the Temporal Squeeze pooling as a layer into off-the-shelf Convolution Neural Networks (CNN), we design a new video classification model, named Temporal Squeeze Network (TeSNet). The resulting Squeezed Images contain the essential movement information from the video frames, corresponding to the optimization of the video classification task. We evaluate our architecture on two video classification benchmarks, and the results achieved are compared to the state-of-the-art.


title: Topological superconductivity in the doped chiral spin liquid on the triangular lattice

abstract: It has long been proposed that doping a chiral spin liquid (CSL) or fractional quantum Hall state can give rise to topological superconductivity. Despite of intensive effort, definitive evidences still remain lacking. We address this problem by studying the $t$-$J$ model supplemented by time-reversal symmetry breaking chiral interaction $J_\chi$ on the triangular lattice using density-matrix renormalization group with a finite concentration $\delta$ of doped holes. It has been established that the undoped, i.e., $\delta$=0, system has a CSL ground state in the parameter region $0.32\le J_\chi/J \le 0.56$. Upon light doping, we find that the ground state of the system is consistent with a Luther-Emery liquid with power-law superconducting and charge-density-wave correlations but short-range spin-spin correlations. In particular, the superconducting correlations, whose pairing symmetry is consistent with $d\pm id$-wave, are dominant at all hole doping concentrations. Our results provide direct evidences that doping the CSL on the triangular lattice can naturally give rise to topological superconductivity.


title: Think Global, Act Local: Relating DNN generalisation and node-level SNR

abstract: The reasons behind good DNN generalisation remain an open question. In this paper we explore the problem by looking at the Signal-to-Noise Ratio of nodes in the network. Starting from information theory principles, it is possible to derive an expression for the SNR of a DNN node output. Using this expression we construct figures-of-merit that quantify how well the weights of a node optimise SNR (or, equivalently, information rate). Applying these figures-of-merit, we give examples indicating that weight sets that promote good SNR performance also exhibit good generalisation. In addition, we are able to identify the qualities of weight sets that exhibit good SNR behaviour and hence promote good generalisation. This leads to a discussion of how these results relate to network training and regularisation. Finally, we identify some ways that these observations can be used in training design.


title: fastai: A Layered API for Deep Learning

abstract: fastai is a deep learning library which provides practitioners with high-level components that can quickly and easily provide state-of-the-art results in standard deep learning domains, and provides researchers with low-level components that can be mixed and matched to build new approaches. It aims to do both things without substantial compromises in ease of use, flexibility, or performance. This is possible thanks to a carefully layered architecture, which expresses common underlying patterns of many deep learning and data processing techniques in terms of decoupled abstractions. These abstractions can be expressed concisely and clearly by leveraging the dynamism of the underlying Python language and the flexibility of the PyTorch library. fastai includes: a new type dispatch system for Python along with a semantic type hierarchy for tensors; a GPU-optimized computer vision library which can be extended in pure Python; an optimizer which refactors out the common functionality of modern optimizers into two basic pieces, allowing optimization algorithms to be implemented in 4-5 lines of code; a novel 2-way callback system that can access any part of the data, model, or optimizer and change it at any point during training; a new data block API; and much more. We have used this library to successfully create a complete deep learning course, which we were able to write more quickly than using previous approaches, and the code was more clear. The library is already in wide use in research, industry, and teaching. NB: This paper covers fastai v2, which is currently in pre-release at http://dev.fast.ai/


title: Two Huge Title and Keyword Generation Corpora of Research Articles

abstract: Recent developments in sequence-to-sequence learning with neural networks have considerably improved the quality of automatically generated text summaries and document keywords, stipulating the need for even bigger training corpora. Metadata of research articles are usually easy to find online and can be used to perform research on various tasks. In this paper, we introduce two huge datasets for text summarization (OAGSX) and keyword generation (OAGKX) research, containing 34 million and 23 million records, respectively. The data were retrieved from the Open Academic Graph which is a network of research profiles and publications. We carefully processed each record and also tried several extractive and abstractive methods of both tasks to create performance baselines for other researchers. We further illustrate the performance of those methods previewing their outputs. In the near future, we would like to apply topic modeling on the two sets to derive subsets of research articles from more specific disciplines.


title: Plasmon Oscillations and de Broglie's Matter Waves Instabilities

abstract: In this research we study the effect of matter-wave instability on electron beam transport with arbitrary degree of degeneracy. Particular class of solutions of the Schr\"{o}dinger-Poisson system is used to model the electron-beam transport at constant speed. It is shown that such electron-beam is described by a coupled driven pseudoforce system solution of which leads to plasmon excitations with dual wave-particle character. The fundamental quantum mechanical de Broglie relation is found to be due to the resonant interaction of particle-like plasmon excitation branch with the electron beam drift. We further obtain a generalized double lengthscale de Broglie wave-particle relation through which various beam-plasmon instability is studied in this model. The quantum charge screening and the chemical potential effects on the matter-wave formation and instabilities are discussed in detail and the well-known Aharonov-Bohm effect is revisited in current quantum hydrodynamic model. Current research may further illuminate the origin of matter-wave in quantum mechanics and lead to clear understanding of novel wave-particle interactions.


title: Computationally efficient univariate filtering for massive data

abstract: The vast availability of large scale, massive and big data has increased the computational cost of data analysis. One such case is the computational cost of the univariate filtering which typically involves fitting many univariate regression models and is essential for numerous variable selection algorithms to reduce the number of predictor variables. The paper manifests how to dramatically reduce that computational cost by employing the score test or the simple Pearson correlation (or the t-test for binary responses). Extensive Monte Carlo simulation studies will demonstrate their advantages and disadvantages compared to the likelihood ratio test and examples with real data will illustrate the performance of the score test and the log-likelihood ratio test under realistic scenarios. Depending on the regression model used, the score test is 30 - 60,000 times faster than the log-likelihood ratio test and produces nearly the same results. Hence this paper strongly recommends to substitute the log-likelihood ratio test with the score test when coping with large scale data, massive data, big data, or even with data whose sample size is in the order of a few tens of thousands or higher.


title: Invariant Risk Minimization Games

abstract: The standard risk minimization paradigm of machine learning is brittle when operating in environments whose test distributions are different from the training distribution due to spurious correlations. Training on data from many environments and finding invariant predictors reduces the effect of spurious features by concentrating models on features that have a causal relationship with the outcome. In this work, we pose such invariant risk minimization as finding the Nash equilibrium of an ensemble game among several environments. By doing so, we develop a simple training algorithm that uses best response dynamics and, in our experiments, yields similar or better empirical accuracy with much lower variance than the challenging bi-level optimization problem of Arjovsky et al. (2019). One key theoretical contribution is showing that the set of Nash equilibria for the proposed game are equivalent to the set of invariant predictors for any finite number of environments, even with nonlinear classifiers and transformations. As a result, our method also retains the generalization guarantees to a large set of environments shown in Arjovsky et al. (2019). The proposed algorithm adds to the collection of successful game-theoretic machine learning algorithms such as generative adversarial networks.


title: Active suspensions of bacteria and passive objects: a model for the near field pair dynamics

abstract: Near field hydrodynamic interactions are essential to determine many important emergent behaviors observed in active suspensions, but have not been successfully modeled so far. In this work we propose an effective model capable of efficiently capturing the essence of the near field hydrodynamic interactions, validated numerically by a pedagogic model system consisting of an E. coli and a spherical tracer. The proposed model effectively captures all the details of near field hydrodynamics through only a tensorial coefficient of resistance, which is fundamentally different from, and thus cannot be replaced by, an effective interaction of conservative nature. In a critical test case that studies the scattering angle of the bacterium-tracer pair dynamics, calculations based on the proposed model reveals a region in parameter space where the bacterium is trapped by the spherical tracer, a phenomenon that is regularly observed in experiments but cannot be explained by any existing model.


title: Adversarial Robustness for Code

abstract: Machine learning and deep learning in particular has been recently used to successfully address many tasks in the domain of code such as finding and fixing bugs, code completion, decompilation, type inference and many others. However, the issue of adversarial robustness of models for code has gone largely unnoticed. In this work, we explore this issue by: (i) instantiating adversarial attacks for code (a domain with discrete and highly structured inputs), (ii) showing that, similar to other domains, neural models for code are vulnerable to adversarial attacks, and (iii) combining existing and novel techniques to improve robustness while preserving high accuracy.


title: Analyzing the Rework Time and Severity of Code Debt: A Case Study Using Technical Debt Catalogs

abstract: This paper presents a case study analyzing Hibernate ecosystem software projects to investigate and demonstrate Code Debt behavior in relation to severity and rework time. The case study carried out revealed that the Code Debt with severity related to impact on software maintenance is the most representative and has the largest rework times to be paid in the Hibernate ecosystem. Besides, it was found that the distributions of rework times of Code Debt for all severities undergo variations in the initial versions of the development.


title: Sigma model on a squashed sphere with a Wess-Zumino term

abstract: A class of two-dimensional sigma models interpolating between $CP^1$ and the $SU(2)$ principal chiral model is discussed. We add the Wess-Zumino-Novikov-Witten term and examine the renormalization group flow of the two coupling constants which characterize the model under consideration. The model flows to the $SU(2)$ WZNW conformal field theory in the IR. There is an ordinary phase in which the model flows from the asymptotically free $CP^1$ model coupled to an extra massless degree of freedom in the UV. At higher loop order we discover that there is also a phase in which the model can flow from non-trivial fixed points in the UV. A non-perturbative confirmation of these extra fixed points would be desirable.


title: Selecting time-series hyperparameters with the artificial jackknife

abstract: This article proposes a generalisation of the delete-$d$ jackknife to solve hyperparameter selection problems for time series. I call it artificial delete-$d$ jackknife to stress that this approach substitutes the classic removal step with a fictitious deletion, wherein observed datapoints are replaced with artificial missing values. This procedure keeps the data order intact and allows plain compatibility with time series. This manuscript justifies the use of this approach asymptotically and shows its finite-sample advantages through simulation studies. Besides, this article describes its real-world advantages by regulating high-dimensional forecasting models for foreign exchange rates.


title: Improving Place Recognition Using Dynamic Object Detection

abstract: We present a novel approach to place recognition well-suited to environments with many dynamic objects--objects that may or may not be present in an agent's subsequent visits. By incorporating an object-detecting preprocessing step, our approach yields high-quality place representations that incorporate object information. Not only does this result in significantly improved place recognition in dynamic environments, it also significantly reduces memory/storage requirements, which may increase the effectiveness of mobile agents with limited resources.