title: A Generalization of Classical Formulas in Numerical Integration and Series Convergence Acceleration

abstract: Summation formulas, such as the Euler-Maclaurin expansion or Gregory's quadrature, have found many applications in mathematics, ranging from accelerating series, to evaluating fractional sums and analyzing asymptotics, among others. We show that these summation formulas actually arise as particular instances of a single series expansion, including Euler's method for alternating series. This new summation formula gives rise to a family of polynomials, which contain both the Bernoulli and Gregory numbers in their coefficients. We prove some properties of those polynomials, such as recurrence identities and symmetries. Lastly, we present one case study, which illustrates one potential application of the new expansion for finite impulse response (FIR) filters.


title: Representative Functional Connectivity Learning for Multiple Clinical groups in Alzheimer's Disease

abstract: Mild cognitive impairment (MCI) is a high-risk dementia condition which progresses to probable Alzheimer's disease (AD) at approximately 10% to 15% per year. Characterization of group-level differences between two subtypes of MCI - stable MCI (sMCI) and progressive MCI (pMCI) is the key step to understand the mechanisms of MCI progression and enable possible delay of transition from MCI to AD. Functional connectivity (FC) is considered as a promising way to study MCI progression since which may show alterations even in preclinical stages and provide substrates for AD progression. However, the representative FC patterns during AD development for different clinical groups, especially for sMCI and pMCI, have been understudied. In this work, we integrated autoencoder and multi-class classification into a single deep model and successfully learned a set of clinical group related feature vectors. Specifically, we trained two non-linear mappings which realized the mutual transformations between original FC space and the feature space. By mapping the learned clinical group related feature vectors to the original FC space, representative FCs were constructed for each group. Moreover, based on these feature vectors, our model achieves a high classification accuracy - 68% for multi-class classification (NC vs SMC vs sMCI vs pMCI vs AD).


title: Inference with generalizable classifier predictions

abstract: This paper addresses the problem of making statistical inference about a population that can only be identified through classifier predictions. The problem is motivated by scientific studies in which human labels of a population are replaced by a classifier. For downstream analysis of the population based on classifier predictions to be sound, the predictions must generalize equally across experimental conditions. In this paper, we formalize the task of statistical inference using classifier predictions, and propose bootstrap procedures to allow inference with a generalizable classifier. We demonstrate the performance of our methods through extensive simulations and a case study with live cell imaging data.


title: Technoeconomic Analysis of Thermal Energy Grid Storage Using Graphite and Tin

abstract: Energy storage is needed to enable dispatchable renewable energy supply and thereby full decarbonization of the grid. However, this can only occur with drastic cost reductions compared to current battery technology, with predicted targets for the cost per unit energy (CPE) below $20/kWh. Notably, for full decarbonization, long duration storage up to 100 hrs will be needed at such low costs, and prior analyses have shown that in such high renewable penetration scenarios, CPE is more critical than other parameters such as round trip efficiency or cost per unit power when comparing the costs of different technologies. Here, we introduce an electricity storage concept that stores electricity as sensible heat in graphite storage blocks and uses multi-junction thermophotovoltaics (TPV) as a heat engine to convert it back to electricity on demand. This design is an outgrowth of the system proposed by Amy et al. in 2019, which has been modified here to use a solid graphite medium and molten tin as a heat transfer fluid rather than silicon as both. The reason for this is two-fold: (1) the CPE of graphite is almost 10X lower than that of silicon, which derives from the lower cost per unit mass (i.e., $0.5/kg vs. $1.5/kg) and the higher heat capacity per unit mass (2000 J/kg-K vs. 950 J/kg-K); and (2) the melting point tin and solubility of tin in graphite are much lower than that of silicon, which lessens the number of issues that have to overcome along the research and development pathway. The usage of graphite also eliminates the need for a second tank, but the main disadvantage of using a solid medium is that one cannot easily provide a steady discharge rate, as the power output from the storage will change with time, as the graphite cools during discharge. Thus, the objective of this work is to examine how these changes in the system design effect the overall technoeconomics.

categories: physics.app-ph physics.flu-dyn physics.pop-ph physics.soc-ph


title: On numerical aspects of parameter identification for the Landau-Lifshitz-Gilbert equation in Magnetic Particle Imaging

abstract: The Landau-Lifshitz-Gilbert equation yields a mathematical model to describe the evolution of the magnetization of a magnetic material, particularly in response to an external applied magnetic field. It allows one to take into account various physical effects, such as the exchange within the magnetic material itself. In particular, the Landau-Lifshitz-Gilbert equation encodes relaxation effects, i.e., it describes the time-delayed alignment of the magnetization field with an external magnetic field. These relaxation effects are an important aspect in magnetic particle imaging, particularly in the calibration process. In this article, we address the data-driven modeling of the system function in magnetic particle imaging, where the Landau-Lifshitz-Gilbert equation serves as the basic tool to include relaxation effects in the model. We formulate the respective parameter identification problem both in the all-at-once and the reduced setting, present reconstruction algorithms that yield a regularized solution and discuss numerical experiments. Apart from that, we propose a practical numerical solver to the nonlinear Landau-Lifshitz-Gilbert equation, not via the classical finite element method, but through solving only linear PDEs in an inverse problem framework.


title: Bayesian hierarchical nonlinear modelling of intra-abdominal volume during pneumoperitoneum for laparoscopic surgery

abstract: Laparoscopy is an operation carried out in the abdomen or pelvis through small incisions with external visual control by a camera. This technique needs the abdomen to be insufflated with carbon dioxide to obtain a working space for surgical instruments' manipulation. Identifying the critical point at which insufflation should be limited is crucial to maximizing surgical working space and minimizing injurious effects. Bayesian nonlinear growth mixed-effects models are applied to data coming from a repeated measures design. This study allows to assess the relationship between the insufflation pressure and the intra--abdominal volume.


title: Toward Automatic Interpretation of 3D Plots

abstract: This paper explores the challenge of teaching a machine how to reverse-engineer the grid-marked surfaces used to represent data in 3D surface plots of two-variable functions. These are common in scientific and economic publications; and humans can often interpret them with ease, quickly gleaning general shape and curvature information from the simple collection of curves. While machines have no such visual intuition, they do have the potential to accurately extract the more detailed quantitative data that guided the surface's construction. We approach this problem by synthesizing a new dataset of 3D grid-marked surfaces (SurfaceGrid) and training a deep neural net to estimate their shape. Our algorithm successfully recovers shape information from synthetic 3D surface plots that have had axes and shading information removed, been rendered with a variety of grid types, and viewed from a range of viewpoints.


title: A multiresolution adaptive wavelet method for nonlinear partial differential equations

abstract: The multiscale complexity of modern problems in computational science and engineering can prohibit the use of traditional numerical methods in multi-dimensional simulations. Therefore, novel algorithms are required in these situations to solve partial differential equations (PDEs) with features evolving on a wide range of spatial and temporal scales. To meet these challenges, we present a multiresolution wavelet algorithm to solve PDEs with significant data compression and explicit error control. We discretize in space by projecting fields and spatial derivative operators onto wavelet basis functions. We provide error estimates for the wavelet representation of fields and their derivatives. Then, our estimates are used to construct a sparse multiresolution discretization which guarantees the prescribed accuracy. Additionally, we embed a predictor-corrector procedure within the temporal integration to dynamically adapt the computational grid and maintain the accuracy of the solution of the PDE as it evolves. We present examples to highlight the accuracy and adaptivity of our approach.


title: Confinement and Renormalization Group Equations in String-inspired Non-local Gauge Theories

abstract: As an extension of the weak perturbation theory obtained in recent analysis on infinite-derivative non-local non-Abelian gauge theories motivated from p-adic string field theory, and postulated as direction of UV-completion in 4-D Quantum Field Theory (QFT), here we investigate the confinement conditions and $\beta-$function in the strong coupling regime. We extend the confinement criterion, previously obtained by Kugo and Ojima for the local theory based on the Becchi-Rouet-Stora-Tyutin (BRST) invariance, to the non-local theory, by using a set of exact solutions of the corresponding local theory. We show that the infinite-derivatives which are active in the UV provides finite contributions also in the infrared (IR) limit and provide a proof of confinement, granted by the absence of the Landau pole. The main difference with the local case is that the IR fixed point is moved to infinity. We also show that in the limit of the energy scale of non-locality $M \rightarrow \infty$ we reproduce the local theory results and see how asymptotic freedom is properly recovered.


title: Hierarchically Regularized Deep Forecasting

abstract: Hierarchical forecasting is a key problem in many practical multivariate forecasting applications - the goal is to simultaneously predict a large number of correlated time series that are arranged in a pre-specified aggregation hierarchy. The main challenge is to exploit the hierarchical correlations to simultaneously obtain good prediction accuracy for time series at different levels of the hierarchy. In this paper, we propose a new approach for hierarchical forecasting which consists of two components. First, decomposing the time series along a global set of basis time series and modeling hierarchical constraints using the coefficients of the basis decomposition. And second, using a linear autoregressive model with coefficients that vary with time. Unlike past methods, our approach is scalable (inference for a specific time series only needs access to its own history) while also modeling the hierarchical structure via (approximate) coherence constraints among the time series forecasts. We experiment on several public datasets and demonstrate significantly improved overall performance on forecasts at different levels of the hierarchy, compared to existing state-of-the-art hierarchical models.


title: Improved Transformer for High-Resolution GANs

abstract: Attention-based models, exemplified by the Transformer, can effectively model long range dependency, but suffer from the quadratic complexity of self-attention operation, making them difficult to be adopted for high-resolution image generation based on Generative Adversarial Networks (GANs). In this paper, we introduce two key ingredients to Transformer to address this challenge. First, in low-resolution stages of the generative process, standard global self-attention is replaced with the proposed multi-axis blocked self-attention which allows efficient mixing of local and global attention. Second, in high-resolution stages, we drop self-attention while only keeping multi-layer perceptrons reminiscent of the implicit neural function. To further improve the performance, we introduce an additional self-modulation component based on cross-attention. The resulting model, denoted as HiT, has a nearly linear computational complexity with respect to the image size and thus directly scales to synthesizing high definition images. We show in the experiments that the proposed HiT achieves state-of-the-art FID scores of 30.83 and 2.95 on unconditional ImageNet $128 \times 128$ and FFHQ $256 \times 256$, respectively, with a reasonable throughput. We believe the proposed HiT is an important milestone for generators in GANs which are completely free of convolutions. Our code is made publicly available at https://github.com/google-research/hit-gan


title: Solutions of the bosonic master-field equation from a supersymmetric matrix model

abstract: It has been argued that the bosonic large-$N$ master field of the IIB matrix model can give rise to an emergent classical spacetime. In a recent paper, we have obtained solutions of a simplified bosonic master-field equation from a related matrix model. In this simplified equation, the effects of dynamic fermions were removed. We now consider the full bosonic master-field equation from a related supersymmetric matrix model for dimensionality $D=3$ and matrix size $N=3$. In this last equation, the effects of dynamic fermions are included. With an explicit realization of the random constants entering this algebraic equation, we establish the existence of nontrivial solutions. The small matrix size, however, does not allow us to make a definitive statement as to the appearance of a diagonal/band-diagonal structure in the obtained matrices.


title: Nonperturbative Renormalization of the Quark Chromoelectric Dipole Moment with the Gradient Flow:Power Divergences

abstract: The CP-violating quark chromoelectric dipole moment (qCEDM) operator, contributing to the electric dipole moment (EDM), mixes under renormalization and particularly on the lattice with the pseudoscalar density. The mixing coefficient is power-divergent with the inverse lattice spacing squared, $1/a^2$, regardless of the lattice action used. We use the gradient flow to define a multiplicatively renormalized qCEDM operator and study its behavior at small flow time. We determine nonperturbatively the linearly divergent coefficient with the flow time, $1/t$, up to subleading logarithmic corrections, and compare it with the 1-loop perturbative expansion in the bare and renormalized strong coupling. We also discuss the O($a$) improvement of the qCEDM defined at positive flow time.


title: Fast algorithm for quantum polar decomposition, pretty-good measurements, and the Procrustes problem

abstract: The polar decomposition of a matrix is a key element in the quantum linear algebra toolbox. We show that the problem of quantum polar decomposition, recently studied in Lloyd et al. [LBP+20], has a simple and concise implementation via the quantum singular value transform (QSVT). We focus on the applications to pretty-good measurements, a close-to-optimal measurement to distinguish quantum states, and the quantum Procrustes problem, the task of learning an optimal unitary mapping between given `input' and `output' quantum states. By transforming the state-preparation unitaries into a block-encoding, a pre-requisite for QSVT, we develop algorithms for these problems whose gate complexity exhibits a polynomial advantage in the size and condition number of the input compared to alternative approaches for the same problem settings [LBP+20, GLMQW20]. For these applications of the polar decomposition, we also obtain an exponential speedup in precision compared to [LBP+20], as the block-encodings remove the need for the costly density matrix exponentiation step. We contribute a rigorous analysis of the approach of [LBP+20].


title: Variational Causal Networks: Approximate Bayesian Inference over Causal Structures

abstract: Learning the causal structure that underlies data is a crucial step towards robust real-world decision making. The majority of existing work in causal inference focuses on determining a single directed acyclic graph (DAG) or a Markov equivalence class thereof. However, a crucial aspect to acting intelligently upon the knowledge about causal structure which has been inferred from finite data demands reasoning about its uncertainty. For instance, planning interventions to find out more about the causal mechanisms that govern our data requires quantifying epistemic uncertainty over DAGs. While Bayesian causal inference allows to do so, the posterior over DAGs becomes intractable even for a small number of variables. Aiming to overcome this issue, we propose a form of variational inference over the graphs of Structural Causal Models (SCMs). To this end, we introduce a parametric variational family modelled by an autoregressive distribution over the space of discrete DAGs. Its number of parameters does not grow exponentially with the number of variables and can be tractably learned by maximising an Evidence Lower Bound (ELBO). In our experiments, we demonstrate that the proposed variational posterior is able to provide a good approximation of the true posterior.


title: Meta Two-Sample Testing: Learning Kernels for Testing with Limited Data

abstract: Modern kernel-based two-sample tests have shown great success in distinguishing complex, high-dimensional distributions with appropriate learned kernels. Previous work has demonstrated that this kernel learning procedure succeeds, assuming a considerable number of observed samples from each distribution. In realistic scenarios with very limited numbers of data samples, however, it can be challenging to identify a kernel powerful enough to distinguish complex distributions. We address this issue by introducing the problem of meta two-sample testing (M2ST), which aims to exploit (abundant) auxiliary data on related tasks to find an algorithm that can quickly identify a powerful test on new target tasks. We propose two specific algorithms for this task: a generic scheme which improves over baselines and a more tailored approach which performs even better. We provide both theoretical justification and empirical evidence that our proposed meta-testing schemes out-perform learning kernel-based tests directly from scarce observations, and identify when such schemes will be successful.


title: On a class of divergence form linear parabolic equations with degenerate coefficients

abstract: We study a class of linear parabolic equations in divergence form with degenerate coefficients on the upper half space. Specifically, the equations are considered in $(-\infty, T) \times \mathbb{R}^d_+$, where $\mathbb{R}^d_+ = \{x \in \mathbb{R}^d\,:\, x_d>0\}$ and $T\in {(-\infty, \infty]}$ is given, and the diffusion matrices are the product of $x_d$ and bounded uniformly elliptic matrices, which are degenerate at $\{x_d=0\}$. As such, our class of equations resembles well the corresponding class of degenerate viscous Hamilton-Jacobi equations. We obtain wellposedness results and regularity type estimates in some appropriate weighted Sobolev spaces for the solutions.


title: Screening vs. gevolution: in chase of a perfect cosmological simulation code

abstract: We compare two competing relativistic approaches to the N-body simulation of the Universe large-scale structure. To this end, employing the corresponding alternative computer codes ("gevolution" and "screening"), we conduct a series of cosmological simulations in boxes of different sizes and calculate the power spectra of the scalar perturbation $\Phi$, the frame-dragging vector potential ${\bf B}$ and the difference between scalar modes $\chi=\Phi-\Psi$. We demonstrate that the corresponding power spectra are in very good agreement between the compared schemes. For example, the relative difference of the power spectra for $\Phi$ is 0.04% maximum. Since the perturbed Einstein equations have much simpler form in the screening approach, the simulation with this code consumes less computational time, saving almost 40% of CPU hours.


title: Grasping force estimation using state-space model and Kalman filter

abstract: The grip force required to handle an object depends on the object's mass and the friction coefficient of its surface. The control of grip force in myoelectric prosthesis is crucial for handling objects adequately. In the current paper we propose a new method for improving the proportional and continuous grasping force estimation to improve control systems for myoelectric prosthesis based on surface electromyography (sEMG) recordings. For this purpose, we develop an approach based on multivariable system identification in the state-space (SS) and continuous force estimation with Kalman Filter (KF). The sEMG recordings of ten healthy individuals performing a grip task were used as data set for model identification. The root mean square (RMS), the mean absolute value (MAV), and the waveform length (WL) extracted from the sEMG signals were used at the model's input while the measured grasping force was the output. The performance of the method was evaluated with the normalized root-mean-squared-error (NRMSE) and the square of the Pearson's correlation coefficient ($ R^2 $). We found the $R^2$ and NRMSE values were 0.92$ \pm $ 0.0319 and 0.723$ \pm $ 0.0563, respectively. The performance of the proposed technique was superior to the results obtained with other regression models, such as the recurrent nonlinear autoregressive exogenous (NARX)-based neural network, the multi-layer perceptron (MLP) network and the linear discriminant analysis (LDA) with a quadratic polynomial fitting (QPF). The results confirm that the method is adequate for real-time applications with myoelectric hand prostheses.


title: Superconductivity in rhombohedral trilayer graphene

abstract: We report the observation of superconductivity in rhombohedral trilayer graphene electrostatically doped with holes. Superconductivity occurs in two distinct regions within the space of gate-tuned charge carrier density and applied electric displacement field, which we denote SC1 and SC2. The high sample quality allows for detailed mapping of the normal state Fermi surfaces by quantum oscillations, which reveal that in both cases superconductivity arises from a normal state described by an annular Fermi sea that is proximal to an isospin symmetry breaking transition where the Fermi surface degeneracy changes. The upper out-of-plane critical field $B_{C\perp}\approx 10 \mathrm{mT}$ for SC1 and $1\mathrm{mT}$ for SC2, implying coherence lengths $\xi$ of 200nm and 600nm, respectively. The simultaneous observation of transverse magnetic electron focusing implies a mean free path $\ell\gtrsim3.5\mathrm{\mu m}$. Superconductivity is thus deep in the clean limit, with the disorder parameter $d=\xi/\ell<0.1$. SC1 emerge from a paramagnetic normal state and is suppressed with in-plane magnetic fields in agreement with the Pauli paramagnetic limit. In contrast, SC2 emerges from a spin-polarized, valley-unpolarized half-metal. Measurements of the in-plane critical field show that this superconductor exceeds the Pauli limit by at least one order of magnitude. We discuss our results in light of several mechanisms including conventional phonon-mediated pairing, pairing due to fluctuations of the proximal isospin order, and intrinsic instabilities of the annular Fermi liquid. Our observation of superconductivity in a clean and structurally simple two-dimensional metal hosting a variety of gate tuned magnetic states may enable a new class of field-effect controlled mesoscopic electronic devices combining correlated electron phenomena.


title: ShapeFit: extracting the power spectrum shape information in galaxy surveys beyond BAO and RSD

abstract: In the standard (classic) approach, galaxy clustering measurements from spectroscopic surveys are compressed into baryon acoustic oscillations and redshift space distortions measurements, which in turn can be compared to cosmological models. Recent works have shown that avoiding this intermediate step and fitting directly the full power spectrum signal (full modelling) leads to much tighter constraints on cosmological parameters. Here we show where this extra information is coming from and extend the classic approach with one additional effective parameter, such that it captures, effectively, the same amount of information as the full modelling approach, but in a model-independent way. We validate this new method (ShapeFit) on mock catalogs, and compare its performance to the full modelling approach finding both to deliver equivalent results. The ShapeFit extension of the classic approach promotes the standard analyses at the level of full modelling ones in terms of information content, with the advantages of i) being more model independent; ii) offering an understanding of the origin of the extra cosmological information; iii) allowing a robust control on the impact of observational systematics.


title: Self-symmetrization of symmetry-breaking dynamics in passive Kerr resonators

abstract: The realization of spontaneous symmetry breaking (SSB) requires a system that exhibits a near perfect symmetry. SSB manifests itself through a pitchfork bifurcation, but that bifurcation is fragile, and perturbed by any asymmetry or imperfections. Consequently, exploiting SSB for real-world applications is challenging and often requires cumbersome stabilization techniques. Here, we reveal a novel method that automatically leads to symmetric conditions, and demonstrate its practical application in coherently-driven, two-mode, passive Kerr resonators. More specifically, we show that introducing a $\pi$-phase defect between the modes of a driven nonlinear resonator makes SSB immune to asymmetries by means of a period-doubled dynamics of the system's modal evolution. The two-roundtrip evolution induces a self-symmetrization of the system through averaging of the parameters, hence enabling the realization of SSB with unprecedented robustness. This mechanism is universal: all symmetry-broken solutions of driven Kerr resonators have a period-doubled counterpart. We experimentally demonstrate this universality by considering the polarization symmetry breaking of several different nonlinear structures found in normal and anomalous dispersion fiber cavities, including homogeneous states, polarization domain walls, and bright vector cavity solitons.


title: Unsupervised Learning of Visual 3D Keypoints for Control

abstract: Learning sensorimotor control policies from high-dimensional images crucially relies on the quality of the underlying visual representations. Prior works show that structured latent space such as visual keypoints often outperforms unstructured representations for robotic control. However, most of these representations, whether structured or unstructured are learned in a 2D space even though the control tasks are usually performed in a 3D environment. In this work, we propose a framework to learn such a 3D geometric structure directly from images in an end-to-end unsupervised manner. The input images are embedded into latent 3D keypoints via a differentiable encoder which is trained to optimize both a multi-view consistency loss and downstream task objective. These discovered 3D keypoints tend to meaningfully capture robot joints as well as object movements in a consistent manner across both time and 3D space. The proposed approach outperforms prior state-of-art methods across a variety of reinforcement learning benchmarks. Code and videos at https://buoyancy99.github.io/unsup-3d-keypoints/


title: A Continuized View on Nesterov Acceleration for Stochastic Gradient Descent and Randomized Gossip

abstract: We introduce the continuized Nesterov acceleration, a close variant of Nesterov acceleration whose variables are indexed by a continuous time parameter. The two variables continuously mix following a linear ordinary differential equation and take gradient steps at random times. This continuized variant benefits from the best of the continuous and the discrete frameworks: as a continuous process, one can use differential calculus to analyze convergence and obtain analytical expressions for the parameters; and a discretization of the continuized process can be computed exactly with convergence rates similar to those of Nesterov original acceleration. We show that the discretization has the same structure as Nesterov acceleration, but with random parameters. We provide continuized Nesterov acceleration under deterministic as well as stochastic gradients, with either additive or multiplicative noise. Finally, using our continuized framework and expressing the gossip averaging problem as the stochastic minimization of a certain energy function, we provide the first rigorous acceleration of asynchronous gossip algorithms.


title: PhyMask: Robust Sensing of Brain Activity and Physiological Signals During Sleep with an All-textile Eye Mask

abstract: Clinical-grade wearable sleep monitoring is a challenging problem since it requires concurrently monitoring brain activity, eye movement, muscle activity, cardio-respiratory features and gross body movements. This requires multiple sensors to be worn at different locations as well as uncomfortable adhesives and discrete electronic components to be placed on the head. As a result, existing wearables either compromise comfort or compromise accuracy in tracking sleep variables. We propose PhyMask, an all-textile sleep monitoring solution that is practical and comfortable for continuous use and that acquires all signals of interest to sleep solely using comfortable textile sensors placed on the head. We show that PhyMask can be used to accurately measure sleep stages and advanced sleep markers such as spindles and k-complexes robustly in the real-world setting. We validate PhyMask against polysomnography and show that it significantly outperforms two commercially-available sleep tracking wearables, Fitbit and Oura Ring.


title: Comments on Sen's restricted preferences and rational choice

abstract: A discussion about conditions on the profiles of voters in order to avoid contradictory decisions that have been designed by Sen and others.


title: A Large-scale Approach to Modelling Molecular Biosignatures: The Diatomics

abstract: This work presents the first steps to modelling synthetic rovibrational spectra for all molecules of astrophysical interest using a new approach implemented in the Prometheus code. The goal is to create a new comprehensive source of first-principles molecular spectra, thus bridging the gap for missing data to help drive future high-resolution studies. Our primary application domain is for molecules identified as signatures of life in planetary atmospheres (biosignatures), but our approach is general and can be applied to other systems. In this work we evaluate the accuracy of our method by studying four diatomic molecules H$_2$, O$_2$, N$_2$ and CO, all of which have well-known spectra. Prometheus uses the Transition-Optimised Shifted Hermite (TOSH) theory to account for anharmonicity for the fundamental $\nu=0 \rightarrow \nu=1$ band, along with thermal profile modeling for the rotational transitions. To this end, we expand TOSH theory to enable the modeling of rotational constants. We show that our simple model achieves results that are a better approximation of the real spectra than those produced through a harmonic approach. We compare our results with high-resolution HITRAN and ExoMol spectral data. We find that modelling accuracy tends to diminish for rovibrational transition away from the band origin, thus highlighting the need for the theory to be further adapted.


title: LBT Reveals Large Dust Particles and a High Mass Loss Rate for K2-22 b

abstract: The disintegrating planet candidate K2-22 b shows periodic and stochastic transits best explained by an escaping debris cloud. However, the mechanism that creates the debris cloud is unknown. The grain size of the debris as well as its sublimation rate can be helpful in understanding the environment that disintegrates the planet. Here, we present simultaneous photometry with the g band at 0.48 microns and KS band at 2.1 microns using the Large Binocular Telescope. During an event with very low dust activity, we put a new upper limit on the size of the planet of 0.71 earth radii or 4500 km. We also detected a medium-depth transit which can be used to constrain the dust particle sizes. We find that the median particle size must be larger than about 0.5 to 1.0 microns, depending on the composition of the debris. This leads to a high mass loss rate of about 3e8 kg/s that is consistent with hydrodynamic escape models. If they are produced by some alternate mechanism such as explosive volcanism, it would require extraordinary geological activity. Combining our upper limits on the planet size with the high mass loss rate, we find a lifetime of the planet of less than 370 Myr. This drops to just 21 Myr when adopting the 0.02 earth masses predicted from hydrodynamical models.


title: Non-Lorentzian Chaos and Cosmological Holography

abstract: We study chaos in non-Lorentzian field theories, specifically Galilean and Carrollian conformal field theories in two dimensions. In a large central charge limit, we find that the Lyapunov exponent saturates the bound on chaos, conjectured originally for relativistic field theories. We recover the same Lyapunov exponent holographically by a shock-wave calculation in three-dimensional flat space cosmologies, providing further evidence for flat space holography.


title: The FL bound and its phenomenological implications

abstract: Demanding that charged Nariai black holes in (quasi-)de Sitter space decay without becoming super-extremal implies a lower bound on the masses of charged particles, known as the Festina Lente (FL) bound. In this paper we fix the $\mathcal{O}(1)$ constant in the bound and elucidate various aspects of it, as well as extensions to $d>4$ and to situations with scalar potentials and dilatonic couplings. We also discuss phenomenological implications of FL including an explanation of why the Higgs potential cannot have a local minimum at the origin, thus explaining why the weak force must be broken. For constructions of meta-stable dS involving anti-brane uplift scenarios, even though the throat region is consistent with FL, the bound implies that we cannot have any light charged matter fields coming from any far away region in the compactified geometry, contrary to the fact that they are typically expected to arise in these scenarios. This strongly suggests that introduction of warped anti-branes in the throat cannot be decoupled from the bulk dynamics as is commonly assumed. Finally, we provide some evidence that in certain situations the FL bound can have implications even with gravity decoupled and illustrate this in the context of non-compact throats.


title: Double Copy Relation for AdS

abstract: We present a double copy relation in AdS$_5$ which relates tree-level four-point amplitudes of supergravity, super Yang-Mills and bi-adjoint scalars.


title: Globular clusters as tracers of the dark matter halo: insights from the E-MOSAICS simulations

abstract: Globular clusters (GCs) are bright objects that span a wide range of galactocentric distances, and are thus probes of the structure of dark matter (DM) haloes. In this work, we explore whether the projected radial profiles of GCs can be used to infer the structural properties of their host DM haloes. We use the simulated GC populations in a sample of 166 central galaxies from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. We find that more massive galaxies host stellar and GC populations with shallower density profiles that are more radially extended. In addition, the metal-poor GC subpopulations tend to have shallower and more extended profiles than the metal-rich subsamples, which we relate to the preferentially accreted origin of the metal-poor GCs. We find strong correlations between the slopes and effective radii of the radial profiles of the GC populations and the structural properties of the DM haloes, such as their power-law slopes, scale radii, and concentration parameters. Accounting for a dependence on the galaxy stellar mass decreases the scatter of the two-dimensional relations. This suggests that the projected number counts of GCs, combined with their galaxy mass, trace the density profile of the DM halo of their host galaxy. When applied to extragalactic GC systems, we recover the scale radii and the extent of the DM haloes of a sample of ETGs with uncertainties smaller than $0.2~\rm dex$. Thus, extragalactic GC systems provide a novel avenue to explore the structure of DM haloes beyond the Local Group.


title: Selection Functions in Astronomical Data Modeling, with the Space Density of White Dwarfs as Worked Example

abstract: Statistical studies of astronomical data sets, in particular of cataloged properties for discrete objects, are central to astrophysics. One cannot model those objects' population properties or incidences without a quantitative understanding of the conditions under which these objects ended up in a catalog or sample, the sample's selection function. As systematic and didactic introductions to this topic are scarce in the astrophysical literature, we aim to provide one, addressing generically the following questions: What is a selection function? What arguments $\vec{q}$ should a selection function depend on? Over what domain must a selection function be defined? What approximations and simplifications can be made? And, how is a selection function used in `modelling'? We argue that volume-complete samples, with the volume drastically curtailed by the faintest objects, reflect a highly sub-optimal selection function that needlessly reduces the number of bright and usually rare objects in the sample. We illustrate these points by a worked example, deriving the space density of white dwarfs (WD) in the Galactic neighbourhood as a function of their luminosity and Gaia color, $\Phi_0(M_G,B-R)$ in [mag$^{-2}$pc$^{-3}$]. We construct a sample of $10^5$ presumed WDs through straightforward selection cuts on the Gaia EDR3 catalog, in magnitude, color, parallax, and astrometric fidelity $\vec{q}=(m_G,B-R,\varpi,p_{af})$. We then combine a simple model for $\Phi_0$ with the effective survey volume derived from this selection function $S_C(\vec{q})$ to derive a detailed and robust estimate of $\Phi_0(M_G,B-R)$. This resulting white dwarf luminosity-color function $\Phi_0(M_G,B-R)$ differs dramatically from the initial number density distribution in the luminosity-color plane: by orders of magnitude in density and by four magnitudes in density peak location.


title: Higher-Form Symmetries and Their Anomalies in M-/F-Theory Duality

abstract: We explore higher-form symmetries of M- and F-theory compactified on elliptic fibrations, determined by the topology of their asymptotic boundaries. The underlying geometric structures are shown to be equivalent to known characterizations of the gauge group topology in F-theory via Mordell--Weil torsion and string junctions. We further study dimensional reductions of the 11d Chern--Simons term in the presence of torsional boundary $G_4$-fluxes, which encode background gauge fields of center 1-form symmetries in the lower-dimensional effective gauge theory. We find contributions that can be interpreted as 't Hooft anomalies involving the 1-form symmetry which originate from a fractionalization of the instanton number of non-Abelian gauge theories in F-/M-theory compactifications to 8d/7d and 6d/5d.


title: StelNet: Hierarchical Neural Network for Automatic Inference in Stellar Characterization

abstract: Characterizing the fundamental parameters of stars from observations is crucial for studying the stars themselves, their planets, and the galaxy as a whole. Stellar evolution theory predicting the properties of stars as a function of stellar age and mass enables translating observables into physical stellar parameters by fitting the observed data to synthetic isochrones. However, the complexity of overlapping evolutionary tracks often makes this task numerically challenging, and with a precision that can be highly variable, depending on the area of the parameter space the observation lies in. This work presents StelNet, a Deep Neural Network trained on stellar evolutionary tracks that quickly and accurately predicts mass and age from absolute luminosity and effective temperature for stars with close to solar metallicity. The underlying model makes no assumption on the evolutionary stage and includes the pre-main sequence phase. We use bootstrapping and train many models to quantify the uncertainty of the model. To break the model's intrinsic degeneracy resulting from overlapping evolutionary paths, we also built a hierarchical model that retrieves realistic posterior probability distributions of the stellar mass and age. We further test and train StelNet using a sample of stars with well-determined masses and ages from the literature.


title: The M-dwarf Ultraviolet Spectroscopic Sample I. Determining Stellar Parameters for Field Stars

abstract: Accurate stellar properties are essential for precise stellar astrophysics and exoplanetary science. In the M-dwarf regime, much effort has gone into defining empirical relations that can use readily accessible observables to assess physical stellar properties. Often, these relations for the quantity of interest are cast as a non-linear function of available data; however, in Bayesian modeling the reverse is needed. In this article, we introduce a new Bayesian framework to self-consistently and simultaneously apply multiple empirical calibrations to fully characterize the mass, luminosity, radius, and effective temperature of a field age M-dwarf. This framework includes a new M-dwarf mass-radius relation with a scatter of 3.1% at fixed mass. We further introduce the M-dwarf Ultraviolet Spectroscopic Sample (MUSS), and apply our methodology to provide consistent stellar parameters for these nearby low-mass stars, selected as having available spectroscopic data in the ultraviolet. These targets are of interest largely as either exoplanet hosts or benchmarks in multi-wavelength stellar activity. We use the field MUSS stars to define a low-mass main sequence in the solar neighborhood through Gaussian Process (GP) regression. These results enable us to empirically measure a feature in the GP derivative at $\mathcal{M}_{b} = 0.337 \pm ^{0.013}_{0.026}$ $M_{\odot}$ that indicates where the M-dwarf Ultraviolet Spectroscopic Sample transitions from fully to partly convective interiors.


title: Multiwavelength evidence for a new flare-mode transitional millisecond pulsar

abstract: We report the discovery of a new low-mass X-ray binary near the center of the unassociated Fermi GeV gamma-ray source 4FGL J0540.0-7552. The source shows the persistent presence of an optical accretion disk and exhibits extreme X-ray and optical variability. It also has an X-ray spectrum well-fit by a hard power law with a Gamma = 1.8 and a high ratio of X-ray to gamma-ray flux. Together, these properties are consistent with the classification of the binary as a transitional millisecond pulsar (tMSP) in the sub-luminous disk state. Uniquely among the candidate tMSPs, 4FGL J0540.0-7552 shows consistent optical, X-ray, and gamma-ray evidence for having undergone a state change, becoming substantially brighter in the optical and X-rays and fainter in GeV gamma-rays sometime in mid-2013. In its current sub-luminous disk state, and like one other candidate tMSP in the Galactic field, 4FGL J0540.0-7552 appears to always be in an X-ray "flare mode", indicating that this could be common phenomenology for tMSPs.


title: Disruption of Hierarchical Clustering in the Vela OB2 Complex and the Cluster Pair Collinder 135 and UBC7 with Gaia EDR3: Evidence of Supernova Quenching

abstract: We identify hierarchical structures in the Vela OB2 complex and the cluster pair Collinder 135 and UBC 7 with Gaia EDR3 using the neural network machine learning algorithm StarGO. Five second-level substructures are disentangled in Vela OB2, which are referred to as Huluwa 1 (Gamma Velorum), Huluwa 2, Huluwa 3, Huluwa 4 and Huluwa 5. For the first time, Collinder 135 and UBC 7 are simultaneously identified as constituent clusters of the pair with minimal manual intervention. We propose an alternative scenario in which Huluwa 1-5 have originated from sequential star formation. The older clusters Huluwa 1-3 with an age of 10-22 Myr, generated stellar feedback to cause turbulence that fostered the formation of the younger-generation Huluwa 4-5 (7-20 Myr). A supernova explosion located inside the Vela IRAS shell quenched star formation in Huluwa 4-5 and rapidly expelled the remaining gas from the clusters. This resulted in global mass stratification across the shell, which is confirmed by the regression discontinuity method. The stellar mass in the lower rim of the shell is $0.32\pm0.14$ $\rm M_\odot$ higher than in the upper rim. Local, cluster-scale mass segregation is observed in the lowest-mass cluster Huluwa 5. Huluwa 1-5 (in Vela OB2) are experiencing significant expansion, while the cluster pair suffers from moderate expansion. The velocity dispersions suggest that all five groups (including Huluwa 1A and Huluwa 1B) in Vela OB2 and the cluster pair are supervirial and are undergoing disruption, and also that Huluwa 1A and Huluwa 1B may be a coeval young cluster pair. N-body simulations predict that Huluwa 1-5 in Vela OB2 and the cluster pair will continue to expand in the future 100 Myr and eventually dissolve.


title: Constraining Saturn's interior with ring seismology: effects of differential rotation and stable stratification

abstract: Normal mode oscillations in Saturn excite density and bending waves in the C Ring, providing a valuable window into the planet's interior. Saturn's fundamental modes (f modes) excite the majority of the observed waves, while gravito-inertial modes (rotationally modified g modes) associated with stable stratification in the deep interior provide a compelling explanation for additional density waves with low azimuthal wavenumbers m. However, multiplets of density waves with nearly degenerate frequencies, including an m=3 triplet, still lack a definitive explanation. We investigate the effects of rapid and differential rotation on Saturn's oscillations, calculating normal modes for independently constrained interior models. We use a non-perturbative treatment of rotation that captures the full effects of the Coriolis and centrifugal forces, and consequently the mixing of sectoral f modes with g modes characterized by very different spherical harmonic degrees. Realistic profiles for differential rotation associated with Saturn's zonal winds can enhance these mode interactions, producing detectable oscillations with frequencies separated by less than 1%. Our calculations demonstrate that a three-mode interaction involving an f mode and two g modes can feasibly explain the finely split m=3 triplet, although the fine-tuning required to produce such an interaction generally worsens agreement with seismological constraints provided by m=2 density waves. Our calculations additionally demonstrate that sectoral f mode frequencies are measurably sensitive to differential rotation in Saturn's convective envelope. Finally, we find that including realistic equatorial antisymmetry in Saturn's differential rotation profile couples modes with even and odd equatorial parity, producing oscillations that could in principle excite both density and bending waves simultaneously.


title: Searching for changing-state AGNs in massive datasets -- I: applying deep learning and anomaly detection techniques to find AGNs with anomalous variability behaviours

abstract: The classic classification scheme for Active Galactic Nuclei (AGNs) was recently challenged by the discovery of the so-called changing-state (changing-look) AGNs (CSAGNs). The physical mechanism behind this phenomenon is still a matter of open debate and the samples are too small and of serendipitous nature to provide robust answers. In order to tackle this problem, we need to design methods that are able to detect AGN right in the act of changing-state. Here we present an anomaly detection (AD) technique designed to identify AGN light curves with anomalous behaviors in massive datasets. The main aim of this technique is to identify CSAGN at different stages of the transition, but it can also be used for more general purposes, such as cleaning massive datasets for AGN variability analyses. We used light curves from the Zwicky Transient Facility data release 5 (ZTF DR5), containing a sample of 230,451 AGNs of different classes. The ZTF DR5 light curves were modeled with a Variational Recurrent Autoencoder (VRAE) architecture, that allowed us to obtain a set of attributes from the VRAE latent space that describes the general behaviour of our sample. These attributes were then used as features for an Isolation Forest (IF) algorithm, that is an anomaly detector for a "one class" kind of problem. We used the VRAE reconstruction errors and the IF anomaly score to select a sample of 8,809 anomalies. These anomalies are dominated by bogus candidates, but we were able to identify 75 promising CSAGN candidates.


title: Evidence for a single-layer van der Waals multiferroic

abstract: Multiferroic materials have garnered wide interest for their exceptional static and dynamical magnetoelectric properties. In particular, type-II multiferroics exhibit an inversion-symmetry-breaking magnetic order which directly induces a ferroelectric polarization through various mechanisms, such as the spin-current or the inverse Dzyaloshinskii-Moriya effect. This intrinsic coupling between the magnetic and dipolar order parameters results in record-strength magnetoelectric effects. Two dimensional materials possessing such intrinsic multiferroic properties have been long sought for harnessing magnetoelectric coupling in nanoelectronic devices. Here, we report the discovery of type-II multiferroic order in a single atomic layer of the transition metal-based van der Waals material NiI2. The multiferroic state of NiI2 is characterized by a proper-screw spin helix with given handedness, which couples to the charge degrees of freedom to produce a chirality-controlled electrical polarization. We use circular dichroic Raman measurements to directly probe the magneto-chiral ground state and its electromagnon modes originating from dynamic magnetoelectric coupling. Using birefringence and second-harmonic generation measurements, we detect a highly anisotropic electronic state simultaneously breaking three-fold rotational and inversion symmetry, and supporting polar order. The evolution of the optical signatures as a function of temperature and layer number surprisingly reveals an ordered magnetic, polar state that persists down to the ultrathin limit of monolayer NiI2. These observations establish NiI2 and transition metal dihalides as a new platform for studying emergent multiferroic phenomena, chiral magnetic textures and ferroelectricity in the two-dimensional limit.


title: Can Thorne-$\dot{\textrm{Z}}$ytkow Objects source GW190814-type events?

abstract: The LIGO-Virgo collaboration reported in their third run the coalescence event GW190814 involving a 2.6 $M_{\odot}$ object with a 23 $M_{\odot}$ black hole. In this letter we study the conditions under which Thorne-$\dot{\textrm{Z}}$ytkow objects (T$\dot{\textrm{Z}}$Os) can be connected to that type of events. We evaluate first the rate of appearance of T$\dot{\textrm{Z}}$Os in the local Universe. Under the assumption that T$\dot{\textrm{Z}}$Os eventually become low mass gap black holes we evaluate how those black holes end up in binaries with other stellar mass black holes and compare to the reported rate for GW190814-type of events (1-23 $\textrm{Gpc}^{-3} \textrm{yr}^{-1}$). We find that T$\dot{\textrm{Z}}$Os in dense stellar clusters can not explain the LIGO-Virgo rate without a T$\dot{\textrm{Z}}$Os population in the field providing a dominant contribution. We also find that T$\dot{\textrm{Z}}$Os formed within hierarchical triple systems in the field with the third more distant star being the progenitor of a stellar mass black hole, may be able to give a rate comparable to that of GW190814-type events. In that case, future observations should discover mergers between stellar mass and low mass gap black holes, with the lower mass spanning the entire low mass gap range.


title: Extending the evolution of the stellar mass-size relation at $z \leq 2$ to low stellar mass galaxies from HFF and CANDELS

abstract: We reliably extend the stellar mass-size relation over $0.2\leq z \leq2$ to low stellar mass galaxies by combining the depth of Hubble Frontier Fields (HFF) with the large volume covered by CANDELS. Galaxies are simultaneously modelled in multiple bands using the tools developed by the MegaMorph project, allowing robust size (i.e., half-light radius) estimates even for small, faint, and high redshift galaxies. We show that above 10$^7$M$_\odot$, star-forming galaxies are well represented by a single power law on the mass-size plane over our entire redshift range. Conversely, the stellar mass-size relation is steep for quiescent galaxies with stellar masses $\geq 10^{10.3}$M$_\odot$ and flattens at lower masses, regardless of whether quiescence is selected based on star-formation activity, rest-frame colours, or structural characteristics. This flattening occurs at sizes of $\sim1$kpc at $z\leq1$. As a result, a double power law is preferred for the stellar mass-size relation of quiescent galaxies, at least above 10$^7$M$_\odot$. We find no strong redshift dependence in the slope of the relation of star-forming galaxies as well as of high mass quiescent galaxies. We also show that star-forming galaxies with stellar masses $\geq$10$^{9.5}$M$_\odot$ and quiescent galaxies with stellar masses $\geq10^{10.3}$M$_\odot$ have undergone significant size growth since $z\sim2$, as expected; however, low mass galaxies have not. Finally, we supplement our data with predominantly quiescent dwarf galaxies from the core of the Fornax cluster, showing that the stellar mass-size relation is continuous below 10$^7$M$_\odot$, but a more complicated functional form is necessary to describe the relation.


title: Fractonic plaquette-dimer liquid beyond renormalization

abstract: We consider close-packed tiling models of geometric objects -- a mixture of hardcore dimers and plaquettes -- as a generalisation of the familiar dimer models. Specifically, on an anisotropic cubic lattice, we demand that each site be covered by either a dimer on a z-link or a plaquette in the x-y plane. The space of such fully packed tilings has an extensive degeneracy. This maps onto a fracton-type `higher-rank electrostatics', which can exhibit a plaquette-dimer liquid and an ordered phase. We analyse this theory in detail, using height representations and T-duality to demonstrate that the concomitant phase transition occurs due to the proliferation of dipoles formed by defect pairs. The resultant critical theory can be considered as a fracton version of the Kosterlitz-Thouless transition. A significant new element is its UV-IR mixing, where the low energy behavior of the liquid phase and the transition out of it is dominated by local (short-wavelength) fluctuations, rendering the critical phenomenon beyond the renormalization group paradigm.


title: Neutrino spin oscillations in conformally gravity coupling models and quintessence surrounding a black hole

abstract: In this paper, we study the spin transitions of neutrinos caused by the interaction with a gravitational field. We consider a model with a scalar field (describing screening effects) conformally coupled to matter and neutrinos. The presence of screening effects suppresses the neutrino spin-flip probability as compared with General Relativity predictions. Such a result could be used, combined with neutrino astronomy, for testing modified theories of gravity and, in turn, screening effects invoked to bypass the solar system and Lab tests. Such an analysis has been also extended to the case of the quintessence field surrounding a black hole. Here we investigate the flavor and spin transitions, showing that also in such a case exists a suppression of the effect compared to General Relativity prediction.


title: Bayesian estimation of our local motion from the Planck-2018 CMB temperature map

abstract: The largest fluctuation in the CMB sky is the CMB dipole, which is believed to be caused by the motion of our observation frame with respect to the CMB rest frame. This motion accounts for the known motion of the Solar System barycentre with a best-fit amplitude of $369$ km/s, in the direction ($\ell= 264^\circ$, $b=48^\circ$) in galactic coordinates. Along with the CMB dipole signal, this motion also causes an inevitable signature of statistical anisotropy in the higher multipoles due to the modulation and aberration of the CMB temperature and polarization fields. This leads to a correlation between adjacent CMB multipoles causing a non-zero value of the off-diagonal terms in the covariance matrix which can be captured in terms of the dipolar spectra of the bipolar spherical harmonics (BipoSH). In our work, we jointly infer the CMB power spectrum and the BipoSH spectrum in a Bayesian framework using the $\textit{Planck}$-2018 $\texttt{SMICA}$ temperature map. We detect amplitude and direction of the local motion consistent with the canonical value $v=369$ km/s inferred from CMB dipole with a statistical significance of $4.54\sigma$, $4.97\sigma$ and $5.23\sigma$ respectively from the masked temperature map with the available sky fraction $40.1\%$, $59.1\%$, and $72.2\%$, confirming the common origin of both the signals. The Bayes factor in favor of the canonical value is between $7$ to $8$ depending on the choice of mask. But it strongly disagrees (by a value of the Bayes factor about $10^{-10}-10^{-11}$) with a higher value of local motion which one can infer from the amplitude of the dipole signal obtained from the CatWISE2020 quasar catalog using the WISE and NEOWISE data set.


title: Identifying the $\nu=\frac{5}{2}$ topological order through charge transport measurements

abstract: We propose an experiment to identify the topological order of the $\nu=\frac{5}{2}$ state through a measurement of the electric conductance of a mesoscopic device. Our setup is based on interfacing $\nu=2, \ \frac{5}{2}$ and $3$ in the same device. Its conductance can unambiguously establish or rule out the particle-hole symmetric Pfaffian topological order, which is supported by recent thermal measurements. Additionally, it distinguishes between the Moore-Read and Anti-Pfaffian topological orders, which are favored by numerical calculations.


title: Disentangling (2+1)d Topological States of Matter with the Entanglement Negativity

abstract: We use the entanglement negativity, a bipartite measure of entanglement in mixed quantum states, to study how multipartite entanglement constrains the real-space structure of the ground state wavefunctions of $(2+1)$-dimensional topological phases. We focus on the (Abelian) Laughlin and (non-Abelian) Moore-Read states at filling fraction $\nu=1/m$. We show that a combination of entanglement negativities, calculated with respect to specific cylinder and torus geometries, determines a necessary condition for when a topological state can be disentangled, i.e., factorized into a tensor product of states defined on cylinder subregions. This condition, which requires the ground state to lie in a definite topological sector, is sufficient for the Laughlin state. On the other hand, we find that a general Moore-Read ground state cannot be disentangled even when the disentangling condition holds.


title: A catalogue of white dwarfs in Gaia EDR3

abstract: We present a catalogue of white dwarf candidates selected from Gaia early data release three (EDR3). We applied several selection criteria in absolute magnitude, colour, and Gaia quality flags to remove objects with unreliable measurements while preserving most stars compatible with the white dwarf locus in the Hertzsprung-Russell diagram. We then used a sample of over 30 000 spectroscopically confirmed white dwarfs and contaminants from the Sloan Digital Sky Survey (SDSS) to map the distribution of these objects in the Gaia absolute magnitude-colour space. Finally, we adopt the same method presented in our previous Gaia DR2 work to calculate a probability of being a white dwarf (Pwd) for $\simeq$1.3 million sources which passed our quality selection. The Pwd values can be used to select a sample of $\simeq$359 000 high-confidence white dwarf candidates in the magnitude range 8< G <21. We calculated stellar parameters (effective temperature, surface gravity, and mass) for all these stars by fitting Gaia astrometry and photometry with synthetic models. We estimate an upper limit of 93 per cent for the overall completeness of our catalogue for white dwarfs with G $\leq$20 mag and effective temperature (Teff)>7000K, at high Galactic latitudes (|b|>20{\deg}). Alongside the main catalogue we include a reduced-proper-motion extension containing $\simeq$10 200 white dwarf candidates with unreliable parallax measurements which could, however be identified on the basis of their proper motion. We also performed a cross-match of our catalogues with SDSS DR16 spectroscopy and provide spectral classification based on visual inspection for all resulting matches.


title: BPS Indices, Modularity and Perturbations in Quantum K-theory

abstract: We study a perturbation family of N=2 3d gauge theories and its relation to quantum K-theory. A 3d version of the Intriligator-Vafa formula is given for the quantum K-theory ring of Grassmannians. The 3d BPS half-index of the gauge theory is connected to the theory of bilateral hypergeometric q-series, and to modular q-characters of a class of conformal field theories in a certain massless limit. Turning on 3d Wilson lines at torsion points leads to mock modular behavior. Perturbed correlators in the IR regime are computed by determining the UV-IR map in the presence of deformations.


title: Aspects of gravitational decoherence in neutrino lensing

abstract: We study decoherence effects in neutrino flavour oscillations in curved spacetime with particular emphasis on the lensing in a Schwarzschild geometry. Assuming Gaussian wave packets for neutrinos, we argue that the decoherence length derived from the exponential suppression of the flavour transition amplitude depends on the proper time of the geodesic connecting the events of the production and detection in general gravitational setting. In the weak gravity limit, the proper time between two events of given proper distance is smaller than that in the flat spacetime. Therefore, in presence of a Schwarzschild object, the neutrino wave packets have to travel relatively more physical distance in space to lapse the same amount of proper time before they decoher. For non-radial propagation applicable to the lensing phenomena, we show that the decoherence, in general, is sensitive to the absolute values of neutrino masses as well as the classical trajectories taken by neutrinos between the source and detector along with the spatial widths of neutrino wave packets. At distances beyond the decoherence length, the probability of neutrino flavour transition due to lensing attains a value which depends only on the leptonic mixing parameters. Hence, the observability of neutrino lensing significantly depends on these parameters and in-turn the lensing can provide useful information about them.


title: Lepton-driven Non-resonant Streaming Instability

abstract: A strong super-Alfv\'{e}nic drift of energetic particles (or cosmic rays, CRs) in a magnetized plasma can amplify the magnetic field significantly through non-resonant streaming instability (NRSI). While the traditional analysis is done for an ion current, here we use kinetic particle-in-cell simulations to study how the NRSI behaves when it is driven by electrons or by a mixture of electrons and positrons. In particular, we characterize growth rate, spectrum, and helicity of the unstable modes, as well the level of magnetic field at saturation. Our results are potentially relevant for several space/astrophysical environments (e.g, electron strahl in the solar wind, at oblique non-relativistic shocks, around pulsar wind nebulae) and also in laboratory experiments.


title: Unpredictability and entanglement in open quantum systems

abstract: We investigate dynamical many-body systems capable of universal computation, which leads to their properties being unpredictable unless the dynamics is simulated from the beginning to the end. Unpredictable behavior can be quantitatively assessed in terms of a data compression of the states occurring during the time evolution, which is closely related to their Kolmogorov complexity. We analyze a master equation embedding of classical cellular automata and demonstrate the existence of a phase transition between predictable and unpredictable behavior as a function of the random noise introduced by the embedding. We then turn to have this dynamics competing with a second process inducing quantum fluctuations and dissipatively driving the system to a highly entangled steady state. Strikingly, for intermediate strength of the quantum fluctuations, we find that both unpredictability and quantum entanglement can coexist even in the long time limit. Finally, we show that the required many-body interactions for the cellular automaton embedding can be efficiently realized within a variational quantum simulator platform based on ultracold Rydberg atoms with high fidelity.


title: Spontaneous Symmetry Breaking, Spectral Statistics, and the Ramp

abstract: Ensembles of quantum chaotic systems are expected to exhibit energy eigenvalues with random-matrix-like level repulsion between pairs of energies separated by less than the inverse Thouless time. Recent research has shown that exact and approximate global symmetries of a system have clear signatures in these spectral statistics, enhancing the spectral form factor or correspondingly weakening level repulsion. This paper extends those results to the case of spontaneous symmetry breaking, and shows that, surprisingly, spontaneously breaking a symmetry further enhances the spectral form factor. For both RMT-inspired toy models and models where the symmetry breaking has a description in terms of fluctuating hydrodynamics, we obtain formulas for this enhancement for arbitrary symmetry breaking patterns, including $Z_n$, $U(1)$, and partially or fully broken non-Abelian symmetries.


title: K-band GRAVITY/VLTI interferometry of "extreme" Herbig Be stars. The size-luminosity relation revisited

abstract: (Abridged:) It has been hypothesized that the location of Herbig Ae/Be stars (HAeBes) within the empirical relation between the inner disk radius (r$_{in}$), inferred from K-band interferometry, and the stellar luminosity (L$_*$), is related to the presence of the innermost gas, the disk-to-star accretion mechanism, the dust disk properties inferred from the spectral energy distributions (SEDs), or a combination of these effects. This work aims to test whether the previously proposed hypotheses do, in fact, serve as a general explanation for the distribution of HAeBes in the size-luminosity diagram. GRAVITY/VLTI spectro-interferometric observations at $\sim $2.2 $ \mu$m have been obtained for five HBes representing two extreme cases concerning the presence of innermost gas and accretion modes. V590 Mon, PDS 281, and HD 94509 show no excess in the near-ultraviolet, Balmer region of the spectra ($\Delta$D$_B$), indicative of a negligible amount of inner gas and disk-to-star accretion, whereas DG Cir and HD 141926 show such strong $\Delta$D$_B$ values that cannot be reproduced from magnetospheric accretion, but probably come from the alternative boundary layer mechanism. Additional data for these and all HAeBes resolved through K-band interferometry have been compiled from the literature and updated using Gaia EDR3 distances, almost doubling previous samples used to analyze the size-luminosity relation. We find no general trend linking the presence of gas inside the dust destruction radius or the accretion mechanism with the location of HAeBes in the size-luminosity diagram. Underlying trends are present and must be taken into account when interpreting the size-luminosity correlation. Still, it is argued that the size-luminosity correlation is most likely to be physically relevant in spite of the previous statistical warning concerning dependencies on distance.


title: Elliptic cohomology is unique up to homotopy

abstract: Homotopy theory folklore tells us that the sheaf defining the cohomology theory Tmf of topological modular forms is unique up to homotopy. Here we provide a proof of this fact, although we claim no originality for the statement. This retroactively reconciles all previous constructions of Tmf.


title: Planning to Fairly Allocate: Probabilistic Fairness in the Restless Bandit Setting

abstract: Restless and collapsing bandits are often used to model budget-constrained resource allocation in settings where arms have action-dependent transition probabilities, such as the allocation of health interventions among patients. However, state-of-the-art Whittle-index-based approaches to this planning problem either do not consider fairness among arms, or incentivize fairness without guaranteeing it. We thus introduce ProbFair, a probabilistically fair policy that maximizes total expected reward and satisfies the budget constraint while ensuring a strictly positive lower bound on the probability of being pulled at each timestep. We evaluate our algorithm on a real-world application, where interventions support continuous positive airway pressure (CPAP) therapy adherence among patients, as well as on a broader class of synthetic transition matrices. We find that ProbFair preserves utility while providing fairness guarantees.


title: Spectroscopic Evidence for the Direct Involvement of Local Moments in the Pairing Process of the Heavy-Fermion Superconductor CeCoIn$_5$

abstract: The microscopic mechanism for electron pairing in heavy-fermion superconductors remains a major challenge in quantum materials. Some form of magnetic mediation is widely accepted with spin fluctuations as a prime candidate. A novel mechanism, 'composite pairing' based on the cooperative two-channel Kondo effect directly involving the f-electron moments has also been proposed for some heavy fermion compounds including CeCoIn$_5$. The origin of the spin resonance peak observed in neutron scattering measurements on CeCoIn$_5$ is still controversial and the corresponding hump-dip structure in the tunneling conductance is missing. This is in contrast to the cuprate and Fe-based high-temperature superconductors, where both characteristic signatures are observed, indicating spin fluctuations are likely involved in the pairing process. Here, we report results from planar tunneling spectroscopy along three major crystallographic orientations of CeCoIn5 over wide ranges of temperature and magnetic field. The pairing gap opens at T$_p$ ~ 5 K, well above the bulk T$_c$ = 2.3 K, and its directional dependence is consistent with d$_{x^2-y^2}$ symmetry. With increasing magnetic field, this pairing gap is suppressed as expected but, intriguingly, a gaplike structure emerges smoothly, increasing linearly up to the highest field applied. This field-induced gaplike feature is only observed below T$_p$. The concomitant appearance of the pairing gap and the field-induced gaplike feature, along with its linear increase with field, indicates that the f-electron local moments are directly involved in the pairing process in CeCoIn$_5$.


title: The Effect of Client Appraisal on the Efficiency of Micro Finance Bank

abstract: One of the major problems confronting financial institutions most especially microfinance institutions is the increasing incidence of loan defaults and consequence loan losses which manifested in their financial performance with huge uncollectible loans and advances. This study assessed the effects of credit management on financial performance on microfinance institutions in Adamawa State, Nigeria. Specifically, we examine the effect of client appraisal on the efficiency of microfinance banks in Adamawa State. The methodology employed in this study is the survey method in which both primary and secondary sources were used in the collection of data. A multi-stage sampling method was adopted in selecting a sample of 21 respondents from a total population of 52 credit officers. Questionnaires were used in the due collection of data from the respondents. Descriptive statistics (simple percentage) and inferential statistics (regression analysis) were used to analyze the data collected and in testing the hypotheses. The study showed that client appraisal has a positive effect on efficiency and productivity.


title: The Effect of Inefficient Accretion on Planetary Differentiation

abstract: Pairwise collisions between terrestrial embryos are the dominant means of accretion during the last stage of planet formation. Hence, their realistic treatment in N-body studies is critical to accurately model the formation of terrestrial planets and to develop interpretations of telescopic and spacecraft observations. In this work, we compare the effects of two collision prescriptions on the core-mantle differentiation of terrestrial planets: a model in which collisions are always completely accretionary (``perfect merging'') and a more realistic model based on neural networks that has been trained on hydrodynamical simulations of giant impacts. The latter model is able to predict the loss of mass due to imperfect accretion and the evolution of non-accreted projectiles in hit-and-run collisions. We find that the results of the neural-network model feature a wider range of final core mass fractions and metal-silicate equilibration pressures, temperatures, and oxygen fugacities than the assumption of perfect merging. When used to model collisions in N-body studies of terrestrial planet formation, the two models provide similar answers for planets more massive than 0.1 Earth's masses. For less massive final bodies, however, the inefficient-accretion model predicts a higher degree of compositional diversity. This phenomenon is not reflected in planet formation models of the solar system that use perfect merging to determine collisional outcomes. Our findings confirm the role of giant impacts as important drivers of planetary diversity and encourage a realistic implementation of inefficient accretion in future accretion studies.


title: An MHD spectral theory approach to Jeans' magnetized gravitational instability

abstract: In this paper, we revisit the governing equations for linear magnetohydrodynamic (MHD) waves and instabilities existing within a magnetized, plane-parallel, self-gravitating slab. Our approach allows for fully non-uniformly magnetized slabs, which deviate from isothermal conditions, such that the well-known Alfv\'en and slow continuous spectra enter the description. We generalize modern MHD textbook treatments, by showing how self-gravity enters the MHD wave equation, beyond the frequently adopted Cowling approximation. This clarifies how Jeans' instability generalizes from hydro to magnetohydrodynamic conditions without assuming the usual Jeans' swindle approach. Our main contribution lies in reformulating the completely general governing wave equations in a number of mathematically equivalent forms, ranging from a coupled Sturm-Liouville formulation, to a Hamiltonian formulation linked to coupled harmonic oscillators, up to a convenient matrix differential form. The latter allows us to derive analytically the eigenfunctions of a magnetized, self-gravitating thin slab. In addition, as an example we give the exact closed form dispersion relations for the hydrodynamical p- and Jeans-unstable modes, with the latter demonstrating how the Cowling approximation modifies due to a proper treatment of self-gravity. The various reformulations of the MHD wave equation open up new avenues for future MHD spectral studies of instabilities as relevant for cosmic filament formation, which can e.g. use modern formal solution strategies tailored to solve coupled Sturm-Liouville or harmonic oscillator problems.


title: Revisiting Model Stitching to Compare Neural Representations

abstract: We revisit and extend model stitching (Lenc & Vedaldi 2015) as a methodology to study the internal representations of neural networks. Given two trained and frozen models $A$ and $B$, we consider a "stitched model'' formed by connecting the bottom-layers of $A$ to the top-layers of $B$, with a simple trainable layer between them. We argue that model stitching is a powerful and perhaps under-appreciated tool, which reveals aspects of representations that measures such as centered kernel alignment (CKA) cannot. Through extensive experiments, we use model stitching to obtain quantitative verifications for intuitive statements such as "good networks learn similar representations'', by demonstrating that good networks of the same architecture, but trained in very different ways (e.g.: supervised vs. self-supervised learning), can be stitched to each other without drop in performance. We also give evidence for the intuition that "more is better'' by showing that representations learnt with (1) more data, (2) bigger width, or (3) more training time can be "plugged in'' to weaker models to improve performance. Finally, our experiments reveal a new structural property of SGD which we call "stitching connectivity'', akin to mode-connectivity: typical minima reached by SGD can all be stitched to each other with minimal change in accuracy.


title: Extracting Global Dynamics of Loss Landscape in Deep Learning Models

abstract: Deep learning models evolve through training to learn the manifold in which the data exists to satisfy an objective. It is well known that evolution leads to different final states which produce inconsistent predictions of the same test data points. This calls for techniques to be able to empirically quantify the difference in the trajectories and highlight problematic regions. While much focus is placed on discovering what models learn, the question of how a model learns is less studied beyond theoretical landscape characterizations and local geometric approximations near optimal conditions. Here, we present a toolkit for the Dynamical Organization Of Deep Learning Loss Landscapes, or DOODL3. DOODL3 formulates the training of neural networks as a dynamical system, analyzes the learning process, and presents an interpretable global view of trajectories in the loss landscape. Our approach uses the coarseness of topology to capture the granularity of geometry to mitigate against states of instability or elongated training. Overall, our analysis presents an empirical framework to extract the global dynamics of a model and to use that information to guide the training of neural networks.


title: Measuring time with stationary quantum clocks

abstract: Time plays a fundamental role in our ability to make sense of the physical laws in the world around us. The nature of time has puzzled people -- from the ancient Greeks to the present day -- resulting in a long running debate between philosophers and physicists alike to whether time needs change to exist (the so-called relatival theory), or whether time flows regardless of change (the so-called substantival theory). One way to decide between the two is to attempt to measure the flow of time with a stationary clock, since if time were substantival, the flow of time would manifest itself in the experiment. Alas, conventional wisdom suggests that in order for a clock to function, it cannot be a static object, thus rendering this experiment seemingly impossible. Here we show, with the aid of counterfactual measurements, the surprising result that a quantum clock can measure the passage of time even while being switched off, thus lending constructive support for the substantival theory of time.


title: A Perspective on Collective Properties of Atoms on 2D Materials

abstract: Atoms deposited on two-dimensional (2D) electronic materials, such as graphene, can exhibit unconventional many-body correlations, not accessible in other settings. All of these are driven by van der Waals forces: between the atoms themselves and atom-material interactions. For example $^4$He atoms on 2D materials can potentially form a variety of exotic quantum states of matter, such as two-dimensional supersolids and superfluids, in addition to solid phases. For the "most quantum" case of a single helium layer we discuss, from a theoretical perspective, how the effective low-energy (Bose-Hubbard) description can take advantage of the extreme sensitivity of this unique system to the interplay between the atomic (helium) and solid-state (graphene) components. Due to the extraordinary variety and tunability of 2D electronic materials, we envisage that a wide range of correlated atomic phases can be realized under favorable conditions. We also outline exciting possibilities in the opposite extreme of many atomic layers forming a liquid on top of graphene -- in this case a so-called "spinodal de-wetting" pattern can form at the liquid-vapor interface which reflects the presence and electronic properties of graphene underneath. Such patterns could be manipulated by choosing different atoms and materials, with potential technological applications.


title: Tri-unitary quantum circuits

abstract: We introduce a novel class of quantum circuits that are unitary along three distinct "arrows of time". These dynamics share some of the analytical tractability of "dual-unitary" circuits, while exhibiting distinctive and richer phenomenology. We find that two-point correlations in these dynamics are strictly confined to three directions in $(1+1)$-dimensional spacetime -- the two light cone edges, $\delta x=\pm v\delta t$, and the static worldline $\delta x=0$. Along these directions, correlation functions are obtained exactly in terms of quantum channels built from the individual gates that make up the circuit. We prove that, for a class of initial states, entanglement grows at the maximum allowed speed up to an entropy density of at least one half of the thermal value, at which point it becomes model-dependent. Finally, we extend our circuit construction to $2+1$ dimensions, where two-point correlation functions are confined to the one-dimensional edges of a tetrahedral light cone -- a subdimensional propagation of information reminiscent of "fractonic" physics.

categories: quant-ph cond-mat.dis-nn cond-mat.stat-mech cond-mat.str-el hep-th


title: Deep neural networks for geometric multigrid methods

abstract: We investigate scaling and efficiency of the deep neural network multigrid method (DNN-MG). DNN-MG is a novel neural network-based technique for the simulation of the Navier-Stokes equations that combines an adaptive geometric multigrid solver, i.e. a highly efficient classical solution scheme, with a recurrent neural network with memory. The neural network replaces in DNN-MG one or multiple finest multigrid layers and provides a correction for the classical solve in the next time step. This leads to little degradation in the solution quality while substantially reducing the overall computational costs. At the same time, the use of the multigrid solver at the coarse scales allows for a compact network that is easy to train, generalizes well, and allows for the incorporation of physical constraints. Previous work on DNN-MG focused on the overall scheme and how to enforce divergence freedom in the solution. In this work, we investigate how the network size affects training and solution quality and the overall runtime of the computations. Our results demonstrate that larger networks are able to capture the flow behavior better while requiring only little additional training time. At runtime, the use of the neural network correction can even reduce the computation time compared to a classical multigrid simulation through a faster convergence of the nonlinear solve that is required at every time step.


title: Next Generation Reservoir Computing

abstract: Reservoir computing is a best-in-class machine learning algorithm for processing information generated by dynamical systems using observed time-series data. Importantly, it requires very small training data sets, uses linear optimization, and thus requires minimal computing resources. However, the algorithm uses randomly sampled matrices to define the underlying recurrent neural network and has a multitude of metaparameters that must be optimized. Recent results demonstrate the equivalence of reservoir computing to nonlinear vector autoregression, which requires no random matrices, fewer metaparameters, and provides interpretable results. Here, we demonstrate that nonlinear vector autoregression excels at reservoir computing benchmark tasks and requires even shorter training data sets and training time, heralding the next generation of reservoir computing.


title: Phase Transitions, Distance Functions, and Implicit Neural Representations

abstract: Representing surfaces as zero level sets of neural networks recently emerged as a powerful modeling paradigm, named Implicit Neural Representations (INRs), serving numerous downstream applications in geometric deep learning and 3D vision. Training INRs previously required choosing between occupancy and distance function representation and different losses with unknown limit behavior and/or bias. In this paper we draw inspiration from the theory of phase transitions of fluids and suggest a loss for training INRs that learns a density function that converges to a proper occupancy function, while its log transform converges to a distance function. Furthermore, we analyze the limit minimizer of this loss showing it satisfies the reconstruction constraints and has minimal surface perimeter, a desirable inductive bias for surface reconstruction. Training INRs with this new loss leads to state-of-the-art reconstructions on a standard benchmark.


title: On the Weyl Asymptotic Formula for Euclidean Domains of Finite Volume

abstract: We give a simple proof of the Weyl asymptotic formula for eigenvalues of the Dirichlet Laplacian, the buckling problem, and the Dirichlet bi-Laplacian in Euclidean domains of finite volume, with no assumptions about the boundary.


title: Improving Paraphrase Detection with the Adversarial Paraphrasing Task

abstract: If two sentences have the same meaning, it should follow that they are equivalent in their inferential properties, i.e., each sentence should textually entail the other. However, many paraphrase datasets currently in widespread use rely on a sense of paraphrase based on word overlap and syntax. Can we teach them instead to identify paraphrases in a way that draws on the inferential properties of the sentences, and is not over-reliant on lexical and syntactic similarities of a sentence pair? We apply the adversarial paradigm to this question, and introduce a new adversarial method of dataset creation for paraphrase identification: the Adversarial Paraphrasing Task (APT), which asks participants to generate semantically equivalent (in the sense of mutually implicative) but lexically and syntactically disparate paraphrases. These sentence pairs can then be used both to test paraphrase identification models (which get barely random accuracy) and then improve their performance. To accelerate dataset generation, we explore automation of APT using T5, and show that the resulting dataset also improves accuracy. We discuss implications for paraphrase detection and release our dataset in the hope of making paraphrase detection models better able to detect sentence-level meaning equivalence.


title: Purity and 2-Calabi-Yau categories

abstract: For various 2-Calabi-Yau categories $\mathscr{C}$ for which the stack of objects $\mathfrak{M}$ has a good moduli space $p\colon\mathfrak{M}\rightarrow \mathcal{M}$, we establish purity of the mixed Hodge module complex $p_{!}\underline{\mathbb{Q}}_{\mathfrak{M}}$. We do this by using formality in 2CY categories, along with \'etale neighbourhood theorems for stacks, to prove that the morphism $p$ is modelled \'etale-locally by the semisimplification morphism from the stack of modules of a preprojective algebra. Then via the integrality theorem in cohomological Donaldson-Thomas theory we prove purity of $p_{!}\underline{\mathbb{Q}}_{\mathfrak{M}}$. It follows that the Beilinson-Bernstein-Deligne-Gabber decomposition theorem for the constant sheaf holds for the morphism $p$, despite the possibly very singular and stacky nature of $\mathfrak{M}$. We use this to define cuspidal cohomology for $\mathfrak{M}$, which is conjecturally a complete space of generators for the BPS algebra associated to $\mathscr{C}$. We prove purity of the Borel-Moore homology of the moduli stack $\mathfrak{M}$, provided its good moduli space $\mathcal{M}$ is projective, or admits a suitable contracting $\mathbb{C}^*$-action. In particular, when $\mathfrak{M}$ is the moduli stack of Gieseker-semistable sheaves on a K3 surface, this proves a conjecture of Halpern-Leistner. We use these results to moreover prove purity for several stacks of coherent sheaves that do not admit a good moduli space. Without the usual assumption that $r$ and $d$ are coprime, we prove that the Borel-Moore homology of the stack of semistable degree $d$ rank $r$ Higgs sheaves is pure and carries a perverse filtration with respect to the Hitchin base, generalising the usual perverse filtration for the Hitchin system to the case of singular stacks of Higgs sheaves.


title: Stochastic resetting: A (very) brief review

abstract: Stochastic processes offer a fundamentally different paradigm of dynamics than deterministic processes, the most prominent example of the latter being Newton's laws of motion. Here, we discuss in a pedagogical manner a simple and illustrative example of stochastic processes in the form of a particle undergoing standard Brownian diffusion, with the additional feature of the particle resetting repeatedly and at random times to its initial condition. Over the years, many different variants of this simple setting have been studied, all of which serve as illustrations of non-trivial and interesting static and dynamic features that characterize stochastic dynamics at long times. We will provide in this work a brief overview of this active and rapidly evolving field by considering the arguably simplest example of Brownian diffusion in one dimension. Along the way, we will learn about some of the general techniques that a physicist employs to study stochastic processes. Relevant to the special issue, we will discuss in detail how introducing resetting in an otherwise diffusive dynamics provides an explicit optimization of the time to locate a target through a special choice of the resetting protocol. We also discuss thermodynamics of resetting, and provide a bird's eye view of some of the recent work in the field of resetting.


title: Modelling Settling-Driven Gravitational Instabilities at the Base of Volcanic Clouds Using the Lattice Boltzmann Method

abstract: Field observations and laboratory experiments have shown that ash sedimentation can be significantly affected by collective settling mechanisms that promote premature ash deposition, with important implications for associated impacts. Among these mechanisms, settling-driven gravitational instabilities result from the formation of a gravitationally-unstable particle boundary layer (PBL) that grows between volcanic ash clouds and the underlying atmosphere. The PBL destabilises once it reaches a critical thickness, triggering the formation of rapid, downward-moving ash fingers that remain poorly characterised. We simulate this process by coupling a Lattice Boltzmann model, which solves the Navier-Stokes equations for the fluid phase, with a Weighted Essentially Non Oscillatory (WENO) finite difference scheme which solves the advection-diffusion-settling equation describing particle transport. Since the physical problem is advection dominated, the use of the WENO scheme reduces numerical diffusivity and ensures accurate tracking of the temporal evolution of the interface between the layers. We have validated the new model by showing that the simulated early-time growth rate of the instability is in very good agreement with that predicted by linear stability analysis, whilst the modelled late-stage behaviour also successfully reproduces quantitative results from published laboratory experiments.


title: Adaptive normalization for IPW estimation

abstract: Inverse probability weighting (IPW) is a general tool in survey sampling and causal inference, used both in Horvitz-Thompson estimators, which normalize by the sample size, and H\'ajek/self-normalized estimators, which normalize by the sum of the inverse probability weights. In this work we study a family of IPW estimators, first proposed by Trotter and Tukey in the context of Monte Carlo problems, that are normalized by an affine combination of these two terms. We show how selecting an estimator from this family in a data-dependent way to minimize asymptotic variance leads to an iterative procedure that converges to an estimator with connections to regression control methods. We refer to this estimator as an adaptively normalized estimator. For mean estimation in survey sampling, this estimator has asymptotic variance that is never worse than the Horvitz--Thompson or H\'ajek estimators, and is smaller except in edge cases. Going further, we show that adaptive normalization can be used to propose improvements of the augmented IPW (AIPW) estimator, average treatment effect (ATE) estimators, and policy learning objectives. Appealingly, these proposals preserve both the asymptotic efficiency of AIPW and the regret bounds for policy learning with IPW objectives, and deliver consistent finite sample improvements in simulations for all three of mean estimation, ATE estimation, and policy learning.


title: Face Age Progression With Attribute Manipulation

abstract: Face is one of the predominant means of person recognition. In the process of ageing, human face is prone to many factors such as time, attributes, weather and other subject specific variations. The impact of these factors were not well studied in the literature of face aging. In this paper, we propose a novel holistic model in this regard viz., ``Face Age progression With Attribute Manipulation (FAWAM)", i.e. generating face images at different ages while simultaneously varying attributes and other subject specific characteristics. We address the task in a bottom-up manner, as two submodules i.e. face age progression and face attribute manipulation. For face aging, we use an attribute-conscious face aging model with a pyramidal generative adversarial network that can model age-specific facial changes while maintaining intrinsic subject specific characteristics. For facial attribute manipulation, the age processed facial image is manipulated with desired attributes while preserving other details unchanged, leveraging an attribute generative adversarial network architecture. We conduct extensive analysis in standard large scale datasets and our model achieves significant performance both quantitatively and qualitatively.


title: Memory effects in quantum dynamics modelled by quantum renewal processes

abstract: Simple, controllable models play an important role to learn how to manipulate and control quantum resources. We focus here on quantum non-Markovianity and model the evolution of open quantum systems by quantum renewal processes. This class of quantum dynamics provides us with a phenomenological approach to characterise dynamics with a variety of non-Markovian behaviours, here described in terms of the trace distance between two reduced states. By adopting a trajectory picture for the open quantum system evolution, we analyse how non-Markovianity is influenced by the constituents defining the quantum renewal process, namely the time-continuous part of the dynamics, the type of jumps and the waiting time distributions. We focus not only on the mere value of the non-Markovianity measure, but also on how different features of the trace distance evolution are altered, including times and number of revivals.


title: In-depth Analysis of Anisotropic Magnetoconductance in Bi$_2$Se$_3$ thin films with electron-electron interaction corrections

abstract: A combination of out-of-plane and in-plane magnetoconductance (MC) study in topological insulators (TI) is often used as an experimental technique to probe weak anti-localization (WAL) response of the topological surface states (TSSs). However, in addition to the above WAL response, weak localization (WL) contribution from conducting bulk states are also known to coexist and contribute to the overall MC; a study that has so far received limited attention. In this article, we accurately extract the above WL contribution by systematically analyzing the temperature and magnetic field dependency of conductivity in Bi$_2$Se$_3$ films. For accurate analysis, we quantify the contribution of electron-electron interactions to the measured MC which is often ignored in recent WAL studies. Moreover, we show that the WAL effect arising from the TSSs with finite penetration depth, for out-of-plane and in-plane magnetic field can together explain the anisotropic magnetoconductance (AMC) and, thus, the investigated AMC study can serve as a useful technique to probe the parameters like phase coherence length and penetration depth that characterise the TSSs in 3D TIs. We also demonstrate that increase in bulk-disorder, achieved by growing the films on amorphous SiO$_2$ substrate rather than on crystalline Al$_2$O$_3$(0001), can lead to stronger decoupling between the top and bottom surface states of the film.