title: On cyclic and nontransitive probabilities

abstract: Motivated by classical nontransitivity paradoxes, we call an $n$-tuple $(x_1,\dots,x_n) \in[0,1]^n$ \textit{cyclic} if there exist independent random variables $U_1,\dots, U_n$ with $P(U_i=U_j)=0$ for $i\not=j$ such that $P(U_{i+1}>U_i)=x_i$ for $i=1,\dots,n-1$ and $P(U_1>U_n)=x_n$. We call the tuple $(x_1,\dots,x_n)$ \textit{nontransitive} if it is cyclic and in addition satisfies $x_i>1/2$ for all $i$. Let $p_n$ (resp.~$p_n^*$) denote the probability that a randomly chosen $n$-tuple $(x_1,\dots,x_n)\in[0,1]^n$ is cyclic (resp.~nontransitive). We determine $p_3$ and $p_3^*$ exactly, while for $n\ge4$ we give upper and lower bounds for $p_n$ that show that $p_n$ converges to $1$ as $n\to\infty$. We also determine the distribution of the smallest, middle, and largest elements in a cyclic triple.


title: A Riemannian Block Coordinate Descent Method for Computing the Projection Robust Wasserstein Distance

abstract: The Wasserstein distance has become increasingly important in machine learning and deep learning. Despite its popularity, the Wasserstein distance is hard to approximate because of the curse of dimensionality. A recently proposed approach to alleviate the curse of dimensionality is to project the sampled data from the high dimensional probability distribution onto a lower-dimensional subspace, and then compute the Wasserstein distance between the projected data. However, this approach requires to solve a max-min problem over the Stiefel manifold, which is very challenging in practice. The only existing work that solves this problem directly is the RGAS (Riemannian Gradient Ascent with Sinkhorn Iteration) algorithm, which requires to solve an entropy-regularized optimal transport problem in each iteration, and thus can be costly for large-scale problems. In this paper, we propose a Riemannian block coordinate descent (RBCD) method to solve this problem, which is based on a novel reformulation of the regularized max-min problem over the Stiefel manifold. We show that the complexity of arithmetic operations for RBCD to obtain an $\epsilon$-stationary point is $O(\epsilon^{-3})$. This significantly improves the corresponding complexity of RGAS, which is $O(\epsilon^{-12})$. Moreover, our RBCD has very low per-iteration complexity, and hence is suitable for large-scale problems. Numerical results on both synthetic and real datasets demonstrate that our method is more efficient than existing methods, especially when the number of sampled data is very large.


title: Elongation and percolation of defect motifs in anisotropic packing problems

abstract: We examine the regime between crystalline and amorphous packings of anisotropic objects on surfaces of different genus by continuously varying their size distribution or shape from monodispersed spheres to bidispersed mixtures or monodispersed ellipsoidal particles; we also consider an anisotropic variant of the Thomson problem with a mixture of charges. With increasing anisotropy, we first observe the disruption of translational order with an intermediate orientationally ordered hexatic phase as proposed by Nelson, Rubinstein and Spaepen, and then a transition to amorphous state. By analyzing the structure of the disclination motifs induced, we show that the hexatic-amorphous transition is caused by the growth and connection of disclination grain boundaries, suggesting this transition lies in the percolation universality class in the scenarios considered.


title: Long-time behaviour of a model for p62-ubiquitin aggregation in cellular autophagy

abstract: The qualitative behavior of a recently formulated ODE model for the dynamics of heterogenous aggregates is analyzed. Aggregates contain two types of particles, oligomers and cross-linkers. The motivation is a preparatory step of cellular autophagy, the aggregation of oligomers of the protein p62 in the presence of ubiquitin cross-linkers. A combination of explicit computations, formal asymptotics, and numerical simulations has led to conjectures on the bifurcation behavior, certain aspects of which are proven rigorously in this work. In particular, the stability of the zero state, where the model has a smoothness deficit is analyzed by a combination of regularizing transformations and blow-up techniques. On the other hand, in a different parameter regime, the existence of polynomially growing solutions is shown by Poincar\'e compactification, combined with a singular perturbation analysis .


title: Out-of-Equilibrium Dynamics and Excess Volatility in Firm Networks

abstract: We study the conditions under which input-output networks can dynamically attain a competitive equilibrium, where markets clear and profits are zero. We endow a classical firm network model with minimal dynamical rules that reduce supply/demand imbalances and excess profits. We show that the time needed to reach equilibrium diverges to infinity as the system approaches an instability point beyond which the Hawkins-Simons condition is violated and competitive equilibrium is no longer admissible. We argue that such slow dynamics is a source of excess volatility, through accumulation and amplification of exogenous shocks. Factoring in essential physical constraints absent in our minimal model, such as causality or inventory management, we then propose a dynamically consistent model that displays a rich variety of phenomena. Competitive equilibrium can only be reached after some time and within some restricted region of parameter space, outside of which one observes spontaneous periodic and chaotic dynamics, reminiscent of real business cycles. This suggests an alternative explanation of excess volatility in terms of purely endogenous fluctuations. Diminishing return to scale and increased perishability of goods are found to ease convergence towards equilibrium.


title: Quantization and soliton-like solutions for the $\Phi\Psi$-model in a optic fiber

abstract: In the framework of a mesoscopical model for dielectric media we provide an analytical description for the electromagnetic field confined in a cylindrical cavity containing a finite dielectric sample. This system is apted to simulate the electromagnetic field in a optic fiber, in which two different regions, a vacuum region and a dielectric one, appear. A complete description for the scattering basis is introduced, together with field quantization and the two-point function. Furthermore, we also determine soliton-like solutions in the dielectric, propagating in the sample of nonlinear dielectric medium.


abstract: In this paper subdivision schemes, which are used for functions approximation and curves generation, are considered. In classical case, for the functions defined on the real line, the theory of subdivision schemes is widely known due to multiple applications in constructive approximation theory, signal processing as well as for generating fractal curves and surfaces. Subdivision schemes on a dyadic half-line, which is a positive half-line, equipped with the standard Lebesgue measure and the digitwise binary addition operation, where the Walsh functions play the role of exponents, are defined and studied. Necessary and sufficient convergence conditions of the subdivision schemes in terms of spectral properties of matrices and in terms of the smoothness of the solution of the corresponding refinement equation are proved. The problem of the convergence of subdivision schemes with non-negative coefficients is also investigated. Explicit convergence criterion of the subdivision schemes with four coefficients is obtained. As an auxiliary result fractal curves on a dyadic half-line are defined and the formula of their smoothness is proved. The paper contains various illustrations and numerical results.


title: Tactile Object Pose Estimation from the First Touch with Geometric Contact Rendering

abstract: In this paper, we present an approach to tactile pose estimation from the first touch for known objects. First, we create an object-agnostic map from real tactile observations to contact shapes. Next, for a new object with known geometry, we learn a tailored perception model completely in simulation. To do so, we simulate the contact shapes that a dense set of object poses would produce on the sensor. Then, given a new contact shape obtained from the sensor output, we match it against the pre-computed set using the object-specific embedding learned purely in simulation using contrastive learning. This results in a perception model that can localize objects from a single tactile observation. It also allows reasoning over pose distributions and including additional pose constraints coming from other perception systems or multiple contacts. We provide quantitative results for four objects. Our approach provides high accuracy pose estimations from distinctive tactile observations while regressing pose distributions to account for those contact shapes that could result from different object poses. We further extend and test our approach in multi-contact scenarios where several tactile sensors are simultaneously in contact with the object. Website: http://mcube.mit.edu/research/tactile_loc_first_touch.html


title: Disk-jet coupling changes as a possible indicator for outbursts from GX 339-4 remaining within the X-ray hard state

abstract: We present quasi-simultaneous radio, (sub-)millimetre, and X-ray observations of the Galactic black hole X-ray binary GX 339-4, taken during its 2017--2018 outburst, where the source remained in the hard X-ray spectral state. During this outburst, GX 339-4 showed no atypical X-ray behaviour that may act as a indicator for an outburst remaining within the hard state. However, quasi-simultaneous radio and X-ray observations showed a flatter than expected coupling between the radio and X-ray luminosities (with a best fit relation of $L_{\rm radio} \propto L_{\rm X}^{0.39 \pm 0.06}$), when compared to successful outbursts from this system ($L_{\rm radio} \propto L_{\rm X}^{0.62 \pm 0.02}$). While our 2017--2018 outburst data only span a limited radio and X-ray luminosity range ($\sim$1 order of magnitude in both, where more than 2-orders of magnitude in $L_{\rm X}$ is desired), including data from other hard-only outbursts from GX 339-4 extends the luminosity range to $\sim$1.2 and $\sim$2.8 orders of magnitude, respectively, and also results in a flatter correlation (where $L_{\rm radio} \propto L_{\rm X}^{0.46 \pm 0.04}$). This result is suggestive that for GX 339-4 a flatter radio -- X-ray correlation, implying a more inefficient coupling between the jet and accretion flow, could act as an indicator for a hard-only outburst. However, further monitoring of both successful and hard-only outbursts over larger luminosity ranges with strictly simultaneous radio and X-ray observations is required from different, single sources, to explore if this applies generally to the population of black hole X-ray binaries, or even GX 339-4 at higher hard-state luminosities.


title: Uncertainty Intervals for Graph-based Spatio-Temporal Traffic Prediction

abstract: Many traffic prediction applications rely on uncertainty estimates instead of the mean prediction. Statistical traffic prediction literature has a complete subfield devoted to uncertainty modelling, but recent deep learning traffic prediction models either lack this feature or make specific assumptions that restrict its practicality. We propose Quantile Graph Wavenet, a Spatio-Temporal neural network that is trained to estimate a density given the measurements of previous timesteps, conditioned on a quantile. Our method of density estimation is fully parameterised by our neural network and does not use a likelihood approximation internally. The quantile loss function is asymmetric and this makes it possible to model skewed densities. This approach produces uncertainty estimates without the need to sample during inference, such as in Monte Carlo Dropout, which makes our method also efficient.


title: On the Binding Problem in Artificial Neural Networks

abstract: Contemporary neural networks still fall short of human-level generalization, which extends far beyond our direct experiences. In this paper, we argue that the underlying cause for this shortcoming is their inability to dynamically and flexibly bind information that is distributed throughout the network. This binding problem affects their capacity to acquire a compositional understanding of the world in terms of symbol-like entities (like objects), which is crucial for generalizing in predictable and systematic ways. To address this issue, we propose a unifying framework that revolves around forming meaningful entities from unstructured sensory inputs (segregation), maintaining this separation of information at a representational level (representation), and using these entities to construct new inferences, predictions, and behaviors (composition). Our analysis draws inspiration from a wealth of research in neuroscience and cognitive psychology, and surveys relevant mechanisms from the machine learning literature, to help identify a combination of inductive biases that allow symbolic information processing to emerge naturally in neural networks. We believe that a compositional approach to AI, in terms of grounded symbol-like representations, is of fundamental importance for realizing human-level generalization, and we hope that this paper may contribute towards that goal as a reference and inspiration.


title: The spaces of dyadic distributions

abstract: In this paper subdivision schemes, which are used for functions approximation and curves generation, are considered. In classical case, for the functions defined on the real line, the theory of subdivision schemes is widely known due to multiple applications in constructive approximation theory, signal processing as well as for generating fractal curves and surfaces. Subdivision schemes on a dyadic half-line, which is the positive half-line, equipped with the standard Lebesgue measure and the digitwise binary addition operation, where the Walsh functions play the role of exponents, are defined and studied. Necessary and sufficient convergence conditions of the subdivision schemes in terms of spectral properties of matrices and in terms of the smoothness of the solution of the corresponding refinement equation are proved. The problem of the convergence of subdivision schemes with non-negative coefficients is also investigated. Explicit convergence criterion of the subdivision schemes with four coefficients is obtained. As an auxiliary result fractal curves on a dyadic half-line are defined and the formula of their smoothness is proved. The paper contains various illustrations and numerical results.


title: Data embedding and prediction by sparse tropical matrix factorization

abstract: Matrix factorization methods are linear models, with limited capability to model complex relations. In our work, we use tropical semiring to introduce non-linearity into matrix factorization models. We propose a method called Sparse Tropical Matrix Factorization (STMF) for the estimation of missing (unknown) values. We evaluate the efficiency of the STMF method on both synthetic data and biological data in the form of gene expression measurements downloaded from The Cancer Genome Atlas (TCGA) database. Tests on unique synthetic data showed that STMF approximation achieves a higher correlation than non-negative matrix factorization (NMF), which is unable to recover patterns effectively. On real data, STMF outperforms NMF on six out of nine gene expression datasets. While NMF assumes normal distribution and tends toward the mean value, STMF can better fit to extreme values and distributions. STMF is the first work that uses tropical semiring on sparse data. We show that in certain cases semirings are useful because they consider the structure, which is different and simpler to understand than it is with standard linear algebra.


title: Synthesis to Deployment: Cyber-Physical Control Architectures

abstract: We consider the problem of how to deploy a controller to a (networked) cyber-physical system (CPS). Controlling a CPS is an involved task, and synthesizing a controller to respect sensing, actuation, and communication constraints is only part of the challenge. In addition to controller synthesis, one should also consider how the controller will be incorporated within the CPS. Put another way, the cyber layer and its interaction with the physical layer need to be taken into account. In this work, we aim to bridge the gap between theoretical controller synthesis and practical CPS deployment. We adopt the system level synthesis (SLS) framework to synthesize a controller, either state-feedback or output-feedback, and provide deployment architectures for the standard SLS controllers. Furthermore, we derive new controller realizations for open-loop stable systems and introduce different state-feedback and output-feedback architectures for deployment, ranging from fully centralized to fully distributed. Finally, we compare the trade-offs among them in terms of robustness, memory, computation, and communication overhead.


title: Towards a Probabilistic Foundation of Relativistic Quantum Theory: The One-Body Born Rule in Curved Spacetime

abstract: In this work we establish a novel approach to the foundations of relativistic quantum theory, which is based on generalizing the quantum-mechanical Born rule for determining particle position probabilities to curved spacetime. A principal motivator for this research has been to overcome internal mathematical problems of quantum field theory (QFT) such as the `problem of infinities' (renormalization), which axiomatic approaches to QFT have shown to be not only of mathematical but also of conceptual nature. The approach presented here is probabilistic by construction, can accommodate a wide array of dynamical models, does not rely on the symmetries of Minkowski spacetime, and respects the general principle of relativity. In the analytical part of this work we consider the $1$-body case under the assumption of smoothness of the mathematical quantities involved. This is identified as a special case of the theory of the general-relativistic continuity equation. While related approaches to the relativistic generalization of the Born rule assume the hypersurfaces of interest to be spacelike and the spacetime to be globally hyperbolic, we employ prior contributions by C. Eckart and J. Ehlers to show that the former condition is naturally replaced by a transversality condition and that the latter one is obsolete. We discuss two distinct formulations of the $1$-body case, which, borrowing terminology from the non-relativistic analog, we term the Lagrangian and Eulerian pictures. We provide a comprehensive treatment of both. The main contribution of this work to the mathematical physics literature is the development of the Lagrangian picture. The Langrangian picture shows how one can address the `problem of time' in this approach and therefore serves as a blueprint for the generalization to many bodies and the case that the number of bodies is not conserved (example given for the latter).


title: Participatory Budgeting with Project Groups

abstract: We study a generalization of the standard approval-based model of participatory budgeting (PB), in which voters are providing approval ballots over a set of predefined projects and -- in addition to a global budget limit, there are several groupings of the projects, each group with its own budget limit. We study the computational complexity of identifying project bundles that maximize voter satisfaction while respecting all budget limits. We show that the problem is generally intractable and describe efficient exact algorithms for several special cases, including instances with only few groups and instances where the group structure is close to be hierarchical, as well as efficient approximation algorithms. Our results could allow, e.g., municipalities to hold richer PB processes that are thematically and geographically inclusive.


title: E3D: Event-Based 3D Shape Reconstruction

abstract: 3D shape reconstruction is a primary component of augmented/virtual reality. Despite being highly advanced, existing solutions based on RGB, RGB-D and Lidar sensors are power and data intensive, which introduces challenges for deployment in edge devices. We approach 3D reconstruction with an event camera, a sensor with significantly lower power, latency and data expense while enabling high dynamic range. While previous event-based 3D reconstruction methods are primarily based on stereo vision, we cast the problem as multi-view shape from silhouette using a monocular event camera. The output from a moving event camera is a sparse point set of space-time gradients, largely sketching scene/object edges and contours. We first introduce an event-to-silhouette (E2S) neural network module to transform a stack of event frames to the corresponding silhouettes, with additional neural branches for camera pose regression. Second, we introduce E3D, which employs a 3D differentiable renderer (PyTorch3D) to enforce cross-view 3D mesh consistency and fine-tune the E2S and pose network. Lastly, we introduce a 3D-to-events simulation pipeline and apply it to publicly available object datasets and generate synthetic event/silhouette training pairs for supervised learning.


title: Electron-beam interaction with emission-line clouds in blazars

abstract: Context: An electron-positron beam escaping from the magnetospheric vacuum gap of an accreting black hole interacts with recombination-line photons from surrounding gas clouds. Inverse-Compton scattering and subsequent pair production initiate unsaturated electromagnetic cascades exhibiting a characteristic spectral energy distribution. Aims: By modelling the interactions of beam electrons (positrons) with hydrogen and helium recombination-line photons, we seek to describe the spectral signature of beam-driven cascades in the broad emission-line region of blazar jets. Methods: Employing coupled kinetic equations for electrons (positrons) and photons including an escape term, we numerically obtain their steady-state distributions, and the escaping photon spectrum. Results: We find that cascade emission resulting from beam interactions can produce a narrow spectral feature at TeV energies. Indications of such an intermittent feature, which defies an explanation in the standard shock-in-jet scenario, have been found at $\approx\,4\,\sigma$ confidence level at an energy of $\approx$ 3 TeV in the spectrum of the blazar Mrk 501. Conclusions: The energetic requirements for explaining the intermittent 3 TeV bump with the beam-interaction model are plausible: Gap discharges that lead to multi-TeV beam electrons (positrons) carrying $\approx$ 0.1 % of the Blandford-Znajek luminosity, which interact with recombination-line photons from gas clouds that reprocess $\approx$ 1 % of the similar accretion luminosity are required.


title: Extracting the signed backbone of intrinsically dense weighted networks

abstract: Networks provide useful tools for analyzing diverse complex systems from natural, social, and technological domains. Growing size and variety of data such as more nodes and links and associated weights, directions, and signs can provide accessory information. Link and weight abundance, on the other hand, results in denser networks with noisy, insignificant, or otherwise redundant data. Moreover, typical network analysis and visualization techniques presuppose sparsity and are not appropriate or scalable for dense and weighted networks. As a remedy, network backbone extraction methods aim to retain only the important links while preserving the useful and elucidative structure of the original networks for further analyses. Here, we provide the first methods for extracting signed network backbones from intrinsically dense unsigned unipartite weighted networks. Utilizing a null model based on statistical techniques, the proposed significance filter and vigor filter allow inferring edge signs. Empirical analysis on migration, voting, temporal interaction, and species similarity networks reveals that the proposed filters extract meaningful and sparse signed backbones while preserving the multiscale nature of the network. The resulting backbones exhibit characteristics typically associated with signed networks such as reciprocity, structural balance, and community structure. The developed tool is provided as a free, open-source software package.


title: Positional Encoding as Spatial Inductive Bias in GANs

abstract: SinGAN shows impressive capability in learning internal patch distribution despite its limited effective receptive field. We are interested in knowing how such a translation-invariant convolutional generator could capture the global structure with just a spatially i.i.d. input. In this work, taking SinGAN and StyleGAN2 as examples, we show that such capability, to a large extent, is brought by the implicit positional encoding when using zero padding in the generators. Such positional encoding is indispensable for generating images with high fidelity. The same phenomenon is observed in other generative architectures such as DCGAN and PGGAN. We further show that zero padding leads to an unbalanced spatial bias with a vague relation between locations. To offer a better spatial inductive bias, we investigate alternative positional encodings and analyze their effects. Based on a more flexible positional encoding explicitly, we propose a new multi-scale training strategy and demonstrate its effectiveness in the state-of-the-art unconditional generator StyleGAN2. Besides, the explicit spatial inductive bias substantially improve SinGAN for more versatile image manipulation.


title: Highly polar molecules consisting of a copper or silver atom interacting with an alkali-metal or alkaline-earth-metal atom

abstract: We theoretically investigate the properties of highly polar diatomic molecules containing $^2S$-state transition-metal atoms. We calculate potential energy curves, permanent electric dipole moments, spectroscopic constants, and leading long-range dispersion-interaction coefficients for molecules consisting of either a Cu or Ag atom interacting with an alkali-metal (Li, Na, K, Rb, Cs, Fr) or alkaline-earth-metal (Be, Mg, Ca, Sr, Ba, Ra) atom. We use ab initio electronic structure methods, such as the coupled cluster and configuration interaction ones, with large Gaussian basis sets and small-core relativistic energy-consistent pseudopotentials. We predict that the studied molecules in the ground electronic state are strongly bound with highly polarized covalent or ionic bonds resulting in very large permanent electric dipole moments. We find that highly excited vibrational levels have maximal electric dipole moments, e.g., exceeding 13$\,$debye for CsAg and 6$\,$debye for BaAg. Results for Cu$_2$, Ag$_2$, and CuAg are also reported. The studied molecules may find application in ultracold dipolar many-body physics, controlled chemistry, or precision measurement experiments.


title: Parametric measures of variability induced by risk measures

abstract: We present a general framework for a comparative theory of variability measures, with a particular focus on the recently introduced one-parameter families of inter-Expected Shortfall differences and inter-expectile differences, that are explored in detail and compared with the widely known and applied inter-quantile differences. From the mathematical point of view, our main result is a characterization of symmetric and comonotonic variability measures as mixtures of inter-Expected Shortfall differences, under a few additional technical conditions. Further, we study the stochastic orders induced by the pointwise comparison of inter-Expected Shortfall and inter-expectile differences, and discuss their relationship with the dilation order. From the statistical point of view, we establish asymptotic consistency and normality of the natural estimators and provide a rule of the thumb for cross-comparisons. Finally, we study the empirical behaviour of the considered classes of variability measures on the S&P 500 Index under various economic regimes, and explore the comparability of different time series according to the introduced stochastic orders.


title: Estimating distances from parallaxes. V: Geometric and photogeometric distances to 1.47 billion stars in Gaia Early Data Release 3

abstract: Stellar distances constitute a foundational pillar of astrophysics. The publication of 1.47 billion stellar parallaxes from Gaia is a major contribution to this. Yet despite Gaia's precision, the majority of these stars are so distant or faint that their fractional parallax uncertainties are large, thereby precluding a simple inversion of parallax to provide a distance. Here we take a probabilistic approach to estimating stellar distances that uses a prior constructed from a three-dimensional model of our Galaxy. This model includes interstellar extinction and Gaia's variable magnitude limit. We infer two types of distance. The first, geometric, uses the parallax together with a direction-dependent prior on distance. The second, photogeometric, additionally uses the colour and apparent magnitude of a star, by exploiting the fact that stars of a given colour have a restricted range of probable absolute magnitudes (plus extinction). Tests on simulated data and external validations show that the photogeometric estimates generally have higher accuracy and precision for stars with poor parallaxes. We provide a catalogue of 1.47 billion geometric and 1.35 billion photogeometric distances together with asymmetric uncertainty measures. Our estimates are quantiles of a posterior probability distribution, so they transform invariably and can therefore also be used directly in the distance modulus (5log10(r)-5). The catalogue may be downloaded or queried using ADQL at various sites (see http://www.mpia.de/homes/calj/gedr3_distances.html) where it can also be cross-matched with the Gaia catalogue.


title: First measurement using a nuclear emulsion detector of the $\nu_{\mu}$ charged-current cross section on iron around the 1$\,$GeV energy region

abstract: We have carried out $\nu_{\mu}$ charged-current interaction measurement on iron using an emulsion detector exposed to the T2K neutrino beam in the J-PARC neutrino facility. The data samples correspond to 4.0$\times$10$^{19}$ protons on target, and the neutrino mean energy is 1.49$\,$GeV. The emulsion detector is suitable for precision measurements of charged particles produced in neutrino-iron interactions with a low momentum threshold thanks to thin-layered structure and sub-$\mu$m spatial resolution. The charged particles are successfully detected, and their multiplicities are measured using the emulsion detector. The cross section was measured to be $\sigma^{\mathrm{Fe}}_{\mathrm{CC}} = (1.28 \pm 0.11({\mathrm{stat.}})^{+0.12}_{-0.11}({\mathrm{syst.}})) \times 10^{-38} \, {\mathrm{cm}}^{2}/{\mathrm{nucleon}}$. The cross section in a limited kinematic phase space of induced muons, $\theta_{\mu} < 45^{\circ}$ and $p_{\mu} > 400 \, {\rm MeV}/c$, on iron was $\sigma^{\mathrm{Fe}}_{\mathrm{CC \hspace{1mm} phase \hspace{0.5mm} space}} = (0.84 \pm 0.07({\mathrm{stat.}})^{+0.07}_{-0.06}({\mathrm{syst.}})) \times 10^{-38} \, {\mathrm{cm}}^{2}/{\mathrm{nucleon}}$. The cross-section results are consistent with previous values obtained via different techniques using the same beamline, and they are well reproduced by current neutrino interaction models. These results demonstrate the capability of the detector towards the detailed measurements of the neutrino-nucleus interactions around the 1$\,$GeV energy region.


title: Isomorphic Bisections of Cubic Graphs

abstract: Graph partitioning, or the dividing of a graph into two or more parts based on certain conditions, arises naturally throughout discrete mathematics, and problems of this kind have been studied extensively. In the 1990s, Ando conjectured that the vertices of every cubic graph can be partitioned into two parts that induce isomorphic subgraphs. Using probabilistic methods together with delicate recolouring arguments, we prove Ando's conjecture for large connected graphs.


title: Atomistic-to-Meso Multi-Scale Data-Driven Graph Surrogate Modeling of Dislocation Glide

abstract: From their birth in the manufacturing process, materials inherently contain defects that affect the mechanical behavior across multiple length and time-scales, including vacancies, dislocations, voids and cracks. Understanding, modeling, and real-time simulation of the underlying stochastic micro-structure defect evolution is therefore vital towards multi-scale coupling and propagating numerous sources of uncertainty from atomistic to eventually aging continuum mechanics. We develop a graph-based surrogate model of dislocation glide for computation of dislocation mobility. We model an edge dislocation as a random walker, jumping between neighboring nodes of a graph following a Poisson stochastic process. The network representation functions as a coarse-graining of a molecular dynamics simulation that provides dislocation trajectories for an empirical computation of jump rates. With this construction, we recover the original atomistic mobility estimates, with remarkable computational speed-up and accuracy. Furthermore, the underlying stochastic process provides the statistics of dislocation mobility associated to the original molecular dynamics simulation, allowing an efficient propagation of material parameters and uncertainties across the scales.


title: Matrix integrals $\&$ finite holography

abstract: We explore the conjectured duality between a class of large $N$ matrix integrals, known as multicritical matrix integrals (MMI), and the series $(2m-1,2)$ of non-unitary minimal models on a fluctuating background. We match the critical exponents of the leading order planar expansion of MMI, to those of the continuum theory on an $S^2$ topology. From the MMI perspective this is done both through a multi-vertex diagrammatic expansion, thereby revealing novel combinatorial expressions, as well as through a systematic saddle point evaluation of the matrix integral as a function of its parameters. From the continuum point of view the corresponding critical exponents are obtained upon computing the partition function in the presence of a given conformal primary. Further to this, we elaborate on a Hilbert space of the continuum theory, and the putative finiteness thereof, on both an $S^2$ and a $T^2$ topology using BRST cohomology considerations. Matrix integrals support this finiteness.


title: MorphGAN: One-Shot Face Synthesis GAN for Detecting Recognition Bias

abstract: To detect bias in face recognition networks, it can be useful to probe a network under test using samples in which only specific attributes vary in some controlled way. However, capturing a sufficiently large dataset with specific control over the attributes of interest is difficult. In this work, we describe a simulator that applies specific head pose and facial expression adjustments to images of previously unseen people. The simulator first fits a 3D morphable model to a provided image, applies the desired head pose and facial expression controls, then renders the model into an image. Next, a conditional Generative Adversarial Network (GAN) conditioned on the original image and the rendered morphable model is used to produce the image of the original person with the new facial expression and head pose. We call this conditional GAN -- MorphGAN. Images generated using MorphGAN conserve the identity of the person in the original image, and the provided control over head pose and facial expression allows test sets to be created to identify robustness issues of a facial recognition deep network with respect to pose and expression. Images generated by MorphGAN can also serve as data augmentation when training data are scarce. We show that by augmenting small datasets of faces with new poses and expressions improves the recognition performance by up to 9% depending on the augmentation and data scarcity.


title: Isotope effects on molecular structures and electronic properties of liquid water via deep potential molecular dynamics based on SCAN functional

abstract: Feynman path-integral deep potential molecular dynamics (PI-DPMD) calculations have been employed to study both light (H$_2$O) and heavy water (D$_2$O) within the isothermal-isobaric ensemble. In particular, the deep neural network is trained based on ab initio data obtained from the strongly constrained and appropriately normed (SCAN) exchange-correlation functional. Because of the lighter mass of hydrogen than deuteron, the properties of light water is more influenced by nuclear quantum effect than those of heavy water. Clear isotope effects are observed and analyzed in terms of hydrogen-bond structure and electronic properties of water that are closely associated with experimental observables. The molecular structures of both liquid H$_2$O and D$_2$O agree well with the data extracted from scattering experiments. The delicate isotope effects on radial distribution functions and angular distribution functions are well reproduced as well. Our approach demonstrates that deep neural network combined with SCAN functional based ab initio molecular dynamics provides an accurate theoretical tool for modeling water and its isotope effects.


title: Conformal Conserved Currents in Embedding Space

abstract: We study conformal conserved currents in arbitrary irreducible representations of the Lorentz group using the embedding space formalism. With the help of the operator product expansion, we first show that conservation conditions can be fully investigated by considering only two- and three-point correlation functions. We then find an explicitly conformally-covariant differential operator in embedding space that implements conservation based on the standard position space operator product expansion differential operator $\partial_\mu$, although the latter does not uplift to embedding space covariantly. The differential operator in embedding space that imposes conservation is the same differential operator $\mathcal{D}_{ijA}$ used in the operator product expansion in embedding space. We provide several examples including conserved currents in irreducible representations that are not symmetric and traceless. With an eye on four-point conformal bootstrap equations for four conserved vector currents $\langle JJJJ\rangle$ and four energy-momentum tensors $\langle TTTT\rangle$, we mostly focus on conservation conditions for $\langle JJ\mathcal{O}\rangle$ and $\langle TT\mathcal{O}\rangle$. Finally, we reproduce and extend the consequences of conformal Ward identities at coincident points by determining three-point coefficients in terms of charges.


title: Video Deblurring by Fitting to Test Data

abstract: Motion blur in videos captured by autonomous vehicles and robots can degrade their perception capability. In this work, we present a novel approach to video deblurring by fitting a deep network to the test video. Our key observation is that some frames in a video with motion blur are much sharper than others, and thus we can transfer the texture information in those sharp frames to blurry frames. Our approach heuristically selects sharp frames from a video and then trains a convolutional neural network on these sharp frames. The trained network often absorbs enough details in the scene to perform deblurring on all the video frames. As an internal learning method, our approach has no domain gap between training and test data, which is a problematic issue for existing video deblurring approaches. The conducted experiments on real-world video data show that our model can reconstruct clearer and sharper videos than state-of-the-art video deblurring approaches. Code and data are available at https://github.com/xrenaa/Deblur-by-Fitting.


title: Quantum Buzzwords

abstract: Many scientists seeking to understand the quantum mechanics of measurement situations (Copenhagen quantum theory) agree on its overwhelmingly successful algorithms to predict the outcomes of laboratory measurements but disagree on what these algorithms mean and how they are to be interpreted. Some of these problems are briefly described and resolutions suggested from the decoherent (or consistent) histories quantum mechanics of closed systems like the Universe.


title: Disconnection and entropic repulsion for the harmonic crystal with random conductances

abstract: We study level-set percolation for the harmonic crystal on $\mathbb{Z}^d$, $d \geq 3$, with uniformly elliptic random conductances. We prove that this model undergoes a non-trivial phase transition at a critical level that is almost surely constant under the environment measure. Moreover, we study the disconnection event that the level-set of this field below a level $\alpha$ disconnects the discrete blow-up of a compact set $A \subseteq \mathbb{R}^d$ from the boundary of an enclosing box. We obtain quenched asymptotic upper and lower bounds on its probability in terms of the homogenized capacity of $A$, utilizing results from Neukamm, Sch\"affner and Schl\"omerkemper, see arXiv:1606.06533. Furthermore, we give upper bounds on the probability that a local average of the field deviates from some profile function depending on $A$, when disconnection occurs. The upper and lower bounds concerning disconnection that we derive are plausibly matching at leading order. In this case, this work shows that conditioning on disconnection leads to an entropic push-down of the field. The results in this article generalize the findings of arXiv:1802.02518 and arXiv:1808.09947 by the authors which treat the case of constant conductances. Our proofs involve novel "solidification estimates" for random walks, which are similar in nature to the corresponding estimates for Brownian motion derived by Sznitman and the second author in arXiv:1706.07229.


title: The Fulib Solution to the TTC 2020 Migration Case

abstract: At Kassel University we are working on a solution for bidirectional transformations based on event sourcing for about a year, now. It turned out, that the TTC 2020 migration case is a special case of a bidirectional transformation and that our approach provides a reasonable solution for it.


title: Gauge equivariant neural networks for quantum lattice gauge theories

abstract: Gauge symmetries play a key role in physics appearing in areas such as quantum field theories of the fundamental particles and emergent degrees of freedom in quantum materials. Motivated by the desire to efficiently simulate many-body quantum systems with exact local gauge invariance, gauge equivariant neural-network quantum states are introduced, which exactly satisfy the local Hilbert space constraints necessary for the description of quantum lattice gauge theory with Zd gauge group on different geometries. Focusing on the special case of Z2 gauge group on a periodically identified square lattice, the equivariant architecture is analytically shown to contain the loop-gas solution as a special case. Gauge equivariant neural-network quantum states are used in combination with variational quantum Monte Carlo to obtain compact descriptions of the ground state wavefunction for the Z2 theory away from the exactly solvable limit, and to demonstrate the confining/deconfining phase transition of the Wilson loop order parameter.

categories: cond-mat.str-el cond-mat.dis-nn cs.LG hep-lat quant-ph


title: The Role of Symmetry in Quantum Query-to-Communication Simulation

abstract: Buhrman, Cleve and Wigderson (STOC'98) showed that for every Boolean function f : {-1,1}^n to {-1,1} and G in {AND_2, XOR_2}, the bounded-error quantum communication complexity of the composed function f o G equals O(Q(f) log n), where Q(f) denotes the bounded-error quantum query complexity of f. This is achieved by Alice running the optimal quantum query algorithm for f, using a round of O(log n) qubits of communication to implement each query. This is in contrast with the classical setting, where it is easy to show that R^{cc}(f o G) is at most 2R(f), where R^{cc} and R denote bounded-error communication and query complexity, respectively. We show that the O(log n) overhead is required for some functions in the quantum setting, and thus the BCW simulation is tight. We note here that prior to our work, the possibility of Q^{cc}(f o G) = O(Q(f)), for all f and all G in {AND_2, XOR_2}, had not been ruled out. More specifically, we show the following. - We show that the log n overhead is *not* required when f is symmetric, generalizing a result of Aaronson and Ambainis for the Set-Disjointness function (Theory of Computing'05). - In order to prove the above, we design an efficient distributed version of noisy amplitude amplification that allows us to prove the result when f is the OR function. - In view of our first result above, one may ask whether the log n overhead in the BCW simulation can be avoided even when f is transitive, which is a weaker notion of symmetry. We give a strong negative answer by showing that the log n overhead is still necessary for some transitive functions even when we allow the quantum communication protocol an error probability that can be arbitrarily close to 1/2. - We also give, among other things, a general recipe to construct functions for which the log n overhead is required in the BCW simulation in the bounded-error communication model.


title: High Fermi velocities and small cyclotron masses in LaAlGe

abstract: We report quantum oscillation measurements of LaAlGe, a Lorentz-violating type-II Weyl semimetal with tilted Weyl cones. Very small quasiparticle masses and very high Fermi velocities were detected at the Fermi surface. Whereas three main frequencies have been observed, angular dependence of two Fermi surface sheets indicates possible two-dimensional (2D) character despite the absence of the 2D structural features such as van der Waals bonds. Such conducting states may offer a good platform for low-dimensional polarized spin current in magnetic RAlGe (R = Ce, Pr) materials.


title: $\mathbb{Z}_2$ lattice gauge theories and Kitaev's toric code: A scheme for analog quantum simulation

abstract: Kitaev's toric code is an exactly solvable model with $\mathbb{Z}_2$-topological order, which has potential applications in quantum computation and error correction. However, a direct experimental realization remains an open challenge. Here, we propose a building block for $\mathbb{Z}_2$ lattice gauge theories coupled to dynamical matter and demonstrate how it allows for an implementation of the toric-code ground state and its topological excitations. This is achieved by introducing separate matter excitations on individual plaquettes, whose motion induce the required plaquette terms. The proposed building block is realized in the second-order coupling regime and is well suited for implementations with superconducting qubits. Furthermore, we propose a pathway to prepare topologically non-trivial initial states during which a large gap on the order of the underlying coupling strength is present. This is verified by both analytical arguments and numerical studies. Moreover, we outline experimental signatures of the ground-state wavefunction and introduce a minimal braiding protocol. Detecting a $\pi$-phase shift between Ramsey fringes in this protocol reveals the anyonic excitations of the toric-code Hamiltonian in a system with only three triangular plaquettes. Our work paves the way for realizing non-Abelian anyons in analog quantum simulators.


title: A lower bound limiting solutions in the hyperbolic case of the generalized Fermat equation

abstract: We find a lower bound for $\chi = 1/p+1/q+1/r$ limiting any solution in the hyperbolic case of the Generalized Fermat Equation $x^p + y^q = z^r$.


title: Applications of Mean Field Games in Financial Engineering and Economic Theory

abstract: This is an expanded version of the lecture given at the AMS Short Course on Mean Field Games, on January 13, 2020 in Denver CO. The assignment was to discuss applications of Mean Field Games in finance and economics. I need to admit upfront that several of the examples reviewed in this chapter were already discussed in book form. Still, they are here accompanied with discussions of, and references to, works which appeared over the last three years. Moreover, several completely new sections are added to show how recent developments in financial engineering and economics can benefit from being viewed through the lens of the Mean Field Game paradigm. The new financial engineering applications deal with bitcoin mining and the energy markets, while the new economic applications concern models offering a smooth transition between macro-economics and finance, and contract theory.


title: Determining the Variations of Ca-K index and Features using a Century Long Equal Contrast Images from Kodaikanal Observatory

abstract: In the earlier analysis of Ca-K spectroheliograms obtained at Kodaikanal Observatory, the "Good" images were used to investigate variations in the chromosphere. Still, the contrast of the images varied on a day-to-day basis. We developed a new methodology to generate images to form a uniform time series. We adjusted each image's contrast until the FWHM of the normalized intensity distribution attained a value between 0.10 and 0.11. This methodology of the "Equal Contrast technique" is expected to compensate for the change of emulsion, development, contrast of the images due to centering of Ca-K line on the exit slit and sky transparency. Besides, this procedure will correct variations in density-to-intensity conversion for different images. We find that the correlation between sunspot and Ca-K line data improves by a large amount. For example, correlation coefficient (CC) between monthly averaged sunspots and Ca-K plage areas for the equal contrast data improves to 0.9 compared to 0.75 for the "Good" data with unequal contrast. The CC for equal contrast images improves to $\sim$0.78 from $\sim$0.46 for the "Okay" data with unequal contrast. Even the CC between the plage area and the daily sunspot number is 0.85 for 100 years of data. This methodology also permits us to study the variations in Enhanced, Active, and Quiet networks with time along with good accuracy for about a century, for the first time. Further, this procedure can be used to combine data from different observatories to make a long time series.


title: INetCEP: In-Network Complex Event Processing for Information-Centric Networking

abstract: Emerging network architectures like Information-centric Networking (ICN) offer simplicity in the data plane by addressing named data. Such flexibility opens up the possibility to move data processing inside network elements for high-performance computation, known as in-network processing. However, existing ICN architectures are limited in terms of data plane programmability due to the lack of (i) in-network processing and (ii) data plane programming abstractions. Such architectures can benefit from Complex Event Processing (CEP), an in-network processing paradigm to efficiently process data inside the data plane. Yet, it is extremely challenging to integrate CEP because the current communication model of ICN is limited to consumer-initiated interaction that comes with significant overhead in a number of requests to process continuous data streams. In contrast, a change to producer-initiated interaction, as favored by CEP, imposes severe limitations for request-reply interactions. In this paper, we propose an in-network CEP architecture, INetCEP that supports unified interaction patterns (consumer- and producer-initiated). In addition, we provide a CEP query language and facilitate CEP operations while increasing the range of applications that can be supported by ICN. We provide an open-source implementation and evaluation of INetCEP over an ICN architecture, Named Function Networking, and two applications: energy forecasting in smart homes and a disaster scenario.


title: Bayesian analysis of Juno/JIRAM's NIR observations of Europa

abstract: Juno spacecraft's spectrometer JIRAM recently observed the moon Europa in the 2-5 {\mu}m wavelength region. Here we present analysis of the average spectrum of a set of observations near 20{\deg}N and 40{\deg}W, focusing on the two forms of water-ice - amorphous and crystalline. We also take this as an opportunity to present a novel Bayesian spectral inversion framework for reflectance spectroscopy. We first validate this framework using simulated spectra of amorphous and crystalline ice mixtures and a laboratory spectrum of crystalline ice. We next analyze the JIRAM data and, through Bayesian model comparisons, find that a two-component intimately mixed model (TC-IM model) of amorphous and crystalline ice is strongly preferred (at 26{\sigma} confidence) over a two-component model of the same species but where their spectra are areally/linearly mixed. We also find that the TC-IM model is strongly preferred (at > 30{\sigma} confidence) over single-component models with only amorphous or crystalline ice, indicating the presence of both these phases of water ice in the data. For the highest SNR estimates of the JIRAM data, the TC-IM model solution corresponds to a mixture with a very large number density fraction (99.952 +/- 0.001 \%) of small (23.12 +/- 1.01 microns) amorphous ice grains, and a very small fraction (0.048 +/- 0.001 \%) of large (565.34 +/- 1.01 microns) crystalline ice grains. The overabundance of small amorphous ice grains we find is consistent with previous studies. The maximum-likelihood spectrum of the TC-IM model, however, is in tension with the data in the regions around 2.5 and 3.6 {\mu}m, and indicates the presence of non-ice components not currently included in our model, primarily due to the limited availability of cryogenic optical constants.


title: On-sky performance and recent results from the Subaru coronagraphic extreme adaptive optics system

abstract: We describe the current on-sky performance of the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument on the Subaru telescope on Maunakea, Hawaii. SCExAO is continuing to advance its AO performance, delivering H band Strehl ratios in excess of 0.9 for bright stars. We describe new advances with SCExAO's wavefront control that lead to a more stable corrected wavefront and diffraction-limited imaging in the optical, modifications to code that better handle read noise suppression within CHARIS, and tests of the spectrophotometric precision and accuracy within CHARIS. We outline steps in the CHARIS Data Processing Pipeline that output publication-grade data products. Finally, we note recent and upcoming science results, including the discovery of new directly-imaged systems and multiwavelength, deeper characterization of planet-forming disks, and upcoming technical advances that will improve SCExAO's sciencec capabilities.


title: Spiral Arm Pattern Motion in the SAO 206462 Protoplanetary Disk

abstract: Spiral arms have been observed in more than a dozen protoplanetary disks, yet the origin of nearly all systems is under debate. Multi-epoch monitoring of spiral arm morphology offers a dynamical way in distinguishing two leading arm formation mechanisms: companion-driven, and gravitational instability induction, since these mechanisms predict distinct motion patterns. By analyzing multi-epoch J-band observations of the SAO 206462 system using the SPHERE instrument on the Very Large Telescope (VLT) in 2015 and 2016, we measure the pattern motion for its two prominent spiral arms in polarized light. On one hand, if both arms are comoving, they can be driven by a planet at $86_{-13}^{+18}$ au on a circular orbit, with gravitational instability motion ruled out. On the other hand, they can be driven by two planets at $120_{-30}^{+30}$ au and $49_{-5}^{+6}$ au, offering a tentative evidence (3.0$\sigma$) that the two spirals are moving independently. The independent arm motion is possibly supported by our analysis of a re-reduction of archival observations using the NICMOS instrument onboard the Hubble Space Telescope (HST) in 1998 and 2005, yet artifacts including shadows can manifest spurious arm motion in HST observations. We expect future re-observations to better constrain the motion mechanism for the SAO 206462 spiral arms.


title: Polarization images of accretion flow around supermassive black holes: imprints of toroidal field structure

abstract: With unprecedented angular resolution, the Event Horizon Telescope (EHT) has opened a new era of black holes. We have previously calculated the expected polarization images of M87* with the EHT observations in mind. There, we demonstrated that circular polarization (CP) images, as well as the linear polarization (LP) maps, can convey quite useful information to us, such as the flow structure and magnetic field configuration around the black hole. In this paper, we make new predictions for the cases in which disk emission dominates over jet emission, bearing Sgr A* in mind. Here we set the proton-to-electron temperature ratio of the disk component to be Tp/Te $\sim$ 2 so as to suppress jet emission relative to emission from accretion flow. As a result, we obtain ring-like images and triple-forked images around the black hole for face-on and edge-on cases, respectively. We also find significant CP components in the images ($\gtrsim 10\%$ in fraction), both with positive and negative signs, amplified through the Faraday conversion, not sensitively depending on the inclination angles. Furthermore, we find a "separatrix" in the CP images, across which the sign of CP is reversed and on which the LP flux is brightest, that can be attributed to the helical magnetic field structure in the disk. These results indicate that future full polarization EHT images are a quite useful tracer of the magnetic field structure. We also discuss to what extent we will be able to extract information regarding the magnetic field configurations, under the scattering in the interstellar plasma, in future EHT polarimetric observations of Sgr A*.


title: Boundary topological entanglement entropy in two and three dimensions

abstract: The topological entanglement entropy is used to measure long-range quantum correlations in the ground state of topological phases. Here we obtain closed form expressions for topological entropy of (2+1)- and (3+1)-dimensional loop gas models, both in the bulk and at their boundaries, in terms of the data of their input fusion categories and algebra objects. Central to the formulation of our results are generalized $\mathcal{S}$-matrices. We conjecture a general property of these $\mathcal{S}$-matrices, with proofs provided in many special cases. This includes constructive proofs for categories up to rank 5.

categories: quant-ph cond-mat.str-el hep-th math-ph math.MP math.QA


title: Charting the Galactic acceleration field I. A search for stellar streams with Gaia DR2 and EDR3 with follow-up from ESPaDOnS and UVES

abstract: We present maps of the stellar streams detected in the Gaia Data Release 2 (DR2) and Early Data Release 3 (EDR3) catalogs using the STREAMFINDER algorithm. We also report the spectroscopic follow-up of the brighter DR2 stream members obtained with the high-resolution CFHT/ESPaDOnS and VLT/UVES spectrographs as well as with the medium-resolution NTT/EFOSC2 spectrograph. Two new stellar streams that do not have a clear progenitor are detected in DR2 (named Hrid and Gunnthra), and seven are detected in EDR3 (named Gaia-6 to Gaia-12). Several candidate streams are also identified. The software also finds very long tidal tails associated with the 15 globular clusters NGC 288, NGC 1261, NGC 1851, NGC 2298, NGC 2808, NGC 3201, M 68, $\omega$Cen, NGC 5466, Palomar 5, M 5, NGC 6101, M 92, NGC 6397 and NGC 7089. These stellar streams will be used in subsequent contributions in this series to chart the properties of the Galactic acceleration field on $\sim$ 100 pc to $\sim$ 100 kpc scales.


title: Cohomology of Lie Superalgebras: Forms, Integral Forms and Coset Superspaces

abstract: We study Chevalley-Eilenberg cohomology of physically relevant Lie superalgebras related to supersymmetric theories, providing explicit expressions for their cocycles in terms of their Maurer-Cartan forms. We then include integral forms in the picture by defining a notion of integral forms related to a Lie superalgebra. We develop a suitable generalization of Chevalley-Eilenberg cohomology extended to integral forms and we prove that it is isomorphic to the ordinary Chevalley-Eilenberg cohomology of the Lie superalgebra. Next we study equivariant Chevalley-Eilenberg cohomology for coset superspaces, which plays a crucial role in supergravity and superstring models. Again, we treat explicitly several examples, providing cocycles' expressions and revealing a characteristic infinite dimensional cohomology.


title: Building Tensor Networks for Holographic States

abstract: We discuss a one-parameter family of states in two-dimensional holographic conformal field theories which are constructed via the Euclidean path integral of an effective theory on a family of hyperbolic slices in the dual bulk geometry. The effective theory in question is the CFT flowed under a $T\overline{T}$ deformation, which "folds" the boundary CFT towards the bulk time-reflection symmetric slice. We propose that these novel Euclidean path integral states in the CFT can be interpreted as continuous tensor network (CTN) states. We argue that these CTN states satisfy a Ryu-Takayanagi-like minimal area upper bound on the entanglement entropies of boundary intervals, with the coefficient being equal to $\frac{1}{4G_N}$; the CTN corresponding to the bulk time-reflection symmetric slice saturates this bound. We also argue that the original state in the CFT can be written as a superposition of such CTN states, with the corresponding wavefunction being the bulk Hartle-Hawking wavefunction.


title: Entanglement Hamiltonian of Interacting Systems: Local Temperature Approximation and Beyond

abstract: We investigate the second quantization form of the entanglement Hamiltonian (EH) of various subregions for the ground-state of several interacting lattice fermions and spin models. The relation between the EH and the model Hamiltonian itself is an unsolved problem for the ground-state of generic local Hamiltonians. In this letter, we demonstrate that the EH is practically local and its dominant components are related to the terms present in the model Hamiltonian up to a smooth spatially varying temperature even for (a) discrete lattice systems, (b) systems with no emergent conformal or Lorentz symmetry, and (c) for subsystems with non-flat boundaries, up to relatively strong interactions. We show that the mentioned local temperature at a given point decays inversely proportional to its distance from the boundary between the subsystem and the environment. We find the subdominant terms in the EH as well and show that they are severely suppressed away from the boundaries of subsystem and are relatively small near them.

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


title: Kaluza-Klein fermion mass matrices from Exceptional Field Theory and ${\cal N} =1$ spectra

abstract: Using Exceptional Field Theory, we determine the infinite-dimensional mass matrices for the gravitino and spin-$1/2$ Kaluza-Klein perturbations above a class of anti-de Sitter solutions of M-theory and massive type IIA string theory with topologically-spherical internal spaces. We then use these mass matrices to compute the spectrum of Kaluza-Klein fermions about some solutions in this class with internal symmetry groups containing SU(3). Combining these results with previously known bosonic sectors of the spectra, we give the complete spectrum about some ${\cal N}=1$ and some non-supersymmetric solutions in this class. The complete spectra are shown to enjoy certain generic features.


title: The ACTIONFINDER: An unsupervised deep learning algorithm for calculating actions and the acceleration field from a set of orbit segments

abstract: We introduce the "ACTIONFINDER", a deep learning algorithm designed to transform a sample of phase-space measurements along orbits in a static potential into action and angle coordinates. The algorithm finds the mapping from positions and velocities to actions and angles in an unsupervised way, by using the fact that points along the same orbit have identical actions. Here we present the workings of the method, and test it on simple axisymmetric models, comparing the derived actions to those generated with the Torus Mapping technique. We show that it recovers the Torus actions for halo-type orbits in a realistic model of the Milky Way to $\sim 0.6$% accuracy with as few as 1024 input phase-space measurements. These actions are much better conserved along orbits than those estimated with the St\"ackel fudge. In our case, the reciprocal mapping from actions and angles to positions and velocities can also be learned. One of the advantages of the ACTIONFINDER is that it does not require the underlying potential to be known in advance, indeed it is designed to return the acceleration field. We expect the algorithm to be useful for analysing the properties of dynamical systems in numerical simulations. However, our ultimate goal with this effort will be to apply it to real stellar streams to recover the Galactic acceleration field in a way that is relatively agnostic about the underlying dark matter properties or the behavior of gravity.


title: Classification of Fermi-LAT sources with deep learning using energy and time spectra

abstract: Despite the growing number of gamma-ray sources detected by Fermi-LAT, about one third of the sources in each survey remains of uncertain type. We present a new deep neural network approach for the classification of unidentified or unassociated gamma-ray sources in the last release of the Fermi-LAT catalogue (4FGL-DR2) obtained with 10 years of data. In contrast to previous work, our method directly uses the measurements of the photon energy spectrum and time series as input for the classification, instead of specific, human-crafted features. Dense neural networks, and for the first time in the context of gamma-ray source classification recurrent neural networks, are studied in depth. We focus on the separation between extragalactic sources, i.e.\ Active Galactic Nuclei, and Galactic pulsars, and on the further classification of pulsars into young and millisecond pulsars. Our neural network architectures provide powerful classifiers, with a performance that is comparable to previous analyses based on human-crafted features. Our benchmark neural network predicts that of the sources of uncertain type in the 4FGL-DR2, 1050 are Active Galactic Nuclei and 78 are Galactic pulsars, with both classes following the expected sky distribution and the clustering in the variability-curvature plane. We investigate the problem of sample selection bias by testing our architectures against a cross-match test data set using an older catalogue, and propose a feature selection algorithm using autoencoders. Our list of high-confidence candidate sources labelled by the neural networks provides a set of targets for further multiwavelength observations addressed to identify their nature. The deep neural network architectures we develop can be easily extended to include specific features, as well as multiwavelength data on the source photon energy and time spectra coming from different instruments.


title: Thermal Instability and Multiphase Gas in the Simulated Interstellar Medium with Conduction, Viscosity and Magnetic Fields

abstract: Thermal instability (TI) plays a crucial role in the formation of multiphase structures and their dynamics in the Interstellar Medium (ISM) and is a leading theory for cold cloud creation in various astrophysical environments. In this paper we use two-dimensional (2D) simulations to investigate thermal instability under the influence of various initial conditions and physical processes. We experiment with Gaussian random field (GRF) density perturbations of different initial power spectra. We also enroll thermal conduction and physical viscosity in isotropic hydrodynamic and anisotropic magnetohydrodynamic (MHD) simulations. We find that the initial GRF spectral index $\alpha$ has a dramatic impact on the pure hydrodynamic development of thermal instability, influencing the size, number and motions of clouds. Cloud fragmentation happens due to two mechanisms: tearing and contraction rebound. In the runs with isotropic conduction and viscosity, the structures and dynamics of the clouds are dominated by evaporation and condensation flows in the non-linear regime, and the flow speed is regulated by viscosity. Cloud disruptions happen as a result of the Darrieus--Landau instability (DLI). Although at very late times, all individual clouds merge into one cold structure in all hydrodynamic runs. In the MHD case, the cloud structure is determined by both the initial perturbations and the initial magnetic field strength. In high $\beta$ runs, anisotropic conduction causes dense filaments to align with the local magnetic fields and the field direction can become reoriented. Strong magnetic fields suppress cross-field contraction and cold filaments can form along or perpendicular to the initial fields.


title: Speckle Imaging Characterization of Radial Velocity Exoplanet Systems

abstract: We conducted speckle imaging observations of 53 stellar systems that were members of long-term radial velocity (RV) monitoring campaigns and exhibited substantial accelerations indicative of planetary or stellar companions in wide orbits. Our observations were made with blue and red filters using the Differential Speckle Survey Instrument at Gemini-South and the NN-Explore Exoplanet Stellar Speckle Imager at the WIYN telescope. The speckle imaging identifies eight luminous companions within two arcseconds of the primary stars. In three of these systems (HD 1388, HD 87359, and HD 104304), the properties of the imaged companion are consistent with the RV measurements, suggesting that these companions may be associated with the primary and the cause of the RV variation. For all 53 stellar systems, we derive differential magnitude limits (i.e., contrast curves) from the imaging. We extend this analysis to include upper limits on companion mass in systems without imaging detections. In 25 systems, we rule out companions with mass greater than 0.2 $M_{\odot}$, suggesting that the observed RV signals are caused by late M dwarfs or substellar (potentially planetary) objects. On the other hand, the joint RV and imaging analysis almost entirely rules out planetary explanations of the RV signal for HD 19522 and suggests that the companion must have an angular separation below a few tenths of an arcsecond. This work highlights the importance of combined RV and imaging observations for characterizing the outer regions of nearby planetary systems.


title: A Geometric View of Closure Phases in Interferometry

abstract: Closure phase is the phase of a closed-loop product of correlations in a $\ge 3$-element interferometer array. Its invariance to element-based phase corruption makes it invaluable for interferometric applications that otherwise require high-accuracy phase calibration. However, its understanding has remained mainly mathematical and limited to the aperture plane (Fourier dual of image plane). Here, we lay the foundations for a geometrical insight. we show that closure phase and its invariance to element-based corruption and to translation are intricately related to the conserved properties (shape, orientation, and size, or SOS) of the principal triangle enclosed by the three fringes formed by a closed triad of array elements, which is referred herein as the "SOS conservation principle". When element-based amplitude calibration is not needed, as is typical in optical interferometry, the 3-element interference image formed from phase-uncalibrated correlations is a true and uncorrupted representation of the source object's morphology, except for a possible shift. Based on this SOS conservation principle, we present two geometric methods to measure the closure phase directly from a 3-element interference image (without requiring an aperture-plane view): (i) the closure phase is directly measurable from any one of the triangle's heights, and (ii) the squared closure phase is proportional to the product of the areas enclosed by the triad of array elements and the principal triangle in the aperture and image planes, respectively. We validate this geometric understanding across a wide range range of interferometric conditions using data from the Very Large Array and the Event Horizon Telescope. This geometric insight can be potentially valuable to other interferometric applications such as optical interferometry. These geometric relationships are generalised for an $N$-element interferometer.

categories: astro-ph.IM physics.class-ph physics.optics quant-ph


title: Theory of competing excitonic orders in insulating WTe$_2$ monolayers

abstract: We develop a theory of the excitonic phase recently proposed as the zero-field insulating state observed near charge neutrality in monolayer WTe$_2$. Using a Hartree-Fock approximation, we numerically identify two distinct gapped excitonic phases: a spin density wave state for weak non-zero interaction strength and spin spiral order at stronger interactions, separated by a narrow window of non-excitonic quantum spin Hall insulator. We introduce a simplified model capturing key features of the WTe$_2$ band structure, in which these phases appear as distinct valley ferromagnetic orders. We link the competition between the excitonic phases to the orbital structure of electronic wavefunctions at the Fermi surface and hence its proximity to the underlying gapped Dirac point in WTe$_2$. We briefly discuss collective modes of the two excitonic states, and comment on implications for experiments.


title: Tidal effects for spinning particles

abstract: Expanding on the recent derivation of tidal actions for scalar particles, we present here the action for a tidally deformed spin-$1/2$ particle. Focusing on operators containing two powers of the Weyl tensor, we combine the Hilbert series with an on-shell amplitude basis to construct the tidal action. With the tidal action in hand, we compute the leading-post-Minkowskian tidal contributions to the spin-1/2 -- spin-1/2 amplitude, arising at $\mathcal{O}(G^{2})$. Our amplitudes provide evidence that the observed long range spin-universality for the scattering of two point particles extends to the scattering of tidally deformed objects. From the scattering amplitude we find the conservative two-body Hamiltonian, linear and angular impulses, eikonal phase, spin kick, and aligned-spin scattering angle. We present analogous results in the electromagnetic case along the way.


title: Quantum Algorithms for Open Lattice Field Theory

abstract: Certain aspects of some unitary quantum systems are well-described by evolution via a non-Hermitian effective Hamiltonian, as in the Wigner-Weisskopf theory for spontaneous decay. Conversely, any non-Hermitian Hamiltonian evolution can be accommodated in a corresponding unitary system + environment model via a generalization of Wigner-Weisskopf theory. This demonstrates the physical relevance of novel features such as exceptional points in quantum dynamics, and opens up avenues for studying many body systems in the complex plane of coupling constants. In the case of lattice field theory, sparsity lends these channels the promise of efficient simulation on standardized quantum hardware. We thus consider quantum operations that correspond to Suzuki-Lee-Trotter approximation of lattice field theories undergoing non-Hermitian time evolution, with potential applicability to studies of spin or gauge models at finite chemical potential, with topological terms, to quantum phase transitions - a range of models with sign problems. We develop non-Hermitian quantum circuits and explore their promise on a benchmark, the quantum one-dimensional Ising model with complex longitudinal magnetic field, showing that observables can probe the Lee-Yang edge singularity. The development of attractors past critical points in the space of complex couplings indicates a potential for study on near-term noisy hardware.


title: ViP-DeepLab: Learning Visual Perception with Depth-aware Video Panoptic Segmentation

abstract: In this paper, we present ViP-DeepLab, a unified model attempting to tackle the long-standing and challenging inverse projection problem in vision, which we model as restoring the point clouds from perspective image sequences while providing each point with instance-level semantic interpretations. Solving this problem requires the vision models to predict the spatial location, semantic class, and temporally consistent instance label for each 3D point. ViP-DeepLab approaches it by jointly performing monocular depth estimation and video panoptic segmentation. We name this joint task as Depth-aware Video Panoptic Segmentation, and propose a new evaluation metric along with two derived datasets for it, which will be made available to the public. On the individual sub-tasks, ViP-DeepLab also achieves state-of-the-art results, outperforming previous methods by 5.1% VPQ on Cityscapes-VPS, ranking 1st on the KITTI monocular depth estimation benchmark, and 1st on KITTI MOTS pedestrian. The datasets and the evaluation codes are made publicly available.


title: Improved spectral gaps for random quantum circuits: large local dimensions and all-to-all interactions

abstract: Random quantum circuits are a central concept in quantum information theory with applications ranging from demonstrations of quantum computational advantage to descriptions of scrambling in strongly-interacting systems and black holes. The utility of random quantum circuits in these settings stems from their ability to rapidly generate quantum pseudo-randomness. In a seminal paper by Brand\~ao, Harrow, and Horodecki, it was proven that the $t$-th moment operator of local random quantum circuits on $n$ qudits with local dimension $q$ has a spectral gap of at least $\Omega(n^{-1}t^{-5-3.1/\log(q)})$, which implies that they are efficient constructions of approximate unitary designs. As a first result, we use Knabe bounds for the spectral gaps of frustration-free Hamiltonians to show that $1D$ random quantum circuits have a spectral gap scaling as $\Omega(n^{-1})$, provided that $t$ is small compared to the local dimension: $t^2\leq O(q)$. This implies a (nearly) linear scaling of the circuit depth in the design order $t$. Our second result is an unconditional spectral gap bounded below by $\Omega(n^{-1}\log^{-1}(n) t^{-\alpha(q)})$ for random quantum circuits with all-to-all interactions. This improves both the $n$ and $t$ scaling in design depth for the non-local model. We show this by proving a recursion relation for the spectral gaps involving an auxiliary random walk. Lastly, we solve the smallest non-trivial case exactly and combine with numerics and Knabe bounds to improve the constants involved in the spectral gap for small values of $t$.


title: Universal Fault-Tolerant Quantum Computing with Stabiliser Codes

abstract: The quantum logic gates used in the design of a quantum computer should be both universal, meaning arbitrary quantum computations can be performed, and fault-tolerant, meaning the gates keep errors from cascading out of control. A number of no-go theorems constrain the ways in which a set of fault-tolerant logic gates can be universal. These theorems are very restrictive, and conventional wisdom holds that a universal fault-tolerant logic gate set cannot be implemented natively, requiring us to use costly distillation procedures for quantum computation. Here, we present a general framework for universal fault-tolerant logic with stabiliser codes, together with a no-go theorem that reveals the very broad conditions constraining such gate sets. Our theorem applies to a wide range of stabiliser code families, including concatenated codes and conventional topological stabiliser codes such as the surface code. The broad applicability of our no-go theorem provides a new perspective on how the constraints on universal fault-tolerant gate sets can be overcome. In particular, we show how non-unitary implementations of logic gates provide a general approach to circumvent the no-go theorem, and we present a rich landscape of constructions for logic gate sets that are both universal and fault-tolerant. That is, rather than restricting what is possible, our no-go theorem provides a signpost to guide us to new, efficient architectures for fault-tolerant quantum computing.


title: Direct measurement of ferroelectric polarization in a tunable semimetal

abstract: Ferroelectricity, the electrostatic counterpart to ferromagnetism, has long been thought to be incompatible with metallicity due to screening of electric dipoles and external electric fields by itinerant charges. Recent measurements, however, demonstrated signatures of ferroelectric switching in the electrical conductance of bilayers and trilayers of WTe$_2$, a semimetallic transition metal dichalcogenide with broken inversion symmetry. An especially promising aspect of this system is that the density of electrons and holes can be continuously tuned by an external gate voltage. This degree of freedom enables investigation of the interplay between ferroelectricity and free carriers, a previously unexplored regime. Here, we employ capacitive sensing in dual-gated mesoscopic devices of bilayer WTe$_2$ to directly measure the spontaneous polarization in the metallic state and quantify the effect of free carriers on the polarization in the conduction and valence bands, separately. We compare our results to a low-energy model for the electronic bands and identify the layer-polarized states that contribute to transport and polarization simultaneously. Bilayer WTe$_2$ is thus shown to be a canonical example of a ferroelectric metal and an ideal platform for exploring polar ordering, ferroelectric transitions, and applications in the presence of free carriers.


title: Calibrating Core Overshooting Parameters With Two-dimensional Hydrodynamical Simulations

abstract: The extent of mixed regions around convective zones is one of the biggest uncertainties in stellar evolution. 1D overshooting descriptions introduce a free parameter ($f_{ov}$) that is in general not well constrained from observations. Especially in small central convective regions the value is highly uncertain due to its tight connection to the pressure scale height. Long-term multi-dimensional hydrodynamic simulations can be used to study the size of the overshooting region and the involved mixing processes. Here we show how one can calibrate an overshooting parameter by performing 2D Maestro simulations of Zero-Age-Main-Sequence stars ranging from $1.3$ to $3.5 M_\odot$. The simulations cover the convective cores of the stars and a large fraction of the surrounding radiative envelope. We follow the convective flow for at least 20 convective turnover times, while the longest simulation covers 430 turnover time scales. This allows us to study how the mixing as well as the convective boundary evolve with time, and how the resulting entrainment can be interpreted in terms of overshooting parameters. We find that increasing the overshooting parameter $f_{ov}$ beyond a certain value in the initial model of our simulations, changes the mixing behaviour completely. This result can be used to put limits on the overshooting parameter. We find $0.010 < f_{ov} < 0.017$ to be in good agreement with our simulations of a $3.5 M_\odot$ mass star. We also identify a diffusive mixing component due to internal gravity waves (IGW) that is active throughout the convectively stable layer, but likely overestimated in our simulations. Furthermore, applying our calibration method to simulations of less massive stars suggests a need for a mass-dependent overshooting description where the mixing in terms of the pressure scale height is reduced for small convective cores.


title: Dual diffeomorphisms and finite distance asymptotic symmetries in 3d gravity

abstract: We study the finite distance boundary symmetry current algebra of the most general first order theory of 3d gravity. We show that the space of quadratic generators contains diffeomorphisms but also a notion of dual diffeomorphisms, which together form either a double Witt or centreless BMS$_3$ algebra. The relationship with the usual asymptotic symmetry algebra relies on a duality between the null and angular directions, which is possible thanks to the existence of the dual diffeomorphisms.


title: Sparse data-driven quadrature rules via $\ell^p$-quasi-norm minimization

abstract: In this paper we show the use of the focal underdetermined system solver to recover sparse empirical quadrature rules for parametrized integrals from existing data, consisting of the values of given parametric functions sampled on a discrete set of points. This algorithm, originally proposed for image and signal reconstruction, relies on an approximated $\ell^p$-quasi-norm minimization. The choice of $0<p<1$ fits the nature of the constraints to which quadrature rules are subject, thus providing a more natural formulation for sparse quadrature recovery compared to the one based on $\ell^1$-norm minimization. We also extend an a priori error estimate available for the $\ell^1$-norm formulation by considering the error resulting from data compression. Finally, we present two numerical examples to illustrate some practical applications. The first concerns the fundamental solution of the linear 1D Schr\"odinger equation, the second example deals with the hyper-reduction of a partial differential equation modelling a nonlinear diffusion process in the framework of the reduced basis method. For both the examples we compare our method with the one based on $\ell^1$-norm minimization and the one relaying on the use of the non-negative least square method. Matlab codes related to the numerical examples and the algorithms described are provided.


title: Tensor network study of the $m=1/2$ magnetization plateau in the Shastry-Sutherland model at finite temperature

abstract: The two-dimensional infinite projected entangled pair states tensor network is evolved in imaginary time with the full update (FU) algorithm to simulate the Shastry-Sutherland model in a magnetic field at finite temperature directly in the thermodynamic limit. We focus on the phase transition into the $m=1/2$ magnetization plateau, which was observed in experiments on SrCu$_2$(BO$_3$)$_2$. For the largest simulated bond dimension, the early evolution in the high-temperature regime is simulated with the simple update (SU) scheme and then, as the correlation length increases, continued with the FU scheme towards the critical regime. We apply a small symmetry-breaking bias field and then extrapolate towards zero bias using a simple scaling theory in the bias field. The combined SU + FU scheme provides an accurate estimate of the critical temperature, even though the results could not be fully converged in the bond dimension in the vicinity of the transition. The critical temperature estimate is improved with a generalized scaling theory that combines two divergent length scales: One due to the bias, and the other due to the finite bond dimension. The obtained results are consistent with the transition being in the universality class of the two-dimensional classical Ising model. The estimated critical temperature is $3.5(2)$ K, which is well above the temperature $2.1$ K used in the experiments.


title: Optimising cost vs accuracy of decentralised analytics in fog computing environments

abstract: The exponential growth of devices and data at the edges of the Internet is rising scalability and privacy concerns on approaches based exclusively on remote cloud platforms. Data gravity, a fundamental concept in Fog Computing, points towards decentralisation of computation for data analysis, as a viable alternative to address those concerns. Decentralising AI tasks on several cooperative devices means identifying the optimal set of locations or Collection Points (CP for short) to use, in the continuum between full centralisation (i.e., all data on a single device) and full decentralisation (i.e., data on source locations). We propose an analytical framework able to find the optimal operating point in this continuum, linking the accuracy of the learning task with the corresponding network and computational cost for moving data and running the distributed training at the CPs. We show through simulations that the model accurately predicts the optimal trade-off, quite often an intermediate point between full centralisation and full decentralisation, showing also a significant cost saving w.r.t. both of them. Finally, the analytical model admits closed-form or numeric solutions, making it not only a performance evaluation instrument but also a design tool to configure a given distributed learning task optimally before its deployment.


title: On the Lattice of Conceptual Measurements

abstract: We present a novel approach for data set scaling based on scale-measures from formal concept analysis, i.e., continuous maps between closure systems, and derive a canonical representation. Moreover, we prove said scale-measures are lattice ordered with respect to the closure systems. This enables exploring the set of scale-measures through by the use of meet and join operations. Furthermore we show that the lattice of scale-measures is isomorphic to the lattice of sub-closure systems that arises from the original data. Finally, we provide another representation of scale-measures using propositional logic in terms of data set features. Our theoretical findings are discussed by means of examples.


title: Revisiting the Water Quality Sensor Placement Problem: Optimizing Network Observability and State Estimation Metrics

abstract: Real-time water quality (WQ) sensors in water distribution networks (WDN) have the potential to enable network-wide observability of water quality indicators, contamination event detection, and closed-loop feedback control of WQ dynamics. To that end, prior research has investigated a wide range of methods that guide the geographic placement of WQ sensors. These methods assign a metric for fixed sensor placement (SP) followed by \textit{metric-optimization} to obtain optimal SP. These metrics include minimizing intrusion detection time, minimizing the expected population and amount of contaminated water affected by an intrusion event. In contrast to the literature, the objective of this paper is to provide a computational method that considers the overlooked metric of state estimation and network-wide observability of the WQ dynamics. This metric finds the optimal WQ sensor placement that minimizes the state estimation error via the Kalman filter for noisy WQ dynamics -- a metric that quantifies WDN observability. To that end, the state-space dynamics of WQ states for an entire WDN are given and the observability-driven sensor placement algorithm is presented. The algorithm takes into account the time-varying nature of WQ dynamics due to changes in the hydraulic profile -- a collection of hydraulic states including heads (pressures) at nodes and flow rates in links which are caused by a demand profile over a certain period of time. Thorough case studies are given, highlighting key findings, observations, and recommendations for WDN operators. Github codes are included for reproducibility.


title: Hard and Soft EM in Bayesian Network Learning from Incomplete Data

abstract: Incomplete data are a common feature in many domains, from clinical trials to industrial applications. Bayesian networks (BNs) are often used in these domains because of their graphical and causal interpretations. BN parameter learning from incomplete data is usually implemented with the Expectation-Maximisation algorithm (EM), which computes the relevant sufficient statistics ("soft EM") using belief propagation. Similarly, the Structural Expectation-Maximisation algorithm (Structural EM) learns the network structure of the BN from those sufficient statistics using algorithms designed for complete data. However, practical implementations of parameter and structure learning often impute missing data ("hard EM") to compute sufficient statistics instead of using belief propagation, for both ease of implementation and computational speed. In this paper, we investigate the question: what is the impact of using imputation instead of belief propagation on the quality of the resulting BNs? From a simulation study using synthetic data and reference BNs, we find that it is possible to recommend one approach over the other in several scenarios based on the characteristics of the data. We then use this information to build a simple decision tree to guide practitioners in choosing the EM algorithm best suited to their problem.


title: MLComp: A Methodology for Machine Learning-based Performance Estimation and Adaptive Selection of Pareto-Optimal Compiler Optimization Sequences

abstract: Embedded systems have proliferated in various consumer and industrial applications with the evolution of Cyber-Physical Systems and the Internet of Things. These systems are subjected to stringent constraints so that embedded software must be optimized for multiple objectives simultaneously, namely reduced energy consumption, execution time, and code size. Compilers offer optimization phases to improve these metrics. However, proper selection and ordering of them depends on multiple factors and typically requires expert knowledge. State-of-the-art optimizers facilitate different platforms and applications case by case, and they are limited by optimizing one metric at a time, as well as requiring a time-consuming adaptation for different targets through dynamic profiling. To address these problems, we propose the novel MLComp methodology, in which optimization phases are sequenced by a Reinforcement Learning-based policy. Training of the policy is supported by Machine Learning-based analytical models for quick performance estimation, thereby drastically reducing the time spent for dynamic profiling. In our framework, different Machine Learning models are automatically tested to choose the best-fitting one. The trained Performance Estimator model is leveraged to efficiently devise Reinforcement Learning-based multi-objective policies for creating quasi-optimal phase sequences. Compared to state-of-the-art estimation models, our Performance Estimator model achieves lower relative error (<2%) with up to 50x faster training time over multiple platforms and application domains. Our Phase Selection Policy improves execution time and energy consumption of a given code by up to 12% and 6%, respectively. The Performance Estimator and the Phase Selection Policy can be trained efficiently for any target platform and application domain.


title: Measuring the Sun's velocity using Gaia EDR3 observations of Stellar Streams

abstract: We measure the Sun's velocity with respect to the Galactic halo using Gaia Early Data Release 3 (EDR3) observations of stellar streams. Our method relies on the fact that, in low-mass streams, the proper motion of stars should be directed along the stream structure in a non-rotating rest frame of the Galaxy, but the observed deviation arises due to the Sun's own reflex motion. This principle allows us to implement a simple geometrical procedure, which we use to analyse 17 streams over a $\sim 3-30$ kpc range. Our constraint on the Sun's motion is independent of any Galactic potential model, and it is also uncorrelated with the Sun's galactocentric distance. We infer the Sun's velocity as $V_{R,\odot}=8.88^{+1.20}_{-1.22}\,\rm{kms^{-1}}$ (radially towards the Galactic centre), $V_{\phi,\odot}=241.91^{+1.61}_{-1.73}\,\rm{kms^{-1}}$ (in the direction of Galactic rotation) and $V_{z,\odot}=3.08^{+1.06}_{-1.10}\,\rm{kms^{-1}}$ (vertically upwards), in global agreement with past measurements through other techniques; although we do note a small but significant difference in the $V_{z,\odot}$ component. Some of these parameters show significant correlation and we provide our MCMC output so it can be used by the reader as an input to future works. The comparison between our Sun's velocity inference and previous results, using other reference frames, indicates that the inner Galaxy is not moving with respect to the inertial frame defined by the halo streams.


title: Off-shell Higgs Couplings in $H^*\to ZZ\to \ell\ell\nu\nu$

abstract: We explore the new physics reach for the off-shell Higgs boson measurement in the ${pp \to H^* \rightarrow Z(\ell^{+}\ell^{-})Z(\nu\bar{\nu})}$ channel at the high-luminosity LHC. The new physics sensitivity is parametrized in terms of the Higgs boson width, effective field theory framework, and a non-local Higgs-top coupling form factor. Adopting Machine-learning techniques, we demonstrate that the combination of a large signal rate and a precise phenomenological probe for the process energy scale, due to the transverse $ZZ$ mass, leads to significant sensitivities beyond the existing results in the literature for the new physics scenarios considered.


title: MetaInfoNet: Learning Task-Guided Information for Sample Reweighting

abstract: Deep neural networks have been shown to easily overfit to biased training data with label noise or class imbalance. Meta-learning algorithms are commonly designed to alleviate this issue in the form of sample reweighting, by learning a meta weighting network that takes training losses as inputs to generate sample weights. In this paper, we advocate that choosing proper inputs for the meta weighting network is crucial for desired sample weights in a specific task, while training loss is not always the correct answer. In view of this, we propose a novel meta-learning algorithm, MetaInfoNet, which automatically learns effective representations as inputs for the meta weighting network by emphasizing task-related information with an information bottleneck strategy. Extensive experimental results on benchmark datasets with label noise or class imbalance validate that MetaInfoNet is superior to many state-of-the-art methods.