title: Overturning negative construal of quantum superposition

abstract: Construal of observable facts or events, that is, the manner in which we understand reality, is based not only on mathematical formulas of a theory suggested as a reasonable explanation for physical phenomena (like general relativity or quantum mechanics), but also on a mathematical model of reasoning used to analyze and appraise statements regarding the objective world (for example, logic of one type or the other). Hence, every time that a certain construal of reality encounters a problem, there is a choice between a modification to the mathematical formalism of the physical theory and a change in the model of reasoning. A case in point is negative construal of quantum superposition causing the problem of definite outcomes. To be sure, according to the said construal, it is not the case that a system being in a superposition of states is exclusively in one of the states constituting the superposition, which in turn implies that macroscopically differing outcomes of observation may appear all at once. The usual approach to the problem of definite outcomes is to modify the quantum mathematical formalism by adding to it some extra postulates (for instance, the postulate of wave function collapse). However, since none of the extra postulates proposed so far has gained broad acceptance, one may try another avenue to resolve the problem, namely, to replace logic with an alternative mathematical model of reasoning. This possibility is studied in the present paper.


title: New magnetic chemically peculiar stars and candidates in the ATLAS First Catalog of Variable Stars

abstract: The number of known variable stars has increased by several magnitudes over the last decade, and automated classification routines are becoming increasingly important to cope with this development. Here we show that the "upside-down CBH variables", which were proposed as a potentially new class of variable stars by Heinze et al. (2018) in the ATLAS First Catalogue of Variable Stars, are, at least to a high percentage, made up of alpha2 Canum Venaticorum (ACV) variables - that is, photometrically variable magnetic chemically peculiar (CP2/He-peculiar) stars - with distinct double-wave light curves. Using suitable selection criteria, we identified 264 candidate ACV variables in the ATLAS variable star catalogue. 62 of these objects were spectroscopically confirmed with spectra from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (all new discoveries except for nine stars) and classified on the MK system. The other 202 stars are here presented as ACV star candidates that require spectroscopic confirmation. The vast majority of our sample of stars are main-sequence objects. Derived masses range from 1.4M(Sun) to 5M(Sun), with half our sample stars being situated in the range from 2 M(Sun) to 2.4 M(Sun), in good agreement with the spectral classifications. Most stars belong to the thin or thick disk; four objects, however, classify as members of the halo population. With a peak magnitude distribution at around 14th magnitude, the here presented stars are situated at the faint end of the known Galactic mCP star population. Our study highlights the need to consider rare variability classes, like ACV variables, in automated classification routines.


title: The Survival of Multiphase Dusty Clouds in Hot Winds

abstract: Much progress has been made recently in the acceleration of $\sim10^{4}$\,K clouds to explain absorption-line measurements of the circumgalactic medium and the warm, atomic phase of galactic winds. However, the origin of the cold, molecular phase in galactic winds has received relatively little theoretical attention. Studies of the survival of $\sim10^{4}$\,K clouds suggest efficient radiative cooling may enable the survival of expelled material from galactic disks. Alternatively, gas colder than 10$^4$\,K may form within the outflow, including molecules if dust survives the acceleration process. We explore the survival of dusty clouds in a hot wind with three-dimensional hydrodynamic simulations including radiative cooling and dust modeled as tracer particles. We find that cold $\sim10^{3}$\,K gas can be destroyed, survive, or transformed entirely to $\sim 10^4\,$K gas. We establish analytic criteria distinguishing these three outcomes which compare characteristic cooling times to the system's `cloud crushing' time. In contrast to typically studied $\sim10^{4}$\,K clouds, colder clouds are entrained faster than the drag time as a result of efficient mixing. We find that while dust can in principle survive embedded in the accelerated clouds, the survival fraction depends critically on the time dust spends in the hot phase and on the effective threshold temperature for destruction. We discuss our results in the context of polluting the circumgalactic medium with dust and metals, as well as understanding observations suggesting rapid acceleration of molecular galactic winds and ram pressure stripped tails of jellyfish galaxies.


title: Exploring the limits of the generation of non-classical states of spins coupled to a cavity by optimal control

abstract: We investigate the generation of non-classical states of spins coupled to a common cavity by means of a collective driving of the spins. We propose a control strategy using specifically designed series of short coherent and squeezing pulses, which have the key advantage of being experimentally implementable with the state-of-the art techniques. The parameters of the control sequence are found by means of optimization algorithms. We consider the cases of two and four spins, the goal being either to reach a well-defined target state or a state maximizing a measure of non-classicality. We discuss the influence of cavity damping and spin offset on the generation of non-classical states. We also explore to which extent squeezing fields help enhancing the efficiency of the control process.


title: On topological motives

abstract: Following Eilenberg-Steenrod axiomatic approach we construct the universal ordinary homology theory for any homological structure on a given category by representing ordinary theories with values in abelian categories. For a convenient category of spaces we then obtain a universal abelian category which can be actually described for CW-complexes as the category of hieratic modules.


title: Property-Aware Relation Networks for Few-Shot Molecular Property Prediction

abstract: Molecular property prediction plays a fundamental role in drug discovery to identify candidate molecules with target properties. However, molecular property prediction is essentially a few-shot problem which makes it hard to use regular machine learning models. In this paper, we propose Property-Aware Relation networks (PAR) to handle this problem. In comparison to existing works, we leverage the fact that both relevant substructures and relationships among molecules change across different molecular properties. We first introduce a property-aware embedding function to transform the generic molecular embeddings to substructure-aware space relevant to the target property. Further, we design an adaptive relation graph learning module to jointly estimate molecular relation graph and refine molecular embeddings w.r.t. the target property, such that the limited labels can be effectively propagated among similar molecules. We adopt a meta-learning strategy where the parameters are selectively updated within tasks in order to model generic and property-aware knowledge separately. Extensive experiments on benchmark molecular property prediction datasets show that PAR consistently outperforms existing methods and can obtain property-aware molecular embeddings and model molecular relation graph properly.


title: Hausdorff dimension of union of lines that cover a curve

abstract: We construct a continuously differentiable curve in the plane that can be covered by a collection of lines such that every line intersects the curve at a single point and the union of the lines has Hausdorff dimension 1. We show that for twice differentiable curves this is impossible. In that case, the union of the lines must have Hausdorff dimension 2. If we use only tangent lines then the differentiability of the curve already implies that the union of the lines must have Hausdorff dimension 2, unless the curve is a line. We also construct a continuous curve, which is in fact the graph of a strictly convex function, such that the union of (one sided) tangent lines has Hausdorff dimension 1.


title: The Excited Spin State of Dimorphos Resulting from the DART Impact

abstract: The NASA Double Asteroid Redirection Test (DART) mission is a planetary defense-driven test of a kinetic impactor on Dimorphos, the satellite of the binary asteroid 65803 Didymos. DART will intercept Dimorphos at a relative speed of ${\sim}6.5 \text{ km s}^{-1}$, perturbing Dimorphos's orbital velocity and changing the binary orbital period. We present three independent methods (one analytic and two numerical) to investigate the post-impact attitude stability of Dimorphos as a function of its axial ratios, $a/b$ and $b/c$ ($a \ge b \ge c$), and the momentum transfer efficiency $\beta$. The first method uses a novel analytic approach in which we assume a circular orbit and a point-mass primary that identifies four fundamental frequencies of motion corresponding to the secondary's mean motion, libration, precession, and nutation frequencies. At resonance locations among these four frequencies, we find that attitude instabilities are possible. Using two independent numerical codes, we recover many of the resonances predicted by the analytic model and indeed show attitude instability. With one code, we use fast Lyapunov indicators to show that the secondary's attitude can evolve chaotically near the resonance locations. Then, using a high-fidelity numerical model, we find that Dimorphos enters a chaotic tumbling state near the resonance locations and is especially prone to unstable rotation about its long axis, which can be confirmed by ESA's Hera mission arriving at Didymos in late 2026. We also show that a fully coupled treatment of the spin and orbital evolution of both bodies is crucial to accurately model the long-term evolution of the secondary's spin state and libration amplitude. Finally, we discuss the implications of a post-impact tumbling or rolling state, including the possibility of terminating BYORP evolution if Dimorphos is no longer in synchronous rotation.


title: Uncertainty Prediction for Machine Learning Models of Material Properties

abstract: Uncertainty quantification in Artificial Intelligence (AI)-based predictions of material properties is of immense importance for the success and reliability of AI applications in material science. While confidence intervals are commonly reported for machine learning (ML) models, prediction intervals, i.e., the evaluation of the uncertainty on each prediction, are seldomly available. In this work we compare 3 different approaches to obtain such individual uncertainty, testing them on 12 ML-physical properties. Specifically, we investigated using the Quantile loss function, machine learning the prediction intervals directly and using Gaussian Processes. We identify each approachs advantages and disadvantages and end up slightly favoring the modeling of the individual uncertainties directly, as it is the easiest to fit and, in most cases, minimizes over-and under-estimation of the predicted errors. All data for training and testing were taken from the publicly available JARVIS-DFT database, and the codes developed for computing the prediction intervals are available through JARVIS-Tools.


title: Spinning test particle motion around a traversable wormhole

abstract: The motion of spinning test particles around a traversable wormhole is investigated using the Mathisson Papapetrous Dixon equations, which couple the Riemann tensor with the antisymmetric tensor $S^{ab}$, related to the spin of the particle. Hence, we study the effective potential, circular orbits, and innermost stable circular orbit ISCO of spinning particles. We found that the spin affects significantly the location of the ISCO, in contrast with the motion of nonspinning particles, where the ISCO is the same in both the upper and lower universes. On the other hand, since the dynamical fourmomentum and kinematical fourvelocity of the spinning particle are not always parallel, we also consider a superluminal bound on the particles motion. In the case of circular orbits at the ISCO, we found that the motion of particles with an adimensional spin parameter lower greater than $s=-1.5$ $(1.5)$ is forbidden. The spin interaction becomes important for Kerr black hole orbiting super massive wormholes SMWH.


title: From block-Toeplitz matrices to differential equations on graphs: towards a general theory for scalable masked Transformers

abstract: In this paper we provide, to the best of our knowledge, the first comprehensive approach for incorporating various masking mechanisms into Transformers architectures in a scalable way. We show that recent results on linear causal attention (Choromanski et al., 2021) and log-linear RPE-attention (Luo et al., 2021) are special cases of this general mechanism. However by casting the problem as a topological (graph-based) modulation of unmasked attention, we obtain several results unknown before, including efficient d-dimensional RPE-masking and graph-kernel masking. We leverage many mathematical techniques ranging from spectral analysis through dynamic programming and random walks to new algorithms for solving Markov processes on graphs. We provide a corresponding empirical evaluation.


title: All the attention you need: Global-local, spatial-channel attention for image retrieval

abstract: We address representation learning for large-scale instance-level image retrieval. Apart from backbone, training pipelines and loss functions, popular approaches have focused on different spatial pooling and attention mechanisms, which are at the core of learning a powerful global image representation. There are different forms of attention according to the interaction of elements of the feature tensor (local and global) and the dimensions where it is applied (spatial and channel). Unfortunately, each study addresses only one or two forms of attention and applies it to different problems like classification, detection or retrieval. We present global-local attention module (GLAM), which is attached at the end of a backbone network and incorporates all four forms of attention: local and global, spatial and channel. We obtain a new feature tensor and, by spatial pooling, we learn a powerful embedding for image retrieval. Focusing on global descriptors, we provide empirical evidence of the interaction of all forms of attention and improve the state of the art on standard benchmarks.


title: Efficient Bayesian Sampling Using Normalizing Flows to Assist Markov Chain Monte Carlo Methods

abstract: Normalizing flows can generate complex target distributions and thus show promise in many applications in Bayesian statistics as an alternative or complement to MCMC for sampling posteriors. Since no data set from the target posterior distribution is available beforehand, the flow is typically trained using the reverse Kullback-Leibler (KL) divergence that only requires samples from a base distribution. This strategy may perform poorly when the posterior is complicated and hard to sample with an untrained normalizing flow. Here we explore a distinct training strategy, using the direct KL divergence as loss, in which samples from the posterior are generated by (i) assisting a local MCMC algorithm on the posterior with a normalizing flow to accelerate its mixing rate and (ii) using the data generated this way to train the flow. The method only requires a limited amount of \textit{a~priori} input about the posterior, and can be used to estimate the evidence required for model validation, as we illustrate on examples.


title: Power-law and log-normal avalanche size statistics in random growth processes

abstract: We study the avalanche statistics observed in a minimal random growth model. The growth is governed by a reproduction rate obeying a probability distribution with finite mean a and variance va. These two control parameters determine if the avalanche size tends to a stationary distribution, (Finite Scale statistics with finite mean and variance or Power-Law tailed statistics with exponent in (1, 3]), or instead to a non-stationary regime with Log-Normal statistics. Numerical results and their statistical analysis are presented for a uniformly distributed growth rate, which are corroborated and generalized by analytical results. The latter show that the numerically observed avalanche regimes exist for a wide family of growth rate distributions and provide a precise definition of the boundaries between the three regimes.


title: Some examples of calculation of massless and massive Feynman integrals

abstract: We show some examples of calculations of massless and massive Feynman integrals.


title: Weakly special test configurations of log canonical Fano varieties

abstract: Let $X$ be a strictly log canonical Fano variety, we show that every lc place of complements is dreamy, and there exists a correspondence between weakly special test configurations of $(X,-K_X)$ and lc places of complements.


title: Light Like Segment Compactification and Braneworlds with Dynamical String Tension

abstract: Braneworld scenarios appear in dynamical string tension theories, where string tension is positive between two surfaces separated by a short distance and at the two surfaces themselves the string tensions become infinite, therefore producing an effective confinement of the strings and therefore of all matter and gravity to the space between these to surfaces. The specific model studied is in the context of the modified measure formulation the string where tension appear as an additional dynamical degree of freedom and these tensions are not universal, but rather each string generates its own tension, which can have a different value for each string. We consider a new background field that can couple to these strings, the tension scalar is capable then of changing locally along the world sheet and then the value of the tension of the string object changes accordingly along the world sheet. Sting tension appears dynamically and each string tension has a characteristic constant of integration that appears as a consequence of the dynamical tension generation . When many types of strings probing the same region of space are considered this tension scalar is constrained by the requirement of quantum conformal invariance. For the case of two types of strings probing the same region of space with different dynamically generated tensions (characterized by different constants of integration) , there are two different metrics, associated to the different strings, that have to satisfy vacuum Einsteins equations and the consistency of these two Einsteins equations determine that the tension field behaves. The two metrics are flat space in Minkowshi space and flat space after a special conformal transformation, then the strings are confined inside a light like Segment or alternatively as expanding Braneworlds where the strings are confined between two expanding bubbles.


title: Geometry of Information: classical and quantum aspects

abstract: In this article, we describe various aspects of categorification of the structures appearing in information theory. These aspects include probabilistic models both of classical and quantum physics, emergence of F-manifolds, and motivic enrichments.


title: Assessing the polymer coil-globule state from the very first spectral modes

abstract: The determination of the coil-globule transition of a polymer is generally based on the reconstruction of scaling laws, implying the need for samples from a rather wide range of different polymer lengths $N$. The spectral point of view developed in this work allows for a very parsimonious description of all the aspects of the finite-size coil-globule transition on the basis of the first two Rouse (cosine) modes only, shedding new light on polymer theory. Capturing the relevant configuration path features, the proposed approach enables to determine the state of a polymer without the need of any information about the polymer length or interaction strength. Importantly, we propose an experimental implementation of our analysis that can be easily performed with modern fluorescent imaging techniques, and would allow differentiation of coil or globule conformations by simply recording the positions of three discernible loci on the polymer.


title: Geometric Optimal Controls for Flapping Wing UAV on a Lie Group

abstract: Inspired by flight characteristics captured from live Monarch butterflies, an optimal control problem is presented while accounting the effects of low-frequency flapping and abdomen undulation. A flapping-wing aerial vehicle is modeled as an articulated rigid body, and its dynamics are developed according to Lagrangian mechanics on an abstract Lie group. This provides an elegant, global formulation of the dynamics for flapping-wing aerial vehicles, avoiding complexities and singularities associated with local coordinates. This is utilized to identify an optimal periodic motion that minimizes energy variations, and an optimal control is formulated to stabilize the periodic motion. Furthermore, the outcome of this paper can be applied to optimal control for any Lagrangian system on a Lie group with a configuration-dependent inertia.


title: One-loop matrix elements of effective superstring interactions: $\alpha'$-expanding loop integrands

abstract: In the low-energy effective action of string theories, non-abelian gauge interactions and supergravity are augmented by infinite towers of higher-mass-dimension operators. We propose a new method to construct one-loop matrix elements with insertions of operators $D^{2k} F^n$ and $D^{2k} R^n$ in the tree-level effective action of type-I and type-II superstrings. Inspired by ambitwistor string theories, our method is based on forward limits of moduli-space integrals using string tree-level amplitudes with two extra points, expanded in powers of the inverse string tension $\alpha'$. Similar to one-loop ambitwistor computations, intermediate steps feature non-standard linearized Feynman propagators which eventually recombine to conventional quadratic propagators. With linearized propagators the loop integrand of the matrix elements obey one-loop versions of the monodromy and KLT relations. We express a variety of four- and five-point examples in terms of quadratic propagators and formulate a criterion on the underlying genus-one correlation functions that should make this recombination possible at all orders in $\alpha'$. The ultraviolet divergences of the one-loop matrix elements are crosschecked against the non-separating degeneration of genus-one integrals in string amplitudes. Conversely, our results can be used as a constructive method to determine degenerations of elliptic multiple zeta values and modular graph forms at arbitrary weight.


title: Strong, recent selective sweeps reshape genetic diversity in freshwater bivalve Megalonaias nervosa

abstract: Freshwater Unionid bivalves have recently faced ecological upheaval through pollution, barriers to dispersal, human harvesting, and changes in fish-host prevalence. Currently, over 70% of species are threatened, endangered or extinct. To characterize the genetic response to these recent selective pressures, we collected population genetic data for one successful bivalve species, Megalonaias nervosa. We identify megabase sized regions that are nearly monomorphic across the population, a signal of strong, recent selection reshaping genetic diversity. These signatures of selection encompass a total of 73Mb, greater response to selection than is commonly seen in population genetic models. We observe 102 duplicate genes with high dN/dS on terminal branches among regions with sweeps, suggesting that gene duplication is a causative mechanism of recent adaptation in M. nervosa. Genes in sweeps reflect functional classes known to be important for Unionid survival, including anticoagulation genes important for fish host parasitization, detox genes, mitochondria management, and shell formation. We identify selective sweeps in regions with no known functional impacts, suggesting mechanisms of adaptation that deserve greater attention in future work on species survival. In contrast, polymorphic transposable element insertions appear to be detrimental and underrepresented among regions with sweeps. TE site frequency spectra are skewed toward singleton variants, and TEs among regions with sweeps are present only at low frequency. Our work suggests that duplicate genes are an essential source of genetic novelty that has helped this successful species succeed in environments where others have struggled. These results suggest that gene duplications deserve greater attention in non-model population genomics, especially in species that have recently faced sudden environmental challenges.


title: Adaptive first-order methods revisited: Convex optimization without Lipschitz requirements

abstract: We propose a new family of adaptive first-order methods for a class of convex minimization problems that may fail to be Lipschitz continuous or smooth in the standard sense. Specifically, motivated by a recent flurry of activity on non-Lipschitz (NoLips) optimization, we consider problems that are continuous or smooth relative to a reference Bregman function - as opposed to a global, ambient norm (Euclidean or otherwise). These conditions encompass a wide range of problems with singular objectives, such as Fisher markets, Poisson tomography, D-design, and the like. In this setting, the application of existing order-optimal adaptive methods - like UnixGrad or AcceleGrad - is not possible, especially in the presence of randomness and uncertainty. The proposed method - which we call adaptive mirror descent (AdaMir) - aims to close this gap by concurrently achieving min-max optimal rates in problems that are relatively continuous or smooth, including stochastic ones.


title: Dynamical vortices in electron-phonon superconductors

abstract: We analyze the structure of an $s-$wave superconducting gap in systems with electron-phonon attraction and electron-electron repulsion. Earlier works have found that superconductivity develops despite strong repulsion, but the gap, $\Delta (\omega_m)$, necessarily changes sign along the Matsubara axis. We analyze the sign-changing gap function from a topological perspective using the knowledge that a nodal point of $\Delta (\omega_m)$ is the center of dynamical vortex. We consider two models with different cutoffs for the repulsive interaction and trace the vortex positions along the Matsubara axis and in the upper frequency half plane upon changing the relative strength of the attractive and repulsive parts of the interaction. We discuss how the presence of dynamical vortices affects the gap structure along the real axis, detectable in ARPES experiments.


title: Machine learning of Kondo physics using variational autoencoders and symbolic regression

abstract: We employ variational autoencoders to extract physical insight from a dataset of one-particle Anderson impurity model spectral functions. Autoencoders are trained to find a low-dimensional, latent space representation that faithfully characterizes each element of the training set, as measured by a reconstruction error. Variational autoencoders, a probabilistic generalization of standard autoencoders, further condition the learned latent space to promote highly interpretable features. In our study, we find that the learned latent variables strongly correlate with well known, but nontrivial, parameters that characterize emergent behaviors in the Anderson impurity model. In particular, one latent variable correlates with particle-hole asymmetry, while another is in near one-to-one correspondence with the Kondo temperature, a dynamically generated low-energy scale in the impurity model. Using symbolic regression, we model this variable as a function of the known bare physical input parameters and "rediscover" the non-perturbative formula for the Kondo temperature. The machine learning pipeline we develop suggests a general purpose approach which opens opportunities to discover new domain knowledge in other physical systems.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.stat-mech cs.LG


title: A Scenario-oriented Approach to Energy-Reserve Joint Procurement and Pricing

abstract: We consider some crucial problems related to the secure and reliable operation of power systems with high renewable penetrations: how much reserve should we procure, how should reserve resources distribute among different locations, and how should we price reserve and charge uncertainty sources. These issues have so far been largely addressed empirically. In this paper, we first develop a scenario-oriented energy-reserve co-optimization model, which directly connects reserve procurement with possible outages and load/renewable power fluctuations without the need for empirical reserve requirements. Accordingly, reserve can be optimally procured system-wide to handle all possible future uncertainties with the minimum expected system total cost. Based on the proposed model, marginal pricing approaches are developed for energy and reserve, respectively. Locational uniform pricing is established for energy, and the similar property is also established for the combination of reserve and re-dispatch. In addition, properties of cost recovery for generators and revenue adequacy for the system operator are also proven.


title: A study of events with photoelectric emission in the DarkSide-50 liquid argon Time Projection Chamber

abstract: Finding unequivocal evidence of dark matter interactions in a particle detector is a major objective of physics research. Liquid argon time projection chambers offer a path to probe Weakly Interacting Massive Particles scattering cross sections on nucleus down to the so-called neutrino floor, in a mass range from few GeV's to hundredths of TeV's. Based on the successful operation of the DarkSide-50 detector at LNGS, a new and more sensitive experiment, DarkSide-20k, has been designed and is now under construction. A thorough understanding of the DarkSide-50 detector response and, therefore, of all kind of observed events, is essential for an optimal design of the new experiment. In this paper, we report on a particular set of events, which were not used for dark matter searches. Namely, standard two-pulse scintillation-ionization signals accompanied by a small amplitude third pulse, originating from single or few electrons, in a time window of less than a maximum drift time. We compare our findings to those of a recent paper of the LUX Collaboration (D.S.Akerib et al. Phys.Rev.D 102, 092004). Indeed, both experiments observe events related to photoionization of the cathode. From the measured rate of these events, we estimate for the first time the quantum efficiency of the tetraphenyl butadiene deposited on the DarkSide-50 cathode at wavelengths around 128 nm, in liquid argon. Also, both experiments observe events likely related to photoionization of impurities in the liquid. The probability of photoelectron emission per unit length turns out to be one order of magnitude smaller in DarkSide-50 than in LUX. This result, together with the much larger measured electron lifetime, coherently hints toward a lower concentration of contaminants in DarkSide-50 than in LUX.


title: Investigation the Effect of Velocity inlet and Carrying Fluid on the Flow inside Coronary Artery

abstract: In this study OpenFOAM 4.4.2 was used to investigate flow inside the coronary artery of the heart. This step is the first step of our future project, which is to include conjugate heat transfer of the heart with three main coronary arteries. Three different velocities were used as inlet boundary conditions to see the effect of velocity increase on velocity, pressure, and wall shear of coronary artery. Also, three different fluids namely the University of Wisconsin solution, Gelatin, and blood was used to investigate the effect of different fluids on flow inside coronary artery. A code based on Reynolds Stress Navier Stokes (RANS) equations was written and was implemented with the real boundary condition that was calculated based on MRI images. In order to improve the accuracy of the current numerical scheme, hex dominant mesh is utilized. When the inlet velocity increases to 0.5 m/s, velocity, wall shear stress, and pressure increase at the narrower parts.


title: Plasmonic Quantum Dots in Twisted Bilayer Graphene

abstract: We derive a material-realistic real-space many-body Hamiltonian for twisted bilayer graphene from first principles, including both single-particle hopping terms for $p_z$ electrons and long-range Coulomb interactions. By disentangling low- and high-energy subspaces of the electronic dispersion, we are able to utilize state-of-the-art constrained Random Phase Approximation calculations to reliably describe the non-local background screening from the high-energy $s$, $p_x$, and $p_y$ electron states for arbitrary twist angles. The twist-dependent low-energy screening from $p_z$ states is subsequently added to obtain a full screening model. We use this approach to study real-space plasmonic patterns in electron-doped twisted bilayer graphene supercells and find, next to classical dipole-like modes, also twist-angle-dependent plasmonic quantum-dot-like excitations with $s$ and $p$ symmetries. Based on their inter-layer charge modulations and their footprints in the electron energy loss spectrum, we can classify these modes into "bright" and "dark" states, which show different dependencies on the twist angle.


title: Long-Range Superconducting Proximity Effect in Nickel Nanowires

abstract: When a ferromagnet is placed in contact with a superconductor, owing to incompatible spin order, the Cooper pairs from the superconductor cannot survive more than one or two nanometers inside the ferromagnet. This is confirmed in the measurements of ferromagnetic nickel (Ni) nanowires contacted by superconducting niobium (Nb) leads. However, when a 3 nm thick copper oxide (CuO) buffer layer made by exposing an evaporated or a sputtered 3 nm Cu film to air, is inserted between the Nb electrodes and the Ni wire, the spatial extent of the superconducting proximity range is dramatically increased from 2 to a few tens of nanometers. Scanning transmission electron microscope study confirms the formation of a 3 nm thick CuO layer when an evaporated Cu film is exposed to air. Magnetization measurements of such a 3 nm CuO film on a SiO2/Si substrate and on Nb/SiO2/Si show clear evidence of ferromagnetism. One way to understand the long-range proximity effect in the Ni nanowire is that the CuO buffer layer with ferromagnetism facilitates the conversion of singlet superconductivity in Nb into triplet supercurrent along the Ni nanowires.


title: Bootstrapping Through Discrete Convolutional Methods

abstract: Bootstrapping was designed to randomly resample data from a fixed sample using Monte Carlo techniques. However, the original sample itself defines a discrete distribution. Convolutional methods are well suited for discrete distributions, and we show the advantages of utilizing these techniques for bootstrapping. The discrete convolutional approach can provide exact numerical solutions for bootstrap quantities, or at least mathematical error bounds. In contrast, Monte Carlo bootstrap methods can only provide confidence intervals which converge slowly. Additionally, for some problems the computation time of the convolutional approach can be dramatically less than that of Monte Carlo resampling. This article provides several examples of bootstrapping using the proposed convolutional technique and compares the results to those of the Monte Carlo bootstrap, and to those of the competing saddlepoint method.


title: Online Graph Topology Learning from Matrix-valued Time Series

abstract: This paper is concerned with the statistical analysis of matrix-valued time series. These are data collected over a network of sensors (typically a set of spatial locations) along time, where a vector of features is observed per time instant per sensor. Thus each sensor is characterized by a vectorial time series. We would like to identify the dependency structure among these sensors and represent it by a graph. When there is only one feature per sensor, the vector auto-regressive models have been widely adapted to infer the structure of Granger causality. The resulting graph is referred to as causal graph. Our first contribution is then extending VAR models to matrix-variate models to serve the purpose of graph learning. Secondly, we propose two online procedures respectively in low and high dimensions, which can update quickly the estimates of coefficients when new samples arrive. In particular in high dimensional regime, a novel Lasso-type is introduced and we develop its homotopy algorithms for the online learning. We also provide an adaptive tuning procedure for the regularization parameter. Lastly, we consider that, the application of AR models onto data usually requires detrending the raw data, however, this step is forbidden in online context. Therefore, we augment the proposed AR models by incorporating trend as extra parameter, and then adapt the online algorithms to the augmented data models, which allow us to simultaneously learn the graph and trend from streaming samples. In this work, we consider primarily the periodic trend. Numerical experiments using both synthetic and real data are performed, whose results support the effectiveness of the proposed methods.


title: Asteroseismology of evolved stars with EGGMiMoSA I. Theoretical mixed-mode patterns from the subgiant to the RGB phase

abstract: This study is the first of a series of papers that provide a technique to analyse the mixed-modes frequency spectra and characterise the structure of stars on the subgiant and red-giant branches. We define seismic indicators, relevant of the stellar structure and study their evolution on a grid of models. The proposed method, EGGMiMoSA, relies on the asymptotic description of mixed modes, defines initial guesses for the parameters, and uses a Levenberg-Marquardt technique to adjust the mixed-modes pattern efficiently. We follow the evolution of the mixed-modes parameters along a grid of models from the subgiant phase to the RGB bump and extend past works. We show the impact of the mass and composition on their evolution. The evolution of the period spacing $\Delta\pi_1$, pressure offset $\epsilon_p$, gravity offset $\epsilon_g$, and coupling factor $q$ as a function of $\Delta\nu$ is little affected by the chemical composition and it follows two different regimes depending on the evolutionary stage. On the subgiant branch, the models display a moderate core-envelope density contrast. The evolution of $\Delta \pi_1$, $\epsilon_p$, $\epsilon_g$, and $q$ thus significantly changes with the mass. Also, we demonstrate that, at fixed Z/X and with proper measurements of $\Delta\pi_1$ and $\Delta\nu$, we may unambiguously constrain the mass, radius and age of a subgiant star. Conversely, on the red-giant branch, the core-envelope density contrast becomes very large. Consequently, the evolution of $\epsilon_p$, $\epsilon_g$ and $q$ as a function of $\Delta\nu$ becomes independent of the mass. This is also true for $\Delta \pi_1$ in stars with masses $\lesssim 1.8M_\odot$ because of core electron degeneracy. This degeneracy is lifted for higher masses, again allowing for a precise measurement of the age. Overall, our computations qualitatively agree with past observed and theoretical studies.


title: Four-year Cosmology Large Angular Scale Surveyor (CLASS) Observations: On-sky Receiver Performance at 40, 90, 150, and 220 GHz Frequency Bands

abstract: The Cosmology Large Angular Scale Surveyor (CLASS) observes the polarized cosmic microwave background (CMB) over the angular scales of 1$^\circ \lesssim \theta \leq$ 90$^\circ$ with the aim of characterizing primordial gravitational waves and cosmic reionization. We report on the on-sky performance of the CLASS Q-band (40 GHz), W-band (90 GHz), and dichroic G-band (150/220 GHz) receivers that have been operational at the CLASS site in the Atacama desert since June 2016, May 2018, and September 2019, respectively. We show that the noise-equivalent power measured by the detectors matches the expected noise model based on on-sky optical loading and lab-measured detector parameters. Using Moon, Venus, and Jupiter observations, we obtain power-to-antenna-temperature calibrations and optical efficiencies for the telescopes. From the CMB survey data, we compute instantaneous array noise-equivalent-temperature sensitivities of 22, 19, 23, and 71 $\mathrm{\mu K}_\mathrm{cmb}\sqrt{\mathrm{s}}$ for the 40, 90, 150, and 220 GHz frequency bands, respectively. These noise temperatures refer to white noise amplitudes, which contribute to sky maps at all angular scales. Future papers will assess additional noise sources impacting larger angular scales.


title: A study on reliability of a $k$-out-of-$n$ system equipped with a cold standby component based on copula

abstract: A $k$-out-of-$n$ system consisting of $n$ exchangeable components equipped with a single cold standby component is considered. All the components including standby are assumed to be dependent and this dependence is modeled by a copula function. An exact expression for the reliability function of the considered system has been obtained. We also compute three different mean residual life functions of the system. Finally, some numerical examples are provided to illustrate the theoretical results. Our studies subsume some of the initial studies in the literature.


title: Port-Hamiltonian Neural Networks for Learning Explicit Time-Dependent Dynamical Systems

abstract: Accurately learning the temporal behavior of dynamical systems requires models with well-chosen learning biases. Recent innovations embed the Hamiltonian and Lagrangian formalisms into neural networks and demonstrate a significant improvement over other approaches in predicting trajectories of physical systems. These methods generally tackle autonomous systems that depend implicitly on time or systems for which a control signal is known apriori. Despite this success, many real world dynamical systems are non-autonomous, driven by time-dependent forces and experience energy dissipation. In this study, we address the challenge of learning from such non-autonomous systems by embedding the port-Hamiltonian formalism into neural networks, a versatile framework that can capture energy dissipation and time-dependent control forces. We show that the proposed \emph{port-Hamiltonian neural network} can efficiently learn the dynamics of nonlinear physical systems of practical interest and accurately recover the underlying stationary Hamiltonian, time-dependent force, and dissipative coefficient. A promising outcome of our network is its ability to learn and predict chaotic systems such as the Duffing equation, for which the trajectories are typically hard to learn.


title: Cosection localization and the Quot scheme $\mathrm{Quot}^{l}_{S}(\mathcal{E})$

abstract: Let $\mathcal{E}$ be a locally free sheaf of rank $r$ on a smooth projective surface $S$. The Quot scheme $\mathrm{Quot}^{l}_{S}(\mathcal{E})$ of length $l$ coherent sheaf quotients of $\mathcal{E}$ is a natural higher rank generalization of the Hilbert scheme of $l$ points of $S$. We study the virtual intersection theory of this scheme. If $C\subset S$ is a smooth canonical curve, we use cosection localization to show that the virtual fundamental class of $\mathrm{Quot}^{l}_{S}(\mathcal{E})$ is $(-1)^{l}$ times the fundamental class of the smooth subscheme $\mathrm{Quot}^{l}_{C}(\mathcal{E}\vert_{C})\subset\mathrm{Quot}^{l}_{S}(\mathcal{E})$. We then prove a structure theorem for virtual tautological integrals over $\mathrm{Quot}^{l}_{S}(\mathcal{E})$. From this we deduce, among other things, the equality of virtual Euler characteristics $\chi^{\mathrm{vir}}(\mathrm{Quot}^{l}_{S}(\mathcal{E}))=\chi^{\mathrm{vir}}(\mathrm{Quot}^{l}_{S}(\mathcal{O}^{\oplus r}))$.


title: Fast suppression of classification error in variational quantum circuits

abstract: Variational quantum circuits (VQCs) have shown great potential in near-term applications. However, the discriminative power of a VQC, in connection to its circuit architecture and depth, is not understood. To unleash the genuine discriminative power of a VQC, we propose a VQC system with the optimal classical post-processing -- maximum-likelihood estimation on measuring all VQC output qubits. Via extensive numerical simulations, we find that the error of VQC quantum data classification typically decay exponentially with the circuit depth, when the VQC architecture is extensive -- the number of gates does not shrink with the circuit depth. This fast error suppression ends at the saturation towards the ultimate Helstrom limit of quantum state discrimination. On the other hand, non-extensive VQCs such as quantum convolutional neural networks are sub-optimal and fail to achieve the Helstrom limit. To achieve the best performance for a given VQC, the optimal classical post-processing is crucial even for a binary classification problem. To simplify VQCs for near-term implementations, we find that utilizing the symmetry of the input properly can improve the performance, while oversimplification can lead to degradation.


title: Seeing and Believing: Evaluating the Trustworthiness of Twitter Users

abstract: Social networking and micro-blogging services, such as Twitter, play an important role in sharing digital information. Despite the popularity and usefulness of social media, there have been many instances where corrupted users found ways to abuse it, as for instance, through raising or lowering user's credibility. As a result, while social media facilitates an unprecedented ease of access to information, it also introduces a new challenge - that of ascertaining the credibility of shared information. Currently, there is no automated way of determining which news or users are credible and which are not. Hence, establishing a system that can measure the social media user's credibility has become an issue of great importance. Assigning a credibility score to a user has piqued the interest of not only the research community but also most of the big players on both sides - such as Facebook, on the side of industry, and political parties on the societal one. In this work, we created a model which, we hope, will ultimately facilitate and support the increase of trust in the social network communities. Our model collected data and analysed the behaviour of~50,000 politicians on Twitter. Influence score, based on several chosen features, was assigned to each evaluated user. Further, we classified the political Twitter users as either trusted or untrusted using random forest, multilayer perceptron, and support vector machine. An active learning model was used to classify any unlabelled ambiguous records from our dataset. Finally, to measure the performance of the proposed model, we used precision, recall, F1 score, and accuracy as the main evaluation metrics.


title: Continual Learning for Automated Audio Captioning Using The Learning Without Forgetting Approach

abstract: Automated audio captioning (AAC) is the task of automatically creating textual descriptions (i.e. captions) for the contents of a general audio signal. Most AAC methods are using existing datasets to optimize and/or evaluate upon. Given the limited information held by the AAC datasets, it is very likely that AAC methods learn only the information contained in the utilized datasets. In this paper we present a first approach for continuously adapting an AAC method to new information, using a continual learning method. In our scenario, a pre-optimized AAC method is used for some unseen general audio signals and can update its parameters in order to adapt to the new information, given a new reference caption. We evaluate our method using a freely available, pre-optimized AAC method and two freely available AAC datasets. We compare our proposed method with three scenarios, two of training on one of the datasets and evaluating on the other and a third of training on one dataset and fine-tuning on the other. Obtained results show that our method achieves a good balance between distilling new knowledge and not forgetting the previous one.


title: A sub-2% Distance to M31 from Photometrically Homogeneous Near-Infrared Cepheid Period-Luminosity Relations Measured with the Hubble Space Telescope

abstract: We present Period-Luminosity Relations (PLRs) for 55 Cepheids in M31 with periods ranging from 4 to 78 days observed with the Hubble Space Telescope (HST) using the same three-band photometric system recently used to calibrate their luminosities. Images were taken with the Wide Field Camera 3 in two optical filters (F555W and F814W) and one near-infrared filter (F160W) using the Drift And SHift (DASH) mode of operation to significantly reduce overheads and observe widely-separated Cepheids in a single orbit. We include additional F160W epochs for each Cepheid from the Panchromatic Hubble Andromeda Treasury (PHAT) and use light curves from the Panoramic Survey Telescope and Rapid Response System of the Andromeda galaxy (PAndromeda) project to determine mean magnitudes. Combined with a 1.28$\%$ absolute calibration of Cepheid PLRs in the Large Magellanic Cloud from Riess et al. (2019a) in the same three filters, we find a distance modulus to M31 of $\mu_0$ = 24.407 $\pm$ 0.032, corresponding to 761 $\pm$ 11 kpc and 1.49$\%$ uncertainty including all error sources, the most precise determination of its distance to date. We compare our results to past measurements using Cepheids and the Tip of the Red Giant Branch (TRGB). This study also provides the groundwork for turning M31 into a precision anchor galaxy in the cosmic distance ladder to measure the Hubble constant together with efforts to measure a fully geometric distance to M31.


title: A New Robust Multivariate Mode Estimator for Eye-tracking Calibration

abstract: We propose in this work a new method for estimating the main mode of multivariate distributions, with application to eye-tracking calibrations. When performing eye-tracking experiments with poorly cooperative subjects, such as infants or monkeys, the calibration data generally suffer from high contamination. Outliers are typically organized in clusters, corresponding to the time intervals when subjects were not looking at the calibration points. In this type of multimodal distributions, most central tendency measures fail at estimating the principal fixation coordinates (the first mode), resulting in errors and inaccuracies when mapping the gaze to the screen coordinates. Here, we developed a new algorithm to identify the first mode of multivariate distributions, named BRIL, which rely on recursive depth-based filtering. This novel approach was tested on artificial mixtures of Gaussian and Uniform distributions, and compared to existing methods (conventional depth medians, robust estimators of location and scatter, and clustering-based approaches). We obtained outstanding performances, even for distributions containing very high proportions of outliers, both grouped in clusters and randomly distributed. Finally, we demonstrate the strength of our method in a real-world scenario using experimental data from eye-tracking calibrations with Capuchin monkeys, especially for distributions where other algorithms typically lack accuracy.


title: Is attention to bounding boxes all you need for pedestrian action prediction?

abstract: The human driver is no longer the only one concerned with the complexity of the driving scenarios. Autonomous vehicles (AV) are similarly becoming involved in the process. Nowadays, the development of AVs in urban places raises essential safety concerns for vulnerable road users (VRUs) such as pedestrians. Therefore, to make the roads safer, it is critical to classify and predict the pedestrians' future behavior. In this paper, we present a framework based on multiple variations of the Transformer models able to infer predict the pedestrian street-crossing decision-making based on the dynamics of its initiated trajectory. We showed that using solely bounding boxes as input features can outperform the previous state-of-the-art results by reaching a prediction accuracy of 91\% and an F1-score of 0.83 on the PIE dataset. In addition, we introduced a large-size simulated dataset (CP2A) using CARLA for action prediction. Our model has similarly reached high accuracy (91\%) and F1-score (0.91) on this dataset. Interestingly, we showed that pre-training our Transformer model on the CP2A dataset and then fine-tuning it on the PIE dataset is beneficial for the action prediction task. Finally, our model's results are successfully supported by the "human attention to bounding boxes" experiment which we created to test humans ability for pedestrian action prediction without the need for environmental context. The code for the dataset and the models is available at: https://github.com/linaashaji/Action_Anticipation


title: First-Order Trotter Error from a Second-Order Perspective

abstract: Simulating quantum dynamics beyond the reach of classical computers is one of the main envisioned applications of quantum computers. The most promising quantum algorithms to this end in the near-term are the simplest, which use the Trotter formula and its higher-order variants to approximate the dynamics of interest. The approximation error of these algorithms is often poorly understood, even in the most basic cases, which are particularly relevant for experiments. Recent studies have reported anomalously low approximation error with unexpected scaling in such cases, which they attribute to quantum interference between the errors from different steps of the algorithm. Here we provide a simpler picture of these effects by relating the Trotter formula to its second-order variant. Our method generalizes state-of-the-art error bounds without the technical caveats of prior studies, and elucidates how each part of the total error arises from the underlying quantum circuit. We compare our bound to the true error numerically, and find a close match over many orders of magnitude in the simulation parameters. Our findings reduce the required circuit depth for the most basic quantum simulation algorithms, and illustrate a useful method for bounding simulation error more broadly.


title: Degenerate linear parabolic equations in divergence form on the upper half space

abstract: We study a class of second-order degenerate linear parabolic equations in divergence form in $(-\infty, T) \times \mathbb R^d_+$ with homogeneous Dirichlet boundary condition on $(-\infty, T) \times \partial \mathbb R^d_+$, where $\mathbb R^d_+ = \{x \in \mathbb R^d\,:\, x_d>0\}$ and $T\in {(-\infty, \infty]}$ is given. The coefficient matrices of the equations are the product of $\mu(x_d)$ and bounded uniformly elliptic matrices, where $\mu(x_d)$ behaves like $x_d^\alpha$ for some given $\alpha \in (0,2)$, which are degenerate on the boundary $\{x_d=0\}$ of the domain. Under a partially VMO assumption on the coefficients, we obtain the wellposedness and regularity of solutions in weighted Sobolev spaces. Our results can be readily extended to systems.


title: Pilot Study Suggests Online Media Literacy Programming Reduces Belief in False News in Indonesia

abstract: Amidst the threat of digital misinformation, we offer a pilot study regarding the efficacy of an online social media literacy campaign aimed at empowering individuals in Indonesia with skills to help them identify misinformation. We found that users who engaged with our online training materials and educational videos were more likely to identify misinformation than those in our control group (total $N$=1000). Given the promising results of our preliminary study, we plan to expand efforts in this area, and build upon lessons learned from this pilot study.


title: Testing Surrogate-Based Optimization with the Fortified Branin-Hoo Extended to Four Dimensions

abstract: Some popular functions used to test global optimization algorithms have multiple local optima, all with the same value, making them all global optima. It is easy to make them more challenging by fortifying them via adding a localized bump at the location of one of the optima. In previous work the authors illustrated this for the Branin-Hoo function and the popular differential evolution algorithm, showing that the fortified Branin-Hoo required an order of magnitude more function evaluations. This paper examines the effect of fortifying the Branin-Hoo function on surrogate-based optimization, which usually proceeds by adaptive sampling. Two algorithms are considered. The EGO algorithm, which is based on a Gaussian process (GP) and an algorithm based on radial basis functions (RBF). EGO is found to be more frugal in terms of the number of required function evaluations required to identify the correct basin, but it is expensive to run on a desktop, limiting the number of times the runs could be repeated to establish sound statistics on the number of required function evaluations. The RBF algorithm was cheaper to run, providing more sound statistics on performance. A four-dimensional version of the Branin-Hoo function was introduced in order to assess the effect of dimensionality. It was found that the difference between the ordinary function and the fortified one was much more pronounced for the four-dimensional function compared to the two dimensional one.


title: The strong converse exponent of discriminating infinite-dimensional quantum states

abstract: The sandwiched R\'enyi divergences of two finite-dimensional density operators quantify their asymptotic distinguishability in the strong converse domain. This establishes the sandwiched R\'enyi divergences as the operationally relevant ones among the infinitely many quantum extensions of the classical R\'enyi divergences for R\'enyi parameter $\alpha>1$. The known proof of this goes by showing that the sandwiched R\'enyi divergence coincides with the regularized measured R\'enyi divergence, which in turn is proved by asymptotic pinching, a fundamentally finite-dimensional technique. Thus, while the notion of the sandwiched R\'enyi divergences was extended recently to density operators on an infinite-dimensional Hilbert space (even for states of a von Neumann algebra), these quantities were so far lacking an operational interpretation similar to the finite-dimensional case, and it has also been open whether they coincide with the regularized measured R\'enyi divergences. In this paper we fill this gap by answering both questions in the positive for density operators on an infinite-dimensional Hilbert space, using a simple finite-dimensional approximation technique. We also initiate the study of the sandwiched R\'enyi divergences, and the related problem of the strong converse exponent, for pairs of positive semi-definite operators that are not necessarily trace-class. This is interesting from the purely mathematical point of view of extending the concept of R\'enyi (and other) divergences to settings beyond the standard one of positive trace-class operators (positive normal functionals in the von Neumann algebra setting). In this spirit, we also discuss the definition and some properties of the more general family of R\'enyi $(\alpha,z)$-divergences of positive semi-definite operators on an infinite-dimensional separable Hilbert space.


title: CCVS: Context-aware Controllable Video Synthesis

abstract: This presentation introduces a self-supervised learning approach to the synthesis of new video clips from old ones, with several new key elements for improved spatial resolution and realism: It conditions the synthesis process on contextual information for temporal continuity and ancillary information for fine control. The prediction model is doubly autoregressive, in the latent space of an autoencoder for forecasting, and in image space for updating contextual information, which is also used to enforce spatio-temporal consistency through a learnable optical flow module. Adversarial training of the autoencoder in the appearance and temporal domains is used to further improve the realism of its output. A quantizer inserted between the encoder and the transformer in charge of forecasting future frames in latent space (and its inverse inserted between the transformer and the decoder) adds even more flexibility by affording simple mechanisms for handling multimodal ancillary information for controlling the synthesis process (eg, a few sample frames, an audio track, a trajectory in image space) and taking into account the intrinsically uncertain nature of the future by allowing multiple predictions. Experiments with an implementation of the proposed approach give very good qualitative and quantitative results on multiple tasks and standard benchmarks.


title: Towards a Benchmark Set for Program Repair Based on Partial Fixes

abstract: Software bugs significantly contribute to software cost and increase the risk of system malfunctioning. In recent years, many automated program-repair approaches have been proposed to automatically fix undesired program behavior. Despite of their great success, specific problems such as fixing bugs with partial fixes still remain unresolved. A partial fix to a known software issue is a programmer's failed attempt to fix the issue the first time. Even though it fails, this fix attempt still conveys important information such as the suspicious software region and the bug type. In this work we do not propose an approach for program repair with partial fixes, but instead answer a preliminary question: Do partial fixes occur often enough, in general, to be relevant for the research area of automated program repair? We crawled 1500 open-source C repositories on GitHub for partial fixes. The result is a benchmark set of 2204 benchmark tasks for automated program repair based on partial fixes. The benchmark set is available open source and open to further contributions and improvement.


title: Representation Consolidation for Training Expert Students

abstract: Traditionally, distillation has been used to train a student model to emulate the input/output functionality of a teacher. A more useful goal than emulation, yet under-explored, is for the student to learn feature representations that transfer well to future tasks. However, we observe that standard distillation of task-specific teachers actually *reduces* the transferability of student representations to downstream tasks. We show that a multi-head, multi-task distillation method using an unlabeled proxy dataset and a generalist teacher is sufficient to consolidate representations from task-specific teacher(s) and improve downstream performance, outperforming the teacher(s) and the strong baseline of ImageNet pretrained features. Our method can also combine the representational knowledge of multiple teachers trained on one or multiple domains into a single model, whose representation is improved on all teachers' domain(s).


title: The Fourier formalism for relativistic axion-photon conversion, with astrophysical applications

abstract: We study the weak mixing of photons and relativistic axion-like particles (axions) in plasmas with background magnetic fields, ${\bf B}$. We show that, to leading order in the axion-photon coupling, the conversion probability, $P_{\gamma \to a}$, is given by the one-dimensional power spectrum of the magnetic field components perpendicular to the particle trajectory. Equivalently, we express $P_{\gamma \to a}$ as the Fourier transform of the magnetic field autocorrelation function, and establish a dictionary between properties of the real-space magnetic field and the energy-dependent conversion probability. For axions more massive than the plasma frequency, ($m_a>\omega_{\rm pl}$), we use this formalism to analytically solve the problem of perturbative axion-photon mixing in a general magnetic field. In the general case where $\omega_{\rm pl}/m_a$ varies arbitrarily along the trajectory, we show that a naive application of the standard formalism for 'resonant' conversion can give highly inaccurate results, and that a careful calculation generically gives non-resonant contributions at least as large as the resonant contribution. Furthermore, we demonstrate how techniques based on the Fast Fourier Transform provide a new, highly efficient numerical method for calculating axion-photon mixing. We briefly discuss magnetic field modelling in galaxy clusters in the light of our results and argue, in particular, that a recently proposed 'regular' model used for studying axion-photon mixing (specifically applied to the Perseus cluster) is inconsistent with observations. Our formalism suggest new methods to search for imprints of axions, and will be important for spectrographs with percent level sensitivity, which includes existing X-ray observations by Chandra as well as the upcoming Athena mission.


title: Cosmic shear beyond 2-point statistics: Accounting for galaxy intrinsic alignment with projected tidal fields

abstract: Developing analysis pipelines based on statistics beyond two-point functions is critical for extracting a maximal amount of cosmological information from current and upcoming weak lensing surveys. In this paper, we study the impact of the intrinsic alignment of galaxies (IA) on three promising probes measured from aperture mass maps -- the lensing peaks, minima and full PDF, in comparison and in combination with the shear two-point correlation functions ($\gamma$-2PCFs). Our two-dimensional IA infusion method converts the light-cone-projected mass sheets into projected tidal tensors, which are then linearly coupled to an intrinsic ellipticity component with a strength controlled by the coupling parameter $A_{\rm IA}$. We validate our method with the $\gamma$-2PCFs statistics, recovering well the analytical calculations from the linear alignment model of Bridle \& King in a full tomographic setting, and for different $A_{\rm IA}$ values. We next use our method to infuse at the galaxy catalogue level a non-linear IA model that includes the density-weighting term introduced in \citet{Blazek2015}, and compute the impact on the three aperture mass map statistics. We find that large \snr peaks are maximally affected, with deviations reaching 30\% (10\%) for a {\it Euclid}-like (KiDS-like) survey. Modelling the signal in a $w$CDM cosmology universe with $N$-body simulations, we forecast the cosmological bias caused by unmodelled IA for 100 deg$^2$ of {\it Euclid}-like data, finding very large offsets in $w_0$ (5-10$\sigma_{\rm stat}$), $\Omega_{\rm m}$ (4-6$\sigma_{\rm stat}$), and $S_8 \equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3}$ ($\sim$3$\sigma_{\rm stat}$). The method presented in this paper offers a compelling avenue to account for IA in beyond-two-point weak lensing statistics, with a flexibility comparable to that of current $\gamma$-2PCFs IA analytical models.


title: Local Control of Supercurrent Density in Epitaxial Planar Josephson Junctions

abstract: The critical current response to an applied out-of-plane magnetic field in a Josephson junction provides insight into the uniformity of its current distribution. In Josephson junctions with semiconducting weak links, the carrier density and, therefore the overall current distribution could be modified electrostatically via metallic gates. Here, we show local control of the current distribution in an epitaxial Al-InAs Josephson junction equipped with five mini-gates. We demonstrate that not only can the junction width be electrostatically defined but also we can locally adjust the current profile to form superconducting quantum interference devices. Our studies show enhanced edge conduction in such long junctions, which can be eliminated by mini-gates to create a uniform current distribution.


title: Exploration of doped quantum magnets with ultracold atoms

abstract: In the last decade, quantum simulators, and in particular cold atoms in optical lattices, have emerged as a valuable tool to study strongly correlated quantum matter. These experiments are now reaching regimes that are numerically difficult or impossible to access. In particular they have started to fulfill a promise which has contributed significantly to defining and shaping the field of cold atom quantum simulations, namely the exploration of doped and frustrated quantum magnets and the search for the origins of high-temperature superconductivity in the fermionic Hubbard model. Despite many future challenges lying ahead, such as the need to further lower the experimentally accessible temperatures, remarkable studies have already emerged. Among them, spin-charge separation in one-dimensional systems has been demonstrated, extended-range antiferromagnetism in two-dimensional systems has been observed, connections to modern day large-scale numerical simulations were made, and unprecedented comparisons with microscopic trial wavefunctions have been carried out at finite doping. In many regards, the field has acquired new realms, putting old ideas to a new test and producing new insights and inspiration for the next generation of physicists. In the first part of this paper, we review the results achieved in cold atom realizations of the Fermi-Hubbard model in recent years. In the second part of this paper, with the stage set and the current state of the field in mind, we propose a new direction for cold atoms to explore: namely mixed-dimensional bilayer systems, where the charge motion is restricted to individual layers which remain coupled through spin-exchange. We propose a novel, strong pairing mechanism in these systems, which puts the formation of hole pairs at experimentally accessible, elevated temperatures within reach.

categories: cond-mat.quant-gas cond-mat.str-el cond-mat.supr-con quant-ph


title: Quantum nature of molecular vibrational quenching: Water - molecular hydrogen collisions

abstract: Rates of conversions of molecular internal energy to and from kinetic energy by means of molecular collision allows to compute collisional line shapes and transport properties of gases. Knowledge of ro-vibrational quenching rates is necessary to connect spectral observations to physical properties of warm astrophysical gasses, including exo-atmospheres. For a system of paramount importance in this context, the vibrational bending mode quenching of H2O by H2, we show here that exchange of vibrational to rotational and kinetic energy remains a quantum process, despite the large numbers of quantum levels involved and the large vibrational energy transfer. The excitation of the quantized rotor of the projectile is by far the most effective ro-vibrational quenching path of water. To do so, we use a fully quantum first principle computation, potential and dynamics, converging it at all stages, in a full coupled channel formalisms. We present here rates for the quenching of the first bendingmode of ortho-H2O by ortho H2, up to 500K, in a fully converged coupled channels formalism.

categories: physics.chem-ph astro-ph.GA astro-ph.SR quant-ph


title: Desiderata for Explainable AI in statistical production systems of the European Central Bank

abstract: Explainable AI constitutes a fundamental step towards establishing fairness and addressing bias in algorithmic decision-making. Despite the large body of work on the topic, the benefit of solutions is mostly evaluated from a conceptual or theoretical point of view and the usefulness for real-world use cases remains uncertain. In this work, we aim to state clear user-centric desiderata for explainable AI reflecting common explainability needs experienced in statistical production systems of the European Central Bank. We link the desiderata to archetypical user roles and give examples of techniques and methods which can be used to address the user's needs. To this end, we provide two concrete use cases from the domain of statistical data production in central banks: the detection of outliers in the Centralised Securities Database and the data-driven identification of data quality checks for the Supervisory Banking data system.


title: An Open-channel Microfluidic Membrane Device for In situ Hyperspectral mapping of enzymatic cellulose hydrolysis

abstract: Synchrotron infrared hyperspectral microscopy is a label-free and non-invasive technique well suited for imaging of chemical events in situ. It can track the spatial and temporal distributions of molecules of interests in a specimen in its native state by the molecule's characteristic vibrational modes. Despite tremendous progress made in recent years, IR hyperspectral imaging of chemical events in biomaterials in liquids remains challenging because of the demanding requirements on environmental control and strong infrared absorption of water. Here we report a multi-phase capillary-driven membrane device for label-free and real-time investigation of enzymatic deconstruction of algal cellulose purified from Cladophora aegagropila.


title: Quantum computations (course of lectures)

abstract: This course of lectures has been taught for several years at the Lomonosov Moscow State University; its modified version in 2021 is read in the Zhejiang University (Hangzhou), in the framework of summer school on quantum computing. The course is devoted to a new type of computations based on quantum mechanics. Quantum computations are fundamentally different from classical ones in that they occur in the space of so-called quantum states, and not in ordinary binary strings. Quantum computing is a real process in which the mathematical description is inextricably linked with quantum physics. Various forms of quantum computing are considered: the Feynman gate model, fermionic and adiabatic computations. A class of problems is described in which quantum computing is not only more efficient than classical ones, but also cannot be replaced by them. The course is designed for students of physics and mathematics and natural science specialties as well as all those interested in this subject. It requires familiarity with the basics of linear algebra and mathematical analysis in the first two courses of the university.


title: Primordial black holes formation and secondary gravitational waves in nonminimal derivative coupling inflation

abstract: We study the possibility of the Primordial Black Holes (PBHs) formation with the aim of finding a considerable fraction of Dark Matter (DM), using the gravitationally enhanced friction mechanism which arises from a nonminimal derivative coupling between the scalar field and the gravity. Assuming the nonminimal coupling parameter as a special function of the scalar field and considering the potential of natural inflation, we find three parameter sets that produce a period of ultra slow-roll inflation. This leads to sufficiently large enhancement in the curvature power spectra to form PBHs. We show that under the gravitationally enhanced friction mechanism, PBHs with a mass around ${\cal O}\big(10^{-12}\big)M_\odot$ can constitute around $96\%$ of the total DM and so this class of PBHs can be taken as a great candidate for DM. We further study the secondary Gravitational Waves (GWs) in our setting and show that our model predicts the peak of the present fractional energy density as $\Omega_{GW0} \sim 10^{-8}$ at different frequencies for all the three parameter sets. This value lies well inside the sensitivity region of some GWs detectors at some frequencies, and therefore the observational compatibility of our model can be appraised by the data from these detectors.


title: 3D Topological Quantum Computing

abstract: In this paper we will present some ideas to use 3D topology for quantum computing extending ideas from a previous paper. Topological quantum computing used \textquotedblleft knotted\textquotedblright{} quantum states of topological phases of matter, called anyons. But anyons are connected with surface topology. But surfaces have (usually) abelian fundamental groups and therefore one needs non-abelian anyons to use it for quantum computing. But usual materials are 3D objects which can admit more complicated topologies. Here, complements of knots do play a prominent role and are in principle the main parts to understand 3-manifold topology. For that purpose, we will construct a quantum system on the complements of a knot in the 3-sphere (see arXiv:2102.04452 for previous work). The whole system is designed as knotted superconductor where every crossing is a Josephson junction and the qubit is realized as flux qubit. We discuss the properties of this systems in particular the fluxion quantization by using the A-polynomial of the knot. Furthermore we showed that 2-qubit operations can be realized by linked (knotted) superconductors again coupled via a Josephson junction.


title: Rotation and Lithium Confirmation of a 500 Parsec Halo for the Open Cluster NGC 2516

abstract: Recent analyses of the Gaia data have identified diffuse stellar populations surrounding nearby open clusters. It is important to verify that these "halos", "tails", and "strings" are of similar ages and compositions as stars in the denser part of the cluster. We present an analysis of NGC 2516 ($\approx$150 Myr), which has a classical tidal radius of 10 pc and an apparent halo of stars spanning 500 pc ($20^\circ$ on-sky). Combining photometry from Gaia, rotation periods from TESS, and lithium measurements from Gaia-ESO and GALAH, we find that the halo of NGC 2516 is the same age as the cluster's core. Two thirds of kinematically selected halo members out to 250 pc from the cluster center have rotation periods consistent with a gyrochronological age of 150 Myr. A comparison sample of field stars shows no such trend. The lithium abundances of stars in the halo are higher than in the field, and are correlated with the stellar rotation rate and binarity fraction, as has been noted in other young open clusters. Broadly speaking, this work supports a new paradigm wherein the halos of open clusters are often more populous than their cores. We highlight implications for spectroscopic survey targeting, open cluster dispersal, and planet searches around young stars.


title: Ultracold spin-balanced fermionic quantum liquids with renormalized $P$-wave interactions

abstract: We consider a spin-balanced degenerate gas of spin-1/2 fermions whose dynamics is governed by low-energy $P$-wave interactions, characterized by the scattering volume $a_1$ and effective momentum $r_1$. The energy per particle $\bar{\cal{E}}$ in the many-body system is calculated by resumming the ladder diagrams comprising both particle-particle and hole-hole intermediate states, following the novel advances recently developed by us in Ann.Phys. 437,168741(2022). This allows to obtain a renormalized result for $\bar{\cal{E}}$ within generic cutoff regularization schemes, with $\bar{\cal{E}}$ directly expressed in terms of the scattering parameters $a_1$ and $r_1$, once the cut off is sent to infinity. The whole set of possible values of $a_1$ and $r_1$ is explored for the first time in the literature looking for minima in the energy per particle with $\bar{\cal{E}}$ given as described. They are actually found, but a further inspection reveals that the associated scattering parameters give rise to resonance poles in the complex momentum-plane with positive imaginary part, which is at odds with the Hermiticity of the Hamiltonian. We also determine that these conflictive poles, with a pole-position momentum that is smaller in absolute value than the Fermi momentum of the system, clearly impact the calculation of $\bar{\cal{E}}$. As a result, we conclude that unpolarized spin-1/2 fermionic normal matter interacting in $P$-wave is not stable. We also study three universal parameters around the unitary limit. Finally, the whole set of values for the parameters $a_1$, $r_1$ is characterized according to whether they give rise to unallowed poles and, if so, by attending to their pole positions relative to the Fermi momentum of the system explored.

categories: cond-mat.quant-gas hep-ph nucl-th quant-ph


title: Long-Range Vector Models at Large N

abstract: We calculate various CFT data for the $O(N)$ vector model with the long-range interaction, working at the next-to-leading order in the $1/N$ expansion. Our results provide additional evidence for the existence of conformal symmetry at the long-range fixed point, as well as the continuity of the CFT data at the long-range to short-range crossover point $s_\star$ of the exponent parameter $s$. We also develop the $N>1$ generalization of the recently proposed IR duality between the long-range and the deformed short-range models, providing further evidence for its non-perturbative validity in the entire region $d/2<s<s_\star$.


title: Two steps forward and one step sideways: the propagation of relativistic jets in realistic binary neutron star merger ejecta

abstract: The association of GRB170817A with GW170817 has confirmed the long-standing hypothesis that binary neutron star (BNS) mergers are the progenitors of at least some short gamma-ray bursts (SGRBs). This connection has ushered in an era in which broadband observations of SGRBs, together with measurements of the time delay between the gravitational waves and the electromagnetic radiation, allow to probe the properties of the emitting outflow and its engine to an unprecedented detail. Since the structure of the radiating outflow is molded by the interaction of a relativistic jet with the binary ejecta, it is of paramount importance to study the system in a realistic setting. Here we present a three-dimensional hydrodynamic simulation of a relativistic jet propagating in the ejecta of a BNS merger, which were computed with a general relativistic magnetohydrodynamic simulation. We find that the jet's centroid oscillates around the axis of the system, due to inhomogeneities encountered in the propagation. These oscillations allow the jet to find the path of least resistance and travel faster than an identical jet in smooth ejecta. In our setup the breakout time is ~0.6 sec, comparable to the expected central engine duration in SGRBs and possibly a non-negligible fraction of the total delay between the gravitational and gamma-ray signals. Our simulation also shows that energy is carried in roughly equal amounts by the jet and by the cocoon, and that about 20 per cent of the injected energy is transferred to the ejecta via mechanical work.


title: The reach of next-to-leading-order perturbation theory for the matter bispectrum

abstract: We provide a comparison between the matter bispectrum derived with different flavours of perturbation theory at next-to-leading order and measurements from an unprecedentedly large suite of $N$-body simulations. We use the $\chi^2$ goodness-of-fit test to determine the range of accuracy of the models as a function of the volume covered by subsets of the simulations. We find that models based on the effective-field-theory (EFT) approach have the largest reach, standard perturbation theory has the shortest, and `classical' resummed schemes lie in between. The gain from EFT, however, is less than in previous studies. We show that the estimated range of accuracy of the EFT predictions is heavily influenced by the procedure adopted to fit the amplitude of the counterterms. For the volumes probed by galaxy redshift surveys, our results indicate that it is advantageous to set three counterterms of the EFT bispectrum to zero and measure the fourth from the power spectrum. We also find that large fluctuations in the estimated reach occur between different realisations. We conclude that it is difficult to unequivocally define a range of accuracy for the models containing free parameters. Finally, we approximately account for systematic effects introduced by the $N$-body technique either in terms of a scale- and shape-dependent bias or by boosting the statistical error bars of the measurements (as routinely done in the literature). We find that the latter approach artificially inflates the reach of EFT models due to the presence of tunable parameters.


title: Exploring the Sgr-Milky-Way-disc interaction using high resolution N-body simulations

abstract: The ongoing merger of the Sagittarius (Sgr) dwarf galaxy with the Milky Way is believed to strongly affect the dynamics of the Milky Way's disc. We present a suite of 13 $N$-body simulations, with 500 million to 1 billion particles, modelling the interaction between the Sagittarius dwarf galaxy (Sgr) and the Galactic disc. To quantify the perturbation to the disc's structure and dynamics in the simulation, we compute the number count asymmetry and the mean vertical velocity in a solar-neighbourhood-like volume. We find that overall the trends in the simulations match those seen in a simple one-dimensional model of the interaction. We explore the effects of changing the mass model of Sgr, the orbital kinematics of Sgr, and the mass of the Milky Way halo. We find that none of the simulations match the observations of the vertical perturbation using Gaia Data Release 2. In the simulation which is the most similar, we find that the final mass of Sgr far exceeds the observed mass of the Sgr remnant, the asymmetry wavelength is too large, and the shape of the asymmetry doesn't match past $z \approx 0.7$ kpc. We therefore conclude that our simulations support the conclusion that Sgr alone could not have caused the observed perturbation to the solar neighbourhood.


title: Sgr A* X-ray flares from non-thermal particle acceleration in a magnetically arrested disc

abstract: Sgr A* exhibits flares in the near-infrared and X-ray bands, with the luminosity in these bands increasing by factors of $10-100$ for $\approx 60$ minutes. One of the models proposed to explain these flares is synchrotron emission of non-thermal particles accelerated by magnetic reconnection events in the accretion flow. We use the results from PIC simulations of magnetic reconnection to post-process 3D two-temperature GRMHD simulations of a magnetically arrested disc (MAD). We identify current sheets, retrieve their properties, estimate their potential to accelerate non-thermal particles and compute the expected non-thermal synchrotron emission. We find that the flux eruptions of MADs can provide suitable conditions for accelerating non-thermal particles to energies $\gamma_e \lesssim 10^6$ and producing simultaneous X-ray and near-infrared flares. For a suitable choice of current-sheet parameters and a simpified synchrotron cooling prescription, the model can simultaneously reproduce the quiescent and flaring X-ray luminosities as well as the X-ray spectral shape. While the near-infrared flares are mainly due to an increase in the temperature near the black hole during the MAD flux eruptions, the X-ray emission comes from narrow current sheets bordering highly magnetized, low-density regions near the black hole and equatorial current sheets where the flux on the black hole reconnects. As a result, not all infrared flares are accompanied by X-ray ones. The non-thermal flaring emission can extend to very hard ($\lesssim 100$ keV) X-ray energies.


title: Unveiling the star formation history of the Upper Scorpius association through its kinematics

abstract: Stellar associations can be discerned as overdensities of sources not only in the physical space but also in the velocity space. The common motion of their members, gradually eroded by the galactic tidal field, is partially reminiscent of the initial kinematic structure. Using recent data from Gaia EDR3, combined with radial velocities from GALAH and APOGEE, we traced back the present positions of stars belonging to Upper Scorpius, a subgroup of Scorpius-Centaurus, the nearest OB association. About one half of the subgroup (the "clustered" population) appears composed of many smaller entities, which were in a more compact configuration in the past. The presence of a kinematic duality is reflected into an age spread between this younger clustered population and an older diffuse population, in turn confirmed by a different fraction $f_D$ of disc-bearing stars ($f_D = 0.24\pm0.02$ vs $f_D = 0.10\pm 0.01$). Star formation in Upper Scorpius appears to have lasted more than 10 Myr and proceeded in small groups that, after a few Myr, dissolve in the field of the older population but retain for some time memory of their initial structure. The difference of ages inferred through isochrones and kinematics, in this regard, could provide a powerful tool to quantify the timescale of gas removal.


title: Spin-based removal of instrumental systematics in 21cm intensity mapping surveys

abstract: Upcoming cosmological intensity mapping surveys will open new windows on the Universe, but they must first overcome a number of significant systematic effects, including polarization leakage. We present a formalism that uses scan strategy information to model the effect of different instrumental systematics on the recovered cosmological intensity signal for `single-dish' (autocorrelation) surveys. This modelling classifies different systematics according to their spin symmetry, making it particularly relevant for dealing with polarization leakage. We show how to use this formalism to calculate the expected contamination from different systematics as a function of the scanning strategy. Most importantly, we show how systematics can be disentangled from the intensity signal based on their spin properties via map-making. We illustrate this, using a set of toy models, for some simple instrumental systematics, demonstrating the ability to significantly reduce the contamination to the observed intensity signal. Crucially, unlike existing foreground removal techniques, this approach works for signals that are non-smooth in frequency, e.g. polarized foregrounds. These map-making approaches are simple to apply and represent an orthogonal and complementary approach to existing techniques for removing systematics from upcoming 21cm intensity mapping surveys.


title: Renormalizable Models of Flavor-Specific Scalars

abstract: New light singlet scalars with flavor-specific couplings represent a phenomenologically distinctive and flavor-safe alternative to the well-studied possibility of Higgs-portal scalars. However, in contrast to the Higgs portal, flavor-specific couplings require an ultraviolet completion involving new heavy states charged under the Standard Model gauge symmetries, leading to a host of additional novel phenomena. Focusing for concreteness on a scenario with up quark-specific couplings, we investigate two simple renormalizable completions, one with an additional vector-like quark and another featuring an extra scalar doublet. We consider the implications of naturalness, flavor- and CP-violation, electroweak precision observables, and direct searches for the new states at the LHC. These bounds, while being model-dependent, are shown to probe interesting regions in the parameter space of the scalar mass and its low-energy effective coupling, complementing the essential phenomenology of the low-energy effective theory at a variety of low and medium energy experiments.


title: X-ray bounds on cooling, composition, and magnetic field of the Cassiopeia A neutron star and young central compact objects

abstract: We present analysis of multiple Chandra and XMM-Newton spectra, separated by 9-19 years, of four of the youngest central compact objects (CCOs) with ages < 2500 yr: CXOU J232327.9+584842 (Cassiopeia A), CXOU J160103.1-513353 (G330.2+1.0), 1WGA J1713.4-3949 (G347.3-0.5), and XMMU J172054.5-372652 (G350.1-0.3). By fitting these spectra with thermal models, we attempt to constrain each CCO's long-term cooling rate, composition, and magnetic field. For the CCO in Cassiopeia A, 14 measurements over 19 years indicate a decreasing temperature at a ten-year rate of 2.2+/-0.2 or 2.8+/-0.3 percent (1sigma error) for a constant or changing X-ray absorption, respectively. We obtain cooling rate upper limits of 17 percent for CXOU J160103.1-513353 and 6 percent for XMMU J172054.5-372652. For the oldest CCO, 1WGA J1713.4-3949, its temperature seems to have increased by 4+/-2 percent over a ten year period. Assuming each CCO's preferred distance and an emission area that is a large fraction of the total stellar surface, a non-magnetic carbon atmosphere spectrum is a good fit to spectra of all four CCOs. If distances are larger and emission areas are somewhat smaller, then equally good spectral fits are obtained using a hydrogen atmosphere with B <= 7x10^10 G or B >= 10^12 G for CXOU J160103.1-513353, B <= 10^10 G or B >= 10^12 G for XMMU J172054.5-372652, and non-magnetic hydrogen atmosphere for 1WGA J1713.4-3949. In a unified picture of CCO evolution, our results suggest most CCOs, and hence a sizable fraction of young neutron stars, have a surface magnetic field that is low early in their life but builds up over several thousand years.


title: Odd-parity stability of black holes in Einstein-Aether gravity

abstract: In Einstein-Aether theory, we study the stability of black holes against odd-parity perturbations on a spherically symmetric and static background. For odd-parity modes, there are two dynamical degrees of freedom arising from the tensor gravitational sector and Aether vector field. We derive general conditions under which neither ghosts nor Laplacian instabilities are present for these dynamical fields. We apply these results to concrete black hole solutions known in the literature and show that some of those solutions can be excluded by the violation of stability conditions. The exact Schwarzschild solution present for $c_{13} = c_{14} = 0$, where $c_i$'s are the four coupling constants of the theory with $c_{ij}=c_i + c_j$, is prone to Laplacian instabilities along the angular direction throughout the horizon exterior. However, we find that the odd-parity instability of high radial and angular momentum modes is absent for black hole solutions with $c_{13} = c_4 = 0$ and $c_1 \geq 0$.


title: Using saturated count models for user-friendly synthesis of categorical data

abstract: Over the past three decades, synthetic data methods for statistical disclosure control have continually evolved, but mainly within the domain of survey data sets. There are certain characteristics of administrative databases, such as their size, which present challenges from a synthesis perspective and require special attention. This paper, through the fitting of saturated count models, presents a synthesis method that is suitable for administrative databases that is tuned by two parameters. The method allows large categorical data sets to be synthesized quickly and allows risk and utility metrics to be satisfied a priori, that is, prior to synthetic data generation. The paper explores how the flexibility afforded by two-parameter count models (the negative binomial and Poisson-inverse Gaussian) can be utilised to protect respondents' - especially uniques' - privacy in synthetic data. Finally, an empirical example is carried out through the synthesis of a database which can be viewed as a good substitute to the English School Census.


title: Studying Bioluminescence Flashes with the ANTARES Deep Sea Neutrino Telescope

abstract: We develop a novel technique to exploit the extensive data sets provided by underwater neutrino telescopes to gain information on bioluminescence in the deep sea. The passive nature of the telescopes gives us the unique opportunity to infer information on bioluminescent organisms without actively interfering with them. We propose a statistical method that allows us to reconstruct the light emission of individual organisms, as well as their location and movement. A mathematical model is built to describe the measurement process of underwater neutrino telescopes and the signal generation of the biological organisms. The Metric Gaussian Variational Inference algorithm is used to reconstruct the model parameters using photon counts recorded by the neutrino detectors. We apply this method to synthetic data sets and data collected by the ANTARES neutrino telescope. The telescope is located 40 km off the French coast and fixed to the sea floor at a depth of 2475 m. The runs with synthetic data reveal that we can reliably model the emitted bioluminescent flashes of the organisms. Furthermore, we find that the spatial resolution of the localization of light sources highly depends on the configuration of the telescope. Precise measurements of the efficiencies of the detectors and the attenuation length of the water are crucial to reconstruct the light emission. Finally, the application to ANTARES data reveals the first precise localizations of bioluminescent organisms using neutrino telescope data.


title: Frequency analysis of Delta Scuti stars towards the Galactic bulge

abstract: We have performed a frequency analysis of 10,092 Delta Scuti-type stars detected in the fourth phase of the Optical Gravitational Lensing Experiment (OGLE) towards the Galactic bulge, which is the most numerous homogeneous sample of Delta Scuti stars observed so far. The main goal was to search for stars pulsating in at least two radial modes simultaneously. We have found 3083 candidates for such stars, which is the largest set obtained to date. Among them, 2655 stars pulsate in two radial modes, 414 stars pulsate in three radial modes, and 14 stars pulsate in four radial modes at the same time. We report the identification of 221 Delta Scuti stars pulsating in the fundamental mode, first overtone, and third overtone simultaneously. We show the most populated Petersen and Bailey diagrams and discuss statistical properties of the identified frequencies based on this numerous sample. Additionally, we present theoretical predictions of period ratios for Delta Scuti stars pulsating in overtones from the fourth to the seventh.