title: Crossing phantom divide in $f(Q)$ gravity

abstract: We investigate the possibility of crossing a phantom divide line in the extension of symmetric teleparallel gravity or the $f(Q)$ gravity, where $Q$ is the non-metricity. We study the cosmic evolution of the effective equation of state parameter for dark energy considering exponential, logarithmic, and combined $f(Q)$ theories. Moreover, the exponential model behaves like the $\Lambda$CDM at high redshifts before deviating to $\omega_{eff}<-1$ or $\omega_{eff}>-1$, respectively, depending on the value of model parameter. It also approaches a de-sitter phase asymptotically. However, the crossing of the phantom divide line, i.e., $\omega= -1$, is realized in the combined $f(Q)$ theory. Furthermore, statefinder diagnostics are studied in order to differentiate between several dark energy models. To ensure the three model's stability, we employ the stability analysis using linear perturbations. We demonstrate how to reassemble $f(Q)$ via a numerical inversion approach based on existing observational constraints on cosmographic parameters and the potential of bridging the phantom divide in the resulting model. It explicitly demonstrates that future crossings of the phantom dividing line are a generic feature of feasible $f(Q)$ gravity models.


title: PAVI: Plate-Amortized Variational Inference

abstract: Given some observed data and a probabilistic generative model, Bayesian inference aims at obtaining the distribution of a model's latent parameters that could have yielded the data. This task is challenging for large population studies where thousands of measurements are performed over a cohort of hundreds of subjects, resulting in a massive latent parameter space. This large cardinality renders off-the-shelf Variational Inference (VI) computationally impractical. In this work, we design structured VI families that can efficiently tackle large population studies. To this end, our main idea is to share the parameterization and learning across the different i.i.d. variables in a generative model -symbolized by the model's plates. We name this concept plate amortization, and illustrate the powerful synergies it entitles, resulting in expressive, parsimoniously parameterized and orders of magnitude faster to train large scale hierarchical variational distributions. We illustrate the practical utility of PAVI through a challenging Neuroimaging example featuring a million latent parameters, demonstrating a significant step towards scalable and expressive Variational Inference.


title: The Z3RO Family of Precoders Cancelling Nonlinear Power Amplification Distortion in Large Array Systems

abstract: Large array systems use a massive number of antenna elements and clever precoder designs to achieve an array gain at the user location. These precoders require linear front-ends, and more specifically linear power amplifiers (PAs), to avoid distortion. This reduces the energy efficiency since PAs are most efficient close to saturation, where they generate most nonlinear distortion. Moreover, the use of conventional precoders can induce a coherent combining of distortion at the user locations, degrading the signal quality. In this work, novel linear precoders, simple to compute and to implement, are proposed that allow working close to saturation, while cancelling the third-order nonlinearity of the PA without prior knowledge of the signal statistics and PA model. Their design consists in saturating a single or a few antennas on purpose together with an negative gain with respect to all other antennas to compensate for the overall nonlinear distortion at the user location. The performance gains of the designs are significant for PAs working close to saturation, as compared to maximum ratio transmission (MRT) precoding and perfect per-antenna digital pre-distortion (DPD) compensation.


title: A Comprehensive Review on Power System Risk-Based Transient Stability

abstract: Power systems are getting more complex than ever and are consequently operating close to their limit of stability. Moreover, with the increasing demand of renewable wind generation, and the requirement to maintain a secure power system, the importance of transient stability cannot be overestimated. Considering its significance in power system security, it is important to suggest a different methodology for enhancing the transient stability, considering uncertainties. Current deterministic industry practices of transient stability assessment ignore the probabilistic nature of variables (fault type, fault location, fault clearing time, etc.). These approaches typically provide a cautious principle and can result in high-priced expansion projects or operational limits. With the increasing system uncertainties and widespread electricity market deregulation, there is a strong inevitability to incorporate risk in the traditional transient stability analysis. Accurate assessment of transient stability in a modern power network is becoming a strict requirement both in planning and in real-time operation, due to the increasingly intricate dynamics of a power system. Further, increasing sources of uncertainty in forecast state and in the reaction to faults highly implies the implementation of risk-based approach in assessing transient stability. Thus, this paper aims to provide a comprehensive review of risk-based transient stability in power networks and the accompanying research. It is believed that this review can be an inception for researchers in the domain of power system planning and security.


title: Energy Conditions in Extended $f(R,G,T)$ Gravity

abstract: In this paper, we consider the flat Friedmann Lematre Robertson-Walke metric in the presence of perfect fluid models and extended $f(R,G,T)$ gravity (where $R$ is the Ricci scalar, $G$ is the Gauss Bonnet invariant and $T$ stands for trace of energy momentum tensor). In this context, we assume some specific realistic $f(R,G,T)$ models configuration that could be used to explore the finite-time future singularities that arise in late-time cosmic accelerating phases. In this scenario, we choose the most recent estimated values for the Hubble, deceleration, snap and jerk parameters to develop the viability and bounds on the models parameters induced by different energy conditions.


title: Surface engineering for ultrathin metal anodes enabling high-performance Zn-ion batteries

abstract: Zn metal battery has been considered a promising alternative energy storage technology in renewable energy storage and grid storage. It is well-known that the surface orientation of a Zn metal anode is vital to the reversibility of a Zn metal battery. Herein, the (101)-oriented thin Zn metal anode (down to 2 {\mu}m) is electrodeposited on a Cu surface by adding dimethyl sulfoxide (DMSO) electrolyte additive in ZnSO4 aqueous solution. Scanning electron microscope (SEM) observation indicates the formation of flat terrace-like compact (101)-oriented surfaces. Insitu optical observation confirms that the (101)-oriented surfaces can be reversibly plated and stripped. DFT calculations reveal two mechanisms for the nucleation and growth of the Zn-(101) surface: (1) formation of Zn(101)//Cu(001) could lower the interface energy as compared to Zn(002)//Cu(001); (2) large reconstruction of the Zn (101) surface with DMSO and H2O absorption. Raman, XPS, and ToF-SIMS characterizations indicate that adding DMSO in ZnCl2 could facilitate the formation of ZnO-based SEI on Zn metal surface, while OH- and S-based SEI can be obtained with DMSO in ZnSO4. The electrochemical testings are performed, which demonstrates a higher cyclability for the (101)-oriented Zn in the half cell as well as a lower charge transfer barrier with respect to the (002)-dominated surface of the same electrode thickness. Zn||V2O5 full cells are further assembled, showing better capacity retention for the (101)-Zn as compared to the (002)-Zn with the same thickness (5 {\mu}m, 3 {\mu}m, and 2 {\mu}m). We hope this study to spur further interest in the control of Zn metal surface crystallographic orientation towards ultrathin Zn metal anodes.


title: Characterization of a kg-scale archaeological lead-based cryogenic detectors for the RES-NOVA experiment

abstract: One of the most energetic events in the Universe are core-collapse Supernovae (SNe), where almost all the star's binding energy is released as neutrinos. These particles are direct probes of the processes occurring in the stellar core and provide unique insights into the gravitational collapse. RES-NOVA will revolutionize how we detect neutrinos from astrophysical sources, by deploying the first ton-scale array of cryogenic detectors made from archaeological lead. Pb offers the highest neutrino interaction cross-section via coherent elastic neutrino-nucleus scattering (CE$\nu$NS). Such process will enable RES-NOVA to be equally sensitive to all neutrino flavors. For the first time, we propose to use archaeological Pb as sensitive target material in order to achieve an ultra-low background level in the region of interest (\textit{O}(1keV)). All these features make possible the deployment of the first cm-scale neutrino telescope for the investigation of astrophysical sources. In this contribution, we will characterize the radiopurity level and the performance of a small-scale proof-of-principle detector of RES-NOVA, consisting in a PbWO$_4$ crystal made from archaeological-Pb operated as cryogenic detector.


title: Optimal time decay estimation for large-solution about 3D compressible MHD equations

abstract: This paper mainly focus on optimal time decay estimation for large-solution about compressible magnetohydrodynamic equations in 3D whole space, provided that $(\sigma_{0}-1,u_{0},M_{0})\in L^1\cap H^2$. In [2](Chen et al.,2019), they proved time decay estimation of $\|(\sigma-1,u,M)\|_{H^1}$ being $(1+t)^{-\frac{3}{4}}$. Based on it, we obtained that of $\|\nabla(\sigma-1,u,M)\|_{H^1}$ being $(1+t)^{-\frac{5}{4}}$ in [24]. Therefore, we are committed to improving that of $\|\nabla^2 (\sigma-1,u,M)\|_{L^2}$ in this paper. Thanks to the method adopted in [25] (Wang and Wen, 2021), we get the optimal time decay estimation to the highest-order derivative for space of solution, which means that time decay estimation of $\|\nabla^2 (\sigma-1,u,M)\|_{L^2}$ is $(1+t)^{-\frac{7}{4}}$.


title: Asymptotics of the number of primes and sums of functions of primes in a subset of natural series

abstract: The paper solves the problems of determining the asymptotics of the number of primes and the sums of functions of primes in a subset of the natural series that satisfies the conditions that the asymptotic density of the number of primes in this subset is constant and not equal to zero.


title: X-ray technological irradiation for TID studies on silicon sensor and electronic devices in a medical facility

abstract: Technological tests of Total Ionizing Dose effects are required not only for silicon particle sensors developed in high energy physics experiments, but also for electronic devices and semiconductor elements used in commercial, automotive and space applications. Using x-ray irradiators, these technological tests and studies can be performed not only in facilities explicitly built for this mission, but also in medical or biological research facilities when some minima requirements are satisfied. Generally this irradiation can be performed without interfering with the medical and biological activities of the facility. In this article will be described the minimum instrumentation required for this type of studies and will be given a detailed description of the preparation and realization of a SiPM x-ray irradiation campaign for TID characterization realized in the Italian TIFPA-INFN Trento Center Laboratory in May 2021, using instruments originally realized for medical and biological irradiation.


title: Active information, missing data and prevalence estimation

abstract: The topic of this paper is prevalence estimation from the perspective of active information. Prevalence among tested individuals has an upward bias under the assumption that individuals' willingness to be tested for the disease increases with the strength of their symptoms. Active information due to testing bias quantifies the degree at which the willingness to be tested correlates with infection status. Interpreting incomplete testing as a missing data problem, the missingness mechanism impacts the degree at which the bias of the original prevalence estimate can be removed. The reduction in prevalence, when testing bias is adjusted for, translates into an active information due to bias correction, with opposite sign to active information due to testing bias. Prevalence and active information estimates are asymptotically normal, a behavior also illustrated through simulations.

categories: stat.ME cs.IT math.IT physics.data-an q-bio.PE q-bio.QM


title: Topological spiral magnetism in the Weyl semimetal SmAlSi

abstract: Weyl electrons are intensely studied due to novel charge transport phenomena such as chiral anomaly, Fermi arcs, and photogalvanic effect. Recent theoretical works suggest that Weyl electrons can also participate in magnetic interactions, and the Weyl-mediated indirect exchange coupling between local moments is proposed as a new mechanism of spiral magnetism that involves chiral electrons. Despite reports of incommensurate and non-collinear magnetic ordering in Weyl semimetals, an actual spiral order has remained hitherto undetected. Here, we present evidence of Weyl-mediated spiral magnetism in SmAlSi from neutron diffraction, transport, and thermodynamic data. We show that the spiral order in SmAlSi results from the nesting between topologically non-trivial Fermi pockets and weak magnetocrystalline anisotropy, unlike related materials (Ce,Pr,Nd)AlSi, where a strong anisotropy prevents the spins from freely rotating. We map the magnetic phase diagram of SmAlSi and reveal an A-phase where topological magnetic excitations may exist. This is corroborated by the observation of a topological Hall effect within the A-phase.


title: On the Rayleigh-Ritz variational method

abstract: We give a simple proof of the well known fact that the approximate eigenvalues provided by the Rayleigh-Ritz variational method are increasingly accurate upper bounds to the exact ones. To this end, we resort to the variational principle, mentioned in most textbooks on quantum chemistry, and to a well known set of projection operators. We think that present approach may be suitable for an advanced course on quantum mechanics or quantum chemistry.


title: REKnow: Enhanced Knowledge for Joint Entity and Relation Extraction

abstract: Relation extraction is an important but challenging task that aims to extract all hidden relational facts from the text. With the development of deep language models, relation extraction methods have achieved good performance on various benchmarks. However, we observe two shortcomings of previous methods: first, there is no unified framework that works well under various relation extraction settings; second, effectively utilizing external knowledge as background information is absent. In this work, we propose a knowledge-enhanced generative model to mitigate these two issues. Our generative model is a unified framework to sequentially generate relational triplets under various relation extraction settings and explicitly utilizes relevant knowledge from Knowledge Graph (KG) to resolve ambiguities. Our model achieves superior performance on multiple benchmarks and settings, including WebNLG, NYT10, and TACRED.


title: Stochastic Zeroth order Descent with Structured Directions

abstract: We introduce and analyze Structured Stochastic Zeroth order Descent (S-SZD), a finite difference approach which approximates a stochastic gradient on a set of $l\leq d$ orthogonal directions, where $d$ is the dimension of the ambient space. These directions are randomly chosen, and may change at each step. For smooth convex functions we prove almost sure convergence of the iterates and a convergence rate on the function values of the form $O(d/l k^{-c})$ for every $c<1/2$, which is arbitrarily close to the one of Stochastic Gradient Descent (SGD) in terms of number of iterations. Our bound also shows the benefits of using $l$ multiple directions instead of one. For non-convex functions satisfying the Polyak-{\L}ojasiewicz condition, we establish the first convergence rates for stochastic zeroth order algorithms under such an assumption. We corroborate our theoretical findings in numerical simulations where assumptions are satisfied and on the real-world problem of hyper-parameter optimization, observing that S-SZD has very good practical performances.


title: Boosting unstable particles

abstract: In relativity, there is no absolute notion of simultaneity, because two clocks that are in different places can always be desynchronized by a Lorentz boost. Here, we explore the implications of this effect for the quantum theory of unstable particles. We show that, when a wavefunction is boosted, its tails travel one to the past and the other to the future. As a consequence, in the new frame of reference, the particle is in a quantum superposition "decayed + non decayed", where the property "decayed-ness" is entangled with the position. Since a particle cannot be localised in a region smaller than the Compton wavelength, there is a non-zero lower bound on this effect, which is fundamental in nature. The surprising implication is that, in a quantum world, decay probabilities can never be Lorentz-invariant. We show that this insight was the missing ingredient to reconcile the seemingly conflicting views about time dilation in relativistic quantum mechanics and quantum field theory.


title: Quasi-Likelihood Analysis of Fractional Brownian Motion with Constant Drift under High-Frequency Observations

abstract: Consider an estimation of the Hurst parameter $H\in(0,1)$ and the volatility parameter $\sigma>0$ for a fractional Brownian motion with a drift term under high-frequency observations with a finite time interval. In the present paper, we propose a consistent estimator of the parameter $\theta=(H,\sigma)$ combining the ideas of a quasi-likelihood function based on a local Gaussian approximation of a high-frequently observed time series and its frequency-domain approximation. Moreover, we prove an asymptotic normality property of the proposed estimator for all $H\in(0,1)$ when the drift process is constant.


title: Globally-Optimal Contrast Maximisation for Event Cameras

abstract: Event cameras are bio-inspired sensors that perform well in challenging illumination conditions and have high temporal resolution. However, their concept is fundamentally different from traditional frame-based cameras. The pixels of an event camera operate independently and asynchronously. They measure changes of the logarithmic brightness and return them in the highly discretised form of time-stamped events indicating a relative change of a certain quantity since the last event. New models and algorithms are needed to process this kind of measurements. The present work looks at several motion estimation problems with event cameras. The flow of the events is modelled by a general homographic warping in a space-time volume, and the objective is formulated as a maximisation of contrast within the image of warped events. Our core contribution consists of deriving globally optimal solutions to these generally non-convex problems, which removes the dependency on a good initial guess plaguing existing methods. Our methods rely on branch-and-bound optimisation and employ novel and efficient, recursive upper and lower bounds derived for six different contrast estimation functions. The practical validity of our approach is demonstrated by a successful application to three different event camera motion estimation problems.


title: Real-time Hyper-Dimensional Reconfiguration at the Edge using Hardware Accelerators

abstract: In this paper we present Hyper-Dimensional Reconfigurable Analytics at the Tactical Edge (HyDRATE) using low-SWaP embedded hardware that can perform real-time reconfiguration at the edge leveraging non-MAC (free of floating-point MultiplyACcumulate operations) deep neural nets (DNN) combined with hyperdimensional (HD) computing accelerators. We describe the algorithm, trained quantized model generation, and simulated performance of a feature extractor free of multiply-accumulates feeding a hyperdimensional logic-based classifier. Then we show how performance increases with the number of hyperdimensions. We describe the realized low-SWaP FPGA hardware and embedded software system compared to traditional DNNs and detail the implemented hardware accelerators. We discuss the measured system latency and power, noise robustness due to use of learnable quantization and HD computing, actual versus simulated system performance for a video activity classification task and demonstration of reconfiguration on this same dataset. We show that reconfigurability in the field is achieved by retraining only the feed-forward HD classifier without gradient descent backpropagation (gradient-free), using few-shot learning of new classes at the edge. Initial work performed used LRCN DNN and is currently extended to use Two-stream DNN with improved performance.


title: Inverting Incomplete Fourier Transforms by a Sparse Regularization Model and Applications in Seismic Wavefield Modeling

abstract: We propose a sparse regularization model for inversion of incomplete Fourier transforms and apply it to seismic wavefield modeling. The objective function of the proposed model employs the Moreau envelope of the $\ell_0$ norm under a tight framelet system as a regularization to promote sparsity. This model leads to a non-smooth, non-convex optimization problem for which traditional iteration schemes are inefficient or even divergent. By exploiting special structures of the $\ell_0$ norm, we identify a local minimizer of the proposed non-convex optimization problem with a global minimizer of a convex optimization problem, which provides us insights for the development of efficient and convergence guaranteed algorithms to solve it. We characterize the solution of the regularization model in terms of a fixed-point of a map defined by the proximity operator of the $\ell_0$ norm and develop a fixed-point iteration algorithm to solve it. By connecting the map with an $\alpha$-averaged nonexpansive operator, we prove that the sequence generated by the proposed fixed-point proximity algorithm converges to a local minimizer of the proposed model. Our numerical examples confirm that the proposed model outperforms significantly the existing model based on the $\ell_1$-norm. The seismic wavefield modeling in the frequency domain requires solving a series of the Helmholtz equation with large wave numbers, which is a computationally intensive task. Applying the proposed sparse regularization model to the seismic wavefield modeling requires data of only a few low frequencies, avoiding solving the Helmholtz equation with large wave numbers. Numerical results show that the proposed method performs better than the existing method based on the $\ell_1$ norm in terms of the SNR values and visual quality of the restored synthetic seismograms.


title: Final Result on the Neutrinoless Double Beta Decay of $^{82}$Se with CUPID-0

abstract: CUPID-0, an array of Zn$^{82}$Se cryogenic calorimeters, was the first medium-scale demonstrator of the scintillating bolometers technology. The first project phase (March 2017 - December 2018) allowed the most stringent limit on the neutrinoless double beta decay half-life of the isotope of interest, $^{82}$Se, to be set. After a six months long detector upgrade, CUPID-0 began its second and last phase (June 2019 - February 2020). In this letter, we describe the search for neutrinoless double beta decay of $^{82}$Se with a total exposure (phase I + II) of 8.82 kg$\times$yr of isotope. We set a limit on the half-life of $^{82}$Se to the ground state of $^{82}$Kr of T$^{0\nu}_{1/2}$($^{82}$Se)$>$ 4.6$\times \mathrm{10}^{24}$ yr (90\% credible interval), corresponding to an effective Majorana neutrino mass m$_{\beta\beta} <$ (263 -- 545) meV. We also set the most stringent lower limits on the neutrinoless decays of $^{82}$Se to the 0$_1^+$, 2$_1^+$ and 2$_2^+$ excited states of $^{82}$Kr, finding 1.8$\times$10$^{23}$ yr, 3.0$\times$10$^{23}$ yr, 3.2$\times$10$^{23}$ yr (90$\%$ credible interval) respectively.


title: Single-parameter aging in the weakly nonlinear limit

abstract: Physical aging deals with slow property changes over time caused by molecular rearrangements. This is relevant for non-crystalline materials like polymers and inorganic glasses, both in production and during subsequent use. The Narayanaswamy theory from 1971 describes physical aging - an inherently nonlinear phenomenon - in terms of a linear convolution integral over the so-called material time $\xi$. The resulting "Tool-Narayanaswamy (TN) formalism" is generally recognized to provide an excellent description of physical aging for small, but still highly nonlinear temperature variations. The simplest version of the TN formalism is single-parameter aging according to which the clock rate $d\xi/dt$ is an exponential function of the property monitored [T. Hecksher et al., J. Chem. Phys. 142, 241103 (2015)]. For temperature jumps starting from thermal equilibrium, this leads to a first-order differential equation for property monitored, involving a system-specific function. The present paper shows analytically that the solution to this equation to first order in the temperature variation has a universal expression in terms of the zeroth-order solution, $R_0(t)$. Numerical data for a binary Lennard-Jones glass former probing the potential energy confirm that, in the weakly nonlinear limit, the theory predicts aging correctly from $R_0(t)$ (which by the fluctuation-dissipation theorem is the normalized equilibrium potential-energy time-autocorrelation function).


title: Black holes in dRGT massive gravity with the signature of EHT observations of M87*

abstract: The recent Event Horizon Telescope (EHT) observations of the M87* have led to a surge of interest in studying the shadow of black holes. Besides, investigation of time evolution and lifetime of black holes helps us to veto/restrict some theoretical models in gravitating systems. Motivated by such exciting properties, we study optical features of black holes, such as the shadow geometrical shape and the energy emission rate in modified gravity. We consider a charged AdS black hole in dRGT massive gravity and look for criteria to restrict the free parameters of the theory. The main goal of this paper is to compare the shadow of the mentioned black hole in a rotating case with the EHT data to obtain the allowed regions of the model parameters. Therefore, we employ the Newman-Janis algorithm to build the rotating counterpart of static solution in dRGT massive gravity. We also calculate the energy emission rate for the rotating case and discuss how the rotation factor and other parameters affect the emission of particles around the black holes.


title: On the square-root approximation finite volume scheme for nonlinear drift-diffusion equations

abstract: We study a finite volume scheme for the approximation of the solution to convection diffusion equations with nonlinear convection and Robin boundary conditions. The scheme builds on the interpretation of such a continuous equation as the hydrodynamic limit of some simple exclusion jump process. We show that the scheme admits a unique discrete solution, that the natural bounds on the solution are preserved, and that it encodes the second principle of thermodynamics in the sense that some free energy is dissipated along time. The convergence of the scheme is then rigorously established thanks to compactness arguments. Numerical simulations are finally provided, highlighting the overall good behavior of the scheme.


title: Distributionally Robust End-to-End Portfolio Construction

abstract: We propose an end-to-end distributionally robust system for portfolio construction that integrates the asset return prediction model with a distributionally robust portfolio optimization model. We also show how to learn the risk-tolerance parameter and the degree of robustness directly from data. End-to-end systems have an advantage in that information can be communicated between the prediction and decision layers during training, allowing the parameters to be trained for the final task rather than solely for predictive performance. However, existing end-to-end systems are not able to quantify and correct for the impact of model risk on the decision layer. Our proposed distributionally robust end-to-end portfolio selection system explicitly accounts for the impact of model risk. The decision layer chooses portfolios by solving a minimax problem where the distribution of the asset returns is assumed to belong to an ambiguity set centered around a nominal distribution. Using convex duality, we recast the minimax problem in a form that allows for efficient training of the end-to-end system.


title: On some properties of random and pseudorandom codes

abstract: We describe some pseudorandom properties of binary linear codes achieving capacity on the binary erasure channel under bit-MAP decoding (as shown in Kudekar et al this includes doubly transitive codes and, in particular, Reed-Muller codes). We show that for all integer $q \ge 2$ the $\ell_q$ norm of the characteristic function of such 'pseudorandom' code decreases as fast as that of any code of the same rate (and equally fast as that of a random code) under the action of the noise operator. In information-theoretic terms this means that the $q^{th}$ R\'enyi entropy of this code increases as fast as possible over the binary symmetric channel. In particular (taking $q = \infty$) this shows that such codes have the smallest asymptotic undetected error probability (equal to that of a random code) over the BSC, for a certain range of parameters. We also study the number of times a certain local pattern, a 'rhombic' $4$-tuple of codewords, appears in a linear code, and show that for a certain range of parameters this number for pseudorandom codes is similar to that for a random code.


title: Fast Deep Autoencoder for Federated learning

abstract: This paper presents a novel, fast and privacy preserving implementation of deep autoencoders. DAEF (Deep Autoencoder for Federated learning), unlike traditional neural networks, trains a deep autoencoder network in a non-iterative way, which drastically reduces its training time. Its training can be carried out in a distributed way (several partitions of the dataset in parallel) and incrementally (aggregation of partial models), and due to its mathematical formulation, the data that is exchanged does not endanger the privacy of the users. This makes DAEF a valid method for edge computing and federated learning scenarios. The method has been evaluated and compared to traditional (iterative) deep autoencoders using seven real anomaly detection datasets, and their performance have been shown to be similar despite DAEF's faster training.


title: Harder-Narasimhan polygons and Laws of Large Numbers

abstract: We build on the recent techniques of Codogni and Patakfalvi, from \cite{Codogni:Patakfalvi:2021}, which were used to establish theorems about semi-positivity of the Chow Mumford line bundles for families of $\K$-semistable Fano varieties. Here we apply the Central Limit Theorem to ascertain the asymptotic probabilistic nature of the vertices of the \emph{Harder and Narasimhan polygons}. As an application of our main result, we use it to establish a filtered vector space analogue of the main technical result of \cite{Codogni:Patakfalvi:2021}. In doing so, we expand upon the slope stability theory, for filtered vector spaces, that was initiated by Faltings and W\"{u}stholz \cite{Faltings:Wustholz}. One source of inspiration for our abstract study of \emph{Harder and Narasimhan data}, which is a concept that we define here, is the lattice reduction methods of Grayson \cite{Grayson:1984}. Another is the work of Faltings and W\"{u}stholz, \cite{Faltings:Wustholz}, and Evertse and Ferretti, \cite{Evertse:Ferretti:2013}, which is within the context of Diophantine approximation for projective varieties.


title: Functional limit theorems for P\'olya urns with growing initial compositions

abstract: In this paper, we prove functional limit theorems for P\'olya urn processes whose number of draws and initial number of balls tend to infinity together. This is motivated by recent work of Borovkov [5], where they prove a functional limit theorem for this model when the urn has identity replacement rule. We generalize this result to arbitrary balanced replacement rules (the total number of balls added to the urn is deterministic). Three asymptotic regimes are possible depending on how one lets the number of initial balls scale with the number of draws of the urn. In each regime, we show a first order deterministic limit and Gaussian second order fluctuations, where the behaviour of these limit processes depend on the regime, the initial composition of the urn, and the urns replacement rule. To prove our main results, we embed the process in continuous-time and use martingale theory. Although these methods are classical since the works of Athreya & Karlin [1] and Janson [11], our setting with initial growing composition necessitates many new ideas. The main difference in proving limiting results for the continuous time embedding in our setting is that, when the initial composition is large compared to the number of draws, the branching process does not have time to reach equilibrium. Because of this, translating the results back to discrete-time is also much harder than in Janson [11]. Interestingly, our continuous-time results hold under weaker assumptions on the replacement structure than classical results for multi-type branching processes; in particular, we do not need any "irreducibility" assumption.


title: Finite electro-elasticity with physics-augmented neural networks

abstract: In the present work, a machine learning based constitutive model for electro-mechanically coupled material behavior at finite deformations is proposed. Using different sets of invariants as inputs, an internal energy density is formulated as a convex neural network. In this way, the model fulfills the polyconvexity condition which ensures material stability, as well as thermodynamic consistency, objectivity, material symmetry, and growth conditions. Depending on the considered invariants, this physics-augmented machine learning model can either be applied for compressible or nearly incompressible material behavior, as well as for arbitrary material symmetry classes. The applicability and versatility of the approach is demonstrated by calibrating it on transversely isotropic data generated with an analytical potential, as well as for the effective constitutive modeling of an analytically homogenized, transversely isotropic rank-one laminate composite and a numerically homogenized cubic metamaterial. These examinations show the excellent generalization properties that physics-augmented neural networks offer also for multi-physical material modeling such as nonlinear electro-elasticity.


title: Evolution of a Relativistic Outflow and X-ray Corona in the Extreme Changing-Look AGN 1ES 1927+654

abstract: 1ES 1927+654 is a paradigm-defying AGN and one of the most peculiar X-ray nuclear transients. In early 2018, this well-known AGN underwent a changing-look event, in which broad optical emission lines appeared and the optical flux increased. Yet, by July 2018, the X-ray flux had dropped by over two orders of magnitude, indicating a dramatic change to the inner accretion flow. With three years of observations with NICER, XMM-Newton, and NuSTAR, we present the X-ray evolution of 1ES 1927+654, which can be broken into three phases-(1) an early super-Eddington phase with rapid variability in X-ray luminosity and spectral parameters, (2) a stable super-Eddington phase at the peak X-ray luminosity, and (3) a steady decline back to the pre-outburst luminosity and spectral parameters. For the first time, we witnessed the formation of the X-ray corona, as the X-ray spectrum transitioned from thermally-dominated to primarily Comptonized. We also track the evolution of the prominent, broad 1 keV feature in the early X-ray spectra and show that this feature can be modeled with blueshifted reflection (z = -0.33) from a single-temperature blackbody irradiating spectrum using xillverTDE, a new flavor of the xillver models. Thus, we propose that the 1 keV feature could arise from reflected emission off the base of an optically thick outflow from a geometrically thick, super-Eddington inner accretion flow, connecting the inner accretion flow with outflows launched during extreme accretion events (e.g. tidal disruption events). Lastly, we compare 1ES 1927+654 to other nuclear transients and discuss applications of xillverTDE to super-Eddington accretors.


title: A simulation study of the estimation quality in the double-Cox model with shared frailty for non-proportional hazards survival analysis

abstract: The Cox regression, a semi-parametric method of survival analysis, is extremely popular in biomedical applications. The proportional hazards assumption is a key requirement in the Cox model. To accommodate non-proportional hazards, we propose to parameterise the shape parameter of the baseline hazard function using the additional, separate Cox-regression term which depends on the vector of the covariates. We call this model the double-Cox model. The R programs for fitting the double-Cox model are available on Github. We formally introduce the double-Cox model with shared frailty and investigate, by simulation, the estimation bias and the coverage of the proposed point and interval estimation methods for the Gompertz and the Weibull baseline hazards. In applications with low frailty variance and a large number of clusters, the marginal likelihood estimation is almost unbiased and the profile likelihood-based confidence intervals provide good coverage for all model parameters. We also compare the results from the over-fitted double-Cox model to those from the standard Cox model with frailty in the case of the scale-only proportional hazards. Results of our simulations on the bias and coverage of the model parameters are provided in 12 Tables and in 145 A4 Figures, 178 pages in total.


title: Pluricanonical Maps and the Fujita Conjecture

abstract: We describe examples showing the sharpness of Fujita's conjecture on adjoint bundles also in the general type case, and use these examples to formulate related bold conjectures on pluricanonical maps of varieties of general type.


title: On stable solutions to the Euler equations in convex planar domains

abstract: In convex planar domains, given an initial vorticity with one sign, we study the regularity and geometric properties of the dynamically stable solutions to the Euler equations in the coadjoint orbit of the initial vorticity. These flows have elliptic stagnation points. Under some nondegeneracy conditions on the data, we show they are Holder continuous and have convex level curves. We also give a detailed description for the set of stagnation points. If the initial vorticity has nice level set topology, these stable solutions are in the L^\infty-strong closure of the coadjoint orbit. We also demonstrate the sharpness of most assumptions we made.


title: Pulse-level Scheduling of Quantum Circuits for Neutral-Atom Devices

abstract: We show how a pulse-level implementation of the multi-qubit gates in neutral-atom device architectures allows for the simultaneous execution of single- and multi-qubit gates acting on overlapping sets of qubits, in a mechanism we name absorption. With absorption as a foundation, we present an algorithm to schedule the execution of a quantum circuit as a pulse sequence on a neutral-atom device with a single channel for single- and multi-qubit gate execution. For any quantum circuit of practical relevance, we observe that the algorithm results in an optimal utilization of the available resources that cannot be surpassed by a different scheduling strategy. Our benchmarks against a custom scheduler attempting to maximize parallelization at the gate level show the time gained by the pulse-level scheduler is proportional to the depth and is most pronounced for quantum circuits with fewer qubits.


title: Practical limits for large-momentum-transfer clock atom interferometers

abstract: Atom interferometry on optical clock transitions is being pursued for numerous long-baseline experiments both terrestrially and for future space missions. Crucial to meeting these experiments' required sensitivities is the implementation of large momentum transfer ($>$10$^3\hbar k$). Here we show that to sequentially apply such a large momentum via $\pi$ pulses places stringent requirements on the frequency noise of the interferometry laser, finding that the linewidth is required to be considerably lower than has previously been suggested. This is due to imperfect pulse fidelity in the presence of noise and is apparent even for an atom at rest interacting with resonant light, making this a fundamental constraint on operational fidelity for a given laser and pulse sequence. Within this framework, we further present and analyse two high-power, frequency-stabilised laser sources designed to perform interferometry on the $^1$S$_0$ - $^3$P$_0$ clock transitions of cadmium and strontium, respectively operating at 332 nm and 698 nm.


title: Global Internet public peering capacity of interconnection: a complex network analysis

abstract: A massive and growing part of Autonomous System (AS)-level traffic exchanges takes place at Internet Exchange Points (IXPs). This paper leverages PeeringDB, a database providing a partial but reasonable view of the global interconnection of ASes at IXPs, to model a complex graph enabling the characterization of the key Internet peering players and their interactions over time. We model a PeeringDB snapshot as a weighted directed bipartite graph, called the pDB c-graph, that captures the port size ASes possess at IXPs using available metadata. This novel model of the Internet is shown to picture relevant features of a complex network that groups ASes and IXPs in geographical areas of influence. From this model, we extract central players of public peering such as hypergiant AS content providers and major regional traffic receivers. Most importantly, this graph model opens the way to apply spectral analysis using reduced Google matrix in order to retrieve the intensity of possible interactions between ASes on the basis of pure connectivity information. As an illustration, we retrieve the timely evolution of the peering network to show how the central content and cloud providers have increased their reach to eyeball networks during Covid-19 pandemic.


title: Accurate Modeling of Lyman-alpha Profiles and their Impact on Photolysis of Terrestrial Planet Atmospheres

abstract: Accurately measuring and modeling the Lyman-$\alpha$ (Ly$\alpha$; $\lambda$1215.67 \AA) emission line from low mass stars is vital for our ability to build predictive high energy stellar spectra, yet interstellar medium (ISM) absorption of this line typically prevents model-measurement comparisons. Ly$\alpha$ also controls the photodissociation of important molecules, like water and methane, in exoplanet atmospheres such that any photochemical models assessing potential biosignatures or atmospheric abundances require accurate Ly$\alpha$ host star flux estimates. Recent observations of three early M and K stars (K3, M0, M1) with exceptionally high radial velocities (>100 km s$^{-1}$) reveal the intrinsic profiles of these types of stars as most of their Ly$\alpha$ flux is shifted away from the geocoronal line core and contamination from the ISM. These observations indicate that previous stellar spectra computed with the PHOENIX atmosphere code have underpredicted the core of Ly$\alpha$ in these types of stars. With these observations, we have been able to better understand the microphysics in the upper atmosphere and improve the predictive capabilities of the PHOENIX atmosphere code. Since these wavelengths drive the photolysis of key molecular species, we also present results analyzing the impact of the resulting changes to the synthetic stellar spectra on observable chemistry in terrestrial planet atmospheres.


title: Weakly-supervised segmentation using inherently-explainable classification models and their application to brain tumour classification

abstract: Deep learning models have shown their potential for several applications. However, most of the models are opaque and difficult to trust due to their complex reasoning - commonly known as the black-box problem. Some fields, such as medicine, require a high degree of transparency to accept and adopt such technologies. Consequently, creating explainable/interpretable models or applying post-hoc methods on classifiers to build trust in deep learning models are required. Moreover, deep learning methods can be used for segmentation tasks, which typically require hard-to-obtain, time-consuming manually-annotated segmentation labels for training. This paper introduces three inherently-explainable classifiers to tackle both of these problems as one. The localisation heatmaps provided by the networks -- representing the models' focus areas and being used in classification decision-making -- can be directly interpreted, without requiring any post-hoc methods to derive information for model explanation. The models are trained by using the input image and only the classification labels as ground-truth in a supervised fashion - without using any information about the location of the region of interest (i.e. the segmentation labels), making the segmentation training of the models weakly-supervised through classification labels. The final segmentation is obtained by thresholding these heatmaps. The models were employed for the task of multi-class brain tumour classification using two different datasets, resulting in the best F1-score of 0.93 for the supervised classification task while securing a median Dice score of 0.67$\pm$0.08 for the weakly-supervised segmentation task. Furthermore, the obtained accuracy on a subset of tumour-only images outperformed the state-of-the-art glioma tumour grading binary classifiers with the best model achieving 98.7\% accuracy.


title: Referring Image Matting

abstract: Different from conventional image matting, which either requires user-defined scribbles/trimap to extract a specific foreground object or directly extracts all the foreground objects in the image indiscriminately, we introduce a new task named Referring Image Matting (RIM) in this paper, which aims to extract the meticulous alpha matte of the specific object that best matches the given natural language description, thus enabling a more natural and simpler instruction for image matting. First, we establish a large-scale challenging dataset RefMatte by designing a comprehensive image composition and expression generation engine to automatically produce high-quality images along with diverse text attributes based on public datasets. RefMatte consists of 230 object categories, 47,500 images, 118,749 expression-region entities, and 474,996 expressions. Additionally, we construct a real-world test set with 100 high-resolution natural images and manually annotate complex phrases to evaluate the out-of-domain generalization abilities of RIM methods. Furthermore, we present a novel baseline method CLIPMat for RIM, including a context-embedded prompt, a text-driven semantic pop-up, and a multi-level details extractor. Extensive experiments on RefMatte in both keyword and expression settings validate the superiority of CLIPMat over representative methods. We hope this work could provide novel insights into image matting and encourage more follow-up studies. The dataset, code and models are available at https://github.com/JizhiziLi/RIM.


title: Curvature Properties of 3-$(\alpha,\delta)$-Sasaki Manifolds

abstract: We investigate curvature properties of 3-$(\alpha,\delta)$-Sasaki manifolds, a special class of almost 3-contact metric manifolds generalizing 3-Sasaki manifolds (corresponding to $\alpha = \delta = 1$) that admit a canonical metric connection with skew torsion and define a Riemannian submersion over a quaternionic K\"ahler manifold with vanishing, positive or negative scalar curvature, according to $\delta = 0$, $\alpha\delta > 0$ or $\alpha\delta < 0$. We shall investigate both the Riemannian curvature and the curvature of the canonical connection, with particular focus on their curvature operators, regarded as symmetric endomorphisms of the space of 2-forms. We describe their spectrum, find distinguished eigenforms, and study the conditions of strongly definite curvature in the sense of Thorpe.


title: EOS-ESTM: A flexible climate model for habitable exoplanets

abstract: Rocky planets with temperate conditions provide the best chance for discovering habitable worlds and life outside the Solar System. In the last decades, new instrumental facilities and large observational campaigns have been driven by the quest for habitable worlds. Climate models aimed at studying the habitability of rocky planets are essential tools to pay off these technological and observational endeavours. In this context, we present EOS-ESTM, a fast and flexible model aimed at exploring the impact on habitability of multiple climate factors, including those unconstrained by observations. EOS-ESTM is built on ESTM, a seasonal-latitudinal energy balance model featuring an advanced treatment of the meridional and vertical transport. The novel features of EOS-ESTM include: (1) parameterizations for simulating the climate impact of oceans, land, ice, and clouds as a function of temperature and stellar zenith distance; (2) a procedure (EOS) for calculating the radiative transfer in atmospheres with terrestrial and non-terrestrial compositions illuminated by solar- and non-solar-type stars. By feeding EOS-ESTM with Earth's stellar, orbital and planetary parameters we derive a reference model that satisfies a large number of observational constraints of the Earth's climate system. Validation tests of non-terrestrial conditions yield predictions that are in line with comparable results obtained with a hierarchy of climate models. The application of EOS-ESTM to planetary atmospheres in maximum greenhouse conditions demonstrates the possibility of tracking the snowball transition at the outer edge of the HZ for a variety of planetary parameters, paving the road for multi-parametric studies of the HZ.


title: Toward a new theory of the fractional quantum Hall effect

abstract: The fractional quantum Hall effect was experimentally discovered in 1982. It was observed that the Hall conductivity $\sigma_{yx}$ of a two-dimensional electron system is quantized, $\sigma_{yx}=e^2/3h$, in the vicinity of the Landau level filling factor $\nu=1/3$. In 1983, Laughlin proposed a trial many-body wave function, which he claimed described a ``new state of matter'' -- a homogeneous incompressible liquid with fractionally charged quasiparticles. Here I develop an exact diagonalization theory that allows calculation of the energy and other physical properties of the ground and excited states of a system of $N$ two-dimensional Coulomb interacting electrons in a strong magnetic field. I analyze the energies, electron densities, and other physical properties of the systems with $N\le 7$ electrons, continuously as a function of magnetic field in the range $1/4\lesssim\nu<1$. The results show that both the ground and excited states of the system resemble a sliding Wigner crystal, whose parameters are influenced by the magnetic field. Energy gaps in the many-particle spectra appear and disappear as the magnetic field changes. I also calculate the physical properties of the $\nu=1/3$ Laughlin state for $N\le 8$ and show that neither this state nor its fractionally charged excitations describe the physical reality. The results obtained shed new light on the nature of the ground and excited states in the fractional quantum Hall effect.


title: Preconditioned infinite GMRES for parameterized linear systems

abstract: We are interested in obtaining approximate solutions to parameterized linear systems of the form $A(\mu) x(\mu) = b$ for many values of the parameter $\mu$. Here $A(\mu)$ is large, sparse, and nonsingular, with a nonlinear analytic dependence on $\mu$. Our approach is based on a companion linearization for parameterized linear systems. The companion matrix is similar to the operator in the infinite Arnoldi method, and we use this to adapt the flexible GMRES setting. In this way, our method returns a function $\tilde{x}(\mu)$ which is cheap to evaluate for different $\mu$, and the preconditioner is applied only approximately. This novel approach leads to increased freedom to carry out the action of the operation inexactly, which provides performance improvement over the method infinite GMRES, without a loss of accuracy in general. We show that the error of our method is estimated based on the magnitude of the parameter $\mu$, the inexactness of the preconditioning, and the spectrum of the linear companion matrix. Numerical examples from a finite element discretization of a Helmholtz equation with a parameterized material coefficient illustrate the competitiveness of our approach. The simulations are reproducible and publicly available online.


title: Teacher Perception of Automatically Extracted Grammar Concepts for L2 Language Learning

abstract: One of the challenges of language teaching is how to organize the rules regarding syntax, semantics, or phonology of the language in a meaningful manner. This not only requires pedagogical skills, but also requires a deep understanding of that language. While comprehensive materials to develop such curricula are available in English and some broadly spoken languages, for many other languages, teachers need to manually create them in response to their students' needs. This process is challenging because i) it requires that such experts be accessible and have the necessary resources, and ii) even if there are such experts, describing all the intricacies of a language is time-consuming and prone to omission. In this article, we present an automatic framework that aims to facilitate this process by automatically discovering and visualizing descriptions of different aspects of grammar. Specifically, we extract descriptions from a natural text corpus that answer questions about morphosyntax (learning of word order, agreement, case marking, or word formation) and semantics (learning of vocabulary) and show illustrative examples. We apply this method for teaching the Indian languages, Kannada and Marathi, which, unlike English, do not have well-developed pedagogical resources and, therefore, are likely to benefit from this exercise. To assess the perceived utility of the extracted material, we enlist the help of language educators from schools in North America who teach these languages to perform a manual evaluation. Overall, teachers find the materials to be interesting as a reference material for their own lesson preparation or even for learner evaluation.


title: Local regularity for the space-homogeneous Landau equation with very soft potentials

abstract: This paper deals with the space-homogenous Landau equation with very soft potentials, including the Coulomb case. This nonlinear equation is of parabolic type with diffusion matrix given by the convolution product of the solution with the matrix $a_{ij} (z)=|z|^\gamma (|z|^2 \delta_{ij} - z_iz_j)$ for $\gamma \in [-3,-2)$. We derive local truncated entropy estimates and use them to establish two facts. Firstly, we prove that the set of singular points (in time and velocity) for the weak solutions constructed as in [C. Villani, Arch. Rational Mech. Anal. 143 (1998), 273-307] has zero $\mathscr{P}^{m_\ast}$ parabolic Hausdorff measure with $m_\ast:= \frac72 |2+\gamma|$. Secondly, we prove that if such a weak solution is axisymmetric, then it is smooth away from the symmetry axis. In particular, radially symmetric weak solutions are smooth away from the origin.


title: Nonparametric Identification of Kronecker Networks

abstract: We address the problem to estimate a dynamic network whose edges describe Granger causality relations and whose topology has a Kronecker structure. Such a structure arises in many real networks and allows to understand the organization of complex networks. We proposed a kernel-based PEM method to learn such networks. Numerical examples show the effectiveness of the proposed method.


title: Ground-state representation for fractional Laplacian on half-line

abstract: We give ground-state representation for the fractional Laplacian with Dirichlet condition on the half-line


title: MEAT: Maneuver Extraction from Agent Trajectories

abstract: Advances in learning-based trajectory prediction are enabled by large-scale datasets. However, in-depth analysis of such datasets is limited. Moreover, the evaluation of prediction models is limited to metrics averaged over all samples in the dataset. We propose an automated methodology that allows to extract maneuvers (e.g., left turn, lane change) from agent trajectories in such datasets. The methodology considers information about the agent dynamics and information about the lane segments the agent traveled along. Although it is possible to use the resulting maneuvers for training classification networks, we exemplary use them for extensive trajectory dataset analysis and maneuver-specific evaluation of multiple state-of-the-art trajectory prediction models. Additionally, an analysis of the datasets and an evaluation of the prediction models based on the agent dynamics is provided.


title: An Image Processing Pipeline for Camera Trap Time-Lapse Recordings

abstract: A new open-source image processing pipeline for analyzing camera trap time-lapse recordings is described. This pipeline includes machine learning models to assist human-in-the-loop video segmentation and animal re-identification. We present some performance results and observations on the utility of this pipeline after using it in a year-long project studying the spatial ecology and social behavior of the gopher tortoise.


title: Classes of Aggregation Rules for Ethical Decision Making in Automated Systems

abstract: We study a class of {\em aggregation rules} that could be applied to ethical AI decision-making. These rules yield the decisions to be made by automated systems based on the information of profiles of preferences over possible choices. We consider two different but very intuitive notions of preferences of an alternative over another one, namely {\it pairwise majority} and {\it position} dominance. Preferences are represented by permutation processes over alternatives and aggregation rules are applied to obtain results that are socially considered to be ethically correct. In this setting, we find many aggregation rules that satisfy desirable properties for an autonomous system. We also address the problem of the stability of the aggregation process, which is important when the information is variable. These results are a contribution for an AI designer that wants to justify the decisions made by an autonomous system.\\ \textit{Keywords:} Aggregation Operators; Permutation Process; Decision Analysis.


title: Approximating optimal SMC proposal distributions in individual-based epidemic models

abstract: Many epidemic models are naturally defined as individual-based models: where we track the state of each individual within a susceptible population. Inference for individual-based models is challenging due to the high-dimensional state-space of such models, which increases exponentially with population size. We consider sequential Monte Carlo algorithms for inference for individual-based epidemic models where we make direct observations of the state of a sample of individuals. Standard implementations, such as the bootstrap filter or the auxiliary particle filter are inefficient due to mismatch between the proposal distribution of the state and future observations. We develop new efficient proposal distributions that take account of future observations, leveraging the properties that (i) we can analytically calculate the optimal proposal distribution for a single individual given future observations and the future infection rate of that individual; and (ii) the dynamics of individuals are independent if we condition on their infection rates. Thus we construct estimates of the future infection rate for each individual, and then use an independent proposal for the state of each individual given this estimate. Empirical results show order of magnitude improvement in efficiency of the sequential Monte Carlo sampler for both SIS and SEIR models.


title: The Tur\'an number for the edge blow-up of trees: the missing case

abstract: The edge blow-up of a graph is the graph obtained from replacing each edge of it by a clique of the same size where the new vertices of the cliques are all different. Wang, Hou, Liu and Ma determined the Tur\'{a}n number of the edge blow-up of trees except one particular case. Answering an problem posed by them, we determined the Tur\'{a}n number of this particular case.


title: Trading Payoffs to Enlarged Neighborhoods? A New Evidence from Evolutionary Game Theory

abstract: Population diversity is an important aspect of Prisoner's Dilemma Game (PDG) research. However, the studies on dynamic diversity and its associated cost still need further investigation. Based on a framework comprising 2-dimensional spatial evolutionary PDG, this work examines the change in a player's neighborhood by enabling each player to pay for an upgrade of their neighborhood to switch from the von Neumann to Moore neighborhood. The upgrade cost (i.e., the cost of the advanced neighborhood) plays a vital role in cooperation promotion and serves as an entry-level to screen players. The results show that a reasonable price (entry-level) supports the cooperators' survival in an environment with high dilemma strength since it allows the formation of "normal-edge-advantage-core" clusters. On the low entry-level side, the privilege of having a larger neighborhood supports cooperation if it is accessible to all the players. On the high entry-level side, encirclements of advantage defectors appear out of the cooperative clusters. To break the encirclement and enable the expansion of the advantage clusters, the entry-level should be increased to interrupt the advantage defectors. The encirclement can be observed only in the deterministic models. Stochastic simulations are provided as robustness benchmarks.


title: On Scaling Laws for Multi-Well Nucleation Problems without Gauge Invariances

abstract: In this article we study scaling laws for simplified multi-well nucleation problems without gauge invariances which are motivated by models for shape-memory alloys. Seeking to explore the role of the order of lamination on the energy scaling for nucleation processes, we provide scaling laws for various model problems in two and three dimensions. In particular, we discuss (optimal) scaling results in the volume and the singular perturbation parameter for settings in which the surrounding parent phase is in the first, the second and the third order lamination convex hull of the wells of the nucleating phase. Furthermore, we provide a corresponding result for the setting of an infinite order laminate which arises in the context of the Tartar square. In particular, our results provide isoperimetric estimates in situations in which strong nonlocal anisotropies are present.


title: Multifidelity Reinforcement Learning with Control Variates

abstract: In many computational science and engineering applications, the output of a system of interest corresponding to a given input can be queried at different levels of fidelity with different costs. Typically, low-fidelity data is cheap and abundant, while high-fidelity data is expensive and scarce. In this work we study the reinforcement learning (RL) problem in the presence of multiple environments with different levels of fidelity for a given control task. We focus on improving the RL agent's performance with multifidelity data. Specifically, a multifidelity estimator that exploits the cross-correlations between the low- and high-fidelity returns is proposed to reduce the variance in the estimation of the state-action value function. The proposed estimator, which is based on the method of control variates, is used to design a multifidelity Monte Carlo RL (MFMCRL) algorithm that improves the learning of the agent in the high-fidelity environment. The impacts of variance reduction on policy evaluation and policy improvement are theoretically analyzed by using probability bounds. Our theoretical analysis and numerical experiments demonstrate that for a finite budget of high-fidelity data samples, our proposed MFMCRL agent attains superior performance compared with that of a standard RL agent that uses only the high-fidelity environment data for learning the optimal policy.


title: Multi-dimensional dual-blind deconvolution approach toward joint radar-communications

abstract: We consider a joint multiple-antenna radar-communications system in a co-existence scenario. Contrary to conventional applications, wherein at least the radar waveform and communications channel are known or estimated \textit{a priori}, we investigate the case when the channels and transmit signals of both systems are unknown. In radar applications, this problem arises in multistatic or passive systems, where transmit signal is not known. Similarly, highly dynamic vehicular or mobile communications may render prior estimates of wireless channel unhelpful. In particular, the radar signal reflected-off multiple targets is overlaid with the multi-carrier communications signal. In order to extract the unknown continuous-valued target parameters (range, Doppler velocity, and direction-of-arrival) and communications messages, we formulate the problem as a sparse dual-blind deconvolution and solve it using atomic norm minimization. Numerical experiments validate our proposed approach and show that precise estimation of continuous-valued channel parameters, radar waveform, and communications messages is possible up to scaling ambiguities.


title: A physical model and a Monte Carlo estimate for the specific intensity spatial derivative, angular derivative and geometric sensitivity

abstract: Starting from the radiative transfer equation and its usual boundary conditions, the objective of this work is to design Monte Carlo algorithms estimating the specific intensity spatial and angular derivatives as well as its geometric sensitivity. The present document is structured in three parts, each of them dedicated to a specific derivative of the intensity. Although they are all assembled here in one document each derivative is of interest independently whether it be for radiative transfers analysis or engineering conception. Therefore, they are thought to be written as three different papers and are presented here as such. Estimating derivatives of the specific intensity when solving radiative transfers using a Monte-Carlo algorithm is challenging. Finite differences are often not sufficiently accurate and directly estimating the derivative from a specific Monte-Carlo algorithm can lead to arduous formal or numerical developments. The proposition here is to work from the radiative transfer equation and its boundary conditions to design a physical model for each derivatives. Only then Monte-Carlo algorithms are built from the derivatives differential equations using the usual equivalent path integral. Since the same methodology is applied to the specific intensity spatial derivative, angular derivative and geometric sensitivity we chose to keep the same writing structure for all three parts so that all common ideas and developments appears exactly the same. We believe this choice to be coherent to facilitate the reader's understanding. Finally, these are preliminary versions of the final papers: for each parts the theory is fully described, but, although they have been implemented, the examples and algorithms sections are not always complete. This will be mentioned in the introductions of the concerned sections.


title: Multi-faceted Graph Attention Network for Radar Target Recognition in Heterogeneous Radar Network

abstract: Radar target recognition (RTR), as a key technology of intelligent radar systems, has been well investigated. Accurate RTR at low signal-to-noise ratios (SNRs) still remains an open challenge. Most existing methods are based on a single radar or the homogeneous radar network, which do not fully exploit frequency-dimensional information. In this paper, a two-stream semantic feature fusion model, termed Multi-faceted Graph Attention Network (MF-GAT), is proposed to greatly improve the accuracy in the low SNR region of the heterogeneous radar network. By fusing the features extracted from the source domain and transform domain via a graph attention network model, the MF-GAT model distills higher-level semantic features before classification in a unified framework. Extensive experiments are presented to demonstrate that the proposed model can greatly improve the RTR performance at low SNRs.


title: Bayesian calibration of coupled computational mechanics models under uncertainty based on interface deformation

abstract: Calibration or parameter identification is used with computational mechanics models related to observed data of the modeled process to find model parameters such that good similarity between model prediction and observation is achieved. We present a Bayesian calibration approach for surface coupled problems in computational mechanics based on measured deformation of an interface when no displacement data of material points is available. The interpretation of such a calibration problem as a statistical inference problem, in contrast to deterministic model calibration, is computationally more robust and allows the analyst to find a posterior distribution over possible solutions rather than a single point estimate. The proposed framework also enables the consideration of unavoidable uncertainties that are present in every experiment and are expected to play an important role in the model calibration process. To mitigate the computational costs of expensive forward model evaluations, we propose to learn the log-likelihood function from a controllable amount of parallel simulation runs using Gaussian process regression. We introduce and specifically study the effect of three different discrepancy measures for deformed interfaces between reference data and simulation. We show that a statistically based discrepancy measure results in the most expressive posterior distribution. We further apply the approach to numerical examples in higher model parameter dimensions and interpret the resulting posterior under uncertainty. In the examples, we investigate coupled multi-physics models of fluid-structure interaction effects in biofilms and find that the model parameters affect the results in a coupled manner.


title: Non-negative forms, volumes of sublevel sets, complete monotonicity and moment matrices

abstract: Let $\mathcal{C}_{d,n}$ be the convex cone consisting of real $n$-variate degree $d$ forms that are strictly positive on $\mathbb{R}^n\setminus \{\mathbf{0}\}$. We prove that the Lebesgue volume of the sublevel set $\{g\leq 1\}$ of $g\in \mathcal{C}_{d,n}$ is a completely monotone function on $\mathcal{C}_{d,n}$ and investigate the related properties. Furthermore, we provide (partial) characterization of forms, whose sublevel sets have finite Lebesgue volume. Finally, we discover an interesting property of a centered Gaussian distribution, establishing a connection between the matrix of its degree $d$ moments and the quadratic form given by the inverse of its covariance matrix.


title: Spectral analysis and fast methods for large matrices arising from PDE approximation

abstract: The main goal of this thesis is to show the crucial role that plays the symbol in analysing the spectrum the sequence of matrices resulting from PDE approximation and in designing a fast method to solve the associated linear problem. In the first part, we study the spectral properties of the matrices arising from $\mathbb{P}_k$ Lagrangian Finite Elements approximation of second order elliptic differential problem with Dirichlet boundary conditions and where the operator is $\mathrm{div} \left(-a(\mathbf{x}) \nabla\cdot\right)$, with $a$ continuous and positive over $\overline \Omega$, $\Omega$ being an open and bounded subset of $\mathbb{R}^d$, $d\ge 1$. We investigate the spectral distribution in the Weyl sense, with a concise overview on localization, clustering, extremal eigenvalues, and asymptotic conditioning. We study in detail the case of constant coefficients on $\Omega=(0,1)^2$ and we give a brief account in the case of variable coefficients and more general domains. While in the second part, we design a fast method of multigrid type for the resolution of linear systems arising from the $\mathbb{Q}_k$ Finite Elements approximation of the same considered problem in one and higher dimensional. The analysis is performed in one dimension, while the numerics are carried out also in higher dimension $d\ge 2$, demonstrating an optimal behavior in terms of the dependency on the matrix size and a robustness with respect to the dimensionality $d$ and to the polynomial degree $k$.


title: An Approach to Quantum 2D Gravity

abstract: We consider a model of 2D gravity with the action quadratic in curvature and represent path integrals as integrals over the SL(2, R) invariant Gaussian functional measure. We reduce these path integrals to the products of Wiener path integrals and calculate the correlation function of the metric in the first perturbative order.


title: How Much is Enough? A Study on Diffusion Times in Score-based Generative Models

abstract: Score-based diffusion models are a class of generative models whose dynamics is described by stochastic differential equations that map noise into data. While recent works have started to lay down a theoretical foundation for these models, an analytical understanding of the role of the diffusion time T is still lacking. Current best practice advocates for a large T to ensure that the forward dynamics brings the diffusion sufficiently close to a known and simple noise distribution; however, a smaller value of T should be preferred for a better approximation of the score-matching objective and higher computational efficiency. Starting from a variational interpretation of diffusion models, in this work we quantify this trade-off, and suggest a new method to improve quality and efficiency of both training and sampling, by adopting smaller diffusion times. Indeed, we show how an auxiliary model can be used to bridge the gap between the ideal and the simulated forward dynamics, followed by a standard reverse diffusion process. Empirical results support our analysis; for image data, our method is competitive w.r.t. the state-of-the-art, according to standard sample quality metrics and log-likelihood.


title: Near-Optimal Distributed Dominating Set in Bounded Arboricity Graphs

abstract: We describe a simple deterministic $O( \varepsilon^{-1} \log \Delta)$ round distributed algorithm for $(2\alpha+1)(1 + \varepsilon)$ approximation of minimum weighted dominating set on graphs with arboricity at most $\alpha$. Here $\Delta$ denotes the maximum degree. We also show a lower bound proving that this round complexity is nearly optimal even for the unweighted case, via a reduction from the celebrated KMW lower bound on distributed vertex cover approximation [Kuhn, Moscibroda, and Wattenhofer JACM'16]. Our algorithm improves on all the previous results (that work only for unweighted graphs) including a randomized $O(\alpha^2)$ approximation in $O(\log n)$ rounds [Lenzen and Wattenhofer DISC'10], a deterministic $O(\alpha \log \Delta)$ approximation in $O(\log \Delta)$ rounds [Lenzen and Wattenhofer DISC'10], a deterministic $O(\alpha)$ approximation in $O(\log^2 \Delta)$ rounds [implicit in Bansal and Umboh IPL'17 and Kuhn, Moscibroda, and Wattenhofer SODA'06], and a randomized $O(\alpha)$ approximation in $O(\alpha\log n)$ rounds [Morgan, Solomon and Wein DISC'21]. We also provide a randomized $O(\alpha \log\Delta)$ round distributed algorithm that sharpens the approximation factor to $\alpha(1+o(1))$. If each node is restricted to do polynomial-time computations, our approximation factor is tight in the first order as it is NP-hard to achieve $\alpha - 1 - \varepsilon$ approximation [Bansal and Umboh IPL'17].


title: A Causal Research Pipeline and Tutorial for Psychologists and Social Scientists

abstract: Causality is a fundamental part of the scientific endeavour to understand the world. Unfortunately, causality is still taboo in much of psychology and social science. Motivated by a growing number of recommendations for the importance of adopting causal approaches to research, we reformulate the typical approach to research in psychology to harmonize inevitably causal theories with the rest of the research pipeline. We present a new process which begins with the incorporation of techniques from the confluence of causal discovery and machine learning for the development, validation, and transparent formal specification of theories. We then present methods for reducing the complexity of the fully specified theoretical model into the fundamental submodel relevant to a given target hypothesis. From here, we establish whether or not the quantity of interest is estimable from the data, and if so, propose the use of semi-parametric machine learning methods for the estimation of causal effects. The overall goal is the presentation of a new research pipeline which can (a) facilitate scientific inquiry compatible with the desire to test causal theories (b) encourage transparent representation of our theories as unambiguous mathematical objects, (c) to tie our statistical models to specific attributes of the theory, thus reducing under-specification problems frequently resulting from the theory-to-model gap, and (d) to yield results and estimates which are causally meaningful and reproducible. The process is demonstrated through didactic examples with real-world data, and we conclude with a summary and discussion of limitations.


title: Phenomenology of an in-host model of hepatitis C

abstract: This paper carries out an analysis of the global properties of solutions of an in-host model of hepatitis C for general values of its parameters. A previously unknown stable steady state on the boundary of the positive orthant is exhibited. It is proved that the model exhibits Hopf bifurcations and hence periodic solutions. A general parametrization of positive steady states is given and it is determined when the number of steady states is odd or even, according to the value of a certain basic reproductive ratio. This implies, in particular, that when this reproductive ratio is greater than one there always exists at least one positive steady state. A positive steady state which bifurcates from an infection-free state when the reproductive ratio passes through one is always stable, i.e. no backward bifurcation occurs in this model.


title: $m$-Symmetric functions, non-symmetric Macdonald polynomials and positivity conjectures

abstract: We study the space, $R_m$, of $m$-symmetric functions consisting of polynomials that are symmetric in the variables $x_{m+1},x_{m+2},x_{m+3},\dots$ but have no special symmetry in the variables $x_1,\dots,x_m$. We obtain $m$-symmetric Macdonald polynomials by $t$-symmetrizing non-symmetric Macdonald polynomials, and show that they form a basis of $R_m$. We define $m$-symmetric Schur functions through a somewhat complicated process involving their dual basis, tableaux combinatorics, and the Hecke algebra generators, and then prove some of their most elementary properties. We conjecture that the $m$-symmetric Macdonald polynomials (suitably normalized and plethystically modified) expand positively in terms of $m$-symmetric Schur functions. We obtain relations on the $(q,t)$-Koska coefficients $K_{\Omega \Lambda}(q,t)$ in the $m$-symmetric world, and show in particular that the usual $(q,t)$-Koska coefficients are special cases of the $K_{\Omega \Lambda}(q,t)$'s. Finally, we show that when $m$ is large, the positivity conjecture, modulo a certain subspace, becomes a positivity conjecture on the expansion of non-symmetric Macdonald polynomials in terms of non-symmetric Hall-Littlewood polynomials.


title: Observing high-k magnons with Mie-resonance-enhanced Brillouin light scattering

abstract: Magnonics is a prospective beyond CMOS technology which uses magnons, the quanta of spin waves, for low-power information processing. Many magnonic concepts and devices were recently demonstrated at macro- and microscale, and now these concepts need to be realized at nanoscale. Brillouin light scattering spectroscopy and microscopy (BLS) has become a standard technique for spin wave visualization and characterization, and enabled many pioneering magnonic experiments. However, due to its fundamental limit in maximum detectable magnon momentum, the conventional BLS cannot be used to detect nanoscale spin waves. Here we show that optically induced Mie resonances in dielectric nanoparticles can be used to extend the range of accessible spin wave wavevectors beyond the BLS fundamental limit. The method is universal and can be used in many magnonic experiments dealing with thermally excited as well as coherently excited high-momentum, short-wavelength spin waves. This discovery significantly extends the usability and relevance of the BLS technique for nanoscale magnonic research.

categories: cond-mat.mes-hall cond-mat.other physics.app-ph physics.optics


title: Using stellar rotation to identify tidally stripped members of the Praesepe open cluster

abstract: As an open cluster orbits the Milky Way, gravitational fields distort it, stripping stars from the core and forming tidal tails. Recent work has identified tidal tails of the Praesepe cluster; we explore rotation periods as a way to confirm these candidate members. In open clusters, the rotation period distribution evolves over time due to magnetic braking. Since tidally stripped stars originally formed within the cluster, they should follow the same period distribution as in the cluster core. We analyze 96 candidate members observed by NASA's Transiting Exoplanet Survey Satellite (TESS) mission. We measure reliable rotation periods for 32 stars, while 64 light curves are noise-dominated. The 32 newly identified rotators are consistent with the period distribution in the core, and with past membership in Praesepe. We therefore suggest that for nearby open clusters, stellar rotation offers a quick and inexpensive method for confirming past members dispersed into tidal tails.


title: A Stochastic Approach to Maxwell Velocity Distribution via Central Limit Theorem and Boltzmann's Entropy Formula

abstract: Maxwell's velocity distribution is known to be universally valid across systems and phases. Here we present a new and general derivation that uses the central limit theorem (CLT) of the probability theory. This essentially uses the idea that repeated intermolecular collisions introduce randomness in the velocity change in the individual components of the velocity vector, leading to, by the CLT, a Gaussian distribution. To complete the derivation, we next show that the mean-square velocity or the standard deviation follows exactly from Boltzmann's formula relating entropy to the density of states, thus avoiding the use of the ideal gas equation of state explicitly. We furthermore prove the Maxwell velocity distribution for a system of harmonic oscillators. This derivation provides a further insight into the origin of Boltzmann's constant in the Maxwell velocity distribution and also in the equipartition theorem. We propose that this derivation provides an approach that explains the universality of Maxwell Gaussian distribution.


title: Lightweight Conditional Model Extrapolation for Streaming Data under Class-Prior Shift

abstract: We introduce LIMES, a new method for learning with non-stationary streaming data, inspired by the recent success of meta-learning. The main idea is not to attempt to learn a single classifier that would have to work well across all occurring data distributions, nor many separate classifiers, but to exploit a hybrid strategy: we learn a single set of model parameters from which a specific classifier for any specific data distribution is derived via classifier adaptation. Assuming a multi-class classification setting with class-prior shift, the adaptation step can be performed analytically with only the classifier's bias terms being affected. Another contribution of our work is an extrapolation step that predicts suitable adaptation parameters for future time steps based on the previous data. In combination, we obtain a lightweight procedure for learning from streaming data with varying class distribution that adds no trainable parameters and almost no memory or computational overhead compared to training a single model. Experiments on a set of exemplary tasks using Twitter data show that LIMES achieves higher accuracy than alternative approaches, especially with respect to the relevant real-world metric of lowest within-day accuracy.


title: Human-AI Interaction Design in Machine Teaching

abstract: Machine Teaching (MT) is an interactive process where a human and a machine interact with the goal of training a machine learning model (ML) for a specified task. The human teacher communicates their task expertise and the machine student gathers the required data and knowledge to produce an ML model. MT systems are developed to jointly minimize the time spent on teaching and the learner's error rate. The design of human-AI interaction in an MT system not only impacts the teaching efficiency, but also indirectly influences the ML performance by affecting the teaching quality. In this paper, we build upon our previous work where we proposed an MT framework with three components, viz., the teaching interface, the machine learner, and the knowledge base, and focus on the human-AI interaction design involved in realizing the teaching interface. We outline design decisions that need to be addressed in developing an MT system beginning from an ML task. The paper follows the Socratic method entailing a dialogue between a curious student and a wise teacher.


title: GD-VAEs: Geometric Dynamic Variational Autoencoders for Learning Nonlinear Dynamics and Dimension Reductions

abstract: We develop data-driven methods incorporating geometric and topological information to learn parsimonious representations of nonlinear dynamics from observations. We develop approaches for learning nonlinear state space models of the dynamics for general manifold latent spaces using training strategies related to Variational Autoencoders (VAEs). Our methods are referred to as Geometric Dynamic (GD) Variational Autoencoders (GD-VAEs). We learn encoders and decoders for the system states and evolution based on deep neural network architectures that include general Multilayer Perceptrons (MLPs), Convolutional Neural Networks (CNNs), and Transpose CNNs (T-CNNs). Motivated by problems arising in parameterized PDEs and physics, we investigate the performance of our methods on tasks for learning low dimensional representations of the nonlinear Burgers equations, constrained mechanical systems, and spatial fields of reaction-diffusion systems. GD-VAEs provide methods for obtaining representations for use in diverse learning tasks involving dynamics.

categories: cs.LG cs.NA math.DS math.NA physics.data-an stat.ML


title: SERE: Exploring Feature Self-relation for Self-supervised Transformer

abstract: Learning representations with self-supervision for convolutional networks (CNN) has been validated to be effective for vision tasks. As an alternative to CNN, vision transformers (ViT) have strong representation ability with spatial self-attention and channel-level feedforward networks. Recent works reveal that self-supervised learning helps unleash the great potential of ViT. Still, most works follow self-supervised strategies designed for CNN, e.g., instance-level discrimination of samples, but they ignore the properties of ViT. We observe that relational modeling on spatial and channel dimensions distinguishes ViT from other networks. To enforce this property, we explore the feature SElf-RElation (SERE) for training self-supervised ViT. Specifically, instead of conducting self-supervised learning solely on feature embeddings from multiple views, we utilize the feature self-relations, i.e., spatial/channel self-relations, for self-supervised learning. Self-relation based learning further enhances the relation modeling ability of ViT, resulting in stronger representations that stably improve performance on multiple downstream tasks. Our source code is publicly available at: https://github.com/MCG-NKU/SERE.


title: Decision-making under uncertainty: a quantum value operator approach

abstract: We propose a quantum expected value theory for decision-making under uncertainty. Quantum density operator as value operator is proposed to simulate people's subjective beliefs. Value operator guides people to choose corresponding actions based on their subjective beliefs through objective world. The value operator can be constructed from quantum gates and logic operations as a quantum decision tree. The genetic programming is used to optimize and auto-generate quantum decision trees.


title: Improved cosmological bounds for a fine-tuned see-saw mechanism of keV sterile neutrinos

abstract: Considering the gauge model SU(2)_L x U(1) with the extension of the lepton sector by right-handed sterile Majorana neutrinos, we distinguish a "minimal parametric mixing scenario" for a complete matrix of 6 x 6 light (active) and sterile neutrinos, which depends, in addition to experimentally measured physical parameters, only on the masses of sterile neutrinos. For such a scenario, improved cosmological bounds are provided, resulting from the lifetime of sterile neutrinos and the fraction of energy carried by sterile neutrino dark matter.


title: On Convergence of FedProx: Local Dissimilarity Invariant Bounds, Non-smoothness and Beyond

abstract: The FedProx algorithm is a simple yet powerful distributed proximal point optimization method widely used for federated learning (FL) over heterogeneous data. Despite its popularity and remarkable success witnessed in practice, the theoretical understanding of FedProx is largely underinvestigated: the appealing convergence behavior of FedProx is so far characterized under certain non-standard and unrealistic dissimilarity assumptions of local functions, and the results are limited to smooth optimization problems. In order to remedy these deficiencies, we develop a novel local dissimilarity invariant convergence theory for FedProx and its minibatch stochastic extension through the lens of algorithmic stability. As a result, we contribute to derive several new and deeper insights into FedProx for non-convex federated optimization including: 1) convergence guarantees independent on local dissimilarity type conditions; 2) convergence guarantees for non-smooth FL problems; and 3) linear speedup with respect to size of minibatch and number of sampled devices. Our theory for the first time reveals that local dissimilarity and smoothness are not must-have for FedProx to get favorable complexity bounds. Preliminary experimental results on a series of benchmark FL datasets are reported to demonstrate the benefit of minibatching for improving the sample efficiency of FedProx.


title: Splitted Levenberg-Marquardt Method for Large-Scale Sparse Problems

abstract: We consider large-scale nonlinear least squares problems with sparse residuals, each of them depending on a small number of variables. A decoupling procedure which results in a splitting of the original problems into a sequence of independent problems of smaller sizes is proposed and analysed. The smaller size problems are modified in a way that offsets the error made by disregarding dependencies that allow us to split the original problem. The resulting method is a modification of the Levenberg-Marquardt method with smaller computational costs. Global convergence is proved as well as local linear convergence under suitable assumptions on sparsity. The method is tested on the network localization simulated problems with up to one million variables and its efficiency is demonstrated.


title: Alternative computation of the Seidel aberration coefficients using the Lie algebraic method

abstract: We give a brief introduction to Hamiltonian optics and Lie algebraic methods. We use these methods to describe the operators governing light propagation, refraction and reflection in phase space. The method offers a systematic way to find aberration coefficients of any order for arbitrary rotationally symmetric optical systems. The coefficients from the Lie method are linked to the Seidel aberration coefficients. Furthermore, the property of summing individual surface contributions is preserved by the Lie algebraic theory. Three examples are given to validate the proposed methodology with good results.


title: Heterotic de-Sitter Solutions

abstract: We classify all warped product de-Sitter backgrounds in heterotic supergravity, up to two loops. We find that warped $dS_n$ backgrounds, with $n\ge 3$, are $\mathbb{R}^{1,n}\times M_{9-n}$, where $M_{9-n}$ is a $(9-n)$-dimensional Riemannian manifold. Moreover, we establish that warped $dS_2$ backgrounds are $AdS_3\times M_7$, where $M_7$ is a 7-dimensional Riemannian manifold.


title: Fitting covariance matrix models to simulations

abstract: Data analysis in cosmology requires reliable covariance matrices. Covariance matrices derived from numerical simulations often require a very large number of realizations to be accurate. When a theoretical model for the covariance matrix exists, the parameters of the model can often be fit with many fewer simulations. We write a likelihood-based method for performing such a fit. We demonstrate how a model covariance matrix can be tested by examining the appropriate $\chi^2$ distributions from simulations. We show that if model covariance has amplitude freedom, the expectation value of second moment of $\chi^2$ distribution with a wrong covariance matrix will always be larger than one using the true covariance matrix. By combining these steps together, we provide a way of producing reliable covariances without ever requiring running a large number of simulations. We demonstrate our method on two examples. First, we measure the two-point correlation function of halos from a large set of $10000$ mock halo catalogs. We build a model covariance with $2$ free parameters, which we fit using our procedure. The resulting best-fit model covariance obtained from just $100$ simulation realizations proves to be as reliable as the numerical covariance matrix built from the full $10000$ set. We also test our method on a setup where the covariance matrix is large by measuring the halo bispectrum for thousands of triangles for the same set of mocks. We build a block diagonal model covariance with $2$ free parameters as an improvement over the diagonal Gaussian covariance. Our model covariance passes the $\chi^2$ test only partially in this case, signaling that the model is insufficient even using free parameters, but significantly improves over the Gaussian one.


title: Observation of $\tau$ lepton pair production in ultraperipheral lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV

abstract: We present an observation of photon-photon production of $\tau$ lepton pairs in ultraperipheral lead-lead collisions. The measurement is based on a data sample with an integrated luminosity of 404 $\mu$b$^{-1}$ collected by the CMS experiment at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The $\gamma\gamma$ $\to$ $\tau^+\tau^-$ process is observed for $\tau\tau$ events with a muon and three charged hadrons in the final state. The measured fiducial cross section is $\sigma(\gamma\gamma$ $\to$ $\tau^+\tau^-)$ = 4.8 $\pm$ 0.6 (stat) $\pm$ 0.5 (syst) $\mu$b, in agreement with leading-order QED predictions. Using $\sigma(\gamma\gamma$ $\to$ $\tau^+\tau^-)$, we estimate a model-dependent value of the anomalous magnetic moment of the $\tau$ lepton of $a_\tau$ = 0.001 $^{+0.055}_{-0.089}$.


title: A Monte Carlo approach to the conformal bootstrap

abstract: We introduce an approach to find approximate numerical solutions of truncated bootstrap equations for Conformal Field Theories (CFTs) in arbitrary dimensions. The method is based on a stochastic search via a Metropolis algorithm guided by an action $S$ which is the logarithm of the truncated bootstrap equations for a single scalar field correlator. While numerical conformal bootstrap methods based on semi-definite programming put rigorous exclusion bounds on CFTs, this method looks for approximate solutions, which correspond to local minima of $S$, when present, and can be even far from the extremality region. By this protocol we find that if no constraint on the operator scaling dimensions is imposed, $S$ has a single minimum, corresponding to the Free Theory. If we fix the external operator dimension, however, we encounter minima that can be studied with our approach. Imposing a conserved stress-tensor, a $\mathbf{Z}_2$ symmetry and one relevant scalar, we identify two regions where local minima of $S$ are present. When projected in the $(\Delta_\sigma, \Delta_{\epsilon})$-plane, $\sigma$ and $\epsilon$ being the external and the lightest exchanged operators, one of these regions essentially coincides with the extremality line found in previous bootstrap studies. The other region is along the generalized free theories in $d = 2$ and below that in both $d = 3$ and $d = 4$. We empirically prove that some of the minima found are associated to known theories, including the $2d$ and $3d$ Ising theories and the $2d$ Yang-Lee model.