title: Why Normalizing Flows Fail to Detect Out-of-Distribution Data

abstract: Detecting out-of-distribution (OOD) data is crucial for robust machine learning systems. Normalizing flows are flexible deep generative models that often surprisingly fail to distinguish between in- and out-of-distribution data: a flow trained on pictures of clothing assigns higher likelihood to handwritten digits. We investigate why normalizing flows perform poorly for OOD detection. We demonstrate that flows learn local pixel correlations and generic image-to-latent-space transformations which are not specific to the target image dataset. We show that by modifying the architecture of flow coupling layers we can bias the flow towards learning the semantic structure of the target data, improving OOD detection. Our investigation reveals that properties that enable flows to generate high-fidelity images can have a detrimental effect on OOD detection.


title: Derivative couplings in gravitational production in the early universe

abstract: Gravitational particle production in the early universe is due to the coupling of matter fields to curvature. This coupling may include derivative terms that modify the kinetic term. The most general first order action contains derivative couplings to the curvature scalar and to the traceless Ricci tensor, which can be dominant in the case of (pseudo-)Nambu-Goldstone bosons or disformal scalars, such as branons. In the presence of these derivative couplings, the density of produced particles for the adiabatic regime in the de Sitter phase (which mimics inflation) is constant in time and decays with the inverse effective mass (which in turn depends on the coupling to the curvature scalar). In the reheating phase following inflation, the presence of derivative couplings to the background curvature modifies in a nontrivial way the gravitational production even in the perturbative regime. We also show that the two couplings -- to the curvature scalar and to the traceless Ricci tensor -- are drastically different, specially for large masses. In this regime, the production becomes highly sensitive to the former coupling while it becomes independent of the latter.


title: Visibility Guided NMS: Efficient Boosting of Amodal Object Detection in Crowded Traffic Scenes

abstract: Object detection is an important task in environment perception for autonomous driving. Modern 2D object detection frameworks such as Yolo, SSD or Faster R-CNN predict multiple bounding boxes per object that are refined using Non-Maximum-Suppression (NMS) to suppress all but one bounding box. While object detection itself is fully end-to-end learnable and does not require any manual parameter selection, standard NMS is parametrized by an overlap threshold that has to be chosen by hand. In practice, this often leads to an inability of standard NMS strategies to distinguish different objects in crowded scenes in the presence of high mutual occlusion, e.g. for parked cars or crowds of pedestrians. Our novel Visibility Guided NMS (vg-NMS) leverages both pixel-based as well as amodal object detection paradigms and improves the detection performance especially for highly occluded objects with little computational overhead. We evaluate vg-NMS using KITTI, VIPER as well as the Synscapes dataset and show that it outperforms current state-of-the-art NMS.


title: A Note on Nesterov's Accelerated Method in Nonconvex Optimization: a Weak Estimate Sequence Approach

abstract: We present a variant of accelerated gradient descent algorithms, adapted from Nesterov's optimal first-order methods, for weakly-quasi-convex and weakly-quasi-strongly-convex functions. We show that by tweaking the so-called estimate sequence method, the derived algorithm achieves optimal convergence rate for weakly-quasi-convex and weakly-quasi-strongly-convex in terms of oracle complexity. In particular, for a weakly-quasi-convex function with Lipschitz continuous gradient, we require $O(\frac{1}{\sqrt{\varepsilon}})$ iterations to acquire an $\varepsilon$-solution; for weakly-quasi-strongly-convex functions, the iteration complexity is $O\left( \ln\left(\frac{1}{\varepsilon}\right) \right)$. Furthermore, we discuss the implications of these algorithms for linear quadratic optimal control problem.


title: Second Order Unconditional Positive Preserving Schemes for Non-equilibrium Reactive Flows with Mass and Mole Balance

abstract: In this study, a family of second order process based modified Patankar Runge-Kutta schemes is proposed with both the mass and mole maintained in balance while preserving the positivity of density and pressure with the time step determined by convection terms. The accuracy analysis is conducted to derive the sufficient and necessary conditions for the Runge-Kutta and Patankar coefficients. Coupling with the finite volume method, the proposed schemes are extended to Euler equations with non-equilibrium reacting source terms. Benchmark tests are given to prove the prior order of accuracy and validate the positive-preserving property for both density and pressure.


title: Optimization and Generalization Analysis of Transduction through Gradient Boosting and Application to Multi-scale Graph Neural Networks

abstract: It is known that the current graph neural networks (GNNs) are difficult to make themselves deep due to the problem known as over-smoothing. Multi-scale GNNs are a promising approach for mitigating the over-smoothing problem. However, there is little explanation of why it works empirically from the viewpoint of learning theory. In this study, we derive the optimization and generalization guarantees of transductive learning algorithms that include multi-scale GNNs. Using the boosting theory, we prove the convergence of the training error under weak learning-type conditions. By combining it with generalization gap bounds in terms of transductive Rademacher complexity, we show that a test error bound of a specific type of multi-scale GNNs that decreases corresponding to the number of node aggregations under some conditions. Our results offer theoretical explanations for the effectiveness of the multi-scale structure against the over-smoothing problem. We apply boosting algorithms to the training of multi-scale GNNs for real-world node prediction tasks. We confirm that its performance is comparable to existing GNNs, and the practical behaviors are consistent with theoretical observations. Code is available at https://github.com/delta2323/GB-GNN.


title: VLA imaging of the XMM-LSS / VIDEO deep field at 1-2 GHz

abstract: Modern radio telescopes are routinely reaching depths where normal starforming galaxies are the dominant observed population. Realising the potential of radio as a tracer of star formation and black hole activity over cosmic time involves achieving such depths over representative volumes, with radio forming part of a larger multiwavelength campaign. In pursuit of this we used the Karl G. Jansky Very Large Array (VLA) to image $\sim$5 deg$^{2}$ of the VIDEO/XMM-LSS extragalactic deep field at 1--2 GHz. We achieve a median depth of 16 $\mu$Jy beam$^{-1}$ with an angular resolution of 4.5\arcsec. Comparisons with existing radio observations of XMM-LSS showcase the improved survey speed of the upgraded VLA: we cover 2.5 times the area and increase the depth by $\sim$20\% in 40\% of the time. Direction-dependent calibration and wide-field imaging were required to suppress the error patterns from off-axis sources of even modest brightness. We derive a catalogue containing 5,762 sources from the final mosaic. Sub-band imaging provides in-band spectral indices for 3,458 (60\%) sources, with the average spectrum becoming flatter than the canonical synchrotron slope below 1 mJy. Positional and flux-density accuracy of the observations, and the differential source counts are in excellent agreement with those of existing measurements. A public release of the images and catalogue accompanies this article.


title: Universal programmable on-chip metasurface building blocks for arbitrary high-order mode manipulation

abstract: On-chip mode-division multiplexing (MDM) has been emerging as a promising technology to further enhance the link capacity and bandwidth of data communications with multiple mode channels. Both mode converters and mode exchangers are indispensable fundamental components for flexible mode operations. While several configurations have been developed previously, it is still very challenging to efficiently manipulate arbitrary high-order modes in a versatile way to reduce the R&D and prototyping costs. Here we initiate a breakthrough with a simple yet universal generic mode operator building block concept utilizing metasurface structures. The programmable arbitrary high-order mode operators can realize mode conversion and exchange simultaneously with competitive and uniform performance, high stability and compact footprints, offering a quintessential step change for on-chip multimode optical interconnections.


title: tmleCommunity: A R Package Implementing Target Maximum Likelihood Estimation for Community-level Data

abstract: Over the past years, many applications aim to assess the causal effect of treatments assigned at the community level, while data are still collected at the individual level among individuals of the community. In many cases, one wants to evaluate the effect of a stochastic intervention on the community, where all communities in the target population receive probabilistically assigned treatments based on a known specified mechanism (e.g., implementing a community-level intervention policy that target stochastic changes in the behavior of a target population of communities). The tmleCommunity package is recently developed to implement targeted minimum loss-based estimation (TMLE) of the effect of community-level intervention(s) at a single time point on an individual-based outcome of interest, including the average causal effect. Implementations of the inverse-probability-of-treatment-weighting (IPTW) and the G-computation formula (GCOMP) are also available. The package supports multivariate arbitrary (i.e., static, dynamic or stochastic) interventions with a binary or continuous outcome. Besides, it allows user-specified data-adaptive machine learning algorithms through SuperLearner, sl3 and h2oEnsemble packages. The usage of the tmleCommunity package, along with a few examples, will be described in this paper.


title: Now that I can see, I can improve: Enabling data-driven finetuning of CNNs on the edge

abstract: In today's world, a vast amount of data is being generated by edge devices that can be used as valuable training data to improve the performance of machine learning algorithms in terms of the achieved accuracy or to reduce the compute requirements of the model. However, due to user data privacy concerns as well as storage and communication bandwidth limitations, this data cannot be moved from the device to the data centre for further improvement of the model and subsequent deployment. As such there is a need for increased edge intelligence, where the deployed models can be fine-tuned on the edge, leading to improved accuracy and/or reducing the model's workload as well as its memory and power footprint. In the case of Convolutional Neural Networks (CNNs), both the weights of the network as well as its topology can be tuned to adapt to the data that it processes. This paper provides a first step towards enabling CNN finetuning on an edge device based on structured pruning. It explores the performance gains and costs of doing so and presents an extensible open-source framework that allows the deployment of such approaches on a wide range of network architectures and devices. The results show that on average, data-aware pruning with retraining can provide 10.2pp increased accuracy over a wide range of subsets, networks and pruning levels with a maximum improvement of 42.0pp over pruning and retraining in a manner agnostic to the data being processed by the network.


title: Pipeline PSRO: A Scalable Approach for Finding Approximate Nash Equilibria in Large Games

abstract: Finding approximate Nash equilibria in zero-sum imperfect-information games is challenging when the number of information states is large. Policy Space Response Oracles (PSRO) is a deep reinforcement learning algorithm grounded in game theory that is guaranteed to converge to an approximate Nash equilibrium. However, PSRO requires training a reinforcement learning policy at each iteration, making it too slow for large games. We show through counterexamples and experiments that DCH and Rectified PSRO, two existing approaches to scaling up PSRO, fail to converge even in small games. We introduce Pipeline PSRO (P2SRO), the first scalable general method for finding approximate Nash equilibria in large zero-sum imperfect-information games. P2SRO is able to parallelize PSRO with convergence guarantees by maintaining a hierarchical pipeline of reinforcement learning workers, each training against the policies generated by lower levels in the hierarchy. We show that unlike existing methods, P2SRO converges to an approximate Nash equilibrium, and does so faster as the number of parallel workers increases, across a variety of imperfect information games. We also introduce an open-source environment for Barrage Stratego, a variant of Stratego with an approximate game tree complexity of $10^{50}$. P2SRO is able to achieve state-of-the-art performance on Barrage Stratego and beats all existing bots. Experiment code is available athttps://github.com/JBLanier/pipeline-psro.


title: Wavefront dislocations reveal the topology of quasi-1D photonic insulators

abstract: Phase singularities appear ubiquitously in wavefields, regardless of the wave equation. Such topological defects can lead to wavefront dislocations, as observed in a humongous number of classical wave experiments. Phase singularities of wave functions are also at the heart of the topological classification of the gapped phases of matter. Despite identical singular features, topological insulators and topological defects in waves remain two distinct fields. Realising 1D microwave insulators, we experimentally observe a wavefront dislocation - a 2D phase singularity - in the local density of states when the systems undergo a topological phase transition. We show theoretically that the change in the number of interference fringes at the transition reveals the topological index that characterises the band topology in the insulator.


title: Correspondence Modules and Persistence Sheaves: A Unifying Perspective on One-Parameter Persistent Homology

abstract: We develop a unifying framework for the treatment of various persistent homology architectures using the notion of correspondence modules. In this formulation, morphisms between vector spaces are given by partial linear relations, as opposed to linear mappings. In the one-dimensional case, among other things, this allows us to: (i) treat persistence modules and zigzag modules as algebraic objects of the same type; (ii) give a categorical formulation of zigzag structures over a continuous parameter; and (iii) construct barcodes associated with spaces and mappings that are richer in geometric information. A structural analysis of one-parameter persistence is carried out at the level of sections of correspondence modules that yield sheaf-like structures, termed persistence sheaves. Under some tameness hypotheses, we prove interval decomposition theorems for persistence sheaves and correspondence modules, as well as an isometry theorem for persistence diagrams obtained from interval decompositions. Applications include: (a) a Mayer-Vietoris sequence that relates the persistent homology of sublevelset filtrations and superlevelset filtrations to the levelset homology module of a real-valued function and (b) the construction of slices of 2-parameter persistence modules along negatively sloped lines.


title: Learning Diverse and Discriminative Representations via the Principle of Maximal Coding Rate Reduction

abstract: To learn intrinsic low-dimensional structures from high-dimensional data that most discriminate between classes, we propose the principle of Maximal Coding Rate Reduction ($\text{MCR}^2$), an information-theoretic measure that maximizes the coding rate difference between the whole dataset and the sum of each individual class. We clarify its relationships with most existing frameworks such as cross-entropy, information bottleneck, information gain, contractive and contrastive learning, and provide theoretical guarantees for learning diverse and discriminative features. The coding rate can be accurately computed from finite samples of degenerate subspace-like distributions and can learn intrinsic representations in supervised, self-supervised, and unsupervised settings in a unified manner. Empirically, the representations learned using this principle alone are significantly more robust to label corruptions in classification than those using cross-entropy, and can lead to state-of-the-art results in clustering mixed data from self-learned invariant features.


title: MCR-TRGB: A Multiwavelength-Covariant, Robust Tip of the Red Giant Branch Measurement Method

abstract: We present a new method to measure colors and magnitudes of the tip of the red giant branch in multiple bandpasses simultaneously by fitting an n-dimensional Gaussian to photometry of candidate tip stars. We demonstrate that this method has several advantages over traditional edge detection, particularly in regimes where the TRGB magnitude is strongly color-dependent, as is the case in the near-infrared. We apply this method to a re-reduction of a set of optical and near-IR HST data originally presented in Dalcanton et al. (2012). The re-reduction takes advantage of the increased depth and accuracy in the NIR photometry enabled by simultaneous reduction with higher resolution optical data in crowded fields (Williams et al. 2014). We compare three possible absolute calibrations of the resulting apparent TRGB measurements, one adopting the same distance moduli as in Dalcanton et al. (2012), and two based on predicted TRGB absolute magnitudes from two widely-used, modern sets of model isochrones. We find systematic offsets among the model absolute calibrations at the ~0.1 mag level, in line with previous investigations. The models also have difficulty reproducing the optical-NIR color-magnitude behavior of our measurements, making these observations a useful benchmark for future improvements.


title: Blocking defector invasion by focusing on the most successful partner

abstract: According to the standard protocol of spatial public goods game, a cooperator player invests not only into his own game but also into the games organized by neighboring partners. In this work, we relax this assumption by allowing cooperators to decide which neighboring group to prefer instead of supporting them uniformly. In particular, we assume that they select their most successful neighbor and focus external investments exclusively into the related group. We show that this very simple alteration of the dynamical rule results in a surprisingly positive evolutionary outcome -- cooperators prevail even in harsh environment represented by small values of the synergy factor in the game. The microscopic mechanism behind the reported success of the cooperator strategy can be explained by a blocking mechanism which affects the propagations of competing strategies in a biased way. Our results, which remain intact by using different interaction topologies, reveal that it could be beneficial to concentrate individual efforts to reach a higher global wellbeing.


title: A new approach to observational cosmology using the scattering transform

abstract: Parameter estimation with non-Gaussian stochastic fields is a common challenge in astrophysics and cosmology. In this paper, we advocate performing this task using the scattering transform, a statistical tool sharing ideas with convolutional neural networks (CNNs) but requiring no training nor tuning. It generates a compact set of coefficients, which can be used as robust summary statistics for non-Gaussian information. It is especially suited for fields presenting localized structures and hierarchical clustering, such as the cosmological density field. To demonstrate its power, we apply this estimator to a cosmological parameter inference problem in the context of weak lensing. On simulated convergence maps with realistic noise, the scattering transform outperforms classic estimators and is on a par with state-of-the-art CNN. It retains the advantages of traditional statistical descriptors, has provable stability properties, allows to check for systematics, and importantly, the scattering coefficients are interpretable. It is a powerful and attractive estimator for observational cosmology and the study of physical fields in general.


title: Recent Performance Studies of the GEM-based TPC Readout (DESY Module)

abstract: For the International Large Detector (ILD) at the planned International Linear Collider (ILC) a Time Projection Chamber (TPC) is foreseen as the main tracking detector. To achieve the required point resolution, Micro-Pattern Gaseous Detectors (MPGD) will be used in the amplification stage. A readout module using a stack of three Gas Electron Multipliers (GEM) for gas amplification was developed at DESY and tested at the DESY II Test Beam Facility. After introducing the readout module and the infrastructure at the test beam facility, the performance related to single point and double-hit resolution of three of these modules is presented. This is followed by results on the particle identification capabilities of the system, using the specific energy loss dE/dx, and simulation studies, aimed to investigate and quantify the impact of high granularity on dE/dx resolution. In addition, a new and improved TPC field cage and the LYCORIS Large-Area Silicon-Strip Telescope for the test beam are described. The LYCORIS beam telescope is foreseen to provide a precise reference of the particle trajectory to validate the momentum resolution measured with a large TPC prototype. For this purpose, it is being installed and tested at the test beam facility within the so-called PCMAG (Persistent Current Magnet).


title: Counting points of bounded height in monoid orbits

abstract: Given a set of endomorphisms on $\mathbb{P}^N$, we establish an upper bound on the number of points of bounded height in the associated monoid orbits. Moreover, we give a more refined estimate with an associated lower bound when the monoid is free. Finally, we show that most sets of rational functions in one variable satisfy these more refined bounds.


title: Intra-Processing Methods for Debiasing Neural Networks

abstract: As deep learning models become tasked with more and more decisions that impact human lives, such as criminal recidivism, loan repayment, and face recognition for law enforcement, bias is becoming a growing concern. Debiasing algorithms are typically split into three paradigms: pre-processing, in-processing, and post-processing. However, in computer vision or natural language applications, it is common to start with a large generic model and then fine-tune to a specific use-case. Pre- or in-processing methods would require retraining the entire model from scratch, while post-processing methods only have black-box access to the model, so they do not leverage the weights of the trained model. Creating debiasing algorithms specifically for this fine-tuning use-case has largely been neglected. In this work, we initiate the study of a new paradigm in debiasing research, intra-processing, which sits between in-processing and post-processing methods. Intra-processing methods are designed specifically to debias large models which have been trained on a generic dataset and fine-tuned on a more specific task. We show how to repurpose existing in-processing methods for this use-case, and we also propose three baseline algorithms: random perturbation, layerwise optimization, and adversarial fine-tuning. All of our techniques can be used for all popular group fairness measures such as equalized odds or statistical parity difference. We evaluate these methods across three popular datasets from the AIF360 toolkit, as well as on the CelebA faces dataset. Our code is available at https://github.com/abacusai/intraprocessing_debiasing.


title: Spectral DiffuserCam: lensless snapshot hyperspectral imaging with a spectral filter array

abstract: Hyperspectral imaging is useful for applications ranging from medical diagnostics to agricultural crop monitoring; however, traditional scanning hyperspectral imagers are prohibitively slow and expensive for widespread adoption. Snapshot techniques exist but are often confined to bulky benchtop setups or have low spatio-spectral resolution. In this paper, we propose a novel, compact, and inexpensive computational camera for snapshot hyperspectral imaging. Our system consists of a tiled spectral filter array placed directly on the image sensor and a diffuser placed close to the sensor. Each point in the world maps to a unique pseudorandom pattern on the spectral filter array, which encodes multiplexed spatio-spectral information. By solving a sparsity-constrained inverse problem, we recover the hyperspectral volume with sub-super-pixel resolution. Our hyperspectral imaging framework is flexible and can be designed with contiguous or non-contiguous spectral filters that can be chosen for a given application. We provide theory for system design, demonstrate a prototype device, and present experimental results with high spatio-spectral resolution.


title: Characterizing the radio continuum nature of sources in the massive star-forming region W75N (B)

abstract: The massive star-forming region W75N~(B) is thought to host a cluster of massive protostars (VLA~1, VLA~2, and VLA~3) undergoing different evolutionary stages. In this work, we present radio continuum data with the highest sensitivity and angular resolution obtained to date in this region, using the VLA-A and covering a wide range of frequencies (4-48~GHz), which allowed us to study the morphology and the nature of the emission of the different radio continuum sources. We also performed complementary studies with multi-epoch VLA data and ALMA archive data at 1.3 mm wavelength. We find that VLA~1 is driving a thermal radio jet at scales of $\approx$0.1 arcsec ($\approx$130 au), but also shows signs of an incipient hyper-compact HII region at scales of $\lesssim$ 1 arcsec ($\lesssim$ 1300~au). VLA~3 is also driving a thermal radio jet at scales of a few tenths of arcsec (few hundred of au). We conclude that this jet is shock-exciting the radio continuum sources Bc and VLA~4 (obscured HH objects), which show proper motions moving outward from VLA~3 at velocities of $\approx$112--118~km/s. We have also detected three new weak radio continuum sources, two of them associated with millimeter continuum cores observed with ALMA, suggesting that these two sources are also embedded YSOs in this massive star-forming region.


title: Ergodic cocycles of IDPFT systems and nonsingular Gaussian actions

abstract: It is proved that each Gaussian cocycle over a mildly mixing Gaussian transformation is either a Gaussian coboundary or sharply weak mixing. The class of nonsingular infinite direct products $T$ of transformations $T_n$, $n\in\Bbb N$, of finite type (IDPFT) is studied. It is shown that if $T_n$ is mildly mixing, $n\in\Bbb N$, the sequence of the Radon-Nikodym derivatives of $T_n$ is asymptotically translation quasi-invariant and $T$ is conservative then the Maharam extension of $T$ is sharply weak mixing. This techniques provides a new approach to the nonsingular Gaussian transformations studied recently by Arano, Isono and Marrakchi.


title: PrivyTRAC: Privacy and Security Preserving Contact Tracing System

abstract: Smartphone location-based methods have been proposed and implemented as an effective alternative to traditional labor intensive contact tracing methods. However, there are serious privacy and security concerns that may impede wide-spread adoption in many societies. Furthermore, these methods rely solely on proximity to patients, based on Bluetooth or GPS signal for example, ignoring lingering effects of virus, including COVID-19, present in the environment. This results in inaccurate risk assessment and incomplete contact tracing. A new system concept, called PrivyTRAC, preserves user privacy, increases security and improves accuracy of smartphone contact tracing. PrivyTRAC enhances users' and patients' privacy by letting users conduct self-evaluation based on the risk maps download to their smartphones. No user information is transmitted to external locations or devices, and no personally identifiable patient information is embedded in the risk maps as they are processed anonymized and aggregated locations of confirmed patients. The risk maps consider both spatial proximity and temporal effects to improve the accuracy of the infection risk estimation. Experiments conducted in the paper illustrate improvement of PrivyTRAC over proximity based methods in terms of true and false positives. An approach to further improve infection risk estimation by incorporating both positive and negative local test results from contacts of confirmed cases is also described.


title: Strongly local p-norm-cut algorithms for semi-supervised learning and local graph clustering

abstract: Graph based semi-supervised learning is the problem of learning a labeling function for the graph nodes given a few example nodes, often called seeds, usually under the assumption that the graph's edges indicate similarity of labels. This is closely related to the local graph clustering or community detection problem of finding a cluster or community of nodes around a given seed. For this problem, we propose a novel generalization of random walk, diffusion, or smooth function methods in the literature to a convex p-norm cut function. The need for our p-norm methods is that, in our study of existing methods, we find those principled methods based on eigenvector, spectral, random walk, or linear system often have difficulty capturing the correct boundary of a target label or target cluster. In contrast, 1-norm or maxflow-mincut based methods capture the boundary, but cannot grow from small seed set; hybrid procedures that use both have many hard to set parameters. In this paper, we propose a generalization of the objective function behind these methods involving p-norms. To solve the p-norm cut problem we give a strongly local algorithm -- one whose runtime depends on the size of the output rather than the size of the graph. Our method can be thought as a nonlinear generalization of the Anderson-Chung-Lang push procedure to approximate a personalized PageRank vector efficiently. Our procedure is general and can solve other types of nonlinear objective functions, such as p-norm variants of Huber losses. We provide a theoretical analysis of finding planted target clusters with our method and show that the p-norm cut functions improve on the standard Cheeger inequalities for random walk and spectral methods. Finally, we demonstrate the speed and accuracy of our new method in synthetic and real world datasets. Our code is available at http://github.com/MengLiuPurdue/SLQ.


title: Cross-temporal forecast reconciliation: Optimal combination method and heuristic alternatives

abstract: Forecast reconciliation is a post-forecasting process aimed to improve the quality of the base forecasts for a system of hierarchical/grouped time series (Hyndman et al., 2011). Contemporaneous (cross-sectional) and temporal hierarchies have been considered in the literature, but - except for Kourentzes and Athanasopoulos (2019) - generally these two features have not been fully considered together. Adopting a notation able to simultaneously deal with both forecast reconciliation dimensions, the paper shows two new results: (i) an iterative cross-temporal forecast reconciliation procedure which extends, and overcomes some weaknesses of, the two-step procedure by Kourentzes and Athanasopoulos (2019), and (ii) the closed-form expression of the optimal (in least squares sense) point forecasts which fulfill both contemporaneous and temporal constraints. The feasibility of the proposed procedures, along with first evaluations of their performance as compared to the most performing `single dimension' (either cross-sectional or temporal) forecast reconciliation procedures, is studied through a forecasting experiment on the 95 quarterly time series of the Australian GDP from Income and Expenditure sides considered by Athanasopoulos et al. (2019).


title: COT-GAN: Generating Sequential Data via Causal Optimal Transport

abstract: We introduce COT-GAN, an adversarial algorithm to train implicit generative models optimized for producing sequential data. The loss function of this algorithm is formulated using ideas from Causal Optimal Transport (COT), which combines classic optimal transport methods with an additional temporal causality constraint. Remarkably, we find that this causality condition provides a natural framework to parameterize the cost function that is learned by the discriminator as a robust (worst-case) distance, and an ideal mechanism for learning time dependent data distributions. Following Genevay et al.\ (2018), we also include an entropic penalization term which allows for the use of the Sinkhorn algorithm when computing the optimal transport cost. Our experiments show effectiveness and stability of COT-GAN when generating both low- and high-dimensional time series data. The success of the algorithm also relies on a new, improved version of the Sinkhorn divergence which demonstrates less bias in learning.


title: Flexible Dataset Distillation: Learn Labels Instead of Images

abstract: We study the problem of dataset distillation - creating a small set of synthetic examples capable of training a good model. In particular, we study the problem of label distillation - creating synthetic labels for a small set of real images, and show it to be more effective than the prior image-based approach to dataset distillation. Methodologically, we introduce a more robust and flexible meta-learning algorithm for distillation, as well as an effective first-order strategy based on convex optimization layers. Distilling labels with our new algorithm leads to improved results over prior image-based distillation. More importantly, it leads to clear improvements in flexibility of the distilled dataset in terms of compatibility with off-the-shelf optimizers and diverse neural architectures. Interestingly, label distillation can also be applied across datasets, for example enabling learning Japanese character recognition by training only on synthetically labeled English letters.


title: Neural Ensemble Search for Uncertainty Estimation and Dataset Shift

abstract: Ensembles of neural networks achieve superior performance compared to stand-alone networks in terms of accuracy, uncertainty calibration and robustness to dataset shift. \emph{Deep ensembles}, a state-of-the-art method for uncertainty estimation, only ensemble random initializations of a \emph{fixed} architecture. Instead, we propose two methods for automatically constructing ensembles with \emph{varying} architectures, which implicitly trade-off individual architectures' strengths against the ensemble's diversity and exploit architectural variation as a source of diversity. On a variety of classification tasks and modern architecture search spaces, we show that the resulting ensembles outperform deep ensembles not only in terms of accuracy but also uncertainty calibration and robustness to dataset shift. Our further analysis and ablation studies provide evidence of higher ensemble diversity due to architectural variation, resulting in ensembles that can outperform deep ensembles, even when having weaker average base learners. To foster reproducibility, our code is available: \url{https://github.com/automl/nes}


title: Large deviations of multichordal SLE$_{0+}$, real rational functions, and zeta-regularized determinants of Laplacians

abstract: We prove a strong large deviation principle (LDP) for multiple chordal SLE$_{0+}$ curves with respect to the Hausdorff metric. In the single-chord case, this result strengthens an earlier partial result by the second author. We also introduce a Loewner potential, which in the smooth case has a simple expression in terms of zeta-regularized determinants of Laplacians. This potential differs from the LDP rate function by an additive constant depending only on the boundary data, that satisfies PDEs arising as a semiclassical limit of the Belavin-Polyakov-Zamolodchikov equations of level two in conformal field theory with central charge $c \to -\infty$. Furthermore, we show that every multichord minimizing the potential in the upper half-plane for given boundary data is the real locus of a rational function and is unique, thus coinciding with the $\kappa \to 0+$ limit of the multiple SLE$_\kappa$. As a by-product, we provide an analytic proof of the Shapiro conjecture in real enumerative geometry, first proved by Eremenko and Gabrielov: if all critical points of a rational function are real, then the function is real up to post-composition by a M\"obius transformation.


title: The role of optimization geometry in single neuron learning

abstract: Recent numerical experiments have demonstrated that the choice of optimization geometry used during training can impact generalization performance when learning expressive nonlinear model classes such as deep neural networks. These observations have important implications for modern deep learning but remain poorly understood due to the difficulty of the associated nonconvex optimization problem. Towards an understanding of this phenomenon, we analyze a family of pseudogradient methods for learning generalized linear models under the square loss - a simplified problem containing both nonlinearity in the model parameters and nonconvexity of the optimization which admits a single neuron as a special case. We prove non-asymptotic bounds on the generalization error that sharply characterize how the interplay between the optimization geometry and the feature space geometry sets the out-of-sample performance of the learned model. Experimentally, selecting the optimization geometry as suggested by our theory leads to improved performance in generalized linear model estimation problems such as nonlinear and nonconvex variants of sparse vector recovery and low-rank matrix sensing.


title: Activity of the Eta-Aquariid and Orionid meteor showers

abstract: We present a multi-instrumental, multidecadal analysis of the activity of the Eta-Aquariid and Orionid meteor showers for the purpose of constraining models of 1P/Halley's meteoroid streams. The interannual variability of the showers' peak activity and period of duration is investigated through the compilation of published visual and radar observations prior to 1985 and more recent measurements reported in the International Meteor Organization (IMO) Visual Meteor DataBase, by the IMO Video Meteor Network and by the Canadian Meteor Orbit Radar (CMOR). These techniques probe the range of meteoroid masses from submilligrams to grams. The Eta-Aquariids and Orionids activity duration, shape, maximum zenithal hourly rates (ZHR) values, and the solar longitude of annual peaks since 1985 are analyzed. When available, annual activity profiles recorded by each detection network were measured and are compared. Observations from the three detection methods show generally good agreement in the showers' shape, activity levels, and annual intensity variations. Both showers display several activity peaks of variable location and strength with time. The Eta-Aquariids are usually two to three times stronger than the Orionids, but the two showers display occasional outbursts with peaks two to four times their usual activity level. CMOR observations since 2002 seem to support the existence of an ~12 year cycle in Orionids activity variations; however, additional and longer term radar and optical observations of the shower are required to confirm such periodicity.


title: Extension of the Generalized Hydrodynamics to the Dimensional Crossover Regime

abstract: In an effort to address integrability breaking in cold gas experiments, we extend the integrable hydrodynamics of the 1d Lieb-Liniger model with two additional components representing the population of atoms in the first and second transverse excited states, thus enabling a description of quasi-1d condensates. Collisions between different components are accounted for through the inclusion of a Boltzmann-type collision integral in the hydrodynamic equation. Contrary to standard generalized hydrodynamics, our extended model captures thermalization of the condensate at a rate consistent with experimental observations from a quantum Newton's cradle setup.


title: Quantum invariants of hyperbolic knots and extreme values of trigonometric products

abstract: In this paper we study the relation between the function $J_{4_1,0}$, which arises from a quantum invariant of the figure-eight knot, and Sudler's trigonometric product. We find $J_{4_1,0}$ up to a constant factor along continued fraction convergents to a quadratic irrational, and we show that its asymptotics deviates from the universal limiting behavior that has been found by Bettin and Drappeau in the case of large partial quotients. We relate the value of $J_{4_1,0}$ to that of Sudler's trigonometric product, and establish asymptotic upper and lower bounds for such Sudler products in response to a question of Lubinsky.


title: Influence of pairing and deformation on charge exchange transitions

abstract: We describe the importance of charge-exchange reactions, and in particular Gamow-Teller transitions, in astrophysical processes and double beta decay, and in understanding of nuclear structure. We first provide an overview of the central role played by the isovector pairing and the quadrupole-quadrupole channels in the description of energy spectra and in the manifestation of collective modes, some associated with deformation of the nuclear shape. We then turned our focus to Gamow-Teller (GT) transitions in relatively light nuclei, especially in the 2p1f shell, where isoscalar pairing may be playing a role in competition with the isovector pairing that dominates in heavier regions. Following a summary of the progress made in recent years on this subject, we report a systematic shell model study aimed at providing further clarification as to how these pairing modes compete. In this study, we use a schematic Hamiltonian that contains a quadrupole-quadrupole interaction as well as both isoscalar and isovector pairing interactions. We first find an optimal set of Hamiltonian parameters for the model, to provide a starting point from which to vary the relevant pairing strengths and thus assess how this impacts the behavior of GT transitions and the corresponding energy spectra and rotational properties of the various nuclei involved in the decays. The analysis includes as an important theme a comparison with experimental data. The need to suppress the isoscalar pairing mode when treating nuclei with a neutron excess to avoid producing spurious results for the ground state spin and parity with the simplified Hamiltonian is highlighted. Varying the strength parameters for the two pairing modes is found to exhibit different but systematic effects on GT transition properties and on the corresponding energy spectra, which are detailed. (abridged)


title: Uncertainty quantification for nonconvex tensor completion: Confidence intervals, heteroscedasticity and optimality

abstract: We study the distribution and uncertainty of nonconvex optimization for noisy tensor completion -- the problem of estimating a low-rank tensor given incomplete and corrupted observations of its entries. Focusing on a two-stage estimation algorithm proposed by Cai et al. (2019), we characterize the distribution of this nonconvex estimator down to fine scales. This distributional theory in turn allows one to construct valid and short confidence intervals for both the unseen tensor entries and the unknown tensor factors. The proposed inferential procedure enjoys several important features: (1) it is fully adaptive to noise heteroscedasticity, and (2) it is data-driven and automatically adapts to unknown noise distributions. Furthermore, our findings unveil the statistical optimality of nonconvex tensor completion: it attains un-improvable $\ell_{2}$ accuracy -- including both the rates and the pre-constants -- when estimating both the unknown tensor and the underlying tensor factors.

categories: stat.ML cs.IT cs.LG math.IT math.OC math.ST stat.TH


title: Is Working From Home The New Norm? An Observational Study Based on a Large Geo-tagged COVID-19 Twitter Dataset

abstract: As the COVID-19 pandemic swept over the world, people discussed facts, expressed opinions, and shared sentiments on social media. Since the reaction to COVID-19 in different locations may be tied to local cases, government regulations, healthcare resources and socioeconomic factors, we curated a large geo-tagged Twitter dataset and performed exploratory analysis by location. Specifically, we collected 650,563 unique geo-tagged tweets across the United States (50 states and Washington, D.C.) covering the date range from January 25 to May 10, 2020. Tweet locations enabled us to conduct region-specific studies such as tweeting volumes and sentiment, sometimes in response to local regulations and reported COVID-19 cases. During this period, many people started working from home. The gap between workdays and weekends in hourly tweet volumes inspired us to propose algorithms to estimate work engagement during the COVID-19 crisis. This paper also summarizes themes and topics of tweets in our dataset using both social media exclusive tools (i.e., #hashtags, @mentions) and the latent Dirichlet allocation model. We welcome requests for data sharing and conversations for more insights. Dataset link: http://covid19research.site/geo-tagged_twitter_datasets/


title: Blowup relations on $\mathbb{C}^2/\mathbb{Z}_2$ from Nakajima-Yoshioka blowup relations

abstract: We obtain bilinear relations on Nekrasov partition functions, arising from study of tau functions of quantum $q$-Painlev\'e equations, from Nakajima-Yoshioka blowup relations by an elementary algebraic approach. Additionaly, using this approach, we prove certain relations on Nekrasov partition functions, modified by Chern-Simons term.


title: Dynamical formation and stability of fermion-boson stars

abstract: Gravitationally bound structures composed by fermions and scalar particles known as fermion-boson stars are regular and static configurations obtained by solving the coupled Einstein-Klein-Gordon-Euler (EKGE) system. In this work, we discuss one possible scenario through which these fermion-boson stars may form by solving numerically the EKGE system under the simplifying assumption of spherical symmetry. Our initial configurations assume an already existing neutron star surrounded by an accreting cloud of a massive and complex scalar field. The results of our simulations show that once part of the initial scalar field is expelled via gravitational cooling the system gradually oscillates around an equilibrium configuration that is asymptotically consistent with a static solution of the system. The formation of fermion-boson stars for large positive values of the coupling constant in the self-interaction term of the scalar-field potential reveal the presence of a node in the scalar field. This suggests that a fermionic core may help stabilize configurations with nodes in the bosonic sector, as happens for purely boson stars in which the ground state and the first excited state coexist.


title: Sums of even powers of k-regulous functions

abstract: We provide an example of a nonnegative $k$-regulous function on $\mathbb{R}^n$ for $k\geq 1$ and $n \geq 2$ which cannot be written as a sum of squares of $k$-regulous functions. We then obtain lower bounds for Pythagoras numbers $p_{2d}(\mathcal{R}^k(\mathbb{R}^n))$ of $k$-regulous functions on $\mathbb{R}^n$ for $k\geq 1$ and $n\geq 2$. We also prove that the second Pythagoras number of the ring of $0$-regulous functions $\mathcal{R}^0(X)$ on an irreducible $0$-regulous affine variety $X$ is finite and bounded from above by $2^{\dim X}$.


title: Topological Dirac Semimetal Phase in Bismuth Based Anode Materials for Sodium-Ion Batteries

abstract: Bismuth has recently attracted interest in connection with Na-ion battery anodes due to its high volumetric capacity. It reacts with Na to form Na$_3$Bi which is a prototypical Dirac semimetal with a nontrivial electronic structure. Density-functional-theory based first-principles calculations are playing a key role in understanding the fascinating electronic structure of Na$_3$Bi and other topological materials. In particular, the strongly-constrained-and-appropriately-normed (SCAN) meta-generalized-gradient-approximation (meta-GGA) has shown significant improvement over the widely used generalized-gradient-approximation (GGA) scheme in capturing energetic, structural, and electronic properties of many classes of materials. Here, we discuss the electronic structure of Na$_3$Bi within the SCAN framework and show that the resulting Fermi velocities and {\it s}-band shift around the $\Gamma$ point are in better agreement with experiments than the corresponding GGA predictions. SCAN yields a purely spin-orbit-coupling (SOC) driven Dirac semimetal state in Na$_3$Bi in contrast with the earlier GGA results. Our analysis reveals the presence of a topological phase transition from the Dirac semimetal to a trivial band insulator phase in Na$_{3}$Bi$_{x}$Sb$_{1-x}$ alloys as the strength of the SOC varies with Sb content, and gives insight into the role of the SOC in modulating conduction properties of Na$_3$Bi.


title: Coherent Reconstruction of Multiple Humans from a Single Image

abstract: In this work, we address the problem of multi-person 3D pose estimation from a single image. A typical regression approach in the top-down setting of this problem would first detect all humans and then reconstruct each one of them independently. However, this type of prediction suffers from incoherent results, e.g., interpenetration and inconsistent depth ordering between the people in the scene. Our goal is to train a single network that learns to avoid these problems and generate a coherent 3D reconstruction of all the humans in the scene. To this end, a key design choice is the incorporation of the SMPL parametric body model in our top-down framework, which enables the use of two novel losses. First, a distance field-based collision loss penalizes interpenetration among the reconstructed people. Second, a depth ordering-aware loss reasons about occlusions and promotes a depth ordering of people that leads to a rendering which is consistent with the annotated instance segmentation. This provides depth supervision signals to the network, even if the image has no explicit 3D annotations. The experiments show that our approach outperforms previous methods on standard 3D pose benchmarks, while our proposed losses enable more coherent reconstruction in natural images. The project website with videos, results, and code can be found at: https://jiangwenpl.github.io/multiperson


title: On general-relativistic hydrogen and hydrogenic ions

abstract: This paper studies how the static non-linear electromagnetic-vacuum spacetime of a point nucleus with negative bare mass affects the self-adjointness of the general-relativistic Dirac Hamiltonian for a test electron, without and with an anomalous magnetic moment. The study interpolates between the previously studied extreme cases of a test electron in (a) the Reissner--Weyl--Nordstr\"om spacetime (Maxwell's electromagnetic vacuum), which supports a very strong curvature singularity with negative infinite bare mass, and (b) the Hoffmann spacetime (Born or Born--Infeld's electromagnetic vacuum) with vanishing bare mass, which features the mildest possible curvature singularity. The main conclusion reached is: {on electrostatic spacetimes of a point nucleus with a strictly negative bare mass} (which may be $-\infty$) essential self-adjointness fails unless the radial electric field diverges sufficiently fast at the nucleus and the anomalous magnetic moment of the electron is taken into account. Thus on the Hoffmann spacetime with (strictly) negative bare mass the Dirac Hamiltonian of a test electron, with or without anomalous magnetic moment, is not essentially self-adjoint. All these operators have self-adjoint extensions, though, with the usual essential spectrum $(-\infty,-\mEL c^2]\cup[\mEL c^2,\infty)$ and an infinite discrete spectrum located in the gap $(-\mEL c^2,\mEL c^2)$


title: A computing strategy and programs to resolve the Gerstenhaber Problem for commuting triples of matrices

abstract: We describe a MATLAB program that could produce a negative answer to the Gerstenhaber Problem by the construction of three commuting $n \times n$ matrices $A,B,C$ over a field $F$ such that the subalgebra $F[A,B,C]$ they generate has dimension greater than $n$. This problem has remained open for nearly 60 years, following Gerstenhaber's surprising result (Annals Math.) that $\dim F[A,B] \le n$ for any two commuting matrices $A,B$. The property fails for four or more commuting matrices. We also make the MATLAB files freely available.


title: The role of metadata in reproducible computational research

abstract: Reproducible computational research (RCR) is the keystone of the scientific method for in silico analyses, packaging the transformation of raw data to published results. In addition to its role in research integrity, RCR has the capacity to significantly accelerate evaluation and reuse. This potential and wide-support for the FAIR principles have motivated interest in metadata standards supporting RCR. Metadata provides context and provenance to raw data and methods and is essential to both discovery and validation. Despite this shared connection with scientific data, few studies have explicitly described the relationship between metadata and RCR. This article employs a functional content analysis to identify metadata standards that support RCR functions across an analytic stack consisting of input data, tools, notebooks, pipelines, and publications. Our article provides background context, explores gaps, and discovers component trends of embeddedness and methodology weight from which we derive recommendations for future work.


title: An EAGLE's View of Ex-situ Galaxy Growth

abstract: Modern observational and analytic techniques now enable the direct measurement of star formation histories and the inference of galaxy assembly histories. However, current theoretical predictions of assembly are not ideally suited for direct comparison with such observational data. We therefore extend the work of prior examinations of the contribution of ex-situ stars to the stellar mass budget of simulated galaxies. Our predictions are specifically tailored for direct testing with a new generation of observational techniques by calculating ex-situ fractions as functions of galaxy mass and morphological type, for a range of surface brightnesses. These enable comparison with results from large FoV IFU spectrographs, and increasingly accurate spectral fitting, providing a look-up method for the estimated accreted fraction. We furthermore provide predictions of ex-situ mass fractions as functions of galaxy mass, galactocentric radius and environment. Using $z=0$ snapshots from the 100cMpc$^3$ and 25cMpc$^3$ EAGLE simulations we corroborate the findings of prior studies, finding that ex-situ fraction increases with stellar mass for central and satellite galaxies in a stellar mass range of 2$\times$10$^{7}$ - 1.9$\times$10$^{12}$ M$_{\odot}$. For those galaxies of mass M$_*$>5$\times$10$^{8}$M$_{\odot}$, we find that the total ex-situ mass fraction is greater for more extended galaxies at fixed mass. When categorising satellite galaxies by their parent group/cluster halo mass we find that the ex-situ fraction decreases with increasing parent halo mass at fixed galaxy mass. This apparently counter-intuitive result may be the result of high passing velocities within large cluster halos inhibiting efficient accretion onto individual galaxies.


title: Monotone operator equilibrium networks

abstract: Implicit-depth models such as Deep Equilibrium Networks have recently been shown to match or exceed the performance of traditional deep networks while being much more memory efficient. However, these models suffer from unstable convergence to a solution and lack guarantees that a solution exists. On the other hand, Neural ODEs, another class of implicit-depth models, do guarantee existence of a unique solution but perform poorly compared with traditional networks. In this paper, we develop a new class of implicit-depth model based on the theory of monotone operators, the Monotone Operator Equilibrium Network (monDEQ). We show the close connection between finding the equilibrium point of an implicit network and solving a form of monotone operator splitting problem, which admits efficient solvers with guaranteed, stable convergence. We then develop a parameterization of the network which ensures that all operators remain monotone, which guarantees the existence of a unique equilibrium point. Finally, we show how to instantiate several versions of these models, and implement the resulting iterative solvers, for structured linear operators such as multi-scale convolutions. The resulting models vastly outperform the Neural ODE-based models while also being more computationally efficient. Code is available at http://github.com/locuslab/monotone_op_net.


title: Six textbook mistakes in computational physics

abstract: This article discusses several erroneous claims which appear in textbooks on numerical methods and computational physics. These are not typos or mistakes an individual author has made, but widespread misconceptions. In an attempt to stop these issues from further propagating, we discuss them here, along with some background comments. In each case, we also provide a correction, which is aimed at summarizing material that is known to experts but is frequently mishandled in the introductory literature. To make the mistakes and corrections easy to understand, we bring up specific examples drawn from elementary physics and math. We also take this opportunity to provide pointers to the specialist literature for readers who wish to delve into a given topic in more detail.


title: Strategies for particle resampling in PIC simulations

abstract: In particle-in-cell simulations, excessive or even unfeasible computational demands can be caused by the growth of the number of particles in the course of prolific ionization or cascaded pair production due to the effects of quantum electrodynamics. Here we discuss how one can organize a dynamic rearrangement of the ensemble to reduce the number of macroparticles, while maintaining acceptable sampling of an arbitrary particle distribution. The approaches of merging and thinning as well as their variants are discussed and the aspects of use are considered.


title: Feynman Lectures on the Strong Interactions

abstract: These twenty-two lectures, with exercises, comprise the extent of what was meant to be a full-year graduate-level course on the strong interactions and QCD, given at Caltech in 1987-88. The course was cut short by the illness that led to Feynman's death. Several of the lectures were finalized in collaboration with Feynman for an anticipated monograph based on the course. The others, while retaining Feynman's idiosyncrasies, are revised similarly to those he was able to check. His distinctive approach and manner of presentation are manifest throughout. Near the end he suggests a novel, nonperturbative formulation of quantum field theory in $D$ dimensions. Supplementary material is provided in appendices and ancillary files, including verbatim transcriptions of three lectures and the corresponding audiotaped recordings.


title: The challenge of measuring the phase function of debris disks. Application to HR\,4796

abstract: Abridged: Debris disks are valuable systems to study dust properties. Because they are optically thin at all wavelengths, we have direct access to the properties of dust grains. One very promising technique to study them is to measure their phase function. Disks that are highly inclined are promising targets as a wider range of scattering angles can be probed. The phase function is usually either inferred by comparing the observations to synthetic disk models assuming a parametrized phase function, or estimating it from the surface brightness of the disk. We argue here that the latter approach can be biased due to projection effects leading to an increase in column density along the major axis of a non flat disk. We present a novel approach to account for those column density effects. The method remains model dependent, as one still requires a disk model to estimate the density variations as a function of the scattering angle. This method allows us however to estimate the shape of the phase function without having to invoke any parametrized form. We apply our method to SPHERE/ZIMPOL observations of HR\,4796 and highlight the differences with previous measurements. Our modelling results suggest that the disk is not vertically flat at optical wavelengths. We discuss some of the caveats of the approach, mostly that our method remains blind to real local increase of the dust density, and that it cannot yet be readily applied to angular differential imaging observations. Similarly to previous studies on HR\,4796, we still cannot reconcile the full picture using a given scattering theory to explain the shape of the phase function, a long lasting problem for debris disks. Nonetheless, we argue that similar effects as the ones highlighted in this study can also bias the determination of the phase function in total intensity.


title: Constraining the cross-section of dark matter with giant radial arcs in galaxy clusters

abstract: We compare the statistics and morphology of giant arcs in galaxy clusters using N-body and non-radiative SPH simulations within the standard cold dark matter model and simulations where dark matter has a non-negligible probability of interaction (parametrized by its cross-section), i.e self-interacting dark matter (SIDM). We use a ray-tracing technique to produce a statistically large number of arcs around six simulated galaxy clusters at different redshifts. Since dark matter is more likely to interact in colliding clusters than in relaxed clusters, and this probability of interaction is largest in denser regions, we focus our analysis on radial arcs (which trace the lensing potential in the central region better than tangential arcs) in galaxy clusters which underwent (or are undergoing) a major merger. We find that self-interacting dark matter produces fewer radial arcs than standard cold dark matter but they are on average more magnified. We also appreciate differences in the arc morphology that could be used to statistically favor one model versus the other.


title: On the analytic representation of Newtonian systems

abstract: We show that the theory of self-adjoint differential equations can be used to provide a satisfactory solution of the inverse variational problem in classical mechanics. A Newtonian equation when transformed to the self-adjoint form allows one to find an appropriate Lagrangian representation (direct analytic representation) for it. On the other hand, the same Newtonian equation in conjunction with its adjoint provides a basis to construct a different Lagrangian representation (indirect analytic representation) for the system. We obtain the time-dependent Lagrangian of the damped Harmonic oscillator from the self-adjoint form of the equation of motion and at the same time identify the adjoint of the equation with the so called Bateman image equation with a view to construct a time-independent indirect Lagrangian representation. We provide a number of case studies to demonstrate the usefulness of the approach derived by us. We also present similar results for a number of nonlinear differential equations by using an integral representation of the Lagrangian function and make some useful comments.

categories: physics.class-ph math-ph math.MP nlin.SI


title: Towards practical differentially private causal graph discovery

abstract: Causal graph discovery refers to the process of discovering causal relation graphs from purely observational data. Like other statistical data, a causal graph might leak sensitive information about participants in the dataset. In this paper, we present a differentially private causal graph discovery algorithm, Priv-PC, which improves both utility and running time compared to the state-of-the-art. The design of Priv-PC follows a novel paradigm called sieve-and-examine which uses a small amount of privacy budget to filter out "insignificant" queries, and leverages the remaining budget to obtain highly accurate answers for the "significant" queries. We also conducted the first sensitivity analysis for conditional independence tests including conditional Kendall's tau and conditional Spearman's rho. We evaluated Priv-PC on 4 public datasets and compared with the state-of-the-art. The results show that Priv-PC achieves 10.61 to 32.85 times speedup and better utility.


title: "Notic My Speech" -- Blending Speech Patterns With Multimedia

abstract: Speech as a natural signal is composed of three parts - visemes (visual part of speech), phonemes (spoken part of speech), and language (the imposed structure). However, video as a medium for the delivery of speech and a multimedia construct has mostly ignored the cognitive aspects of speech delivery. For example, video applications like transcoding and compression have till now ignored the fact how speech is delivered and heard. To close the gap between speech understanding and multimedia video applications, in this paper, we show the initial experiments by modelling the perception on visual speech and showing its use case on video compression. On the other hand, in the visual speech recognition domain, existing studies have mostly modeled it as a classification problem, while ignoring the correlations between views, phonemes, visemes, and speech perception. This results in solutions which are further away from how human perception works. To bridge this gap, we propose a view-temporal attention mechanism to model both the view dependence and the visemic importance in speech recognition and understanding. We conduct experiments on three public visual speech recognition datasets. The experimental results show that our proposed method outperformed the existing work by 4.99% in terms of the viseme error rate. Moreover, we show that there is a strong correlation between our model's understanding of multi-view speech and the human perception. This characteristic benefits downstream applications such as video compression and streaming where a significant number of less important frames can be compressed or eliminated while being able to maximally preserve human speech understanding with good user experience.


title: Temporal Phenotyping using Deep Predictive Clustering of Disease Progression

abstract: Due to the wider availability of modern electronic health records, patient care data is often being stored in the form of time-series. Clustering such time-series data is crucial for patient phenotyping, anticipating patients' prognoses by identifying "similar" patients, and designing treatment guidelines that are tailored to homogeneous patient subgroups. In this paper, we develop a deep learning approach for clustering time-series data, where each cluster comprises patients who share similar future outcomes of interest (e.g., adverse events, the onset of comorbidities). To encourage each cluster to have homogeneous future outcomes, the clustering is carried out by learning discrete representations that best describe the future outcome distribution based on novel loss functions. Experiments on two real-world datasets show that our model achieves superior clustering performance over state-of-the-art benchmarks and identifies meaningful clusters that can be translated into actionable information for clinical decision-making.


title: Explaining Local, Global, And Higher-Order Interactions In Deep Learning

abstract: We present a simple yet highly generalizable method for explaining interacting parts within a neural network's reasoning process. First, we design an algorithm based on cross derivatives for computing statistical interaction effects between individual features, which is generalized to both 2-way and higher-order (3-way or more) interactions. We present results side by side with a weight-based attribution technique, corroborating that cross derivatives are a superior metric for both 2-way and higher-order interaction detection. Moreover, we extend the use of cross derivatives as an explanatory device in neural networks to the computer vision setting by expanding Grad-CAM, a popular gradient-based explanatory tool for CNNs, to the higher order. While Grad-CAM can only explain the importance of individual objects in images, our method, which we call Taylor-CAM, can explain a neural network's relational reasoning across multiple objects. We show the success of our explanations both qualitatively and quantitatively, including with a user study. We will release all code as a tool package to facilitate explainable deep learning.


title: Targeted Adversarial Perturbations for Monocular Depth Prediction

abstract: We study the effect of adversarial perturbations on the task of monocular depth prediction. Specifically, we explore the ability of small, imperceptible additive perturbations to selectively alter the perceived geometry of the scene. We show that such perturbations can not only globally re-scale the predicted distances from the camera, but also alter the prediction to match a different target scene. We also show that, when given semantic or instance information, perturbations can fool the network to alter the depth of specific categories or instances in the scene, and even remove them while preserving the rest of the scene. To understand the effect of targeted perturbations, we conduct experiments on state-of-the-art monocular depth prediction methods. Our experiments reveal vulnerabilities in monocular depth prediction networks, and shed light on the biases and context learned by them.


title: Neutrino (Antineutrino)-Nucleus Interactions in the Shallow- and Deep-Inelastic Scattering Regions

abstract: In $\nu/\bar{\nu}$-N/A interactions SIS is technically defined in terms of the four-momentum transfer to the hadronic system as non-resonant meson production with $Q^2 \lessapprox 1~GeV^2$. This non-resonant meson production intermixes with resonant meson production in a regime of similar effective hadronic mass W of the interaction. As $Q^2$ grows and surpasses this $\approx 1~GeV^2$ limit, non-resonant interactions begin to take place with quarks within the nucleon indicating the start of DIS region. SIS and DIS regions have received varying degrees of attention from the community. While the theoretical / phenomenological study of $\nu$-nucleon and $\nu$-nucleus DIS scattering is advanced, such studies of a large portion of the SIS region, particularly the SIS to DIS transition region, have hardly begun. Experimentally, the SIS and the DIS regions for $\nu$-nucleon scattering have minimal results and only in the experimental study of the $\nu$-nucleus DIS region are there significant results for some nuclei. Since current and future neutrino oscillation experiments have contributions from both higher W SIS and DIS kinematic regions and these regions are in need of both considerable theoretical and experimental study, this review will concentrate on these SIS to DIS transition and DIS kinematic regions surveying our knowledge and the current challenges.


title: Vulnerability Coverage for Secure Configuration

abstract: We present a novel idea on adequacy testing called ``{vulnerability coverage}.'' The introduced coverage measure examines the underlying software for the presence of certain classes of vulnerabilities often found in the National Vulnerability Database (NVD) website. The thoroughness of the test input generation procedure is performed through the adaptation of evolutionary algorithms namely Genetic Algorithms (GA) and Particle Swarm Optimization (PSO). The methodology utilizes the Common Vulnerability Scoring System (CVSS), a free and open industry standard for assessing the severity of computer system security vulnerabilities, as a fitness measure for test inputs generation. The outcomes of these evolutionary algorithms are then evaluated in order to identify the vulnerabilities that match a class of vulnerability patterns for testing purposes.


title: Detection of Coincidentally Correct Test Cases through Random Forests

abstract: The performance of coverage-based fault localization greatly depends on the quality of test cases being executed. These test cases execute some lines of the given program and determine whether the underlying tests are passed or failed. In particular, some test cases may be well-behaved (i.e., passed) while executing faulty statements. These test cases, also known as coincidentally correct test cases, may negatively influence the performance of the spectra-based fault localization and thus be less helpful as a tool for the purpose of automated debugging. In other words, the involvement of these coincidentally correct test cases may introduce noises to the fault localization computation and thus cause in divergence of effectively localizing the location of possible bugs in the given code. In this paper, we propose a hybrid approach of ensemble learning combined with a supervised learning algorithm namely, Random Forests (RF) for the purpose of correctly identifying test cases that are mislabeled to be the passing test cases. A cost-effective analysis of flipping the test status or trimming (i.e., eliminating from the computation) the coincidental correct test cases is also reported.


title: Vulnerability Coverage as an Adequacy Testing Criterion

abstract: Mainstream software applications and tools are the configurable platforms with an enormous number of parameters along with their values. Certain settings and possible interactions between these parameters may harden (or soften) the security and robustness of these applications against some known vulnerabilities. However, the large number of vulnerabilities reported and associated with these tools make the exhaustive testing of these tools infeasible against these vulnerabilities infeasible. As an instance of general software testing problem, the research question to address is whether the system under test is robust and secure against these vulnerabilities. This paper introduces the idea of ``vulnerability coverage,'' a concept to adequately test a given application for a certain classes of vulnerabilities, as reported by the National Vulnerability Database (NVD). The deriving idea is to utilize the Common Vulnerability Scoring System (CVSS) as a means to measure the fitness of test inputs generated by evolutionary algorithms and then through pattern matching identify vulnerabilities that match the generated vulnerability vectors and then test the system under test for those identified vulnerabilities. We report the performance of two evolutionary algorithms (i.e., Genetic Algorithms and Particle Swarm Optimization) in generating the vulnerability pattern vectors.


title: Are Bell-tests only about local incompatibility?

abstract: The view exists that Bell-tests would only be about local incompatibility of quantum observables and that quantum non-locality would be an unnecessary concept in physics. In this note, we emphasize that it is not incompatibility at the local level that is important for the violation of Bell-CHSH inequality, but incompatibility at the non-local level of the joint measurements. Hence, non-locality remains a necessary concept to properly interpret the outcomes of certain joint quantum measurements.


title: Is SARS-CoV-2 a new Frankenstein monster virus?

abstract: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a novel virus of beta coronavirus genus originated in Wuhan, China in December 2019, resulted in the pandemic spread of Corona Virus Disease 2019 (COVID-19) worldwide. This genus also contains SARS-CoV (originated in China 2002-2003) and MERS-CoV (found in Saudi Arabia 2012). The nucleotide sequences of SARS-CoV-2 is closer to SARS-CoV (with approximately 80% identity) than to MERS-CoV. Despite these similarities, SARS-CoV-2 has two main features over the other coronaviruses. The first is the high contagious rate and the second is the immune response evasion. The higher transmission ability makes this virus quickly spread worldwide with a high mortality rate and more economic losses. This review will provide an overview of the current knowledge on the role of some viral and host cell factors in higher transmission and contagious spread of SARS-CoV-2.


title: Bell's Theorem, Quantum Probabilities, and Superdeterminism

abstract: In this short survey article, I discuss Bell's theorem and some strategies that attempt to avoid the conclusion of non-locality. I focus on two that intersect with the philosophy of probability: (1) quantum probabilities and (2) superdeterminism. The issues they raised not only apply to a wide class of no-go theorems about quantum mechanics but are also of general philosophical interest.


title: Autofocusing technologies for whole slide imaging and automated microscopy

abstract: Whole slide imaging (WSI) has moved digital pathology closer to diagnostic practice in recent years. Due to the inherent tissue topography variability, accurate autofocusing remains a critical challenge for WSI and automated microscopy systems. The traditional focus map surveying method is limited in its ability to acquire a high degree of focus points while still maintaining high throughput. Real-time approaches decouple image acquisition from focusing, thus allowing for rapid scanning while maintaining continuous accurate focus. This work reviews the traditional focus map approach and discusses the choice of focus measure for focal plane determination. It also discusses various real-time autofocusing approaches including reflective-based triangulation, confocal pinhole detection, low-coherence interferometry, tilted sensor approach, independent dual sensor scanning, beam splitter array, phase detection, dual-LED illumination, and deep-learning approaches. The technical concepts, merits, and limitations of these methods are explained and compared to those of a traditional WSI system. This review may provide new insights for the development of high-throughput automated microscopy imaging systems that can be made broadly available and utilizable without loss of capacity.


title: Exact solutions of a damped harmonic oscillator in a time dependent noncommutative space

abstract: In this paper we have obtained the exact eigenstates of a two dimensional damped harmonic oscillator in time dependent noncommutative space. It has been observed that for some specific choices of the damping factor and the time dependent frequency of the oscillator, there exists interesting solutions of the time dependent noncommutative parameters following from the solutions of the Ermakov-Pinney equation. Further, these solutions enable us to get exact analytic forms for the phase which relates the eigenstates of the Hamiltonian with the eigenstates of the Lewis invariant. We then obtain expressions for the matrix elements of the coordinate operators raised to a finite arbitrary power. From these general results we then compute the expectation value of the Hamiltonian. The expectation values of the energy are found to vary with time for different solutions of the Ermakov-Pinney equation corresponding to different choices of the damping factor and the time dependent frequency of the oscillator.


title: Optomechanical non-reciprocity and its equivalence to antiresonance: Control of isolation frequency using mechanical drive

abstract: We demonstrate that optomechanical non-reciprocity is equivalent to the anti-resonance, often discussed in the context of coupled driven harmonic oscillators. We show that that suitable phase-difference between the cavity driving fields make the relevant optomechanical couplings complex, which leads to non-reciprocity in the field fluctuations and anti-resonance in average field amplitudes. This analogy with anti-resonance demonstrates that only for a particular frequency (the so-called isolation frequency) of input signal, maximum non-reciprocity can be achieved. In contrast to the previous studies, we here show that one can dynamically control this isolation frequency by applying a mechanical drive of suitable frequency to the membrane in the optomechanical setup.


title: Self-Supervised Domain Mismatch Estimation for Autonomous Perception

abstract: Autonomous driving requires self awareness of its perception functions. Technically spoken, this can be realized by observers, which monitor the performance indicators of various perception modules. In this work we choose, exemplarily, a semantic segmentation to be monitored, and propose an autoencoder, trained in a self-supervised fashion on the very same training data as the semantic segmentation to be monitored. While the autoencoder's image reconstruction performance (PSNR) during online inference shows already a good predictive power w.r.t. semantic segmentation performance, we propose a novel domain mismatch metric DM as the earth mover's distance between a pre-stored PSNR distribution on training (source) data, and an online-acquired PSNR distribution on any inference (target) data. We are able to show by experiments that the DM metric has a strong rank order correlation with the semantic segmentation within its functional scope. We also propose a training domain-dependent threshold for the DM metric to define this functional scope.


title: Learning to map source code to software vulnerability using code-as-a-graph

abstract: We explore the applicability of Graph Neural Networks in learning the nuances of source code from a security perspective. Specifically, whether signatures of vulnerabilities in source code can be learned from its graph representation, in terms of relationships between nodes and edges. We create a pipeline we call AI4VA, which first encodes a sample source code into a Code Property Graph. The extracted graph is then vectorized in a manner which preserves its semantic information. A Gated Graph Neural Network is then trained using several such graphs to automatically extract templates differentiating the graph of a vulnerable sample from a healthy one. Our model outperforms static analyzers, classic machine learning, as well as CNN and RNN-based deep learning models on two of the three datasets we experiment with. We thus show that a code-as-graph encoding is more meaningful for vulnerability detection than existing code-as-photo and linear sequence encoding approaches. (Submitted Oct 2019, Paper #28, ICST)


title: Using The Baryonic Tully-Fisher Relation to Measure $H_o$

abstract: We explore the use of the baryonic Tully-Fisher relation (bTFR) as a new distance indicator. Advances in near-IR imaging and stellar population models, plus precise rotation curves, have reduced the scatter in the bTFR such that distance is the dominant source of uncertainty. Using 50 galaxies with accurate distances from Cepheids or tip magnitude of the red giant branch, we calibrate the bTFR on a scale independent of $H_o$. We then apply this calibrated bTFR to 95 independent galaxies from the SPARC sample, using CosmicFlows-3 velocities, to deduce the local value of $H_o$. We find $H_o$ = 75.1 +/- 2.3 (stat) +/- 1.5 (sys) km s$^{-1}$ Mpc$^{-1}$.


title: Weinberg dimension-5 operator by vector-like lepton doublets

abstract: It is well known that a Weinberg dimension-5 operator for small neutrino masses can be realized at tree level in three types of renormalizable models: (i) the type-I seesaw mediated by fermion singlets, (ii) the type-II seesaw mediated by Higgs triplets, (iii) the type-III seesaw mediated by fermion triplets. We here point out such operator can be also induced at tree level by vector-like lepton doublets in association with unusual fermion singlets, Higgs triplets or fermion triplets. If these unusual fermion singlets, Higgs triplets or fermion triplets are heavy enough, their decays can generate a lepton asymmetry to explain the cosmic baryon asymmetry, meanwhile, the vector-like lepton doublets can lead to a novel inverse or linear seesaw with rich observable phenomena. We further specify our scenario can be naturally embedded into a grand unification theory without the conventional type-I, type-II or type-III seesaw.


title: Fe III emission in quasars: evidence for a dense turbulent medium

abstract: Recent improvements to atomic energy-level data allow, for the first time, accurate predictions to be made for the Fe III line emission strengths in the spectra of luminous, $L_\text{bol}=10^{46}-10^{48}$ erg/s, Active Galactic Nuclei. The Fe III emitting gas must be primarily photoionized, consistent with observations of line reverberation. We use CLOUDY models exploring a wide range of parameter space, together with 26,500 rest-frame ultraviolet spectra from the Sloan Digital Sky Survey, to constrain the physical conditions of the line emitting gas. The observed Fe III emission is best accounted for by dense ($n_H=10^{14}$ cm$^{-3}$) gas which is microturbulent, leading to smaller line optical depths and fluorescent excitation. Such high density gas appears to be present in the central regions of the majority of luminous quasars. Using our favoured model, we present theoretical predictions for the relative strengths of the Fe III UV34 $\lambda\lambda$1895,1914,1926 multiplet. This multiplet is blended with the Si III] $\lambda$1892 and C III] $\lambda$1909 emission lines and an accurate subtraction of UV34 is essential when using these lines to infer information about the physics of the broad line region in quasars.


title: The Origin of the Bimodal Distribution of Magnetic Fields in Early-type Stars

abstract: In early-type stars a fossil magnetic field may be generated during the star formation process or be the result of a stellar merger event. Surface magnetic fields are thought to be erased by (sub)surface convection layers, which typically leave behind weak disordered fields. However, if the fossil field is strong enough it can prevent the onset of (sub)surface convection and so be preserved onto the main sequence. We calculate the critical field strength at which this occurs, and find that it corresponds well with the lower limit amplitude of observed fields in strongly magnetised Ap/Bp stars ($\approx$ 300 G). The critical field strength is predicted to increase slightly during the main sequence evolution, which could also explain the observed decline in the fraction of magnetic stars. This supports the conclusion that the bimodal distribution of observed magnetic fields in early-type stars reflects two different field origin stories: strongly magnetic fields are fossils fields inherited from star formation or a merger event, and weak fields are the product of on-going dynamo action.


title: Time Delay Lens Modelling Challenge

abstract: In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.


title: Universal scattering laws for quiescent bouncing cosmology

abstract: Cosmological bounces occur in many gravity theories. We define singularity scattering maps relating large scale geometries before and after the bounce (assuming no BKL oscillations) and encoding microscopic details of the theory. By classifying all suitably local maps we uncover three universal laws: scaling of Kasner exponents, canonical transformation of matter, directional metric scaling. These are indeed obeyed by Bianchi I bounces in string theory, loop quantum cosmology and modified matter models; our classification then determines how inhomogeneities and anisotropies traverse bounces, and precisely extracts model-dependent degrees of freedom.


title: Quantifying the Stellar Halo's Response to the LMC's Infall with Spherical Harmonics

abstract: The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall.


title: The Formation of a Stellar Association in the NGC 7000/IC 5070 Complex: Results from Kinematic Analysis of Stars and Gas

abstract: We examine the clustering and kinematics of young stellar objects (YSOs) in the North America/Pelican Nebulae, as revealed by Gaia astrometry, in relation to the structure and motions of the molecular gas, as indicated in molecular line maps. The Gaia parallaxes and proper motions allow us to significantly refine previously published lists of YSOs, demonstrating that many of the objects previously thought to form a distributed population turn out to be non-members. The members are subdivided into at least 6 spatio-kinematic groups, each of which is associated with its own molecular cloud component or components. Three of the groups are expanding, with velocity gradients of 0.3-0.5 km/s/pc, up to maximum velocities of ~8 km/s away from the groups' centers. The two known O-type stars associated with the region, 2MASS J20555125+4352246 and HD 199579, are rapidly escaping one of these groups, following the same position-velocity relation as the low-mass stars. We calculate that a combination of gas expulsion and tidal forces from the clumpy distribution of molecular gas could impart the observed velocity gradients within the groups. However, on a global scale, the relative motions of the groups do not appear either divergent or convergent. The velocity dispersion of the whole system is consistent with the kinetic energy gained due to gravitational collapse of the complex. Most of the stellar population has ages similar to the free-fall timescales for the natal clouds. Thus, we suggest the nearly free-fall collapse of a turbulent molecular cloud as the most likely scenario for star formation in this complex.


title: Temperature and Metallicity Gradients in the Hot Gas Outflows of M82

abstract: We utilize deep Chandra X-ray Observatory imaging and spectra of M82, the prototype of a starbursting galaxy with a multiphase wind, to map the hot plasma properties along the minor axis of the galaxy. We extract spectra from 11 regions up to 2.5 kpc from the starbursting midplane and model the data as a multi-temperature, optically thin thermal plasma with contributions from a non-thermal (power-law) component and from charge exchange (CX). We examine the gradients in best-fit parameters, including the intrinsic column density, plasma temperature, metal abundances, and number density of the hot gas as a function of distance from the M82 nucleus. We find that the temperatures and number densities of the warm-hot and hot plasma peak at the starbursting ridge and decreases along the minor axis. The temperature and density profiles are inconsistent with spherical adiabatic expansion of a super-heated wind and suggest mass loading and mixing of the hot phase with colder material. Non-thermal emission is detected in all of the regions considered, and CX comprises 8-25% of the total absorption-corrected, broad-band (0.5-7 keV) X-ray flux. We show that the abundances of O, Ne, Mg, and Fe are roughly constant across the regions considered, while Si and S peak within 500 pc of the central starburst. These findings support a direct connection between the M82 superwind and the warm-hot, metal-rich circumgalactic medium (CGM).


title: Determining the nuclear neutron distribution from Coherent Elastic neutrino-Nucleus Scattering: current results and future prospects

abstract: Coherent elastic neutrino-nucleus scattering (CE$\nu$NS), a process recently measured for the first time at ORNL's Spallation Neutron Source, is directly sensitive to the weak form factor of the nucleus. The European Spallation Source (ESS), presently under construction, will generate the most intense pulsed neutrino flux suitable for the detection of CE$\nu$NS. In this paper we quantify its potential to determine the root mean square radius of the point-neutron distribution, for a variety of target nuclei and a suite of detectors. To put our results in context we also derive, for the first time, a constraint on this parameter from the analysis of the energy and timing data of the CsI detector at the COHERENT experiment.


title: Evidence from the H3 Survey that the Stellar Halo is Entirely Comprised of Substructure

abstract: In the $\Lambda$CDM paradigm the Galactic stellar halo is predicted to harbor the accreted debris of smaller systems. To identify these systems, the H3 Spectroscopic Survey, combined with $Gaia$, is gathering 6D phase-space and chemical information in the distant Galaxy. Here we present a comprehensive inventory of structure within 50 kpc from the Galactic center using a sample of 5684 giants at $|b|>40^{\circ}$ and $|Z|>2$ kpc. We identify known structures including the high-$\alpha$ disk, the in-situ halo (disk stars heated to eccentric orbits), Sagittarius (Sgr), $Gaia$-Sausage-Enceladus (GSE), the Helmi Streams, Sequoia, and Thamnos. Additionally, we identify the following new structures: (i) Aleph ([Fe/H]$=-0.5$), a low eccentricity structure that rises a surprising 10 kpc off the plane, (ii, iii) Arjuna ([Fe/H]$=-1.2$) and I'itoi ([Fe/H]$<-2$), which comprise the high-energy retrograde halo along with Sequoia, and (iv) Wukong ([Fe/H]$=-1.6$), a prograde phase-space overdensity chemically distinct from GSE. For each structure we provide [Fe/H], [$\alpha$/Fe], and orbital parameters. Stars born within the Galaxy are a major component at $|Z|\sim$2 kpc ($\approx$60$\%$), but their relative fraction declines sharply to $\lesssim$5$\%$ past 15 kpc. Beyond 15 kpc, $>$80$\%$ of the halo is built by two massive ($M_{\star}\sim10^{8}-10^{9}M_{\odot}$) accreted dwarfs: GSE ([Fe/H]$=-1.2$) within 25 kpc, and Sgr ([Fe/H]$=-1.0$) beyond 25 kpc. This explains the relatively high overall metallicity of the halo ([Fe/H]$\approx-1.2$). We attribute $\gtrsim$95$\%$ of the sample to one of the listed structures, pointing to a halo built entirely from accreted dwarfs and heating of the disk.