title: Cluster Analysis of a Symbolic Regression Search Space

abstract: In this chapter we take a closer look at the distribution of symbolic regression models generated by genetic programming in the search space. The motivation for this work is to improve the search for well-fitting symbolic regression models by using information about the similarity of models that can be precomputed independently from the target function. For our analysis, we use a restricted grammar for uni-variate symbolic regression models and generate all possible models up to a fixed length limit. We identify unique models and cluster them based on phenotypic as well as genotypic similarity. We find that phenotypic similarity leads to well-defined clusters while genotypic similarity does not produce a clear clustering. By mapping solution candidates visited by GP to the enumerated search space we find that GP initially explores the whole search space and later converges to the subspace of highest quality expressions in a run for a simple benchmark problem.


title: A Contrastive Learning Approach to Auroral Identification and Classification

abstract: Unsupervised learning algorithms are beginning to achieve accuracies comparable to their supervised counterparts on benchmark computer vision tasks, but their utility for practical applications has not yet been demonstrated. In this work, we present a novel application of unsupervised learning to the task of auroral image classification. Specifically, we modify and adapt the Simple framework for Contrastive Learning of Representations (SimCLR) algorithm to learn representations of auroral images in a recently released auroral image dataset constructed using image data from Time History of Events and Macroscale Interactions during Substorms (THEMIS) all-sky imagers. We demonstrate that (a) simple linear classifiers fit to the learned representations of the images achieve state-of-the-art classification performance, improving the classification accuracy by almost 10 percentage points over the current benchmark; and (b) the learned representations naturally cluster into more clusters than exist manually assigned categories, suggesting that existing categorizations are overly coarse and may obscure important connections between auroral types, near-earth solar wind conditions, and geomagnetic disturbances at the earth's surface. Moreover, our model is much lighter than the previous benchmark on this dataset, requiring in the area of fewer than 25\% of the number of parameters. Our approach exceeds an established threshold for operational purposes, demonstrating readiness for deployment and utilization.


title: Notes on Worldsheet-Like Variables for Cluster Configuration Spaces

abstract: We continue the exploration of various appearances of cluster algebras in scattering amplitudes and related topics in physics. The cluster configuration spaces generalize the familiar moduli space ${\mathcal M}_{0,n}$ to finite-type cluster algebras. We study worldsheet-like variables, which for classical types have also appeared in the study of the symbol alphabet of Feynman integrals. We provide a systematic derivation of these variables from $Y$-systems, which allows us to express the dihedral coordinates in terms of them and to write the corresponding cluster string integrals in compact forms. We mainly focus on the $D_n$ type and show how to reach the boundaries of the configuration space, and write the saddle-point equations in terms of these variables. Moreover, these variables make it easier to study various topological properties of the space using a finite-field method. We propose conjectures about quasi-polynomial point count, dimensions of cohomology, and the number of saddle points for the $D_n$ space up to $n=10$, which greatly extend earlier results.


title: Physics-Augmented Learning: A New Paradigm Beyond Physics-Informed Learning

abstract: Integrating physical inductive biases into machine learning can improve model generalizability. We generalize the successful paradigm of physics-informed learning (PIL) into a more general framework that also includes what we term physics-augmented learning (PAL). PIL and PAL complement each other by handling discriminative and generative properties, respectively. In numerical experiments, we show that PAL performs well on examples where PIL is inapplicable or inefficient.

categories: cs.LG physics.class-ph physics.comp-ph physics.data-an


title: Transport anomalies in the layered compound BaPt4Se6

abstract: We report a layered ternary selenide BaPt4Se6 featuring sesqui-selenide Pt2Se3 layers sandwiched by Ba atoms. The Pt2Se3 layers in this compound can be derived from the Dirac-semimetal PtSe2 phase with Se vacancies that form a honeycomb structure. This structure results in a Pt (VI) and Pt (II) mixed-valence compound with both PtSe6 octahedra and PtSe4 square net coordination configurations. Temperature dependent electrical transport measurements suggest two distinct anomalies: a resistivity crossover, mimic to the metal-insulator (M-I) transition at ~150K, and a resistivity plateau at temperatures below 10K. The resistivity crossover is not associated with any structural, magnetic or charge order modulated phase transitions. Magnetoresistivity, Hall and heat capacity measurements concurrently suggest an existing hidden state below 5K in this system. Angle-resolved photoemission spectroscopy measurements reveal a metallic state and no dramatic reconstruction of the electronic structure up to 200K.


title: The Impact of Sub-Resolution Porosity on Numerical Simulations of Multiphase Flow

abstract: Sub-resolution porosity (SRP) is an ubiquitous, yet often ignored, feature in Digital Rock Physics. It embodies the trade-off between image resolution and field-of-view, and it is a direct result of choosing an imaging resolution that is larger than the smallest pores in a heterogeneous rock sample. In this study, we investigate the impacts of SRP on multiphase flow in porous rocks. To do so, we use our newly developed Multiphase Micro-Continuum model to perform first-of-a-kind direct numerical simulations of two-phase flow in porous samples containing SRP. We show that SRP properties (porosity, permeability, wettability) can impact predicted absolute permeabilities, fluid breakthrough times, residual saturations, and relative permeabilities by factors of up to 2, 1.5, 3, and 20, respectively. In particular, our results reveal that SRP can function as a persistent connector preventing the formation of isolated wetting fluid domains during drainage, thus dramatically increasing relative permeabilities to both fluids at low saturations. Overall, our study confirms previous evidence that the influence of SRP cannot be disregarded without incurring significant errors in numerical predictions or experimental analyses of multiphase flow in heterogeneous porous media.


title: Asymmetric scattering between kinks and wobblers

abstract: The asymmetric scattering between wobblers and kinks in the standard $\phi^4$ model is numerically investigated in two different scenarios. First, the collision between wobblers with opposite phase is analyzed. Here, a destructive interference between the shape modes of the colliding wobblers takes place at the impact time. The second scenario involves the scattering between a wobbler and an (unexcited) kink. In this case the energy transfer from the wobbler to the kink can be examined. The dependence of the final velocities and wobbling amplitudes of the scattered wobblers on the collision velocity and on the initial wobbling amplitude is discussed. Both situations lead to very different fractal structures in the velocity diagrams.


title: Intra-Day Price Simulation with Generative Adversarial Modelling of the Order Flow

abstract: Intra-day price variations in financial markets are driven by the sequence of orders, called the order flow, that is submitted at high frequency by traders. This paper introduces a novel application of the Sequence Generative Adversarial Networks framework to model the order flow, such that random sequences of the order flow can then be generated to simulate the intra-day variation of prices. As a benchmark, a well-known parametric model from the quantitative finance literature is selected. The models are fitted, and then multiple random paths of the order flow sequences are sampled from each model. Model performances are then evaluated by using the generated sequences to simulate price variations, and we compare the empirical regularities between the price variations produced by the generated and real sequences. The empirical regularities considered include the distribution of the price log-returns, the price volatility, and the heavy-tail of the log-returns distributions. The results show that the order sequences from the generative model are better able to reproduce the statistical behaviour of real price variations than the sequences from the benchmark.


title: Parallel spinor flows on three-dimensional Cauchy hypersurfaces

abstract: The three-dimensional parallel spinor flow is the evolution flow defined by a parallel spinor on a globally hyperbolic Lorentzian four-manifold. We prove that, despite the fact that Lorentzian metrics admitting parallel spinors are not necessarily Ricci flat, the parallel spinor flow preserves the vacuum momentum and Hamiltonian constraints and therefore the Einstein and parallel spinor flows coincide on common initial data. Using this result, we provide an initial data characterization of parallel spinors on Ricci flat Lorentzian four-manifolds, which in turn yields the first initial data characterization of Ricci-flat pp-waves. Furthermore, we explicitly solve the left-invariant parallel spinor flow on simply connected Lie groups, obtaining along the way necessary and sufficient conditions for the flow to be immortal. These are, to the best of our knowledge, the first non-trivial examples of evolution flows of parallel spinors. Finally, we use some of these examples to construct families of $\eta\,$-Einstein cosymplectic structures and to produce solutions to the left-invariant Ricci flow in three dimensions. This suggests the intriguing possibility of using first-order hyperbolic spinorial flows to construct special solutions of curvature flows in Riemannian signature.


title: The Role of Communication Technology Across the Life Course: A Field Guide to Social Support in East York

abstract: We examine how Canadians living in the East York section of Toronto exchange social support. Just as we have had to deconstruct social support to understand its component parts, we now deconstruct how different types of communication technologies play socially supportive roles. We draw on 101 in-depth interviews conducted in 2013-2014 to shed light on the support networks of a sample of East York residents and discern the role of communication technologies in the exchange of different types of social support across age groups. Our findings show that not much has changed since the 1960s in terms of the social ties that our sample of East Yorkers have, and the types of support mobilized via social networks: companionship, small and large services, emotional aid, and financial support. What has changed is how communication technologies interweave in complex ways with different types of social ties (partners, siblings, friends, etc.) to mobilize social support. We found that with siblings and extended kin communication technologies could boost the frequency of interaction and help exchange support at a distance. With friendship ties, communication technologies provide a continuous, constant flow of interaction. We draw implications for theories of social support and for social policy linked to interventions aimed at helping vulnerable groups during the COVID-19 pandemic.


title: The eDiscovery Medicine Show

abstract: The practice of bloodletting gradually fell into disfavor as a growing body of scientific evidence showed its ineffectiveness and demonstrated the effectiveness of various pharmaceuticals for the prevention and treatment of certain diseases. At the same time, the patent medicine industry promoted ineffective remedies at medicine shows featuring entertainment, testimonials, and pseudo-scientific claims with all the trappings--but none of the methodology--of science. Today, many producing parties and eDiscovery vendors similarly promote obsolete technology as well as unvetted tools labeled "artificial intelligence" or "technology-assisted review," along with unsound validation protocols. This situation will end only when eDiscovery technologies and tools are subject to testing using the methods of information retrieval.


title: Quantum anomalous Hall effect from inverted charge transfer gap

abstract: A general mechanism is presented by which topological physics arises in strongly correlated systems without flat bands. Starting from a charge transfer insulator, topology emerges when the charge transfer energy between the cation and anion is reduced to invert the lower Hubbard band and the spin-degenerate charge transfer band. A universal low-energy theory is developed for the inversion of charge transfer gap in a quantum antiferromagnet. The inverted state is found to be a quantum anomalous Hall (QAH) insulator with non-coplanar magnetism. Interactions play two essential roles in this mechanism: producing the insulating gap and quasiparticle bands prior to the band inversion, and causing the change of magnetic order necessary for the QAH effect after inversion. Our theory explains the electric field induced transition from correlated insulator to QAH state in AB-stacked TMD bilayer MoTe$_2$/WSe$_2$.


title: Online Object Model Reconstruction and Reuse for Lifelong Improvement of Robot Manipulation

abstract: This work proposes a robotic pipeline for picking and constrained placement of objects without geometric shape priors. Compared to recent efforts developed for similar tasks, where every object was assumed to be novel, the proposed system recognizes previously manipulated objects and performs online model reconstruction and reuse. Over a lifelong manipulation process, the system keeps learning features of objects it has interacted with and updates their reconstructed models. Whenever an instance of a previously manipulated object reappears, the system aims to first recognize it and then register its previously reconstructed model given the current observation. This step greatly reduces object shape uncertainty allowing the system to even reason for parts of objects, which are currently not observable. This also results in better manipulation efficiency as it reduces the need for active perception of the target object during manipulation. To get a reusable reconstructed model, the proposed pipeline adopts: i) TSDF for object representation, and ii) a variant of the standard particle filter algorithm for pose estimation and tracking of the partial object model. Furthermore, an effective way to construct and maintain a dataset of manipulated objects is presented. A sequence of real-world manipulation experiments is performed. They show how future manipulation tasks become more effective and efficient by reusing reconstructed models of previously manipulated objects, which were generated during their prior manipulation, instead of treating objects as novel every time.


title: CMB lensing reconstruction biases from masking extragalactic sources

abstract: Observed Cosmic Microwave Background (CMB) maps are contaminated by foregrounds, some of which are usually masked to perform cosmological analyses. If masks are correlated to the lensing signal, such as those removing extragalactic emissions located in matter overdensities, measurements over the unmasked sky may give biased estimates. We quantify the impact of these mask-induced biases for the reconstructed CMB lensing auto- and cross-correlation power spectra with external matter tracers. We show that they arise both from changes in the lensing power, and via modifications to the reconstruction power spectrum corrections, $N_L^{(0)}$, $N_L^{(1)}$ and $N_L^{(3/2)}$). For direct masking of the CMB lensing field, we derive simple analytic models of the masking effect and show that it is potentially large. We show that mask-induced biases are significantly reduced by optimal filtering of the CMB maps in the lensing reconstruction. We test the resulting lensing power spectrum biases on numerical simulations, masking radio sources, and peaks of thermal Sunyaev-Zeldovich (tSZ) and cosmic infrared background (CIB) emission. For the lensing auto spectrum, masking biases can only be measured with a statistical significance $\lesssim 3\sigma$ for future data sets. The same applies to the cross-correlation power spectra between CMB lensing and tSZ and CIB even though biases are larger (up to ~30%). We find that masking tSZ-selected galaxy clusters leads to the largest mask biases, potentially detectable with high significance. We find that the calibration of cluster masses using CMB lensing, in particular for objects at $z\lesssim 0.6$, might be significantly affected by mask biases for near-future observations if the lensing signal recovered inside the mask holes is used without further corrections. Conversely, mass calibration of high redshift objects will still deliver unbiased results.


title: PDC-Net+: Enhanced Probabilistic Dense Correspondence Network

abstract: Establishing robust and accurate correspondences between a pair of images is a long-standing computer vision problem with numerous applications. While classically dominated by sparse methods, emerging dense approaches offer a compelling alternative paradigm that avoids the keypoint detection step. However, dense flow estimation is often inaccurate in the case of large displacements, occlusions, or homogeneous regions. In order to apply dense methods to real-world applications, such as pose estimation, image manipulation, or 3D reconstruction, it is therefore crucial to estimate the confidence of the predicted matches. We propose the Enhanced Probabilistic Dense Correspondence Network, PDC-Net+, capable of estimating accurate dense correspondences along with a reliable confidence map. We develop a flexible probabilistic approach that jointly learns the flow prediction and its uncertainty. In particular, we parametrize the predictive distribution as a constrained mixture model, ensuring better modelling of both accurate flow predictions and outliers. Moreover, we develop an architecture and an enhanced training strategy tailored for robust and generalizable uncertainty prediction in the context of self-supervised training. Our approach obtains state-of-the-art results on multiple challenging geometric matching and optical flow datasets. We further validate the usefulness of our probabilistic confidence estimation for the tasks of pose estimation, 3D reconstruction, image-based localization, and image retrieval. Code and models are available at https://github.com/PruneTruong/DenseMatching.


title: $f$-Cal: Calibrated aleatoric uncertainty estimation from neural networks for robot perception

abstract: While modern deep neural networks are performant perception modules, performance (accuracy) alone is insufficient, particularly for safety-critical robotic applications such as self-driving vehicles. Robot autonomy stacks also require these otherwise blackbox models to produce reliable and calibrated measures of confidence on their predictions. Existing approaches estimate uncertainty from these neural network perception stacks by modifying network architectures, inference procedure, or loss functions. However, in general, these methods lack calibration, meaning that the predictive uncertainties do not faithfully represent the true underlying uncertainties (process noise). Our key insight is that calibration is only achieved by imposing constraints across multiple examples, such as those in a mini-batch; as opposed to existing approaches which only impose constraints per-sample, often leading to overconfident (thus miscalibrated) uncertainty estimates. By enforcing the distribution of outputs of a neural network to resemble a target distribution by minimizing an $f$-divergence, we obtain significantly better-calibrated models compared to prior approaches. Our approach, $f$-Cal, outperforms existing uncertainty calibration approaches on robot perception tasks such as object detection and monocular depth estimation over multiple real-world benchmarks.


title: Extinction in the Large Magellanic Cloud Bar around NGC1854, NGC1856, and NGC1858

abstract: We report on the extinction properties in the fields around the clusters NGC 1854, NGC 1856, and NGC 1858 in the bar of the Large Magellanic Cloud. The colour-magnitude diagrams of the stars in all these regions show an elongated red giant clump that reveals a variable amount of extinction across these fields, ranging from Av~0.2 to Av~1.9, including Galactic foreground extinction. The extinction properties nonetheless are remarkably uniform. The slope of the reddening vectors measured in the (V-I,V) and (B-I,B) colour-magnitude planes is fully in line with the Av/E(B-V)~5.5 value found in the outskirts of 30 Dor. This indicates the presence of an additional grey extinction component in the optical requiring big grains to be about twice as abundant as in the diffuse Galactic interstellar medium (ISM). Areas of higher extinction appear to be systematically associated with regions of more intense star formation, as measured by the larger number of stars more massive than 8 Msun, thus making injection of big grains into the ISM by SNII explosion the likely mechanism at the origin of the observed grey extinction component.


title: Modeling Ephraim Chambers' Knowledge Structure from a Naive Standpoint

abstract: In the preface to his Cyclopaedia published in 1728 Ephraim Chambers offers readers a systematized structure of his attempt to produce a universal repository of human knowledge. Divided into an interconnected taxonomic tree and domain vocabulary, this structure forms the basis of one effort from the Metadata Research Center to study historical ontologies. The knowledge structure is being encoded into a Simple Knowledge Organization System (SKOS) form as well as a Web Ontology Language (OWL) version. This paper explores the expressive and functional differences between these SKOS and OWL versions of Chambers' knowledge structure. As part of this goal, the paper research focused on the construction and application of rules in each system to produce a more computationally ready representation of Chambers' structure in SKOS, which is more thesaurus-like, and OWL, which represents additional ontological nuances. First, studying the various textual aspects at the semantic, syntactic, and typographic levels allowed for the relationships between terms to manifest from which rules governing expression of the connections between elements developed. Second, because each language, SKOS and OWL, functionally expresses different logical relationships, their possibilities and limitations offer a ground for further analyzing the resultant knowledge structures; although, each stemmed from the same basic source of Chambers' text. Lastly this paper will examine rule making and expression in light of Paul Grice's theory of conversational implicature to understand how a naive agent formulates and applies these rules to a knowledge structure.


title: Unsolved Problems in ML Safety

abstract: Machine learning (ML) systems are rapidly increasing in size, are acquiring new capabilities, and are increasingly deployed in high-stakes settings. As with other powerful technologies, safety for ML should be a leading research priority. In response to emerging safety challenges in ML, such as those introduced by recent large-scale models, we provide a new roadmap for ML Safety and refine the technical problems that the field needs to address. We present four problems ready for research, namely withstanding hazards ("Robustness"), identifying hazards ("Monitoring"), reducing inherent model hazards ("Alignment"), and reducing systemic hazards ("Systemic Safety"). Throughout, we clarify each problem's motivation and provide concrete research directions.


title: Separating Circuits : Switching Lemmas and Random Restrictions

abstract: This was submitted as a final project for CS254B, taught by Li Yang Tan and Tom Knowles. The field of Circuit Complexity utilises careful analysis of Boolean Circuit Functions in order to extract meaningful information about a range of complexity classes. In particular, the complexity class $P / \text{Poly}$ has played a central role in much of the historical attempts to tackle the problem of whether solution and verification are equivalent i.e. the central $P$ versus $NP$ problem. Whilst circuits can potentially be easier to analyse than Turing Machines due to their non-uniform nature of computation (program size is allowed to depend on the input size), it is notoriously hard to establish lower bounds for them. In this report, we will touch upon several results published by Hastad, Sipser and Razborov that will highlight a dynamic interplay between circuit complexity and many of the central ideas of modern-day complexity theory, and in particular the central importance of Hastad's Switching Lemma.


title: Toxicity of Carbon Nanomaterials

abstract: The outstanding multidisciplinary applicability of nanomaterials has paved the path for the rapid advancement of nanoscience during the last few decades. Such technological progress subsequently results in an inevitable environmental exposure of nanomaterials. Presently, nanomaterials are employed in an extensive range of commercial products. Safe and sustainable incorporation of nanomaterials in industrial products requires a profound and comprehensive understanding of their potential toxicity. Among different nanomaterials, carbon nanomaterials marked its notable superiority toward the development of state-of-the-art nanotechnology due to the significant contribution of each of the carbon allotropes with varied dimensionality. The zero-dimensional fullerene, one-dimensional carbon nanotube, and two-dimensional graphene possess an exclusive combination of distinctive properties that are utilized in most of the nanotechnology-based products nowadays. However, potential risk factors are associated with the production and the use of carbon nanomaterials. Consequently, the number of studies regarding the assessment of the toxicity of these nanomaterials has increased rapidly in the past decade. This chapter will summarize the recent scientific efforts on the toxicity evaluation of different carbon nanomaterials.


title: A Remark on Mathieu's Series

abstract: We establish a new lower bound for Mathieu's series and present a new derivation of its expansions in terms of Riemann Zeta functions.


title: Strategies to reach ultra-high capacitance values for supercapacitors: materials design

abstract: This review paper highlights the recent developments in supercapacitors by pointing out the significance of appropriate electrode and device designs. We reported ten extremely high-performance supercapacitors with specific capacitance values among the highest available until now to the best of our knowledge. These state-of-the-art designs employing innovative electrode materials have been discussed along with their short descriptions. The supercapacitors collected here possess the most promising potential for facilitating next-generation systems in energy harvesting and storage. This review is just the surface that can help provide a pathway for supercapacitor research, which is still wide open for exploring and developing new advanced materials for energy applications of the future.


title: Click-through Rate Prediction with Auto-Quantized Contrastive Learning

abstract: Click-through rate (CTR) prediction becomes indispensable in ubiquitous web recommendation applications. Nevertheless, the current methods are struggling under the cold-start scenarios where the user interactions are extremely sparse. We consider this problem as an automatic identification about whether the user behaviors are rich enough to capture the interests for prediction, and propose an Auto-Quantized Contrastive Learning (AQCL) loss to regularize the model. Different from previous methods, AQCL explores both the instance-instance and the instance-cluster similarity to robustify the latent representation, and automatically reduces the information loss to the active users due to the quantization. The proposed framework is agnostic to different model architectures and can be trained in an end-to-end fashion. Extensive results show that it consistently improves the current state-of-the-art CTR models.


title: New Hybrid Techniques for Business Recommender Systems

abstract: Besides the typical applications of recommender systems in B2C scenarios such as movie or shopping platforms, there is a rising interest in transforming the human-driven advice provided e.g. in consultancy via the use of recommender systems. We explore the special characteristics of such knowledge-based B2B services and propose a process that allows to incorporate recommender systems into them. We suggest and compare several recommender techniques that allow to incorporate the necessary contextual knowledge (e.g. company demographics). These techniques are evaluated in isolation on a test set of business intelligence consultancy cases. We then identify the respective strengths of the different techniques and propose a new hybridisation strategy to combine these strengths. Our results show that the hybridisation leads to a substantial performance improvement over the individual methods.


title: Time modulated media with digital meta-atoms

abstract: We develop the theory of acoustic wave propagation in a waveguide containing an array of time modulated digital meta-atoms, showing the equivalence between this array and a homogeneous, time varying, dispersive material. In the limit of an adiabatic time variation we find some choices of meta-atom coupling strength lead to exceptional points, where one mode is exponentially amplified and the adiabatic approximation breaks down. In the highly non--adiabatic limit we derive the analogue of reflection coefficients for an abrupt change in the meta-atom coupling. Due to dispersion we have both conversion between different types of modes, and different directions of propagation, which can be tuned through modifying the programmed response of the digital meta--atoms.


title: Improved Beckner's inequality for axially symmetric functions on $\mathbb{S}^n$

abstract: In this article we present various uniqueness and existence results for Q-curvature type equations with a Paneitz operator on $\s^n$ in axially symmetric function spaces. In particular, we show uniqueness results for $n=6, 8$ and improve the best constant of Beckner's inequality in these dimensions for axially symmetric functions under the constraint that their centers of mass are at the origin. As a consequence, the associated first Szeg\"o limit theorem is also proven for axially symmetric functions.


title: Fine-tuning Vision Transformers for the Prediction of State Variables in Ising Models

abstract: Transformers are state-of-the-art deep learning models that are composed of stacked attention and point-wise, fully connected layers designed for handling sequential data. Transformers are not only ubiquitous throughout Natural Language Processing (NLP), but, recently, they have inspired a new wave of Computer Vision (CV) applications research. In this work, a Vision Transformer (ViT) is applied to predict the state variables of 2-dimensional Ising model simulations. Our experiments show that ViT outperform state-of-the-art Convolutional Neural Networks (CNN) when using a small number of microstate images from the Ising model corresponding to various boundary conditions and temperatures. This work opens the possibility of applying ViT to other simulations, and raises interesting research directions on how attention maps can learn about the underlying physics governing different phenomena.

categories: cs.CV cond-mat.stat-mech cs.LG physics.comp-ph


title: PSI: Constructing ad-hoc Simplices to Interpolate High-Dimensional Unstructured Data

abstract: Interpolating unstructured data using barycentric coordinates becomes infeasible at high dimensions due to the prohibitive memory requirements of building a Delaunay triangulation. We present a new algorithm to construct ad-hoc simplices that are empirically guaranteed to contain the target coordinates, based on a nearest neighbor heuristic and an iterative dimensionality reduction through projection. We use these simplices to interpolate the astrophysical cooling function $\Lambda$ and show that this new approach produces good results with just a fraction of the previously required memory.


title: Optical excitations of Skyrmions, knotted solitons, and defects in atoms

abstract: Analogies between non-trivial topologies of matter and light have inspired numerous studies, including defect formation in structured light and topological photonic band-structures. Three-dimensional topological objects of localized particle-like nature attract broad interest across discipline boundaries from elementary particle physics and cosmology to condensed matter physics. Here we show how simple structured light beams can be transformed into optical excitations of atoms with considerably more complex topologies representing three-dimensional particle-like Skyrmions. This construction can also be described in terms of linked Hopf maps, analogous to knotted solitons of the Skyrme-Faddeev model. We identify the transverse polarization density current as the effective magnetic gauge potential for the Chern-Simons helicity term. While we prepare simpler two-dimensional baby-Skyrmions and singular defects using the traditional Stokes vectors on the Poincar\'e sphere for light, particle-like topologies can only be achieved in the full optical hypersphere description that no longer discards the variation of the total electromagnetic phase of vibration.

categories: cond-mat.quant-gas hep-th physics.atom-ph physics.optics quant-ph


title: Lobbying Influence -- The Role of Money, Strategies and Measurements

abstract: Comparing the results for preference attainment, self-perceived influence and reputational influence, this paper analyzes the relationship between financial resources and lobbying influence. The empirical analysis builds on data from an original survey with 312 Swiss energy policy stakeholders combined with document data from multiple policy consultation submission processes. The results show that the distribution of influence varies substantially depending on the measure. While financial resources for political purposes predict influence across all measures, the relationship is positive only for some. An analysis of indirect effects sheds light on the potential mechanisms that translate financial resources into influence.


title: Distributed UAV-enabled zero-forcing cooperative jamming scheme for safeguarding future wireless networks

abstract: In this work, we investigate the impact of two cooperative unmanned aerial vehicle (UAV)-based jammers on the secrecy performance of a ground wireless network in the presence of an eavesdropper. For that purpose, we investigate the secrecy-area related metrics, Jamming Coverage and Jamming Efficiency. Moreover, we propose a hybrid metric, the so-called Weighted Secrecy Coverage (WSC) and a virtual distributed multiple-input-multiple-output (MIMO)-based zero-forcing precoding scheme to avoid the jamming effects on the legitimate receiver. For evaluating these metrics, we derive a closed-form position-based metric, the secrecy improvement. Our mathematical derivations and comparative simulations show that the proposed zero-forcing scheme leads to an improvement on the secrecy performance in terms of the WSC, and provides conditions for improvement of Jamming Efficiency. They also show positioning trends on the UAVs over a fixed orbit around the legitimate transmitter as well as power allocation trends for optimal secrecy.


title: All-Around Real Label Supervision: Cyclic Prototype Consistency Learning for Semi-supervised Medical Image Segmentation

abstract: Semi-supervised learning has substantially advanced medical image segmentation since it alleviates the heavy burden of acquiring the costly expert-examined annotations. Especially, the consistency-based approaches have attracted more attention for their superior performance, wherein the real labels are only utilized to supervise their paired images via supervised loss while the unlabeled images are exploited by enforcing the perturbation-based \textit{"unsupervised"} consistency without explicit guidance from those real labels. However, intuitively, the expert-examined real labels contain more reliable supervision signals. Observing this, we ask an unexplored but interesting question: can we exploit the unlabeled data via explicit real label supervision for semi-supervised training? To this end, we discard the previous perturbation-based consistency but absorb the essence of non-parametric prototype learning. Based on the prototypical network, we then propose a novel cyclic prototype consistency learning (CPCL) framework, which is constructed by a labeled-to-unlabeled (L2U) prototypical forward process and an unlabeled-to-labeled (U2L) backward process. Such two processes synergistically enhance the segmentation network by encouraging more discriminative and compact features. In this way, our framework turns previous \textit{"unsupervised"} consistency into new \textit{"supervised"} consistency, obtaining the \textit{"all-around real label supervision"} property of our method. Extensive experiments on brain tumor segmentation from MRI and kidney segmentation from CT images show that our CPCL can effectively exploit the unlabeled data and outperform other state-of-the-art semi-supervised medical image segmentation methods.


title: A framework for quantitative analysis of Computed Tomography images of viral pneumonitis: radiomic features in COVID and non-COVID patients

abstract: Purpose: to optimize a pipeline of clinical data gathering and CT images processing implemented during the COVID-19 pandemic crisis and to develop artificial intelligence model for different of viral pneumonia. Methods: 1028 chest CT image of patients with positive swab were segmented automatically for lung extraction. A Gaussian model developed in Python language was applied to calculate quantitative metrics (QM) describing well-aerated and ill portions of the lungs from the histogram distribution of lung CT numbers in both lungs of each image and in four geometrical subdivision. Furthermore, radiomic features (RF) of first and second order were extracted from bilateral lungs using PyRadiomic tools. QM and RF were used to develop 4 different Multi-Layer Perceptron (MLP) classifier to discriminate images of patients with COVID (n=646) and non-COVID (n=382) viral pneumonia. Results: The Gaussian model applied to lung CT histogram correctly described healthy parenchyma 94% of the patients. The resulting accuracy of the models for COVID diagnosis were in the range 0.76-0.87, as the integral of the receiver operating curve. The best diagnostic performances were associated to the model based on RF of first and second order, with 21 relevant features after LASSO regression and an accuracy of 0.81$\pm$0.02 after 4-fold cross validation Conclusions: Despite these results were obtained with CT images from a single center, a platform for extracting useful quantitative metrics from CT images was developed and optimized. Four artificial intelligence-based models for classifying patients with COVID and non-COVID viral pneumonia were developed and compared showing overall good diagnostic performances


title: In-phase oscillations from the cooperation of cellular and network positive feedback in synaptically-coupled oscillators

abstract: We study the emergent dynamics of a network of synaptically coupled slow-fast oscillators. Synaptic coupling provides a network-level positive feedback mechanism that cooperates with cellular-level positive feedback to ignite in-phase network oscillations. Using analytical bifurcation analysis, we prove that the Perron-Frobenius eigenvector of the network adjacency matrix fully controls the oscillation pattern locally in a neighborhood of a Hopf bifurcation. Besides shifting the focus from the spectral properties of the network Laplacian matrix to the network adjacency matrix, we discuss other key differences between synaptic and diffusive coupling.


title: Rational Approximations of Quasi-Periodic Problems via Projected Green's Functions

abstract: We introduce the projected Green's function technique to study quasi-periodic systems such as the Andre-Aubry-Harper (AAH) model and beyond. In particular, we use projected Green's functions to construct a "rational approximate" sequence of transfer matrix equations consistent with quasi-periodic topology, where convergence of these sequences corresponds to the existence of extended eigenfunctions. We motivate this framework by applying it to a few well studied cases such as the almost-Mathieu operator (AAH model), as well as more generic non-dual models that challenge standard routines. The technique is flexible and can be used to extract both analytic and numerical results, e.g. we analytically extract a modified phase diagram for Liouville irrationals. As a numerical tool, it does not require the fixing of boundary conditions and circumvents a primary failing of numerical techniques in quasi-periodic systems, extrapolation from finite size. Instead, it uses finite size scaling to define convergence bounds on the full irrational limit.

categories: cond-mat.mes-hall cond-mat.dis-nn math-ph math.MP quant-ph


title: Multiband gravitational wave cosmology with stellar origin black hole binaries

abstract: Massive stellar origin black hole binaries (SBHBs), originating from stars above the pair-instability mass gap, are primary candidates for multiband gravitational wave (GW) observations. Here we study the possibility to use them as effective dark standard sirens to constrain cosmological parameters. The long lasting inspiral signal emitted by these systems is accessible by the future $Laser \; Interferometer \; Space \; Antenna$ (LISA), while the late inspiral and merger are eventually detected by third generation ground-based telescopes such as the $Einstein \; Telescope$ (ET). The direct measurement of the luminosity distance and the sky position to the source, together with the inhomogeneous redshift distribution of possible host galaxies, allow us to infer cosmological parameters by probabilistic means. The efficiency of this statistical method relies in high parameter estimation performances. We show that this multiband approach allows a precise determination of the Hubble constant H$_0$ with just ${\cal O}(10)$ detected sources. For selected SBHB population models, assuming $4$ ($10$) years of LISA observations, we find that H$_0$ is typically determined at $\sim 2\%$ ($\sim 1.5\%$), whereas $\Omega_m$ is only mildly constrained with a typical precision of $30\%$ ($20\%$). We discuss the origin of some outliers in our final estimates and we comment on ways to reduce their presence.


title: Evolution of the Ultraviolet Upturn at $0.3<z<1$: exploring helium rich stellar populations

abstract: We measure the near-UV (rest-frame $\sim 2400$\AA) to optical color for early-type galaxies in 12 clusters at $0.3 < z < 1.0$. We show that this is a suitable proxy for the more common far-ultraviolet bandpass used to measure the ultraviolet upturn and find that the upturn is detected to $z=0.6$ in these data, in agreement with previous work. We find evidence that the strength of the upturn starts to wane beyond this redshift and largely disappears at $z=1$. Our data is most consistent with models where early-type galaxies contain minority stellar populations with non-cosmological helium abundances, up to around 46\%, formed at $z \geq 3$, resembling multiple stellar population globular clusters in our Galaxy. This suggests that elliptical galaxies and globular clusters share similar chemical evolution and star formation histories. The vast majority of the stellar mass in these galaxies also must have been in place at $z > 3$.


title: Diffusivity in force-free simulations of global magnetospheres

abstract: Abstract: Assuming that the numerical diffusivity triggered by violations of the force-free electrodynamics constraints is a proxy for the physical resistivity, we examine its impact on the overall dynamics of force-free aligned pulsar magnetospheres endowed with an equatorial current sheet. We assess the constraint violations as a diffusivity source. The effects of modifications on electric fields used to restore force-free conditions are not confined to the equatorial current sheet, but modify the magnetospheric dynamics on timescales shorter than the pulsar rotational period. These corrections propagate especially via a channel that was unexplored, namely, changes induced to the electric charge density, $\rho$. We quantify the global consequences of diffusivity by comparing different techniques to model $\rho$. By default, we combine a conservative $\rho$-evolution with hyperbolic/parabolic cleaning of inaccuracies in the Maxwell equations. As an alternative, we enforce a constrained evolution where $\rho$ is directly computed as the electric field divergence. The conservative approach reduces the Poynting flux dissipated in the equatorial current sheet by an order of magnitude, along with an increase of the pulsar luminosity driven by a shift of the Y-point location. The luminosity changes according to $L_{\rm Y}\propto \alpha^{0.11}$, where $\alpha$ is the ratio of diffusion to advection timescales, controlling the amount of (numerical) diffusivity. Our models suggest interpreting the luminosity dependence on the Y-point location as differences in resistivities encountered at the equatorial current sheet. Alternatively, they could be interpreted in terms of the pair formation multiplicity, $\kappa$, smaller diffusion being consistent with $\kappa\gg 1$.


title: Causality, Unitarity, and the Weak Gravity Conjecture

abstract: We consider the shift of charge-to-mass ratio for extremal black holes in the context of effective field theory, motivated by the Weak Gravity Conjecture. We constrain extremality corrections in different regimes subject to unitarity and causality constraints. In the asymptotic IR, we demonstrate that for any supersymmetric theory in flat space, and for all minimally coupled theories, logarithmic running at one loop pushes the Wilson coefficient of certain four-derivative operators to be larger at lower energies, guaranteeing the existence of sufficiently large black holes with $Q>M$. We identify two exceptional cases of nonsupersymmetric theories involving large numbers of light states and Planck-scale nonminimal couplings, in which the sign of the running is reversed, leading to black holes with negative corrections to $Q/M$ in the deep IR, but argue that these do not rule out extremal black holes as the requisite charged states for the WGC. We separately show that causality and unitarity imply that the leading threshold corrections to the effective action from integrating out massive states, in any weakly coupled theory, can be written as a sum of squares and is manifestly positive for black hole backgrounds. Quite beautifully, the shift in the extremal $Q/M$ ratio is directly proportional to the shift in the on-shell action, guaranteeing that these threshold corrections push $Q>M$ in compliance with the WGC. Our results apply for black holes with or without dilatonic coupling and charged under any number of ${\rm U}(1)$s.


title: Cosmology with the kinetic Sunyaev-Zeldovich effect: Independent of the optical depth and $\sigma_8$

abstract: The cosmological constraints from the kinetic Sunyaev-Zeldovich experiments are degenerate with the optical depth measurement, which is commonly known as the optical-depth degeneracy. In this work, we introduce a new statistic based on the first moment of relative velocity between pairs in a triplet, which is capable of constraining cosmological parameters independent of the optical depth, and $\sigma_8$. Using 22,000 $N$-body simulations from the Quijote suite, we quantify the information content in the new statistic using Fisher matrix forecast. We find that it is able to obtain strong constraints on the cosmological parameters, particularly on the summed neutrino mass. The constraints have a factor of 6.2-12.9, and 2.3-5.7 improvement on all cosmological model parameters when compared to those obtained from the mean pairwise velocity, and the redshift-space halo power spectrum, respectively. Thus the new statistic paves a way forward to constrain cosmological parameters independent of the optical depth and $\sigma_8$ using data from future kinetic Sunyaev-Zeldovich experiments alone.


title: Paradoxes before the paradox: surface gravity and the information loss problem

abstract: The information loss paradox is widely regarded as one of the biggest open problems in theoretical physics. Several classical and quantum features must be present to enable its formulation. First, an event horizon is needed to justify the objective status of tracing out degrees of freedom inside the black hole. Second, evaporation must be completed (or nearly completed) in finite time according to a distant observer, and thus the formation of the black hole should also occur in finite time. In spherical symmetry these requirements constrain the possible metrics strongly enough to obtain a unique black hole formation scenario and match their parameters with the semiclassical results. However, the two principal generalizations of surface gravity, the quantity that determines the Hawking temperature, do not agree with each other on the dynamical background. Neither can correspond to the emission of nearly-thermal radiation. We infer from this that the information loss problem cannot be consistently posed in its standard form.


title: Line-graph-lattice crystal structures of stoichiometric materials

abstract: The origin of many quantum-material phenomena is intimately related to the presence of flat electronic bands. In quantum simulation, such bands have been realized through line-graph lattices, a class of lattices known to exhibit flat bands. Based on that work, we conduct a high-throughput screening for line-graph lattices among the crystalline structures of the Materials Flatband Database and report on new candidates for line-graph materials and lattice models. In particular, we find materials with line-graph-lattice structures beyond the two most commonly known examples, the kagom\'e and pyrochlore lattices. We also identify materials which may exhibit flat topological bands. Finally, we examine the various line-graph lattices detected and highlight those with gapped flat bands and those most frequently represented among this set of materials. With the identification of real stoichiometric materials and theoretical lattice geometries, the results of this work may inform future studies of flat-band many-body physics in both condensed matter experiment and theory.


title: Annihilogenesis

abstract: We investigate a novel interplay between the decay and annihilation of a particle whose mass undergoes a large shift during a first order phase transition, leading to the particles becoming trapped in the false vacuum and enhancing their annihilation rates as the bubbles of true vacuum expand. This opens up a large region of the parameter space where annihilations can be important. We apply this scenario to baryogenesis, where we find that annihilations can be enhanced enough to generate the requires baryon asymmetry even for relatively tiny annihilation cross sections with modest CP asymmetries.


title: Probing Early Super-massive Black Hole Growth and Quasar Evolution with Near-infrared Spectroscopy of 37 Reionization-era Quasars at 6.3 < z <= 7.64

abstract: We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the MgII emission lines are in the range $(0.3-3.6)\times10^{9}\,M_{\odot}$, which requires massive seed black holes with masses $\gtrsim10^{3-4}\,M_{\odot}$, assuming Eddington accretion since $z=30$. The Eddington ratio distribution peaks at $\lambda_{\rm Edd}\sim0.8$ and has a mean of 1.08, suggesting high accretion rates for these quasars. The CIV - MgII emission line velocity differences in our sample show an increase of CIV blueshift towards higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The FeII/MgII flux ratios derived for these quasars up to $z=7.6$, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for $z>6.5$ quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the CIV line. Our sample yields a strong broad absorption line quasar fraction of $\sim$24%, higher than the fractions in lower redshift quasar samples, although this could be affected by small sample statistics and selection effects.


title: ArtPop: A Stellar Population and Image Simulation Python Package

abstract: We present Artificial Stellar Populations (ArtPop), an open-source Python package for synthesizing stellar populations and generating artificial images of fully populated stellar systems. The code is designed to be intuitive to use and as modular as possible, making it possible to use each of its functionalities independently or together. ArtPop has a wide range of scientific and pedagogical use cases, including the measurement of detection efficiencies in current and future imaging surveys, the calculation of integrated stellar population parameters, quantitative comparisons of isochrone models, and the development and validation of astronomical image processing algorithms. In this paper, we give an overview of the ArtPop package, provide simple coding examples to demonstrate its implementation, and present results from some potential applications of the code. We provide links to the source code that created each example and figure throughout the paper. ArtPop is under active development, and we welcome bug reports, feature requests, and code contributions from the community.


title: Evidence for unconventional superconductivity in twisted bilayer graphene

abstract: The emergence of superconductivity with doping from correlated insulators in magic-angle twisted bilayer graphene (MATBG) has raised the intriguing possibility that its pairing mechanism is distinct from that of conventional superconductors, as described by the Bardeen-Cooper-Schrieffer (BCS) theory. While there is now ample evidence for strong electronic correlations in MATBG, recent studies have claimed that unlike correlated insulators, superconductivity persists even when these interactions are partially screened. This suggests that the pairing in MATBG might be conventional in nature, a consequence of the large density of states of its nearly flat bands, perhaps phonon-mediated as in BCS superconductors. Here we combine tunneling and Andreev reflection spectroscopy with the scanning tunneling microscope (STM) to observe several key experimental signatures for unconventional superconductivity in MATBG. We show that the tunneling spectra below the transition temperature $T_c$ are inconsistent with that of a conventional s-wave superconductor, but rather resemble that of a nodal superconductor with an anisotropic pairing mechanism. We observe a large discrepancy between the tunneling energy gap $\Delta_T$, which far exceeds the mean-field BCS ratio (with $2\Delta_T/k_BT_c \sim 25$) and the energy gap $\Delta_{AR}$ extracted from Andreev reflection spectroscopy ($2\Delta_{AR}/k_BT_c \sim 6$). The tunneling gap persists even when superconductivity is suppressed, indicating its emergence from a pseudogap phase, with a suppressed density of states at the Fermi level. Moreover, the pseudogap state and superconductivity are both absent when MATBG is aligned with the hexagonal boron nitride (hBN) underneath. These findings and other observations reported here provide a preponderance of evidence for a non-BCS mechanism for superconductivity in MATBG.


title: Types of Transients in the Centers of Post-Starburst and Quiescent Balmer-Strong Galaxies

abstract: Tidal Disruption Events (TDEs) have been found to show a preference for post-starburst (PS) and quiescent Balmer-strong (QBS) galaxies. This preference can be used to help find TDEs in transient surveys. But what other transients might "contaminate" such a search, and by how much? We examine all reported transients coincident with the centers of galaxies in the French & Zabludoff (2018) catalog of spectroscopically confirmed PS and QBS galaxies and photometrically identified PS and QBS galaxy candidates. We find that TDEs and Type Ia supernovae (SNe) are the only types of transients classified in the centers of these galaxies (aside from one active galactic nucleus flare), with Type Ia SNe being 8.3+-0.2 times more prevalent than TDEs (1-sigma confidence bounds). This factor is ~2.7 times lower than in a control sample of quiescent galaxies. Narrowing the sample to spectroscopically confirmed QBS galaxies does not change these statistics much. In spectroscopically confirmed PS galaxies, however, TDEs are the ones that outnumber Type Ia SNe 2+-0.6 to 1. Unfortunately, there are few such galaxies in the catalog. By classifying transients from the entire catalog, three times more TDEs are expected to be found, but with a ~16-times larger Type Ia SN contamination. We use the public ZTF photometric archive to search for possibly missed TDEs in the French & Zabludoff (2018) galaxies. We find three unclassified clear transients -- none of which are likely missed TDEs based on their light-curve colors.


title: Dense matter equation of state and phase transitions from a Generalized Skyrme model

abstract: Skyrmion crystals are the field configurations which minimize the energy per baryon in the infinitely large topological charge sector of the Skyrme model, at least for sufficiently high density. They are, therefore, an important tool to describe the ground state of cold, symmetric nuclear matter at high density regimes. In this work, we analyze different crystalline phases and the existence of phase transitions between them within the generalized Skyrme model, with the ultimate goal of describing symmetric nuclear matter in a wide regime of densities. Furthermore, we propose a new energy-minimizing phase for densities lower than the nuclear saturation point ($n_0$) which also presents a good qualitative behavior in the zero density limit, thereby improving the description of strongly interacting matter in the region $n<n_0$.


title: Multi-gap topology of the Wilson loop operator in mirror symmetric insulators

abstract: We study the multi-gap topology of the periodic spectra of Wilson loop operators (WLOs) in mirror symmetric insulators. We develop two topological invariants each associated with a mirror-invariant gap in the Wilson loop spectrum. We propose that both topological invariants in combination determine the general higher-order bulk-boundary correspondence in 2D mirror symmetric, boundary-obstructed topological insulators. Finally, we demonstrate that these new multi-gap topological invariants apply to anomalous cases beyond those captured by the nested Wilson loop, and we subsequently develop an understanding of the correlation between WLOs along two orthogonal directions.


title: A fully Lagrangian, non-parametric bias model for dark-matter halos

abstract: We present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at the field level using the mass-weighted halo field in the AbacusSummit simulations at $z=0.5$. Unlike the perturbative halo bias model widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity $\delta_1$ in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts, and may give rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial of $\delta_1$ and other quantities. Especially for massive halos with $6\times10^{12}\ h^{-1}\ M_\odot$, the halo-to-mass ratio starts soaring up at $\delta_1>0$, substantially different from the predictions of the bias expansion. We show that for the halo masses we consider ($M>3\times10^{11}\ h^{-1}\ M_\odot$) a non-parametric halo-to-mass ratio as a function of $\delta_1$ and its local derivative $\nabla^2\delta_1$ can recover the halo power spectra to sub-percent level for wavenumbers $k=0.01-0.1\ h\ {\rm Mpc}^{-1}$ given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At $k<0.01\ h\ {\rm Mpc}^{-1}$ and for massive halos with $M>6\times10^{12}\ h^{-1}\ M_\odot$, our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions.


title: Eigenvalue repulsions in the quasinormal spectra of the Kerr-Newman black hole

abstract: We study the gravito-electromagnetic perturbations of the Kerr-Newman (KN) black hole metric and identify the two $-$ photon sphere and near-horizon $-$ families of quasinormal modes (QNMs) of the KN black hole, computing the frequency spectra (for all the KN parameter space) of the modes with the slowest decay rate. We uncover a novel phenomenon for QNMs that is unique to the KN system, namely eigenvalue repulsion between QNM families. Such a feature is common in solid state physics where \eg it is responsible for energy bands/gaps in the spectra of electrons moving in certain Schr\"odinger potentials. Exploiting the enhanced symmetries of the near-horizon limit of the near-extremal KN geometry we also develop a matching asymptotic expansion that allows us to solve the perturbation problem using separation of variables and provides an excellent approximation to the KN QNM spectra near extremality. The KN QNM spectra here derived are required not only to account for the gravitational emission in astrophysical environments, such as the ones probed by LIGO, Virgo and LISA, but also allow to extract observational implications on several new physics scenarios, such as mini-charged dark-matter or certain modified theories of gravity, degenerate with the KN solution at the scales of binary mergers.


title: Physical Characterization of Metal-rich Near-Earth Asteroids 6178 (1986 DA) and 2016 ED85

abstract: Metal-rich near-Earth asteroids (NEAs) represent a small fraction of the NEA population that is mostly dominated by S- and C-type asteroids. Because of this, their identification and study provide us with a unique opportunity to learn more about the formation and evolution of this particular type of bodies, as well as their relationship with meteorites found on Earth. We present near-infrared (NIR) spectroscopic data of NEAs 6178 (1986 DA) and 2016 ED85. We found that the spectral characteristics of these objects are consistent with those of metal-rich asteroids, showing red slopes, convex shapes, and a weak pyroxene absorption band at $\sim$0.93 $\mu$m. The compositional analysis showed that they have a pyroxene chemistry of Fs$_{40.6\pm3.3}$Wo$_{8.9\pm1.1}$ and a mineral abundance of $\sim$15% pyroxene and 85% metal. We determined that these objects were likely transported to the near-Earth space via the 5:2 mean motion resonance with Jupiter. Asteroid spectra were compared with the spectra of mesosiderites and bencubbinites. Differences in the NIR spectra and pyroxene chemistry suggest that bencubbinites are not good meteorite analogs. Mesosiderites were found to have a similar pyroxene chemistry and produced a good spectral match when metal was added to the silicate component. We estimated that the amounts of Fe, Ni, Co, and the platinum group metals present in 1986 DA could exceed the reserves worldwide.


title: Lithium in the lower red giant branch of 5 Galactic globular clusters

abstract: Lithium is one of the few elements produced during the Big Bang Nucleosynthesis in the early universe. Moreover, its fragility makes it useful as a proxy for stellar environmental conditions. As such, the lithium abundance in old systems is at the core of different astrophysical problems. Stars in the lower red giant branch allow studying globular clusters where main sequence stars are too faint to be observed. We use these stars to analyze the initial Li content of the clusters and compare it to cosmological predictions, to measure spreads in Li between different stellar populations, and to study signs of extra depletion in these giants. We use GIRAFFE spectra to measure the lithium and sodium abundances of lower red giant branch stars in 5 globular clusters. These cover an extensive range in metallicity, from [Fe/H]$\sim-0.7$ to [Fe/H]$\sim-2.3$ dex. We find that the lithium abundance in these lower red giant branch stars forms a plateau, with values from $\mathrm{A(Li)_{NLTE}}=0.84$ to $1.03$ dex, showing no clear correlation with metallicity. When using stellar evolutionary models to calculate the primordial abundance of these clusters, we recover values $\mathrm{A(Li)_{NLTE}}=2.1-2.3$ dex, consistent with the constant value observed in warm metal-poor halo stars, the Spite plateau. Additionally, we find no difference in the lithium abundance of first and second population stars in each cluster. We also report the discovery of a Li-rich giant in the cluster NGC3201, with $\mathrm{A(Li)_{NLTE}}=1.63\pm0.18$ dex, where the enrichment mechanism is probably pollution from external sources.


title: Electrical Networks, Lagrangian Grassmannians and Symplectic Groups

abstract: We refine the result of T. Lam \cite{L} on embedding the space $E_n$ of electrical networks on a planar graph with $n$ boundary points into the totally non-negative Grassmannian $\mathrm{Gr}_{\geq 0}(n-1,2n)$ by proving first that the image lands in $\mathrm{Gr}(n-1,V)\subset \mathrm{Gr}(n-1,2n)$ where $V\subset \mathbb{R}^{2n}$ is a certain subspace of dimension $2n-2$. The role of this reduction in the dimension of the ambient space is crucial for us. We show next that the image lands in fact inside the Lagrangian Grassmannian $\mathrm{LG}(n-1,V)\subset \mathrm{Gr}(n-1,V)$. As it is well known $\mathrm{LG}(n-1)$ can be identified with $\mathrm{Gr}(n-1,2n-2)\cap \mathbb{P} L$ where $L\subset \bigwedge^{n-1}\mathbb R^{2n-2}$ is a subspace of dimension equal to the Catalan number $C_n$, moreover it is the space of the fundamental representation of the symplectic group $Sp(2n-2)$ which corresponds to the last vertex of the Dynkin diagram. We show further that the linear relations cutting the image of $E_n$ out of $\mathrm{Gr}(n-1,2n)$ found in \cite{L} define that space $L$. This connects the combinatorial description of $E_n$ discovered in \cite{L} and representation theory of the symplectic group.


title: Modelling of ionising feedback with Smoothed Particle Hydrodynamics and Monte Carlo Radiative Transfer on a Voronoi grid

abstract: The ionising feedback of young massive stars is well known to influence the dynamics of the birth environment and hence plays an important role in regulating the star formation process in molecular clouds. For this reason, modern hydrodynamics codes adopt a variety of techniques accounting for these radiative effects. A key problem hampering these efforts is that the hydrodynamics are often solved using smoothed particle hydrodynamics (SPH), whereas radiative transfer is typically solved on a grid. Here we present a radiation-hydrodynamics (RHD) scheme combining the SPH code Phantom and the Monte Carlo Radiative Transfer (MCRT) code CMacIonize, using the particle distribution to construct a Voronoi grid on which the MCRT is performed. We demonstrate that the scheme successfully reproduces the well-studied problem of D-type H II region expansion in a uniform density medium. Furthermore, we use this simulation setup to study the robustness of the RHD code with varying choice of grid structure, density mapping method, and mass and temporal resolution. To test the scheme under more realistic conditions, we apply it to a simulated star-forming cloud reminiscing those in the Central Molecular Zone of our galaxy, in order to estimate the amount of ionised material that a single source could create. We find that a stellar population of several $10^3~\rm{M_{\odot}}$ is needed to noticeably ionise the cloud. Based on our results, we formulate a set of recommendations to guide the numerical setup of future and more complex simulations of star forming clouds.


title: Nanotrappy: An open-source versatile package for cold-atom trapping close to nanostructures

abstract: Trapping cold neutral atoms in close proximity to nanostructures has raised a large interest in recent years, pushing the frontiers of cavity-QED and boosting the emergence of the waveguide-QED field of research. The design of efficient dipole trapping schemes in evanescent fields is a crucial requirement and a difficult task. Here we present an open-source Python package for calculating optical trapping potentials for neutral atoms, especially in the vicinity of nanostructures. Given field distributions and for a variety of trap configurations, nanotrappy computes the three-dimensional trapping potentials as well as the trap properties, ranging from trap positions to trap frequencies and state-dependent light shifts. We demonstrate the versatility for various seminal structures in the field, e.g., optical nanofiber, alligator slow-mode photonic-crystal waveguide and microtoroid. This versatile package facilitates the systematic design of structures and provides a full characterization of trapping potentials with applications to coherent manipulation of atoms and quantum information science.

categories: quant-ph physics.atom-ph physics.comp-ph physics.optics


title: Mass Matters: No Evidence for Ubiquitous Lithium Production in Low-Mass Clump Giants

abstract: Known sources of lithium (Li) in the universe include the big bang, novae, asymptotic giant branch stars, and cosmic ray spallation. During their longer-lived evolutionary phases, stars are not expected to add to the Li budget of the Galaxy, but to largely deplete it. In this context, recent analyses of Li data from GALAH and LAMOST for field red clump (RC) stars have concluded that there is the need for a new production channel of Li, ubiquitous among low-mass stars, and that would be triggered on the upper red giant branch (RGB) or at helium ignition. This is distinct from the "Li-rich giant" problem and reflects bulk RC star properties. We provide an analysis of the GALAH Li data that accounts for the distribution of progenitor masses of field RC stars observed today. Such progenitors are different than today's field RGB stars. Using standard post main-sequence stellar evolution, we show that the distribution of Li among field RC giants as observed by GALAH is consistent with standard model predictions, and does not require new Li production mechanisms. Our model predicts a large fraction of very low Li abundances from low mass progenitors, with higher abundances from higher mass ones. Moreover, there should be a large number of upper limits for RC giants, and higher abundances should correspond to higher masses. The most recent GALAH data indeed confirm the presence of large numbers of upper limits, and a much lower mean Li abundance in RC stars, in concordance with our interpretation.


title: Bit Complexity of Jordan Normal Form and Spectral Factorization

abstract: We study the bit complexity of two related fundamental computational problems in linear algebra and control theory. Our results are: (1) An $\tilde{O}(n^{\omega+3}a+n^4a^2+n^\omega\log(1/\epsilon))$ time algorithm for finding an $\epsilon-$approximation to the Jordan Normal form of an integer matrix with $a-$bit entries, where $\omega$ is the exponent of matrix multiplication. (2) An $\tilde{O}(n^6d^6a+n^4d^4a^2+n^3d^3\log(1/\epsilon))$ time algorithm for $\epsilon$-approximately computing the spectral factorization $P(x)=Q^*(x)Q(x)$ of a given monic $n\times n$ rational matrix polynomial of degree $2d$ with rational $a-$bit coefficients having $a-$bit common denominators, which satisfies $P(x)\succeq 0$ for all real $x$. The first algorithm is used as a subroutine in the second one. Despite its being of central importance, polynomial complexity bounds were not previously known for spectral factorization, and for Jordan form the best previous best running time was an unspecified polynomial in $n$ of degree at least twelve \cite{cai1994computing}. Our algorithms are simple and judiciously combine techniques from numerical and symbolic computation, yielding significant advantages over either approach by itself.


title: Computing Ground State Properties with Early Fault-Tolerant Quantum Computers

abstract: Significant effort in applied quantum computing has been devoted to the problem of ground state energy estimation for molecules and materials. Yet, for many applications of practical value, additional properties of the ground state must be estimated. These include Green's functions used to compute electron transport in materials and the one-particle reduced density matrices used to compute electric dipoles of molecules. In this paper, we propose a quantum-classical hybrid algorithm to efficiently estimate such ground state properties with high accuracy using low-depth quantum circuits. We provide an analysis of various costs (circuit repetitions, maximal evolution time, and expected total runtime) as a function of target accuracy, spectral gap, and initial ground state overlap. This algorithm suggests a concrete approach to using early fault tolerant quantum computers for carrying out industry-relevant molecular and materials calculations.


title: High contrast imaging at the photon noise limit with self-calibrating WFS/C systems

abstract: High contrast imaging (HCI) systems rely on active wavefront control (WFC) to deliver deep raw contrast in the focal plane, and on calibration techniques to further enhance contrast by identifying planet light within the residual speckle halo. Both functions can be combined in an HCI system and we discuss a path toward designing HCI systems capable of calibrating residual starlight at the fundamental contrast limit imposed by photon noise. We highlight the value of deploying multiple high-efficiency wavefront sensors (WFSs) covering a wide spectral range and spanning multiple optical locations. We show how their combined information can be leveraged to simultaneously improve WFS sensitivity and residual starlight calibration, ideally making it impossible for an image plane speckle to hide from WFS telemetry. We demonstrate residual starlight calibration in the laboratory and on-sky, using both a coronagraphic setup, and a nulling spectro-interferometer. In both case, we show that bright starlight can calibrate residual starlight.


title: Stellar Population Astrophysics (SPA) with the TNG: Stock 2, a little-studied open cluster with an eMSTO

abstract: Stock 2 is a little-studied open cluster that shows an extended main-sequence turnoff (eMSTO). In order to investigate this phenomenon and characterise the cluster itself we performed high-resolution spectroscopy in the framework of the Stellar Population Astrophysics (SPA) project. We employed the High Accuracy Radial velocity Planet Searcher in North hemisphere spectrograph (HARPS-N) at the Telescopio Nazionale Galileo (TNG). We completed our observations with additional spectra taken with the Catania Astrophysical Observatory Spectrograph (CAOS). In total we observed 46 stars (dwarfs and giants), which represent, by far, the largest sample collected for this cluster to date. We provide the stellar parameters, extinction, radial and projected rotational velocities for most of the stars. Chemical abundances for 21 species with atomic numbers up to 56 have also been derived. We notice a differential reddening in the cluster field whose average value is 0.27 mag. It seems to be the main responsible for the observed eMSTO, since it cannot be explained as the result of different rotational velocities, as found in other clusters. We estimate an age for Stock 2 of 450$\pm$150 Ma which corresponds to a MSTO stellar mass of $\approx$2.8 M$_{\odot}$. The cluster mean radial velocity is around 8.0 km s$^{-1}$. We find a solar-like metallicity for the cluster, [Fe/H]=$-$0.07$\pm$0.06, compatible with its Galactocentric distance. MS stars and giants show chemical abundances compatible within the errors, with the exceptions of Barium and Strontium, which are clearly overabundant in giants, and Cobalt, which is only marginally overabundant. Finally, Stock 2 presents a chemical composition fully compatible with that observed in other open clusters of the Galactic thin disc.


title: Nuclear fusion catalyzed by doubly charged scalars: Implications for energy production

abstract: A number of popular extensions of the Standard Model of particle physics predict the existence of doubly charged scalar particles $X^{\pm\pm}$. Such particles may be long-lived or even stable. If exist, $X^{--}$ could form atomic bound states with light nuclei and catalyze their fusion by essentially eliminating the Coulomb barrier between them. Such an $X$-catalyzed fusion ($X$CF) process does not require high temperatures or pressure and may have important applications for energy production. A similar process of muon-catalyzed fusion ($\mu$CF) has been shown not to be a viable source of energy because of the sticking of negative muons to helium nuclei produced in the fusion of hydrogen isotopes, which stops the catalytic process. We analyze $X$CF in deuterium environments and show that the $X$-particles can only stick to $^6$Li nuclei, which are produced in the third-stage reactions downstream in the catalytic cycle. The corresponding sticking probability is very low, and, before getting bound to $^6$Li, each $X$-particle can catalyze $\sim 3.5\cdot 10^{9}$ fusion cycles, producing $\sim 7\cdot 10^{4}$ TeV of energy. We also discuss the ways of reactivating the $X$-particles from the Coulomb-bound (${\rm ^6Li}X$) states, which would allow re-using them in $X$CF reactions.

categories: hep-ph astro-ph.SR hep-ex nucl-ex nucl-th


title: Constraints on Kerr-Newman black holes from merger-ringdown gravitational-wave observations

abstract: We construct a template to model the post-merger phase of a binary black hole coalescence in the presence of a remnant $U(1)$ charge. We include the quasi-normal modes typically dominant during a binary black hole coalescence, $(\ell,m,n) = \{(2,2,0), (2,2,1)\}$ and also present analytical fits for the quasinormal mode frequencies of a Kerr-Newman black hole in terms of its spin and charge, here also including the $(3,3,0)$ mode. Aside from astrophysical electric charge, our template can accommodate extensions of the Standard Model, such as a dark photon. Applying the model to LIGO-Virgo detections, we find that we are unable to distinguish between the charged and uncharged hypotheses from a purely post-merger analysis of the current events. However, restricting the mass and spin to values compatible with the analysis of the full signal, we obtain a 90th percentile bound $\bar{q} < 0.33$ on the black hole charge-to-mass ratio, for the most favorable case of GW150914. Under similar assumptions, by simulating a typical loud signal observed by the LIGO-Virgo network at its design sensitivity, we assess that this model can provide a robust measurement of the charge-to-mass ratio only for values $\bar{q} \gtrsim 0.5$; here we also assume that the mode amplitudes are similar to the uncharged case in creating our simulated signal. Lower values, down to $\bar{q} \sim 0.3$, could instead be detected when evaluating the consistency of the pre-merger and post-merger emission.


title: Neutral Exotica at FASER$\nu$ and SND@LHC

abstract: The $(g-2)_\mu$ anomaly indicates that the second generation of leptons should have new interactions beyond the standard model. The high flux of $\nu_\mu$ and $\bar{\nu}_\mu$ at the forward experiments such as FASER$\nu$ and SND@LHC makes them suitable setups to search for new interactions of the second generation leptons. In this paper, we build a model in which the second generation left-handed leptons couple to a new right-handed neutrino, $N$ and a new Higgs doublet which also couples to the quarks. The scattering of high energy $\nu_\mu$ off nuclei can produce $N$. We investigate how forward experiments can test this model by looking for the $N$ production vertex followed by the displaced vertex of the $N$ decay. Discovering even a single such event can be a harbinger to look for the spectacular signals of the new Higgs doublet production at the LHC. We discuss the possibility of explaining the $(g-2)_\mu$ anomaly by adding more generations of $N$ which will lead to chain decays of $N$ and multiple leptons with distinct signals both at forward experiments and at the CMS and ATLAS detectors. Finally, we show that by adding a new light singlet scalar mixed with the neutral component of the new Higgs doublet ({\it i.e.,} 2HDM+$S$ model), the statistics of the data sample can be dramatically increased.


title: Smart at what cost? Characterising Mobile Deep Neural Networks in the wild

abstract: With smartphones' omnipresence in people's pockets, Machine Learning (ML) on mobile is gaining traction as devices become more powerful. With applications ranging from visual filters to voice assistants, intelligence on mobile comes in many forms and facets. However, Deep Neural Network (DNN) inference remains a compute intensive workload, with devices struggling to support intelligence at the cost of responsiveness.On the one hand, there is significant research on reducing model runtime requirements and supporting deployment on embedded devices. On the other hand, the strive to maximise the accuracy of a task is supported by deeper and wider neural networks, making mobile deployment of state-of-the-art DNNs a moving target. In this paper, we perform the first holistic study of DNN usage in the wild in an attempt to track deployed models and match how these run on widely deployed devices. To this end, we analyse over 16k of the most popular apps in the Google Play Store to characterise their DNN usage and performance across devices of different capabilities, both across tiers and generations. Simultaneously, we measure the models' energy footprint, as a core cost dimension of any mobile deployment. To streamline the process, we have developed gaugeNN, a tool that automates the deployment, measurement and analysis of DNNs on devices, with support for different frameworks and platforms. Results from our experience study paint the landscape of deep learning deployments on smartphones and indicate their popularity across app developers. Furthermore, our study shows the gap between bespoke techniques and real-world deployments and the need for optimised deployment of deep learning models in a highly dynamic and heterogeneous ecosystem.


title: An Accelerated Stochastic Gradient for Canonical Polyadic Decomposition

abstract: We consider the problem of structured canonical polyadic decomposition. If the size of the problem is very big, then stochastic gradient approaches are viable alternatives to classical methods, such as Alternating Optimization and All-At-Once optimization. We extend a recent stochastic gradient approach by employing an acceleration step (Nesterov momentum) in each iteration. We compare our approach with state-of-the-art alternatives, using both synthetic and real-world data, and find it to be very competitive.


title: A nonstandard-analytic proof of a theorem regarding noncommutative ergodic optimizations

abstract: In a previous article, we extended the notion of ergodic optimization to the setting of C*-dynamical systems of countable discrete groups. Among the key results of that paper was that given an action $G \stackrel{\Xi}{\curvearrowright} \mathfrak{M}$ of a countable discrete amenable group $G$ on a W*-probability space $(\mathfrak{M}, \rho)$ by $\rho$-preserving $*$-automorphisms of $\mathfrak{M}$, a positive element $x \in \mathfrak{M}$, and a right F{\o}lner sequence $\mathcal{F} = (F_k)_{k \in \mathbb{N} }$ for $G$, the sequence $$\left( \left\| \frac{1}{|F_k|} \sum_{g \in F_k} \Xi_g x \right\| \right)_{ k \in \mathbb{N} }$$ converges to a value $\Gamma(x)$ which can be described in the language of ergodic optimization. We provide here an alternate, more direct proof of that theorem using the tools of nonstandard analysis.


title: Guided Probabilistic Simulation of Complex Systems Toward Rare and Extreme Events

abstract: Simulation based or dynamic probabilistic risk assessment methodologies were primarily developed for proving a more realistic and complete representation of complex systems accident response. Such simulation based methodologies have proven to be particularly powerful for systems with control loops and complex interactions between its elements, be they hardware, software, or human, as they provide a natural probabilistic environment to include physical models of system behavior (e.g., coupled neutronics and thermal hydraulic codes for nuclear power plants), mechanistic models of materials or hardware systems to predict failure, and those of natural hazards. Despite the advancements in simulation based methodologies, the fundamental challenge still persists as the space of possible probabilistic system trajectories is nearly infinite in size in simulating even systems of relatively low complexity. In this paper, a framework is developed to identify rare and extreme events and enabling the use of reverse trajectories to trace failures (or other system states) to causes for potential mitigation actions. This framework consists of an Intelligent Guidance module, Trajectory Generation module and Physical Simulation module. The Intelligent Guidance module provides planning information to the Trajectory Generation module that creates scenarios by interacting with the Physical Simulation in its environment. In turn, system trajectories or scenarios are created and post processed to provide updating information to the Intelligent Guidance module or aggregate the results when stopping criteria are met. The objective of guided simulation is to control the growth of the scenario tree and to efficiently identify important scenarios that meet single or multiple criteria. We present several solution strategies, both qualitative and data driven for each module.


title: One-shot Key Information Extraction from Document with Deep Partial Graph Matching

abstract: Automating the Key Information Extraction (KIE) from documents improves efficiency, productivity, and security in many industrial scenarios such as rapid indexing and archiving. Many existing supervised learning methods for the KIE task need to feed a large number of labeled samples and learn separate models for different types of documents. However, collecting and labeling a large dataset is time-consuming and is not a user-friendly requirement for many cloud platforms. To overcome these challenges, we propose a deep end-to-end trainable network for one-shot KIE using partial graph matching. Contrary to previous methods that the learning of similarity and solving are optimized separately, our method enables the learning of the two processes in an end-to-end framework. Existing one-shot KIE methods are either template or simple attention-based learning approach that struggle to handle texts that are shifted beyond their desired positions caused by printers, as illustrated in Fig.1. To solve this problem, we add one-to-(at most)-one constraint such that we will find the globally optimized solution even if some texts are drifted. Further, we design a multimodal context ensemble block to boost the performance through fusing features of spatial, textual, and aspect representations. To promote research of KIE, we collected and annotated a one-shot document KIE dataset named DKIE with diverse types of images. The DKIE dataset consists of 2.5K document images captured by mobile phones in natural scenes, and it is the largest available one-shot KIE dataset up to now. The results of experiments on DKIE show that our method achieved state-of-the-art performance compared with recent one-shot and supervised learning approaches. The dataset and proposed one-shot KIE model will be released soo


title: Systematical studies of the E1 photon strength functions combining Skyrme-HFB+QRPA model and experimental giant dipole resonance properties

abstract: Valuable theoretical predictions of nuclear dipole excitations in the whole nuclear chart are of great interest for different applications, including in particular nuclear astrophysics. We present here the systematic study of the electric dipole (E1) photon strength functions (PSFs) combining the microscopic Hartree-Fock-Bogoliubov plus Quasiparticle Random Phase Approximation (HFB+QRPA) model and the parametrizations constrained by the available experimental giant dipole resonance (GDR) data. For about 10000 nuclei with 8<Z<124 lying between the proton and the neutron drip-lines on nuclear chart, the particle-hole strength distributions are computed using the HFB+QRPA model under the assumption of spherical symmetry and making use of the BSk27 Skyrme effective interaction derived from the most accurate HFB mass model (HFB-27) so far achieved. Large-scale calculations of the BSk27+QRPA E1 PSFs are performed in the framework of a specific folding procedure, in which three phenomenological improvements are considered. First, two interference factors are introduced and adjusted to reproduce at best the available experimental GDR data. Second, an empirical expression accounting for the deformation effect is applied to describe the peak splitting of the strength function. Third, the width of the strength function is corrected by a temperature-dependent term, which effectively increases the de-excitation photon strength function at low-energy. The E1 PSFs as well as the extracted GDR peaks and widths are compared with available experimental data. A relatively good agreement with data indicates the reliability of the calculations. Eventually, the astrophysical (n,g) rates for all the 10000 nuclei with 8<Z<124 are estimated using the present E1 PSFs. The resulting reaction rates are compared with previous BSk7+QRPA results and Gogny-HFB+QRPA predictions based on the D1M interaction.


title: Multicolor Ramsey numbers for Berge cycles

abstract: In this paper, for small uniformities, we determine the order of magnitude of the multicolor Ramsey numbers for Berge cycles of length $4$, $5$, $6$, $7$, $10$, or $11$. Our result follows from a more general setup which can be applied to other hypergraph Ramsey problems. Using this, we additionally determine the order of magnitude of the multicolor Ramsey number for Berge-$K_{a,b}$ for certain $a$, $b$, and uniformities.


title: Constructing Prediction Intervals Using the Likelihood Ratio Statistic

abstract: Statistical prediction plays an important role in many decision processes such as university budgeting (depending on the number of students who will enroll), capital budgeting (depending on the remaining lifetime of a fleet of systems), the needed amount of cash reserves for warranty expenses (depending on the number of warranty returns), and whether a product recall is needed (depending on the number of potentially life-threatening product failures). In statistical inference, likelihood ratios have a long history of use for decision making relating to model parameters (e.g., in evidence-based medicine and forensics). We propose a general prediction method, based on a likelihood ratio (LR) involving both the data and a future random variable. This general approach provides a way to identify prediction interval methods that have excellent statistical properties. For example, if a prediction method can be based on a pivotal quantity, our LR-based method will often identify it. For applications where a pivotal quantity does not exist, the LR-based method provides a procedure with good coverage properties for both continuous or discrete-data prediction applications.


title: Forecasting the COVID-19 vaccine uptake rate: An infodemiological study in the US

abstract: A year following the initial COVID-19 outbreak in China, many countries have approved emergency vaccines. Public-health practitioners and policymakers must understand the predicted populational willingness for vaccines and implement relevant stimulation measures. This study developed a framework for predicting vaccination uptake rate based on traditional clinical data-involving an autoregressive model with autoregressive integrated moving average (ARIMA)- and innovative web search queries-involving a linear regression with ordinary least squares/least absolute shrinkage and selection operator, and machine-learning with boost and random forest. For accuracy, we implemented a stacking regression for the clinical data and web search queries. The stacked regression of ARIMA (1,0,8) for clinical data and boost with support vector machine for web data formed the best model for forecasting vaccination speed in the US. The stacked regression provided a more accurate forecast. These results can help governments and policymakers predict vaccine demand and finance relevant programs.


title: Gaussian process error modeling for chiral effective-field-theory calculations of $np\leftrightarrow d\gamma$ at low energies

abstract: We calculate the energy-dependent cross section of the $np\leftrightarrow d\gamma$ process in chiral effective field theory and apply state-of-the-art tools for quantification of theory uncertainty. We focus on the low-energy regime, where the magnetic dipole and the electric dipole transitions cross over, including the range relevant for big-bang nucleosynthesis. Working with the leading one- and two-body electromagnetic currents, we study the order-by-order convergence of this observable in the chiral expansion of the nuclear potential. We find that the Gaussian process error model describes the observed convergence very well, allowing us to present Bayesian credible intervals for the truncation error with correlations between the cross sections at different energies taken into account. We obtain a 1$\sigma$ estimate of about 0.2\% for the uncertainty from the truncation of the nuclear potential. This is an important step towards calculations with statistically interpretable uncertainties for astrophysical reactions involving light nuclei.


title: Lorentz symmetry and high-energy neutrino astronomy

abstract: The search of violation of Lorentz symmetry, or Lorentz violation (LV), is an active research field. The effects of LV are expected to be very small and special systems are often used to search it. High-energy astrophysical neutrinos offer a unique system to search signatures of LV due to the three factors: high neutrino energy, long propagation distance, and the presence of quantum mechanical interference. In this brief review, we introduce tests of LV and summarize existing searches of LV using atmospheric and astrophysical neutrinos.


title: Competence-Aware Path Planning via Introspective Perception

abstract: Robots deployed in the real world over extended periods of time need to reason about unexpected failures, learn to predict them, and to proactively take actions to avoid future failures. Existing approaches for competence-aware planning are either model-based, requiring explicit enumeration of known failure modes, or purely statistical, using state- and location-specific failure statistics to infer competence. We instead propose a structured model-free approach to competence-aware planning by reasoning about plan execution failures due to errors in perception, without requiring a priori enumeration of failure sources or requiring location-specific failure statistics. We introduce competence-aware path planning via introspective perception (CPIP), a Bayesian framework to iteratively learn and exploit task-level competence in novel deployment environments. CPIP factorizes the competence-aware planning problem into two components. First, perception errors are learned in a model-free and location-agnostic setting via introspective perception prior to deployment in novel environments. Second, during actual deployments, the prediction of task-level failures is learned in a context-aware setting. Experiments in a simulation show that the proposed CPIP approach outperforms the frequentist baseline in multiple mobile robot tasks, and is further validated via real robot experiments in an environment with perceptually challenging obstacles and terrain.