title: Resonant Interaction of Modulation-correlated Quantum Electron Wavepackets with Bound Electron States

abstract: Free-Electron Bound-Electron Resonant Interaction (FEBERI) is the resonant inelastic interaction of periodically density-bunched free electrons with a quantum two level system. We present a comprehensive relativistic quantum mechanical theory for this interaction in a model in which the electrons are represented as quantum electron wavepackets (QEW). The analysis reveals the wave-particle duality nature of the QEW, delineating the point-particle-like and wave-like interaction regimes, and manifesting the physical reality of the wavefunction dimensions and its density modulation characteristics in interaction with matter. The analysis comprehends the case of laser-beam-modulated multiple QEWs that are modulation-phase correlated. Based on the Born interpretation of the electron wavefunction we predict quantum transitions enhancement proportional to the number of electrons squared, analogous to superradiance.


title: A deep neural network algorithm for semilinear elliptic PDEs with applications in insurance mathematics

abstract: In insurance mathematics optimal control problems over an infinite time horizon arise when computing risk measures. Their solutions correspond to solutions of deterministic semilinear (degenerate) elliptic partial differential equations. In this paper we propose a deep neural network algorithm for solving such partial differential equations in high dimensions. The algorithm is based on the correspondence of elliptic partial differential equations to backward stochastic differential equations with random terminal time.


title: Diameters of Graphs on Reduced Words of 12 and 21-Inflations

abstract: It is a classical result that any permutation in the symmetric group can be generated by a sequence of adjacent transpositions. The sequences of minimal length are called reduced words, and in this paper we study the graphs of these reduced words, with edges determined by relations in the underlying Coxeter group. Recently, the diameter has been calculated for the longest permutation $n\ldots 21$ by Reiner and Roichman as well as Assaf. In this paper we find inductive formulas for the diameter of the graphs of 12-inflations and many 21-inflations. These results extend to the associated graphs on commutation and long braid classes. Also, these results give a recursive formula for the diameter of the longest permutation, which matches that of Reiner, Roichman and Assaf. Lastly, We make progress on conjectured bounds of the diameter by Reiner and Roichman, which are based on the underlying hyperplane arrangement, and find families of permutations that achieve the upper bound and lower bound of the conjecture. In particular permutations that avoid 312 or 231 have graphs that achieve the upper bound.


title: Entanglement Entropy for $T \bar T$, $J \bar T$, $T \bar J$ deformed holographic CFT

abstract: We derive the geodesic equation for determining the Ryu-Takayanagi surface in $AdS_3$ deformed by single trace $\mu T \bar T + \varepsilon_+ J \bar T + \varepsilon_- T \bar J$ deformation for generic values of $(\mu, \varepsilon_+, \varepsilon_-)$ for which the background is free of singularities. For generic values of $\varepsilon_\pm$, Lorentz invariance is broken, and the Ryu-Takayanagi surface embeds non-trivially in time as well as spatial coordinates. We solve the geodesic equation and characterize the UV and IR behavior of the entanglement entropy and the Casini-Huerta $c$-function. We comment on various features of these observables in the $(\mu, \varepsilon_+, \varepsilon_-)$ parameter space. We discuss the matching at leading order in small $(\mu, \varepsilon_+, \varepsilon_-)$ expansion of the entanglement entropy between the single trace deformed holographic system and a class of double trace deformed theories where a strictly field theoretic analysis is possible. We also comment on expectation value of a large rectangular Wilson loop-like observable.


title: Identifying Transition States of Chemical Kinetic Systems using Network Embedding Techniques

abstract: Using random walk sampling methods for feature learning on networks, we develop a method for generating low-dimensional node embeddings for directed graphs and identifying transition states of stochastic chemical reacting systems. We modified objective functions adopted in existing random walk based network embedding methods to handle directed graphs and neighbors of different degrees. Through optimization via gradient ascent, we embed the weighted graph vertices into a low-dimensional vector space Rd while preserving the neighborhood of each node. We then demonstrate the effectiveness of the method on dimension reduction through several examples regarding identification of transition states of chemical reactions, especially for entropic systems.


title: A Helmholtz equation solver using unsupervised learning: Application to transcranial ultrasound

abstract: Transcranial ultrasound therapy is increasingly used for the non-invasive treatment of brain disorders. However, conventional numerical wave solvers are currently too computationally expensive to be used online during treatments to predict the acoustic field passing through the skull (e.g., to account for subject-specific dose and targeting variations). As a step towards real-time predictions, in the current work, a fast iterative solver for the heterogeneous Helmholtz equation in 2D is developed using a fully-learned optimizer. The lightweight network architecture is based on a modified UNet that includes a learned hidden state. The network is trained using a physics-based loss function and a set of idealized sound speed distributions with fully unsupervised training (no knowledge of the true solution is required). The learned optimizer shows excellent performance on the test set, and is capable of generalization well outside the training examples, including to much larger computational domains, and more complex source and sound speed distributions, for example, those derived from x-ray computed tomography images of the skull.

categories: physics.comp-ph cs.LG eess.IV physics.med-ph


title: The intensity and evolution of the extreme storms in January 1938

abstract: Major solar eruptions occasionally direct interplanetary coronal mass ejections (ICMEs) to Earth and cause significant geomagnetic storms and low-latitude aurorae. While single extreme storms are of significant threats to the modern civilization, storms occasionally appear in sequence and, acting synergistically, cause 'perfect storms' at Earth. The stormy interval in January 1938 was one of such cases. Here, we analyze the contemporary records to reveal its time series on their source active regions, solar eruptions, ICMEs, geomagnetic storms, low-latitude aurorae, and cosmic-ray (CR) variations. Geomagnetic records show that three storms occurred successively on 17/18 January (Dcx ~ -171 nT) on 21/22 January (Dcx ~ -328 nT) and on 25/26 January (Dcx ~ -336 nT). The amplitudes of the cosmic-ray variations and sudden storm commencements show the impact of the first ICME as the largest (~ 6% decrease in CR and 72 nT in SSC) and the ICME associated with the storms that followed as more moderate (~ 3% decrease in CR and 63 nT in SSC; ~ 2% decrease in CR and 63 nT in SSC). Interestingly, a significant solar proton event occurred on 16/17 January and the Cheltenham ionization chamber showed a possible ground level enhancement. During the first storm, aurorae were less visible at mid-latitudes, whereas during the second and third storms, the equatorward boundaries of the auroral oval were extended down to 40.3{\deg} and 40.0{\deg} in invariant latitude. This contrast shows that the initial ICME was probably faster, with a higher total magnitude but a smaller southward component.


title: Theory of optically detected spin noise in nanosystems

abstract: Theory of spin noise in low dimensional systems and bulk semiconductors is reviewed. Spin noise is usually detected by optical means, continuously measuring the rotation angle of the polarization plane of the probe beam passing through the sample. Spin noise spectra yield rich information about the spin properties of the system including, for example, $g$-factors of the charge carriers, spin relaxation times, parameters of the hyperfine interaction, spin-orbit interaction constants, frequencies and widths of the optical resonances. The review describes basic models of spin noise, methods of its theoretical description, and their relation with the experimental results. We also discuss the relation between the spin noise spectroscopy, the strong and weak quantum measurements and the spin flip Raman scattering, and analyze similar effects including manifestations of the charge, current and valley polarization fluctuations in the optical response. Possible directions for further development of the spin noise spectroscopy are outlined.


title: Domain adaptation under structural causal models

abstract: Domain adaptation (DA) arises as an important problem in statistical machine learning when the source data used to train a model is different from the target data used to test the model. Recent advances in DA have mainly been application-driven and have largely relied on the idea of a common subspace for source and target data. To understand the empirical successes and failures of DA methods, we propose a theoretical framework via structural causal models that enables analysis and comparison of the prediction performance of DA methods. This framework also allows us to itemize the assumptions needed for the DA methods to have a low target error. Additionally, with insights from our theory, we propose a new DA method called CIRM that outperforms existing DA methods when both the covariates and label distributions are perturbed in the target data. We complement the theoretical analysis with extensive simulations to show the necessity of the devised assumptions. Reproducible synthetic and real data experiments are also provided to illustrate the strengths and weaknesses of DA methods when parts of the assumptions in our theory are violated.


title: Optimal bounds for the colorful fractional Helly theorem

abstract: The well known fractional Helly theorem and colorful Helly theorem can be merged into the so called colorful fractional Helly theorem. It states: For every $\alpha \in (0, 1]$ and every non-negative integer $d$, there is $\beta_{col} = \beta_{col}(\alpha, d) \in (0, 1]$ with the following property. Let $\mathcal{F}_1, \dots, \mathcal{F}_{d+1}$ be finite nonempty families of convex sets in $\mathbb{R}^d$ of sizes $n_1, \dots, n_{d+1}$ respectively. If at least $\alpha n_1 n_2 \cdots n_{d+1}$ of the colorful $(d+1)$-tuples have a nonempty intersection, then there is $i \in [d+1]$ such that $\mathcal{F}_i$ contains a subfamily of size at least $\beta_{col} n_i$ with a nonempty intersection. (A colorful $(d+1)$-tuple is a $(d+1)$-tuple $(F_1, \dots , F_{d+1})$ such that $F_i$ belongs to $\mathcal{F}_i$ for every $i$.) The colorful fractional Helly theorem was first stated and proved by B\'ar\'any, Fodor, Montejano, Oliveros, and P\'or in 2014 with $\beta_{col} = \alpha/(d+1)$. In 2017 Kim proved the theorem with better function $\beta_{col}$, which in particular tends to $1$ when $\alpha$ tends to $1$. Kim also conjectured what is the optimal bound for $\beta_{col}(\alpha, d)$ and provided the upper bound example for the optimal bound. The conjectured bound coincides with the optimal bounds for the (non-colorful) fractional Helly theorem proved independently by Eckhoff and Kalai around 1984. We verify Kim's conjecture by extending Kalai's approach to the colorful scenario. Moreover, we obtain optimal bounds also in more general setting when we allow several sets of the same color.


title: Global higher integrability for minimisers of convex functionals with (p,q)-growth

abstract: We prove global $W^{1,q}(\Omega,\mathbb{R}^m)$-regularity for minimisers of convex functionals of the form $\mathscr{F}(u)=\int_\Omega F(x,Du)\mathrm{d} x$. $W^{1,q}(\Omega,\mathbb{R}^m)$ regularity is also proven for minimisers of the associated relaxed functional. Our main assumptions on $F(x,z)$ are a uniform $\alpha$-H\"older continuity assumption in $x$ and controlled $(p,q)$-growth conditions in $z$ with $q<\frac{(n+\alpha)p}{n}$.


title: Cosmological accumulation of conserved numbers in Kaluza-Klein theories

abstract: We develop a new mechanism for the accumulation of conserved numbers in the early Universe in Kaluza-Klein-like theories. The relaxation of the primordial extra space perturbations existing in the early Universe leads to the establishment of a symmetric final state and the appearance of Killing vectors. As a result, the initial non-zero value of symmetry associated numbers occurs after the inflation. We show this conceptual idea on a toy model of 2-dimensional apple-like extra space with U(1) symmetry. This mechanism naturally arises in the Kaluza-Klein theories and can be used to explain the observed cosmological baryon asymmetry.


title: A Single-Loop Smoothed Gradient Descent-Ascent Algorithm for Nonconvex-Concave Min-Max Problems

abstract: Nonconvex-concave min-max problem arises in many machine learning applications including minimizing a pointwise maximum of a set of nonconvex functions and robust adversarial training of neural networks. A popular approach to solve this problem is the gradient descent-ascent (GDA) algorithm which unfortunately can exhibit oscillation in case of nonconvexity. In this paper, we introduce a "smoothing" scheme which can be combined with GDA to stabilize the oscillation and ensure convergence to a stationary solution. We prove that the stabilized GDA algorithm can achieve an $O(1/\epsilon^2)$ iteration complexity for minimizing the pointwise maximum of a finite collection of nonconvex functions. Moreover, the smoothed GDA algorithm achieves an $O(1/\epsilon^4)$ iteration complexity for general nonconvex-concave problems. Extensions of this stabilized GDA algorithm to multi-block cases are presented. To the best of our knowledge, this is the first algorithm to achieve $O(1/\epsilon^2)$ for a class of nonconvex-concave problem. We illustrate the practical efficiency of the stabilized GDA algorithm on robust training.


title: An Algebraic Foundation of Amended Dimensional Analysis

abstract: We present an innovative approach to dimensional analysis, based on a general representation theorem for complete quantity functions admitting a covariant scalar representation; this theorem is in turn grounded in a purely algebraic theory of quantity spaces. Examples of dimensional analysis based on this approach are given, showing that it allows results obtained by traditional dimensional analysis to be strengthened. For example, the orbital period of a two-body system can be derived without use of equations of motion.


title: Recursive Random Contraction Revisited

abstract: In this note, we revisit the recursive random contraction algorithm of Karger and Stein for finding a minimum cut in a graph. Our revisit is occasioned by a paper of Fox, Panigrahi, and Zhang which gives an extension of the Karger-Stein algorithm to minimum cuts and minimum $k$-cuts in hypergraphs. When specialized to the case of graphs, the algorithm is somewhat different than the original Karger-Stein algorithm. We show that the analysis becomes particularly clean in this case: we can prove that the probability that a fixed minimum cut in an $n$ node graph is returned by the algorithm is bounded below by $1/(2H_n-2)$, where $H_n$ is the $n$th harmonic number. We also consider other similar variants of the algorithm, and show that no such algorithm can achieve an asymptotically better probability of finding a fixed minimum cut.


title: Dynamics of gravitational collapse in the axisymmetric Einstein-Vlasov system

abstract: We numerically investigate the dynamics near black hole formation of solutions to the Einstein--Vlasov system in axisymmetry. Our results are obtained using a particle-in-cell and finite difference code based on the $(2+1)+1$ formulation of the Einstein field equations in axisymmetry. Solutions are launched from non-stationary initial data and exhibit type I critical behaviour. In particular, we find lifetime scaling in solutions containing black holes, and support that the critical solutions are stationary. Our results contain examples of solutions that form black holes, perform damped oscillations, and appear to disperse. We prove that complete dispersal of the solution implies that it has nonpositive binding energy.


title: GANs & Reels: Creating Irish Music using a Generative Adversarial Network

abstract: In this paper we present a method for algorithmic melody generation using a generative adversarial network without recurrent components. Music generation has been successfully done using recurrent neural networks, where the model learns sequence information that can help create authentic sounding melodies. Here, we use DC-GAN architecture with dilated convolutions and towers to capture sequential information as spatial image information, and learn long-range dependencies in fixed-length melody forms such as Irish traditional reel.


title: WaveTransform: Crafting Adversarial Examples via Input Decomposition

abstract: Frequency spectrum has played a significant role in learning unique and discriminating features for object recognition. Both low and high frequency information present in images have been extracted and learnt by a host of representation learning techniques, including deep learning. Inspired by this observation, we introduce a novel class of adversarial attacks, namely `WaveTransform', that creates adversarial noise corresponding to low-frequency and high-frequency subbands, separately (or in combination). The frequency subbands are analyzed using wavelet decomposition; the subbands are corrupted and then used to construct an adversarial example. Experiments are performed using multiple databases and CNN models to establish the effectiveness of the proposed WaveTransform attack and analyze the importance of a particular frequency component. The robustness of the proposed attack is also evaluated through its transferability and resiliency against a recent adversarial defense algorithm. Experiments show that the proposed attack is effective against the defense algorithm and is also transferable across CNNs.


title: High-pressure synthesis of Dirac materials: layered van der Waals bonded BeN$_4$ polymorph

abstract: High pressure chemistry is known to inspire the creation of unexpected new classes of compounds with exceptional properties. Here we report the synthesis at ~90 GPa of novel beryllium polynitrides, monoclinic and triclinic BeN4. The triclinic phase, upon decompression to ambient conditions, transforms into a compound with atomic-thick BeN4 layers interconnected via weak van der Waals bonds consisting of polyacetylene-like nitrogen chains with conjugated {\pi}-systems and Be atoms in square-planar coordination. Theoretical calculations for a single BeN4 layer show that its electronic lattice is described by a slightly distorted honeycomb structure reminiscent of the graphene lattice and the presence of Dirac points in the electronic band structure at the Fermi level. The BeN4 layer, i.e. beryllonitrene, represents a qualitatively new class of 2D materials that can be built of a metal atom and polymeric nitrogen chains and host anisotropic Dirac fermions.


title: Understanding the Failure Modes of Out-of-Distribution Generalization

abstract: Empirical studies suggest that machine learning models often rely on features, such as the background, that may be spuriously correlated with the label only during training time, resulting in poor accuracy during test-time. In this work, we identify the fundamental factors that give rise to this behavior, by explaining why models fail this way {\em even} in easy-to-learn tasks where one would expect these models to succeed. In particular, through a theoretical study of gradient-descent-trained linear classifiers on some easy-to-learn tasks, we uncover two complementary failure modes. These modes arise from how spurious correlations induce two kinds of skews in the data: one geometric in nature, and another, statistical in nature. Finally, we construct natural modifications of image classification datasets to understand when these failure modes can arise in practice. We also design experiments to isolate the two failure modes when training modern neural networks on these datasets.


title: Quantum advantage for differential equation analysis

abstract: Quantum algorithms for both differential equation solving and for machine learning potentially offer an exponential speedup over all known classical algorithms. However, there also exist obstacles to obtaining this potential speedup in useful problem instances. The essential obstacle for quantum differential equation solving is that outputting useful information may require difficult post-processing, and the essential obstacle for quantum machine learning is that inputting the training set is a difficult task just by itself. In this paper, we demonstrate, when combined, these difficulties solve one another. We show how the output of quantum differential equation solving can serve as the input for quantum machine learning, allowing dynamical analysis in terms of principal components, power spectra, and wavelet decompositions. To illustrate this, we consider continuous time Markov processes on epidemiological and social networks. These quantum algorithms provide an exponential advantage over existing classical Monte Carlo methods.

categories: quant-ph cs.DS cs.NA math-ph math.MP math.NA


title: COVID-19 incidences and its association with environmental quality: A country-level assessment in India

abstract: This study explored the association between the five key air pollutants (Nitrogen Dioxide (NO2), Sulphur Dioxide (SO2), Particulate Matter (PM2.5, PM10), and Carbon Monoxide (CO)) and COVID-19 incidences in India. The COVID-19 confirmed cases, air pollution concentration and meteorological variables (temperature, wind speed, surface pressure) for district and city scale were obtained for 2019 and 2020. The location-based air pollution observations were converted to a raster surface using interpolation. The deaths and positive cases are reported so far were found highest in Mumbai (436 and 11394), followed by Ahmedabad (321 and 4991), Pune (129 and 2129), Kolkata (99 and 783), Indore (83 and 1699), Jaipur (53 and 1111), Ujjain (42 and 201), Surat (37 and 799), Vadodara (31 and 400), Chennai (23 and 2647), Bhopal (22 and 652), Thane (21 and 1889), respectively. Unlike the other studies, this study has not found any substantial association between air pollution and COVID-19 incidences at the district level. Considering the number of confirmed cases, the coefficient of determination (R2) values estimated as 0.003 for PM2.5, 0.002 for PM10 and SO2, 0.001 for CO, and 0.0002 for NO2, respectively. This suggests an absolute no significant association between air pollution and COVID-19 incidences (both confirmed cases and death) in India. The same association was observed for the number of deaths as well. For COVID-19 confirmed cases, none of the five pollutants has exhibited any statistically significant association. Additionally, except the wind speed, the climate variables have no produced any statistically significant association with the COVID-19 incidences.


title: Contextual BERT: Conditioning the Language Model Using a Global State

abstract: BERT is a popular language model whose main pre-training task is to fill in the blank, i.e., predicting a word that was masked out of a sentence, based on the remaining words. In some applications, however, having an additional context can help the model make the right prediction, e.g., by taking the domain or the time of writing into account. This motivates us to advance the BERT architecture by adding a global state for conditioning on a fixed-sized context. We present our two novel approaches and apply them to an industry use-case, where we complete fashion outfits with missing articles, conditioned on a specific customer. An experimental comparison to other methods from the literature shows that our methods improve personalization significantly.


title: Discrete-time portfolio optimization under maximum drawdown constraint with partial information and deep learning resolution

abstract: We study a discrete-time portfolio selection problem with partial information and maxi\-mum drawdown constraint. Drift uncertainty in the multidimensional framework is modeled by a prior probability distribution. In this Bayesian framework, we derive the dynamic programming equation using an appropriate change of measure, and obtain semi-explicit results in the Gaussian case. The latter case, with a CRRA utility function is completely solved numerically using recent deep learning techniques for stochastic optimal control problems. We emphasize the informative value of the learning strategy versus the non-learning one by providing empirical performance and sensitivity analysis with respect to the uncertainty of the drift. Furthermore, we show numerical evidence of the close relationship between the non-learning strategy and a no short-sale constrained Merton problem, by illustrating the convergence of the former towards the latter as the maximum drawdown constraint vanishes.


title: Storing short single-photon-level optical pulses in Bose-Einstein condensates for high-performance quantum memory

abstract: Large-scale quantum networks require quantum memories featuring long-lived storage of non-classical light together with efficient, high-speed and reliable operation. The concurrent realization of these features is challenging due to inherent limitations of matter platforms and light-matter interaction protocols. Here, we propose an approach to overcome this obstacle, based on the implementation of the Autler-Townes-splitting (ATS) quantum-memory protocol on a Bose-Einstein condensate (BEC) platform. We demonstrate a proof-of-principle of this approach by storing short pulses of single-photon-level light as a collective spin-excitation in a rubidium BEC. For 20 ns long-pulses, we achieve an ultra-low-noise memory with an efficiency of 30% and lifetime of 15 $\mu$s. The non-adiabatic character of the ATS protocol (leading to high-speed and low-noise operation) in combination with the intrinsically large atomic densities and ultra-low temperatures of the BEC platform (offering highly efficient and long-lived storage) opens up a new avenue towards high-performance quantum memories.


title: Glueballs and Strings in $Sp(2N)$ Yang-Mills theories

abstract: Motivated in part by the pseudo-Nambu Goldstone Boson mechanism of electroweak symmetry breaking in Composite Higgs Models, in part by dark matter scenarios with strongly coupled origin, as well as by general theoretical considerations related to the large-N extrapolation, we perform lattice studies of the Yang-Mills theories with $Sp(2N)$ gauge groups. We measure the string tension and the mass spectrum of glueballs, extracted from appropriate 2-point correlation functions of operators organised as irreducible representations of the octahedral symmetry group. We perform the continuum extrapolation and study the magnitude of finite-size effects, showing that they are negligible in our calculation. We present new numerical results for $N=1$, $2$, $3$, $4$, combine them with data previously obtained for $N=2$, and extrapolate towards $N\rightarrow \infty$. We confirm explicitly the expectation that, as already known for $N=1,2$ also for $N=3,4$ a confining potential rising linearly with the distance binds a static quark to its antiquark. We compare our results to the existing literature on other gauge groups, with particular attention devoted to the large-$N$ limit. We find agreement with the known values of the mass of the $0^{++}$, $0^{++*}$ and $2^{++}$ glueballs obtained taking the large-$N$ limit in the $SU(N)$ groups. In addition, we determine for the first time the mass of some heavier glueball states at finite $N$ in $Sp(2N)$ and extrapolate the results towards $N \rightarrow +\infty$ taking the limit in the latter groups. Since the large-$N$ limit of $Sp(2N)$ is the same as in $SU(N)$, our results are relevant also for the study of QCD-like theories.


title: Effective field theory for hydrodynamics without boosts

abstract: We formulate the Schwinger-Keldysh effective field theory of hydrodynamics without boost symmetry. This includes a spacetime covariant formulation of classical hydrodynamics without boosts with an additional conserved particle/charge current coupled to Aristotelian background sources. We find that, up to first order in derivatives, the theory is characterised by the thermodynamic equation of state and a total of 29 independent transport coefficients, in particular, 3 hydrostatic, 9 non-hydrostatic non-dissipative, and 17 dissipative. Furthermore, we study the spectrum of linearised fluctuations around anisotropic equilibrium states with non-vanishing fluid velocity. This analysis reveals a pair of sound modes that propagate at different speeds along and opposite to the fluid flow, one charge diffusion mode, and two distinct shear modes along and perpendicular to the fluid velocity. We present these results in a new hydrodynamic frame that is linearly stable irrespective of the boost symmetry in place. This provides a unified covariant stable approach for simultaneously treating Lorentzian, Galilean, and Lifshitz fluids within an effective field theory framework and sets the stage for future studies of non-relativistic intertwined patterns of symmetry breaking.


title: On construction of finite averaging sets for $SL(2, \mathbb{C})$ via its Cartan decomposition

abstract: Averaging physical quantities over Lie groups appears in many contexts across the rapidly developing branches of physics like quantum information science or quantum optics. Such an averaging process can be always represented as averaging with respect to a finite number of elements of the group, called a finite averaging set. In the previous research such sets, known as $t$-designs, were constructed only for the case of averaging over unitary groups (hence the name unitary $t$-designs). In this work we investigate the problem of constructing finite averaging sets for averaging over general non-compact matrix Lie groups, which is much more subtle task due to the fact that the the uniform invariant measure on the group manifold (the Haar measure) is infinite. We provide a general construction of such sets based on the Cartan decomposition of the group, which splits the group into its compact and non-compact components. The averaging over the compact part can be done in a uniform way, whereas the averaging over the non-compact one has to be endowed with a suppresing weight function, and can be approached using generalised Gauss quadratures. This leads us to the general form of finite averaging sets for semisimple matrix Lie groups in the product form of finite averaging sets with respect to the compact and non-compact parts. We provide an explicit calculation of such sets for the group $SL(2, \mathbb{C})$, although our construction can be applied to other cases. Possible applications of our results cover finding finite ensambles of random operations in quantum information science and quantum optics, which can be used in constructions of randomised quantum algorithms, including optical interferometric implementations.


title: Stable and efficient Petrov-Galerkin methods for a kinetic Fokker-Planck equation

abstract: We propose a stable Petrov-Galerkin discretization of a kinetic Fokker-Planck equation constructed in such a way that uniform inf-sup stability can be inferred directly from the variational formulation. Inspired by well-posedness results for parabolic equations, we derive a lower bound for the dual inf-sup constant of the Fokker-Planck bilinear form by means of stable pairs of trial and test functions. The trial function of such a pair is constructed by applying the kinetic transport operator and the inverse velocity Laplace-Beltrami operator to a given test function. For the Petrov-Galerkin projection we choose an arbitrary discrete test space and then define the discrete trial space using the same application of transport and inverse Laplace-Beltrami operator. As a result, the spaces replicate the stable pairs of the continuous level and we obtain a well-posed numerical method with a discrete inf-sup constant identical to the inf-sup constant of the continuous problem independently of the mesh size. We show how the specific basis functions can be efficiently computed by low-dimensional elliptic problems, and confirm the practicability and performance of the method with numerical experiments.


title: Quickest Bayesian and non-Bayesian detection of false data injection attack in remote state estimation

abstract: In this paper, quickest detection of false data injection attack on remote state estimation is considered. A set of $N$ sensors make noisy linear observations of a discrete-time linear process with Gaussian noise, and report the observations to a remote estimator. The challenge is the presence of a few potentially malicious sensors which can start strategically manipulating their observations at a random time in order to skew the estimates. The quickest attack detection problem for a known {\em linear} attack scheme in the Bayesian setting with a Geometric prior on the attack initiation instant is posed as a constrained Markov decision process (MDP), in order to minimize the expected detection delay subject to a false alarm constraint, with the state involving the probability belief at the estimator that the system is under attack. State transition probabilities are derived in terms of system parameters, and the structure of the optimal policy is derived analytically. It turns out that the optimal policy amounts to checking whether the probability belief exceeds a threshold. Next, generalized CUSUM based attack detection algorithm is proposed for the non-Bayesian setting where the attacker chooses the attack initiation instant in a particularly adversarial manner. It turns out that computing the statistic for the generalised CUSUM test in this setting relies on the same techniques developed to compute the state transition probabilities of the MDP. Numerical results demonstrate significant performance gain under the proposed algorithms against competing algorithms.


title: Light-Weight DDoS Mitigation at Network Edge with Limited Resources

abstract: The Internet of Things (IoT) has been growing rapidly in recent years. With the appearance of 5G, it is expected to become even more indispensable to people's lives. In accordance with the increase of Distributed Denial-of-Service (DDoS) attacks from IoT devices, DDoS defense has become a hot research topic. DDoS detection mechanisms executed on routers and SDN environments have been intensely studied. However, these methods have the disadvantage of requiring the cost and performance of the devices. In addition, there is no existing DDoS mitigation algorithm on the network edge that can be performed with the low-cost and low performance equipments. Therefore, this paper proposes a light-weight DDoS mitigation scheme at the network edge using limited resources of inexpensive devices such as home gateways. The goal of the proposed scheme is to simply detect and mitigate flooding attacks. It utilizes unused queue resources to detect malicious flows by random shuffling of queue allocation and discard the packets of the detected flows. The performance of the proposed scheme was confirmed via theoretical analysis and computer simulation. The simulation results match the theoretical results and the proposed algorithm can efficiently detect malicious flows using limited resources.


title: Quantum scattering process and information transfer out of a black hole

abstract: We calculate the probability amplitude for tree-level elastic electron-muon scattering in Minkowski spacetime with carefully prepared initial and final wave packets. The obtained nonzero amplitude implies a nonvanishing probability for detecting a recoil electron outside the light-cone of the initial muon. Transposing this Minkowski-spacetime scattering result to a near-horizon spacetime region of a massive Schwarzschild black hole and referring to a previously proposed Gedankenexperiment, we conclude that, in principle, it is possible to have information transfer from inside the black-hole horizon to outside.


title: Generic existence of multiplicity-1 minmax minimal hypersurfaces via Allen--Cahn

abstract: In Guaraco's 2018 work a new proof was given of the existence of a closed minimal hypersurface in a compact Riemannian manifold $N^{n+1}$ with $n\geq 2$. This was achieved by employing an Allen--Cahn approximation scheme and a one-parameter minmax for the Allen--Cahn energy (relying on works by Hutchinson, Tonegawa, Wickramasekera to pass to the limit as the Allen-Cahn parameter tends to $0$). The minimal hypersurface obtained may a priori carry a locally constant integer multiplicity. Here we consider a minmax construction that is a modification of the one in the aforementioned work, by allowing an initial freedom on the choice of the valley points between which the mountain pass construction is carried out, and then optimising over said choice. We prove that, when $2\leq n\leq 6$ and the metric is bumpy, this minmax leads to a (smooth closed) minimal hypersurface with multiplicity $1$. (When $n=2$ this conclusion also follows from Chodosh--Mantoulidis's recent work.) As immediate corollary we obtain that every compact Riemannian manifold of dimension $n+1$, $2\leq n\leq 6$, endowed with a bumpy metric, admits a two-sided smooth closed minimal hypersurface (this existence conclusion also follows from Zhou's recent result for minmax constructions via Almgren--Pitts theory).


title: Virtual Morse-Bott index, moduli spaces of pairs, and applications to topology of smooth four-manifolds

abstract: We previously developed an approach to Bialynicki-Birula theory for holomorphic $\mathbb{C}^*$ actions on complex analytic spaces and the concept of virtual Morse-Bott indices for singular critical points of Hamiltonian functions for the induced circle actions (see Feehan, arXiv:2206.14710). For Hamiltonian functions of circle actions on closed, complex Kaehler manifolds, the virtual Morse-Bott index coincides with the classical Morse-Bott index due to Bott (1954) and Frankel (1959). A key principle in our approach is that positivity of the virtual Morse-Bott index at a critical point of the Hamiltonian function implies that the critical point cannot be a local minimum even when that critical point is a singular point in the moduli space. In this monograph, we consider our method in the context of the moduli space of non-Abelian monopoles over a closed, complex, Kaehler surface. We use the Hirzebruch-Riemann-Roch Theorem to compute virtual Morse-Bott indices of all critical strata (Seiberg-Witten moduli subspaces) and we prove that these indices are positive in a setting motivated by the conjecture that all closed, smooth four-manifolds of Seiberg-Witten simple type obey the Bogomolov-Miyaoka-Yau inequality.


title: Non-Hermitian Yukawa interactions of fermions with axions: potential microscopic origin and dynamical mass generation

abstract: In this mini review, we discuss some recent developments regarding properties of (quantum) field-theory models containing anti-Hermitian Yukawa interactions between pseudoscalar fields (axions) and Dirac (or Majorana) fermions. Specifically, after motivating physically such interactions, in the context of string-inspired low-energy effective field theories, involving right-handed neutrinos and axion fields, we proceed to discuss their formal consistency within the so-called Parity-Time-reversal(PT)-symmetry framework, as well as dynamical mass generation, induced by the Yukawa interactions, for both fermions and axions. The Yukawa couplings are assumed weak, given that they are conjectured to have been generated by non-perturbative effects in the underlying microscopic string theory. The models under discussion contain, in addition to the Yukawa terms, also anti-Hermitian anomalous derivative couplings of the pseudoscalar fields to axial fermion currents, as well as interactions of the fermions with non-Hermitian axial backgrounds. We discuss the role of such additional couplings on the Yukawa-induced dynamically-generated masses. For the case where the fermions are right-handed neutrinos, we compare such masses with the radiative ones induced by both, the anti-Hermitian anomalous terms and the anti-Hermitian Yukawa interactions in phenomenologically relevant models.


title: Demonstration of electron cooling using a pulsed beam from an electrostatic electron cooler

abstract: Cooling of hadron beams is critically important in the next generation of hadron storage rings for delivery of unprecedented performance. One such application is the electron-ion collider presently under development in the US. The desire to develop electron coolers for operation at much higher energies than previously achieved necessitates the use of radio-frequency (RF) fields for acceleration as opposed to the conventional, electrostatic approach. While electron cooling is a mature technology at low energy utilizing a DC beam, RF acceleration requires the cooling beam to be bunched, thus extending the parameter space to an unexplored territory. It is important to experimentally demonstrate the feasibility of cooling with electron bunches and further investigate how the relative time structure of the two beams affects the cooling properties; thus, a set of four pulsed-beam cooling experiments was carried out by a collaboration of Jefferson Lab and Institute of Modern Physics (IMP). The experiments have successfully demonstrated cooling with a beam of electron bunches in both the longitudinal and transverse directions for the first time. We have measured the effect of the electron bunch length and longitudinal ion focusing strength on the temporal evolution of the longitudinal and transverse ion beam profile and demonstrate that if the synchronization can be accurately maintained, the dynamics are not adversely affected by the change in time structure.


title: A Framework for Learning Predator-prey Agents from Simulation to Real World

abstract: In this paper, we propose an evolutionary predatorprey robot system which can be generally implemented from simulation to the real world. We design the closed-loop robot system with camera and infrared sensors as inputs of controller. Both the predators and prey are co-evolved by NeuroEvolution of Augmenting Topologies (NEAT) to learn the expected behaviours. We design a framework that integrate Gym of OpenAI, Robot Operating System (ROS), Gazebo. In such a framework, users only need to focus on algorithms without being worried about the detail of manipulating robots in both simulation and the real world. Combining simulations, real-world evolution, and robustness analysis, it can be applied to develop the solutions for the predator-prey tasks. For the convenience of users, the source code and videos of the simulated and real world are published on Github.


title: A computational periporomechanics model for localized failure in unsaturated porous media

abstract: We implement a computational periporomechanics model for simulating localized failure in unsaturated porous media. The coupled periporomechanics model is based on the peridynamic state concept and the effective force state concept. The coupled governing equations are integral-differential equations without assuming the continuity of solid displacement and fluid pressures. The fluid flow and effective force states are determined by nonlocal fluid pressure and deformation gradients through the recently formulated multiphase constitutive correspondence principle. The coupled peri-poromechanics is implemented numerically for high-performance computing by an implicit multiphase meshfree method utilizing the message passing interface. The numerical implementation is validated by simulating classical poromechanics problems and comparing the numerical results with analytical solutions and experimental data. Numerical examples are presented to demonstrate the robustness of the fully coupled peri-poromechanics in modeling localized failures in unsaturated porous media.


title: Eccentricity queries and beyond using Hub Labels

abstract: Hub labeling schemes are popular methods for computing distances on road networks and other large complex networks, often answering to a query within a few microseconds for graphs with millions of edges. In this work, we study their algorithmic applications beyond distance queries. We focus on eccentricity queries and distance-sum queries, for several versions of these problems on directed weighted graphs, that is in part motivated by their importance in facility location problems. On the negative side, we show conditional lower bounds for these above problems on unweighted undirected sparse graphs, via standard constructions from "Fine-grained" complexity. However, things take a different turn when the hub labels have a sublogarithmic size. Indeed, given a hub labeling of maximum label size $\leq k$, after pre-processing the labels in total $2^{{O}(k)} \cdot |V|^{1+o(1)}$ time, we can compute both the eccentricity and the distance-sum of any vertex in $2^{{O}(k)} \cdot |V|^{o(1)}$ time. It can also be applied to the fast global computation of some topological indices. Finally, as a by-product of our approach, on any fixed class of unweighted graphs with bounded expansion, we can decide whether the diameter of an $n$-vertex graph in the class is at most $k$ in $f(k) \cdot n^{1+o(1)}$ time, for some "explicit" function $f$.


title: Tidal Love Numbers of Kerr Black Holes

abstract: The open question of whether a Kerr black hole can become tidally deformed or not has profound implications for fundamental physics and gravitational-wave astronomy. We consider a Kerr black hole embedded in a weak and slowly varying, but otherwise arbitrary, multipolar tidal environment. By solving the static Teukolsky equation for the gauge-invariant Weyl scalar $\psi_0$, and by reconstructing the corresponding metric perturbation in an ingoing radiation gauge, for a general harmonic index $\ell$, we compute the linear response of a Kerr black hole to the tidal field. This linear response vanishes identically for a Schwarzschild black hole and for an axisymmetric perturbation of a spinning black hole. For a nonaxisymmetric perturbation of a spinning black hole, however, the linear response does not vanish, and it contributes to the Geroch-Hansen multipole moments of the perturbed Kerr geometry. As an application, we compute explicitly the rotational black hole tidal Love numbers that couple the induced quadrupole moments to the quadrupolar tidal fields, to linear order in the black hole spin, and we introduce the corresponding notion of tidal Love tensor. Finally, we show that those induced quadrupole moments are closely related to the well-known physical phenomenon of tidal torquing of a spinning body interacting with a tidal gravitational environment.


title: Momentum Entanglement for Atom Interferometry

abstract: Compared to light interferometers, the flux in cold-atom interferometers is low and the associated shot noise large. Sensitivities beyond these limitations require the preparation of entangled atoms in different momentum modes. Here, we demonstrate a source of entangled atoms that is compatible with state-of-the-art interferometers. Entanglement is transferred from the spin degree of freedom of a Bose-Einstein condensate to well-separated momentum modes, witnessed by a squeezing parameter of -3.1(8) dB. Entanglement-enhanced atom interferometers open up unprecedented sensitivities for quantum gradiometers or gravitational wave detectors.


title: Broadening and sharpening of the Drude peak through antiferromagnetic fluctuations

abstract: Antiferromagnetic or charge density wave fluctuations couple with light through the recently discovered {\pi}-ton contribution to the optical conductivity, and quite generically constitute the dominant vertex corrections in low-dimensional correlated electron systems. Here we study the arguably simplest version of these $\pi$-tons based on the semi-analytical random phase approximation (RPA) ladder in the transversal particle-hole channel. The vertex corrections to the optical conductivity are calculated directly for real frequencies. We validate that the RPA qualitatively reproduces the {\pi}-ton vertex corrections to the Drude peak in the Hubbard model. Depending on the temperature we find vertex corrections to broaden or sharpen the Drude peak.


title: The impact of (n,$\gamma$) reaction rate uncertainties of unstable isotopes on the i-process nucleosynthesis of the elements from Ba to W

abstract: The abundances of n-capture elements in the CEMP-r/s stars agree with predictions of intermediate n-density nucleosynthesis, at $N_\mathrm{n}\sim 10^{13}$-$10^{15} \mathrm{cm}^{-3}$, in rapidly-accreting white dwarfs (RAWDs). We have performed Monte-Carlo simulations of this i-process nucleosynthesis to determine the impact of (n,$\gamma$) reaction rate uncertainties of 164 unstable isotopes, from $^{131}$I to $^{189}$Hf, on the predicted abundances of 18 elements from Ba to W. The impact study is based on two representative one-zone models with constant values of $N_\mathrm{n} = 3.16\times 10^{14}\ \mathrm{cm}^{-3}$ and $N_\mathrm{n} = 3.16\times 10^{13}\ \mathrm{cm}^{-3}$ and on a multi-zone model based on a realistic stellar evolution simulation of He-shell convection entraining H in a RAWD model with [Fe/H]=-2.6. For each of the selected elements, we have identified up to two (n,$\gamma$) reactions having the strongest correlations between their rate variations constrained by Hauser-Feshbach computations and the predicted abundances, with the Pearson product-moment correlation coefficients $|r_\mathrm{P}| > 0.15$. We find that the discrepancies between the predicted and observed abundances of Ba and Pr in the CEMP-i star CS31062-050 are significantly diminished if the rate of $^{137}$Cs(n,$\gamma)^{138}$Cs is reduced and the rates of $^{141}$Ba(n,$\gamma)^{142}$Ba or $^{141}$La(n,$\gamma)^{142}$La increased. The uncertainties of temperature-dependent $\beta$-decay rates of the same unstable isotopes have a negligible effect on the predicted abundances. One-zone Monte-Carlo simulations can be used instead of computationally time-consuming multi-zone Monte-Carlo simulations in reaction rate uncertainty studies if they use comparable values of $N_\mathrm{n}$ (abridged).


title: A geographical study of $\overline{\mathcal M}_{2}(\mathbb{P}^2,4)^{\text{main}}$

abstract: We discuss criteria for a stable map of genus two and degree $4$ to the projective plane to be smoothable, as an application of our modular desingularisation of $\overline{\mathcal M}_{2,n}(\mathbb{P}^r,d)^{\text{main}}$ via logarithmic geometry and Gorenstein singularities.


title: Magnetoresistance oscillation study of the half-quantum vortex in doubly connected mesoscopic superconducting cylinders of Sr2RuO4

abstract: The observation of the highly unusual half-quantum vortex (HQV) in a single crystalline superconductor excludes unequivocally the spin-singlet symmetry of the superconducting order parameter. HQVs were observed previously in mesoscopic samples of Sr2RuO4 in cantilever torque magnetometry measurements, thus providing direct evidence for spin-triplet pairing in the material. In addition, it raised important questions on HQV, including its stability and dynamics. These issues have remained largely unexplored, in particular, experimentally. We report in this paper the detection of HQVs in mesoscopic, doubly connected cylinders of single-crystalline Sr2RuO4 of a mesoscopic size and the examination of the effect of the in-plane magnetic field needed for the observation of the HQV by magnetoresistance (MR) oscillations measurements. Several distinct features found in our data, especially a dip and secondary peaks in the MR oscillations seen only in the presence of a sufficiently large in-plane magnetic field as well as a large measurement current, are linked to the formation of the HQV fluxoid state in and crossing of an Abrikosov HQV through the sample. The conclusion is drawn from the analysis of our data using a model of thermally activated vortex crossing overcoming a free-energy barrier which is modulated by the applied magnetic flux enclosed in the cylinder as well as the measurement current. Evidence for the trapping of an HQV fluxoid state in the sample was also found. Our observation of the HQV in mesoscopic Sr2RuO4 provided not only additional evidence for spin-triplet superconductivity in Sr2RuO4 but also insights into the physics of HQV, including its spontaneous spin polarization, stability, and dynamics. Our study also revealed a possible effect of the measurement current on the magnitude of the spontaneous spin polarization associated with the HQV.


title: Ray-marching Thurston geometries

abstract: We describe algorithms that produce accurate real-time interactive in-space views of the eight Thurston geometries using ray-marching. We give a theoretical framework for our algorithms, independent of the geometry involved. In addition to scenes within a geometry $X$, we also consider scenes within quotient manifolds and orbifolds $X / \Gamma$. We adapt the Phong lighting model to non-euclidean geometries. The most difficult part of this is the calculation of light intensity, which relates to the area density of geodesic spheres. We also give extensive practical details for each geometry.


title: A solution to Erd\H{o}s and Hajnal's odd cycle problem

abstract: In 1981, Erd\H{o}s and Hajnal asked whether the sum of the reciprocals of the odd cycle lengths in a graph with infinite chromatic number is necessarily infinite. Let $\mathcal{C}(G)$ be the set of cycle lengths in a graph $G$ and let $\mathcal{C}_\text{odd}(G)$ be the set of odd numbers in $\mathcal{C}(G)$. We prove that, if $G$ has chromatic number $k$, then $\sum_{\ell\in \mathcal{C}_\text{odd}(G)}1/\ell\geq (1/2-o_k(1))\log k$. This solves Erd\H{o}s and Hajnal's odd cycle problem, and, furthermore, this bound is asymptotically optimal. In 1984, Erd\H{o}s asked whether there is some $d$ such that each graph with chromatic number at least $d$ (or perhaps even only average degree at least $d$) has a cycle whose length is a power of 2. We show that an average degree condition is sufficient for this problem, solving it with methods that apply to a wide range of sequences in addition to the powers of 2. Finally, we use our methods to show that, for every $k$, there is some $d$ so that every graph with average degree at least $d$ has a subdivision of the complete graph $K_k$ in which each edge is subdivided the same number of times. This confirms a conjecture of Thomassen from 1984.


title: Isometric embeddings in trees and their use in the diameter problem

abstract: We prove that given a discrete space with $n$ points which is either embedded in a system of $k$ trees, or the Cartesian product of $k$ trees, we can compute all eccentricities in ${\cal O}(2^{{\cal O}(k\log{k})}(N+n)^{1+o(1)})$ time, where $N$ is the cumulative total order over all these $k$ trees. This is near optimal under the Strong Exponential-Time Hypothesis, even in the very special case of an $n$-vertex graph embedded in a system of $\omega(\log{n})$ spanning trees. However, given such an embedding in the strong product of $k$ trees, there is a much faster ${\cal O}(N + kn)$-time algorithm for this problem. All our positive results can be turned into approximation algorithms for the graphs and finite spaces with a quasi isometric embedding in trees, if such embedding is given as input, where the approximation factor (resp., the approximation constant) depends on the distortion of the embedding (resp., of its stretch). The existence of embeddings in the Cartesian product of finitely many trees has been thoroughly investigated for cube-free median graphs. We give the first-known quasi linear-time algorithm for computing the diameter within this graph class. It does not require an embedding in a product of trees to be given as part of the input. On our way, being given an $n$-node tree $T$, we propose a data structure with ${\cal O}(n\log{n})$ pre-processing time in order to compute in ${\cal O}(k\log^2{n})$ time the eccentricity of any subset of $k$ nodes. We combine the latter technical contribution, of independent interest, with a recent distance-labeling scheme that was designed for cube-free median graphs.


title: Getzler-Kapranov complexes and moduli stacks of curves

abstract: In this paper, we study the so-called Getzler-Kapranov complexes and their relation to the cohomology of moduli stacks of curves.


title: A Local Search Framework for Experimental Design

abstract: We present a local search framework to design and analyze both combinatorial algorithms and rounding algorithms for experimental design problems. This framework provides a unifying approach to match and improve all known results in D/A/E-design and to obtain new results in previously unknown settings. For combinatorial algorithms, we provide a new analysis of the classical Fedorov's exchange method. We prove that this simple local search algorithm works well as long as there exists an almost optimal solution with good condition number. Moreover, we design a new combinatorial local search algorithm for E-design using the regret minimization framework. For rounding algorithms, we provide a unified randomized exchange algorithm to match and improve previous results for D/A/E-design. Furthermore, the algorithm works in the more general setting to approximately satisfy multiple knapsack constraints, which can be used for weighted experimental design and for incorporating fairness constraints into experimental design.


title: Renormalization Group Improvement of the Effective Potential: an EFT Approach

abstract: We apply effective field theory (EFT) methods to compute the renormalization group improved effective potential for theories with a large mass hierarchy. Our method allows one to compute the effective potential in a systematic expansion in powers of the mass ratio, as well as to sum large logarithms of mass ratios using renormalization group evolution. The effective potential is the sum of one-particle irreducible diagrams (1PI) but information about which diagrams are 1PI is lost after matching to the EFT, since heavy lines get shrunk to a point. We therefore introduce a tadpole condition in place of the 1PI condition, and use the renormalization group improved value of the tadpole in computing the effective potential. We explain why the effective potential computed using an EFT is not the same as the effective potential of the EFT. We illustrate our method using the $O(N)$ model, a theory of two scalars in the unbroken and broken phases, and the Higgs-Yukawa model. Our leading-log result, obtained by integrating the one-loop $\beta$-functions, correctly reproduces the log-squared term in explicit two-loop calculations. Our method does not have a Goldstone boson infrared divergence problem.


title: A data processing system for balloon-borne telescopes

abstract: The JEM-EUSO Collaboration aims at studying Ultra High Energy Cosmic Rays (UHECR) from space. To reach this goal, a series of pathfinder missions has been developed to prove the observation principle and to raise the technological readiness level of the instrument. Among these, the EUSO-SPB2 (Extreme Universe Space Observatory on a Super Pressure Balloon, mission two) foresees the launch of two telescopes on an ultra-long duration balloon. One is a fluorescence telescope designed to detect UHECR via the UV fluorescence emission of the showers in the atmosphere. The other one measures direct Cherenkov light emission from lower energy cosmic rays and other optical backgrounds for cosmogenic tau neutrino detection. In this paper, we describe the data processing system which has been designed to perform data management and instrument control for the two telescopes. It is a complex which controls front-end electronics, tags events with arrival time and payload position through a GPS system, provides signals for time synchronization of the event and measures live and dead time of the telescope. In addition, the data processing system manages mass memory for data storage, performs housekeeping monitor, and controls power on and power off sequences. The target flight duration for the NASA super pressure program is 100 days, consequently, the requirements on the electronics and the data handling are quite severe. The system operates at high altitude in unpressurised environment, which introduces a technological challenge for heat dissipation.


title: Learning Bayesian Networks from Ordinal Data

abstract: Bayesian networks are a powerful framework for studying the dependency structure of variables in a complex system. The problem of learning Bayesian networks is tightly associated with the given data type. Ordinal data, such as stages of cancer, rating scale survey questions, and letter grades for exams, are ubiquitous in applied research. However, existing solutions are mainly for continuous and nominal data. In this work, we propose an iterative score-and-search method - called the Ordinal Structural EM (OSEM) algorithm - for learning Bayesian networks from ordinal data. Unlike traditional approaches designed for nominal data, we explicitly respect the ordering amongst the categories. More precisely, we assume that the ordinal variables originate from marginally discretizing a set of Gaussian variables, whose structural dependence in the latent space follows a directed acyclic graph. Then, we adopt the Structural EM algorithm and derive closed-form scoring functions for efficient graph searching. Through simulation studies, we illustrate the superior performance of the OSEM algorithm compared to the alternatives and analyze various factors that may influence the learning accuracy. Finally, we demonstrate the practicality of our method with a real-world application on psychological survey data from 408 patients with co-morbid symptoms of obsessive-compulsive disorder and depression.


title: The ins and outs of speaker recognition: lessons from VoxSRC 2020

abstract: The VoxCeleb Speaker Recognition Challenge (VoxSRC) at Interspeech 2020 offers a challenging evaluation for speaker recognition systems, which includes celebrities playing different parts in movies. The goal of this work is robust speaker recognition of utterances recorded in these challenging environments. We utilise variants of the popular ResNet architecture for speaker recognition and perform extensive experiments using a range of loss functions and training parameters. To this end, we optimise an efficient training framework that allows powerful models to be trained with limited time and resources. Our trained models demonstrate improvements over most existing works with lighter models and a simple pipeline. The paper shares the lessons learned from our participation in the challenge.


title: Naive analytics equilibrium

abstract: We study interactions with uncertainty about demand sensitivity. In our solution concept (1) firms choose seemingly-optimal strategies given the level of sophistication of their data analytics, and (2) the levels of sophistication form best responses to one another. Under the ensuing equilibrium firms underestimate price elasticities and overestimate advertising effectiveness, as observed empirically. The misestimates cause firms to set prices too high and to over-advertise. In games with strategic complements (substitutes), profits Pareto dominate (are dominated by) those of the Nash equilibrium. Applying the model to team production games explains the prevalence of overconfidence among entrepreneurs and salespeople.


title: Algorithmic pure states for the negative spherical perceptron

abstract: We consider the spherical perceptron with Gaussian disorder. This is the set $S$ of points $\sigma \in \mathbb{R}^N$ on the sphere of radius $\sqrt{N}$ satisfying $\langle g_a , \sigma \rangle \ge \kappa\sqrt{N}\,$ for all $1 \le a \le M$, where $(g_a)_{a=1}^M$ are independent standard gaussian vectors and $\kappa \in \mathbb{R}$ is fixed. Various characteristics of $S$ such as its surface measure and the largest $M$ for which it is non-empty, were computed heuristically in statistical physics in the asymptotic regime $N \to \infty$, $M/N \to \alpha$. The case $\kappa<0$ is of special interest as $S$ is conjectured to exhibit a hierarchical tree-like geometry known as "full replica-symmetry breaking" (FRSB) close to the satisfiability threshold $\alpha_{\text{SAT}}(\kappa)$, and whose characteristics are captured by a Parisi variational principle akin to the one appearing in the Sherrington-Kirkpatrick model. In this paper we design an efficient algorithm which, given oracle access to the solution of the Parisi variational principle, exploits this conjectured FRSB structure for $\kappa<0$ and outputs a vector $\hat{\sigma}$ satisfying $\langle g_a , \hat{\sigma}\rangle \ge \kappa \sqrt{N}$ for all $1\le a \le M$ and lying on a sphere of non-trivial radius $\sqrt{\bar{q} N}$, where $\bar{q} \in (0,1)$ is the right-end of the support of the associated Parisi measure. We expect $\hat{\sigma}$ to be approximately the barycenter of a pure state of the spherical perceptron. Moreover we expect that $\bar{q} \to 1$ as $\alpha \to \alpha_{\text{SAT}}(\kappa)$, so that $\big\langle g_a,\hat{\sigma}/|\hat{\sigma}|\big\rangle \geq (\kappa-o(1))\sqrt{N}$ near criticality.


title: Semi-polarized meromorphic Hitchin and Calabi-Yau integrable systems

abstract: It was shown by Diaconescu, Donagi and Pantev that Hitchin systems of type ADE are isomorphic to certain Calabi-Yau integrable systems. In this paper, we prove an analogous result in the setting of meromorphic Hitchin systems of type A which are known to be Poisson integrable systems. We consider a symplectization of the meromorphic Hitchin integrable system, which is a semi-polarized integrable system in the sense of Kontsevich and Soibelman. On the Hitchin side, we show that the moduli space of unordered diagonally framed Higgs bundles forms an integrable system in this sense and recovers the meromorphic Hitchin system as the fiberwise compact quotient. Then we construct a new family of quasi-projective Calabi-Yau threefolds and show that its relative intermediate Jacobian fibration, as semi-polarized integrable systems, is isomorphic to the moduli space of unordered diagonally framed Higgs bundles.


title: Thermal transport in $T\bar{T}$-deformed conformal field theories: from integrability to holography

abstract: In this paper we consider the energy and momentum transport in (1+1)-dimension conformal field theories (CFTs) that are deformed by an irrelevant operator $T\bar{T}$, using the integrability based generalized hydrodynamics, and holography. The two complementary methods allow us to study the energy and momentum transport after the in-homogeneous quench, derive the exact non-equilibrium steady states (NESS) and calculate the Drude weights and the diffusion constants. Our analysis reveals that all of these quantities satisfy universal formulae regardless of the underlying CFT, thereby generalizing the universal formulae for these quantities in pure CFTs. As a sanity check, we also confirm that the exact momentum diffusion constant agrees with the conformal perturbation. These fundamental physical insights have important consequences for our understanding of the $T\bar{T}$-deformed CFTs. First of all, they provide the first check of the $T\bar{T}$-deformed $\mathrm{AdS}_3$/$\mathrm{CFT}_2$ correspondence from the dynamical standpoint. And secondly, we are able to identify a remarkable connection between the $T\bar{T}$-deformed CFTs and reversible cellular automata.


title: Around the diameter of AT-free graphs

abstract: A graph algorithm is truly subquadratic if it runs in ${\cal O}(m^b)$ time on connected $m$-edge graphs, for some positive $b < 2$. Roditty and Vassilevska Williams (STOC'13) proved that under plausible complexity assumptions, there is no truly subquadratic algorithm for computing the diameter of general graphs. In this work, we present positive and negative results on the existence of such algorithms for computing the diameter on some special graph classes. Specifically, three vertices in a graph form an asteroidal triple (AT) if between any two of them there exists a path that avoids the closed neighbourhood of the third one. We call a graph AT-free if it does not contain an AT. We first prove that for all $m$-edge AT-free graphs, one can compute all the eccentricities in truly subquadratic ${\cal O}(m^{3/2})$ time. Then, we extend our study to several subclasses of chordal graphs -- all of them generalizing interval graphs in various ways --, as an attempt to understand which of the properties of AT-free graphs, or natural generalizations of the latter, can help in the design of fast algorithms for the diameter problem on broader graph classes. For instance, for all chordal graphs with a dominating shortest path, there is a linear-time algorithm for computing a diametral pair if the diameter is at least four. However, already for split graphs with a dominating edge, under plausible complexity assumptions, there is no truly subquadratic algorithm for deciding whether the diameter is either $2$ or $3$.


title: On $k$-para-K\"ahler Lie algebras a subclass of $k$-symplectic Lie algebras

abstract: $k$-Para-K\"ahler Lie algebras are a generalization of para-K\"ahler Lie algebras $(k=1)$ and constitute a subclass of $k$-symplectic Lie algebras. In this paper, we show that the characterization of para-K\"ahler Lie algebras as left symmetric bialgebras can be generalized to $k$-para-K\"ahler Lie algebras leading to the introduction of two new structures which are different but both generalize the notion of left symmetric algebra. This permits also the introduction of generalized $S$-matrices. We determine then all the $k$-symplectic Lie algebras of dimension $(k+1)$ and all the six dimensional 2-para-K\"ahler Lie algebras.


title: Post-processing of quantum instruments

abstract: Studying sequential measurements is of the utmost importance to both the foundational aspects of quantum theory and the practical implementations of quantum technologies, with both of these applications being abstractly described by the concatenation of quantum instruments into a sequence of certain length. In general, the choice of instrument at any given step in the sequence can be conditionally chosen based on the classical results of all preceding instruments. For two instruments in a sequence we consider the conditional second instrument as an effective way of post-processing the first instrument into a new one. This is similar to how a measurement described by a positive operator-valued measure (POVM) can be post-processed into another by way of classical randomization of its outcomes using a stochastic matrix. In this work we study the post-processing relation of instruments and the partial order it induces on their equivalence classes. We characterize the greatest and the least element of this order, give examples of post-processings between different types of instruments and draw connections between post-processings of some of these instruments and their induced POVMs.


title: $\sigma$-Ridge: group regularized ridge regression via empirical Bayes noise level cross-validation

abstract: Features in predictive models are not exchangeable, yet common supervised models treat them as such. Here we study ridge regression when the analyst can partition the features into $K$ groups based on external side-information. For example, in high-throughput biology, features may represent gene expression, protein abundance or clinical data and so each feature group represents a distinct modality. The analyst's goal is to choose optimal regularization parameters $\lambda = (\lambda_1, \dotsc, \lambda_K)$ -- one for each group. In this work, we study the impact of $\lambda$ on the predictive risk of group-regularized ridge regression by deriving limiting risk formulae under a high-dimensional random effects model with $p\asymp n$ as $n \to \infty$. Furthermore, we propose a data-driven method for choosing $\lambda$ that attains the optimal asymptotic risk: The key idea is to interpret the residual noise variance $\sigma^2$, as a regularization parameter to be chosen through cross-validation. An empirical Bayes construction maps the one-dimensional parameter $\sigma$ to the $K$-dimensional vector of regularization parameters, i.e., $\sigma \mapsto \widehat{\lambda}(\sigma)$. Beyond its theoretical optimality, the proposed method is practical and runs as fast as cross-validated ridge regression without feature groups ($K=1$).


title: Characteristic formulation of the Regge-Wheeler and Zerilli Green functions

abstract: We present a characteristic initial value approach to calculating the Green function of the Regge-Wheeler and Zerilli equations. We combine well-known numerical methods with newly derived initial data to obtain a scheme which can in principle be generalised to any desired order of convergence. We demonstrate the approach with implementations up to sixth-order in the grid spacing. By combining the results of our numerical code with late-time tail expansions and methods of subtracting the direct part of the Green function, we show that the scalar self-force in Schwarzschild spacetime can be computed to better accuracy than previous Green-function based approaches. We also demonstrate agreement with frequency-domain methods for computing the Green function in the gravitational case. Finally, we apply the Regge-Wheeler and Zerilli Green functions to the computation of the gravitational energy flux.


title: Tensor Completion via Tensor Networks with a Tucker Wrapper

abstract: In recent years, low-rank tensor completion (LRTC) has received considerable attention due to its applications in image/video inpainting, hyperspectral data recovery, etc. With different notions of tensor rank (e.g., CP, Tucker, tensor train/ring, etc.), various optimization based numerical methods are proposed to LRTC. However, tensor network based methods have not been proposed yet. In this paper, we propose to solve LRTC via tensor networks with a Tucker wrapper. Here by "Tucker wrapper" we mean that the outermost factor matrices of the tensor network are all orthonormal. We formulate LRTC as a problem of solving a system of nonlinear equations, rather than a constrained optimization problem. A two-level alternative least square method is then employed to update the unknown factors. The computation of the method is dominated by tensor matrix multiplications and can be efficiently performed. Also, under proper assumptions, it is shown that with high probability, the method converges to the exact solution at a linear rate. Numerical simulations show that the proposed algorithm is comparable with state-of-the-art methods.


title: Down the bot hole: actionable insights from a 1-year analysis of bots activity on Twitter

abstract: Nowadays, social media represent persuasive tools that have been progressively weaponized to affect people's beliefs, spread manipulative narratives, and sow conflicts along divergent factions. Software-controlled accounts (i.e., bots) are one of the main actors associated with manipulation campaigns, especially in the political context. Uncovering the strategies behind bots' activities is of paramount importance to detect and curb such campaigns. In this paper, we present a long term (one year) analysis of bots activity on Twitter in the run-up to the 2018 U.S. Midterm Elections. We identify different classes of accounts based on their nature (bot vs. human) and engagement within the online discussion and we observe that hyperactive bots played a pivotal role in the dissemination of conspiratorial narratives, while dominating the political debate since the year before the election. Our analysis, on the horizon of the upcoming U.S. 2020 Presidential Election, reveals both alarming findings of humans' susceptibility to bots and actionable insights that can contribute to curbing coordinated campaigns.


title: Cream of the Crop: Distilling Prioritized Paths For One-Shot Neural Architecture Search

abstract: One-shot weight sharing methods have recently drawn great attention in neural architecture search due to high efficiency and competitive performance. However, weight sharing across models has an inherent deficiency, i.e., insufficient training of subnetworks in hypernetworks. To alleviate this problem, we present a simple yet effective architecture distillation method. The central idea is that subnetworks can learn collaboratively and teach each other throughout the training process, aiming to boost the convergence of individual models. We introduce the concept of prioritized path, which refers to the architecture candidates exhibiting superior performance during training. Distilling knowledge from the prioritized paths is able to boost the training of subnetworks. Since the prioritized paths are changed on the fly depending on their performance and complexity, the final obtained paths are the cream of the crop. We directly select the most promising one from the prioritized paths as the final architecture, without using other complex search methods, such as reinforcement learning or evolution algorithms. The experiments on ImageNet verify such path distillation method can improve the convergence ratio and performance of the hypernetwork, as well as boosting the training of subnetworks. The discovered architectures achieve superior performance compared to the recent MobileNetV3 and EfficientNet families under aligned settings. Moreover, the experiments on object detection and more challenging search space show the generality and robustness of the proposed method. Code and models are available at https://github.com/microsoft/cream.git.


title: Quintessential $\alpha$-attractor inflation: forecasts for Stage IV galaxy surveys

abstract: Single-field models of $\alpha$-attractor quintessential inflation provide a unified picture of the two periods of early- and late-time cosmic acceleration, where both inflation and dark energy are described by a single scalar degree of freedom rolling down a runaway potential. These theoretically well-motivated models have distinct observational predictions that are in agreement with existing cosmological data. We show that the next generation of large-scale structure surveys, even when no other cosmological data sets are considered, will strongly constrain the parameter space of these models, and test them against the standard cosmological model and more conventional non-quintessential inflation. In particular, we expect $\mathcal{O}(10^{-5}\mathrm{-}10^{-4})$ constraints on the present values of the dark energy equation of state and its time derivative, $w_0$ and $w_a$. We also forecast more than one order of magnitude tighter constraints on the spectral index of primordial curvature perturbations $n_s$ compared to the expectations for the standard model. This demonstrates the powerful synergy between the upcoming large-scale structure probes of inflation and those aiming to measure the tensor-to-scalar ratio $r$ through the observation of $B$-mode polarization of the cosmic microwave background.


title: Black-Box Optimization of Object Detector Scales

abstract: Object detectors have improved considerably in the last years by using advanced CNN architectures. However, many detector hyper-parameters are generally manually tuned, or they are used with values set by the detector authors. Automatic Hyper-parameter optimization has not been explored in improving CNN-based object detectors hyper-parameters. In this work, we propose the use of Black-box optimization methods to tune the prior/default box scales in Faster R-CNN and SSD, using Bayesian Optimization, SMAC, and CMA-ES. We show that by tuning the input image size and prior box anchor scale on Faster R-CNN mAP increases by 2% on PASCAL VOC 2007, and by 3% with SSD. On the COCO dataset with SSD there are mAP improvement in the medium and large objects, but mAP decreases by 1% in small objects. We also perform a regression analysis to find the significant hyper-parameters to tune.


title: Passport-aware Normalization for Deep Model Protection

abstract: Despite tremendous success in many application scenarios, deep learning faces serious intellectual property (IP) infringement threats. Considering the cost of designing and training a good model, infringements will significantly infringe the interests of the original model owner. Recently, many impressive works have emerged for deep model IP protection. However, they either are vulnerable to ambiguity attacks, or require changes in the target network structure by replacing its original normalization layers and hence cause significant performance drops. To this end, we propose a new passport-aware normalization formulation, which is generally applicable to most existing normalization layers and only needs to add another passport-aware branch for IP protection. This new branch is jointly trained with the target model but discarded in the inference stage. Therefore it causes no structure change in the target model. Only when the model IP is suspected to be stolen by someone, the private passport-aware branch is added back for ownership verification. Through extensive experiments, we verify its effectiveness in both image and 3D point recognition models. It is demonstrated to be robust not only to common attack techniques like fine-tuning and model compression, but also to ambiguity attacks. By further combining it with trigger-set based methods, both black-box and white-box verification can be achieved for enhanced security of deep learning models deployed in real systems. Code can be found at https://github.com/ZJZAC/Passport-aware-Normalization.


title: Sensitivity of the Polyakov loop to chiral symmetry restoration

abstract: In the heavy, static quark mass regime of QCD, the Polyakov loop is well known to be an order parameter of the deconfinement phase transition; however, the sensitivity of the Polyakov loop to the deconfinement of light, dynamical quarks is less clear. On the other hand, from the perspective of an effective Lagrangian written in the vicinity of the chiral transition, the Polyakov loop is an energy-like operator and should hence scale as any energy-like operator would. We show here that the Polyakov loop and heavy-quark free energy are sensitive to the chiral transition, i.e. their scaling is consistent with energy-like observables in 3-$d$ $O(N)$ universality classes.


title: AdS one-loop partition functions from bulk and edge characters

abstract: We show that the one-loop partition function of any higher spin field in $(d+1)$-dimensional Anti-de Sitter spacetime can be expressed as an integral transform of an $\text{SO}(2,d)$ bulk character and an $\text{SO}(2,d-2)$ edge character. We apply this character integral formula to various higher-spin Vasiliev gravities and find miraculous (almost) cancellations between bulk and edge characters that lead to agreement with the predictions of HS/CFT holography. We also discuss about the relation between the character integral representation and Rindler-AdS thermal partition function.


title: Black Hole Interior in Unitary Gauge Construction

abstract: A quantum system with a black hole accommodates two widely different, though physically equivalent, descriptions. In one description, based on global spacetime of general relativity, the existence of the interior region is manifest, while understanding unitarity requires nonperturbative quantum gravity effects such as replica wormholes. The other description adopts a manifestly unitary, or holographic, description, in which the interior emerges effectively as a collective phenomenon of fundamental degrees of freedom. In this paper we study the latter approach, which we refer to as the unitary gauge construction. In this picture, the formation of a black hole is signaled by the emergence of a surface (stretched horizon) possessing special dynamical properties: quantum chaos, fast scrambling, and low energy universality. These properties allow for constructing interior operators, as we do explicitly, without relying on details of microscopic physics. A key role is played by certain coarse modes in the zone region (hard modes), which determine the degrees of freedom relevant for the emergence of the interior. We study how the interior operators can or cannot be extended in the space of microstates and analyze irreducible errors associated with such extension. This reveals an intrinsic ambiguity of semiclassical theory formulated with a finite number of degrees of freedom. We provide an explicit prescription of calculating interior correlators in the effective theory, which describes only a finite region of spacetime. We study the issue of state dependence of interior operators in detail and discuss a connection of the resulting picture with the quantum error correction interpretation of holography.


title: Boundary Conditions in Topological AdS$_4$/CFT$_3$

abstract: We revisit the construction in four-dimensional gauged $Spin(4)$ supergravity of the holographic duals to topologically twisted three-dimensional $\mathcal{N}=4$ field theories. Our focus in this paper is to highlight some subtleties related to preserving supersymmetry in AdS/CFT, namely the inclusion of finite counterterms and the necessity of a Legendre transformation to find the dual to the field theory generating functional. Studying the geometry of these supergravity solutions, we conclude that the gravitational free energy is indeed independent from the metric of the boundary, and it vanishes for any smooth solution.


title: Observation of a non-Hermitian phase transition in an optical quantum gas

abstract: Quantum gases of light, as photons or polariton condensates in optical microcavities, are collective quantum systems enabling a tailoring of dissipation from e.g. cavity loss. This makes them a tool to study dissipative phases, an emerging subject in quantum manybody physics. Here we experimentally demonstrate a non-Hermitian phase transition of a photon Bose-Einstein condensate to a new dissipative phase, characterized by a biexponential decay of the condensate's second-order coherence. The phase transition occurs due to the emergence of an exceptional point in the quantum gas. While Bose-Einstein condensation is usually connected to ordinary lasing by a smooth crossover, the observed phase transition separates the novel, biexponential phase from both lasing and an intermediate, oscillatory condensate regime. Our findings pave the way for studies of a wide class of dissipative quantum phases, for instance in topological or lattice systems.


title: Logarithmic corrections to scaling in the four-dimensional uniform spanning tree

abstract: We compute the precise logarithmic corrections to mean-field scaling for various quantities describing the uniform spanning tree of the four-dimensional hypercubic lattice $\mathbb{Z}^4$. We are particularly interested in the distribution of the past of the origin, that is, the finite piece of the tree that is separated from infinity by the origin. We prove that the probability that the past contains a path of length $n$ is of order $(\log n)^{1/3}n^{-1}$, that the probability that the past contains at least $n$ vertices is of order $(\log n)^{1/6} n^{-1/2}$, and that the probability that the past reaches the boundary of the box $[-n,n]^4$ is of order $(\log n)^{2/3+o(1)}n^{-2}$. An important part of our proof is to prove concentration estimates for the capacity of the four-dimensional loop-erased random walk which may be of independent interest. Our results imply that the Abelian sandpile model also exhibits non-trivial polylogarithmic corrections to mean-field scaling in four dimensions, although it remains open to compute the precise order of these corrections.


title: RelationNet++: Bridging Visual Representations for Object Detection via Transformer Decoder

abstract: Existing object detection frameworks are usually built on a single format of object/part representation, i.e., anchor/proposal rectangle boxes in RetinaNet and Faster R-CNN, center points in FCOS and RepPoints, and corner points in CornerNet. While these different representations usually drive the frameworks to perform well in different aspects, e.g., better classification or finer localization, it is in general difficult to combine these representations in a single framework to make good use of each strength, due to the heterogeneous or non-grid feature extraction by different representations. This paper presents an attention-based decoder module similar as that in Transformer~\cite{vaswani2017attention} to bridge other representations into a typical object detector built on a single representation format, in an end-to-end fashion. The other representations act as a set of \emph{key} instances to strengthen the main \emph{query} representation features in the vanilla detectors. Novel techniques are proposed towards efficient computation of the decoder module, including a \emph{key sampling} approach and a \emph{shared location embedding} approach. The proposed module is named \emph{bridging visual representations} (BVR). It can perform in-place and we demonstrate its broad effectiveness in bridging other representations into prevalent object detection frameworks, including RetinaNet, Faster R-CNN, FCOS and ATSS, where about $1.5\sim3.0$ AP improvements are achieved. In particular, we improve a state-of-the-art framework with a strong backbone by about $2.0$ AP, reaching $52.7$ AP on COCO test-dev. The resulting network is named RelationNet++. The code will be available at https://github.com/microsoft/RelationNet2.


title: Proceedings 9th International Workshop on Theorem Proving Components for Educational Software

abstract: The 9th International Workshop on Theorem-Proving Components for Educational Software (ThEdu'20) was scheduled to happen on June 29 as a satellite of the IJCAR-FSCD 2020 joint meeting, in Paris. The COVID-19 pandemic came by surprise, though, and the main conference was virtualised. Fearing that an online meeting would not allow our community to fully reproduce the usual face-to-face networking opportunities of the ThEdu initiative, the Steering Committee of ThEdu decided to cancel our workshop. Given that many of us had already planned and worked for that moment, we decided that ThEdu'20 could still live in the form of an EPTCS volume. The EPTCS concurred with us, recognising this very singular situation, and accepted our proposal of organising a special issue with papers submitted to ThEdu'20. An open call for papers was then issued, and attracted five submissions, all of which have been accepted by our reviewers, who produced three careful reports on each of the contributions. The resulting revised papers are collected in the present volume. We, the volume editors, hope that this collection of papers will help further promoting the development of theorem-proving-based software, and that it will collaborate to improve the mutual understanding between computer mathematicians and stakeholders in education. With some luck, we would actually expect that the very special circumstances set up by the worst sanitary crisis in a century will happen to reinforce the need for the application of certified components and of verification methods for the production of educational software that would be available even when the traditional on-site learning experiences turn out not to be recommendable.


title: The realizability of discs with ribbons on a M\"obius strip

abstract: An hieroglyph on n letters is a cyclic sequence of the letters 1,2, . . . , n of length 2n such that each letter appears in the sequence twice.Take an hieroglyph H. Take a convex polygon with 2n sides. Put the letters in the sequence of letters of the hieroglyph on the sides of the convexpolygon in the same order. For each letter i glue the ends of a ribbon to thepair of sides corresponding to the letter i. Call the resulting surface a disk with ribbons corresponding to the hieroglyph H. An hieroglyph H is weakly realizable on the M\"obius strip if some disk with ribbons corresponding to H can be cut out of the M\"obius strip. We give a criterion for weak realizability, which gives a quadratic (in the number of letters) algorithm. Our criterion is based on the Mohar criterion for realizability of a disk with ribbons in the M\"obius strip.


title: Leading-Color Two-Loop QCD Corrections for Three-Photon Production at Hadron Colliders

abstract: We compute the two-loop helicity amplitudes for the production of three photons at hadron colliders in QCD at leading-color. Using the two-loop numerical unitarity method coupled with analytic reconstruction techniques, we obtain the decomposition of the two-loop amplitudes in terms of master integrals in analytic form. These expressions are valid to all orders in the dimensional regulator. We use them to compute the two-loop finite remainders, which are given in a form that can be efficiently evaluated across the whole physical phase space. We further package these results in a public code which assembles the helicity-summed squared two-loop remainders, whose numerical stability across phase-space is demonstrated. This is the first time that a five-point two-loop process is publicly available for immediate phenomenological applications.


title: Targeting for long-term outcomes

abstract: Decision makers often want to target interventions so as to maximize an outcome that is observed only in the long-term. This typically requires delaying decisions until the outcome is observed or relying on simple short-term proxies for the long-term outcome. Here we build on the statistical surrogacy and policy learning literatures to impute the missing long-term outcomes and then approximate the optimal targeting policy on the imputed outcomes via a doubly-robust approach. We first show that conditions for the validity of average treatment effect estimation with imputed outcomes are also sufficient for valid policy evaluation and optimization; furthermore, these conditions can be somewhat relaxed for policy optimization. We apply our approach in two large-scale proactive churn management experiments at The Boston Globe by targeting optimal discounts to its digital subscribers with the aim of maximizing long-term revenue. Using the first experiment, we evaluate this approach empirically by comparing the policy learned using imputed outcomes with a policy learned on the ground-truth, long-term outcomes. The performance of these two policies is statistically indistinguishable, and we rule out large losses from relying on surrogates. Our approach also outperforms a policy learned on short-term proxies for the long-term outcome. In a second field experiment, we implement the optimal targeting policy with additional randomized exploration, which allows us to update the optimal policy for future subscribers. Over three years, our approach had a net-positive revenue impact in the range of $4-5 million compared to the status quo.


title: Radio-frequency Dark Photon Dark Matter across the Sun

abstract: Dark photon as an ultralight dark matter candidate can interact with the Standard Model particles via kinetic mixing. We propose to search for the ultralight dark photon dark matter using radio telescopes with solar observations. The dark photon dark matter can efficiently convert into photons in the outermost region of the solar atmosphere, the solar corona, where the plasma mass of photons is close to the dark photon rest mass. Due to the strong resonant conversion and benefiting from the short distance between the Sun and the Earth, the radio telescopes can lead the dark photon search sensitivity in the mass range of $4 \times 10^{-8} - 4\times 10^{-6} \, \rm{eV}$, corresponding to the frequency $10 - 1000 \, {\rm MHz}$. As a promising example, the operating radio telescope LOFAR can reach the kinetic mixing $\epsilon \sim 10^{-13}$ ($10^{-14}$) within 1 (100) hour solar observations. The future experiment SKA phase 1 can reach $\epsilon \sim 10^{-16} - 10^{-14}$ with $1$ hour solar observations.


title: Collapsar R-Process Yields Can Reproduce [Eu/Fe] Abundance Scatter in Metal-Poor Stars

abstract: It is unclear if neutron star mergers can explain the observed r-process abundances of metal-poor stars. Collapsars, defined here as rotating massive stars whose collapse results in a rapidly accreting disk around a black hole that can launch jets, are a promising alternative. We find that we can produce a self-consistent model in which a population of collapsars with stochastic europium yields synthesizes all of the r-process material in metal-poor ([Fe/H] < -2.5) stars. Our model reproduces the observed scatter and evolution of scatter of [Eu/Fe] abundances. We find that if collapsars are the dominant r-process site for metal-poor stars, r-process synthesis may be linked to supernovae that produce long gamma-ray bursts. Our results also allow for the possibility that core-collapse supernovae beyond those that launch gamma-ray bursts also produce r-process material (e.g., potentially a subset of Type Ic-BL supernovae). Furthermore, we identify collapsar jet properties (isotropic energy, engine luminosity, or engine time) which may trace r-process yield and verify that the amount of r-process yield produced per collapsar in our model (~0.07 Msun) is consistent with other independent estimates. In the future, achieving 0.05 dex precision on distribution scatter or a reliable selection function would further constrain our probe of r-process production. Our model would also hold for another prompt r-process site with a power-law yield, and work is needed to determine if, for example, fast-merging neutron stars can also explain abundance scatter.