title: Preludes to the Eilenberg-Moore and the Leray-Serre spectral sequences

abstract: The Leray-Serre and the Eilenberg-Moore spectral sequence are fundamental tools for computing the cohomology of a group or, more generally, of a space. We describe the relationship between these two spectral sequences when both of them share the same abutment. There exists a joint tri-graded refinement of the Leray-Serre and the Eilenberg-Moore spectral sequence. This refinement involves two more spectral sequences, the preludes from the title, which abut to the initial terms of the Leray-Serre and the Eilenberg-Moore spectral sequence, respectively. We show that one of these always degenerates from its second page on and that the other one satisfies a local-to-global property: It degenerates for all possible base spaces if and only if it does so when the base space is contractible. When the preludes degenerate early enough, they appear to echo Deligne's d\'ecalage, but in general, this is an illusion. We discuss several principal fibrations to illustrate the possible cases and give applications, in particular, to Lie groups, torus bundles, and generalizations.


title: Combining Pseudo-Point and State Space Approximations for Sum-Separable Gaussian Processes

abstract: Gaussian processes (GPs) are important probabilistic tools for inference and learning in spatio-temporal modelling problems such as those in climate science and epidemiology. However, existing GP approximations do not simultaneously support large numbers of off-the-grid spatial data-points and long time-series which is a hallmark of many applications. Pseudo-point approximations, one of the gold-standard methods for scaling GPs to large data sets, are well suited for handling off-the-grid spatial data. However, they cannot handle long temporal observation horizons effectively reverting to cubic computational scaling in the time dimension. State space GP approximations are well suited to handling temporal data, if the temporal GP prior admits a Markov form, leading to linear complexity in the number of temporal observations, but have a cubic spatial cost and cannot handle off-the-grid spatial data. In this work we show that there is a simple and elegant way to combine pseudo-point methods with the state space GP approximation framework to get the best of both worlds. The approach hinges on a surprising conditional independence property which applies to space--time separable GPs. We demonstrate empirically that the combined approach is more scalable and applicable to a greater range of spatio-temporal problems than either method on its own.


title: Integrated information as a common signature of dynamical and information-processing complexity

abstract: The apparent dichotomy between information-processing and dynamical approaches to complexity science forces researchers to choose between two diverging sets of tools and explanations, creating conflict and often hindering scientific progress. Nonetheless, given the shared theoretical goals between both approaches, it is reasonable to conjecture the existence of underlying common signatures that capture interesting behaviour in both dynamical and information-processing systems. Here we argue that a pragmatic use of Integrated Information Theory (IIT), originally conceived in theoretical neuroscience, can provide a potential unifying framework to study complexity in general multivariate systems. Furthermore, by leveraging metrics put forward by the integrated information decomposition ($\Phi$ID) framework, our results reveal that integrated information can effectively capture surprisingly heterogeneous signatures of complexity -- including metastability and criticality in networks of coupled oscillators as well as distributed computation and emergent stable particles in cellular automata -- without relying on idiosyncratic, ad-hoc criteria. These results show how an agnostic use of IIT can provide important steps towards bridging the gap between informational and dynamical approaches to complex systems.


title: Residual Error: a New Performance Measure for Adversarial Robustness

abstract: Despite the significant advances in deep learning over the past decade, a major challenge that limits the wide-spread adoption of deep learning has been their fragility to adversarial attacks. This sensitivity to making erroneous predictions in the presence of adversarially perturbed data makes deep neural networks difficult to adopt for certain real-world, mission-critical applications. While much of the research focus has revolved around adversarial example creation and adversarial hardening, the area of performance measures for assessing adversarial robustness is not well explored. Motivated by this, this study presents the concept of residual error, a new performance measure for not only assessing the adversarial robustness of a deep neural network at the individual sample level, but also can be used to differentiate between adversarial and non-adversarial examples to facilitate for adversarial example detection. Furthermore, we introduce a hybrid model for approximating the residual error in a tractable manner. Experimental results using the case of image classification demonstrates the effectiveness and efficacy of the proposed residual error metric for assessing several well-known deep neural network architectures. These results thus illustrate that the proposed measure could be a useful tool for not only assessing the robustness of deep neural networks used in mission-critical scenarios, but also in the design of adversarially robust models.


title: A Coarse-to-Fine Instance Segmentation Network with Learning Boundary Representation

abstract: Boundary-based instance segmentation has drawn much attention since of its attractive efficiency. However, existing methods suffer from the difficulty in long-distance regression. In this paper, we propose a coarse-to-fine module to address the problem. Approximate boundary points are generated at the coarse stage and then features of these points are sampled and fed to a refined regressor for fine prediction. It is end-to-end trainable since differential sampling operation is well supported in the module. Furthermore, we design a holistic boundary-aware branch and introduce instance-agnostic supervision to assist regression. Equipped with ResNet-101, our approach achieves 31.7\% mask AP on COCO dataset with single-scale training and testing, outperforming the baseline 1.3\% mask AP with less than 1\% additional parameters and GFLOPs. Experiments also show that our proposed method achieves competitive performance compared to existing boundary-based methods with a lightweight design and a simple pipeline.


title: Direct Evidence for Magnetic Reflection of Heavy Ions from High Mach Number Collisionless Shocks

abstract: Strong shocks in collisionless plasmas, such as supernovae shocks and shocks driven by coronal mass ejections, are known to be a primary source of energetic particles. Due to their different mass per charge ratio, the interaction of heavy ions with the shock layer differs from that of protons, and injection of these ions into acceleration processes is a challenge. Here we show the first direct observational evidence of magnetic reflection of alpha particles from a high Mach number quasi-perpendicular shock using in-situ spacecraft measurements. The intense magnetic amplification at the shock front associated with nonstationarity modulates the trajectory of alpha particles, some of which travel back upstream as they gyrate in the enhanced magnetic field and experience further acceleration in the upstream region. Our results in particular highlight the important role of high magnetic amplification in seeding heavy ions into the energization processes at nonstationary reforming shocks.


title: Arnold web and dynamical tunneling in a four-site Bose-Hubbard model

abstract: We present the quantum and classical study of a four-well trapped Bose-Einstein condensate (BEC) modeled by the Bose-Hubbard Hamiltonian. The model used is fashioned as a minimal bipartite system consisting of a trimer coupled weakly to a monomer. Semiclassical insights into the fragmentation dynamics of the condensate is obtained by mapping out the Arnold web of the high dimensional classical phase space. A remarkable one-to-one correspondence between the quantum state space in occupation number representation and the classical Arnold web in action space is observed. The relevance of dynamical tunneling for the dynamics of Fock states near resonance junctions on the Arnold web is highlighted. We also predict the possibility of manipulating the transport within the trimer subsystem due to the weak coupling to the monomer.


title: Operator estimates for homogenization of the Robin Laplacian in a perforated domain

abstract: Let $\varepsilon>0$ be a small parameter. We consider the domain $\Omega_\varepsilon:=\Omega\setminus D_\varepsilon$, where $\Omega$ is an open domain in $\mathbb{R}^n$, and $D_\varepsilon$ is a family of small balls of the radius $d_\varepsilon=o(\varepsilon)$ distributed periodically with period $\varepsilon$. Let $\Delta_\varepsilon$ be the Laplace operator in $\Omega_\varepsilon$ subject to the Robin condition ${\partial u\over \partial n}+\gamma_\varepsilon u = 0$ with $\gamma_\varepsilon\ge 0$ on the boundary of the holes and the Dirichlet condition on the exterior boundary. Kaizu (1985, 1989) and Brillard (1988) have shown that, under appropriate assumptions on $d_\varepsilon$ and $\gamma_\varepsilon$, the operator $\Delta_\varepsilon$ converges in the strong resolvent sense to the sum of the Dirichlet Laplacian in $\Omega$ and a constant potential. We improve this result deriving estimates on the rate of convergence in terms of $L^2\to L^2$ and $L^2\to H^1$ operator norms. As a byproduct we establish the estimate on the distance between the spectra of the associated operators.


title: Community Detection in Interval-Weighted Networks

abstract: In this paper we introduce and develop the concept of Interval-Weighted Networks (IWN), a novel approach in Social Network Analysis, where the edge weights are represented by closed intervals composed with precise information, comprehending intrinsic variability. We extend IWN for both Newman's modularity and modularity gain and the Louvain algorithm (LA), considering a tabular representation of networks by contingency tables. We apply our methodology in a real-world commuter network in mainland Portugal between the twenty three NUTS 3 regions. The optimal partition of regions is developed and compared using two new different approaches, designated as ``Classic Louvain'' (CL) and ``Hybrid Louvain'' (HL), which allow taking into account the variability observed in the original network, thereby minimizing the loss of information present in the raw data. Our findings suggest the division of the twenty three Portuguese regions in three main communities. However, we find different geographical partitions according to the community detection methodology used. This analysis can be useful in many real-world applications, since it takes into account that the weights may vary within the ranges, rather than being constant.


title: Charged wormholes in spacetimes of embedding class one

abstract: The existence of charged black holes has suggested that wormholes may also be charged. The purpose of this paper is to construct a general model of a charged wormhole that proves to be a natural extension of the original Morris-Thorne wormhole. This goal is achieved by means of the classical embedding theory that has played a major role in the general theory of relativity.


title: Classes of normally and nearly normally torsion-free monomial ideals

abstract: In this paper, our main focus is to explore different classes of nearly normally torsion-free ideals. We first characterize all finite simple connected graphs with nearly normally torsion-free cover ideals. Next, we characterize all normally torsion-free $t$-spread principal Borel ideals that can also be viewed as edge ideals of uniform multipartite hypergraphs.


title: Semantic navigation with domain knowledge

abstract: Several deployment locations of mobile robotic systems are human made (i.e. urban firefighter, building inspection, property security) and the manager may have access to domain-specific knowledge about the place, which can provide semantic contextual information allowing better reasoning and decision making. In this paper we propose a system that allows a mobile robot to operate in a location-aware and operator-friendly way, by leveraging semantic information from the deployment location and integrating it to the robots localization and navigation systems. We integrate Building Information Models (BIM) into the Robotic Operating System (ROS), to generate topological and metric maps fed to an layered path planner (global and local). A map merging algorithm integrates newly discovered obstacles into the metric map, while a UWB-based localization system detects equipment to be registered back into the semantic database. The results are validated in simulation and real-life deployments in buildings and construction sites.


title: Traceable localization enables accurate integration of quantum emitters and photonic structures with high yield

abstract: Traceability to the International System of Units (SI) is fundamental to measurement accuracy and reliability. In this study, we demonstrate subnanometer traceability of localization microscopy, establishing a metrological foundation for the maturation and application of super-resolution methods. To do so, we create a master standard by measuring the positions of submicrometer apertures in an array by traceable atomic-force microscopy. We perform correlative measurements of this master standard by optical microscopy, calibrating scale factor and correcting aberration effects. We introduce the concept of a localization uncertainty field due to optical localization errors and scale factor uncertainty, with regions of position traceability to within a 68 % coverage interval of +/- 1 nm. These results enable localization metrology with high throughput, which we apply to measure working standards that we fabricate by electron-beam lithography, validating the accuracy of mean pitch and closing the loop for disseminating and integrating reference arrays. We then apply our novel methods to calibrate an optical microscope with a sample cryostat, accounting for thermal contraction by use of a submicrometer pillar array in silicon as a reference material and elucidating complex distortion. This new calibration enables the accurate integration of quantum emitters and photonic structures with high yield, as we demonstrate theoretically through simulations of the dependence of the Purcell factor of radiative enhancement on registration errors across a wide field. Our study illuminates the challenges and opportunities of achieving traceable localization for comparison of position data across lithography and microscopy systems, from ambient to cryogenic temperatures.


title: Oscillatory rheotaxis of active droplets in microchannels

abstract: Biological microswimmers are known to navigate upstream of an external flow (positive rheotaxis) in trajectories ranging from linear, spiral to oscillatory. Such rheotaxis stems from the interplay between the motion and complex shapes of the microswimmers, e.g. the chirality of the rotating flagella, the shear flow characteristics, and the hydrodynamic interaction with a confining surface. Here, we show that an isotropic, active droplet microswimmer exhibits a unique oscillatory rheotaxis in a microchannel despite its simple spherical geometry. The swimming velocity, orientation, and the chemical wake of the active droplet undergo periodic variations between the confining walls during the oscillatory navigation. Using a hydrodynamic model and concepts of dynamical systems, we demonstrate that the oscillatory rheotaxis of the active droplet emerges primarily from the interplay between the hydrodynamic interaction of the finite-sized microswimmer with all the microchannel walls, and the shear flow characteristics. Such oscillatory rheotactic behavior is different from the directed motion near a planar wall observed previously for artificial microswimmers in shear flows. Our results provide a realistic understanding of the behaviour of active particles in confined microflows, as will be encountered in majority of the applications like targeted drug delivery.


title: iCISCF: An Iterative Configuration Interaction-based Multiconfigurational Self-consistent Field Theory for Large Active Spaces

abstract: An iterative configuration interaction (iCI)-based multiconfigurational self-consistent field (SCF) theory, iCISCF, is proposed to handle systems that require large complete active spaces (CAS). The success of iCISCF stems from three ingredients: (1) efficient selection of individual configuration state functions spanning the CAS, meanwhile maintaining full spin symmetry; (2) the use of Jacobi rotation for the optimization of active orbitals, in conjunction with a quasi-Newton algorithm for the core/active-virtual and core-active orbital rotations; (3) a second-order perturbative treatment of the residual space left over by the selection procedure (i.e., iCISCF(2)). Just like selected iCI being a very accurate approximation to CASCI, iCISCF(2) is a very accurate approximation to CASSCF. Several examples that go beyond the capability of CASSCF are taken as showcases to reveal the performances of iCISCF and iCISCF(2).


title: A new approach to proper orthogonal decomposition with difference quotients

abstract: In a recent work [B. Koc et al., arXiv:2010.03750, SIAM J. Numer. Anal., to appear], the authors showed that including difference quotients (DQs) is necessary in order to prove optimal pointwise in time error bounds for proper orthogonal decomposition (POD) reduced order models of the heat equation. In this work, we introduce a new approach to including DQs in the POD procedure. Instead of computing the POD modes using all of the snapshot data and DQs, we only use the first snapshot along with all of the DQs and special POD weights. We show that this approach retains all of the numerical analysis benefits of the standard POD DQ approach, while using a POD data set that has half the number of snapshots as the standard POD DQ approach, i.e., the new approach is more computationally efficient. We illustrate our theoretical results with numerical experiments.


title: Magnetic monopole as the shadow side of the electric charge

abstract: It is shown that a magnetic monopole appears as the tension state of the primary electric charge at its motion through each section of the path equal to the particle's de Broglie wavelength. This conclusion is followed from a submicroscopic consideration of particles and their motion in the framework of the theory of physical space in the form of a tessellattice developed by Michel Bounias and the author. The periodical change of the particle's charged state to its monopole state can easily be introduced in the conventional Maxwell equations and the magnetic monopole automatically shows up in the structure of Maxwell's equations. The monopole is also presented in any quark system as a quark obeys dynamics that are also characterized by the appropriate de Broglie wavelength and hence the electric charge changes periodically to the magnetic monopole. A (anti)neutrino emerges as the typical electron's monopole. When the charged particle becomes the monopole, it also loses its mass (the mass passes to the particle's inerton cloud) and thus the neutrino is a massless particle, or more correctly massless monopole.


title: Dual gauge concepts in electrodynamics, and new limitations on gauge invariance

abstract: Gauge invariance, a core principle in electrodynamics, has two separate meanings. One concept treats the photon as the gauge particle for electrodynamics. It is based on symmetries of the Lagrangian, and requires no mention of electric or magnetic fields. The second concept depends directly on the electric and magnetic fields, and how they can be represented by potential functions that are not unique. A general proof that potentials are more fundamental than fields serves to resolve discrepancies. Physical symmetries, however, are altered by gauge transformations and strongly limit gauge freedom. A new constraint on the form of allowable gauge transformations must be introduced that applies to both gauge concepts.


title: $T\bar{T}+\Lambda_2$ Deformed CFT on the Stretched dS$_3$ Horizon

abstract: The $T\bar{T}+\Lambda_2$ deformation of a Holographic CFT on a cylinder (or torus in Euclidean signature) inhabits the stretched cosmological horizon in dS$_3$ when the deformation parameter is tuned to a particular value. I will describe how this insight allows us to compute the entropy of the stretched cosmic horizon including the logarithmic correction.


title: New generating and counting Functions of prime numbers applied to approximate Chebyschev 2nd class function and the least action principle applied to find non-trivial roots of the Zeta function and to Riemann Hypothesis

abstract: We introduce a new set of prime numbers functions including an exact Generating Function and a Discriminating Function of Prime Numbers neither based on prime number tables nor on algorithms. Instead these functions are defined in terms of ordinary elementary functions, therefore having the advantage of being analytic and readily calculable. Also presented are four applications of our new Prime Numbers Generating Function, namely: obtaining a new analytic formula for counting prime numbers, obtaining an approximant to Euler product function, obtaining an approximant to Riemann Zeta (sigma,tau) function based on our primes discriminating function, an accurate approximant to the Chebyshev function of second class in terms of our primes generating function, and the application of this approximant in sharp estimates related to the validity of the Riemann Hypothesis. We also apply the variational calculus of classical mechanics to obtain the non-trivial roots of Riemann zeta function in the complex plane, in an original and novel approach. A variational test function based on the modulus squared of Riemann function is defined, and then Hamilton Principle is applied to analytically obtain the non-trivial roots of Riemann zeta in a completely original way. We optimize our analytical procedure by defining a more general test function that depends explicitly on the abscissa variable sigma, and present a procedure to find non-trivial roots along the critical line, as demanded by Riemann Hypothesis, thus confirming it. Our method even allows us to define a function, that behaves analogous to the Riemann zeta function, and that even admits a functional type equation.


title: A learned conditional prior for the VAE acoustic space of a TTS system

abstract: Many factors influence speech yielding different renditions of a given sentence. Generative models, such as variational autoencoders (VAEs), capture this variability and allow multiple renditions of the same sentence via sampling. The degree of prosodic variability depends heavily on the prior that is used when sampling. In this paper, we propose a novel method to compute an informative prior for the VAE latent space of a neural text-to-speech (TTS) system. By doing so, we aim to sample with more prosodic variability, while gaining controllability over the latent space's structure. By using as prior the posterior distribution of a secondary VAE, which we condition on a speaker vector, we can sample from the primary VAE taking explicitly the conditioning into account and resulting in samples from a specific region of the latent space for each condition (i.e. speaker). A formal preference test demonstrates significant preference of the proposed approach over standard Conditional VAE. We also provide visualisations of the latent space where well-separated condition-specific clusters appear, as well as ablation studies to better understand the behaviour of the system.


title: CT Image Synthesis Using Weakly Supervised Segmentation and Geometric Inter-Label Relations For COVID Image Analysis

abstract: While medical image segmentation is an important task for computer aided diagnosis, the high expertise requirement for pixelwise manual annotations makes it a challenging and time consuming task. Since conventional data augmentations do not fully represent the underlying distribution of the training set, the trained models have varying performance when tested on images captured from different sources. Most prior work on image synthesis for data augmentation ignore the interleaved geometric relationship between different anatomical labels. We propose improvements over previous GAN-based medical image synthesis methods by learning the relationship between different anatomical labels. We use a weakly supervised segmentation method to obtain pixel level semantic label map of images which is used learn the intrinsic relationship of geometry and shape across semantic labels. Latent space variable sampling results in diverse generated images from a base image and improves robustness. We use the synthetic images from our method to train networks for segmenting COVID-19 infected areas from lung CT images. The proposed method outperforms state-of-the-art segmentation methods on a public dataset. Ablation studies also demonstrate benefits of integrating geometry and diversity.


title: Variational approach to the Schr\"odinger equation with a delta-function potential

abstract: We obtain accurate eigenvalues of the one-dimensional Schr\"odinger equation with a Hamiltonian of the form $H_{g}=H+g\delta (x)$, where $\delta (x)$ is the Dirac delta function. We show that the well known Rayleigh-Ritz variational method is a suitable approach provided that the basis set takes into account the effect of the Dirac delta on the wavefunction.


title: A Benchmarks Library for Extended Parametric Timed Automata

abstract: Parametric timed automata are a powerful formalism for reasoning on concurrent real-time systems with unknown or uncertain timing constants. In order to test the efficiency of new algorithms, a fair set of benchmarks is required. We present an extension of the IMITATOR benchmarks library, that accumulated over the years a number of case studies from academic and industrial contexts. We extend here the library with several dozens of new benchmarks; these benchmarks highlight several new features: liveness properties, extensions of (parametric) timed automata (including stopwatches or multi-rate clocks), and unsolvable toy benchmarks. These latter additions help to emphasize the limits of state-of-the-art parameter synthesis techniques, with the hope to develop new dedicated algorithms in the future.


title: True-pairs of Real Linear Operators and Factorization of Real Polynomials

abstract: A linear operator on a finite dimensional nonzero real vector space may not have an eigenvalue. We define a related notion of a true-pair of a linear operator, and then show that each linear operator on a finite dimensional nonzero real vector space has a true-pair. This is usually proved by using the Fundamental theorem of algebra and Cayley-Hamilton theorem. We construct an inductive proof of this fact without using these theorems. From this we deduce that a polynomial with real coefficients can be written as a product of linear factors and quadratic factors with negative discriminant. It thus gives a proof of the latter fact about polynomials with real coefficients, which does not use complex numbers.


title: Dual-view Molecule Pre-training

abstract: Inspired by its success in natural language processing and computer vision, pre-training has attracted substantial attention in cheminformatics and bioinformatics, especially for molecule based tasks. A molecule can be represented by either a graph (where atoms are connected by bonds) or a SMILES sequence (where depth-first-search is applied to the molecular graph with specific rules). Existing works on molecule pre-training use either graph representations only or SMILES representations only. In this work, we propose to leverage both the representations and design a new pre-training algorithm, dual-view molecule pre-training (briefly, DMP), that can effectively combine the strengths of both types of molecule representations. The model of DMP consists of two branches: a Transformer branch that takes the SMILES sequence of a molecule as input, and a GNN branch that takes a molecular graph as input. The training of DMP contains three tasks: (1) predicting masked tokens in a SMILES sequence by the Transformer branch, (2) predicting masked atoms in a molecular graph by the GNN branch, and (3) maximizing the consistency between the two high-level representations output by the Transformer and GNN branches separately. After pre-training, we can use either the Transformer branch (this one is recommended according to empirical results), the GNN branch, or both for downstream tasks. DMP is tested on nine molecular property prediction tasks and achieves state-of-the-art performances on seven of them. Furthermore, we test DMP on three retrosynthesis tasks and achieve state-of-the-art results on them.


title: The Gravitational-Wave Physics II: Progress

abstract: It has been a half-decade since the first direct detection of gravitational waves, which signifies the coming of the era of the gravitational-wave astronomy and gravitational-wave cosmology. The increasing number of the detected gravitational-wave events has revealed the promising capability of constraining various aspects of cosmology, astronomy, and gravity. Due to the limited space in this review article, we will briefly summarize the recent progress over the past five years, but with a special focus on some of our own work for the Key Project ``Physics associated with the gravitational waves'' supported by the National Natural Science Foundation of China. In particular, (1) we have presented the mechanism of the gravitational-wave production during some physical processes of the early Universe, such as inflation, preheating and phase transition, and the cosmological implications of gravitational-wave measurements; (2) we have put constraints on the neutron star maximum mass according to GW170817 observations; (3) we have developed a numerical relativity algorithm based on the finite element method and a waveform model for the binary black hole coalescence along an eccentric orbit.


title: Efficient Black-Box Importance Sampling for VaR and CVaR Estimation

abstract: This paper considers Importance Sampling (IS) for the estimation of tail risks of a loss defined in terms of a sophisticated object such as a machine learning feature map or a mixed integer linear optimisation formulation. Assuming only black-box access to the loss and the distribution of the underlying random vector, the paper presents an efficient IS algorithm for estimating the Value at Risk and Conditional Value at Risk. The key challenge in any IS procedure, namely, identifying an appropriate change-of-measure, is automated with a self-structuring IS transformation that learns and replicates the concentration properties of the conditional excess from less rare samples. The resulting estimators enjoy asymptotically optimal variance reduction when viewed in the logarithmic scale. Simulation experiments highlight the efficacy and practicality of the proposed scheme


title: An estimate of asymptotics of the moments of additive arithmetic functions with a limit distribution defined on a subset of the natural series

abstract: We study the asymptotics of the moments of arithmetic functions that have a limit distribution, not necessarily normal, defined on a subset of the natural series that satisfies certain requirements. Several assertions are proved on estimating the asymptotics of the moments of strongly additive arithmetic functions and also with additive functions of the class H that have a limit distribution and are defined on a subset of the natural series. The first version of the article is devoted to the study of the asymptotics of the moments of arithmetic functions that have a limit distribution on the natural series. The second version of the article is devoted to the study of the asymptotics of the moments of arithmetic functions that have a limit distribution on an arithmetic progression.


title: Nonparametric Hamiltonian Monte Carlo

abstract: Probabilistic programming uses programs to express generative models whose posterior probability is then computed by built-in inference engines. A challenging goal is to develop general purpose inference algorithms that work out-of-the-box for arbitrary programs in a universal probabilistic programming language (PPL). The densities defined by such programs, which may use stochastic branching and recursion, are (in general) nonparametric, in the sense that they correspond to models on an infinite-dimensional parameter space. However standard inference algorithms, such as the Hamiltonian Monte Carlo (HMC) algorithm, target distributions with a fixed number of parameters. This paper introduces the Nonparametric Hamiltonian Monte Carlo (NP-HMC) algorithm which generalises HMC to nonparametric models. Inputs to NP-HMC are a new class of measurable functions called "tree representable", which serve as a language-independent representation of the density functions of probabilistic programs in a universal PPL. We provide a correctness proof of NP-HMC, and empirically demonstrate significant performance improvements over existing approaches on several nonparametric examples.


title: Minimal and characteristic polynomials of symmetric matrices in characteristic two

abstract: Let $k$ be a field of characteristic two. We prove that a non constant monic polynomial $f\in k[X]$ of degree $n$ is the minimal/characteristic polynomial of a symmetric matrix with entries in $k$ if and only if it is not the product of pairwise distinct inseparable irreducible polynomials. In this case, we prove that $f$ is the minimal polynomial of a symmetric matrix of size $n$. We also prove that any element $\alpha\in k_{alg}$ of degree $n\geq 1$ is the eigenvalue of a symmetrix matrix of size $n$ or $n+1$, the first case happening if and only if the minimal polynomial of $\alpha$ is not the product of pairwise distinct inseparable irreducible polynomials.


title: VSAC: Efficient and Accurate Estimator for H and F

abstract: We present VSAC, a RANSAC-type robust estimator with a number of novelties. It benefits from the introduction of the concept of independent inliers that improves significantly the efficacy of the dominant plane handling and, also, allows near error-free rejection of incorrect models, without false positives. The local optimization process and its application is improved so that it is run on average only once. Further technical improvements include adaptive sequential hypothesis verification and efficient model estimation via Gaussian elimination. Experiments on four standard datasets show that VSAC is significantly faster than all its predecessors and runs on average in 1-2 ms, on a CPU. It is two orders of magnitude faster and yet as precise as MAGSAC++, the currently most accurate estimator of two-view geometry. In the repeated runs on EVD, HPatches, PhotoTourism, and Kusvod2 datasets, it never failed.


title: An Analysis of the Deployment of Models Trained on Private Tabular Synthetic Data: Unexpected Surprises

abstract: Diferentially private (DP) synthetic datasets are a powerful approach for training machine learning models while respecting the privacy of individual data providers. The effect of DP on the fairness of the resulting trained models is not yet well understood. In this contribution, we systematically study the effects of differentially private synthetic data generation on classification. We analyze disparities in model utility and bias caused by the synthetic dataset, measured through algorithmic fairness metrics. Our first set of results show that although there seems to be a clear negative correlation between privacy and utility (the more private, the less accurate) across all data synthesizers we evaluated, more privacy does not necessarily imply more bias. Additionally, we assess the effects of utilizing synthetic datasets for model training and model evaluation. We show that results obtained on synthetic data can misestimate the actual model performance when it is deployed on real data. We hence advocate on the need for defining proper testing protocols in scenarios where differentially private synthetic datasets are utilized for model training and evaluation.


title: Topological gravity in 3+1D and a possible origin of dark matte

abstract: Dark matter is one of the deepest mystery of the universe. So far there is no natural explanation why the dark matter should exist and even dominate the universe. In this paper, we begin with a 3+1D topological gravity theory which is super renormalizable with vanishing beta functions, then we argue that Einstein gravity can emerge by condensing loop-like excitation from the underlying topological gravity theory. In the meanwhile, the uncondensed loop-like excitations serves as a natural candidate of dark matter and a generalized Einstein equation can be derived in the presence of loop-source(dark matter) background. Surprisingly, we find that such kind of dark matter will not contribute to scalar curvature, however, it will become a source of torsion. Finally, we derive the generalized Einstein equation in the presence of Dirac field. Very different from the usual Einstein-Carton theory, our theory further predicts that any type of normal matter, including Dirac field will not produce torsion. All these unique predictions can be tested by future experiments. Our framework suggests that topological invariant principle might play a more profound role than the well-known general covariance principle, especially towards understanding the nature of dark matter and quantum gravity in 3+1D.


title: On self shrinkers of medium entropy in $\mathbb{R}^4$

abstract: In this article we study smooth asymptotically conical self shrinkers in $\mathbb{R}^4$ with Colding-Minicozzi entropy bounded above by $\Lambda_{1}$.


title: Study of neutrino-nucleus reactions with CRISP Program (0 < $E_\nu$ < 3 GeV)

abstract: The neutrino-nucleus reactions are studied at energies from 0 to 3 GeV, using the CRISP program. To simulate these reactions, CRISP uses the Monte Carlo method through an intranuclear cascade model. Quase-elastic and baryonic resonance formation channels for the neutrino-nucleon interaction are considered. The total and differential particle emission cross-sections were obtained, obtaining a good agreement with the values reported by the MiniBooNE experiment. The influence of nuclear effects on the studied reactions, such as fermionic motion, the Pauli blocking mechanism, and the nucleonic separation energy, was shown. It was not possible to simultaneously reproduce the $\nu_\mu + D$ and $\nu_\mu + ^{12}C$ reactions using the same axial mass value. For the charged current quasi-elastic channel, $M_A = 0.95 \ GeV$ for the $\nu_\mu + D$ reaction, and $M_A = 1.35 \ GeV$ for the $\nu_\mu + ^{12}C$ reaction. This can be solved if one considers, in addition to the neutrino-nucleon interaction, the neutrino interaction with a pair of nucleons, just as we demonstrate in the last part of this work.


title: Dynamical implications of convexity beyond dynamical convexity

abstract: We establish sharp dynamical implications of convexity on symmetric spheres that do not follow from dynamical convexity. It allows us to show the existence of elliptic and non-hyperbolic periodic orbits and to furnish new examples of dynamically convex contact forms, in any dimension, that are not equivalent to convex ones via contactomorphisms that preserve the symmetry. Moreover, these examples are $C^1$-stable in the sense that they are actually not equivalent to convex ones via contactomorphisms that are $C^1$-close to those preserving the symmetry. We also show the multiplicity of symmetric non-hyperbolic and symmetric (not necessarily non-hyperbolic) closed Reeb orbits under suitable pinching conditions.


title: Coulomb Drag between a Carbon Nanotube and Monolayer Graphene

abstract: We have measured Coulomb drag between an individual single-walled carbon nanotube (SWNT) as a one-dimensional (1D) conductor and the two-dimensional (2D) conductor monolayer graphene, separated by a few-atom-thick boron nitride layer. The graphene carrier density is tuned across the charge neutrality point (CNP) by a gate, while the SWNT remains degenerate. At high temperatures, the drag resistance changes sign across the CNP, as expected for momentum transfer from drive to drag layer, and exhibits layer exchange Onsager reciprocity. We find that layer reciprocity is broken near the graphene CNP at low temperatures due to nonlinear drag response associated with temperature dependent drag and thermoelectric effects. The drag resistance shows power-law dependences on temperature and carrier density characteristic of 1D Fermi liquid-2D Dirac fluid drag. The 2D drag signal at high temperatures decays with distance from the 1D source slower than expected for a diffusive current distribution, suggesting additional interaction effects in the graphene in the hydrodynamic transport regime.


title: Cross-sections of photonuclear reactions on $^{nat}$Mo targets at end-point bremsstrahlung energy up to $E_{{\gamma}max}$ = 100 MeV

abstract: Experiments to determine the yields and bremsstrahlung flux-averaged cross-sections <${\sigma}(E_{{\gamma}max})$> of photonuclear reactions on the natural Mo targets were performed on the beam from the NSC KIPT electron linear accelerator LUE-40 with the use of the ${\gamma}$-activation technique. The bremsstrahlung end-point energies were in the range $E_{{\gamma}max}$ = 35...80 MeV. The bremsstrahlung quantum flux was calculated with the program GEANT4.9.2 and, in addition, was monitored using the $^{100}$Mo$({\gamma},n)^{99}$Mo reaction. Calculations of the yields and average cross-sections <${\sigma}(E_{{\gamma}max})$> for photonuclear reactions on stable Mo isotopes were computed using the ${\sigma}(E)$ cross-sections from the TALYS1.95 code (for the level density model LD1). A comparison of experimental and calculated cross-sections <${\sigma}(E_{{\gamma}max})$> for reactions $^{92}$Mo$({\gamma},2n)^{90}$Mo and $^{92}$Mo$({\gamma},pn)^{90}$Nb was performed.


title: Existence and Uniqueness of Exact WKB Solutions for Second-Order Singularly Perturbed Linear ODEs

abstract: We prove an existence and uniqueness theorem for exact WKB solutions of general singularly perturbed linear second-order ODEs in the complex domain. These include the one-dimensional time-independent complex Schr\"odinger equation. Notably, our results are valid both in the case of generic WKB trajectories as well as closed WKB trajectories. We also explain in what sense exact and formal WKB solutions form a basis. As a corollary of the proof, we establish the Borel summability of formal WKB solutions for a large class of problems, and derive an explicit formula for the Borel transform.


title: Towards hexagonal $C_{2v}$ systems with anisotropic DMI: Characterization of epitaxial Co$(10\bar{1}0)$/Pt$(110)$ multilayer films

abstract: Ferromagnet/heavy metal (FM/HM) multilayer thin films with $C_{2v}$ symmetry have the potential to host antiskyrmions and other chiral spin textures via an anisotropic Dzyaloshinkii-Moriya interaction (DMI). Here, we present a candidate material system that also has a strong uniaxial magnetocrystalline anisotropy aligned in the plane of the film. This system is based on a new Co/Pt epitaxial relationship, which is the central focus of this work: hexagonal closed-packed Co$(10\bar{1}0)[0001]$ $\parallel$ face-centered cubic Pt$(110)[001]$. We characterized the crystal structure and magnetic properties of our films using X-ray diffraction techniques and magnetometry respectively, including q-scans to determine stacking fault densities and their correlation with the measured magnetocrystalline anisotropy constant and thickness of Co. In future ultrathin multilayer films, we expect this epitaxial relationship to further enable an anisotropic DMI while supporting interfacial perpendicular magnetic anisotropy. The anticipated confluence of these properties, along with the tunability of multilayer films, make this material system a promising testbed for unveiling new spin configurations in FM/HM films.


title: Observation of the Critical State to Multiple-Type Dirac Semimetal Phases in KMgBi

abstract: Dirac semimetals (DSMs) are classified into different phases based on the types of the Dirac fermions. Tuning the transition among different types of the Dirac fermions in one system remains challenging. Recently, KMgBi was predicted to be located at a critical state that various types of Dirac fermions can be induced owing to the existence of a flat band. Here, we carried out systematic studies on the electronic structure of KMgBi single crystal by combining angle-resolve photoemission spectroscopy (ARPES) and scanning tunneling microscopy/spectroscopy (STM/STS). The flat band was clearly observed near the Fermi level. We also revealed a small bandgap of ~ 20 meV between the flat band and the conduction band. These results demonstrate the critical state of KMgBi that transitions among various types of Dirac fermions can be tuned in one system.


title: Active Offline Policy Selection

abstract: This paper addresses the problem of policy selection in domains with abundant logged data, but with a restricted interaction budget. Solving this problem would enable safe evaluation and deployment of offline reinforcement learning policies in industry, robotics, and recommendation domains among others. Several off-policy evaluation (OPE) techniques have been proposed to assess the value of policies using only logged data. However, there is still a big gap between the evaluation by OPE and the full online evaluation. Yet, large amounts of online interactions are often not possible in practice. To overcome this problem, we introduce active offline policy selection - a novel sequential decision approach that combines logged data with online interaction to identify the best policy. We use OPE estimates to warm start the online evaluation. Then, in order to utilize the limited environment interactions wisely we decide which policy to evaluate next based on a Bayesian optimization method with a kernel that represents policy similarity. We use multiple benchmarks, including real-world robotics, with a large number of candidate policies to show that the proposed approach improves upon state-of-the-art OPE estimates and pure online policy evaluation.


title: Less is More: Feature Selection for Adversarial Robustness with Compressive Counter-Adversarial Attacks

abstract: A common observation regarding adversarial attacks is that they mostly give rise to false activation at the penultimate layer to fool the classifier. Assuming that these activation values correspond to certain features of the input, the objective becomes choosing the features that are most useful for classification. Hence, we propose a novel approach to identify the important features by employing counter-adversarial attacks, which highlights the consistency at the penultimate layer with respect to perturbations on input samples. First, we empirically show that there exist a subset of features, classification based in which bridge the gap between the clean and robust accuracy. Second, we propose a simple yet efficient mechanism to identify those features by searching the neighborhood of input sample. We then select features by observing the consistency of the activation values at the penultimate layer.


title: Growth and characterization of large (Y,La)TiO$_3$ and (Y,Ca)TiO$_3$ single crystals

abstract: The Mott-insulating rare-earth titanates (RTiO$_3$, R being a rare-earth ion) are an important class of materials that encompasses interesting spin-orbital phases as well as ferromagnet-antiferromagnet and insulator-metal transitions. The growth of these materials has been plagued by difficulties related to overoxidation, which arises from a strong tendency of Ti$^{3+}$ to oxidize to Ti$^{4+}$. We describe our efforts to grow sizable single crystals of YTiO$_3$ and its La-substituted and Ca-doped variants with the optical travelling-solvent floating-zone technique. We present sample characterization $via$ chemical composition analysis, magnetometry, charge transport, neutron scattering, x-ray absorption spectroscopy and x-ray magnetic circular dichroism to understand macroscopic physical property variations associated with overoxidation. Furthermore, we demonstrate a good signal-to-noise ratio in inelastic magnetic neutron scattering measurements of spin-wave excitations. A superconducting impurity phase, found to appear in Ca-doped samples at high doping levels, is identified as TiO.


title: An Empirical Investigation into Deep and Shallow Rule Learning

abstract: Inductive rule learning is arguably among the most traditional paradigms in machine learning. Although we have seen considerable progress over the years in learning rule-based theories, all state-of-the-art learners still learn descriptions that directly relate the input features to the target concept. In the simplest case, concept learning, this is a disjunctive normal form (DNF) description of the positive class. While it is clear that this is sufficient from a logical point of view because every logical expression can be reduced to an equivalent DNF expression, it could nevertheless be the case that more structured representations, which form deep theories by forming intermediate concepts, could be easier to learn, in very much the same way as deep neural networks are able to outperform shallow networks, even though the latter are also universal function approximators. In this paper, we empirically compare deep and shallow rule learning with a uniform general algorithm, which relies on greedy mini-batch based optimization. Our experiments on both artificial and real-world benchmark data indicate that deep rule networks outperform shallow networks.


title: Minimizing capacity among linear images of rotationally invariant conductors

abstract: Logarithmic capacity is shown to be minimal for a planar set having $N$-fold rotational symmetry ($N \geq 3$), among all conductors obtained from the set by area-preserving linear transformations. Newtonian and Riesz capacities obey a similar property in all dimensions, when suitably normalized linear transformations are applied to a set having irreducible symmetry group. A corollary is P\'{o}lya and Schiffer's lower bound on capacity in terms of moment of inertia.


title: Gravitational Bremsstrahlung and Hidden Supersymmetry of Spinning Bodies

abstract: The recently established formalism of a worldline quantum field theory, which describes the classical scattering of massive bodies in Einstein gravity, is generalized up to quadratic order in spin -- for a pair of Kerr black holes revealing a hidden ${\mathcal N}=2$ supersymmetry. The far-field time-domain waveform of the gravitational waves produced in such a spinning encounter is computed at leading order in the post-Minkowskian (weak field, but generic velocity) expansion, and exhibits this supersymmetry. From the waveform we extract the leading-order total radiated angular momentum in a generic reference frame, and the total radiated energy in the center-of-mass frame to leading order in a low-velocity approximation.


title: Rigorous bounds on the Analytic $S$-matrix

abstract: We consider a dual $S$-matrix Bootstrap approach in $d\geq 3$ space-time dimensions which relies solely on the rigorously proven analyticity, crossing, and unitarity properties of the scattering amplitudes. As a proof of principle, we provide rigorous upper and lower numerical bounds on the quartic coupling for the scattering of identical scalar particles in four dimensions.


title: Bridging the Gap Between Object Detection and User Intent via Query-Modulation

abstract: When interacting with objects through cameras, or pictures, users often have a specific intent. For example, they may want to perform a visual search. With most object detection models relying on image pixels as their sole input, undesired results are not uncommon. Most typically: lack of a high-confidence detection on the object of interest, or detection with a wrong class label. The issue is especially severe when operating capacity-constrained mobile object detectors on-device. In this paper we investigate techniques to modulate mobile detectors to explicitly account for the user intent, expressed as an embedding of a simple query. Compared to standard detectors, query-modulated detectors show superior performance at detecting objects for a given user query. Thanks to large-scale training data synthesized from standard object detection annotations, query-modulated detectors also outperform a specialized referring expression recognition system. Query-modulated detectors can also be trained to simultaneously solve for both localizing a user query and standard detection, even outperforming standard mobile detectors at the canonical COCO task.


title: On-Device Personalization of Automatic Speech Recognition Models for Disordered Speech

abstract: While current state-of-the-art Automatic Speech Recognition (ASR) systems achieve high accuracy on typical speech, they suffer from significant performance degradation on disordered speech and other atypical speech patterns. Personalization of ASR models, a commonly applied solution to this problem, is usually performed in a server-based training environment posing problems around data privacy, delayed model-update times, and communication cost for copying data and models between mobile device and server infrastructure. In this paper, we present an approach to on-device based ASR personalization with very small amounts of speaker-specific data. We test our approach on a diverse set of 100 speakers with disordered speech and find median relative word error rate improvement of 71% with only 50 short utterances required per speaker. When tested on a voice-controlled home automation platform, on-device personalized models show a median task success rate of 81%, compared to only 40% of the unadapted models.


title: A proof of conjecture of Li-Yang

abstract: In this paper, we prove a conjecture posed by Li-Yang in \cite{ly3}. We prove the following result: Let $f(z)$ be a nonconstant entire function, and let $a(z)\not\equiv\infty, b(z)\not\equiv\infty$ be two distinct small meromorphic functions of $f(z)$. If $f(z)$ and $f^{(k)}(z)$ share $a(z)$ and $b(z)$ IM. Then $f(z)\equiv f^{(k)}(z)$, which confirms a conjecture due to Li and Yang (in Illinois J. Math. 44:349-362, 2000).


title: Frank-Wolfe and friends: a journey into projection-free first-order optimization methods

abstract: Invented some 65 years ago in a seminal paper by Marguerite Straus-Frank and Philip Wolfe, the Frank-Wolfe method recently enjoys a remarkable revival, fuelled by the need of fast and reliable first-order optimization methods in Data Science and other relevant application areas. This review tries to explain the success of this approach by illustrating versatility and applicability in a wide range of contexts, combined with an account on recent progress in variants, both improving on the speed and efficiency of this surprisingly simple principle of first-order optimization.


title: A Probabilistic Representation of DNNs: Bridging Mutual Information and Generalization

abstract: Recently, Mutual Information (MI) has attracted attention in bounding the generalization error of Deep Neural Networks (DNNs). However, it is intractable to accurately estimate the MI in DNNs, thus most previous works have to relax the MI bound, which in turn weakens the information theoretic explanation for generalization. To address the limitation, this paper introduces a probabilistic representation of DNNs for accurately estimating the MI. Leveraging the proposed MI estimator, we validate the information theoretic explanation for generalization, and derive a tighter generalization bound than the state-of-the-art relaxations.


title: Pattern selection in the Schnakenberg equations: From normal to anomalous diffusion

abstract: Pattern formation in the classical and fractional Schnakenberg equations is studied to understand the nonlocal effects of anomalous diffusion. Starting with linear stability analysis, we find that if the activator and inhibitor have the same diffusion power, the Turing instability space depends only on the ratio of diffusion coefficients $\kappa_1/\kappa_2$. However, the smaller diffusive powers might introduce larger unstable wave numbers with wider band, implying that the patterns may be more chaotic in the fractional cases. We then apply a weakly nonlinear analysis to predict the parameter regimes for spot, stripe, and mixed patterns in the Turing space. Our numerical simulations confirm the analytical results and demonstrate the differences of normal and anomalous diffusion on pattern formation. We find that in the presence of superdiffusion the patterns exhibit multiscale structures. The smaller the diffusion powers, the larger the unstable wave numbers and the smaller the pattern scales.


title: Lagrangian fields, Calabi functions, and local symplectic groupoids

abstract: A Lagrangian field on a symplectic manifold $M$ is a family $\Lambda=\{\Lambda_x|x \in M\}$ of pointed Lagrangian submanifolds of $M$. This notion is a generalization of a real Lagrangian polarization for which each $\Lambda_x$ is the leaf containing $x$. Two Lagrangian fields $\Lambda$ and $\tilde \Lambda$ are called transversal if $\Lambda_x$ intersects $\tilde\Lambda_x$ transversally at $x$ for every $x$. Two transversal Lagrangian fields determine an almost para-K\"ahler structure on $M$. We construct a local symplectic groupoid on a neighborhood of the zero section of $T^\ast M$ from two transversal Lagrangian fields on $M$. The Lagrangian manifold of $n$-cycles of this groupoid in $(T^\ast M)^n$ has a generating function whose germ around the diagonal of $M^n$ is given by the $n$-point cyclic Calabi function of a closed (1,1)-form on a neighborhood of the diagonal of $M^2$ obtained from the symplectic form on $M$.


title: Liberating Confinement from Lagrangians: 1-form Symmetries and Lines in 4d N=1 from 6d N=(2,0)

abstract: We study confinement in 4d N=1 theories obtained by deforming 4d N=2 theories of Class S. We argue that confinement in a vacuum of the N=1 theory is encoded in the 1-cycles of the associated N=1 curve. This curve is the spectral cover associated to a generalized Hitchin system describing the profiles of two Higgs fields over the Riemann surface upon which the 6d (2,0) theory is compactified. Using our method, we reproduce the expected properties of confinement in various classic examples, such as 4d N=1 pure Super-Yang-Mills theory and the Cachazo-Seiberg-Witten setup. More generally, this work can be viewed as providing tools for probing confinement in non-Lagrangian N=1 theories, which we illustrate by constructing an infinite class of non-Lagrangian N=1 theories that contain confining vacua. The simplest model in this class is an N=1 deformation of the N=2 theory obtained by gauging $SU(3)^3$ flavor symmetry of the $E_6$ Minahan-Nemeschansky theory.


title: Perturbed recombination from inhomogeneous photon injection and application to accreting primordial black holes

abstract: Exotic electromagnetic energy injection in the early Universe may alter cosmological recombination, and ultimately cosmic microwave background (CMB) anisotropies. Moreover, if energy injection is inhomogeneous, it may induce a spatially-varying ionization fraction, and non-Gaussianity in the CMB. The observability of these signals, however, is contingent upon how far the injected particles propagate and deposit their energy into the primordial plasma, relative to the characteristic scale of energy injection fluctuations. In this study we inspect the spatial properties of energy deposition and perturbed recombination resulting from an inhomogeneous energy injection of sub-10 MeV photons, relevant to accreting primordial black holes (PBHs). We develop a novel Monte-Carlo radiation transport code accounting for all relevant photon interactions in this energy range, and including secondary electron energy deposition efficiency through a new analytic approximation. For a specified injected photon spectrum, the code outputs an injection-to-deposition Green's function depending on time and distance from the injection point. Combining this output with a linearized solution of the perturbed recombination problem, we derive time- and scale-dependent deposition-to-ionization Green's functions. We apply this general framework to accreting PBHs, whose luminosity is strongly spatially modulated by supersonic relative velocities between cold dark matter and baryons. We find that the resulting spatial fluctuations of the free-electron fraction are of the same magnitude as its mean deviation from standard recombination, from which current CMB power spectra constraints are derived. This work suggests that the sensitivity to accreting PBHs might be substantially improved by propagating these inhomogeneities to CMB anisotropy power spectra and non-Gaussian statistics, which we study in subsequent papers.


title: Location and topology of the fundamental gap in photonic crystals

abstract: The fundamental, or first, band gap is of unmatched importance in the study of photonic crystals. Here, we address precisely where this gap can be opened in the band structure of three-dimensional photonic crystals. Although strongly constrained by symmetry, this problem cannot be addressed directly with conventional band-symmetry analysis due to the existence of a photonic polarization vortex at zero frequency. We develop an approach for overcoming the associated symmetry singularity by incorporating fictitious, auxiliary longitudinal modes. Our strategy also enables us to extend recent developments in symmetry-based topological analysis to the fundamental gap of three-dimensional photonic crystals. Exploiting this, we systematically study the topology of the minimal fundamental gaps. This reveals the existence of topological gap-obstructions that push the fundamental gap higher than what a conventional analysis would suggest. Our work demonstrates that topology can play a crucial role in the opening of the fundamental photonic gap and informs future theoretical and experimental searches for conventional and topological band gaps in three-dimensional photonic crystals.


title: MADE: Exploration via Maximizing Deviation from Explored Regions

abstract: In online reinforcement learning (RL), efficient exploration remains particularly challenging in high-dimensional environments with sparse rewards. In low-dimensional environments, where tabular parameterization is possible, count-based upper confidence bound (UCB) exploration methods achieve minimax near-optimal rates. However, it remains unclear how to efficiently implement UCB in realistic RL tasks that involve non-linear function approximation. To address this, we propose a new exploration approach via \textit{maximizing} the deviation of the occupancy of the next policy from the explored regions. We add this term as an adaptive regularizer to the standard RL objective to balance exploration vs. exploitation. We pair the new objective with a provably convergent algorithm, giving rise to a new intrinsic reward that adjusts existing bonuses. The proposed intrinsic reward is easy to implement and combine with other existing RL algorithms to conduct exploration. As a proof of concept, we evaluate the new intrinsic reward on tabular examples across a variety of model-based and model-free algorithms, showing improvements over count-only exploration strategies. When tested on navigation and locomotion tasks from MiniGrid and DeepMind Control Suite benchmarks, our approach significantly improves sample efficiency over state-of-the-art methods. Our code is available at https://github.com/tianjunz/MADE.


title: Shell model analysis of the $ B(E2,2^+ \rightarrow 0^+)$'s in the A=70 T=1 triplet

abstract: he $B(E2,2^+ \rightarrow 0^+)$ transition strengths of the T=1 isobaric triplet $^{70}$Kr, $^{70}$Br, $^{70}$Se, recently measured at RIKEN/RIBF, are discussed in terms of state of the art large scale shell model calculations using the JUN45 and JUN45+LNPS plus Coulomb interactions. In this letter we argue that, depending on the effective charges used, the calculations are either in line with the experimental data within statistical uncertainties, or the anomaly happens in $^{70}$Br, rather than $^{70}$Kr. In the latter case, we suggest that it can be due to the presence of a hitherto undetected 1$^+$ T=0 state below the yrast 2$^+$ T=1 state. Our results do not support a shape change of $^{70}$Kr with respect to the other members of the isobaric multiplet.


title: How to train your ViT? Data, Augmentation, and Regularization in Vision Transformers

abstract: Vision Transformers (ViT) have been shown to attain highly competitive performance for a wide range of vision applications, such as image classification, object detection and semantic image segmentation. In comparison to convolutional neural networks, the Vision Transformer's weaker inductive bias is generally found to cause an increased reliance on model regularization or data augmentation ("AugReg" for short) when training on smaller training datasets. We conduct a systematic empirical study in order to better understand the interplay between the amount of training data, AugReg, model size and compute budget. As one result of this study we find that the combination of increased compute and AugReg can yield models with the same performance as models trained on an order of magnitude more training data: we train ViT models of various sizes on the public ImageNet-21k dataset which either match or outperform their counterparts trained on the larger, but not publicly available JFT-300M dataset.


title: End-to-end Temporal Action Detection with Transformer

abstract: Temporal action detection (TAD) aims to determine the semantic label and the temporal interval of every action instance in an untrimmed video. It is a fundamental and challenging task in video understanding. Previous methods tackle this task with complicated pipelines. They often need to train multiple networks and involve hand-designed operations, such as non-maximal suppression and anchor generation, which limit the flexibility and prevent end-to-end learning. In this paper, we propose an end-to-end Transformer-based method for TAD, termed TadTR. Given a small set of learnable embeddings called action queries, TadTR adaptively extracts temporal context information from the video for each query and directly predicts action instances with the context. To adapt Transformer to TAD, we propose three improvements to enhance its locality awareness. The core is a temporal deformable attention module that selectively attends to a sparse set of key snippets in a video. A segment refinement mechanism and an actionness regression head are designed to refine the boundaries and confidence of the predicted instances, respectively. With such a simple pipeline, TadTR requires lower computation cost than previous detectors, while preserving remarkable performance. As a self-contained detector, it achieves state-of-the-art performance on THUMOS14 (56.7% mAP) and HACS Segments (32.09% mAP). Combined with an extra action classifier, it obtains 36.75% mAP on ActivityNet-1.3. Code is available at https://github.com/xlliu7/TadTR.


title: Riemannian Convex Potential Maps

abstract: Modeling distributions on Riemannian manifolds is a crucial component in understanding non-Euclidean data that arises, e.g., in physics and geology. The budding approaches in this space are limited by representational and computational tradeoffs. We propose and study a class of flows that uses convex potentials from Riemannian optimal transport. These are universal and can model distributions on any compact Riemannian manifold without requiring domain knowledge of the manifold to be integrated into the architecture. We demonstrate that these flows can model standard distributions on spheres, and tori, on synthetic and geological data. Our source code is freely available online at http://github.com/facebookresearch/rcpm


title: Impact of Image Persistence in the Roman Space Telescope High-Latitude Survey

abstract: The High Latitude Survey of the Nancy Grace Roman Space Telescope is expected to measure the positions and shapes of hundreds of millions of galaxies in an area of 2220 deg$^2$. This survey will provide high-quality weak lensing data with unprecedented systematics control. The Roman Space Telescope will survey the sky in near infrared (NIR) bands using Teledyne H4RG HgCdTe photodiode arrays. These NIR arrays exhibit an effect called persistence: charges that are trapped in the photodiodes during earlier exposures are gradually released into later exposures, leading to contamination of the images and potentially to errors in measured galaxy properties such as fluxes and shapes. In this work, we use image simulations that incorporate the persistence effect to study its impact on galaxy shape measurements and weak lensing signals. No significant spatial correlations are found between the galaxy shape changes induced by persistence. On the scales of interest for weak lensing cosmology, the effect of persistence on the weak lensing correlation function is about two orders of magnitude lower than the Roman Space Telescope additive shear error budget, indicating that the persistence effect is expected to be a subdominant contributor to the systematic error budget for weak lensing with the Roman Space Telescope given its current design.


title: AdS Bulk Locality from Sharp CFT Bounds

abstract: It is a long-standing conjecture that any CFT with a large central charge and a large gap $\Delta_{\text{gap}}$ in the spectrum of higher-spin single-trace operators must be dual to a local effective field theory in AdS. We prove a sharp form of this conjecture by deriving numerical bounds on bulk Wilson coefficients in terms of $\Delta_{\text{gap}}$ using the conformal bootstrap. Our bounds exhibit the scaling in $\Delta_{\text{gap}}$ expected from dimensional analysis in the bulk. Our main tools are dispersive sum rules that provide a dictionary between CFT dispersion relations and S-matrix dispersion relations in appropriate limits. This dictionary allows us to apply recently-developed flat-space methods to construct positive CFT functionals. We show how AdS$_{4}$ naturally resolves the infrared divergences present in 4D flat-space bounds. Our results imply the validity of twice-subtracted dispersion relations for any S-matrix arising from the flat-space limit of AdS/CFT.


title: A Systematic Review of Computational Thinking in Early Ages

abstract: Nowadays, technology has become dominant in the daily lives of most people around the world. From children to older people, technology is present, helping in the most diverse daily tasks and allowing accessibility. However, many times these people are just end-users, without any incentive to the development of computational thinking (CT). With advances in technologies, the abstraction of coding, programming languages, and the hardware resources involved will become a reality. However, while we have not progressed to this stage, it is necessary to encourage the development of CT teaching from an early age. This work will present state of the art concerning teaching initiatives and tools on programming (e.g., ScratchJr), robotics (e.g., KIBO), and other playful tools (e.g., Happy Maps) for the development of CT in the early ages, specifically filling the gap of CT at the kindergarten level. This survey presents a systematic review of the literature, emphasizing computational and robotic tools used in preschool classes to develop the CT. The systematic review evaluated more than 60 papers from 2010 to December 2020, electing 31 papers and adding three papers from the qualitative stage. The paper's amount was classified in taxonomy to show CT's principal tools and initiates applied to children early. To conclude this survey, an extensive discussion about the terms and authors related to this research area is present.


title: Filling-Enforced Obstructed Atomic Insulators

abstract: Topological band theory has achieved great success in the high-throughput search for topological band structures both in paramagnetic and magnetic crystal materials. However, a significant proportion of materials are topologically trivial insulators at the Fermi level. In this paper, we show that, remarkably, for a subset of the topologically trivial insulators, knowing only their electron number and the Wyckoff positions of the atoms we can separate them into two groups: the obstructed atomic insulator (OAI) and the atomic insulator (AI). The interesting group, the OAI, have a center of charge not localized on the atoms. Using the theory of topological quantum chemistry, in this work we first derive the necessary and sufficient conditions for a topologically trivial insulator to be a filling enforced obstructed atomic insulator (feOAI) in the 1651 Shubnikov space groups. Remarkably, the filling enforced criteria enable the identification of obstructed atomic bands without knowing the representations of the band structures. Hence, no ab-initio calculations are needed for the filling enforced criteria, although they are needed to obtain the band gaps. With the help of the Topological Quantum Chemistry website, we have performed a high-throughput search for feOAIs and have found that 957 ICSD entries (638 unique materials) are paramagnetic feOAIs, among which 738 (475) materials have an indirect gap. The metallic obstructed surface states of feOAIs are also showcased by several material examples.


title: GPLA-12: An Acoustic Signal Dataset of Gas Pipeline Leakage

abstract: In this paper, we introduce a new acoustic leakage dataset of gas pipelines, called as GPLA-12, which has 12 categories over 684 training/testing acoustic signals. Unlike massive image and voice datasets, there have relatively few acoustic signal datasets, especially for engineering fault detection. In order to enhance the development of fault diagnosis, we collect acoustic leakage signals on the basis of an intact gas pipe system with external artificial leakages, and then preprocess the collected data with structured tailoring which are turned into GPLA-12. GPLA-12 dedicates to serve as a feature learning dataset for time-series tasks and classifications. To further understand the dataset, we train both shadow and deep learning algorithms to observe the performance. The dataset as well as the pretrained models have been released at both www.daip.club and github.com/Deep-AI-Application-DAIP


title: Efficient Computation of $N$-point Correlation Functions in $D$ Dimensions

abstract: We present efficient algorithms for computing the $N$-point correlation functions (NPCFs) of random fields in arbitrary $D$-dimensional homogeneous and isotropic spaces. Such statistics appear throughout the physical sciences, and provide a natural tool to describe stochastic processes. algorithms for computing the NPCF components have $\mathcal{O}(n^N)$ complexity (for a data set containing $n$ particles); their application is thus computationally infeasible unless $N$ is small. By projecting the statistic onto a suitably-defined angular basis, we show that the estimators can be written in a separable form, with complexity $\mathcal{O}(n^2)$, or $\mathcal{O}(n_{\rm g}\log n_{\rm g})$ if evaluated using a Fast Fourier Transform on a grid of size $n_{\rm g}$. Our decomposition is built upon the $D$-dimensional hyperspherical harmonics; these form a complete basis on the $(D-1)$-sphere and are intrinsically related to angular momentum operators. Concatenation of $(N-1)$ such harmonics gives states of definite combined angular momentum, forming a natural separable basis for the NPCF. As $N$ and $D$ grow, the number of basis components quickly becomes large, providing a practical limitation to this (and all other) approaches: however, the dimensionality is greatly reduced in the presence of symmetries; for example, isotropic correlation functions require only states of zero combined angular momentum. We provide a \textsc{Julia} package implementing our estimators, and show how they can be applied to a variety of scenarios within cosmology and fluid dynamics. The efficiency of such estimators will allow higher-order correlators to become a standard tool in the analysis of random fields.

categories: astro-ph.IM physics.app-ph physics.comp-ph physics.data-an


title: MadFlow: automating Monte Carlo simulation on GPU for particle physics processes

abstract: We present MadFlow, a first general multi-purpose framework for Monte Carlo (MC) event simulation of particle physics processes designed to take full advantage of hardware accelerators, in particular, graphics processing units (GPUs). The automation process of generating all the required components for MC simulation of a generic physics process and its deployment on hardware accelerator is still a big challenge nowadays. In order to solve this challenge, we design a workflow and code library which provides to the user the possibility to simulate custom processes through the MadGraph5_aMC@NLO framework and a plugin for the generation and exporting of specialized code in a GPU-like format. The exported code includes analytic expressions for matrix elements and phase space. The simulation is performed using the VegasFlow and PDFFlow libraries which deploy automatically the full simulation on systems with different hardware acceleration capabilities, such as multi-threading CPU, single-GPU and multi-GPU setups. The package also provides an asynchronous unweighted events procedure to store simulation results. Crucially, although only Leading Order is automatized, the library provides all ingredients necessary to build full complex Monte Carlo simulators in a modern, extensible and maintainable way. We show simulation results at leading-order for multiple processes on different hardware configurations.


title: High star cluster formation efficiency in the strongly lensed Sunburst Lyman-continuum galaxy at z=2.37

abstract: We investigate the strongly lensed (\mu x10-100) Lyman continuum (LyC) galaxy, dubbed Sunburst, at z=2.37, taking advantage of a new accurate model of the lens. A characterization of the intrinsic (delensed) properties of the galaxy yields a size of ~3 sq.kpc, a luminosity Muv=-20.3,and a stellar mass M~10^9 Msun;16% of the ultraviolet light is located in a 3 Myr old gravitationally-bound young massive star cluster (YMC) with an effective radius of Re~8 pc and a dynamical mass of ~10^7 Msun (similar to the stellar mass), from which LyC radiation is detected (\lambda < 912A). The inferred outflowing gas velocity (>300 km/s) exceeds the escape velocity of the star cluster. The resulting escape fraction of the ionizing radiation emerging from the Sunburst galaxy is >6-12%, whilst it is >46-93% if inferred from the YMC. 12 additional likely star clusters with 3<Re<20 pc are identified in the galaxy from which we derive a cluster formation efficiency \Gamma>~30%, which is consistent with the high \Gamma derived in local galaxies experiencing extreme gas physical conditions. The presence of the YMC influences the morphology (nucleation), photometry (photometric jumps) and spectroscopic output (nebular emission) of the entire galaxy. The de-lensed LyC and UV (1600A) magnitudes of the YMC are ~30.6 and ~26.9, whilst the galaxy has m1600~24.8. A relatively large rest-frame equivalent width of EWrest(Hb+[OIII]4959-5007)~450A emerges from the galaxy with the YMC contributing to ~30%. If O-type stars are mainly forged in star clusters, then such engines were the key ionizing agents during reionization and the increasing occurrence of high EW lines (Hb+[OIII]) observed at z>6.5 might be an indirect signature of a high \Gamma at reionization.Future facilities (like VLT/MAVIS or ELT), will probe bound clusters on moderately magnified (\mu<5-10) galaxies across cosmic epochs up to reionization[ABRIDGED]


title: Say Their Names: Resurgence in the collective attention toward Black victims of fatal police violence following the death of George Floyd

abstract: The murder of George Floyd by police in May 2020 sparked international protests and renewed attention in the Black Lives Matter movement. Here, we characterize ways in which the online activity following George Floyd's death was unparalleled in its volume and intensity, including setting records for activity on Twitter, prompting the saddest day in the platform's history, and causing George Floyd's name to appear among the ten most frequently used phrases in a day, where he is the only individual to have ever received that level of attention who was not known to the public earlier that same week. Further, we find this attention extended beyond George Floyd and that more Black victims of fatal police violence received attention following his death than during other past moments in Black Lives Matter's history. We place that attention within the context of prior online racial justice activism by showing how the names of Black victims of police violence have been lifted and memorialized over the last 12 years on Twitter. Our results suggest that the 2020 wave of attention to the Black Lives Matter movement centered past instances of police violence in an unprecedented way, demonstrating the impact of the movement's rhetorical strategy to "say their names."


title: Horizon Acoustics of the GHS Black Hole and the Spectrum of ${\rm AdS}_2$

abstract: We uncover a novel structure in Einstein-Maxwell-dilaton gravity: an ${\rm AdS}_2 \times S^2$ solution in string frame, which can be obtained by a near-horizon limit of the extreme GHS black hole with dilaton coupling $\lambda \neq 1$. Unlike the Bertotti-Robinson spacetime, our solution has independent length scales for the ${\rm AdS}_2$ and $S^2$, with ratio controlled by $\lambda$. We solve the perturbation problem for this solution, finding the independently propagating towers of states in terms of superpositions of gravitons, photons, and dilatons and their associated effective potentials. These potentials describe modes obeying conformal quantum mechanics, with couplings that we compute, and can be recast as giving the spectrum of the effective masses of the modes. By dictating the conformal weights of boundary operators, this spectrum provides crucial data for any future construction of a holographic dual to these ${\rm AdS}_2\times S^2$ configurations.


title: A disk and no signatures of tidal distortion in the galaxy "lacking" dark matter NGC 1052-DF2

abstract: Using ultra-deep imaging ($\mu_g = 30.4$ mag/arcsec$^2$; 3$\sigma$, 10"x10"), we probed the surroundings of the first galaxy "lacking" dark matter KKS2000[04] (NGC 1052-DF2). Signs of tidal stripping in this galaxy would explain its claimed low content of dark matter. However, we find no evidence of tidal tails. In fact, the galaxy remains undisturbed down to a radial distance of 80 arcsec. This radial distance triples previous spatial explorations of the stellar distribution of this galaxy. In addition, the distribution of its globular clusters (GCs) is not extended in relation to the bulk of the galaxy (the radius containing half of the GCs is 21 arcsec). We also found that the surface brightness radial profiles of this galaxy in the g and r bands decline exponentially from 35 to 80 arcsec. That, together with a constant ellipticity and position angle in the outer parts of the galaxy strongly suggests the presence of a low-inclination disk. This is consistent with the evidence of rotation found for this object. This finding implies that the dynamical mass of this galaxy is a factor of 2 higher than previously reported, bringing the dark matter content of this galaxy in line with galaxies of similar stellar mass.


title: The one-loop tadpole in the geoSMEFT

abstract: Making use of the geometric formulation of the Standard Model Effective Field Theory we calculate the one-loop tadpole diagrams to all orders in the Standard Model Effective Field Theory power counting. This work represents the first calculation of a one-loop amplitude beyond leading order in the Standard Model Effective Field Theory, and discusses the potential to extend this methodology to perform similar calculations of observables in the near future.


title: Numerical Analysis of Long-term Variability of AGN Jets through RMHD Simulations

abstract: Relativistic AGN jets exhibit multi-timescale variability and a broadband non-thermal spectrum extending from radio to gamma-rays. These highly magnetized jets are prone to undergo several Magneto-hydrodynamic (MHD) instabilities during their propagation in space and could trigger jet radiation and particle acceleration. This work aims to study the implications of relativistic kink mode instability on the observed long-term variability in the context of the twisting in-homogeneous jet model. To achieve this, we investigate the physical configurations preferable for forming kink mode instability by performing high-resolution 3D relativistic MHD simulations of a portion of highly magnetized jets. In particular, we perform simulations of cylindrical plasma column with Lorentz factor $\geq 5$ and study the effects of magnetization values and axial wave-numbers with decreasing pitch on the onset and growth of kink instability. We have confirmed the impact of axial wave-number on the dynamics of the plasma column including the growth of the instability. In this work, we have further investigated the connection between the dynamics of the plasma column with its time-varying emission features. From our analysis, we find a correlated trend between the growth rate of kink mode instability and the flux variability obtained from the simulated light curve.


title: A strong lensing model of the galaxy cluster PSZ1 G311.65-18.48

abstract: We present a strong lensing analysis of the galaxy cluster PSZ1 G311.65-18.48 (z=0.443) using multi-band observations with Hubble Space Telescope, complemented with VLT/MUSE spectroscopic data. The MUSE observations provide redshift estimates for the lensed sources and help reducing the mis-identification of the multiple images. Spectroscopic data are also used to measure the inner velocity dispersions of 15 cluster galaxies and calibrate the scaling relations to model the subhalo cluster component. The model is based on 62 multiple images grouped in 17 families belonging to 4 different sources. The majority of them are multiple images of compact stellar knots belonging to a single star-forming galaxy at z=2.3702. This source is strongly lensed by the cluster to form the Sunburst Arc system. To accurately reproduce all the multiple images, we build a parametric mass model, which includes both cluster-scale and galaxy-scale components. The resulting model has a r.m.s. separation between the model-predicted and the observed positions of the multiple images of only 0.14''. We conclude that PSZ1 G311.65-18.48 has a relatively round projected shape and a large Einstein radius (29'' for z_s = 2.3702), which could indicate that the cluster is elongated along the line of sight. The Sunburst Arc source is located at the intersection of a complex network of caustics, which explains why parts of the arc are imaged with unprecedented multiplicity (up to 12 times).


title: Free fermion entanglement with a semitransparent interface: the effect of graybody factors on entanglement islands

abstract: We study the entanglement entropy of free fermions in 2d in the presence of a partially transmitting interface that splits Minkowski space into two half-spaces. We focus on the case of a single interval that straddles the defect, and compute its entanglement entropy in three limits: Perturbing away from the fully transmitting and fully reflecting cases, and perturbing in the amount of asymmetry of the interval about the defect. Using these results within the setup of the Poincar\'e patch of AdS$_2$ statically coupled to a zero temperature flat space bath, we calculate the effect of a partially transmitting AdS$_2$ boundary on the location of the entanglement island region. The partially transmitting boundary is a toy model for black hole graybody factors. Our results indicate that the entanglement island region behaves in a monotonic fashion as a function of the transmission/reflection coefficient at the interface.


title: Demonstration of Confinement and Chiral Symmetry Breaking in $SO(N_c)$ Gauge Theories

abstract: We demonstrate that $SO(N_{c})$ gauge theories with matter fields in the vector representation confine due to monopole condensation and break the $SU(N_{F})$ chiral symmetry to $SO(N_{F})$ via the quark bilinear. Our results are obtained by perturbing the ${\cal N}=1$ supersymmetric theory with anomaly-mediated supersymmetry breaking.


title: An investigation of the $\alpha_S$ and heavy quark mass dependence in the MSHT20 global PDF analysis

abstract: We investigate the MSHT20 global PDF sets, demonstrating the effects of varying the strong coupling $\alpha_S(M_Z^2)$ and the masses of the charm and bottom quarks. We determine the preferred value, and accompanying uncertainties, when we allow $\alpha_S(M_Z^2)$ to be a free parameter in the MSHT20 global analyses of deep-inelastic and related hard scattering data, at both NLO and NNLO in QCD perturbation theory. We also study the constraints on $\alpha_S(M_Z^2)$ which come from the individual data sets in the global fit by repeating the NNLO and NLO global analyses at various fixed values of $\alpha_S(M_Z^2)$, spanning the range $\alpha_S(M_Z^2)=0.108$ to $0.130$ in units of $0.001$. We make all resulting PDFs sets available. We find that the best fit values are $\alpha_S(M_Z^2)=0.1203\pm 0.0015$ and $0.1174\pm 0.0013$ at NLO and NNLO respectively. We investigate the relationship between the variations in $\alpha_S(M_Z^2)$ and the uncertainties on the PDFs, and illustrate this by calculating the cross sections for key processes at the LHC. We also perform fits where we allow the heavy quark masses $m_c$ and $m_b$ to vary away from their default values and make PDF sets available in steps of $\Delta m_c =0.05~{\rm GeV}$ and $\Delta m_b =0.25~{\rm GeV}$, using the pole mass definition of the quark masses. As for varying $\alpha_S(M_Z^2)$ values, we present the variation in the PDFs and in the predictions. We examine the comparison to data, particularly the HERA data on charm and bottom cross sections and note that our default values are very largely compatible with best fits to data. We provide PDF sets with 3 and 4 active quark flavours, as well as the standard value of 5 flavours.


title: Classification of singularities of cluster algebras of finite type: the case of trivial coefficients

abstract: We provide a complete classification of the singularities of cluster algebras of finite type with trivial coefficients. Alongside, we develop a constructive desingularization of these singularities via blowups in regular centers over fields of arbitrary characteristic. Furthermore, from the same perspective, we study a family of cluster algebras which are not of finite type and which arise from a star shaped quiver.


title: EFPE: Efficient fully precessing eccentric gravitational waveforms for binaries with long inspirals

abstract: In this paper, we derive a set of equations of motions for binaries on eccentric orbits undergoing spin-induced precession that can efficiently be integrated on the radiation-reaction timescale. We find a family of solutions with a computation cost improved by a factor $10$ - $50$ down to $\sim 10$ ms per waveform evaluation compared to waveforms obtained by directly integrating the precession equations, that maintain a mismatch of the order $10^{-4}$ - $10^{-6}$ for waveforms lasting a million orbital cycles and a thousand spin-induced precession cycles. We express it in terms of parameters that make the solution regular in the equal-mass limit, thus bypassing a problem of previous similar solutions. We point to ways in which the solution presented in this paper can be perturbed to take into account effects such as general quadrupole momenta and post-Newtonian corrections to the precession equations. This new waveform, with its improved efficiency and its accuracy, makes possible Bayesian parameter estimation using the full spin and eccentricity parameter volume for long lasting inspiralling signals such as stellar-origin black hole binaries observed by LISA.


title: Detection of Aerosols at Microbar Pressures in an Exoplanet Atmosphere

abstract: Formation of hazes at microbar pressures has been explored by theoretical models of exoplanet atmospheres to explain Rayleigh scattering and/or featureless transmission spectra, however observational evidence of aerosols in the low pressure formation environments has proved elusive. Here, we show direct evidence of aerosols existing at $\sim$1 microbar pressures in the atmosphere of the warm sub-Saturn WASP-69b using observations taken with Space Telescope Imaging Spectrograph (STIS) and Wide Field Camera 3 (WFC3) instruments on the Hubble Space Telescope. The transmission spectrum shows a wavelength-dependent slope induced by aerosol scattering that covers 11 scale heights of spectral modulation. Drawing on the extensive studies of haze in our Solar System, we model the transmission spectrum based on a scaled version of Jupiter's haze density profile to show that WASP-69b transmission spectrum can be produced by scattering from an approximately constant density of particles extending throughout the atmospheric column from 40 millibar to microbar pressures. These results are consistent with theoretical expectations based on microphysics of the aerosol particles that have suggested haze can exist at microbar pressures in exoplanet atmospheres.


title: Planetary magnetosphere evolution around post-main-sequence stars

abstract: Accompanying the mounting detections of planets orbiting white dwarfs and giant stars are questions about their physical history and evolution, particularly regarding detectability of their atmospheres and potential for habitability. Here we determine how the size of planetary magnetospheres evolve over time from the end of the main sequence through to the white dwarf phase due to the violent winds of red giant and asymptotic giant branch stars. By using a semianalytic prescription, we investigate the entire relevant phase space of planet type, planet orbit and stellar host mass (1-7Msun). We find that a planetary magnetosphere will always be quashed at some point during the giant branch phases unless the planet's magnetic field strength is at least two orders of magnitude higher than Jupiter's current value. We also show that the time variation of the stellar wind and density generates a net increase in wind ram pressure and does not allow a magnetosphere to be maintained at any time for field strengths less than 10^{-5} T (0.1 G). This lack of protection hints that currently potentially habitable planets orbiting white dwarfs would have been previously inhospitable.


title: Irradiation-driven escape of primordial planetary atmospheres I. The ATES photoionization hydrodynamics code

abstract: Intense X-ray and ultraviolet stellar irradiation can heat and inflate the atmospheres of closely orbiting exoplanets, driving mass outflows that may be significant enough to evaporate a sizable fraction of the planet atmosphere over the system lifetime. The recent surge in the number of known exoplanets, together with the imminent deployment of new ground and space-based facilities for exoplanet discovery and characterization, requires a prompt and efficient assessment of the most promising targets for intensive spectroscopic follow-ups. To this purpose, we developed ATES (ATmospheric EScape); a new hydrodynamics code that is specifically designed to compute the temperature, density, velocity and ionization fraction profiles of highly irradiated planetary atmospheres, along with the current, steady-state mass loss rate. ATES solves the one-dimensional Euler, mass and energy conservation equations in radial coordinates through a finite-volume scheme. The hydrodynamics module is paired with a photoionization equilibrium solver that includes cooling via bremsstrahlung, recombination and collisional excitation/ionization for the case of a primordial atmosphere entirely composed of atomic hydrogen and helium, whilst also accounting for advection of the different ion species. Compared against the results of 14 moderately-to-highly irradiated planets simulated with The PLUTO-CLOUDY Interface (TPCI), ATES yields remarkably good agreement at a significantly smaller fraction of the computational time. A convergence study shows that ATES recovers stable, steady-state hydrodynamic solutions for systems with $\log(-\phi_p) \lesssim 12.9 + 0.17\log F_{\rm XUV}$. Incidentally, atmospheres of systems above this threshold are generally thought to be undergoing Jeans escape. The code, which also features a user-friendly graphic interface, is available publicly as an online repository.