title: Recurrent Solar Energetic Particle Flux Enhancements Observed near Earth and Mars

abstract: August 1 to November 15, 2016 period was characterized by the presence of Corotating Interaction Regions (CIRs) and a few weak Coronal Mass Ejections (CMEs) in the heliosphere. In this study we show recurrent energetic electron and proton enhancements observed near Earth (1 AU) and Mars (1.43-1.38 AU) during this period. The observations near Earth are using data from instruments aboard ACE, SOHO, and SDO whereas those near Mars are by the SEP, SWIA, and MAG instruments aboard MAVEN. During this period, the energetic electron fluxes observed near Earth and Mars showed prominent periodic enhancements over four solar rotations, with major periodicities of ~27 days and ~13 days. Periodic radar blackout/weakening of radar signals at Mars are observed by MARSIS/MEX, associated with these solar energetic electron enhancements. During this period, a weak CME and a High Speed Stream (HSS)-related interplanetary shock could interact with the CIR and enhance energetic proton fluxes near 1.43-1.38 AU, and as a result, ~27 day periodicity in proton fluxes is significantly diminished at 1.43-1.38 AU. These events also cause unexpected impact on the Martian topside ionosphere, such as topside ionospheric depletion and compression observed by LPW and NGIMS onboard MAVEN. These observations are unique not only because of the recurring nature of electron enhancements seen at two vantage points, but also because they reveal unexpected impact of the weak CME and interplanetary shock on the Martian ionosphere, which provide new insight into the impact of CME-HSS interactions on Martian plasma environment.

categories: astro-ph.SR astro-ph.EP physics.geo-ph physics.space-ph


title: Statistical complexity is maximized close to criticality in cortical dynamics

abstract: Complex systems are typically characterized as an intermediate situation between a complete regular structure and a random system. Brain signals can be studied as a striking example of such systems: cortical states can range from highly synchronous and ordered neuronal activity (with higher spiking variability) to desynchronized and disordered regimes (with lower spiking variability). It has been recently shown, by testing independent signatures of criticality, that a phase transition occurs in a cortical state of intermediate spiking variability. Here, we use a symbolic information approach to show that, despite the monotonical increase of the Shannon entropy between ordered and disordered regimes, we can determine an intermediate state of maximum complexity based on the Jensen disequilibrium measure. More specifically, we show that statistical complexity is maximized close to criticality for cortical spiking data of urethane-anesthetized rats, as well as for a network model of excitable elements that presents a critical point of a non-equilibrium phase transition.


title: Extending the Hint Factory for the assistance dilemma: A novel, data-driven HelpNeed Predictor for proactive problem-solving help

abstract: Determining when and whether to provide personalized support is a well-known challenge called the assistance dilemma. A core problem in solving the assistance dilemma is the need to discover when students are unproductive so that the tutor can intervene. Such a task is particularly challenging for open-ended domains, even those that are well-structured with defined principles and goals. In this paper, we present a set of data-driven methods to classify, predict, and prevent unproductive problem-solving steps in the well-structured open-ended domain of logic. This approach leverages and extends the Hint Factory, a set of methods that leverages prior student solution attempts to build data-driven intelligent tutors. We present a HelpNeed classification, that uses prior student data to determine when students are likely to be unproductive and need help learning optimal problem-solving strategies. We present a controlled study to determine the impact of an Adaptive pedagogical policy that provides proactive hints at the start of each step based on the outcomes of our HelpNeed predictor: productive vs. unproductive. Our results show that the students in the Adaptive condition exhibited better training behaviors, with lower help avoidance, and higher help appropriateness (a higher chance of receiving help when it was likely to be needed), as measured using the HelpNeed classifier, when compared to the Control. Furthermore, the results show that the students who received Adaptive hints based on HelpNeed predictions during training significantly outperform their Control peers on the posttest, with the former producing shorter, more optimal solutions in less time. We conclude with suggestions on how these HelpNeed methods could be applied in other well-structured open-ended domains.


title: Towards Topic-Guided Conversational Recommender System

abstract: Conversational recommender systems (CRS) aim to recommend high-quality items to users through interactive conversations. To develop an effective CRS, the support of high-quality datasets is essential. Existing CRS datasets mainly focus on immediate requests from users, while lack proactive guidance to the recommendation scenario. In this paper, we contribute a new CRS dataset named \textbf{TG-ReDial} (\textbf{Re}commendation through \textbf{T}opic-\textbf{G}uided \textbf{Dial}og). Our dataset has two major features. First, it incorporates topic threads to enforce natural semantic transitions towards the recommendation scenario. Second, it is created in a semi-automatic way, hence human annotation is more reasonable and controllable. Based on TG-ReDial, we present the task of topic-guided conversational recommendation, and propose an effective approach to this task. Extensive experiments have demonstrated the effectiveness of our approach on three sub-tasks, namely topic prediction, item recommendation and response generation. TG-ReDial is available at https://github.com/RUCAIBox/TG-ReDial.


title: Testing Differential Privacy with Dual Interpreters

abstract: Applying differential privacy at scale requires convenient ways to check that programs computing with sensitive data appropriately preserve privacy. We propose here a fully automated framework for {\em testing} differential privacy, adapting a well-known "pointwise" technique from informal proofs of differential privacy. Our framework, called DPCheck, requires no programmer annotations, handles all previously verified or tested algorithms, and is the first fully automated framework to distinguish correct and buggy implementations of PrivTree, a probabilistically terminating algorithm that has not previously been mechanically checked. We analyze the probability of DPCheck mistakenly accepting a non-private program and prove that, theoretically, the probability of false acceptance can be made exponentially small by suitable choice of test size. We demonstrate DPCheck's utility empirically by implementing all benchmark algorithms from prior work on mechanical verification of differential privacy, plus several others and their incorrect variants, and show DPCheck accepts the correct implementations and rejects the incorrect variants. We also demonstrate how DPCheck can be deployed in a practical workflow to test differentially privacy for the 2020 US Census Disclosure Avoidance System (DAS).


title: Rainbow Solutions to the Sidon Equation in Cyclic Groups

abstract: Given a coloring of group elements, a rainbow solution to an equation is a solution whose every element is assigned a different color. The rainbow number of $\mathbb{Z}_n$ for an equation $eq$, denoted $rb(\mathbb{Z}_n,eq)$, is the smallest number of colors $r$ such that every exact $r$-coloring of $\mathbb{Z}_n$ admits a rainbow solution to the equation $eq$. We prove that for every exact $4$-coloring of $\mathbb{Z}_p$, where $p\geq 3$ is prime, there exists a rainbow solution to the Sidon equation $x_1+x_2=x_3+x_4$. Furthermore, we determine the rainbow number of $\mathbb{Z}_n$ for the Sidon equation.


title: All exact Borel subalgebras and all directed bocses are normal

abstract: Recently, Brzezi\'nski, Koenig and K\"ulshammer have introduced the notion of normal exact Borel subalgebra of a quasihereditary algebra. They have shown that there exists a one-to-one correspondence between normal directed bocses and quasihereditary algebras with a normal and homological exact Borel subalgebra. In this short note, we prove that every exact Borel subalgebra is automatically normal. As a corollary, we conclude that every directed bocs has a group-like element. These results simplify Brzezi\'nski, Koenig and K\"ulshammer's bijection.


title: The English Patient: Evaluating Local Lockdowns Using Real-Time COVID-19 & Consumption Data

abstract: We find UK 'local lockdowns' of cities and small regions, focused on limiting how many people a household can interact with and in what settings, are effective in turning the tide on rising positive COVID-19 cases. Yet, by focusing on household mixing within the home, these local lockdowns have not inflicted the large declines in consumption observed in March 2020 when the first virus wave and first national lockdown occurred. Our study harnesses a new source of real-time, transaction-level consumption data that we show to be highly correlated with official statistics. The effectiveness of local lockdowns are evaluated applying a difference-in-difference approach which exploits nearby localities not subject to local lockdowns as comparison groups. Our findings indicate that policymakers may be able to contain virus outbreaks without killing local economies. However, the ultimate effectiveness of local lockdowns is expected to be highly dependent on co-ordination between regions and an effective system of testing.


title: A Representation of Hyponormal Absolutely Norm Attaining Operators

abstract: In this article, we characterize absolutely norm attaining normal operators in terms of the essential spectrum. Later we prove a structure theorem for hyponormal absolutely norm attaining (or $\mathcal{AN}$-operators in short) and deduce conditions for the normality of the operator.


title: Identifying when Precession can be Measured in Gravitational Waveforms

abstract: In binary-black-hole systems where the black-hole spins are misaligned with the orbital angular momentum, precession effects leave characteristic modulations in the emitted gravitational waveform. Here, we investigate where in the parameter space we will be able to accurately identify precession, for likely observations over coming LIGO-Virgo-KAGRA observing runs. Despite the large number of parameters that characterise a precessing binary, we perform a large scale systematic study to identify the impact of each source parameter on the measurement of precession. We simulate a fiducial binary at moderate mass-ratio, signal-to-noise ratio (SNR), and spins, such that precession will be clearly identifiable, then successively vary each parameter while holding the remaining parameters fixed. As expected, evidence for precession increases with signal-to noise-ratio (SNR), higher in-plane spins, more unequal component masses, and higher inclination, but our study provides a quantitative illustration of each of these effects, and informs our intuition on parameter dependencies that have not yet been studied in detail, for example, the effect of varying the relative strength of the two polarisations, the total mass, and the aligned-spin components. We also measure the "precession SNR" $\rho_p$, which was introduced in Refs[1,2] to quantify the signal power associated with precession. By comparing $\rho_p$ with both Bayes factors and the recovered posterior distributions, we find it is a reliable metric for measurability that accurately predicts when the detected signal contains evidence for precession.


title: Convergence of the Finite Volume Method on Unstructured Meshes for a 3D Phase Field Model of Solidification

abstract: We present a convergence result for the finite volume method applied to a particular phase field problem suitable for simulation of pure substance solidification. The model consists of the heat equation and the phase field equation with a general form of the reaction term which encompasses a variety of existing models governing dendrite growth and elementary interface tracking problems. We apply the well known compact embedding techniques in the context of the finite volume method on admissible unstructured polyhedral meshes. We develop the necessary interpolation theory and derive an a priori estimate to obtain boundedness of the key terms. Based on this estimate, we conclude the convergence of all of the terms in the equation system.


title: A User-Friendly Computational Framework for Robust Structured Regression with the L$_2$ Criterion

abstract: We introduce a user-friendly computational framework for implementing robust versions of a wide variety of structured regression methods with the L$_{2}$ criterion. In addition to introducing an algorithm for performing L$_{2}$E regression, our framework enables robust regression with the L$_{2}$ criterion for additional structural constraints, works without requiring complex tuning procedures on the precision parameter, can be used to identify heterogeneous subpopulations, and can incorporate readily available non-robust structured regression solvers. We provide convergence guarantees for the framework and demonstrate its flexibility with some examples. Supplementary materials for this article are available online.


title: Topological shadows and complexity of islands in multiboundary wormholes

abstract: Recently, remarkable progress in recovering the Page curve of an evaporating black hole (BH) in Jackiw-Teitelboim gravity has been achieved through use of Quantum Extremal surfaces (QES). Multi-boundary Wormhole (MbW) models have been crucial in parallel model building in three dimensions. Motivated by this we here use the latter models to compute the subregion complexity of the Hawking quanta of the evaporating BH in AdS$_{3}$ and obtain the Page curve associated with this information theoretic measure. We use three- and $n$-boundary wormhole constructions to elucidate our computations of volumes below the Hubeny-Rangamani-Takayanagi (HRT) surfaces at different times. Time is represented by the growing length of the throat horizons corresponding to smaller exits of the multi-boundary wormhole and the evaporating bigger exit shrinks with evolving time. We track the change in choice of HRT surfaces with time and plot the volume with time. The smooth transition of Page curve is realized by a discontinuous jump at Page time in volume subregion complexity plots and the usual Page transition is realized as a phase transition due to the inclusion of the island in this context. We discuss mathematical intricacies and physical insights regarding the inclusion of the extra volume at Page time. The analysis is backed by calculations and lessons from kinematic space and tensor networks.


title: Isometric and affine copies of a set in volumetric Helly results

abstract: We show that for any compact convex set $K$ in $\mathbb{R}^d$ and any finite family $\mathcal{F}$ of convex sets in $\mathbb{R}^d$, if the intersection of every sufficiently small subfamily of $\mathcal{F}$ contains an isometric copy of $K$ of volume $1$, then the intersection of the whole family contains an isometric copy of $K$ scaled by a factor of $(1-\varepsilon)$, where $\varepsilon$ is positive and fixed in advance. Unless $K$ is very similar to a disk, the shrinking factor is unavoidable. We prove similar results for affine copies of $K$. We show how our results imply the existence of randomized algorithms that approximate the largest copy of $K$ that fits inside a given polytope $P$ whose expected runtime is linear on the number of facets of $P$.


title: Active volatile drops on liquid baths

abstract: In most experimental studies, active drops propel in a liquid bulk due to self-generated interfacial stresses of solutal origin. Here, we demonstrate the self-propulsion of a volatile drop on the surface of a liquid bath due to stresses of thermal origin. Evaporative heat pumping is converted into directed motion driven by thermocapillary stresses, which emerge on the drop surface as a result of a symmetry breaking of the drop temperature field. The dependence of the drop speed on the activity source, i.e. the evaporation flux, is derived with scaling arguments and captures the experimental data.


title: Analyzing the Weyl construction for dynamical Cartan subalgebras

abstract: When the reduced twisted $C^*$-algebra $C^*_r(\mathcal{G}, c)$ of a non-principal groupoid $\mathcal{G}$ admits a Cartan subalgebra, Renault's work on Cartan subalgebras implies the existence of another groupoid description of $C^*_r(\mathcal{G}, c)$. In an earlier paper, joint with Reznikoff and Wright, we identified situations where such a Cartan subalgebra arises from a subgroupoid $\mathcal{S}$ of $\mathcal{G}$. In this paper, we study the relationship between the original groupoids $\mathcal{S}, \mathcal{G}$ and the Weyl groupoid and twist associated to the Cartan pair. We first identify the spectrum $\mathfrak{B}$ of the Cartan subalgebra $C^*_r(\mathcal{S}, c)$. We then show that the quotient groupoid $\mathcal{G}/\mathcal{S}$ acts on $\mathfrak{B}$, and that the corresponding action groupoid is exactly the Weyl groupoid of the Cartan pair. Lastly we show that, if the quotient map $\mathcal{G}\to\mathcal{G}/\mathcal{S}$ admits a continuous section, then the Weyl twist is also given by an explicit continuous $2$-cocycle on $\mathcal{G}/\mathcal{S} \ltimes \mathfrak{B}$.


title: The universal function of the diffractive process in color dipole picture

abstract: In this presentation, we obtain the corresponding universal function to the diffractive process and show the cross section exhibits the geometrical scaling. It is observed the diffractive theory according to the color dipole approach at small-$x$ is a convenient framework that reveals the color transparency and the saturation phenomena. Also we calculate the contribution of heavy quark productions in the diffractive cross section for high energy that is determined by the small size dipole configuration. The ratio of the diffractive cross section to the total cross section in the electron-proton collision is the other important quantity that is computed in this work.


title: All quasihereditary algebras with a regular exact Borel subalgebra

abstract: Not every quasihereditary algebra $(A,\Phi,\unlhd)$ has an exact Borel subalgebra. A theorem by Koenig, K\"ulshammer and Ovsienko asserts that there always exists a quasihereditary algebra Morita equivalent to $A$ that has a regular exact Borel subalgebra, but a characterisation of such a Morita representative is not directly obtainable from their work. This paper gives a criterion to decide whether a quasihereditary algebra contains a regular exact Borel subalgebra and provides a method to compute all the representatives of $A$ that have a regular exact Borel subalgebra. It is shown that the Cartan matrix of a regular exact Borel subalgebra of a quasihereditary algebra $(A,\Phi,\unlhd)$ only depends on the composition factors of the standard and costandard $A$-modules and on the dimension of the $\operatorname{Hom}$-spaces between standard $A$-modules. We also characterise the basic quasihereditary algebras that admit a regular exact Borel subalgebra.


title: Hierarchical PCA and Modeling Asset Correlations

abstract: Modeling cross-sectional correlations between thousands of stocks, across countries and industries, can be challenging. In this paper, we demonstrate the advantages of using Hierarchical Principal Component Analysis (HPCA) over the classic PCA. We also introduce a statistical clustering algorithm for identifying of homogeneous clusters of stocks, or "synthetic sectors". We apply these methods to study cross-sectional correlations in the US, Europe, China, and Emerging Markets.


title: DART: A Lightweight Quality-Suggestive Data-to-Text Annotation Tool

abstract: We present a lightweight annotation tool, the Data AnnotatoR Tool (DART), for the general task of labeling structured data with textual descriptions. The tool is implemented as an interactive application that reduces human efforts in annotating large quantities of structured data, e.g. in the format of a table or tree structure. By using a backend sequence-to-sequence model, our system iteratively analyzes the annotated labels in order to better sample unlabeled data. In a simulation experiment performed on annotating large quantities of structured data, DART has been shown to reduce the total number of annotations needed with active learning and automatically suggesting relevant labels.


title: Chiral Hinge Magnons in Second-Order Topological Magnon Insulators

abstract: When interacting spins in condensed matter order ferromagnetically, their ground state wave function is topologically trivial. Nonetheless, in two dimensions, the ferromagnetic state can support spin excitations with nontrivial topology, an exotic state known as topological magnon insulator (TMI). Here, we theoretically unveil and numerically confirm a novel ferromagnetic state in three dimensions dubbed second-order TMI, whose hallmarks are excitations at its hinges, where facets intersect. Since ferromagnetism naturally comes with broken time-reversal symmetry, the hinge magnons are chiral, rendering backscattering impossible. Hence, they trace out a three-dimensional path about the sample unimpeded by defects and are topologically protected by the spectral gap. They are remarkably robust against disorder and simultaneously highly tunable by atomic-level engineering of the sample termination. Our findings empower magnonics with the tools of higher-order topology, a promising route to combine low-energy information transfer free of Joule heating with three-dimensional vertical integration.


title: Deep SVBRDF Estimation on Real Materials

abstract: Recent work has demonstrated that deep learning approaches can successfully be used to recover accurate estimates of the spatially-varying BRDF (SVBRDF) of a surface from as little as a single image. Closer inspection reveals, however, that most approaches in the literature are trained purely on synthetic data, which, while diverse and realistic, is often not representative of the richness of the real world. In this paper, we show that training such networks exclusively on synthetic data is insufficient to achieve adequate results when tested on real data. Our analysis leverages a new dataset of real materials obtained with a novel portable multi-light capture apparatus. Through an extensive series of experiments and with the use of a novel deep learning architecture, we explore two strategies for improving results on real data: finetuning, and a per-material optimization procedure. We show that adapting network weights to real data is of critical importance, resulting in an approach which significantly outperforms previous methods for SVBRDF estimation on real materials. Dataset and code are available at https://lvsn.github.io/real-svbrdf


title: Towards Blood Flow in the Virtual Human: Efficient Self-Coupling of HemeLB

abstract: Many scientific and medical researchers are working towards the creation of a virtual human - a personalised digital copy of an individual - that will assist in a patient's diagnosis, treatment and recovery. The complex nature of living systems means that the development of this remains a major challenge. We describe progress in enabling the HemeLB lattice Boltzmann code to simulate 3D macroscopic blood flow on a full human scale. Significant developments in memory management and load balancing allow near linear scaling performance of the code on hundreds of thousands of computer cores. Integral to the construction of a virtual human, we also outline the implementation of a self-coupling strategy for HemeLB. This allows simultaneous simulation of arterial and venous vascular trees based on human-specific geometries.


title: A Survey On Anti-Spoofing Methods For Face Recognition with RGB Cameras of Generic Consumer Devices

abstract: The widespread deployment of face recognition-based biometric systems has made face Presentation Attack Detection (face anti-spoofing) an increasingly critical issue. This survey thoroughly investigates the face Presentation Attack Detection (PAD) methods, that only require RGB cameras of generic consumer devices, over the past two decades. We present an attack scenario-oriented typology of the existing face PAD methods and we provide a review of over 50 of the most recent face PAD methods and their related issues. We adopt a comprehensive presentation of the methods that have most influenced face PAD following the proposed typology, and in chronological order. By doing so, we depict the main challenges, evolutions and current trends in the field of face PAD, and provide insights on its future research. From an experimental point of view, this survey paper provides a summarized overview of the available public databases and extensive comparative experimental results of different PAD methods.


title: Optimal control of vaccination and plasma transfusion with potential usefulness for COVID-19

abstract: The SEIR model is a compartmental model used to simulate the dynamics of an epidemic. In this chapter, we introduce two control functions in the compartmental SEIR model representing vaccination and plasma transfusion. Optimal control problems are proposed to study the effects of these two control measures, on the reduction of infected individuals and increase of recovered ones, with minimal costs. Up to our knowledge, the plasma transfusion treatment has never been considered as a control strategy for epidemics mitigation. The proposed vaccination and treatment strategies may have a real application in the challenging and hard problem of controlling the COVID-19 pandemic.


title: Automatic generation of reviews of scientific papers

abstract: With an ever-increasing number of scientific papers published each year, it becomes more difficult for researchers to explore a field that they are not closely familiar with already. This greatly inhibits the potential for cross-disciplinary research. A traditional introduction into an area may come in the form of a review paper. However, not all areas and sub-areas have a current review. In this paper, we present a method for the automatic generation of a review paper corresponding to a user-defined query. This method consists of two main parts. The first part identifies key papers in the area by their bibliometric parameters, such as a graph of co-citations. The second stage uses a BERT based architecture that we train on existing reviews for extractive summarization of these key papers. We describe the general pipeline of our method and some implementation details and present both automatic and expert evaluations on the PubMed dataset.


title: A novel derivation of rigorous macroscopic limits from a micro-meso description of signal-triggered cell migration in fibrous environments

abstract: In this work we upscale a prototypical kinetic transport equation which models a cell population moving in a fibrous environment with a chemo- or haptotactic signal influencing both the direction and the magnitude of the cell velocity. The presented approach to scaling does not rely on orthogonality and treats parabolic and hyperbolic scalings in a unified manner. It is shown that the steps of the formal limit procedures are mirrored by rigorous operations with finite measures provided that the measure-valued position-direction fiber distribution enjoys some spacial continuity.


title: Steenrod operations on the de Rham cohomology of algebraic stacks

abstract: Building up on work of Epstein, May and Drury, we define and investigate the mod $p$ Steenrod operations on the de Rham cohomology of smooth algebraic stacks over a field of characteristic $p>0$. We then compute the action of the operations on the de Rham cohomology of classifying stacks for finite groups, connected reductive groups for which $p$ is not a torsion prime, and (special) orthogonal groups when $p=2$.


title: Memory in 3D cyclically driven granular material

abstract: We perform experimental and numerical studies of a granular system under cyclic-compression to investigate reversibility and memory effects. We focus on the quasi-static forcing of dense systems, which is most relevant to a wide range of geophysical, industrial, and astrophysical problems. We find that soft-sphere simulations with proper stiffness and friction quantitatively reproduce both the translational and rotational displacements of the grains. We then utilize these simulations to demonstrate that such systems are capable of storing the history of previous compressions. While both mean translational and rotational displacements encode such memory, the response is fundamentally different for translations compared to rotations. For translational displacements, this memory of prior forcing depends on the coefficient of static inter-particle friction, but rotational memory is not altered by the level of friction.


title: On the detectability of ultra-compact binary pulsar systems

abstract: Using neural networks, we integrate the ability to account for Doppler smearing due to a pulsar's orbital motion with the pulsar population synthesis package \psrpoppy\ to develop accurate modeling of the observed binary pulsar population. As a first application, we show that binary neutron star systems where the two components have highly unequal mass are, on average, easier to detect than systems which are symmetric in mass. We then investigate the population of ultra-compact ($1.5 \, {\rm min} \leq P_{\rm b} \leq 15\,\rm min$) neutron star--white dwarf (NS--WD) and double neutron star (DNS) systems which are promising sources for the Laser Interferometer Space Antenna gravitational-wave detector. Given the non-detection of these systems in radio surveys thus far, we estimate a 95\% confidence upper limit of $\sim$1450 and $\sim$1100 ultra-compact NS--WD and DNS systems in the Milky Way that are beaming towards the Earth respectively. We also show that using survey integration times in the range 20~s to 200~s with time-domain resampling will maximize the signal-to-noise ratio as well as the probability of detection of these ultra-compact binary systems. Among all the large scale radio pulsar surveys, those that are currently being carried out at the Arecibo radio telescope have $\sim$50--80\% chance of detecting at least one of these systems using current integration integration times and $\sim$80--95\% using optimal integration times in the next several years.


title: Space Charge Aberrations in Single Shot Time-Resolved Transmission Electron Microscopy

abstract: In this paper, we discuss effects of space charge fields on imaging performance in a single shot time-resolved electron microscope (TEM). Using a Green's function perturbation method, we derive analytical estimates for the effects of space charge nonlinearity on the image formation process and the associated aberration coefficients. The dependence of these coefficients on the initial beam phase space distribution is elucidated. The results are validated by particle tracking simulations and provide fundamental scaling laws for the trade-off between temporal and spatial resolution in single-shot time-resolved TEM.


title: Olympus: a benchmarking framework for noisy optimization and experiment planning

abstract: Research challenges encountered across science, engineering, and economics can frequently be formulated as optimization tasks. In chemistry and materials science, recent growth in laboratory digitization and automation has sparked interest in optimization-guided autonomous discovery and closed-loop experimentation. Experiment planning strategies based on off-the-shelf optimization algorithms can be employed in fully autonomous research platforms to achieve desired experimentation goals with the minimum number of trials. However, the experiment planning strategy that is most suitable to a scientific discovery task is a priori unknown while rigorous comparisons of different strategies are highly time and resource demanding. As optimization algorithms are typically benchmarked on low-dimensional synthetic functions, it is unclear how their performance would translate to noisy, higher-dimensional experimental tasks encountered in chemistry and materials science. We introduce Olympus, a software package that provides a consistent and easy-to-use framework for benchmarking optimization algorithms against realistic experiments emulated via probabilistic deep-learning models. Olympus includes a collection of experimentally derived benchmark sets from chemistry and materials science and a suite of experiment planning strategies that can be easily accessed via a user-friendly python interface. Furthermore, Olympus facilitates the integration, testing, and sharing of custom algorithms and user-defined datasets. In brief, Olympus mitigates the barriers associated with benchmarking optimization algorithms on realistic experimental scenarios, promoting data sharing and the creation of a standard framework for evaluating the performance of experiment planning strategies


title: Quasielastic Electromagnetic Scattering Cross Sections and World Data Comparisons in the {\fontfamily{qcr}\selectfont GENIE} Monte Carlo Event Generator

abstract: The usage of Monte Carlo neutrino event generators (MC$\nu$EGs) is a norm within the high-energy $\nu$ scattering community. The relevance of quasielastic (QE) energy regimes to $\nu$ oscillation experiments implies that accurate calculations of $\nu A$ cross sections in this regime will be a key contributor to reducing the systematic uncertainties affecting the extraction of oscillation parameters. In spite of this, many MC$\nu$EGs utilize highly phenomenological, parameterized models of QE scattering cross sections. Moreover, a culture of validation of MC$\nu$EGs against prolific electron ($e$) scattering data has been historically lacking. In this work, we implement new $e A$ cross sections obtained from nuclear ab initio approaches in GENIE, the primary MC$\nu$EG utilized by the FNAL community. In particular, we utilize results from Quantum MC methods which solve the many-body nuclear problem in the Short-Time Approximation (STA), allowing consistent retention of two-nucleon dynamics which are crucial to explain available nuclear electromagnetic (electroweak) data over a wide range of energy and momentum transfers. This new implementation in GENIE is fully tested against the world QE electromagnetic data, finding agreement with available data below $\sim2\,$GeV of beam energy with the aid of a scaling function formalism. The STA is currently limited to study $A\leq12$ nuclei, however, its semi-inclusive multibody identity components are exportable to other many-body computational techniques such as Auxiliary Field Diffusion MC which can reach $A\leq40$ systems while continuing to realize the factorization contained within the STA's multinucleon dynamics. Together, these developments promise to make future experiments such as DUNE more accurate in their assessment of MC$\nu$EG systematics, $\nu$ properties, and potentially empower the discovery of physics beyond the Standard Model.


title: Strongly lensed supernovae as a self-sufficient probe of the distance duality relation

abstract: The observation of strongly lensed Type Ia supernovae enables both the luminosity and angular diameter distance to a source to be measured simultaneously using a single observation. This feature can be used to measure the distance duality parameter $\eta(z)$ without relying on multiple datasets and cosmological assumptions to reconstruct the relation between angular and luminosity distances. In this paper, we show how this can be achieved by future observations of strongly lensed Type Ia systems. Using simulated datasets, we reconstruct the function $\eta(z)$ using both parametric and non-parametric approaches, focusing on Genetic Algorithms and Gaussian processes for the latter. In the parametric approach, we find that in the realistic scenario of $N_{\rm lens}=20$ observed systems, the parameter $\epsilon_0$ used to describe the trend of $\eta(z)$ can be constrained with the precision achieved by current SNIa and BAO surveys, while in the futuristic case ($N_{\rm lens}=1000$) these observations could be competitive with the forecast precision of upcoming LSS and SN surveys. Using the machine learning approaches of Genetic Algorithms and Gaussian processes, we find that both reconstruction methods are generally well able to correctly recover the underlying fiducial model in the mock data, even in the realistic case of $N_{\rm lens}=20$. Both approaches learn effectively from the features of the mock data points, yielding $1\sigma$ constraints that are in excellent agreement with the parameterised results.


title: Automating Inference of Binary Microlensing Events with Neural Density Estimation

abstract: Automated inference of binary microlensing events with traditional sampling-based algorithms such as MCMC has been hampered by the slowness of the physical forward model and the pathological likelihood surface. Current analysis of such events requires both expert knowledge and large-scale grid searches to locate the approximate solution as a prerequisite to MCMC posterior sampling. As the next generation, space-based microlensing survey with the Roman Space Observatory is expected to yield thousands of binary microlensing events, a new scalable and automated approach is desired. Here, we present an automated inference method based on neural density estimation (NDE). We show that the NDE trained on simulated Roman data not only produces fast, accurate, and precise posteriors but also captures expected posterior degeneracies. A hybrid NDE-MCMC framework can further be applied to produce the exact posterior.

categories: astro-ph.IM astro-ph.EP astro-ph.SR cs.LG physics.data-an


title: Online and Distribution-Free Robustness: Regression and Contextual Bandits with Huber Contamination

abstract: In this work we revisit two classic high-dimensional online learning problems, namely linear regression and contextual bandits, from the perspective of adversarial robustness. Existing works in algorithmic robust statistics make strong distributional assumptions that ensure that the input data is evenly spread out or comes from a nice generative model. Is it possible to achieve strong robustness guarantees even without distributional assumptions altogether, where the sequence of tasks we are asked to solve is adaptively and adversarially chosen? We answer this question in the affirmative for both linear regression and contextual bandits. In fact our algorithms succeed where conventional methods fail. In particular we show strong lower bounds against Huber regression and more generally any convex M-estimator. Our approach is based on a novel alternating minimization scheme that interleaves ordinary least-squares with a simple convex program that finds the optimal reweighting of the distribution under a spectral constraint. Our results obtain essentially optimal dependence on the contamination level $\eta$, reach the optimal breakdown point, and naturally apply to infinite dimensional settings where the feature vectors are represented implicitly via a kernel map.


title: Why reducing the cosmic sound horizon alone can not fully resolve the Hubble tension

abstract: The mismatch between the locally measured expansion rate of the universe and the one inferred from the cosmic microwave background measurements by Planck in the context of the standard $\Lambda$CDM, known as the Hubble tension, has become one of the most pressing problems in cosmology. A large number of amendments to the $\Lambda$CDM model have been proposed in order to solve this tension. Many of them introduce new physics, such as early dark energy, modifications of the standard model neutrino sector, extra radiation, primordial magnetic fields or varying fundamental constants, with the aim of reducing the sound horizon at recombination $r_{\star}$. We demonstrate here that any model which only reduces $r_{\star}$ can never fully resolve the Hubble tension while remaining consistent with other cosmological datasets. We show explicitly that models which achieve a higher Hubble constant with lower values of matter density $\Omega_m h^2$ run into tension with the observations of baryon acoustic oscillations, while models with larger $\Omega_mh^2$ develop tension with galaxy weak lensing data.


title: Deformable DETR: Deformable Transformers for End-to-End Object Detection

abstract: DETR has been recently proposed to eliminate the need for many hand-designed components in object detection while demonstrating good performance. However, it suffers from slow convergence and limited feature spatial resolution, due to the limitation of Transformer attention modules in processing image feature maps. To mitigate these issues, we proposed Deformable DETR, whose attention modules only attend to a small set of key sampling points around a reference. Deformable DETR can achieve better performance than DETR (especially on small objects) with 10 times less training epochs. Extensive experiments on the COCO benchmark demonstrate the effectiveness of our approach. Code is released at https://github.com/fundamentalvision/Deformable-DETR.


title: Non-perturbative gauge transformations of arbitrary fermion correlation functions in quantum electrodynamics

abstract: We study the transformation of the dressed electron propagator and the general $N$-point functions under a change in the covariant gauge of internal photon propagators. We re-establish the well known Landau-Khalatnikov-Fradkin transformation for the propagator and generalise it to arbitrary correlation functions in configuration space, finding that it coincides with the analogous result for scalar fields. We comment on the consequences for perturbative application in momentum-space.


title: Oscillatory magnetic fields for neutron resonance spin-echo spectroscopy

abstract: The generation of high frequency oscillatory magnetic fields represents a fundamental component underlying the successful implementation of neutron resonant spin-echo spectrometers, a class of instrumentation critical for the high-resolution extraction of dynamical excitations (structural and magnetic) in materials. In this paper, the setup of the resonant circuits at the longitudinal resonant spin-echo spectrometer RESEDA is described in comprehensive technical detail. We demonstrate that these circuits are capable of functioning at frequencies up to 3.6 MHz and over a broad bandwidth down to 35 kHz using a combination of signal generators, amplifiers, impedance matching transformers, and a carefully designed cascade of tunable capacitors and customized coils.


title: How to search for multiple messengers -- a general framework beyond two messengers

abstract: Quantification of the significance of a candidate multi-messenger detection of cosmic events is an emerging need in the astrophysics and astronomy communities. In this paper we show that a model-independent optimal search does not exist, and we present a general Bayesian method for the optimal model-dependent search, which is scalable to any number and any kind of messengers, and applicable to any model. In the end, we demonstrate it through an example for a joint gravitational wave, high-energy neutrino, short gamma-ray burst event search; which has not been examined heretofore.


title: Enhanced Lidov-Kozai migration and the formation of the transiting giant planet WD1856+534b

abstract: We investigate the possible origin of the transiting giant planet WD1856+534b, the first strong exoplanet candidate orbiting a white dwarf, through high-eccentricity migration (HEM) driven by the Lidov-Kozai (LK) effect. The host system's overall architecture is an hierarchical quadruple in the '2+2' configuration, owing to the presence of a tertiary companion system of two M-dwarfs. We show that a secular inclination resonance in 2+2 systems can significantly broaden the LK window for extreme eccentricity excitation ($e \gtrsim 0.999$), allowing the giant planet to migrate for a wide range of initial orbital inclinations. Octupole effects can also contribute to the broadening of this 'extreme' LK window. By requiring that perturbations from the companion stars be able to overcome short-range forces and excite the planet's eccentricity to $e \simeq 1$, we obtain an absolute limit of $a_{1} \gtrsim 8 \, {\rm AU} \, (a_{3} / 1500 \, {\rm AU})^{6/7}$ for the planet's semi-major axis just before migration (where $a_{3}$ is the semi-major axis of the 'outer' orbit). We suggest that, to achieve a wide LK window through the 2+2 resonance, WD1856b likely migrated from $30 \, {\rm AU} \lesssim a_{1} \lesssim 60 \, {\rm AU}$, corresponding to $\sim 10$--$20 \, {\rm AU}$ during the host's main-sequence phase. We discuss possible difficulties of all flavours of HEM affecting the occurrence rate of short-period giant planets around white dwarfs.


title: Mitigating the effects of undersampling in weak lensing shear estimation with metacalibration

abstract: Metacalibration is a state-of-the-art technique for measuring weak gravitational lensing shear from well-sampled galaxy images. We investigate the accuracy of shear measured with metacalibration from fitting elliptical Gaussians to undersampled galaxy images. In this case, metacalibration introduces aliasing effects leading to an ensemble multiplicative shear bias about 0.01 for Euclid, and even larger for the Roman Space Telescope, well exceeding the missions' requirements. We find that this aliasing bias can be mitigated by computing shapes from weighted moments with wider Gaussians as weight functions, thereby trading bias for a slight increase in variance of the measurements. We show that this approach is robust to the point-spread function in consideration and meets the stringent requirements of Euclid for galaxies with moderate to high signal-to-noise ratios. We therefore advocate metacalibration as a viable shear measurement option for weak lensing from upcoming space missions.


title: Did Sgr cause the vertical waves in the solar neighbourhood?

abstract: The vertical distribution of stars in the solar neighbourhood is not in equilibrium but contains a wave signature in both density and velocity space originating from a perturbation. With the discovery of the phase-space spiral in Gaia data release 2, determining the origin of this perturbation has become even more urgent. We develop and test a fast method for calculating the perturbation from a passing satellite on the vertical component of a part of a disc galaxy. This fast method allows us to test a large variety of possible perturbations to the vertical disc very quickly. We apply our method to the range of possible perturbations to the solar neighbourhood stemming from the recent passage of the Sagittarius dwarf galaxy (Sgr), varying its mass, mass profile, and present-day position within their observational uncertainties, and its orbit within different realistic models for the Milky Way's gravitational potential. We find that we are unable to reproduce the observed asymmetry in the vertical number counts and its concomitant breathing mode in velocity space for any plausible combination of Sgr and Milky-Way properties. In all cases, either the amplitude or the perturbation wavelength of the number-count asymmetry and of the oscillations in the mean vertical velocity produced by the passage of Sgr is in large disagreement with the observations from Gaia DR2. We conclude that Sgr cannot have caused the observed oscillations in the vertical disc or the Gaia phase-space spiral.


title: Coordinated Assembly of Galaxy Groups and Clusters in the IllustrisTNG Simulations

abstract: Recent stellar population analysis of early-type galaxy spectra has demonstrated that the low-mass galaxies in cluster centers have high [$\alpha/\rm Fe$] and old ages characteristic of massive galaxies and unlike the low-mass galaxy population in the outskirts of clusters and fields. This phenomenon has been termed "coordinated assembly" to highlight the fact that the building blocks of massive cluster central galaxies are drawn from a special subset of the overall low-mass galaxy population. Here we explore this idea in the IllustrisTNG simulations, particularly the TNG300 run, in order to understand how environment, especially cluster centers, shape the star formation histories of quiescent satellite galaxies in groups and clusters ($M_{200c,z=0}\geq10^{13} M_{\odot}$). Tracing histories of quenched satellite galaxies with $M_{\star,z=0}\geq10^{10} M_{\odot}$, we find that those in more massive dark matter halos, and located closer to the primary galaxies, are quenched earlier, have shorter star formation timescales, and older stellar ages. The star formation timescale-$M_{\star}$ and stellar age-$M_{\star}$ scaling relations are in good agreement with observations, and are predicted to vary with halo mass and cluster-centric distance. The dependence on environment arises due to the infall histories of satellite galaxies: galaxies that are located closer to cluster centers in more massive dark matter halos at $z=0$ were accreted earlier on average. The delay between infall and quenching time is shorter for galaxies in more massive halos, and depends on the halo mass at its first accretion, showing that group pre-processing is a crucial aspect in satellite quenching.


title: Statistical Signatures of Nanoflare Activity. II. A Nanoflare Explanation for Periodic Brightenings in Flare Stars observed by NGTS

abstract: Several studies have documented periodic and quasi-periodic signals from the time series of dMe flare stars and other stellar sources. Such periodic signals, observed within quiescent phases (i.e., devoid of larger-scale microflare or flare activity), range in period from $1-1000$ seconds and hence have been tentatively linked to ubiquitous $p$-mode oscillations generated in the convective layers of the star. As such, most interpretations for the observed periodicities have been framed in terms of magneto-hydrodynamic wave behavior. However, we propose that a series of continuous nanoflares, based upon a power-law distribution, can provide a similar periodic signal in the associated time series. Adapting previous statistical analyses of solar nanoflare signals, we find the first statistical evidence for stellar nanoflare signals embedded within the noise envelope of M-type stellar lightcurves. Employing data collected by the Next Generation Transit Survey (NGTS), we find evidence for stellar nanoflare activity demonstrating a flaring power-law index of $3.25 \pm 0.20 $, alongside a decay timescale of $200 \pm 100$ s. We also find that synthetic time series, consistent with the observations of dMe flare star lightcurves, are capable of producing quasi-periodic signals in the same frequency range as $p$-mode signals, despite being purely comprised of impulsive signatures. Phenomena traditionally considered a consequence of wave behaviour may be described by a number of high frequency but discrete nanoflare energy events. This new physical interpretation presents a novel diagnostic capability, by linking observed periodic signals to given nanoflare model conditions.


title: Limits and Security of Free-Space Quantum Communications

abstract: The study of free-space quantum communications requires tools from quantum information theory, optics and turbulence theory. Here we combine these tools to bound the ultimate rates for key and entanglement distribution through a free-space link, where the propagation of quantum systems is generally affected by diffraction, atmospheric extinction, turbulence, pointing errors, and background noise. Besides establishing ultimate limits, we also show that the composable secret-key rate achievable by a suitable (pilot-guided and post-selected) coherent-state protocol is sufficiently close to these limits, therefore showing the suitability of free-space channels for high-rate quantum key distribution. Our work provides analytical tools for assessing the composable finite-size security of coherent-state protocols in general conditions, from the standard assumption of a stable communication channel (as is typical in fiber-based connections) to the more challenging scenario of a fading channel (as is typical in free-space links).

categories: quant-ph cond-mat.other physics.ao-ph physics.app-ph physics.optics


title: The critical dark matter halo mass for Population III star formation: dependence on Lyman-Werner radiation, baryon-dark matter streaming velocity, and redshift

abstract: A critical dark matter halo mass ($M_{\rm crit}$) for Population III (Pop III) stars can be defined as the typical minimum halo mass that hosts sufficient cold dense gas required for the formation of the first stars. The presence of Lyman-Werner (UV) radiation, which can dissociate molecular hydrogen, and the baryon-dark matter streaming velocity both delay the formation of Pop III stars by increasing $M_{\rm crit}$. In this work, we constrain $M_{\rm crit}$ as a function of Lyman-Werner flux (including self-shielding), baryon-dark matter streaming, and redshift using cosmological simulations with a large sample of halos utilizing the adaptive mesh refinement (AMR) code ENZO. We provide a fit for $M_{\rm crit}$ as a function of these quantities which we expect to be particularly useful for semi-analytical models of early galaxy formation. In addition, we find: (i) the measured redshift dependence of $M_{\rm crit}$ in the absence of radiation or streaming is $(1+z)^{-1.58}$, consistent with a constant virial temperature; (ii) increasing the UV background increases $M_{\rm crit}$ while steepening the redshift dependence, up to $(1+z)^{-5.7}$; (iii) baryon-dark matter streaming boosts $M_{\rm crit}$ but flattens the dependence on redshift; (iv) the combination of the two effects is not simply multiplicative.


title: The effects of turbulence on galactic nuclear gas rings

abstract: The gas dynamics in the inner few kiloparsecs of barred galaxies often results in configurations that give rise to nuclear gas rings. However, the generic dynamical description of the formation of galactic nuclear rings does not take into account the effects of thermal pressure and turbulence. Here we perform 3D hydrodynamic simulations of gas in a galactic barred potential out to a radius of $2$~kpc. We include self-gravity and a large-scale turbulence-driving module. We study how the formation of gaseous nuclear rings is affected by changing the bar pattern speed and the strength of the turbulence, and compare the results to simulations with a relatively high sound speed. We use two values for the bar pattern speed that have been discussed in the literature for our Milky Way Galaxy (40 and 63 km s$^{-1}$ kpc$^{-1}$). Our results show that turbulence produces broader and smaller nuclear rings, and enhances gas migration towards the inner few pc of the galaxy, compared to simulations without turbulence.


title: Simultaneous extraction of $\alpha_s$ and $m_t$ from LHC $t\bar{t}$ differential distributions

abstract: We present a joint extraction of the strong coupling $\alpha_s$ and the top-quark pole mass $m_t$ from measurements of top-quark pair production performed by the ATLAS and CMS experiments at the 8 TeV LHC. For the first time, differential NNLO theory predictions for different values of the top-quark mass are utilised for four kinematic distributions: the average transverse momentum of the top-quark, its average rapidity and the pair invariant mass and rapidity. The use of fastNLO tables for these distributions allows rapid evaluation of the differential theory predictions for different PDF sets. We consider the single differential distributions from the experiments both separately and in combination in order to obtain the best fit to theory. Our final values are $\alpha_s=0.1159^{+0.0013}_{-0.0014}$ and $m_t=173.8^{+0.8}_{-0.8}$ GeV which are compatible with previous extractions using top-quark measurements. In the case of $m_t$, our value is also compatible with the world average value collated by the Particle Data Group.


title: Counting D1-D5-P Microstates in Supergravity

abstract: We quantize the D1-D5-P microstate geometries known as superstrata directly in supergravity. We use Rychkov's consistency condition [hep-th/0512053] which was derived for the D1-D5 system; for superstrata, this condition turns out to be strong enough to fix the symplectic form uniquely. For the $(1,0,n)$ superstrata, we further confirm this quantization by a bona-fide explicit computation of the symplectic form using the semi-classical covariant quantization method in supergravity. We use the resulting quantizations to count the known supergravity superstrata states, finding agreement with previous countings that the number of these states grows parametrically smaller than those of the corresponding black hole.


title: Automatic Post-selection by Ancillae Thermalisation

abstract: Tasks such as classification of data and determining the groundstate of a Hamiltonian cannot be carried out through purely unitary quantum evolution. Instead, the inherent non-unitarity of the measurement process must be harnessed. Post-selection and its extensions provide a way to do this. However they make inefficient use of time resources -- a typical computation might require $O(2^m)$ measurements over $m$ qubits to reach a desired accuracy. We propose a method inspired by the eigenstate thermalisation hypothesis, that harnesses the induced non-linearity of measurement on a subsystem. Post-selection on $m$ ancillae qubits is replaced with tracing out $O(\log\epsilon / \log(1-p))$ (where p is the probability of a successful measurement) to attain the same accuracy as the post-selection circuit. We demonstrate this scheme on the quantum perceptron and phase estimation algorithm. This method is particularly advantageous on current quantum computers involving superconducting circuits.


title: Super-sample covariance of the thermal Sunyaev-Zel'dovich effect

abstract: The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful probe of cosmology. The statistical errors in the tSZ power spectrum measurements are dominated by the presence of massive clusters in a survey volume that are easy to identify on individual cluster basis. First, we study the impact of super sample covariance (SSC) on the tSZ power spectrum measurements, and find that the sample variance is dominated by the connected non-Gaussian (cNG) covariance arising mainly from Poisson number fluctuations of massive clusters in the survey volume. Second, we find that removing such individually-detected, massive clusters from the analysis significantly reduces the cNG contribution, thereby leading the SSC to be a leading source of the sample variance. We then show, based on Fisher analysis, that the power spectrum measured from the remaining diffuse tSZ effects can be used to obtain tight constraints on cosmological parameters as well as the hydrostatic mass bias parameter. Our method offers complementary use of individual tSZ cluster counts and the power spectrum measurements of diffuse tSZ signals for cosmology and intracluster gas physics.


title: Device-independent quantum key distribution from computational assumptions

abstract: In device-independent quantum key distribution (DIQKD), an adversary prepares a device consisting of two components, distributed to Alice and Bob, who use the device to generate a secure key. The security of existing DIQKD schemes holds under the assumption that the two components of the device cannot communicate with one another during the protocol execution. This is called the no-communication assumption in DIQKD. Here, we show how to replace this assumption, which can be hard to enforce in practice, by a standard computational assumption from post-quantum cryptography: we give a protocol that produces secure keys even when the components of an adversarial device can exchange arbitrary quantum communication, assuming the device is computationally bounded. Importantly, the computational assumption only needs to hold during the protocol execution -- the keys generated at the end of the protocol are information-theoretically secure as in standard DIQKD protocols.


title: Elucidating Galaxy Assembly Bias in SDSS

abstract: We investigate the level of galaxy assembly bias in the Sloan Digital Sky Survey (SDSS) main galaxy sample using ELUCID, a state-of-the-art constrained simulation that accurately reconstructed the initial density perturbations within the SDSS volume. On top of the ELUCID haloes, we develop an extended HOD model that includes the assembly bias of central and satellite galaxies, parameterized as $\mathcal{Q}_\mathrm{cen}$ and $\mathcal{Q}_\mathrm{sat}$, respectively, to predict a suite of one- and two-point observables. In particular, our fiducial constraint employs the probability distribution of the galaxy number counts measured on $8\,\mathrm{Mpc}\,h^{-1}$ scales $N_8^g$ and the projected cross-correlation functions of quintiles of galaxies selected by $N_8^g$ with our entire galaxy sample. We perform extensive tests of the efficacy of our method by fitting the same observables to mock data using both constrained and non-constrained simulations. We discover that in many cases the level of cosmic variance between the two simulations can produce biased constraints that lead to an erroneous detection of galaxy assembly bias if the non-constrained simulation is used. When applying our method to the SDSS data, the ELUCID reconstruction effectively removes an otherwise strong degeneracy between cosmic variance and galaxy assembly bias in SDSS, enabling us to derive an accurate and stringent constraint on the latter. Our fiducial ELUCID constraint, for galaxies above a stellar mass threshold $M_*{=}10^{10.2}\,h^{-2}\,M_\odot$, is $\mathcal{Q}_\mathrm{cen}{=}{-}0.09\pm{0.05}$ and $\mathcal{Q}_\mathrm{sat}{=}0.09\pm{0.10}$, indicating no evidence for a significant~($>2\sigma$) galaxy assembly bias in the local Universe probed by SDSS. Finally, our method provides a promising path to the robust modelling of the galaxy-halo connection within future surveys like DESI and PFS.


title: Modified virial theorem for highly magnetized white dwarfs

abstract: Generally the virial theorem provides a relation between various components of energy integrated over a system. This helps us to understand the underlying equilibrium. Based on the virial theorem we can estimate, for example, the maximum allowed magnetic field in a star. Recent studies have proposed the existence of highly magnetized white dwarfs, with masses significantly higher than the Chandrasekhar limit. Surface magnetic fields of such white dwarfs could be more than 10^9 G with the central magnitude several orders higher. These white dwarfs could be significantly smaller in size than their ordinary counterparts (with surface fields restricted to about 10^9 G). In this paper we reformulate the virial theorem for non-rotating, highly magnetized white dwarfs (B-WDs) in which, unlike in previous formulations, the contribution of the magnetic pressure to the magnetohydrostatic balance cannot be neglected. Along with the new equation of magnetohydrostatic equilibrium, we approach the problem by invoking magnetic flux conservation and by varying the internal magnetic field with the matter density as a power law. Either of these choices are supported by previous independent work and neither violates any important physics. They are useful while there is no prior knowledge of field profile within a white dwarf. We then compute the modified gravitational, thermal and magnetic energies and examine how the magnetic pressure influences the properties of such white dwarfs. Based on our results we predict important properties of these B-WDs, which turn out to be independent of our chosen field profiles.


title: Accounting for object detection bias in weak gravitational lensing studies

abstract: Weak lensing by large-scale structure is a powerful probe of cosmology if the apparent alignments in the shapes of distant galaxies can be accurately measured. Most studies have therefore focused on improving the fidelity of the shape measurements themselves, but the preceding step of object detection has been largely ignored. In this paper we study the impact of object detection for a Euclid-like survey and show that it leads to biases that exceed requirements for the next generation of cosmic shear surveys. In realistic scenarios, blending of galaxies is an important source of detection bias. We find that MetaDetection is able to account for blending, leading to average multiplicative biases that meet requirements for Stage IV surveys, provided a sufficiently accurate model for the point spread function is available. Further work is needed to estimate the performance for actual surveys. Combined with sufficiently realistic image simulations, this provides a viable way forward towards accurate shear estimates for Stage IV surveys.


title: Nonminimal Couplings and the Forgotten Field of Axion Inflation

abstract: We study the multifield dynamics of axion models nonminimally coupled to gravity. As usual, we consider a canonical $U(1)$ symmetry-breaking model in which the axion is the phase of a complex scalar field. If the complex scalar field has a nonminimal coupling to gravity, then the (oft-forgotten) radial component can drive a phase of inflation prior to an inflationary phase driven by the axion field. In this setup, the mass of the axion field is dependent on the radial field because of the nonminimal coupling, and the axion remains extremely light during the phase of radial inflation. As the radial field approaches the minimum of its potential, there is a transition to a second phase of inflation during which the axion field contributes a substantial fraction of the system's total energy density, yielding a phase of "natural inflation." This system exhibits ultra-light isocurvature perturbations, which are converted to adiabatic perturbations by a turning field-space trajectory. For models wherein the CMB pivot scale exited the horizon well before the turn, this acts to suppress the tensor-to-scalar ratio $r$, without generating CMB non-Gaussianity or observable isocurvature perturbations. Finally, we note that for certain parameters the interaction strength between axion and gauge fields can be suppressed during the first phase of inflation relative to its value during the second phase by several orders of magnitude. This decouples the constraints on the inflationary production of gauge fields (e.g., from primordial black holes) from the constraints on their production during (p)reheating.


title: Improving Helium Abundance Determinations with Leo P as a Case Study

abstract: Currently, the primordial helium abundance is best estimated through spectroscopic observations of H II regions in metal-poor galaxies. However these determinations are limited by several systematic uncertainties which ultimately limit our ability to accurately ascertain the primordial abundance. In this study, we improve the methodologies for solving for the reddening, the emission contributions from collisional excitation of the H I atoms, the effects underlying absorption in the H I and He I emission lines, and the treatment of the blended H I and He I emission at $\lambda$3889 with the aim of lowering the systematic uncertainties in helium abundance determinations. To apply these methods, we have obtained observations of the He I $\lambda$10830 emission line in the brightest H II region in the extremely metal-poor (3$\%$ Z$_{\odot}$) galaxy Leo P with the LUCI1 instrument on the LBT. We combine this measurement with previous MODS/LBT observations to derive an improved helium abundance. In doing so, our present analysis results in a decrease in the uncertainty in the helium abundance of Leo P by approximately 70%. This result is combined with data from other observations to estimate the primordial helium mass fraction, Y$_{p}$ $=$ 0.2453 $\pm$ 0.0034.


title: Adatom generated spatially oscillating strain fields on isotropic and other substrates

abstract: An extended elastic eigenvector approach for adatom interactions is applied to elastically isotropic substrates. The oscillating interactions between adatom monomers or dimers are caused by strain fields in the substrate. Expansive or contractive strain fields are surrounding adatoms in form of quasi circular rings. The exact shape of the ring structure depends on the elastic constants of the substrate and especially on the surface Brillouin zone. Interactions on isotropic substrates are depicted and compared with an anisotropic example. Overlays of 2-D plots can help to identify adatom structures caused by elastic effects and to prepare numerical calculations of their elastic energies. The consequences of different sorts of adatoms with different interaction sign are sketched, large short distance attraction could occur in such cases. Limitations of the method are discussed and possible applications are formulated.


title: Evidence for galaxy assembly bias in BOSS CMASS redshift-space galaxy correlation function

abstract: Building accurate and flexible galaxy-halo connection models is crucial in modeling galaxy clustering on non-linear scales. Recent studies have found that halo concentration by itself cannot capture the full galaxy assembly bias effect and that the local environment of the halo can be an excellent indicator of galaxy assembly bias. In this paper, we propose an extended halo occupation distribution model (HOD) that includes both a concentration-based assembly bias term and an environment-based assembly bias term. We use this model to achieve a good fit (chi^2/DoF = 1.35) on the 2D redshift-space 2-point correlation function (2PCF) of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy sample. We find that the inclusion of both assembly bias terms is strongly favored by the data and the standard 5-parameter HOD is strongly rejected. More interestingly, the redshift-space 2PCF drives the assembly bias parameters in a way that preferentially assigns galaxies to lower mass halos. This results in galaxy-galaxy lensing predictions that are within 1sigma agreement with the observation, alleviating the perceived tension between galaxy clustering and lensing. We also showcase a consistent 3-5sigma preference for a positive environment-based assembly bias that persists over variations in the fit. We speculate that the environmental dependence might be driven by underlying processes such as mergers and feedback, but might also be indicative of a larger halo boundaries such as the splashback radius. Regardless, this work highlights the importance of building flexible galaxy-halo connection models and demonstrates the extra constraining power of the redshift-space 2PCF.


title: New bounds on the size of Nearly Perfect Matchings in almost regular hypergraphs

abstract: Let $H$ be a $k$-uniform $D$-regular simple hypergraph on $N$ vertices. Based on an analysis of the R\"odl nibble, Alon, Kim and Spencer (1997) proved that if $k \ge 3$, then $H$ contains a matching covering all but at most $ND^{-1/(k-1)+o(1)}$ vertices, and asked whether this bound is tight. In this paper we improve their bound by showing that for all $k > 3$, $H$ contains a matching covering all but at most $ND^{-1/(k-1)-\eta}$ vertices for some $\eta = \Theta(k^{-3}) > 0$, when $N$ and $D$ are sufficiently large. Our approach consists of showing that the R\"odl nibble process not only constructs a large matching but it also produces many well-distributed `augmenting stars' which can then be used to significantly improve the matching constructed by the R\"odl nibble process. Based on this, we also improve the results of Kostochka and R\"odl (1998) and Vu (2000) on the size of matchings in almost regular hypergraphs with small codegree. As a consequence, we improve the best known bounds on the size of large matchings in combinatorial designs with general parameters. Finally, we improve the bounds of Molloy and Reed (2000) on the chromatic index of hypergraphs with small codegree (which can be applied to improve the best known bounds on the chromatic index of Steiner triple systems and more general designs).


title: The Effects of Asymmetric Dark Matter on Stellar Evolution I: Spin-Dependent Scattering

abstract: Most of the dark matter (DM) search over the last few decades has focused on WIMPs, but the viable parameter space is quickly shrinking. Asymmetric Dark Matter (ADM) is a WIMP-like DM candidate with slightly smaller masses and no present day annihilation, meaning that stars can capture and build up large quantities. The captured ADM can transport energy through a significant volume of the star. We investigate the effects of spin-dependent ADM energy transport on stellar structure and evolution in stars with $0.9 \leq M_{\star}/\mathrm{M}_{\odot} \leq 5.0$ in varying DM environments. We wrote a MESA module that calculates the capture of DM and the subsequent energy transport within the star. We fix the DM mass to 5 GeV and the cross section to $10^{-37} \mathrm{cm{^2}}$, and study varying environments by scaling the DM capture rate. For stars with radiative cores ($M_{\star} \lesssim 1.3\ \mathrm{M}_{\odot}$), the presence of ADM flattens the temperature and burning profiles in the core and increases MS ($X_c > 10^{-3}$) lifetimes by up to $\sim 20\%$. We find that strict requirements on energy conservation are crucial to the simulation of ADM's effects on these stars. In higher-mass stars, ADM energy transport shuts off core convection, limiting available fuel and shortening MS lifetimes by up to $\sim 40\%$. This may translate to changes in the luminosity and effective temperature of the MS turnoff in population isochrones. The tip of the red giant branch may occur at lower luminosities. The effects are largest in DM environments with high densities and/or low velocity dispersions, making dwarf and early forming galaxies most likely to display the effects.

categories: astro-ph.GA astro-ph.CO astro-ph.SR hep-ph


title: FastVC: Fast Voice Conversion with non-parallel data

abstract: This paper introduces FastVC, an end-to-end model for fast Voice Conversion (VC). The proposed model can convert speech of arbitrary length from multiple source speakers to multiple target speakers. FastVC is based on a conditional AutoEncoder (AE) trained on non-parallel data and requires no annotations at all. This model's latent representation is shown to be speaker-independent and similar to phonemes, which is a desirable feature for VC systems. While the current VC systems primarily focus on achieving the highest overall speech quality, this paper tries to balance the development concerning resources needed to run the systems. Despite the simple structure of the proposed model, it outperforms the VC Challenge 2020 baselines on the cross-lingual task in terms of naturalness.


title: Machine Learning for Gas and Oil Exploration

abstract: Drilling boreholes for gas and oil extraction is an expensive process and profitability strongly depends on characteristics of the subsurface. As profitability is a key success factor, companies in the industry utilise well logs to explore the subsurface beforehand. These well logs contain various characteristics of the rock around the borehole, which allow petrophysicists to determine the expected amount of contained hydrocarbon. However, these logs are often incomplete and, as a consequence, the subsequent analyses cannot exploit the full potential of the well logs. In this paper we demonstrate that Machine Learning can be applied to \emph{fill in the gaps} and estimate missing values. We investigate how the amount of training data influences the accuracy of prediction and how to best design regression models (Gradient Boosting and neural network) to obtain optimal results. We then explore the models' predictions both quantitatively, tracking the prediction error, and qualitatively, capturing the evolution of the measured and predicted values for a given property with depth. Combining the findings has enabled us to develop a predictive model that completes the well logs, increasing their quality and potential commercial value.


title: A new upper bound for the multiple knapsack problem

abstract: In this paper, a new upper bound for the Multiple Knapsack Problem (MKP) is proposed, based on the idea of relaxing MKP to a {\em Bounded Sequential Multiple Knapsack Problem}, i.e., a multiple knapsack problem in which item sizes are divisible. Such a relaxation, called sequential relaxation, is obtained by suitably replacing the items of a MKP instance with items with divisible sizes. Experimental results on benchmark instances show that the upper bound is effective when the ratio between the number of items and the number of knapsacks is small.


title: Ribbonlength of families of folded ribbon knots

abstract: We study Kauffman's model of folded ribbon knots: knots made of a thin strip of paper folded flat in the plane. The folded ribbonlength is the length to width ratio of such a ribbon knot. We give upper bounds on the folded ribbonlength of 2-bridge, $(2,q)$ torus, twist, and pretzel knots, and these upper bounds turn out to be linear in the crossing number. We give a new way to fold $(p,q)$ torus knots and show that their folded ribbonlength is bounded above by $2p$. This means, for example, that the trefoil knot can be constructed with a folded ribbonlength of 6. We then show that any $(p,q)$ torus knot $K$ with $p\geq q>2$ has a constant $c>0$, such that the folded ribbonlength is bounded above by $c\cdot Cr(K)^{1/2}$. This provides an example of an upper bound on folded ribbonlength that is sub-linear in crossing number.


title: Remarkable Predictive Power of the Modified Long Wavelength Approximation

abstract: The modified long-wavelength approximation (MLWA), a next order approximation beyond the Rayleigh limit, has been applied usually only to the dipole $\ell=1$ contribution and for the range of size parameters $x$ not exceeding $x\lesssim 1$ to estimate far- and near-field electromagnetic properties of plasmonic nanoparticles. Provided that the MLWA functional form for the $T$-matrix in a given channel $\ell$ is limited to the ratio $T\sim iR/(F+D-iR)$, where $F$ is the familiar size-independent Fr\"ohlich term and $R\sim {\cal O}(x^{2\ell+1})$ is a radiative reaction term, there is a one-parameter freedom in selecting the dynamic depolarization term $D\sim {\cal O}(x^2)$ which preserves the fundamental feature of the MLWA that its predictions coincide with those of the Mie theory up to the order ${\cal O}(x^2)$. By exploiting this untapped design freedom, we demonstrate on a number of different metals (Ag, Al, Au, Mg), and using real material data, that the MLWA may surprisingly yield very accurate results for plasmonic spheres both for (i) $x$ up to $\gtrsim 1$ and beyond, and (ii) higher order multipoles ($\ell>1$), essentially doubling its expected range of validity. Because the MLWA obviates the need of using spherical Bessel and Hankel functions and allows for an intuitive description of (nano)particle properties in terms of a driven damped harmonic oscillator parameters, a significantly improved analysis and understanding of nanoparticle scattering and near-field properties can be achieved.


title: MatDRAM: A pure-MATLAB Delayed-Rejection Adaptive Metropolis-Hastings Markov Chain Monte Carlo Sampler

abstract: Markov Chain Monte Carlo (MCMC) algorithms are widely used for stochastic optimization, sampling, and integration of mathematical objective functions, in particular, in the context of Bayesian inverse problems and parameter estimation. For decades, the algorithm of choice in MCMC simulations has been the Metropolis-Hastings (MH) algorithm. An advancement over the traditional MH-MCMC sampler is the Delayed-Rejection Adaptive Metropolis (DRAM). In this paper, we present MatDRAM, a stochastic optimization, sampling, and Monte Carlo integration toolbox in MATLAB which implements a variant of the DRAM algorithm for exploring the mathematical objective functions of arbitrary-dimensions, in particular, the posterior distributions of Bayesian models in data science, Machine Learning, and scientific inference. The design goals of MatDRAM include nearly-full automation of MCMC simulations, user-friendliness, fully-deterministic reproducibility, and the restart functionality of simulations. We also discuss the implementation details of a technique to automatically monitor and ensure the diminishing adaptation of the proposal distribution of the DRAM algorithm and a method of efficiently storing the resulting simulated Markov chains. The MatDRAM library is open-source, MIT-licensed, and permanently located and maintained as part of the ParaMonte library at https://github.com/cdslaborg/paramonte.

categories: physics.data-an astro-ph.IM cs.CE q-bio.QM stat.AP


title: PoinT-5: Pointer Network and T-5 based Financial NarrativeSummarisation

abstract: Companies provide annual reports to their shareholders at the end of the financial year that describes their operations and financial conditions. The average length of these reports is 80, and it may extend up to 250 pages long. In this paper, we propose our methodology PoinT-5 (the combination of Pointer Network and T-5 (Test-to-text transfer Transformer) algorithms) that we used in the Financial Narrative Summarisation (FNS) 2020 task. The proposed method uses pointer networks to extract important narrative sentences from the report, and then T-5 is used to paraphrase extracted sentences into a concise yet informative sentence. We evaluate our method using ROUGE-N (1,2), L, and SU4. The proposed method achieves the highest precision scores in all the metrics and highest F1 scores in ROUGE1, and LCS and the only solution to cross the MUSE solution baseline in ROUGE-LCS metrics.


title: From "universal" profiles to "universal" scaling laws in X-ray galaxy clusters

abstract: As the end products of the hierarchical process of cosmic structure formation, galaxy clusters present some predictable properties, like those mostly driven by gravity, and some others, more affected by astrophysical dissipative processes, that can be recovered from observations and that show remarkable "universal" behaviour once rescaled by halo mass and redshift. However, a consistent picture that links these universal radial profiles and the integrated values of the thermodynamical quantities of the intracluster medium, also quantifying the deviations from the standard self-similar gravity-driven scenario, has to be demonstrated. In this work, we use a semi-analytic model based on a universal pressure profile in hydrostatic equilibrium within a cold dark matter halo with a defined relation between mass and concentration to reconstruct the scaling laws between the X-ray properties of galaxy clusters. We also quantify any deviation from the self-similar predictions in terms of temperature dependence of a few physical quantities such as the gas mass fraction, the relation between spectroscopic temperature and its global value, and, if present, the hydrostatic mass bias. This model allows to reconstruct both the observed profiles and the scaling laws between integrated quantities. We use the Planck-selected ESZ sample to calibrate the predicted scaling laws between gas mass, temperature, luminosity and total mass. Our universal model reproduces well the observed thermodynamic properties and provides a way to interpret the observed deviations from the standard self-similar behaviour. By combining these results with the constraints on the observed $Y_{SZ}-T$ relation, we show how we can quantify the level of gas clumping affecting the studied sample, estimate the clumping-free gas mass fraction, and suggest the average level of hydrostatic bias present.


title: Luminous solar neutrinos I: Dipole portals

abstract: Solar neutrinos upscattering inside the Earth can source unstable particles that can decay inside terrestrial detectors. Contrary to naive expectations we show that when the decay length is much shorter than the radius of the \emph{Earth} (rather than the detector), the event rate is independent of the decay length. In this paper we study a transition dipole operator (neutrino dipole portal) and show that Borexino's existing data probes previously untouched parameter space in the 0.5--20 MeV regime, complementing recent cosmological and supernova bounds. We briefly comment on similarities and differences with luminous dark matter and comment on future prospects for analogous signals stemming from atmospheric neutrinos.


title: Anisotropy and multifractal analysis of turbulent velocity and temperature in the roughness sublayer of a forested canopy

abstract: Anisotropy and multifractality in velocity and temperature time series sampled at multiple heights in the roughness sublayer (RSL) over a boreal mixed-coniferous forest are reported. In particular, a turbulent-stress invariant analysis along with a scalewise version of it are conducted to elucidate the nature of relaxation of large-scale anisotropy to quasi-isotropic states at small scales. As the return to isotropy is linked to nonlinear interactions and correlations between different fluctuating velocity components across scales, we study the velocity and temperature time series by using multifractal detrended fluctuation analysis and multiscale multifractal analysis to assess the effects of thermal stratification and surface roughness on turbulence in the RSL. The findings are compared so as to quantify the anisotropy and multifractality ubiquitous to RSL turbulent flow. As we go up in the RSL, (a) the length scale at which return to isotropy commences increases because of the weakening of the surface effects and (b) the largest scales become increasingly anisotropic. The anisotropy in multifractal exponents for the velocity fluctuations is diminished when we use the extended-self-similarity procedure to extract the multifractal-exponent ratios.


title: Eigenmode analysis of membrane stability in inviscid flow

abstract: We study the instability of a thin membrane (of zero bending rigidity) to out-of-plane deflections, when the membrane is immersed in an inviscid fluid flow and sheds a trailing vortex-sheet wake. We solve the nonlinear eigenvalue problem iteratively with large ensembles of initial guesses, for three canonical boundary conditions---both ends fixed, one end fixed and one free, and both free. Over several orders of magnitude of membrane mass density, we find instability by divergence or flutter (particularly at large mass density, or with one or both ends free). The most unstable eigenmodes generally become "wavier" at smaller mass density and smaller tension, but with regions of nonmonotonic behavior. We find good quantitative agreement with unsteady time-stepping simulations at small amplitude, but only qualitative similarities with the eventual steady-state large-amplitude motions.


title: A Fully Tensorized Recurrent Neural Network

abstract: Recurrent neural networks (RNNs) are powerful tools for sequential modeling, but typically require significant overparameterization and regularization to achieve optimal performance. This leads to difficulties in the deployment of large RNNs in resource-limited settings, while also introducing complications in hyperparameter selection and training. To address these issues, we introduce a "fully tensorized" RNN architecture which jointly encodes the separate weight matrices within each recurrent cell using a lightweight tensor-train (TT) factorization. This approach represents a novel form of weight sharing which reduces model size by several orders of magnitude, while still maintaining similar or better performance compared to standard RNNs. Experiments on image classification and speaker verification tasks demonstrate further benefits for reducing inference times and stabilizing model training and hyperparameter selection.


title: Nuclear Spin Assisted Magnetic Field Angle Sensing

abstract: Quantum sensing exploits the strong sensitivity of quantum systems to measure small external signals. The nitrogen-vacancy (NV) center in diamond is one of the most promising platforms for real-world quantum sensing applications, predominantly used as a magnetometer. However, its magnetic field sensitivity vanishes when a bias magnetic field acts perpendicular to the NV axis. Here, we introduce a novel sensing strategy assisted by the nitrogen nuclear spin that uses the entanglement between the electron and nuclear spins to restore the magnetic field sensitivity. This, in turn, allows us to detect small changes in the magnetic field angle relative to the NV axis. Furthermore, based on the same underlying principle, we show that the NV coupling strength to magnetic noise, and hence its coherence time, exhibits a strong asymmetric angle dependence. This allows us to uncover the directional properties of the local magnetic environment and to realize maximal decoupling from anisotropic noise.


title: Connected Generalized Inverse Limits and Intermediate Value Property

abstract: In this paper, we consider inverse limits of $[0,1]$ using upper semicontinuous set-valued functions. We introduce two generalizations of the Intermediate Value Property and prove that inverse limits with upper semicontinuous set-valued bonding functions are connected if the bonding functions are surjective, have connected graphs, and have either generalization of the Intermediate Value Property. Examples are given to demonstrate that if any of the conditions is dropped, the result does not hold in general. An example is given to show that an inverse limit may be connected even if the bonding functions do not have either Intermediate Value Property. Further, we compare the structures of set-valued functions with the two types of the Intermediate Value Property.


title: Multiscale Elliptic PDEs Upscaling and Function Approximation via Subsampled Data

abstract: There is an intimate connection between numerical upscaling of multiscale PDEs and scattered data approximation of heterogeneous functions: the coarse variables selected for deriving an upscaled equation (in the former) correspond to the sampled information used for approximation (in the latter). As such, both problems can be thought of as recovering a target function based on some coarse data that are either artificially chosen by an upscaling algorithm, or determined by some physical measurement process. The purpose of this paper is then to study that, under such a setup and for a specific elliptic problem, how the lengthscale of the coarse data, which we refer to as the subsampled lengthscale, influences the accuracy of recovery, given limited computational budgets. Our analysis and experiments identify that, reducing the subsampling lengthscale may improve the accuracy, implying a guiding criterion for coarse-graining or data acquisition in this computationally constrained scenario, especially leading to direct insights for the implementation of the Gamblets method in the numerical homogenization literature. Moreover, reducing the lengthscale to zero may lead to a blow-up of approximation error if the target function does not have enough regularity, suggesting the need for a stronger prior assumption on the target function to be approximated. We introduce a singular weight function to deal with it, both theoretically and numerically. This work sheds light on the interplay of the lengthscale of coarse data, the computational costs, the regularity of the target function, and the accuracy of approximations and numerical simulations.


title: Memory effects, arches and polar defect ordering at the cross-over from wet to dry active nematics

abstract: We use analytic arguments and numerical solutions of the continuum, active nematohydrodynamic equations to study how friction alters the behaviour of active nematics. Concentrating on the case where there is nematic ordering in the passive limit, we show that, as the friction is increased, memory effects become more prominent and $+1/2$ topological defects leave increasingly persistent trails in the director field as they pass. The trails are preferential sites for defect formation and they tend to impose polar order on any new $+1/2$ defects. In the absence of noise and for high friction, it becomes very difficult to create defects, but trails formed by any defects present at the beginning of the simulations persist and organise into parallel arch-like patterns in the director field. We show aligned arches of equal width are approximate steady state solutions of the equations of motion which co-exist with the nematic state. We compare our results to other models in the literature, in particular dry systems with no hydrodynamics, where trails, arches and polar defect ordering have also been observed.


title: Bicycle paths, elasticae and sub-Riemannian geometry

abstract: We relate the sub-Riemannian geometry on the group of rigid motions of the plane to `bicycling mathematics'. We show that this geometry's geodesics correspond to bike paths whose front tracks are either non-inflectional Euler elasticae or straight lines, and that its infinite minimizing geodesics (or `metric lines') correspond to bike paths whose front tracks are either straight lines or `Euler's solitons' (also known as Syntractrix or Convicts' curves).


title: Orthogonal Decomposition of Tensor Trains

abstract: In this paper we study the problem of decomposing a given tensor into a tensor train such that the tensors at the vertices are orthogonally decomposable. When the tensor train has length two, and the orthogonally decomposable tensors at the two vertices are symmetric, we recover the decomposition by considering random linear combinations of slices. Furthermore, if the tensors at the vertices are symmetric and low-rank but not orthogonally decomposable, we show that a whitening procedure can transform the problem into the orthogonal case. When the tensor network has length three or more and the tensors at the vertices are symmetric and orthogonally decomposable, we provide an algorithm for recovering them subject to some rank conditions. Finally, in the case of tensor trains of length two in which the tensors at the vertices are orthogonally decomposable but not necessarily symmetric, we show that the decomposition problem reduces to the novel problem of decomposing a matrix into an orthogonal matrix multiplied by diagonal matrices on either side. We provide and compare two solutions, one based on Sinkhorn's theorem and one on Procrustes' algorithm. We conclude with a multitude of open problems in linear and multilinear algebra that arose in our study.


title: Efficient Real-Time Radial Distortion Correction for UAVs

abstract: In this paper we present a novel algorithm for onboard radial distortion correction for unmanned aerial vehicles (UAVs) equipped with an inertial measurement unit (IMU), that runs in real-time. This approach makes calibration procedures redundant, thus allowing for exchange of optics extemporaneously. By utilizing the IMU data, the cameras can be aligned with the gravity direction. This allows us to work with fewer degrees of freedom, and opens up for further intrinsic calibration. We propose a fast and robust minimal solver for simultaneously estimating the focal length, radial distortion profile and motion parameters from homographies. The proposed solver is tested on both synthetic and real data, and perform better or on par with state-of-the-art methods relying on pre-calibration procedures.


title: Signed Distance Laplacian Matrices for Signed Graphs

abstract: A signed graph is a graph whose edges are labeled either positive or negative. Corresponding to the two signed distance matrices defined for signed graphs, we define two signed distance laplacian matrices. We characterize balance in signed graphs using these matrices and find signed distance laplacian spectra of some classes of unbalanced signed graphs.


title: Gaussian MRF Covariance Modeling for Efficient Black-Box Adversarial Attacks

abstract: We study the problem of generating adversarial examples in a black-box setting, where we only have access to a zeroth order oracle, providing us with loss function evaluations. Although this setting has been investigated in previous work, most past approaches using zeroth order optimization implicitly assume that the gradients of the loss function with respect to the input images are \emph{unstructured}. In this work, we show that in fact substantial correlations exist within these gradients, and we propose to capture these correlations via a Gaussian Markov random field (GMRF). Given the intractability of the explicit covariance structure of the MRF, we show that the covariance structure can be efficiently represented using the Fast Fourier Transform (FFT), along with low-rank updates to perform exact posterior estimation under this model. We use this modeling technique to find fast one-step adversarial attacks, akin to a black-box version of the Fast Gradient Sign Method~(FGSM), and show that the method uses fewer queries and achieves higher attack success rates than the current state of the art. We also highlight the general applicability of this gradient modeling setup.