title: Topology and Routing Problems: The Circular Frame

abstract: In this work, we solve the problem of finding non-intersecting paths between points on a plane with a new approach by borrowing ideas from geometric topology, in particular, from the study of polygonal schema in mathematics. We use a topological transformation on the 2-dimensional planar routing environment that simplifies the routing problem into a problem of connecting points on a circle with straight line segments that do not intersect in the interior of the circle. These points are either the points that need to be connected by non-intersecting paths or special `reference' points that parametrize the topology of the original environment prior to the transformation. When all the necessary points on the circle are fully connected, the transformation is reversed such that the line segments combine to become the non-intersecting paths that connect the start and end points in the original environment. We interpret the transformed environment in which the routing problem is solved as a new data structure where any routing problem can be solved efficiently. We perform experiments and show that the routing time and success rate of the new routing algorithm outperforms the ones for the A*-algorithm.


title: Effects of Magnetic Fields in Hot White Dwarfs

abstract: In this work, we study the effects of temperature on magnetic white dwarfs. We model their interior as a nuclei lattice surrounded by a relativistic free Fermi gas of electrons, accounting for effects from temperature, Landau levels and anomalous magnetic moment. We find that, at low densities (corresponding to the outer regions of star), both temperature and magnetic field effects play an important role in the calculation of microscopic thermodynamical quantities. To study macroscopic stellar structures within a general-relativistic approach, we solve numerically the coupled Einstein's-Maxwell's equations for fixed entropy per particle configurations and discuss how temperature affects stellar magnetic field profiles, masses and radii.


title: Hierarchical Graph Neural Networks

abstract: Over the recent years, Graph Neural Networks have become increasingly popular in network analytic and beyond. With that, their architecture noticeable diverges from the classical multi-layered hierarchical organization of the traditional neural networks. At the same time, many conventional approaches in network science efficiently utilize the hierarchical approaches to account for the hierarchical organization of the networks, and recent works emphasize their critical importance. This paper aims to connect the dots between the traditional Neural Network and the Graph Neural Network architectures as well as the network science approaches, harnessing the power of the hierarchical network organization. A Hierarchical Graph Neural Network architecture is proposed, supplementing the original input network layer with the hierarchy of auxiliary network layers and organizing the computational scheme updating the node features through both - horizontal network connections within each layer as well as the vertical connection between the layers. It enables simultaneous learning of the individual node features along with the aggregated network features at variable resolution and uses them to improve the convergence and stability of the individual node feature learning. The proposed Hierarchical Graph Neural network architecture is successfully evaluated on the network embedding and modeling as well as network classification, node labeling, and community tasks and demonstrates increased efficiency in those.


title: GADTs, Functoriality, Parametricity: Pick Two

abstract: GADTs can be represented either as their Church encodings a la Atkey, or as fixpoints a la Johann and Polonsky. While a GADT represented as its Church encoding need not support a map function satisfying the functor laws, the fixpoint representation of a GADT must support such a map function even to be well-defined. The two representations of a GADT thus need not be the same in general. This observation forces a choice of representation of data types in languages supporting GADTs. In this paper we show that choosing whether to represent data types as their Church encodings or as fixpoints determines whether or not a language supporting GADTs can have parametric models. This choice thus has important consequences for how we can program with, and reason about, these advanced data types.


title: Deep Coded Aperture Design: An End-to-End Approach for Computational Imaging Tasks

abstract: Covering from photography to depth and spectral estimation, diverse computational imaging (CI) applications benefit from the versatile modulation of coded apertures (CAs). The light wave fields as space, time, or spectral can be modulated to obtain projected encoded information at the sensor that is then decoded by efficient methods, such as the modern deep learning decoders. Despite the CA can be fabricated to produce an analog modulation, a binary CA is mostly preferred since easier calibration, higher speed, and lower storage are achieved. As the performance of the decoder mainly depends on the structure of the CA, several works optimize the CA ensembles by customizing regularizers for a particular application without considering critical physical constraints of the CAs. This work presents an end-to-end (E2E) deep learning-based optimization of CAs for CI tasks. The CA design method aims to cover a wide range of CI problems easily changing the loss function of the deep approach. The designed loss function includes regularizers to fulfill the widely used sensing requirements of the CI applications. Mainly, the regularizers can be selected to optimize the transmittance, the compression ratio, and the correlation between measurements, while a binary CA solution is encouraged, and the performance of the CI task is maximized in applications such as restoration, classification, and semantic segmentation.


title: Meta-validation of bipartite network projections

abstract: Monopartite projections of bipartite networks are useful tools for modeling indirect interactions in complex systems. The standard approach to identify significant links is statistical validation using a suitable null network model, such as the popular configuration model (CM) that constrains node degrees and randomizes everything else. However different CM formulations exist, depending on how the constraints are imposed and for which sets of nodes. Here we systematically investigate the application of these formulations in validating the same network, showing that they lead to different results even when the same significance threshold is used. Instead a much better agreement is obtained for the same density of validated links. We thus propose a meta-validation approach that allows to identify model-specific significance thresholds for which the signal is strongest, and at the same time to obtain results independent of the way in which the null hypothesis is formulated. We illustrate this procedure using data on scientific production of world countries.

categories: physics.soc-ph cond-mat.dis-nn cs.SI physics.data-an


title: Integrated modelling of microvibrations induced by Solar Array Drive Mechanism for worst-case end-to-end analysis and robust disturbance estimation

abstract: Modern and future observation Space missions face increasingly more demanding pointing performance requirements. This is accompanied with the development of larger lightweight flexible structures. This paper outlines a methodology for modeling a generic multi-body flexible spacecraft in an end-to-end fashion within the Two-Inputs Two-Outputs Port framework: from disturbance all the way to pointing performance by tacking into account all the uncertainties of each sub-system. A particular focus is dedicated to a novel generic model for Solar Drive Array Mechanisms and its harmonic disturbances produced by the micro-stepping driver and the gearbox imperfections. The proposed model is validated with the on-board telemetries of a European spacecraft. Moreover this paper shows how the proposed model-based approach can be easily used for control design and closed-loop robust performance analysis. Worst-case scenarios of mechanical parameters combination are investigated to provide robust performance certificate of the reference study case. A Linear Parameter-Varying observer, scheduled by the solar array rotor angle, is finally proposed in order to estimate the disturbance torques induced by the gearbox imperfections just using the measurements coming from a classical attitude control system: a star tracker and a gyrometer.


title: Frequency-explicit approximability estimates for time-harmonic Maxwell's equations

abstract: We consider time-harmonic Maxwell's equations set in an heterogeneous medium with perfectly conducting boundary conditions. Given a divergence-free right-hand side lying in $L^2$, we provide a frequency-explicit approximability estimate measuring the difference between the corresponding solution and its best approximation by high-order N\'ed\'elec finite elements. Such an approximability estimate is crucial in both the a priori and a posteriori error analysis of finite element discretizations of Maxwell's equations, but the derivation is not trivial. Indeed, it is hard to take advantage of high-order polynomials given that the right-hand side only exhibits $L^2$ regularity. We proceed in line with previously obtained results for the simpler setting of the scalar Helmholtz equation, and propose a regularity splitting of the solution. In turn, this splitting yields sharp approximability estimates generalizing known results for the scalar Helmoltz equation and showing the interest of high-order methods.


title: Branch cuts: writing, editing, and ramified complexities

abstract: As I was preparing my tenure application, the University of Wisconsin Board of Regents voted to redefine tenure, removing many of the institution's historical protections. Reevaluating my career priorities in light of these changes and a resurgent two-body problem, I recognized that my fundamental goal was communicating mathematical ideas. I found a new role as an editor at Mathematical Reviews, part of the American Mathematical Society. To my surprise, thinking more about my identity as a writer and editor also changed my perspective on my own sexuality and gender identity, inspiring new approaches to leadership.


title: SERVAS! Secure Enclaves via RISC-V Authenticryption Shield

abstract: Isolation is a long-standing challenge of software security. Traditional privilege rings and virtual memory are more and more augmented with concepts such as capabilities, protection keys, and powerful enclaves. At the same time, we are evidencing an increased need for physical protection, shifting towards full memory encryption schemes. This results in a complex interplay of various security mechanisms, increasing the burden for system architects and security analysts. In this work, we tackle the isolation challenge with a new isolation primitive called authenticryption shield that unifies both traditional and advanced isolation policies while offering the potential for future extensibility. At the core, we build upon an authenticated memory encryption scheme that gives cryptographic isolation guarantees and, thus, streamlines the security reasoning. We showcase the versatility of our approach by designing and prototyping SERVAS -- an innovative enclave architecture for RISC-V. Unlike current enclave systems, SERVAS facilitates efficient and secure enclave memory sharing. While the memory encryption constitutes the main overhead, entering or exiting a SERVAS enclave requires only 3.5x of a simple syscall, instead of 71x for Intel SGX.


title: Double-matched matrix decomposition for multi-view data

abstract: We consider the problem of extracting joint and individual signals from multi-view data, that is data collected from different sources on matched samples. While existing methods for multi-view data decomposition explore single matching of data by samples, we focus on double-matched multi-view data (matched by both samples and source features). Our motivating example is the miRNA data collected from both primary tumor and normal tissues of the same subjects; the measurements from two tissues are thus matched both by subjects and by miRNAs. Our proposed double-matched matrix decomposition allows to simultaneously extract joint and individual signals across subjects, as well as joint and individual signals across miRNAs. Our estimation approach takes advantage of double-matching by formulating a new type of optimization problem with explicit row space and column space constraints, for which we develop an efficient iterative algorithm. Numerical studies indicate that taking advantage of double-matching leads to superior signal estimation performance compared to existing multi-view data decomposition based on single-matching. We apply our method to miRNA data as well as data from the English Premier League soccer matches, and find joint and individual multi-view signals that align with domain specific knowledge.


title: Learning-enhanced robust controller synthesis with rigorous statistical and control-theoretic guarantees

abstract: The combination of machine learning with control offers many opportunities, in particular for robust control. However, due to strong safety and reliability requirements in many real-world applications, providing rigorous statistical and control-theoretic guarantees is of utmost importance, yet difficult to achieve for learning-based control schemes. We present a general framework for learning-enhanced robust control that allows for systematic integration of prior engineering knowledge, is fully compatible with modern robust control and still comes with rigorous and practically meaningful guarantees. Building on the established Linear Fractional Representation and Integral Quadratic Constraints framework, we integrate Gaussian Process Regression as a learning component and state-of-the-art robust controller synthesis. In a concrete robust control example, our approach is demonstrated to yield improved performance with more data, while guarantees are maintained throughout.


title: A method to coarse-grain multi-agent stochastic systems with regions of multistability

abstract: Hybrid multiscale modelling has emerged as a useful framework for modelling complex biological phenomena. However, when accounting for stochasticity in the internal dynamics of agents, these models frequently become computationally expensive. Traditional techniques to reduce the computational intensity of such models can lead to a reduction in the richness of the dynamics observed, compared to the original system. Here we use large deviation theory to decrease the computational cost of a spatially-extended multi-agent stochastic system with a region of multi-stability by coarse-graining it to a continuous time Markov chain on the state space of stable steady states of the original system. Our technique preserves the original description of the stable steady states of the system and accounts for noise-induced transitions between them. We apply the method to a bistable system modelling phenotype specification of cells driven by a lateral inhibition mechanism. For this system, we demonstrate how the method may be used to explore different pattern configurations and unveil robust patterns emerging on longer timescales. We then compare the full stochastic, coarse-grained and mean-field descriptions via pattern quantification metrics and in terms of the numerical cost of each method. Our results show that the coarse-grained system exhibits the lowest computational cost while preserving the rich dynamics of the stochastic system. The method has the potential to reduce the computational complexity of hybrid multiscale models, making them more tractable for analysis, simulation and hypothesis testing.


title: Vector-like Leptons in Light of the Cabibbo-Angle Anomaly

abstract: The first row of the Cabibbo-Cobayashi-Maskawa (CKM) matrix shows a discrepancy of $\sim $$3\,\sigma$ with unitarity, known as the "Cabibbo Angle Anomaly" (CAA). This tension can be explained via modified $W$ couplings to leptons, but in order to consistently assess the validity of this solution, a global fit including the electroweak (EW) precision constraints and the low energy tests of lepton flavour universality (LFU) is necessary. Performing such a fit for gauge-invariant dimension-6 operators, we find that even when assuming LFU, including the CKM elements $V_{us}$ and $V_{ud}$ into the electroweak fit has a relevant impact, shifting the best fit point significantly. Vector-like leptons (VLLs) are prime candidates for a corresponding UV completion since they can affect the $W$ and $Z$ couplings to leptons already at tree-level. Studying each pattern of new physics (NP) given by the six possible representations of VLLs (under the SM gauge group) we find that any single representation describes experimental data at most slighlty better than the SM hypothesis. However, allowing for several representations of VLLs at the same time, we find a simple scenario consisting of a singlet, $N$, coupling to electrons and a triplet, $\Sigma_1$, coupling to muons which not only explains the CAA but also improves the electroweak fit in such a way that its best fit point is preferred by more than $4\,\sigma$ with respect to the SM.


title: Correlation in Neuronal Calcium Spiking: Quantification based on Empirical Mutual Information Rate

abstract: Quantification of neuronal correlations in neuron populations helps us to understand neural coding rules. Such quantification could also reveal how neurons encode information in normal and disease conditions like Alzheimer's and Parkinson's. While neurons communicate with each other by transmitting spikes, there would be a change in calcium concentration within the neurons inherently. Accordingly, there would be correlations in calcium spike trains and they could have heterogeneous memory structures. In this context, estimation of mutual information rate in calcium spike trains assumes primary significance. However, such estimation is difficult with available methods which would consider longer blocks for convergence without noticing that neuronal information changes in short time windows. Against this backdrop, we propose a faster method that exploits the memory structures in pair of calcium spike trains to quantify mutual information shared between them. Our method has shown superior performance with example Markov processes as well as experimental spike trains. Such mutual information rate analysis could be used to identify signatures of neuronal behavior in large populations in normal and abnormal conditions.


title: Gamma-convergence of a gradient-flow structure to a non-gradient-flow structure

abstract: We study the asymptotic behaviour of a gradient system in a regime in which the driving energy becomes singular. For this system gradient-system convergence concepts are ineffective. We characterize the limiting behaviour in a different way, by proving $\Gamma$-convergence of the so-called energy-dissipation functional, which combines the gradient-system components of energy and dissipation in a single functional. The $\Gamma$-limit of these functionals again characterizes a variational evolution, but this limit functional is not the energy-dissipation functional of any gradient system. The system in question describes the diffusion of a particle in a one-dimensional double-well energy landscape, in the limit of small noise. The wells have different depth, and in the small-noise limit the process converges to a Markov process on a two-state system, in which jumps only happen from the higher to the lower well. This transmutation of a gradient system into a variational evolution of non-gradient type is a model for how many one-directional chemical reactions emerge as limit of reversible ones. The $\Gamma$-convergence proved in this paper both identifies the `fate' of the gradient system for these reactions and the variational structure of the limiting irreversible reactions.


title: Reconfiguring Independent Sets on Interval Graphs

abstract: We study reconfiguration of independent sets in interval graphs under the token sliding rule. We show that if two independent sets of size $k$ are reconfigurable in an $n$-vertex interval graph, then there is a reconfiguration sequence of length $\mathcal{O}(k\cdot n^2)$. We also provide a construction in which the shortest reconfiguration sequence is of length $\Omega(k^2\cdot n)$. As a counterpart to these results, we also establish that $\textsf{Independent Set Reconfiguration}$ is PSPACE-hard on incomparability graphs, of which interval graphs are a special case.


title: Influence of bulk mass distribution on orbital precession of S2 star in Yukawa gravity

abstract: In this study we investigate possible applications of observed S2 orbit around Galactic Center for constraining the Yukawa gravity at scales in the range between several tens and several thousands astronomical units (AU) to obtain graviton mass constraints. In our model we suppose that bulk distribution of matter (includes stellar cluster, interstellar gas distribution and dark matter) exists near Supermassive Black Hole (SMBH) in our Galactic Center. We obtain the values of orbital precession angle for different values of mass density of matter and we require that the value of orbital precession is the same like in General Relativity (GR). From that request we determine gravity parameter $\lambda$ and the upper value for graviton mass. We found that in the cases where the density of extended mass is higher, the maximum allowed value for parameter $\lambda$ is smaller and the upper limit for graviton mass is higher. It is due to the fact that the extended mass causes the retrograde orbital precession. We believe that this study is a very efficient tool to evaluate a gravitational potential at the Galactic Center, parameter $\lambda$ of the Yukawa gravity model, and to constrain the graviton mass.


title: ResMLP: Feedforward networks for image classification with data-efficient training

abstract: We present ResMLP, an architecture built entirely upon multi-layer perceptrons for image classification. It is a simple residual network that alternates (i) a linear layer in which image patches interact, independently and identically across channels, and (ii) a two-layer feed-forward network in which channels interact independently per patch. When trained with a modern training strategy using heavy data-augmentation and optionally distillation, it attains surprisingly good accuracy/complexity trade-offs on ImageNet. We also train ResMLP models in a self-supervised setup, to further remove priors from employing a labelled dataset. Finally, by adapting our model to machine translation we achieve surprisingly good results. We share pre-trained models and our code based on the Timm library.


title: Dynamic tariffs-based demand response in retail electricity market under uncertainty

abstract: Demand response (DR) programs play a crucial role in improving system reliability and mitigating price volatility by altering the core profile of electricity consumption. This paper proposes a game-theoretical model that captures the dynamic interplay between retailers (leaders) and consumers (followers) in a tariffs-based electricity market under uncertainty. The proposed procedure offers theoretical and economic insights by analyzing demand flexibility within a hierarchical decision-making framework. In particular, two main market configurations are examined under uncertainty: i) there exists a retailer that exercises market power over consumers, and ii) the retailer and the consumers participate in a perfect competitive game. The former case is formulated as a mathematical program with equilibrium constraints (MPEC), whereas the latter case is recast as a mixed-integer linear program (MILP). These problems are solved by deriving equivalent tractable reformulations based on the Karush-Kuhn-Tucker (KKT) optimality conditions of each agent's problem. Numerical simulations based on real data from the European Energy Exchange platform are used to illustrate the performance of the proposed methodology. The results indicate that the proposed model effectively characterizes the interactions between retailers and flexible consumers in both perfect and imperfect market structures. Under perfect competition, the economic benefits extend not only to consumers but also to overall social welfare. Conversely, in an imperfect market, retailers leverage consumer flexibility to enhance their expected profits, transferring the risk of uncertainty to end-users. Additionally, the degree of consumer flexibility and their valuation of electricity consumption play significant roles in shaping market outcomes.


title: Covariant quantum kernels for data with group structure

abstract: The use of kernel functions is a common technique to extract important features from data sets. A quantum computer can be used to estimate kernel entries as transition amplitudes of unitary circuits. Quantum kernels exist that, subject to computational hardness assumptions, cannot be computed classically. It is an important challenge to find quantum kernels that provide an advantage in the classification of real-world data. We introduce a class of quantum kernels that can be used for data with a group structure. The kernel is defined in terms of a unitary representation of the group and a fiducial state that can be optimized using a technique called kernel alignment. We apply this method to a learning problem on a coset-space that embodies the structure of many essential learning problems on groups. We implement the learning algorithm with $27$ qubits on a superconducting processor.


title: Universal features in lifetime distribution of clusters in hydrogen bonding liquids

abstract: Hydrogen bonding liquids, typically water and alcohols, are known to form labile structures (network, chains, etc...), hence the lifetime of such structures is an important microscopic parameter, which can be calculated in computer simulations. Since these cluster entities are mostly statistical in nature, one would expect that, in the short time regime, their lifetime distribution would be a broad Gaussian-like function of time, with a single maximum representing their mean lifetime, and weakly dependent on criteria such as the bonding distance and angle, much similarly to non-hydrogen bonding simple liquids, while the long time part is known to have some power law dependence. Unexpectedly, all the hydrogen bonding liquids studied herein, namely water and alcohols, display highly hierarchic three types of specific lifetimes, in the sub-picosecond range 0-0.5ps The dominant lifetime very strongly depends on the bonding distance criterion and is related to hydrogen bonded pairs. This mode is absent in non-H-bonding simple liquids. The secondary and tertiary mean lifetimes are related to clusters, and are nearly independent on the bonding criterion. Of these two lifetimes, only the first one can be related to that of simple liquids, which poses the question of the nature of the third life time. The study of acohols reveals that this 3rd lifetime is related to the topology of H-bonded clusters, and that its distribution may be also affected by the alkyl tail surrounding "bath". This study reveals that hydrogen bonding liquids have a universal hierarchy of hydrogen bonding lifetimes with a timescale regularity across very different types, and which depend on the topology of the cluster structures


title: Microscopics of Quantum Annealing in the Disordered Dipolar Ising Ferromagnet LiHo$_{x}$Y$_{1-x}$F$_4$

abstract: Quantum annealing (QA) refers to an optimization process that uses quantum fluctuations to find the global minimum of a rugged energy landscape with many local minima. Conceptually, QA is often framed in the context of the disordered transverse field Ising model, in which a magnetic field applied perpendicular to the Ising axis tunes the quantum fluctuations and enables the system to tunnel through energy barriers, and hence reach the ground state more quickly. A solid state material closely related to this model, LiHo$_{0.45}$Y$_{0.55}$F$_4$, was shown to exhibit faster dynamics after a QA protocol, compared to thermal annealing (TA), but little is known about the actual process of optimization involved or the nature of the spin correlations in the state that is reached. Here, we report on the microscopics of QA in this material using diffuse magnetic neutron scattering. Comparing a QA to a TA protocol that reach the same end-point, we find very similar final diffuse scattering which consists of pinch-point scattering, largely consistent with critical scattering near a phase boundary in the dipolar Ising ferromagnetic model. However, comparing the time evolution at the end of the protocols, we find that the spin correlations evolve more significantly after TA, suggesting that QA produces a state closer to equilibrium. We also observe experimental evidence that the transverse field produces random fields, which had been previously predicted for this material and studied in other contexts. Thus, while the material does exhibit a quantum speedup under quantum annealing conditions, it is not a simple annealing problem; the energy landscape being optimized is changing as the optimization proceeds.


title: A Benchmarking on Cloud based Speech-To-Text Services for French Speech and Background Noise Effect

abstract: This study presents a large scale benchmarking on cloud based Speech-To-Text systems: {Google Cloud Speech-To-Text}, {Microsoft Azure Cognitive Services}, {Amazon Transcribe}, {IBM Watson Speech to Text}. For each systems, 40158 clean and noisy speech files about 101 hours are tested. Effect of background noise on STT quality is also evaluated with 5 different Signal-to-noise ratios from 40dB to 0dB. Results showed that {Microsoft Azure} provided lowest transcription error rate $9.09\%$ on clean speech, with high robustness to noisy environment. {Google Cloud} and {Amazon Transcribe} gave similar performance, but the latter is very limited for time-constraint usage. Though {IBM Watson} could work correctly in quiet conditions, it is highly sensible to noisy speech which could strongly limit its application in real life situations.


title: Interpretable machine learning for high-dimensional trajectories of aging health

abstract: We have built a computational model for individual aging trajectories of health and survival, which contains physical, functional, and biological variables, and is conditioned on demographic, lifestyle, and medical background information. We combine techniques of modern machine learning with an interpretable interaction network, where health variables are coupled by explicit pair-wise interactions within a stochastic dynamical system. Our dynamic joint interpretable network (DJIN) model is scalable to large longitudinal data sets, is predictive of individual high-dimensional health trajectories and survival from baseline health states, and infers an interpretable network of directed interactions between the health variables. The network identifies plausible physiological connections between health variables as well as clusters of strongly connected health variables. We use English Longitudinal Study of Aging (ELSA) data to train our model and show that it performs better than multiple dedicated linear models for health outcomes and survival. We compare our model with flexible lower-dimensional latent-space models to explore the dimensionality required to accurately model aging health outcomes. Our DJIN model can be used to generate synthetic individuals that age realistically, to impute missing data, and to simulate future aging outcomes given arbitrary initial health states.


title: Investigating Degradation Modes in Zn-AgO Aqueous Batteries with $\textit{In-Situ}$ X-ray Micro Computed Tomography

abstract: To meet growing energy demands, degradation mechanisms of energy storage devices must be better understood. As a non-destructive tool, X-ray Computed Tomography (CT) has been increasingly used by the battery community to perform $\textit{in-situ}$ experiments that can investigate dynamic phenomena. However, few have used X-ray CT to study representative battery systems over long cycle lifetimes (>100 cycles). Here, we report the $\textit{in-situ}$ CT study of Zn-Ag batteries and demonstrate the effects of current collector parasitic gassing over long-term storage and cycling. We design performance representative $\textit{in-situ}$ CT cells that can achieve >250 cycles at a high areal capacity of $\mathrm{12.5\;mAh/cm^2}$. Combined with electrochemical experiments, the effects of current collector parasitic gassing are revealed with micro-scale CT (MicroCT). The volume expansion and evolution of ZnO and Zn depletion is quantified with cycling and elevated temperature testing. The experimental insights are then utilized to develop larger form-factor $\mathrm{4\;cm^2}$ cells with electrochemically compatible current collectors. With this, we demonstrate over 325 cycles at a high capacity of $\mathrm{12.5\;mAh/cm^2}$ for a $\mathrm{4\;cm^2}$ form-factor. This work demonstrates that $\textit{in-situ}$ X-ray CT used in long cycle-lifetime studies can be applied to examine a multitude of other battery chemistries to improve their performances.


title: Comment on "Massive electrodynamics and the magnetic monopoles"

abstract: In this paper we correct previous work on magnetic charge plus a photon mass. We show that contrary to previous claims this system has a very simple, closed form solution which is the Dirac string potential multiplied by a exponential decaying part. Interesting features of this solution are discussed, namely, (i) the Dirac string becomes a real feature of the solution, (ii) the breaking of gauge symmetry via the photon mass leads to a breaking of the rotational symmetry of the monopole's magnetic field, (iii) the Dirac quantization condition is potentially altered.


title: Fluctuation-induced quantum friction in nanoscale water flows

abstract: The flow of water in carbon nanochannels has defied understanding thus far, with accumulating experimental evidence for ultra-low friction, exceptionally high water flow rates, and curvature-dependent hydrodynamic slippage. These unique properties have raised considerable interest in carbon-based membranes for desalination, molecular sieving and energy harvesting. However, the mechanism of water-carbon friction remains unknown, with neither current theories, nor classical or ab initio molecular dynamics simulations providing satisfactory rationalisation for its singular behaviour. Here, we develop a quantum theory of the solid-liquid interface, which reveals a new contribution to friction, due to the coupling of charge fluctuations in the liquid to electronic excitations in the solid. We expect that this quantum friction, which is absent in Born-Oppenheimer molecular dynamics, is the dominant friction mechanism for water on carbon-based materials. As a key result, we demonstrate a dramatic difference in quantum friction between the water-graphene and water-graphite interface, due to the coupling of water Debye collective modes with a thermally excited plasmon specific to graphite. This suggests an explanation for the radius-dependent slippage of water in carbon nanotubes, in terms of the nanotubes' electronic excitations. Our findings open the way to quantum engineering of hydrodynamic flows through the confining wall electronic properties.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.soft cond-mat.stat-mech


title: Using reinforcement learning to design an AI assistantfor a satisfying co-op experience

abstract: In this project, we designed an intelligent assistant player for the single-player game Space Invaders with the aim to provide a satisfying co-op experience. The agent behaviour was designed using reinforcement learning techniques and evaluated based on several criteria. We validate the hypothesis that an AI-driven computer player can provide a satisfying co-op experience.


title: General formulas of global characteristic coefficients of Collatz function

abstract: The purpose of this paper is to show three general formulas of three global characteristic coefficients of Collatz function. The Collatz function is defined by the following operation on an arbitrary positive integer if N is odd multiply it by 3 and add 1 then the sum obtained is divided by 2, if N is even divide it by 2. Based on the principle, we define the n-order function denoted by T^n such as the different expressions of that function are results of applying Collatz function n times to a natural numbers which are expressed in well determined forms. Based on these expressions, we can characterize that n-order function by three global characteristic coefficients. In the first, we define these three global coefficients. Secondly, we show that each global characteristic coefficient has a general expression as a function of n.


title: On partitions of integers with restrictions involving squares

abstract: In this paper, we study partitions of positive integers with restrictions involving squares. We mainly establish the following two results (which were conjectured by Sun in 2013): (i) Each positive integer $n$ can be written as $n=x+y+z$ with $x,y,z$ positive integers such that $x^2+y^2+z^2$ is a square, unless $n$ has the form $n=2^{a}3^{b}$ or $2^{a}7$ with $a$ and $b$ nonnegative integers. (ii) Each integer $n>7$ with $n\not=11,14,17$ can be written as $n=x+y+2z$ with $x,y,z$ positive integers such that $x^2+y^2+2z^2$ is a square.


title: Diff-Explainer: Differentiable Convex Optimization for Explainable Multi-hop Inference

abstract: This paper presents Diff-Explainer, the first hybrid framework for explainable multi-hop inference that integrates explicit constraints with neural architectures through differentiable convex optimization. Specifically, Diff-Explainer allows for the fine-tuning of neural representations within a constrained optimization framework to answer and explain multi-hop questions in natural language. To demonstrate the efficacy of the hybrid framework, we combine existing ILP-based solvers for multi-hop Question Answering (QA) with Transformer-based representations. An extensive empirical evaluation on scientific and commonsense QA tasks demonstrates that the integration of explicit constraints in an end-to-end differentiable framework can significantly improve the performance of non-differentiable ILP solvers (8.91% - 13.3%). Moreover, additional analysis reveals that Diff-Explainer is able to achieve strong performance when compared to standalone Transformers and previous multi-hop approaches while still providing structured explanations in support of its predictions.


title: Deep Learning Hamiltonian Monte Carlo

abstract: We generalize the Hamiltonian Monte Carlo algorithm with a stack of neural network layers and evaluate its ability to sample from different topologies in a two dimensional lattice gauge theory. We demonstrate that our model is able to successfully mix between modes of different topologies, significantly reducing the computational cost required to generated independent gauge field configurations. Our implementation is available at https://github.com/saforem2/l2hmc-qcd .

categories: hep-lat cond-mat.stat-mech cs.LG stat.ML


title: Pseudo-Goldstone Dark Matter in $SO(10)$

abstract: We propose a pseudo-Goldstone boson dark matter (pGDM) particle in $SO(10)$ grand unified theory (GUT). Due to its Goldstone nature, this pGDM evades the direct DM detection experiments which, otherwise, severely constrain the parameter space of DM models. In $SO(10)$, the pGDM is embedded as a linear combination of the Standard Model (SM) singlet scalars in ${\bf 16_H}$ and ${\bf 126_H}$ representations. We consider two scenarios for the intermediate route of $SO(10)$ symmetry breaking (SB) to the SM: $SU(5) \times U(1)_X$ and Pati-Salam the $SU(4)_c \times SU(2)_L \times SU(2)_R$ (4-2-2) gauge groups. The vacuum expectation value of ${\bf 126_H}$, which triggers the breaking of $U(1)_X$ and 4-2-2 symmetry in the two scenarios, respectively, determines the pGDM lifetime whose astrophysical lower bound provides one of the most stringent constraints. For the 4-2-2 route to $SO(10)$, the successful SM gauge coupling unification requires the 4-2-2 breaking scale to be ${\cal O} (10^{11})$ GeV, and most of the parameter space is excluded. For the $SU(5) \times U(1)_X$ route, on the other hand, the $U(1)_X$ breaking scale can be significantly higher, and a wide range of the parameter space is allowed. Furthermore, the proton lifetime in the $SU(5)$ case is predicted to be $4.53 \times 10^{34}$ years, which lies well within the sensitivity reach of the Hyper-Kamiokande experiment. We also examine the constraints on the model parameter space from the Large Hadron Collider and the indirect DM search by Fermi-LAT and MAGIC experiments.


title: Compound Arbitrarily Varying Channels

abstract: We propose a communication model, that we call compound arbitrarily varying channels (CAVC), which unifies and generalizes compound channels and arbitrarily varying channels (AVC). A CAVC can be viewed as a noisy channel with a fixed, but unknown, compound-state and an AVC-state which may vary with every channel use. The AVC-state is controlled by an adversary who is aware of the compound-state. We study three problems in this setting: 'communication', 'communication and compound-state identification', and 'communication or compound-state identification'. For these problems, we study conditions for feasibility and capacity under deterministic coding and random coding.


title: Cosmological constraints from unWISE and Planck CMB lensing tomography

abstract: A number of recent, low-redshift, lensing measurements hint at a universe in which the amplitude of lensing is lower than that predicted from the $\Lambda$CDM model fit to the data of the Planck CMB mission. Here we use the auto- and cross-correlation signal of unWISE galaxies and Planck CMB lensing maps to infer cosmological parameters at low redshift. In particular, we consider three unWISE samples (denoted as "blue", "green" and "red") at median redshifts $z \sim 0.6$, $1.1$ and 1.5, which fully cover the Dark Energy dominated era. Our cross-correlation measurements, with combined significance $S/N \sim 80$, are used to infer the amplitude of low-redshift fluctuations, $\sigma_8$; the fraction of matter in the Universe, $\Omega_m$; and the combination $S_8 \equiv \sigma_8 (\Omega_m / 0.3)^{0.5}$ to which these low-redshift lensing measurements are most sensitive. The combination of blue, green and red samples gives a value $S_8=0.784\pm 0.015$, that is fully consistent with other low-redshift lensing measurements and in 2.4$\sigma$ tension with the CMB predictions from Planck. This is noteworthy, because CMB lensing probes the same physics as previous galaxy lensing measurements, but with very different systematics, thus providing an excellent complement to previous measurements.


title: Cosmic Axion Force

abstract: Nambu-Goldstone bosons, or axions, may be ubiquitous. Some of the axions may have small masses and thus serve as mediators of long-range forces. In this paper, we study the force mediated by an extremely light axion, $\phi$, between the visible sector and the dark sector, where dark matter lives. Since nature does not preserve the CP symmetry, the coupling between dark matter and $\phi$ is generically CP-violating. In this case, the induced force is extremely long-range and behaves as an effective magnetic field. If the force acts on electrons or nucleons, the spins of them on Earth precess around a fixed direction towards the galactic center. This provides an experimental opportunity for $\phi$ with mass, $m_\phi$, and decay constant, $f_\phi$, satisfying $m_\phi\lesssim 10^{-25}\,$ eV, $f_\phi\lesssim 10^{14}\,$GeV if the daily modulation of the effective magnetic field signals in magnetometers is measured by using the coherent averaging method. The effective magnetic field induced by an axionic compact object, such as an axion domain wall, is also discussed.


title: Ramification of Volterra-type Rough Paths

abstract: We extend the new approach introduced in arXiv:1912.02064v2 [math.PR] and arXiv:2102.10119v1 [math.PR] for dealing with stochastic Volterra equations using the ideas of Rough Path theory and prove global existence and uniqueness results. The main idea of this approach is simple: Instead of the iterated integrals of a path comprising the data necessary to solve any equation driven by that path, now iterated integral convolutions with the Volterra kernel comprise said data. This leads to the corresponding abstract objects called Volterra-type Rough Paths, as well as the notion of the convolution product, an extension of the natural tensor product used in Rough Path Theory.


title: Brightness Asymmetry of Black Hole Images as a Probe of Observer Inclination

abstract: The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which show a ring-like emission structure with the South side only slightly brighter than the North side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17 deg for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counter examples, we show that low-inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black-hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on.


title: Kernel Two-Sample Tests for Manifold Data

abstract: We present a study of a kernel-based two-sample test statistic related to the Maximum Mean Discrepancy (MMD) in the manifold data setting, assuming that high-dimensional observations are close to a low-dimensional manifold. We characterize the test level and power in relation to the kernel bandwidth, the number of samples, and the intrinsic dimensionality of the manifold. Specifically, when data densities $p$ and $q$ are supported on a $d$-dimensional sub-manifold ${M}$ embedded in an $m$-dimensional space and are H\"older with order $\beta$ (up to 2) on ${M}$, we prove a guarantee of the test power for finite sample size $n$ that exceeds a threshold depending on $d$, $\beta$, and $\Delta_2$ the squared $L^2$-divergence between $p$ and $q$ on the manifold, and with a properly chosen kernel bandwidth $\gamma$. For small density departures, we show that with large $n$ they can be detected by the kernel test when $\Delta_2$ is greater than $n^{- { 2 \beta/( d + 4 \beta ) }}$ up to a certain constant and $\gamma$ scales as $n^{-1/(d+4\beta)}$. The analysis extends to cases where the manifold has a boundary and the data samples contain high-dimensional additive noise. Our results indicate that the kernel two-sample test has no curse-of-dimensionality when the data lie on or near a low-dimensional manifold. We validate our theory and the properties of the kernel test for manifold data through a series of numerical experiments.


title: A stringy perspective on the coincidence problem

abstract: We argue that, for generic string compactifications, dark energy is likely to signal the beginning of the end of our universe, perhaps even through decompactification, with possible implications for the cosmological coincidence problem. Thanks to the scarcity (absence?) of stable de Sitter vacua, dark energy in string theory is assumed to take the form of a quintessence field in slow roll. As it rolls, a tower of heavy states will generically descend, triggering an apocalyptic phase transition in the low energy cosmological dynamics after at most a few hundred Hubble times. As a result, dark energy domination cannot continue indefinitely and there is at least a percentage chance that we find ourselves in the first Hubble epoch. We use a toy model of quintessence coupled to a tower of heavy states to explicitly demonstrate the breakdown in the cosmological dynamics as the tower becomes light. This occurs through a large number of corresponding particles being produced after a certain time, overwhelming quintessence. We also discuss some implications for early universe inflation.


title: Robust output feedback model predictive control using online estimation bounds

abstract: We present a framework to design nonlinear robust output feedback model predictive control (MPC) schemes that ensure constraint satisfaction under noisy output measurements and disturbances. We provide novel estimation methods to bound the magnitude of the estimation error based on: stability properties of the observer; detectability; set-membership estimation; moving horizon estimation (MHE). Robust constraint satisfaction is guaranteed by suitably incorporating these online validated bounds on the estimation error in a homothetic tube based MPC formulation. In addition, we show how the performance can be further improved by combining MHE and MPC in a single optimization problem. The framework is applicable to a general class of detectable and (incrementally) stabilizable nonlinear systems. While standard output feedback MPC schemes use offline computed worst-case bounds on the estimation error, the proposed framework utilizes online validated bounds, thus reducing conservatism and improving performance. We demonstrate the reduced conservatism of the proposed framework using a nonlinear 10-state quadrotor example.


title: Flavour Physics with SuperIso

abstract: We describe SuperIso v4.1, a public code for calculation of flavour observables, in the Standard Model and for New Physics scenarios in a model-independent way as well as for specific beyond the SM models such as two Higgs doublet models (2HDM), minimal supersymmetric Standard Model (MSSM) and next-to-minimal supersymmetric Standard Model (NMSSM).


title: Dark QCD Matters

abstract: We investigate the nightmare scenario of dark sectors that are made of non-abelian gauge theories with fermions, gravitationally coupled to the Standard Model (SM). While testing these scenarios is experimentally challenging, they are strongly motivated by the accidental stability of dark baryons and pions, that explain the cosmological stability of dark matter (DM). We study the production of these sectors which are minimally populated through gravitational freeze-in, leading to a dark sector temperature much lower than the SM, or through inflaton decay, or renormalizable interactions producing warmer DM. Despite having only gravitational couplings with the SM these scenarios turn out to be rather predictive depending roughly on three parameters: the dark sector temperature, the confinement scale and the dark pion mass. In particular, when the initial temperature is comparable to the SM one these scenarios are very constrained by structure formation, $\Delta N_{\rm eff}$ and limits on DM self-interactions. Dark sectors with same temperature or warmer than SM are typically excluded.


title: Magnetohydrodynamic effect on first star formation: prestellar core collapse and protostar formation

abstract: Recent theoretical studies have suggested that a magnetic field may play a crucial role in the first star formation in the universe. However, the influence of the magnetic field on the first star formation has yet to be understood well. In this study, we perform three-dimensional magnetohydrodynamic simulations taking into account all the relevant cooling processes and non-equilibrium chemical reactions up to the protostar density, in order to study the collapse of magnetized primordial gas cores with self-consistent thermal evolution. Our results show that the thermal evolution of the central core is hardly affected by a magnetic field, because magnetic forces do not prevent the contraction along the fields lines. We also find that the magnetic braking extracts the angular momentum from the core and suppresses fragmentation depending on the initial strength of the magnetic field. The angular momentum transport by the magnetic outflows is less effective than that by the magnetic braking because the outflows are launched only in a late phase of the collapse. Our results indicate that the magnetic effects become important for the field strength $B> 10^{-8}(n_{\rm H}/1\ \rm cm^{-3})^{2/3}\ \rm G$, where $n_{\rm H}$ is the number density, during the collapse phase. Finally, we compare our results with simulations using a barotropic approximation and confirm that this approximation is reasonable at least for the collapse phase. Nevertheless, self-consistent treatment of the thermal and chemical processes is essential for extending simulations to the accretion phase, in which radiative feedback by protostars plays a crucial role.


title: DNA Nanotechnology Meets Nanophotonics

abstract: In the past decades, DNA has been intensely studied and exploited in different research areas of nanoscience and nanotechnology. At first glance, DNA-based nanophotonics seems to deviate quite far from the original goal of Nadrian Seeman, the founder of DNA nanotechnology, who hoped to organize biological entities using DNA in high-resolution crystals. As a matter of fact, DNA-based nanophotonics does closely follow his central spirit. That is, apart from being a genetic material for inheritance, DNA is also an ideal material for building molecular devices.


title: Generalising Multilingual Concept-to-Text NLG with Language Agnostic Delexicalisation

abstract: Concept-to-text Natural Language Generation is the task of expressing an input meaning representation in natural language. Previous approaches in this task have been able to generalise to rare or unseen instances by relying on a delexicalisation of the input. However, this often requires that the input appears verbatim in the output text. This poses challenges in multilingual settings, where the task expands to generate the output text in multiple languages given the same input. In this paper, we explore the application of multilingual models in concept-to-text and propose Language Agnostic Delexicalisation, a novel delexicalisation method that uses multilingual pretrained embeddings, and employs a character-level post-editing model to inflect words in their correct form during relexicalisation. Our experiments across five datasets and five languages show that multilingual models outperform monolingual models in concept-to-text and that our framework outperforms previous approaches, especially for low resource languages.


title: Self-collision of a portal wormhole

abstract: We consider the self-collision of portals in classical general relativity. Portals are wormholes supported by a single loop of negative mass cosmic string, and being wormholes, portals have a nontrivial topology. Portals can be constructed so that the curvature is zero everywhere outside the cosmic string, with vanishing ADM mass. The conical singularities of these wormholes can be smoothed, yielding a spatial topology of $S^2 \times S^1$ with a point corresponding to spatial infinity removed. If one attempts to collide the mouths of a smoothed portal to induce self-annihilation, one naively might think that a Euclidean topology is recovered, which would violate the classical no topology change theorems. We consider a particular limit of smoothed portals supported by an anisotropic fluid, and find that while the portal mouths do not experience an acceleration as they are brought close together, a curvature singularity forms in the limit that the separation distance vanishes. We find that in general relativity, the interaction between portal mouths is not primarily gravitational in nature, but depends critically on matter interactions.


title: Lyman Alpha Emitters and the 21cm Power Spectrum as Probes of Density-Ionization Correlation in the Epoch of Reionization

abstract: Due to the large cross section of Ly$\alpha$ photons with hydrogen, Lyman Alpha Emitters (LAEs) are sensitive to the presence of neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR): the period in the Universe's history where neutral hydrogen in the IGM is ionized. The correlation of the ionized regions in the IGM with respect to the underlying intrinsic LAEs has a pronounced effect on the number of observed LAEs and their apparent clustering. As a result, observations of LAEs during the EoR can be used as a probe of the EoR morphology. Here we build on previous works where we parametrize the density-ionisation correlation during the EoR, and study how the observed number density and angular correlation function (ACF) of LAEs depends on this parametrization. Using Subaru measurements of the number density of LAEs and their ACF at z = 6.6, we place constraints on the EoR morphology. We find that measurements of LAEs at z = 6.6 alone cannot distinguish between different density-ionization models at $68\%$ credibility. However, adding information regarding the number density, and ACF, of LAEs at $z = 6.6$ to 21cm power spectrum measurements using the hydrogen Epoch of Reionization Array (HERA) at the midpoint of reionization can rule out uncorrelated and outside-in reionization at $99\%$ credibility.


title: ELVES I: Structures of Dwarf Satellites of MW-like Galaxies; Morphology, Scaling Relations, and Intrinsic Shapes

abstract: The structure of a dwarf galaxy is an important probe into the effects of stellar feedback and environment. Using an unprecedented sample of 223 low-mass satellites from the ongoing Exploration of Local VolumE Satellites (ELVES) Survey, we explore the structures of dwarf satellites in the mass range $10^{5.5}<M_\star<10^{8.5}$M$_\odot$. We survey satellites around $80\%$ of the massive, $M_K<-22.4$ mag, hosts in the Local Volume. Our sample of dwarf satellites is complete to luminosities of $M_V<-9$ mag and surface brightness $\mu_{0,V}<26.5$ mag arcsec$^{-2}$ within at least $\sim200$ projected kpc. We separate the satellites into late- and early-type, finding the mass-size relations are very similar between them, to within $\sim5\%$. This similarity indicates that the quenching and transformation of a late-type dwarf into an early-type involves only very mild size evolution. Considering the distribution of apparent ellipticities, we infer the intrinsic shapes of the early- and late-type samples. Combining with literature samples, we find that both types of dwarfs get thicker at fainter luminosities but early-types are always rounder at fixed luminosity. Finally, we compare the LV satellites with dwarf samples from the cores of the Virgo and Fornax clusters. We find that the cluster satellites show similar scaling relations to the LV early-type dwarfs but are roughly $10\%$ larger at fixed mass, which we interpret as being due to tidal heating in the cluster environments. The dwarf structure results presented here are a useful reference for simulations of dwarf galaxy formation and the transformation of dwarf irregulars into spheroidals.


title: The Weizmann Fast Astronomical Survey Telescope (W-FAST): System Overview

abstract: A relatively unexplored phase space of transients and stellar variability is that of second and sub-second time-scales. We describe a new optical observatory operating in the Negev desert in Israel, with a 55 cm aperture, a field of view of 2.6x2.6 deg (~7deg^2) equipped with a high frame rate, low read noise, CMOS camera. The system can observe at a frame rate of up to 90HZ (full frame), while nominally observations are conducted at 10-25Hz. The data, generated at a rate of over 6Gbits/s at a frame rate of 25Hz, are analyzed in real time. The observatory is fully robotic and capable of autonomously collecting data on a few thousand stars in each field each night. We present the system overview, performance metrics, science objectives, and some first results, e.g., the detection of a high rate of glints from geosynchronous satellites, reported in Nir et al. 2020.


title: CHAOS VI: Direct Abundances in NGC 2403

abstract: We report the direct abundances for the galaxy NGC 2403 as observed by the CHemical Abundances Of Spirals (CHAOS) project. Using the Multi-Object Double Spectrograph on the Large Binocular Telescope, we observe two fields with H II regions that cover an R$_g$/R$_e$ range of 0.18 to 2.31. 32 H II regions contain at least one auroral line detection, and we detect a total of 122 temperature-sensitive auroral lines. Here, for the first time, we use the intrinsic scatter in the T$_e$-T$_e$ diagrams, added in quadrature to the uncertainty on the measured temperature, to determine the uncertainty on an electron temperature inferred for one ionization zone from a measurement in a different ionization zone. We then use all available temperature data within an H II region to obtain a weighted average temperature within each ionization zone. We re-derive the oxygen abundances of all CHAOS galaxies using this new temperature prioritization method, and we find that the gradients are consistent with the results of Berg et al. (2020). For NGC 2403, we measure a direct oxygen abundance gradient of -0.09($\pm$0.03) dex/Re, with an intrinsic dispersion of 0.037($\pm$0.017) dex, and an N/O abundance gradient of -0.17($\pm$0.03) dex/Re with an intrinsic dispersion of 0.060($\pm$0.018) dex. For direct comparison, we use the line intensities from the study of NGC 2403 by Berg et al. (2013), and find their recomputed values for the O/H and N/O gradients are consistent with ours.


title: Optimal observing strategies for velocity-suppressed dark matter annihilation

abstract: Numerous particle models for the cosmological dark matter feature a pair-annihilation rate that scales with powers of the relative velocity between the annihilating particles. As a result, the annihilation rate in the central regions of a dark matter halo can be significantly lower than at the halo's periphery for particular ambient gravitational potentials. While this might be offset by an increasing dark matter pair number density in the inner halo, it raises the question: what angular region for dark matter models with velocity-suppressed annihilation rates optimizes signal-to-noise? Here, we consider simplified background models for galactic and extragalactic targets and demonstrate that the optimal observing strategy varies greatly case-by-case. Generally, a bright central source warrants an annular region of interest, while a flatter background warrants as large as possible an angular region, possibly including the central regions.


title: Hierarchical mergers of stellar-mass black holes and their gravitational-wave signatures

abstract: We review theoretical findings, astrophysical modeling, and current gravitational-wave evidence of hierarchical stellar-mass black-hole mergers. While most of the compact binary mergers detected by LIGO and Virgo are expected to consist of first-generation black holes formed from the collapse of stars, others might instead be of second (or higher) generation, containing the remnants of previous black-hole mergers. Such a subpopulation of hierarchically assembled black holes presents distinctive gravitational-wave signatures, namely higher masses, possibly within the pair-instability mass gap, and dimensionless spins clustered at the characteristic value of $\sim$0.7. In order to produce hierarchical mergers, astrophysical environments need to overcome the relativistic recoils imparted to black-hole merger remnants, a condition which prefers hosts with escape speeds $\gtrsim$ 100 km/s. Promising locations for efficient production of hierarchical mergers include nuclear star clusters and accretion disks surrounding active galactic nuclei, though environments that are less efficient at retaining merger products such as globular clusters may still contribute significantly to the detectable population of repeated mergers. While GW190521 is the single most promising hierarchical-merger candidate to date, constraints coming from large population analyses are becoming increasingly more powerful.

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


title: ELVES II: GCs and Nuclear Star Clusters of Dwarf Galaxies; The Importance of Environment

abstract: We present the properties of the globular clusters (GCs) and nuclear star clusters (NSCs) of low-mass ($10^{5.5}<M_\star<10^{8.5}$ $M_\odot$) early-type satellites of Milky Way-like and small group hosts in the Local Volume (LV) using deep, ground-based data from the ongoing Exploration of Local VolumE Satellites (ELVES) Survey. This sample of 177 dwarfs significantly increases the statistics for studying the star clusters of dwarfs in low-density environments, offering an important comparison to samples from nearby galaxy clusters. The LV dwarfs exhibit significantly lower nucleation fractions at fixed galaxy mass than dwarfs in nearby clusters. The mass of NSCs of LV dwarfs show a similar scaling of $M_{\star,\mathrm{NSC}}\propto M_{\star,\mathrm{gal}}^{0.4}$ as that found in clusters but offset to lower NSC masses. To deal with foreground/background contamination in the GC analysis, we employ both a statistical subtraction and Bayesian approach to infer the average GC system properties from all dwarfs simultaneously. We find that the GC occupation fraction and average abundance are both increasing functions of galaxy stellar mass, and the LV dwarfs show significantly lower average GC abundance at fixed galaxy mass than a comparable sample of Virgo dwarfs analyzed in the same way, demonstrating that GC prevalence also shows an important secondary dependence on the dwarf's environment. This result strengthens the connection between GCs and NSCs in low-mass galaxies. We discuss these observations in the context of modern theories of GC and NSC formation, finding that the environmental dependencies can be well-explained by these models.


title: The Pristine Inner Galaxy Survey (PIGS) III: carbon-enhanced metal-poor stars in the bulge

abstract: The most metal-deficient stars hold important clues about the early build-up and chemical evolution of the Milky Way, and carbon-enhanced metal-poor (CEMP) stars are of special interest. However, little is known about CEMP stars in the Galactic bulge. In this paper, we use the large spectroscopic sample of metal-poor stars from the Pristine Inner Galaxy Survey (PIGS) to identify CEMP stars ([C/Fe] > +0.7) in the bulge region and to derive a CEMP fraction. We identify 96 new CEMP stars in the inner Galaxy, of which 62 are very metal-poor ([Fe/H] < -2.0); this is more than a ten-fold increase compared to the seven previously known bulge CEMP stars. The cumulative fraction of CEMP stars in PIGS is $42^{\,+14\,}_{\,-13} \%$ for stars with [Fe/H] < -3.0, and decreases to $16^{\,+3\,}_{\,-3} \%$ for [Fe/H] < -2.5 and $5.7^{\,+0.6\,}_{\,-0.5} \%$ for [Fe/H] < -2.0. The PIGS inner Galaxy CEMP fraction for [Fe/H] < -3.0 is consistent with the halo fraction found in the literature, but at higher metallicities the PIGS fraction is substantially lower. While this can partly be attributed to a photometric selection bias, such bias is unlikely to fully explain the low CEMP fraction at higher metallicities. Considering the typical carbon excesses and metallicity ranges for halo CEMP-s and CEMP-no stars, our results point to a possible deficiency of both CEMP-s and CEMP-no stars (especially the more metal-rich) in the inner Galaxy. The former is potentially related to a difference in the binary fraction, whereas the latter may be the result of a fast chemical enrichment in the early building blocks of the inner Galaxy.


title: The van der Waals interaction as the starting point for an effective field theory

abstract: We consider the system of three ${}^4$He atoms to assess whether a pure van der Waals potential can be used as a starting point for an effective field theory to describe three-body processes in ultracold atomic systems. Using a long-range van der Waals interaction in combination with short-distance two-body counterterms, we analyze the dependence of two- and three-body observables on the short-distance regulator that is required due to the singular nature of the van der Waals interaction. We benchmark our approach with results obtained with the realistic ${}^4$He-${}^4$He LM2M2 interaction and find good agreement. We furthermore show that in this effective field theory approach no three-body force is required at leading order and that universal van der Waals physics leads to a universal correlation between three-body observables in the absence of an Efimov three-body parameter.


title: Photoevaporation vs. core-powered mass-loss: model comparison with the 3D radius gap

abstract: The EUV/X-ray photoevaporation and core-powered mass-loss models are both capable of reproducing the bimodality in the sizes of small, close-in exoplanets observed by the \textit{Kepler} space mission, often referred to as the "radius gap". However, it is unclear which of these two mechanisms dominates the atmospheric mass-loss which is likely sculpting the radius gap. In this work, we propose a new method of differentiating between the two models, which relies on analysing the radius gap in 3D parameter space. Using models for both mechanisms, and by performing synthetic transit surveys we predict the size and characteristics of a survey capable of discriminating between the two models. We find that a survey of $\gtrsim 5000$ planets, with a wide range in stellar mass and measurement uncertainties at a $\lesssim 5\%$ level is sufficient. Our methodology is robust against moderate false positive contamination of $\lesssim 10\%$. We perform our analysis on two surveys (which do not satisfy our requirements): the \textit{California Kepler Survey} and the \textit{Gaia-Kepler Survey} and find, unsurprisingly, that both data-sets are consistent with either model. We propose a hypothesis test to be performed on future surveys which can robustly ascertain which of the two mechanisms formed the radius gap, provided one dominates over the other.


title: More Exact Results on Chiral Gauge Theories: the Case of the Symmetric Tensor

abstract: We study dynanics of $SU(N-4)$ gauge theories with fermions in rank-2 symmetric tensor and $N$ anti-fundamental representations, by perturbing supersymmetric theories with anomaly-mediated supersymmetry breaking. We find the $SU(N)\times U(1)$ global symmetry is dynamically broken to $SO(N)$ for $N\geq 17$, a different result from conjectures in the literature. For $N<17$, theories flow to infrared fixed points.


title: Entanglement generation in $(1+1)D$ QED scattering processes

abstract: We study real-time meson-meson scattering processes in $(1+1)$-dimensional QED by means of Tensor Networks. We prepare initial meson wave packets with given momentum and position introducing an approximation based on the free fermions model. Then, we compute the dynamics of two initially separated colliding mesons, observing a rich phenomenology as the interaction strength and the initial states are varied in the weak and intermediate coupling regimes. Finally, we consider elastic collisions and measure some scattering amplitudes as well as the entanglement generated by the process. Remarkably, we identify two different regimes for the asymptotic entanglement between the outgoing mesons: it is perturbatively small below a threshold coupling, past which its growth as a function of the coupling abruptly accelerates.

categories: hep-lat cond-mat.stat-mech cond-mat.str-el hep-ph physics.comp-ph quant-ph


title: Long Period Jovian Tilts the Orbits of Two sub-Neptunes Relative to Stellar Spin Axis in Kepler-129

abstract: We present the discovery of Kepler-129 d ($P_{d}=7.2^{+0.4}_{-0.3}$ yr, $m\sin i_{d}=8.3^{+1.1}_{-0.7}\ \rm M_{Jup}$, $ e_{d}=0.15^{+0.07}_{-0.05} $) based on six years of radial velocity (RV) observations from Keck/HIRES. Kepler-129 also hosts two transiting sub-Neptunes: Kepler-129 b ($P_{b}=15.79$ days, $r_{b}=2.40\pm{0.04}\ \rm{R_{\oplus}}$) and Kepler-129 c ($P_{c}=82.20$ days, $r_{c}=2.52\pm{0.07}\ \rm{R_{\oplus}}$) for which we measure masses of $m_{b}<20\ \rm{M_{\oplus}}$ and $m_{c}=43^{+13}_{-12}\ \rm{M_{\oplus}}$. Kepler-129 is an hierarchical system consisting of two tightly-packed inner planets and an external companion whose mass is close to the deuterium burning limit. In such a system, two inner planets precess around the orbital normal of the outer companion, causing their inclinations to oscillate with time. Based on an asteroseismic analysis of Kepler data, we find tentative evidence that Kepler-129 b and c are misaligned with stellar spin axis by $\gtrsim 38$ deg, which could be torqued by Kepler-129 d if it is inclined by $\gtrsim 19$ deg relative to inner planets. Using N-body simulations, we provide additional constraints on the mutual inclination between Kepler-129 d and inner planets by estimating the fraction of time during which two inner planets both transit. The probability that two planets both transit decreases as their misalignment with Kepler-129 d increases. We also find a more massive Kepler-129 c enables the two inner planets to become strongly coupled and more resistant to perturbations from Kepler-129 d. The unusually high mass of Kepler-129 c provides a valuable benchmark for both planetary dynamics and interior structure, since the best-fit mass is consistent with this $\rm{2.5\ R_{\oplus}}$ planet having a rocky surface.


title: Optically driving the radiative Auger transition

abstract: In a radiative Auger process, optical decay is accompanied by simultaneous excitation of other carriers. The radiative Auger process gives rise to weak red-shifted satellite peaks in the optical emission spectrum. These satellite peaks have been observed over a large spectral range: in the X-ray emission of atoms; close to visible frequencies on donors in semiconductors and quantum emitters; and at infrared frequencies as shake-up lines in two-dimensional systems. So far, all the work on the radiative Auger process has focussed on detecting the spontaneous emission. However, the fact that the radiative Auger process leads to photon emission suggests that the transition can also be optically excited. In such an inverted radiative Auger process, excitation would correspond to simultaneous photon absorption and electronic de-excitation. Here, we demonstrate optical driving of the radiative Auger transition on a trion in a semiconductor quantum dot. The radiative Auger and the fundamental transition together form a $\Lambda$-system. On driving both transitions of this $\Lambda$-system simultaneously, we observe a reduction of the fluorescence signal by up to $70\%$. Our results demonstrate a type of optically addressable transition connecting few-body Coulomb interactions to quantum optics. The results open up the possibility of carrying out THz spectroscopy on single quantum emitters with all the benefits of optics: coherent laser sources, efficient and fast single-photon detectors. In analogy to optical control of an electron spin, the $\Lambda$-system between the radiative Auger and the fundamental transitions allows optical control of the emitters' orbital degree of freedom.


title: Testing isomorphism between tuples of subspaces

abstract: Given two tuples of subspaces, can you tell whether the tuples are isomorphic? We develop theory and algorithms to address this fundamental question. We focus on isomorphisms in which the ambient vector space is acted on by either a unitary group or general linear group. If isomorphism also allows permutations of the subspaces, then the problem is at least as hard as graph isomorphism. Otherwise, we provide a variety of polynomial-time algorithms with Matlab implementations to test for isomorphism. Keywords: subspace isomorphism, Grassmannian, Bargmann invariants, $H^\ast$-algebras, quivers, graph isomorphism


title: Smoothness of non-reductive fixed point sets and cohomology of non-reductive GIT quotients

abstract: We establish a method for calculating the Poincar\'e series of moduli spaces constructed as quotients of smooth varieties by suitable non-reductive group actions; examples of such moduli spaces include moduli spaces of unstable vector or Higgs bundles on a smooth projective curve, with a Harder-Narasimhan type of length two. To do so, we first prove a result concerning the smoothness of fixed point sets for suitable non-reductive group actions on smooth varieties. This enables us to prove that quotients of smooth varieties by such non-reductive group actions, which can be constructed using Non-Reductive GIT via a sequence of blow-ups, have at worst finite quotient singularities. We conclude the paper by providing explicit formulae for the Poincar\'e series of these non-reductive GIT quotients.


title: Effect of mass loss due to stellar winds on the formation of supermassive black hole seeds in dense nuclear star clusters

abstract: The observations of high redshifts quasars at $z\gtrsim 6$ have revealed that supermassive black holes (SMBHs) of mass $\sim 10^9\,\mathrm{M_{\odot}}$ were already in place within the first $\sim$ Gyr after the Big Bang. Supermassive stars (SMSs) with masses $10^{3-5}\,\mathrm{M_{\odot}}$ are potential seeds for these observed SMBHs. A possible formation channel of these SMSs is the interplay of gas accretion and runaway stellar collisions inside dense nuclear star clusters (NSCs). However, mass loss due to stellar winds could be an important limitation for the formation of the SMSs and affect the final mass. In this paper, we study the effect of mass loss driven by stellar winds on the formation and evolution of SMSs in dense NSCs using idealised N-body simulations. Considering different accretion scenarios, we have studied the effect of the mass loss rates over a wide range of metallicities $Z_\ast=[.001-1]\mathrm{Z_{\odot}}$ and Eddington factors $f_{\rm Edd}=L_\ast/L_{\mathrm{Edd}}=0.5,0.7,\,\&\, 0.9$. For a high accretion rate of $10^{-4}\,\mathrm{M_{\odot}yr^{-1}}$, SMSs with masses $\gtrsim 10^3\MSun$ could be formed even in a high metallicity environment. For a lower accretion rate of $10^{-5}\,\mathrm{M_{\odot}yr^{-1}}$, SMSs of masses $\sim 10^{3-4}\,\mathrm{M_{\odot}}$ can be formed for all adopted values of $Z_\ast$ and $f_{\rm Edd}$, except for $Z_\ast=\mathrm{Z_{\odot}}$ and $f_{\rm Edd}=0.7$ or 0.9. For Eddington accretion, SMSs of masses $\sim 10^3\,\mathrm{M_{\odot}}$ can be formed in low metallicity environments with $Z_\ast\lesssim 0.01\mathrm{Z_{\odot}}$. The most massive SMSs of masses $\sim 10^5\,\mathrm{M_{\odot}}$ can be formed for Bondi-Hoyle accretion in environments with $Z_\ast \lesssim 0.5\mathrm{Z_{\odot}}$.


title: Gravitational wave signatures of black hole quasi-normal mode instability

abstract: Black hole (BH) spectroscopy has emerged as a powerful approach to extract spacetime information from gravitational wave (GW) observed signals. Yet, quasinormal mode (QNM) spectral instability under high wave-number perturbations has been recently shown to be a common classical general relativistic phenomenon [1]. This requires to assess its impact on the BH QNM spectrum, in particular on BH QNM overtone frequencies. We conclude: i) perturbed BH QNM overtones are indeed potentially observable in the GW waveform, providing information on small-scale environment BH physics, and ii) their detection poses a challenging data analysis problem of singular interest for LISA astrophysics. We adopt a two-fold approach, combining theoretical results from scattering theory with a fine-tuned data analysis on a highly-accurate numerical GW ringdown signal. The former introduces a set of effective parameters (partially lying on a BH Weyl law) to characterise QNM instability physics. The latter provides a proof-of-principle demonstrating that the QNM spectral instability is indeed accessible in the time-domain GW waveform, though certainly requiring large signal-to-noise ratios. Particular attention is devoted to discuss the patterns of isospectrality loss under QNM instability, since the disentanglement between axial and polar GW parities may already occur within the near-future detection range.


title: New Numerical Interface Scheme for the Kurganov-Tadmor second-order Method

abstract: In this paper, we develop a numerical scheme to handle interfaces across computational domains in multi-block schemes for the approximation of systems of conservation laws. We are interested in transmitting shock discontinuities without lowering the overall precision of the method. We want to accomplish this without using information from interior points of adjacent grids, that is, sharing only information from boundary points of those grids. To achieve this, we choose to work with the second-order Kurganov-Tadmor (KT) method at interior points, relaxing it to first order at interfaces. This allows us to keep second-order overall accuracy (in the relevant norm) and at the same time preserve the TVD property of the original scheme. After developing the method we performed several standard one and two-dimensional tests. Among them, we used the one-dimensional advection and Burgers equations to verify the second-order convergence of the method. We also tested the two-dimensional Euler equations with an implosion and a Gresho vortex\cite{liska2003}. In particular, in the two-dimensional implosion test we can see that regardless of the orientation of shocks with respect to the interface, they travel across them without appreciable deformation both in amplitude and front direction.


title: Physics reach of ${D}_{(s)}\to \pi(K) \ell\ell$ and other charming null test opportunities

abstract: We discuss possibilities to test physics beyond the Standard Model in $\vert\Delta c\vert=\vert\Delta u\vert= 1$ semileptonic, hadronic and missing energy decay modes. Clean null test observables such as angular observables, CP-asymmetries and lepton universality tests are presented and model-independent correlations as well as details within flavorful, anomaly-free $Z^\prime$ models are worked out.


title: A Bayesian Gaussian Process for Estimating a Causal Exposure Response Curve in Environmental Epidemiology

abstract: Motivated by environmental policy questions, we address the challenges of estimation, change point detection, and uncertainty quantification of a causal exposure-response function (CERF). Under a potential outcome framework, the CERF describes the relationship between a continuously varying exposure (or treatment) and its causal effect on an outcome. We propose a new Bayesian approach that relies on a Gaussian process (GP) model to estimate the CERF nonparametrically. To achieve the desired separation of design and analysis phases, we parametrize the covariance (kernel) function of the GP to mimic matching via a Generalized Propensity Score (GPS). The hyper-parameters as well as the form of the kernel function of the GP are chosen to optimize covariate balance. Our approach achieves automatic uncertainty evaluation of the CERF with high computational efficiency, and enables change point detection through inference on derivatives of the CERF. We provide theoretical results showing the correspondence between our Bayesian GP framework and traditional approaches in causal inference for estimating causal effects of a continuous exposure. We apply the methods to 520,711 ZIP-code-level observations to estimate the causal effect of long-term exposures to PM2.5, ozone, and NO2 on all-cause mortality among Medicare enrollees in the US. A computationally efficient implementation of the proposed GP models is provided in the GPCERF R package, which is available on CRAN.


title: Lipschitz Chain Approximation of Metric Integral Currents

abstract: Every integral current in a locally compact metric space $X$ can be approximated by a Lipschitz chain with respect to the normal mass, provided that Lipschitz maps into $X$ can be extended slightly.


title: CASTing a Net: Supporting Teachers with Search Technology

abstract: Past and current research has typically focused on ensuring that search technology for the classroom serves children. In this paper, we argue for the need to broaden the research focus to include teachers and how search technology can aid them. In particular, we share how furnishing a behind-the-scenes portal for teachers can empower them by providing a window into the spelling, writing, and concept connection skills of their students.


title: An extension theory for partial groups

abstract: Partial groups are a natural generalization of discrete groups recently introduced by Chermak in connection with the theory of fusion systems. In this paper we develop an extension theory for partial groups based in the classical theory of fibre bundles of simplicial sets.


title: Duplex Sequence-to-Sequence Learning for Reversible Machine Translation

abstract: Sequence-to-sequence learning naturally has two directions. How to effectively utilize supervision signals from both directions? Existing approaches either require two separate models, or a multitask-learned model but with inferior performance. In this paper, we propose REDER (Reversible Duplex Transformer), a parameter-efficient model and apply it to machine translation. Either end of REDER can simultaneously input and output a distinct language. Thus REDER enables reversible machine translation by simply flipping the input and output ends. Experiments verify that REDER achieves the first success of reversible machine translation, which helps outperform its multitask-trained baselines by up to 1.3 BLEU.


title: Anisotropic Frictional Response of Texture Induced Strained Graphene

abstract: Friction-induced energy dissipation impedes the performance of nanoscale devices during their relative motion. Nevertheless, an ingeniously designed structure which utilizes graphene topping can tune the friction force signal by inducing local strain. The present work reports capping of graphene over Si grooved surfaces of different pitch lengths from sub-nanoscale (P=40 nm) to a quarter of a micron (P= 250 nm). The variation in the pitch lengths induces different strains in graphene revealed by scanning probe techniques, Raman spectroscopy and molecular dynamics (MD) simulation. The asymmetric straining of C-C bonds over the groove architecture is exploited through friction force microscopy in different directions of orthogonal and parallel to groove axis. The presence of graphene lubricates the textured surface by a factor of 10 and periodically dissipated friction force, which was found to be stochastic over the bare surface. For the first time, we presented transformation of the lubrication into an ultra-low friction force by a factor of 20 over the crest scanning parallel to the groove axis. Such anisotropy is found to be insignificant at the bare textured system, clearly demonstrating the strain-dependent regulation of friction force. Our results are applicable for graphene, and other 2D materials covered corrugated structures with movable components such as NEMS, nanoscale gears and robotics.


title: The landscape of promising non-supersymmetric string models

abstract: Leptoquarks extending the Standard Model (SM) are attracting an increasing attention in the recent literature. Hence, the identification of 4D SM-like models and the classification of allowed leptoquarks from strings is an important step in the study of string phenomenology. We perform the most extensive search for SM-like models from the non-supersymmetric heterotic string $\mathrm{SO}(16)\times\mathrm{SO}(16)$, resulting in more than 170,000 inequivalent promising string models from 138 Abelian toroidal orbifolds. We explore the 4D massless particle spectra of these models in order to identify all exotics beside the three generations of quarks and leptons. Hereby, we learn which leptoquark can be realized in this string setup. Moreover, we analyze the number of SM Higgs doublets which is generically larger than one. Then, we identify SM-like models with a minimal particle content. These so-called "almost SM" models appear most frequently in the orbifold geometries $\mathbb Z_2\times\mathbb Z_4$ (2,4) and (1,6). Finally, we apply machine learning to our dataset in order to predict the orbifold geometry where a given particle spectrum can be found most likely.


title: Impact of DER Communication Delay in AGC: Cyber-Physical Dynamic Simulation

abstract: Distributed energy resource (DER) frequency regulations are promising technologies for future grid operation. Unlike conventional generators, DERs might require open communication networks to exchange signals with control centers, possibly through DER aggregators; therefore, the impacts of the communication variations on the system stability need to be investigated. This paper develops a cyber-physical dynamic simulation model based on the Hierarchical Engine for Large-Scale Co-Simulation (HELICS) to evaluate the impact of the communication variations, such as delays in DER frequency regulations. The feasible delay range can be obtained under different parameter settings. The results show that the risk of instability generally increases with the communication delay.


title: Necessary and Sufficient Girth Conditions for Tanner Graphs of Quasi-Cyclic LDPC Codes

abstract: This paper revisits the connection between the girth of a protograph-based LDPC code given by a parity-check matrix and the properties of powers of the product between the matrix and its transpose in order to obtain the necessary and sufficient conditions for a code to have given girth between 6 and 12, and to show how these conditions can be incorporated into simple algorithms to construct codes of that girth. To this end, we highlight the role that certain submatrices that appear in these products have in the construction of codes of desired girth. In particular, we show that imposing girth conditions on a parity-check matrix is equivalent to imposing conditions on a square submatrix obtained from it and we show how this equivalence is particularly strong for a protograph based parity-check matrix of variable node degree 2, where the cycles in its Tanner graph correspond one-to-one to the cycles in the Tanner graph of a square submatrix obtained by adding the permutation matrices (or products of these) in the composition of the parity-check matrix. We end the paper with exemplary constructions of codes with various girths and computer simulations. Although, we mostly assume the case of fully connected protographs of variable node degree 2 and 3, the results can be used for any parity-check matrix/protograph-based Tanner graph.


title: Conditional a posteriori error bounds for high order DG time stepping approximations of semilinear heat models with blow-up

abstract: This work is concerned with the development of an adaptive numerical method for semilinear heat flow models featuring a general (possibly) nonlinear reaction term that may cause the solution to blow up in finite time. The fully discrete scheme consists of a high order discontinuous Galerkin (dG) time stepping method and a conforming finite element discretisation (cG) in space. The proposed adaptive procedure is based on rigorously devised conditional a posteriori error bounds in the $L^{\infty}(L^{\infty})$ norm. Numerical experiments complement the theoretical results.


title: Estimating Parkinsonism Severity in Natural Gait Videos of Older Adults with Dementia

abstract: Drug-induced parkinsonism affects many older adults with dementia, often causing gait disturbances. New advances in vision-based human pose-estimation have opened possibilities for frequent and unobtrusive analysis of gait in residential settings. This work leverages novel spatial-temporal graph convolutional network (ST-GCN) architectures and training procedures to predict clinical scores of parkinsonism in gait from video of individuals with dementia. We propose a two-stage training approach consisting of a self-supervised pretraining stage that encourages the ST-GCN model to learn about gait patterns before predicting clinical scores in the finetuning stage. The proposed ST-GCN models are evaluated on joint trajectories extracted from video and are compared against traditional (ordinal, linear, random forest) regression models and temporal convolutional network baselines. Three 2D human pose-estimation libraries (OpenPose, Detectron, AlphaPose) and the Microsoft Kinect (2D and 3D) are used to extract joint trajectories of 4787 natural walking bouts from 53 older adults with dementia. A subset of 399 walks from 14 participants is annotated with scores of parkinsonism severity on the gait criteria of the Unified Parkinson's Disease Rating Scale (UPDRS) and the Simpson-Angus Scale (SAS). Our results demonstrate that ST-GCN models operating on 3D joint trajectories extracted from the Kinect consistently outperform all other models and feature sets. Prediction of parkinsonism scores in natural walking bouts of unseen participants remains a challenging task, with the best models achieving macro-averaged F1-scores of 0.53 +/- 0.03 and 0.40 +/- 0.02 for UPDRS-gait and SAS-gait, respectively. Pre-trained model and demo code for this work is available: https://github.com/TaatiTeam/stgcn_parkinsonism_prediction.


title: Front speed and pattern selection of a propagating chemical front in an active fluid

abstract: Spontaneous pattern formation in living systems is driven by reaction-diffusion chemistry and active mechanics. The feedback between chemical and mechanical forces is often essential to robust pattern formation, yet it remains poorly understood in general. In this analytical and numerical paper, we study an experimentally-motivated minimal model of coupling between reaction-diffusion and active matter: a propagating front of an autocatalytic and stress-generating species. In the absence of activity, the front is described by the the well-studied KPP equation. We find that front propagation is maintained even in active systems, with crucial differences: an extensile stress increases the front speed beyond a critical magnitude of the stress, while a contractile stress has no effect on the front speed but can generate a periodic instability in the high-concentration region behind the front. We expect our results to be useful in interpreting pattern formation in active systems with mechano-chemical coupling in vivo and in vitro.


title: Perron value and moment of rooted trees

abstract: The Perron value $\rho(T)$ of a rooted tree $T$ has a central role in the study of the algebraic connectivity and characteristic set, and it can be considered a weight of spectral nature for $T$. A different, combinatorial weight notion for $T$ - the moment $\mu(T)$ - emerges from the analysis of Kemeny's constant in the context of random walks on graphs. In the present work, we compare these two weight concepts showing that $\mu(T)$ is "almost" an upper bound for $\rho(T)$ and the ratio $\mu(T)/\rho(T)$ is unbounded but at most linear in the order of $T$. To achieve these primary goals, we introduce two new objects associated with $T$ - the Perron entropy and the neckbottle matrix - and we investigate how different operations on the set of rooted trees affect the Perron value and the moment.


title: Proof of CMB-driven X-ray brightening of high-z radio galaxies

abstract: We present a definitive assessment of the role of Inverse Compton scattering of Cosmic Microwave Background photons (IC/CMB) in the context of radio galaxies. Owing to the steep increase of the CMB radiation energy density, IC/CMB is supposed to become progressively more important with respect to radio synchrotron cooling as the redshift increases. For typical energies at play, this process will up-scatter the CMB photons into the X-ray band, and is thus expected to yield a redshift-dependent, concurrent X-ray brightening and radio dimming of the jet-powered structures. Here we show how a conclusive proof of this effect hinges on high-resolution imaging data in which the extended lobes can be distinguished from the compact hot spots where synchrotron-self-Compton dominates the X-ray emission regardless of redshift. We analyze Chandra and Very Large Array data of 11 radio galaxies between 1.3<z<4.3, and demonstrate that the emission from their lobes is fully consistent with the expectations from IC/CMB in equipartition. Once the dependence on size and radio luminosity are properly accounted for, the measured lobe X-ray luminosities bear the characteristic (1+z)**4 proportionality expected of a CMB seed radiation field. Whereas this effect can effectively quench the (rest-frame) GHz radio emission from z>3 radio galaxies below <1 mJy, IC/CMB alone can not be responsible for a deficit in high-z, radio-loud AGN if--as we argue--such AGN typically have bright, compact hot spots.


title: zeus: A Python implementation of Ensemble Slice Sampling for efficient Bayesian parameter inference

abstract: We introduce zeus, a well-tested Python implementation of the Ensemble Slice Sampling (ESS) method for Bayesian parameter inference. ESS is a novel Markov chain Monte Carlo (MCMC) algorithm specifically designed to tackle the computational challenges posed by modern astronomical and cosmological analyses. In particular, the method requires only minimal hand--tuning of 1-2 hyper-parameters that are often trivial to set; its performance is insensitive to linear correlations and it can scale up to 1000s of CPUs without any extra effort. Furthermore, its locally adaptive nature allows to sample efficiently even when strong non-linear correlations are present. Lastly, the method achieves a high performance even in strongly multimodal distributions in high dimensions. Compared to emcee, a popular MCMC sampler, zeus performs 9 and 29 times better in a cosmological and an exoplanet application respectively.

categories: astro-ph.IM astro-ph.CO astro-ph.EP physics.comp-ph