title: Democratizing the Edge: A Pervasive Edge Computing Framework

abstract: The needs of emerging applications, such as augmented and virtual reality, federated machine learning, and autonomous driving, have motivated edge computing--the push of computation capabilities to the edge. Various edge computing architectures have emerged, including multi-access edge computing and edge-cloud, all with the premise of reducing communication latency and augmenting privacy. However, these architectures rely on static and pre-deployed infrastructure, falling short in harnessing the abundant resources at the network's edge. In this paper, we discuss the design of Pervasive Edge Computing (PEC)--a democratized edge computing framework, which enables end-user devices (e.g., smartphones, IoT devices, and vehicles) to dynamically participate in a large-scale computing ecosystem. Our vision of the democratized edge involves the real-time composition of services using available edge resources like data, software, and compute-hardware from multiple stakeholders. We discuss how the novel Named-Data Networking architecture can facilitate service deployment, discovery, invocation, and migration. We also discuss the economic models critical to the adoption of PEC and the outstanding challenges for its full realization.


title: All in the Exponential Family: Bregman Duality in Thermodynamic Variational Inference

abstract: The recently proposed Thermodynamic Variational Objective (TVO) leverages thermodynamic integration to provide a family of variational inference objectives, which both tighten and generalize the ubiquitous Evidence Lower Bound (ELBO). However, the tightness of TVO bounds was not previously known, an expensive grid search was used to choose a "schedule" of intermediate distributions, and model learning suffered with ostensibly tighter bounds. In this work, we propose an exponential family interpretation of the geometric mixture curve underlying the TVO and various path sampling methods, which allows us to characterize the gap in TVO likelihood bounds as a sum of KL divergences. We propose to choose intermediate distributions using equal spacing in the moment parameters of our exponential family, which matches grid search performance and allows the schedule to adaptively update over the course of training. Finally, we derive a doubly reparameterized gradient estimator which improves model learning and allows the TVO to benefit from more refined bounds. To further contextualize our contributions, we provide a unified framework for understanding thermodynamic integration and the TVO using Taylor series remainders.


title: Object Goal Navigation using Goal-Oriented Semantic Exploration

abstract: This work studies the problem of object goal navigation which involves navigating to an instance of the given object category in unseen environments. End-to-end learning-based navigation methods struggle at this task as they are ineffective at exploration and long-term planning. We propose a modular system called, `Goal-Oriented Semantic Exploration' which builds an episodic semantic map and uses it to explore the environment efficiently based on the goal object category. Empirical results in visually realistic simulation environments show that the proposed model outperforms a wide range of baselines including end-to-end learning-based methods as well as modular map-based methods and led to the winning entry of the CVPR-2020 Habitat ObjectNav Challenge. Ablation analysis indicates that the proposed model learns semantic priors of the relative arrangement of objects in a scene, and uses them to explore efficiently. Domain-agnostic module design allow us to transfer our model to a mobile robot platform and achieve similar performance for object goal navigation in the real-world.


title: Measuring Robustness to Natural Distribution Shifts in Image Classification

abstract: We study how robust current ImageNet models are to distribution shifts arising from natural variations in datasets. Most research on robustness focuses on synthetic image perturbations (noise, simulated weather artifacts, adversarial examples, etc.), which leaves open how robustness on synthetic distribution shift relates to distribution shift arising in real data. Informed by an evaluation of 204 ImageNet models in 213 different test conditions, we find that there is often little to no transfer of robustness from current synthetic to natural distribution shift. Moreover, most current techniques provide no robustness to the natural distribution shifts in our testbed. The main exception is training on larger and more diverse datasets, which in multiple cases increases robustness, but is still far from closing the performance gaps. Our results indicate that distribution shifts arising in real data are currently an open research problem. We provide our testbed and data as a resource for future work at https://modestyachts.github.io/imagenet-testbed/ .


title: On the Dimensions of the Realization Spaces of Polytopes

abstract: Robertson (1988) suggested a model for the realization space of a convex d-dimensional polytope and an approach via the implicit function theorem, which -- in the case of a full rank Jacobian -- proves that the realization space is a manifold of dimension NG(P):=d(f_0+f_{d-1})-f_{0,d-1}, which is the natural guess for the dimension given by the number of variables minus the number of quadratic equations that are used in the definition of the realization space. While this indeed holds for many natural classes of polytopes (including simple and simplicial polytopes, as well as all polytopes of dimension at most 3),and Robertson claimed this to be true for all polytopes, Mnev's (1986/1988) Universality Theorem implies that it is not true in general: Indeed, (1) the centered realization space is not a smoothly embedded manifold in general, and (2) it does not have the dimension NG(P) in general. In this paper we develop Jacobian criteria for the analysis of realization spaces. From these we get easily that for various large and natural classes of polytopes the realization spaces are indeed manifolds, whose dimensions are given by NG(P). However, we also identify the smallest polytopes where the dimension count (2) and thus Robertson's claim fails, among them the bipyramid over a triangular prism. For the property (1), we analyze the classical 24-cell: We show that the realization space has at least the dimension 48, and it has points where it is a manifold of this dimension, but it is not smoothly embedded as a manifold everywhere.


title: Symbol Alphabets from Plabic Graphs

abstract: Symbol alphabets of n-particle amplitudes in N=4 super-Yang-Mills theory are known to contain certain cluster variables of Gr(4,n) as well as certain algebraic functions of cluster variables. In this paper we suggest an algorithm for computing these symbol alphabets from plabic graphs by solving matrix equations of the form C Z = 0 to associate functions on Gr(m,n) to parameterizations of certain cells of Gr(k,n) indexed by plabic graphs. For m=4 and n=8 we show that this association precisely reproduces the 18 algebraic symbol letters of the two-loop NMHV eight-point amplitude from four plabic graphs.


title: New N=2 superconformal field theories from S-folds

abstract: We study the four-dimensional N=2 superconformal field theories that describe D3-branes probing the recently constructed N=2 S-folds in F-theory. We introduce a novel, infinite class of superconformal field theories related to S-fold theories via partial Higgsing. We determine several properties of both the S-fold models and this new class of theories, including their central charges, Coulomb branch spectrum, and moduli spaces of vacua, by bringing to bear an array of field-theoretical techniques, to wit, torus-compactifications of six-dimensional N=(1,0) theories, class S technology, and the SCFT/VOA correspondence.


title: Traffic Performance Score for Measuring the Impact of COVID-19 on Urban Mobility

abstract: Measuring traffic performance is critical for public agencies who manage traffic and individuals who plan trips, especially when special events happen. The COVID-19 pandemic has significantly influenced almost every aspect of daily life, including urban traffic patterns. Thus, it is important to measure the impact of COVID-19 on transportation to further guide agencies and residents to properly respond to changes in traffic patterns. However, most existing traffic performance metrics incorporate only a single traffic parameter and measure only the performance of individual corridors. To overcome these challenges, in this study, a Traffic Performance Score (TPS) is proposed that incorporates multiple parameters for measuring network-wide traffic performance. An interactive web-based TPS platform that provides real-time and historical spatial-temporal traffic performance analysis is developed by the STAR Lab at the University of Washington. Based on data from this platform, this study analyzes the impact of COVID-19 on different road segments and the traffic network as a whole. Considering this pandemic has greatly reshaped social and economic operations, this study also evaluates how COVID-19 is changing the urban mobility from both travel demand and driving behavior perspectives.


title: Group Ensemble: Learning an Ensemble of ConvNets in a single ConvNet

abstract: Ensemble learning is a general technique to improve accuracy in machine learning. However, the heavy computation of a ConvNets ensemble limits its usage in deep learning. In this paper, we present Group Ensemble Network (GENet), an architecture incorporating an ensemble of ConvNets in a single ConvNet. Through a shared-base and multi-head structure, GENet is divided into several groups to make explicit ensemble learning possible in a single ConvNet. Owing to group convolution and the shared-base, GENet can fully leverage the advantage of explicit ensemble learning while retaining the same computation as a single ConvNet. Additionally, we present Group Averaging, Group Wagging and Group Boosting as three different strategies to aggregate these ensemble members. Finally, GENet outperforms larger single networks, standard ensembles of smaller networks, and other recent state-of-the-art methods on CIFAR and ImageNet. Specifically, group ensemble reduces the top-1 error by 1.83% for ResNeXt-50 on ImageNet. We also demonstrate its effectiveness on action recognition and object detection tasks.


title: Missing in Axion: where are XENON1T's big black holes?

abstract: We pioneer the black hole mass gap as a powerful new tool for constraining new particles. A new particle that couples to the Standard Model---such as an axion---acts as an additional source of loss in the cores of population-III stars, suppressing mass lost due to winds and quenching the pair-instability. This results in heavier astrophysical black holes. As an example, using stellar simulations we show that the solar axion explanation of the recent XENON1T excess implies astrophysical black holes of ~ 56 MS, squarely within the black hole mass gap predicted by the Standard Model.

categories: hep-ph astro-ph.CO astro-ph.GA astro-ph.HE gr-qc


title: Revisiting Comparative Performance of DNS Resolvers in the IPv6 and ECS Era

abstract: This paper revisits the issue of the performance of DNS resolution services available to Internet users. While several prior studies addressed this important issue, significant developments, namely, the IPv6 finally getting traction and the adoption of the ECS extension to DNS by major DNS resolution services, warrant a reassessment under these new realities. We find that DNS resolution services differ drastically -- by an order of magnitude in some locations -- in their query response time. We also find established resolvers (Google DNS and OpenDNS) to lag far behind relative newcomers (Cloudflair and Quad9) in terms of DNS latency, and trace the cause to drastically lower cache hit rates, which we further trace to less cache sharing within the resolver platform. In addition, we find that public resolvers have largely closed the gap with ISP resolvers in the quality of CDNs' client-to-edge-server mappings as measured by latency. Finally, in most locations, we observe IPv6 penalty in the latency of client-to-CDN-edge-server mappings produced by the resolvers. Moreover, this penalty, while often significant, still does not rise above typical thresholds employed by the Happy Eyeballs algorithm for preferring IPv4 communication. resolvers. Thus, dual-stacked clients in these locations may experience suboptimal performance.


title: Methanol formation chemistry with revised reactions scheme

abstract: The aim of the presented work is to analyze the impact of experimentally evaluated reactions of hydrogen abstraction on surfaces of interstellar grains on the chemical evolution of methanol and its precursors on grains and in the gas phase under conditions of cold dark cloud and during the collapse of the translucent cloud into the dark cloud. Analysis of simulation results shows that those reactions are highly efficient destruction channels for HCO and H2CO on grain surfaces, and significantly impact the abundances of almost all molecules participating in the formation of CH3OH. Next, in models with those reactions maximum abundances of methanol in gas and on grain surface decrease by more than 2-3 orders of magnitude in comparison to models without surface abstraction reactions of hydrogen. Finally, we study the impact of binding energies of CH2OH and CH3O radicals on methanol chemistry.


title: Swapping Autoencoder for Deep Image Manipulation

abstract: Deep generative models have become increasingly effective at producing realistic images from randomly sampled seeds, but using such models for controllable manipulation of existing images remains challenging. We propose the Swapping Autoencoder, a deep model designed specifically for image manipulation, rather than random sampling. The key idea is to encode an image with two independent components and enforce that any swapped combination maps to a realistic image. In particular, we encourage the components to represent structure and texture, by enforcing one component to encode co-occurrent patch statistics across different parts of an image. As our method is trained with an encoder, finding the latent codes for a new input image becomes trivial, rather than cumbersome. As a result, it can be used to manipulate real input images in various ways, including texture swapping, local and global editing, and latent code vector arithmetic. Experiments on multiple datasets show that our model produces better results and is substantially more efficient compared to recent generative models.


title: Scale-invariant dynamics in the Solar System

abstract: The covariant scale-invariant dynamics (SID) theory has recently been proposed as a possible explanation for the observed dynamical discrepancies in galaxies (Maeder & Gueorguiev 2020). SID implies that these discrepancies $-$ commonly attributed to dark matter $-$ arise instead from a non-standard velocity-dependent force that causes two-body near-Keplerian orbits to expand. We show that the predicted expansion of the Earth-Moon orbit is incompatible with lunar laser ranging data at $>200\sigma$. Moreover, SID predicts that the gravitating mass of any object was much smaller in the past. If true, a low-mass red giant star must be significantly older than in standard theory. This would make it much older than the conventional age of the Universe, which however is expected to be similarly old in SID. Moreover, it is not completely clear whether SID truly contains new physics beyond General Relativity, with several previous works arguing that the extra degree of freedom is purely mathematical. We conclude that the SID model is falsified at high significance by observations across a range of scales, even if it is theoretically well formulated.


title: Knowledge-Aware Language Model Pretraining

abstract: How much knowledge do pretrained language models hold? Recent research observed that pretrained transformers are adept at modeling semantics but it is unclear to what degree they grasp human knowledge, or how to ensure they do so. In this paper we incorporate knowledge-awareness in language model pretraining without changing the transformer architecture, inserting explicit knowledge layers, or adding external storage of semantic information. Rather, we simply signal the existence of entities to the input of the transformer in pretraining, with an entity-extended tokenizer; and at the output, with an additional entity prediction task. Our experiments show that solely by adding these entity signals in pretraining, significantly more knowledge is packed into the transformer parameters: we observe improved language modeling accuracy, factual correctness in LAMA knowledge probing tasks, and semantics in the hidden representations through edge probing.We also show that our knowledge-aware language model (KALM) can serve as a drop-in replacement for GPT-2 models, significantly improving downstream tasks like zero-shot question-answering with no task-related training.


title: Throughput Analysis for Virtual MIMO WSNs over Measured MIMO Channels

abstract: A recently conducted indoor-to-outdoor measurement campaign for investigating the propagation characteristics of an $8 \times 8$ virtual multiple-input-multiple-output (MIMO) based wireless sensor network (WSN) is presented in this paper. The campaign is conducted in an instrumentation room devoid of windows, but filled with different noisy electrical and measuring units. The channel impulse responses are reported when a 20 MHz wide signal is transmitted at 2.53 GHz. Measurements are collected for 15 different spatial combinations of the transmit antennas. After analyzing the collected data, system capacity and achievable transmission rates are calculated for each measurement scenario. Using these values, we examined the best configuration for positioning the sensors that can maximize overall network throughput. Results demonstrated that distributing sensors on all 4 walls of the room can achieve the highest possible information rate.


title: Constructing Basis Path Set by Eliminating Path Dependency

abstract: The way the basis path set works in neural network remains mysterious, and the generalization of newly appeared G-SGD algorithm to more practical network is hindered. The Basis Path Set Searching problem is formulated from the perspective of graph theory, to find the basis path set in a regular complicated neural network. Our paper aims to discover the underlying cause of path dependency between two independent substructures. Algorithm DEAH is designed to solve the Basis Path Set Searching problem by eliminating such path dependency. The path subdivision chain is proposed to effectively eliminate the path dependency inside the chain and between chains. The theoretical proofs and analysis of polynomial time complexity are presented. The paper therefore provides one methodology to find the basis path set in a more general neural network, which offers theoretical and algorithmic support for the application of G-SGD algorithm in more practical scenarios.


title: Evolution of spin freezing transition and structural, magnetic phase diagram of Dy$_{2-x}$La$_x$Zr$_2$O$_7$; $x$ = 0-2.0

abstract: Dy$_{2}$Zr$_{2}$O$_{7}$ a disordered pyrochlore system, exhibits the spin ice freezing under the application of magnetic field. Our studies suggest the stabilization of pyrochlore phase in Dy$_{2-x}$La$_{x}$Zr$_{2}$O$_{7}$ with the substitution of nonmagnetic La, along with the biphasic mixture for the intermediate compositions. We observed that the higher La compositions (1.5 $\leq$ x $\leq$ 1.9), show spin freezing (T $\sim$ 17 K) similar to the field induced spin ice freezing for low La compositions (0 $\leq$ x $\leq$ 0.5), and the well known spin ice systems Dy$_{2}$Ti$_{2}$O$_{7}$ and Ho$_{2}$Ti$_{2}$O$_{7}$. The low temperature magnetic state for higher La compositions (1.5 $\leq$ x $\leq$ 1.9) culminates into spin glass state below 6 K. The Cole-Cole plot and Casimir-du Pr$\acute{e}$ fit shows narrow distribution of spin relaxation time in these compounds.


title: LSTM and GPT-2 Synthetic Speech Transfer Learning for Speaker Recognition to Overcome Data Scarcity

abstract: In speech recognition problems, data scarcity often poses an issue due to the willingness of humans to provide large amounts of data for learning and classification. In this work, we take a set of 5 spoken Harvard sentences from 7 subjects and consider their MFCC attributes. Using character level LSTMs (supervised learning) and OpenAI's attention-based GPT-2 models, synthetic MFCCs are generated by learning from the data provided on a per-subject basis. A neural network is trained to classify the data against a large dataset of Flickr8k speakers and is then compared to a transfer learning network performing the same task but with an initial weight distribution dictated by learning from the synthetic data generated by the two models. The best result for all of the 7 subjects were networks that had been exposed to synthetic data, the model pre-trained with LSTM-produced data achieved the best result 3 times and the GPT-2 equivalent 5 times (since one subject had their best result from both models at a draw). Through these results, we argue that speaker classification can be improved by utilising a small amount of user data but with exposure to synthetically-generated MFCCs which then allow the networks to achieve near maximum classification scores.


title: On O+ ion heating by BBELF waves at low altitude: Test particle simulations

abstract: We investigate mechanisms of wave-particle heating of ionospheric O$^{+}$ ions resulting from broadband extremely low frequency (BBELF) waves using numerical test particle simulations that take into account ion-neutral collisions, in order to explain observations from the Enhanced Polar Outflow Probe (e-POP) satellite at low altitudes ($\sim$400 km)[Shen et al., 2018]. We argue that in order to reproduce ion temperatures observed at e-POP altitudes, the most effective ion heating mechanism is through cyclotron acceleration by short-scale electrostatic ion cyclotron (EIC) waves with perpendicular wavelengths $\lambda_{\perp}$$\leq$200 m. The interplay between finite perpendicular wavelengths, wave amplitudes, and ion-neutral collision frequencies collectively determine the ionospheric ion heating limit, which begins to decrease sharply with decreasing altitude below approximately 500 km, where the ratio $\frac{\nu_c}{f_{ci}}$ becomes larger than 10$^{-3}$, $\nu_c$ and $f_{ci}$ denoting the O$^+$--O collision frequency and ion cyclotron frequency. We derive, both numerically and analytically, the ion gyroradius limit from heating by an EIC wave at half the cyclotron frequency. The limit is 0.28$\lambda_{\perp}$. The ion gyroradius limit from an EIC wave can be surpassed either through adding waves with different $\lambda_{\perp}$, or through stochastic "breakout", meaning ions diffuse in energy beyond the gyroradius limit due to stochastic heating from large-amplitude waves. Our two-dimensional simulations indicate that small-scale ($<$1 km) Alfv\'en waves cannot account for the observed ion heating through trapping or stochastic heating.


title: Hubbard ladders at small $U$ revisited

abstract: We re-examine the zero temperature phase diagram of the two-leg Hubbard ladder in the small $U$ limit, both analytically and using density-matrix renormalization group (DMRG). We find a ubiquitous Luther-Emery phase, but with a crossover in behavior at a characteristic interaction strength, $U^\star$; for $U \gtrsim U^\star$, there is a single emergent correlation length $\log[\xi] \sim 1/U$, characterizing the gapped modes of the system, but for $U\lesssim U^\star$ there is a hierarchy of length scales, differing parametrically in powers of $U$, reflecting a two-step renormalization group flow to the ultimate fixed point. Finally, to illustrate the versatility of the approach developed here, we sketch its implications for a half-filled triangular lattice Hubbard model on a cylinder, and find results in conflict with inferences concerning the small $U$ phase from recent DMRG studies of the same problem.


title: Implementing a Fast Unbounded Quantum Fanout Gate Using Power-Law Interactions

abstract: The standard circuit model for quantum computation presumes the ability to directly perform gates between arbitrary pairs of qubits, which is unlikely to be practical for large-scale experiments. Power-law interactions with strength decaying as $1/r^\alpha$ in the distance $r$ provide an experimentally realizable resource for information processing, whilst still retaining long-range connectivity. We leverage the power of these interactions to implement a fast quantum fanout gate with an arbitrary number of targets. Our implementation allows the quantum Fourier transform (QFT) and Shor's algorithm to be performed on a $D$-dimensional lattice in time logarithmic in the number of qubits for interactions with $\alpha \le D$. As a corollary, we show that power-law systems with $\alpha \le D$ are difficult to simulate classically even for short times, under a standard assumption that factoring is classically intractable. Complementarily, we develop a new technique to give a general lower bound, linear in the size of the system, on the time required to implement the QFT and the fanout gate in systems that are constrained by a linear light cone. This allows us to prove an asymptotically tighter lower bound for long-range systems than is possible with previously available techniques.


title: Inverse stellar population age gradients of post-starburst galaxies at z=0.8 with LEGA-C

abstract: We use deep, spatially resolved spectroscopy from the LEGA-C Survey to study radial variations in the stellar population of 17 spectroscopically-selected post-starburst (PSB) galaxies. We use spectral fitting to measure two Lick indices, $H{\delta}_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $H{\delta}_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $H{\delta}_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ \r{A}/$R_e$), mainly due to seeing convolution. A two-SSP model suggests intrinsic gradients are as strong as observed in local PSB galaxies ($\approx -0.8$ \r{A}$/R_e$). We interpret these results in terms of inside-out growth (for the bulk of the quiescent population) vs star formation occurring last in the centre (for PSB galaxies). At $z\approx0.8$, central starbursts are often the result of gas-rich mergers, as evidenced by the high fraction of PSB galaxies with disturbed morphologies and tidal features (40%). Our results provide additional evidence for multiple paths to quiescence: a standard path, associated with inside-out disc formation and with gradually decreasing star-formation activity, without fundamental structural transformation, and a fast path, associated with centrally-concentrated starbursts, leaving an inverse age gradient and smaller half-light radius.


title: Precision tests of fundamental physics with $\eta$ and $\eta^\prime$ mesons

abstract: Decays of the neutral and long-lived $\eta$ and $\eta'$ mesons provide a unique, flavor-conserving laboratory to test low-energy Quantum Chromodynamics and search for new physics beyond the Standard Model. They have drawn world-wide attention in recent years and have inspired broad experimental programs in different high-intensity-frontier centers. New experimental data will offer critical inputs to precisely determine the light quark mass ratios, $\eta$-$\eta'$ mixing parameters, and hadronic contributions to the anomalous magnetic moment of the muon. At the same time, it will provide a sensitive probe to test potential new physics. This includes searches for hidden photons, light Higgs scalars, and axion-like particles that are complementary to worldwide efforts to detect new light particles below the GeV mass scale, as well as tests of discrete symmetry violation. In this review, we give an update on theoretical developments, discuss the experimental opportunities, and identify future research needed in this field.


title: Compact objects in general relativity: From Buchdahl stars to quasiblack holes

abstract: A Buchdahl star is a highly compact star for which the boundary radius $R$ obeys $R=\frac98 r_+$, where $r_+$ is the gravitational radius of the star itself. A quasiblack hole is a maximum compact star, or more generically a maximum compact object, for which the boundary radius $R$ obeys $R=r_+$. Quasiblack holes are objects on the verge of becoming black holes. Continued gravitational collapse ends in black holes and has to be handled with the Oppenheimer-Snyder formalism. Quasistatic contraction ends in a quasiblack hole and should be treated with appropriate techniques. Quasiblack holes, not black holes, are the real descendants of Mitchell and Laplace dark stars. Quasiblack holes have many interesting properties. We develop the concept of a quasiblack hole, give several examples of such an object, define what it is, draw its Carter-Penrose diagram, study its pressure properties, obtain its mass formula, derive the entropy of a nonextremal quasiblack hole, and through an extremal quasiblack hole give a solution to the puzzling entropy of extremal black holes.


title: Non-Analytic Non-Equilibrium Field Theory: Stochastic Reheating of the Ising Model

abstract: Many-body non-equilibrium steady states can be described by a Landau-Ginzburg theory if one allows non-analytic terms in the potential. We substantiate this claim by working out the case of the Ising magnet in contact with a thermal bath and undergoing stochastic reheating: It is reset to a paramagnet at random times. By a combination of stochastic field theory and Monte Carlo simulations, we unveil how the usual $\varphi^4$ potential is deformed by non-analytic operators of intrinsic non-equilibrium nature. We demonstrate their infrared relevance at low temperatures by a renormalization-group analysis of the non-equilibrium steady state. The equilibrium ferromagnetic fixed point is thus destabilized by stochastic reheating, and we identify the new non-equilibrium fixed point.


title: Long range dark matter self-interactions and plasma instabilities

abstract: So far, the observed effects of dark matter are compatible with it having purely gravitational interactions. However, in many models, dark matter has additional interactions with itself, with the Standard Model, and/or with additional hidden sector states. In this paper, we discuss models in which dark matter interacts through a light vector mediator, giving rise to a long-ranged force between dark matter particles. Coherent scattering effects through this force can lead to the exponential growth of small perturbations, in analogy to electromagnetic plasma instabilities. These instabilities can be significant at couplings many orders of magnitude below those for which the usual particle-by-particle constraints on dark matter self-interactions apply. While this possibility has been noted in the literature, we provide the first systematic study of such instabilities, including the case where the mediator has finite mass. The latter is relevant for models of kinetically mixed `dark photon' mediators, which represent an important target for proposed dark matter detection experiments. Our analyses are of the growth of small perturbations, so do not immediately provide observational constraints on dark matter models - however, they do motivate further study of large regions of parameter space.


title: Gauge-Symmetry Protection Using Single-Body Terms

abstract: Quantum-simulator hardware promises new insights into problems from particle and nuclear physics. A major challenge is to reproduce gauge invariance, as violations of this quintessential property of lattice gauge theories can have dramatic consequences, e.g., the generation of a photon mass in quantum electrodynamics. Here, we introduce an experimentally friendly method to protect gauge invariance in $\mathrm{U}(1)$ lattice gauge theories against coherent errors in a controllable way. Our method employs only single-body energy-penalty terms, thus enabling practical implementations. As we derive analytically, some sets of penalty coefficients render undesired gauge sectors inaccessible by unitary dynamics for exponentially long times, and, for few-body error terms, with resources independent of system size. These findings constitute an exponential improvement over previously known results from energy-gap protection or perturbative treatments. In our method, the gauge-invariant subspace is protected by an emergent global symmetry, meaning it can be immediately applied to other symmetries. In our numerical benchmarks for continuous-time and digital quantum simulations, gauge protection holds for all calculated evolution times (up to $t>10^{10}/J$ for continuous time, with $J$ the relevant energy scale). Crucially, our gauge-protection technique is simpler to realize than the associated ideal gauge theory, and can thus be readily implemented in current ultracold-atom analog simulators as well as digital noisy intermediate scale quantum (NISQ) devices.

categories: quant-ph cond-mat.quant-gas cond-mat.str-el hep-ex hep-lat


title: Constraining the Charge-Sign and Rigidity-Dependence of Solar Modulation

abstract: Our ability to identify the sources of cosmic rays and understand how these particles propagate through the interstellar medium is hindered by the combined effects of the solar wind and its embedded magnetic field, collectively known as solar modulation. In this paper, we build upon our previous work to model and constrain the effects of solar modulation on the cosmic-ray spectrum, using data from AMS-02 and BESS Polar II collected between 2007 and 2012, during which the heliospheric magnetic field was in a state of negative polarity. Our model uses measurements of the heliospheric magnetic field and the tilt angle of the heliospheric current sheet to accurately predict the effects of solar modulation as a function of time, charge, and rigidity. By incorporating data from a period of negative polarity, we have been able to robustly observe and constrain the charge-dependent effects of solar modulation.

categories: astro-ph.HE astro-ph.SR hep-ph physics.space-ph


title: The $F$-Symbols for Transparent Haagerup-Izumi Categories with $G = \mathbb{Z}_{2n+1}$

abstract: A fusion category is called transparent if the associator involving any invertible object is the identity map. For the Haagerup-Izumi fusion rings with $G = \mathbb{Z}_{2n+1}$ (the $\mathbb{Z}_3$ case is the Haagerup fusion ring with six simple objects), the transparent ansatz reduces the number of independent $F$-symbols from order $\mathcal{O}(n^6)$ to $\mathcal{O}(n^2)$, rendering the pentagon identity practically solvable. Transparent Haagerup-Izumi fusion categories are thereby constructively classified up to $G = \mathbb{Z}_9$, recovering all known Haagerup-Izumi fusion categories to this order, and producing new ones. Transparent Haagerup-Izumi fusion categories additionally satisfying $S_4$ tetrahedral invariance are further classified up to $G = \mathbb{Z}_{15}$, and the explicit $F$-symbols for the unitary ones, including the Haagerup $\mathcal{H}_3$ fusion category, are compactly presented. The $F$-symbols for the Haagerup $\mathcal{H}_2$ fusion category are also presented. Going beyond, the transparent ansatz offers a viable course towards constructing novel fusion categories for new fusion rings.

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


title: X-Ray Spectroscopy of Newly Identified ULXs Associated With M87's Globular Cluster Population

abstract: We have identified seven ultraluminous X-ray sources (ULXs) which are coincident with globular cluster candidates (GCs) associated with M87. ULXs in the old GC environment represent a new population of ULXs, and ones likely to be black holes. In this study we perform detailed X-ray spectroscopic followup to seven GC ULXs across a wealth of archival Chandra observations and long time baseline of 16 years. This study brings the total known sample of GC ULXs to 17. Two of these sources show variability in their X-ray luminosity of an order of magnitude over many years, and one of these sources shows intra-observational variability on the scale of hours. While the X-ray spectra of the majority of globular cluster ULXs are best fit by single component models, one of the sources studied in this paper is the second GC ULX to be best fit by a two component model. We compare this new sample of GC ULXs to the previously studied sample, and compare the X-ray and optical properties counterparts across the samples. We find that the clusters that host ULXs in M87 have metallicities from g-z=1.01 to g-z=1.70. The best-fit power-law indices of the X-ray spectra range from Gamma=1.37-2.21, and the best fit inner black-body disk temperatures range from kT=0.56-1.90 keV.


title: Systematics of Type IIA moduli stabilisation

abstract: We analyse the flux-induced scalar potential for type IIA orientifolds in the presence of $p$-form, geometric and non-geometric fluxes. Just like in the Calabi-Yau case, the potential presents a bilinear structure, with a factorised dependence on axions and saxions. This feature allows one to perform a systematic search for vacua, which we implement for the case of geometric backgrounds. Guided by stability criteria, we consider configurations with a particular on-shell F-term pattern, for which we derive a no-go result for de Sitter extrema. We classify branches of supersymmetric and non-supersymmetric vacua, and argue that the latter are perturbatively stable for a large subset of them. Our solutions reproduce and generalise previous results in the literature, obtained either from the 4d or 10d viewpoint.


title: The effective field theory of low scale see-saw at colliders

abstract: We study the Standard Model effective field theory ($\nu$SMEFT) extended with operators involving right-handed neutrinos, focussing on the regime where the right-handed neutrinos decay promptly on collider scales to a photon and a Standard Model neutrino. This scenario arises naturally for right-handed neutrinos with masses of the order $m_N \sim 0.1 \dots 10\, \text{GeV}$. We limit the relevant dimension-six operator coefficients using LEP and LHC searches with photons and missing energy in the final state as well as pion and tau decays. While bounds on new physics contributions are generally in the TeV scale for order one operator coefficients, some coefficients, however, remain very poorly constrained or even entirely evade bounds from current data. Consequently, we identify such weakly constrained scenarios and propose new searches for rare top and tau decays involving photons to probe potential new physics in the $\nu$SMEFT parameter space. Our analysis highlights the importance of performing dedicated searches for new rare tau and top decays.


title: Sliced Iterative Normalizing Flows

abstract: We develop an iterative (greedy) deep learning (DL) algorithm which is able to transform an arbitrary probability distribution function (PDF) into the target PDF. The model is based on iterative Optimal Transport of a series of 1D slices, matching on each slice the marginal PDF to the target. The axes of the orthogonal slices are chosen to maximize the PDF difference using Wasserstein distance at each iteration, which enables the algorithm to scale well to high dimensions. As special cases of this algorithm, we introduce two sliced iterative Normalizing Flow (SINF) models, which map from the data to the latent space (GIS) and vice versa (SIG). We show that SIG is able to generate high quality samples of image datasets, which match the GAN benchmarks, while GIS obtains competitive results on density estimation tasks compared to the density trained NFs, and is more stable, faster, and achieves higher $p(x)$ when trained on small training sets. SINF approach deviates significantly from the current DL paradigm, as it is greedy and does not use concepts such as mini-batching, stochastic gradient descent and gradient back-propagation through deep layers.


title: x-cut Cosmic Shear: Optimally Removing Sensitivity to Baryonic and Nonlinear Physics with an Application to the Dark Energy Survey Year 1 Shear Data

abstract: We present a new method, called $x$-cut cosmic shear, which optimally removes sensitivity to poorly modeled scales from the two-point cosmic shear signal. We show that the $x$-cut cosmic shear covariance matrix can be computed from the correlation function covariance matrix in a few minutes, enabling a likelihood analysis at virtually no additional computational cost. Further we show how to generalize $x$-cut cosmic shear to galaxy-galaxy lensing. Performing an $x$-cut cosmic shear analysis of the Dark Energy Survey Year 1 (DESY1) shear data, we reduce the error on $S_8= \sigma_8 (\Omega_m / 0.3) ^ {0.5}$ by $32 \%$ relative to a correlation function analysis with the same priors and angular scale cut criterion, while showing our constraints are robust to different baryonic feedback models. Largely driven by information at small angular scales, our result, $S_8= 0.734 \pm 0.026$, yields a $2.6 \sigma$ tension with the Planck Legacy analysis of the cosmic microwave background. As well as alleviating baryonic modelling uncertainties, our method can be used to optimally constrain a large number of theories of modified gravity where computational limitations make it infeasible to model the power spectrum down to extremely small scales. The key parts of our code are made publicly available.


title: Heavy $Z^\prime$ Bosons in the Secluded $U(1)^\prime$ Model at Hadron Colliders

abstract: We study $Z'$ phenomenology at hadron colliders in an $U(1)'$ extended MSSM. We choose a $U(1)'$ model with a secluded sector, where the tension between the electroweak scale and developing a large enough mass for $Z'$ is resolved by incorporating three additional singlet superfields into the model. We perform a detailed analysis of the production, followed by decays, including into supersymmetric particles, of a $Z'$ boson with mass between 4 and 5.2 TeV, with particular emphasis on its possible discovery. We select three different scenarios consistent with the latest available experimental data and relic density constraints, and concentrate on final signals with two leptons, four leptons and six leptons. Including the SM background from processes with two, three or four vector bosons, we show the likelihood of observing a $Z^\prime$ boson is not promising for the HL-LHC at 14 TeV. While at 27 and 100 TeV, the situation is more optimistic, and we devise specific benchmark scenarios which could be observed.


title: Dispelling Myths on Superposition Attacks: Formal Security Model and Attack Analyses

abstract: It is of folkloric belief that the security of classical cryptographic protocols is automatically broken if the Adversary is allowed to perform superposition queries and the honest players forced to perform actions coherently on quantum states. Another widely held intuition is that enforcing measurements on the exchanged messages is enough to protect protocols from these attacks. However, the reality is much more complex. Security models dealing with superposition attacks only consider unconditional security. Conversely, security models considering computational security assume that all supposedly classical messages are measured, which forbids by construction the analysis of superposition attacks. Boneh and Zhandry have started to study the quantum computational security for classical primitives in their seminal work at Crypto'13, but only in the single-party setting. To the best of our knowledge, an equivalent model in the multiparty setting is still missing. In this work, we propose the first computational security model considering superposition attacks for multiparty protocols. We show that our new security model is satisfiable by proving the security of the well-known One-Time-Pad protocol and give an attack on a variant of the equally reputable Yao Protocol for Secure Two-Party Computations. The post-mortem of this attack reveals the precise points of failure, yielding highly counter-intuitive results: Adding extra classical communication, which is harmless for classical security, can make the protocol become subject to superposition attacks. We use this newly imparted knowledge to construct the first concrete protocol for Secure Two-Party Computation that is resistant to superposition attacks. Our results show that there is no straightforward answer to provide for either the vulnerabilities of classical protocols to superposition attacks or the adapted countermeasures.


title: The K2 and TESS Synergy. I. Updated Ephemerides and Parameters for K2-114, K2-167, K2-237, and K2-261

abstract: Although the Transiting Exoplanet Survey Satellite (TESS) primary mission observed the northern and southern ecliptic hemispheres, generally avoiding the ecliptic, and the Kepler space telescope during the K2 mission could only observe near the ecliptic, many of the K2 fields extend far enough from the ecliptic plane that sections overlap with TESS fields. Using photometric observations from both K2 and TESS, combined with archival spectroscopic observations, we globally modeled four known planetary systems discovered by K2 that were observed in the first year of the primary TESS mission. Specifically, we provide updated ephemerides and system parameters for K2-114 b, K2-167 b, K2-237 b, and K2-261 b. These were some of the first K2 planets to be observed by TESS in the first year and include three Jovian sized planets and a sub-Neptune with orbital periods less than 12 days. In each case, the updated ephemeris significantly reduces the uncertainty in prediction of future times of transit, which is valuable for planning observations with the James Webb Space Telescope and other future facilities. The TESS extended mission is expected to observe about half of the K2 fields, providing the opportunity to perform this type of analysis on a larger number of systems.


title: The Velocity Dispersion Function for Quiescent Galaxies in Nine Strong-Lensing Clusters

abstract: We measure the central stellar velocity dispersion function for quiescent galaxies in a set of nine northern clusters in the redshift range $0.18 < z < 0.29$ and with strong lensing arcs in Hubble Space Telescope images. The velocity dispersion function links galaxies directly to their dark matter halos. From dense SDSS and MMT/Hectospec spectroscopy we identify $231 - 479$ spectroscopic members in each cluster. We derive physical properties of cluster members including redshift, $D_{n}4000$, and central stellar velocity dispersion and we include a table of these measurements for 3419 cluster members. We construct the velocity dispersion functions for quiescent galaxies with $D_{n}4000 > 1.5$ and within $R_{200}$. The cluster velocity dispersion functions all show excesses at $\sigma \gtrsim 250 km s^{-1}$ compared to the field velocity dispersion function. The velocity dispersion function slope at large velocity dispersion ($\sigma > 160 km s^{-1}$) is steeper for more massive clusters, consistent with the trend observed for cluster luminosity functions. The spatial distribution of galaxies with large velocity dispersion at radii larger than $R_{200}$ further underscores the probable major role of dry mergers in the growth of massive cluster galaxies during cluster assembly.


title: Products of positive operators

abstract: On finite dimensional spaces, it is apparent that an operator is the product of two positive operators if and only if it is similar to a positive operator. Here, the class ${\mathcal L}^{+2}$ of bounded operators on separable infinite dimensional Hilbert spaces which can be written as the product of two bounded positive operators is studied. The structure is much richer, and connects (but is not equivalent to) quasi-similarity and quasi-affinity to a positive operator. The spectral properties of operators in ${\mathcal L}^{+2}$ are developed, and membership in ${\mathcal L}^{+2}$ among special classes, including algebraic and compact operators, is examined.


title: Distributed Safe Learning using an Invariance-based Safety Framework

abstract: In large-scale networks of uncertain dynamical systems, where communication is limited and there is a strong interaction among subsystems, learning local models and control policies offers great potential for designing high-performance controllers. At the same time, the lack of safety guarantees, here considered in the form of constraint satisfaction, prevents the use of data-driven techniques to safety-critical distributed systems. This paper presents a safety framework that guarantees constraint satisfaction for uncertain distributed systems while learning. The framework considers linear systems with coupling in the dynamics and subject to bounded parametric uncertainty, and makes use of robust invariance to guarantee safety. In particular, a robust non-convex invariant set, given by the union of multiple ellipsoidal invariant sets, and a nonlinear backup control law, given by the combination of multiple stabilizing linear feedbacks, are computed offline. In presence of unsafe inputs, the safety framework applies the backup control law, preventing the system to violate the constraints. As the robust invariant set and the backup stabilizing controller are computed offline, the online operations reduce to simple function evaluations, which enables the use of the proposed framework on systems with limited computational resources. The capabilities of the safety framework are illustrated by three numerical examples.


title: Parkinson's Disease Detection Using Ensemble Architecture from MR Images

abstract: Parkinson's Disease(PD) is one of the major nervous system disorders that affect people over 60. PD can cause cognitive impairments. In this work, we explore various approaches to identify Parkinson's using Magnetic Resonance (MR) T1 images of the brain. We experiment with ensemble architectures combining some winning Convolutional Neural Network models of ImageNet Large Scale Visual Recognition Challenge (ILSVRC) and propose two architectures. We find that detection accuracy increases drastically when we focus on the Gray Matter (GM) and White Matter (WM) regions from the MR images instead of using whole MR images. We achieved an average accuracy of 94.7\% using smoothed GM and WM extracts and one of our proposed architectures. We also perform occlusion analysis and determine which brain areas are relevant in the architecture decision making process.


title: X-ray Writing of Fiber Bragg Gratings (FBGs)

abstract: Fabrication of nanoscale optical Fiber Bragg Gratings (FBGs) is one of the key manufacturing processes for fiber optics, which has many important applications in data communication and distributed remote sensing at a distance up to hundreds of kms. However, the fabrication of the FBGs is challenging. To meet such need, we propose a novel method to write the FBGs using high-flux synchrotron x-ray with a nanoscale gold/Si gratings to modulate the x-ray flux and thus change the refractive index contrast of the optical fiber to form the FBGs. The nanoscale gold/Si with high aspect ratio has been fabricated at the Center for Nanoscale Materials (CNM) at Argonne National Laboratory (ANL) and the preliminary experiment of x-ray writing of the optical fiber FBGs has been carried out at the Advanced Photon Source (APS) of ANL. The preliminary experiment shows promising result of the x-ray writing method of the optical fiber FBGs.


title: Three-dimensional quantum gravity according to ST modular bootstrap

abstract: We combine the large-$c$ ST modular bootstrap equations with the Cardy formula for the asymptotic growth of the density of states to prove that any $2d$ unitary, compact, conformal field theory (CFT) with no higher spin conserved currents leads to conflicting inequalities whenever the entire spectrum of non-trivial primaries lies above the BTZ threshold. As a consequence, the holographic dual of $3d$ pure gravity, if it exists, cannot be a $2d$ CFT. Consistent solutions of ST bootstrap equations require additional primaries lying below the BTZ threshold. The lowest non-trivial primary necessarily has odd spin


title: Coloring linear hypergraphs: the Erd\H{o}s-Faber-Lov\'asz conjecture and the Combinatorial Nullstellensatz

abstract: The long-standing Erd\H{o}s-Faber-Lov\'asz conjecture states that every $n$-uniform linear hypergaph with $n$ edges has a proper vertex-coloring using $n$ colors. In this paper we propose an algebraic framework to the problem and formulate a corresponding stronger conjecture. Using the Combinatorial Nullstellensatz, we reduce the Erd\H{o}s-Faber-Lov\'asz conjecture to the existence of non-zero coefficients in certain polynomials. These coefficients are in turn related to the number of orientations with prescribed in-degree sequences of some auxiliary graphs. We prove the existence of certain orientations, which verifies a necessary condition for our algebraic approach to work.


title: Typical structure of hereditary graph families. I. Apex-free families

abstract: A family of graphs $\mathcal{F}$ is hereditary if $\mathcal{F}$ is closed under isomorphism and taking induced subgraphs. The speed of $\mathcal{F}$ is the sequence $\{|\mathcal{F}^n|\}_{n \in \mathbb{N}}$, where $\mathcal{F}^n$ denotes the set of graphs in $\mathcal{F}$ with the vertex set $[n]$. Alon, Balogh, Bollob\'{a}s and Morris [The structure of almost all graphs in a hereditary property, JCTB 2011] gave a rough description of typical graphs in a hereditary family and used it to show for every proper hereditary family $\mathcal{F}$ there exist $\varepsilon>0$ and an integer $l \geq 1$ such that $$|\mathcal{F}^n| = 2^{(1-1/l)n^2/2+o(n^{2-\varepsilon})}.$$ The main result of this paper gives a more precise description of typical structure for a restricted class of hereditary families. As a consequence we characterize hereditary families with the speed just above the threshold $2^{(1-1/l)n^2/2}$, generalizing a result of Balogh and Butterfield [Excluding induced subgraphs: Critical graphs, RSA 2011].


title: On Beltrami states near black hole event horizon

abstract: In this article, we study Beltrami equilibria for plasmas in near the horizon of a spinning black hole, and develop a framework for constructing the magnetic field profile in the near horizon limit for Clebsch flows in the single-fluid approximation. We find that the horizon profile for the magnetic field is shown to satisfy a system of first-order coupled ODEs dependent on a boundary condition for the magnetic field. For states in which the generalized vorticity vanishes (the generalized `superconducting' plasma state), the horizon profile becomes independent of the boundary condition, and depend only on the thermal properties of the plasma. Our analysis makes use of the full form for the time-independent Ampere's law in the 3+1 formalism, generalizing earlier conclusions for the case of vanishing vorticity, namely the complete magnetic field expulsion near the equator of an axisymmetric black horizon assuming that the thermal properties of the plasma are symmetric about the equatorial plane. For the general case, we find and discuss additional conditions required for the expulsion of magnetic fields at given points on the black hole horizon. We perform a length scale analysis which indicates the emergence of two distinct length scales characterizing the magnetic field variation and strength of the Beltrami term, respectively.


title: Typical structure of hereditary graph families. II. Exotic examples

abstract: A graph $G$ is $H$-free if it does not contain an induced subgraph isomorphic to $H$. The study of the typical structure of $H$-free graphs was initiated by Erd\H{o}s, Kleitman and Rothschild, who have shown that almost all $C_3$-free graphs are bipartite. Since then the typical structure of $H$-free graphs has been determined for several families of graphs $H$, including complete graphs, trees and cycles. Recently, Reed and Scott proposed a conjectural description of the typical structure of $H$-free graphs for all graphs $H$, which extends all previously known results in the area. We construct an infinite family of graphs for which the Reed-Scott conjecture fails, and use the methods we developed in the prequel paper to describe the typical structure of $H$-free graphs for graphs $H$ in this family. Using similar techniques, we construct an infinite family of graphs $H$ for which the maximum size of a homogenous set in a typical $H$-free graph is sublinear in the number of vertices, answering a question of Loebl et al. and Kang et al.


title: Rethink Maximum Mean Discrepancy for Domain Adaptation

abstract: Existing domain adaptation methods aim to reduce the distributional difference between the source and target domains and respect their specific discriminative information, by establishing the Maximum Mean Discrepancy (MMD) and the discriminative distances. However, they usually accumulate to consider those statistics and deal with their relationships by estimating parameters blindly. This paper theoretically proves two essential facts: 1) minimizing the MMD equals to maximize the source and target intra-class distances respectively but jointly minimize their variance with some implicit weights, so that the feature discriminability degrades; 2) the relationship between the intra-class and inter-class distances is as one falls, another rises. Based on this, we propose a novel discriminative MMD. On one hand, we consider the intra-class and inter-class distances alone to remove a redundant parameter, and the revealed weights provide their approximate optimal ranges. On the other hand, we design two different strategies to boost the feature discriminability: 1) we directly impose a trade-off parameter on the implicit intra-class distance in MMD to regulate its change; 2) we impose the similar weights revealed in MMD on inter-class distance and maximize it, then a balanced factor could be introduced to quantitatively leverage the relative importance between the feature transferability and its discriminability. The experiments on several benchmark datasets not only prove the validity of theoretical results but also demonstrate that our approach could perform better than the comparative state-of-art methods substantially.


title: Nonequilibrium many-body quantum engine driven by time-translation symmetry breaking

abstract: Quantum many-body systems out of equilibrium can host intriguing phenomena such as transitions to exotic dynamical states. Although this emergent behaviour can be observed in experiments, its potential for technological applications is largely unexplored. Here, we investigate the impact of collective effects on quantum engines that extract mechanical work from a many-body system. Using an opto-mechanical cavity setup with an interacting atomic gas as a working fluid, we demonstrate theoretically that such engines produce work under periodic driving. The stationary cycle of the working fluid features nonequilibrium phase transitions, resulting in abrupt changes of the work output. Remarkably, we find that our many-body quantum engine operates even without periodic driving. This phenomenon occurs when its working fluid enters a phase that breaks continuous time-translation symmetry: the emergent time-crystalline phase can sustain the motion of a load generating mechanical work. Our findings pave the way for designing novel nonequilibrium quantum machines.


title: Falsification-Based Robust Adversarial Reinforcement Learning

abstract: Reinforcement learning (RL) has achieved enormous progress in solving various sequential decision-making problems, such as control tasks in robotics. Since policies are overfitted to training environments, RL methods have often failed to be generalized to safety-critical test scenarios. Robust adversarial RL (RARL) was previously proposed to train an adversarial network that applies disturbances to a system, which improves the robustness in test scenarios. However, an issue of neural network-based adversaries is that integrating system requirements without handcrafting sophisticated reward signals are difficult. Safety falsification methods allow one to find a set of initial conditions and an input sequence, such that the system violates a given property formulated in temporal logic. In this paper, we propose falsification-based RARL (FRARL): this is the first generic framework for integrating temporal logic falsification in adversarial learning to improve policy robustness. By applying our falsification method, we do not need to construct an extra reward function for the adversary. Moreover, we evaluate our approach on a braking assistance system and an adaptive cruise control system of autonomous vehicles. Our experimental results demonstrate that policies trained with a falsification-based adversary generalize better and show less violation of the safety specification in test scenarios than those trained without an adversary or with an adversarial network.


title: Scattering of rod-like swimmers in low Reynolds number environments

abstract: In their search for metabolic resources microbes swim through viscous environments that present physical anisotropies, including steric obstacles across a wide range of sizes. Hydrodynamic forces are known to significantly alter swimmer trajectories near flat and low-curvature surfaces. In this work, we imaged hundreds-of-thousands of high-curvature scattering interactions between swimming bacteria and micro-fabricated pillars with radii from ~1 to ~10 cell lengths. As a function of impact parameter, cell-pillar interactions produced distinct chiral distributions for scattering angle -- including unexpected 'counter-rotator' trajectories -- well-described by a sterics-only model. Our data and model suggest that alteration of swimmer trajectories is subject to distinct mechanisms when interacting with objects of different size; primarily steric for objects below ~10 cell lengths and requiring incorporation of hydrodynamics at larger scales. These alterations in trajectory impact swim dynamics and may affect microbial populations in ways that depend on the shape and placement of obstacles within an environment.


title: On a spectral pattern of the Von-Neumann probes

abstract: In this paper we have extended our previous work \citep{probe1} and considered spectral characteristics of interstellar non-relativistic Type-II and Type-III Von-Neumann probes. It has been shown that by means of the proton capture, required for replication, strong bremstrahlung emission will be generated. We have found that for both types of civilizations the probes might be visible mainly in the infrared spectral band, but as it has been found the probes might be visible in the ultraviolet as well. For both cases it was shown that differential power has dips for equally spaced frequencies, which might be a significant fingerprint to identify such exotic interstellar objects.


title: Coherent Unipolar Cherenkov and Diffraction Radiation Generated by Relativistic Electrons

abstract: For usual radiation the integral of the electric field strength over the time is zero. First time the possibility of unipolar radiation generation has been considered theoretically in the paper [EG.Bessonov. Preprint No 76 of the Lebedev Physical Institute of the Academy of Science of the USSR, Moscow, 1990]. According to this work the unipolar radiation is radiation for which the integral of the electric field strength over the time differs significantly from zero. Futher a number of theoretical articles have been devoted to this problem mainly in respect to a synchrotron radiation. However up to now there are no experimental investigations of this phenomenon. In this paper we present the results of experimental observation of the unipolar Cherenkov and diffraction radiation generated by relativistic electrons in millimeter wavelength region. For this purpose the detector sensitive to the selected direction of the electric field strength has been elaborated. We observed the coherent Cherenkov radiation and backward diffraction radiation of bunched electron beam when the electrons moving near the targets. The partial unipolarity effect has been registered for the Cherenkov radiation and the almost full unipolarity has been observed for the backward diffraction radiation.


title: Experience Report: Smuggling a Little Bit of Coq Inside a CAD Development Context (Extended Abstract)

abstract: While the use of formal verification techniques is well established in the development of mission-critical software, it is still rare in the production of most other kinds of software. We share our experience that a formal verification tool such as Coq can be very useful and practical in the context of off-the-shelf software development -- CAD in particular -- at least in some occasions. The emphasis is on 3 main areas: factors that can enable the use of Coq in an industrial context; some typical examples of tasks, where Coq can offer an advantage; examples of issues to overcome - and some non-issues - when integrating Coq in a standard development process.


title: Role of cosmic rays in the early stages of galactic outflows

abstract: Using an idealized set-up, we investigate the dynamical role of cosmic rays (CRs) in the early stages of galactic outflows for galaxies of halo masses $10^{8}$, $10^{11}$ and $10^{12}$ $M_\odot$. The outflow is launched from a central region in the galactic disk where we consider three different constant star formation rates ($0.1$, $1$, and $10$ $M_{\odot}$ yr$^{-1}$) over a dynamical timescale of $50$ Myr. We determine the temperature distribution of the gas and find that CRs can reduce the temperature of the shocked gas, which is consistent with previous results. However, we show that CRs do not have any noticeable effect on the mass loading by the outflow. We find that CRs can reduce the size of the outflow, which contradicts previous claims of efficient dynamical impact of CRs; however, it is consistent with earlier theoretical models of cosmic ray driven blastwave as well as stellar wind. We discuss the dependence of our results on CR injection prescriptions and compare them with earlier studies. We conclude that in the early stages of galactic outflows the dynamical role of CRs is not important.


title: Canonical forms of two-qubit states under local operations

abstract: Canonical forms of two-qubits under the action of stochastic local operations and classical communications (SLOCC) offer great insight for understanding non-locality and entanglement shared by them. They also enable geometric picture of two-qubit states within the Bloch ball. It has been shown (Verstraete et.al. {Phys. Rev. A, 64, 010101(R) (2001)) that an arbitrary two-qubit state gets transformed under SLOCC into one of the {\em two} different canonical forms. One of these happens to be the Bell diagonal form of two-qubit states and the other non-diagonal canonical form is obtained for a family of rank deficient two-qubit states. The method employed by Verstraete et.al. required highly non-trivial results on matrix decompositions in $n$ dimensional spaces with indefinite metric. Here we employ an entirely different approach -- inspired by the methods developed by Rao et. al., (J. Mod. Opt. 45, 955 (1998)) in classical polarization optics -- which leads naturally towards the identification of two inequivalent SLOCC invariant canonical forms for two-qubit states. In addition, our approach results in a simple geometric visualization of two-qubit states in terms of their SLOCC canonical forms.


title: Comment on "A tight distance-dependent estimator for screening three-center Coulomb integrals over Gaussian basis functions" [J. Chem. Phys. 142, 154106 (2015)]

abstract: We extend the tight distance-dependent estimator proposed by Hollman et al. [J. Chem. Phys. 142, 154106 (2015)] for the three-center Coulomb integrals over Gaussian atomic orbitals to handle the two-center case. We also propose minor modifications of the original three-center estimator for the case of contracted ket Gaussians and concentric bra Gaussians.


title: Accelerated Message Passing for Entropy-Regularized MAP Inference

abstract: Maximum a posteriori (MAP) inference in discrete-valued Markov random fields is a fundamental problem in machine learning that involves identifying the most likely configuration of random variables given a distribution. Due to the difficulty of this combinatorial problem, linear programming (LP) relaxations are commonly used to derive specialized message passing algorithms that are often interpreted as coordinate descent on the dual LP. To achieve more desirable computational properties, a number of methods regularize the LP with an entropy term, leading to a class of smooth message passing algorithms with convergence guarantees. In this paper, we present randomized methods for accelerating these algorithms by leveraging techniques that underlie classical accelerated gradient methods. The proposed algorithms incorporate the familiar steps of standard smooth message passing algorithms, which can be viewed as coordinate minimization steps. We show that these accelerated variants achieve faster rates for finding $\epsilon$-optimal points of the unregularized problem, and, when the LP is tight, we prove that the proposed algorithms recover the true MAP solution in fewer iterations than standard message passing algorithms.


title: Response by the Montreal AI Ethics Institute to the Santa Clara Principles on Transparency and Accountability in Online Content Moderation

abstract: In April 2020, the Electronic Frontier Foundation (EFF) publicly called for comments on expanding and improving the Santa Clara Principles on Transparency and Accountability (SCP), originally published in May 2018. The Montreal AI Ethics Institute (MAIEI) responded to this call by drafting a set of recommendations based on insights and analysis by the MAIEI staff and supplemented by workshop contributions from the AI Ethics community convened during two online public consultation workshops. In its submission, MAIEI provides 12 overarching recommendations for the SCP, these include: 1) ensure there is more diversity in the content moderation process; 2) increase transparency into how platforms guide content-ranking; 3) disclose anonymized data on the training and/or cultural background of the content moderators for a platform; 4) tailor content moderation tools for specific issues; 5) draft specific guidelines for messaging applications with regards to data protection in content moderation; 6) take into account cultural differences relevant to what constitutes acceptable behavior online; 7) ensure platforms are transparent in regards to political advertising; 8) ensure greater transparency into the user-generated flagging/reporting systems deployed by a platform; 9) clarify if user content is flagged or reported through an automated system; 10) provide more data on the types of content removed from platforms; 11) provide clear guidelines on the appeal process, as well as data on prior appeals; 12) create a system for periodically revisiting the SCP so it reflects various technological advancements, modifications in law and policy, as well as changing trends or movements in content moderation.


title: Nonlocal Effective Electromagnetic Wave Characteristics of Composite Media: Beyond the Quasistatic Regime

abstract: We derive exact nonlocal expressions for the effective dielectric constant tensor ${\boldsymbol \varepsilon}_e({\bf k}_I, \omega)$ of disordered two-phase composites and metamaterials from first principles. This formalism extends the long-wavelength limitations of conventional homogenization estimates of ${\boldsymbol \varepsilon}_e({\bf k}_I, \omega)$ for arbitrary microstructures so that it can capture spatial dispersion well beyond the quasistatic regime (where $\omega$ and ${\bf k}_I$ are frequency and wavevector of the incident radiation). This is done by deriving nonlocal strong-contrast expansions that exactly account for multiple scattering for the range of wavenumbers for which our extended homogenization theory applies, i.e., $0 \le |{\bf k}_I| \ell \lesssim 1$ (where $\ell$ is a characteristic heterogeneity length scale). Due to the fast-convergence properties of such expansions, their lower-order truncations yield accurate closed-form approximate formulas for ${\varepsilon}_e({\bf k}_I,\omega)$ that incorporate microstructural information via the spectral density, which is easy to compute for any composite. The accuracy of these microstructure-dependent approximations is validated by comparison to full-waveform simulation methods for both 2D and 3D ordered and disordered models of composite media. Thus, our closed-form formulas enable one to predict accurately and efficiently the effective wave characteristics well beyond the quasistatic regime without having to perform full-blown simulations. Among other results, we show that certain disordered hyperuniform particulate composites exhibit novel wave characteristics. Our results demonstrate that one can design the effective wave characteristics of a disordered composite by engineering the microstructure to possess tailored spatial correlations at prescribed length scales.


title: The magnetized disk-halo transition region of M51

abstract: The grand-design face-on spiral galaxy M51 is an excellent laboratory for studying magnetic fields in galaxies. We present new observations of M51 using the VLA at the frequency range of S-band (2-4GHz), to shed new light on the transition region between the disk and halo. We present images of the distributions of the total intensity, polarized intensity, degree of polarization, and rotation measure (RM). The RM distribution in S-band shows a fluctuating pattern without any apparent large-scale structure. We discuss a model of the depolarization of synchrotron radiation in a multi-layer magneto-ionic medium and compare the model predictions to the polarization data of M51 between 1-8GHz. Since the model predictions strongly differ within the wavelength range of the S-band, the new data are essential. The parameters of the model are adjusted to fit to the data of polarization fractions in a few selected regions. In three spiral arm regions, the turbulent field in the disk dominates with strengths between 18muG and 24muG, while the regular field strengths are 8-16muG. In one inter-arm region, the regular field strength of 18muG exceeds that of the turbulent field of 11muG. The regular field strengths in the halo are 3-5muG. The observed RMs in the disk-halo transition region are probably dominated by tangled regular fields, as predicted from models of evolving dynamos, and/or vertical fields, as predicted from numerical simulations of Parker instabilities or galactic winds. Both types of magnetic fields have frequent reversals on scales similar to or larger than the beam size (550pc) that contribute to an increase of the RM dispersion and to distortions of any large-scale pattern of the regular field. Our study devises new ways of analyzing and interpreting broadband multi-frequency polarization data that will be applicable to future data from, for example, the Square Kilometre Array.


title: Implications of a possible TeV break in the cosmic-ray electron and positron flux

abstract: A TeV spectral break in the total flux of cosmic-ray electrons and positrons (CREs) at which the spectral power index softens from $\sim 3$ to $\sim 4$ has been observed by H.E.S.S. and recently confirmed by DAMPE with a high significance of $6.6~\sigma$. Such an observation is apparently inconsistent with the data from other experiments such as Fermi-LAT, AMS-02 and CALET. We perform a global analysis to the latest CRE data including Fermi-LAT, AMS-02, CALET, DAMPE and H.E.S.S. with energy scale uncertainties taken into account to improve the consistency between the data sets. The fit result strongly favors the existence of the break at $\sim 1$ TeV with an even higher statistical significance of $13.3~\sigma$. In view of the tentative CRE break, we revisit a number of models of nearby sources, such as a single generic Pulsar Wind Nebula (PWN), known multiple PWNe from the ATNF catalog, and their combinations with either an additional Dark Matter (DM) component or a Supernova Remnant (SNR). We show that the CRE break at $\sim 1$ TeV, together with the known CR positron excess points towards the possibility that the nearby sources should be highly charge asymmetric. Among the models under consideration, the one with a PWN plus SNR is most favored by the current data. The favoured distance and age of the PWN and SNR sources are both within $0.6$ kpc and around $10^{5}$ yr respectively. Possible candidate sources include PSR J0954-5430, Vela and Monogem ring, etc. We find that for the models under consideration, the additional DM component is either unnecessary, or predicts too much photons in tension with the H.E.S.S. data of $\gamma$-rays from the direction of the Galactic Center. We also show that the current measurement of the anisotropies in the arrival direction of the CRE can be useful in determining the property of the sources.


title: Modeling features in the redshift-space halo power spectrum with perturbation theory

abstract: We study the ability of perturbative models with effective field theory contributions and infra-red resummation to model the redshift space clustering of biased tracers in models where the linear power spectrum has "features" -- either imprinted during inflation or induced by non-standard expansion histories. We show that both Eulerian and Lagrangian perturbation theory are capable of reproducing the Fourier space two-point functions of halos up to the non-linear scale from a suite of $4096^3$ particle N-body simulations. This is the first demonstration that perturbative models can accurately fit the redshift-space clustering of biased tracers in N-body simulations of such theories. By comparing different theoretical models and IR resummation schemes we assess the current theoretical uncertainty in predicting power spectra for models with features. Our results suggest that future surveys will be able to detect or tightly constrain features in the primordial spectrum below the one percent level across a wide range of scales.


title: Performance analysis of Zero Black-Derman-Toy interest rate model in catastrophic events: COVID-19 case study

abstract: In this paper we continue the research of our recent interest rate tree model called Zero Black-Derman-Toy (ZBDT) model, which includes the possibility of a jump at each step to a practically zero interest rate. This approach allows to better match to risk of financial slowdown caused by catastrophic events. We present how to valuate a wide range of financial derivatives for such a model. The classical Black-Derman-Toy (BDT) model and novel ZBDT model are described and analogies in their calibration methodology are established. Finally two cases of applications of the novel ZBDT model are introduced. The first of them is the hypothetical case of an S-shape term structure and decreasing volatility of yields. The second case is an application of the ZBDT model in the structure of United States sovereign bonds in the current $2020$ economic slowdown caused by the Coronavirus pandemic. The objective of this study is to understand the differences presented by the valuation in both models for different derivatives.


title: DPCP-p: A Distributed Locking Protocol for Parallel Real-Time Tasks

abstract: Real-time scheduling and locking protocols are fundamental facilities to construct time-critical systems. For parallel real-time tasks, predictable locking protocols are required when concurrent sub-jobs mutually exclusive access to shared resources. This paper for the first time studies the distributed synchronization framework of parallel real-time tasks, where both tasks and global resources are partitioned to designated processors, and requests to each global resource are conducted on the processor on which the resource is partitioned. We extend the Distributed Priority Ceiling Protocol (DPCP) for parallel tasks under federated scheduling, with which we proved that a request can be blocked by at most one lower-priority request. We develop task and resource partitioning heuristics and propose analysis techniques to safely bound the task response times. Numerical evaluation (with heavy tasks on 8-, 16-, and 32-core processors) indicates that the proposed methods improve the schedulability significantly compared to the state-of-the-art locking protocols under federated scheduling.


title: An integrable family of torqued, damped, rigid rotors

abstract: Expositions of the Euler equations for the rotation of a rigid body often invoke the idea of a specially damped system whose energy dissipates while its angular momentum magnitude is conserved in the body frame. An attempt to explicitly construct such a damping function leads to a more general, but still integrable, system of cubic equations whose trajectories are confined to nested sets of quadric surfaces in angular momentum space. For some choices of parameters, the lines of fixed points along both the largest and smallest moment of inertia axes can be simultaneously attracting. Limiting cases are those that conserve either the energy or the magnitude of the angular momentum. Parallels with rod mechanics, micromagnetics, and particles with effective mass are briefly discussed.


title: Learning Search Space Partition for Black-box Optimization using Monte Carlo Tree Search

abstract: High dimensional black-box optimization has broad applications but remains a challenging problem to solve. Given a set of samples $\{\vx_i, y_i\}$, building a global model (like Bayesian Optimization (BO)) suffers from the curse of dimensionality in the high-dimensional search space, while a greedy search may lead to sub-optimality. By recursively splitting the search space into regions with high/low function values, recent works like LaNAS shows good performance in Neural Architecture Search (NAS), reducing the sample complexity empirically. In this paper, we coin LA-MCTS that extends LaNAS to other domains. Unlike previous approaches, LA-MCTS learns the partition of the search space using a few samples and their function values in an online fashion. While LaNAS uses linear partition and performs uniform sampling in each region, our LA-MCTS adopts a nonlinear decision boundary and learns a local model to pick good candidates. If the nonlinear partition function and the local model fits well with ground-truth black-box function, then good partitions and candidates can be reached with much fewer samples. LA-MCTS serves as a \emph{meta-algorithm} by using existing black-box optimizers (e.g., BO, TuRBO) as its local models, achieving strong performance in general black-box optimization and reinforcement learning benchmarks, in particular for high-dimensional problems.


title: Computing Conceptual Distances between Breast Cancer Screening Guidelines: An Implementation of a Near-Peer Epistemic Model of Medical Disagreement

abstract: Using natural language processing tools, we investigate the differences of recommendations in medical guidelines for the same decision problem -- breast cancer screening. We show that these differences arise from knowledge brought to the problem by different medical societies, as reflected in the conceptual vocabularies used by the different groups of authors.The computational models we build and analyze agree with the near-peer epistemic model of expert disagreement proposed by Garbayo. Even though the article is a case study focused on one set of guidelines, the proposed methodology is broadly applicable. In addition to proposing a novel graph-based similarity model for comparing collections of documents, we perform an extensive analysis of the model performance. In a series of a few dozen experiments, in three broad categories, we show, at a very high statistical significance level of 3-4 standard deviations for our best models, that the high similarity between expert annotated model and our concept based, automatically created, computational models is not accidental. Our best model achieves roughly 70% similarity. We also describe possible extensions of this work.


title: "QM19 summary talk": Outlook and future of heavy-ion collisions

abstract: A summary of the QM19 conference is given by highlighting a few selected results. These are discussed as examples to illustrate the exciting future of heavy-ion collisions and the need for further instrumentation. (The arXiv version is significantly longer than the printed proceedings, with more figures.)


title: ConFoc: Content-Focus Protection Against Trojan Attacks on Neural Networks

abstract: Deep Neural Networks (DNNs) have been applied successfully in computer vision. However, their wide adoption in image-related applications is threatened by their vulnerability to trojan attacks. These attacks insert some misbehavior at training using samples with a mark or trigger, which is exploited at inference or testing time. In this work, we analyze the composition of the features learned by DNNs at training. We identify that they, including those related to the inserted triggers, contain both content (semantic information) and style (texture information), which are recognized as a whole by DNNs at testing time. We then propose a novel defensive technique against trojan attacks, in which DNNs are taught to disregard the styles of inputs and focus on their content only to mitigate the effect of triggers during the classification. The generic applicability of the approach is demonstrated in the context of a traffic sign and a face recognition application. Each of them is exposed to a different attack with a variety of triggers. Results show that the method reduces the attack success rate significantly to values < 1% in all the tested attacks while keeping as well as improving the initial accuracy of the models when processing both benign and adversarial data.


title: High Resolution, Wide Field, Narrow Band, Snapshot Imaging

abstract: We investigate the imaging performance of an interferometric array in the case of wide field, high resolution, narrow band, snapshot imaging. We find that, when uv-cell sizes are sufficiently small (ie. image sizes are sufficiently large), each instantaneous visibility record is gridded into its own uv-cell. This holds even for dense arrays, like the core of the next generation VLA. In this particular, application, Uniform weighting of the gridded visibilities approaches Natural weighting, with its often deleterious consequences on the resulting synthesized beam. For a core-dominated array, we show that the resulting image noise is highly correlated on scales comparable to the spatial frequencies of the core baselines. In general, this study accentuates the fact that, for imaging applications that require high resolution (Plains array and greater), many of the core antennas can be employed as a separate subarray for low resolution science, without sacrificing the quality of the high resolution science.


title: Embeddings of Function Spaces via the Caffarelli-Silvestre Extension, Capacities and Wolff potentials

abstract: Let $P_{\alpha} f(x,t)$ be the Caffarelli-Silvestre extension of a smooth function $f(x): \mathbb{R}^n \rightarrow \mathbb{R}^{n+1}_+:=\mathbb{R}^n\times (0,\infty).$ The purpose of this article is twofold. Firstly, we want to characterize a nonnegative measure $\mu$ on $\mathbb{R}^{n+1}_+$ such that $f(x)\rightarrow P_{\alpha} f(x,t)$ induces bounded embeddings from the Lebesgue spaces $L^p(\mathbb{R}^n)$ to the $L^q(\mathbb{R}^{n+1}_+,\mu).$ On one hand, these embeddings will be characterized by using a newly introduced $L^p-$capacity associated with the Caffarelli-Silvestre extension. In doing so, the mixed norm estimates of $P_{\alpha} f(x,t),$ the dual form of the $L^p-$capacity, the $L^p-$capacity of general balls, and a capacitary strong type inequality will be established, respectively. On the other hand, when $p>q>1,$ these embeddings will also be characterized in terms of the Hedberg-Wolff potential of $\mu.$ Secondly, we characterize a nonnegative measure $\mu$ on $\mathbb{R}^{n+1}_+$ such that $f(x)\rightarrow P_{\alpha} f(x,t)$ induces bounded embeddings from the homogeneous Sobolev spaces $\dot{W}^{\beta,p}(\mathbb{R}^n)$ to the $L^q(\mathbb{R}^{n+1}_+,\mu)$ in terms of the fractional perimeter of open sets for endpoint cases and the fractional capacity for general cases.


title: Quantifying intrinsic causal contributions via structure preserving interventions

abstract: We propose a notion of causal influence that describes the `intrinsic' part of the contribution of a node on a target node in a DAG. By recursively writing each node as a function of the upstream noise terms, we separate the intrinsic information added by each node from the one obtained from its ancestors. To interpret the intrinsic information as a {\it causal} contribution, we consider `structure-preserving interventions' that randomize each node in a way that mimics the usual dependence on the parents and does not perturb the observed joint distribution. To get a measure that is invariant with respect to relabelling nodes we use Shapley based symmetrization and show that it reduces in the linear case to simple ANOVA after resolving the target node into noise variables. We describe our contribution analysis for variance and entropy, but contributions for other target metrics can be defined analogously. The code is available in the package gcm of the open source library DoWhy.


title: Bayesian Coresets: Revisiting the Nonconvex Optimization Perspective

abstract: Bayesian coresets have emerged as a promising approach for implementing scalable Bayesian inference. The Bayesian coreset problem involves selecting a (weighted) subset of the data samples, such that the posterior inference using the selected subset closely approximates the posterior inference using the full dataset. This manuscript revisits Bayesian coresets through the lens of sparsity constrained optimization. Leveraging recent advances in accelerated optimization methods, we propose and analyze a novel algorithm for coreset selection. We provide explicit convergence rate guarantees and present an empirical evaluation on a variety of benchmark datasets to highlight our proposed algorithm's superior performance compared to state-of-the-art on speed and accuracy.


title: Monte-Carlo Imaging for Optical Interferometry

abstract: We present a flexible code created for imaging from the bispectrum and visibility-squared. By using a simulated annealing method, we limit the probability of converging to local chi-squared minima as can occur when traditional imaging methods are used on data sets with limited phase information. We present the results of our code used on a simulated data set utilizing a number of regularization schemes including maximum entropy. Using the statistical properties from Monte-Carlo Markov chains of images, we show how this code can place statistical limits on image features such as unseen binary companions.


title: Adaptive Discretization for Model-Based Reinforcement Learning

abstract: We introduce the technique of adaptive discretization to design an efficient model-based episodic reinforcement learning algorithm in large (potentially continuous) state-action spaces. Our algorithm is based on optimistic one-step value iteration extended to maintain an adaptive discretization of the space. From a theoretical perspective we provide worst-case regret bounds for our algorithm which are competitive compared to the state-of-the-art model-based algorithms. Moreover, our bounds are obtained via a modular proof technique which can potentially extend to incorporate additional structure on the problem. From an implementation standpoint, our algorithm has much lower storage and computational requirements due to maintaining a more efficient partition of the state and action spaces. We illustrate this via experiments on several canonical control problems, which shows that our algorithm empirically performs significantly better than fixed discretization in terms of both faster convergence and lower memory usage. Interestingly, we observe empirically that while fixed-discretization model-based algorithms vastly outperform their model-free counterparts, the two achieve comparable performance with adaptive discretization.


title: Modular Fuss-Catalan numbers

abstract: The modular Catalan numbers $C_{k,n}$, introduced by Hein and Huang in 2016 count equivalence classes of parenthesizations of $x_0 * x_1 * \dots *x_n$ where $*$ is a binary $k$-associative operation and $k$ is a positive integer. The classical notion of associativity is just 1-associativity, in which case $C_{1,n} = 1$ and the size of the unique class is given by the Catalan number $C_n$. In this paper we introduce modular Fuss-Catalan numbers $C_{k,n}^{m}$ which count equivalence classes of parenthesizations of $x_0 * x_1 * \dots *x_n$ where $*$ is an $m$-ary $k$-associative operation for $m \geq 2$. Our main results are a closed formula for $C_{k,n}^{m}$ and a characterisation of $k$-associativity.


title: Open Questions on Nonequilibrium Thermodynamics of Chemical Reaction Networks

abstract: Chemical reaction networks (CRNs) are prototypical complex systems because reactions are nonlinear and connected in intricate ways, and they are also essential to understand living systems. Here, I discuss how recent developments in nonequilibrium thermodynamics provide new insight on how CRNs process energy and perform sophisticated tasks, and describe open challenges in the field.


title: Adversarial Example Games

abstract: The existence of adversarial examples capable of fooling trained neural network classifiers calls for a much better understanding of possible attacks to guide the development of safeguards against them. This includes attack methods in the challenging non-interactive blackbox setting, where adversarial attacks are generated without any access, including queries, to the target model. Prior attacks in this setting have relied mainly on algorithmic innovations derived from empirical observations (e.g., that momentum helps), lacking principled transferability guarantees. In this work, we provide a theoretical foundation for crafting transferable adversarial examples to entire hypothesis classes. We introduce Adversarial Example Games (AEG), a framework that models the crafting of adversarial examples as a min-max game between a generator of attacks and a classifier. AEG provides a new way to design adversarial examples by adversarially training a generator and a classifier from a given hypothesis class (e.g., architecture). We prove that this game has an equilibrium, and that the optimal generator is able to craft adversarial examples that can attack any classifier from the corresponding hypothesis class. We demonstrate the efficacy of AEG on the MNIST and CIFAR-10 datasets, outperforming prior state-of-the-art approaches with an average relative improvement of $29.9\%$ and $47.2\%$ against undefended and robust models (Table 2 & 3) respectively.


title: Some i-Mark games

abstract: Let $S$ be a set of positive integers, and let $D$ be a set of integers larger than $1$. The game $i$-Mark$(S,D)$ is an impartial combinatorial game introduced by Sopena (2016), which is played with a single pile of tokens. In each turn, a player can subtract $s \in S$ from the pile, or divide the size of the pile by $d \in D$, if the pile size is divisible by $d$. Sopena partially analyzed the games with $S=[1, t-1]$ and $D=\{d\}$ for $d \not\equiv 1 \pmod t$, but left the case $d \equiv 1 \pmod t$ open. We solve this problem by calculating the Sprague-Grundy function of $i$-Mark$([1,t-1],\{d\})$ for $d \equiv 1 \pmod t$, for all $t,d \geq 2$. We also calculate the Sprague-Grundy function of $i$-Mark$(\{2\},\{2k + 1\})$ for all $k$, and show that it exhibits similar behavior. Finally, following Sopena's suggestion to look at games with $|D|>1$, we derive some partial results for the game $i$-Mark$(\{1\}, \{2, 3\})$, whose Sprague-Grundy function seems to behave erratically and does not show any clean pattern. We prove that each value $0,1,2$ occurs infinitely often in its SG sequence, with a maximum gap length between consecutive appearances.


title: Sequential Transfer in Reinforcement Learning with a Generative Model

abstract: We are interested in how to design reinforcement learning agents that provably reduce the sample complexity for learning new tasks by transferring knowledge from previously-solved ones. The availability of solutions to related problems poses a fundamental trade-off: whether to seek policies that are expected to achieve high (yet sub-optimal) performance in the new task immediately or whether to seek information to quickly identify an optimal solution, potentially at the cost of poor initial behavior. In this work, we focus on the second objective when the agent has access to a generative model of state-action pairs. First, given a set of solved tasks containing an approximation of the target one, we design an algorithm that quickly identifies an accurate solution by seeking the state-action pairs that are most informative for this purpose. We derive PAC bounds on its sample complexity which clearly demonstrate the benefits of using this kind of prior knowledge. Then, we show how to learn these approximate tasks sequentially by reducing our transfer setting to a hidden Markov model and employing spectral methods to recover its parameters. Finally, we empirically verify our theoretical findings in simple simulated domains.


title: Scalable Monte Carlo Inference and Rescaled Local Asymptotic Normality

abstract: In this paper, we generalize the property of local asymptotic normality (LAN) to an enlarged neighborhood, under the name of rescaled local asymptotic normality (RLAN). We obtain sufficient conditions for a regular parametric model to satisfy RLAN. We show that RLAN supports the construction of a statistically efficient estimator which maximizes a cubic approximation to the log-likelihood on this enlarged neighborhood. In the context of Monte Carlo inference, we find that this maximum cubic likelihood estimator can maintain its statistical efficiency in the presence of asymptotically increasing Monte Carlo error in likelihood evaluation.