Showing 1–50 of 441 results for author: Allen, B

Free Electron Theory for Thin Metal Films

Authors: Philip B. Allen

Abstract: Quantum free electrons, i.e. plane waves, with wavevector k, and occupancy constrained by the Pauli exclusion principle, are explained in all introductory texts about solids. A free-electron description works surprisingly well for many properties of "simple" metals. It is assumed that the interior of the metal is essentially infinite, and surfaces are presumed irrelevant. Over the past 30 years, e… ▽ More Quantum free electrons, i.e. plane waves, with wavevector k, and occupancy constrained by the Pauli exclusion principle, are explained in all introductory texts about solids. A free-electron description works surprisingly well for many properties of "simple" metals. It is assumed that the interior of the metal is essentially infinite, and surfaces are presumed irrelevant. Over the past 30 years, experiments that visualize surfaces have revolutionized solid state physics, stimulating new theory and applications. Therefore, a basic question is, how can the free electron picture be applied to properties of solids where surfaces play a prominent role? Various versions of an extended free-electron theory are used, but not always explained pedagogically. This paper focusses on idealized metallic films. Three versions (an oversimplified one and two stages of improvement) of a free-electron description of metal films are given. These versions are illustrated in detail for the specific example of a slab of aluminum with six layers of atoms. △ Less

Submitted 23 June, 2024; originally announced June 2024.

Comments: 12 pages, 11 figures; to be submitted to the American Journal of Physics

Source anisotropies and pulsar timing arrays

Authors: Bruce Allen, Deepali Agarwal, Joseph D. Romano, Serena Valtolina

Abstract: Pulsar timing arrays (PTA) hunt for gravitational waves (GW) by searching for the correlations that GWs induce in the time-of-arrival residuals from different pulsars. If the GW sources are of astrophysical origin, then they are located at discrete points on the sky. However, PTA data are often modeled, and subsequently analyzed, via a "standard Gaussian ensemble". That ensemble is obtained in the… ▽ More Pulsar timing arrays (PTA) hunt for gravitational waves (GW) by searching for the correlations that GWs induce in the time-of-arrival residuals from different pulsars. If the GW sources are of astrophysical origin, then they are located at discrete points on the sky. However, PTA data are often modeled, and subsequently analyzed, via a "standard Gaussian ensemble". That ensemble is obtained in the limit of an infinite density of vanishingly weak, Poisson-distributed sources. In this paper, we move away from that ensemble, to study the effects of two types of "source anisotropy". The first (a), which is often called "shot noise", arises because there are $N$ discrete GW sources at specific sky locations. The second (b) arises because the GW source positions are not a Poisson process, for example, because galaxy locations are clustered. Here, we quantify the impact of (a) and (b) on the mean and variance of the pulsar-averaged Hellings and Downs correlation. For conventional PTA sources, we show that the effects of shot noise (a) are much larger than the effects of clustering (b). △ Less

Submitted 23 June, 2024; originally announced June 2024.

Comments: 16 pages, 3 figures

Volume Above Distance Below with Boundary II

Authors: Brian Allen, Edward Bryden

Abstract: It was shown by B. Allen, R. Perales, and C. Sormani that on a closed manifold where the diameter of a sequence of Riemannian metrics is bounded, if the volume converges to the volume of a limit manifold, and the sequence of Riemannian metrics are $C^0$ converging from below then one can conclude volume preserving Sormani-Wenger Intrinsic Flat convergence. The result was extended to manifolds with… ▽ More It was shown by B. Allen, R. Perales, and C. Sormani that on a closed manifold where the diameter of a sequence of Riemannian metrics is bounded, if the volume converges to the volume of a limit manifold, and the sequence of Riemannian metrics are $C^0$ converging from below then one can conclude volume preserving Sormani-Wenger Intrinsic Flat convergence. The result was extended to manifolds with boundary by B. Allen and R. Perales by a doubling with necks procedure which produced a closed manifold and reduced the case with boundary to the case without boundary. The consequence of the doubling with necks procedure was requiring a stronger condition than necessary on the boundary. Using the estimates for the Sormani-Wenger Intrinsic Flat distance on manifolds with boundary developed by B. Allen and R. Perales, we show that only a bound on the area of the boundary is needed in order to conclude volume preserving intrinsic flat convergence for manifolds with boundary. We also provide an example which shows that one should not expect convergence without a bound on area. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 12 pages, comments welcome!

On the Scalar Curvature Compactness Conjecture in the Conformal Case

Authors: Brian Allen, Wenchuan Tian, Changliang Wang

Abstract: Is a sequence of Riemannian manifolds with positive scalar curvature, satisfying some conditions to keep the sequence reasonable, compact? What topology should one use for the convergence and what is the regularity of the limit space? In this paper we explore these questions by studying the case of a sequence of Riemannian manifolds which are conformal to the $n$-dimensional round sphere. We are a… ▽ More Is a sequence of Riemannian manifolds with positive scalar curvature, satisfying some conditions to keep the sequence reasonable, compact? What topology should one use for the convergence and what is the regularity of the limit space? In this paper we explore these questions by studying the case of a sequence of Riemannian manifolds which are conformal to the $n$-dimensional round sphere. We are able to show that the sequence of conformal factors are compact in several analytic senses and are able to establish $C^0$ convergence away from a singular set of small volume in a similar fashion as C. Dong. Under a bound on the total scalar curvature we are able to show that the limit conformal factor has weak positive scalar curvature in the sense of weakly solving the conformal positive scalar curvature equation. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 30 pages, comments welcome

Evaluating Class Membership Relations in Knowledge Graphs using Large Language Models

Authors: Bradley P. Allen, Paul T. Groth

Abstract: A backbone of knowledge graphs are their class membership relations, which assign entities to a given class. As part of the knowledge engineering process, we propose a new method for evaluating the quality of these relations by processing descriptions of a given entity and class using a zero-shot chain-of-thought classifier that uses a natural language intensional definition of a class. We evaluat… ▽ More A backbone of knowledge graphs are their class membership relations, which assign entities to a given class. As part of the knowledge engineering process, we propose a new method for evaluating the quality of these relations by processing descriptions of a given entity and class using a zero-shot chain-of-thought classifier that uses a natural language intensional definition of a class. We evaluate the method using two publicly available knowledge graphs, Wikidata and CaLiGraph, and 7 large language models. Using the gpt-4-0125-preview large language model, the method's classification performance achieves a macro-averaged F1-score of 0.830 on data from Wikidata and 0.893 on data from CaLiGraph. Moreover, a manual analysis of the classification errors shows that 40.9% of errors were due to the knowledge graphs, with 16.0% due to missing relations and 24.9% due to incorrectly asserted relations. These results show how large language models can assist knowledge engineers in the process of knowledge graph refinement. The code and data are available on Github. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 11 pages, 1 figure, 2 tables, accepted at the European Semantic Web Conference Special Track on Large Language Models for Knowledge Engineering, Hersonissos, Crete, GR, May 2024, for associated code and data, see https://github.com/bradleypallen/evaluating-kg-class-memberships-using-llms

ACM Class: I.2.7; I.2.4

Simulating Chemistry on Bosonic Quantum Devices

Authors: Rishab Dutta, Delmar G. A. Cabral, Ningyi Lyu, Nam P. Vu, Yuchen Wang, Brandon Allen, Xiaohan Dan, Rodrigo G. Cortiñas, Pouya Khazaei, Max Schäfer, Alejandro C. C. d. Albornoz, Scott E. Smart, Scott Nie, Michel H. Devoret, David A. Mazziotti, Prineha Narang, Chen Wang, James D. Whitfield, Angela K. Wilson, Heidi P. Hendrickson, Daniel A. Lidar, Francisco Pérez-Bernal, Lea F. Santos, Sabre Kais, Eitan Geva , et al. (1 additional authors not shown)

Abstract: Bosonic quantum devices offer a novel approach to realize quantum computations, where the quantum two-level system (qubit) is replaced with the quantum (an)harmonic oscillator (qumode) as the fundamental building block of the quantum simulator. The simulation of chemical structure and dynamics can then be achieved by representing or mapping the system Hamiltonians in terms of bosonic operators. In… ▽ More Bosonic quantum devices offer a novel approach to realize quantum computations, where the quantum two-level system (qubit) is replaced with the quantum (an)harmonic oscillator (qumode) as the fundamental building block of the quantum simulator. The simulation of chemical structure and dynamics can then be achieved by representing or mapping the system Hamiltonians in terms of bosonic operators. In this perspective, we review recent progress and future potential of using bosonic quantum devices for addressing a wide range of challenging chemical problems, including the calculation of molecular vibronic spectra, the simulation of gas-phase and solution-phase adiabatic and nonadiabatic chemical dynamics, the efficient solution of molecular graph theory problems, and the calculations of electronic structure. △ Less

Submitted 5 July, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: 40 pages including references, 13 figures, revised

Pulsar Timing Array Harmonic Analysis and Source Angular Correlations

Authors: Bruce Allen

Abstract: Gravitational waves (GWs) influence the arrival times of radio signals coming from pulsars. Here, we investigate the harmonic space approach to describing a pulsar's response to GWs. We derive and discuss the "diagonalized form" of the response, which is a sum of spin-2-weighted spherical harmonics of the GW direction multiplied by normal (spin-weight 0) spherical harmonics of the pulsar direction… ▽ More Gravitational waves (GWs) influence the arrival times of radio signals coming from pulsars. Here, we investigate the harmonic space approach to describing a pulsar's response to GWs. We derive and discuss the "diagonalized form" of the response, which is a sum of spin-2-weighted spherical harmonics of the GW direction multiplied by normal (spin-weight 0) spherical harmonics of the pulsar direction. We show how this allows many useful objects, for example, the Hellings and Downs two-point function, to be easily calculated. The approach also provides a clear description of the gauge dependence. We then employ this harmonic approach to model the effects of angular correlations in the sky locations of GW sources (sometimes called "statistical isotropy"). To do this, we construct rotationally invariant ensembles made up of many Gaussian subensembles, each of which breaks rotational invariance. Using harmonic techniques, we compute the cosmic covariance and the total covariance of the Hellings and Downs correlation in these models. The results may be used to assess the impact of angular source correlations on the Hellings and Downs correlation, and for optimal reconstruction of the Hellings and Downs curve in models where GW sources have correlated sky locations. △ Less

Submitted 29 April, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: Updated with citations to recent work by others; shifted to covariance C(O.O') to unify notation with [51]. 22 pages, 1 figure

SHROOM-INDElab at SemEval-2024 Task 6: Zero- and Few-Shot LLM-Based Classification for Hallucination Detection

Authors: Bradley P. Allen, Fina Polat, Paul Groth

Abstract: We describe the University of Amsterdam Intelligent Data Engineering Lab team's entry for the SemEval-2024 Task 6 competition. The SHROOM-INDElab system builds on previous work on using prompt programming and in-context learning with large language models (LLMs) to build classifiers for hallucination detection, and extends that work through the incorporation of context-specific definition of task,… ▽ More We describe the University of Amsterdam Intelligent Data Engineering Lab team's entry for the SemEval-2024 Task 6 competition. The SHROOM-INDElab system builds on previous work on using prompt programming and in-context learning with large language models (LLMs) to build classifiers for hallucination detection, and extends that work through the incorporation of context-specific definition of task, role, and target concept, and automated generation of examples for use in a few-shot prompting approach. The resulting system achieved fourth-best and sixth-best performance in the model-agnostic track and model-aware tracks for Task 6, respectively, and evaluation using the validation sets showed that the system's classification decisions were consistent with those of the crowd-sourced human labellers. We further found that a zero-shot approach provided better accuracy than a few-shot approach using automatically generated examples. Code for the system described in this paper is available on Github. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 6 pages, 6 figures, 4 tables, camera-ready copy, accepted to the 18th International Workshop on Semantic Evaluation (SemEval-2024), for associated code and data see https://github.com/bradleypallen/shroom

Standardizing Knowledge Engineering Practices with a Reference Architecture

Authors: Bradley P. Allen, Filip Ilievski

Abstract: Knowledge engineering is the process of creating and maintaining knowledge-producing systems. Throughout the history of computer science and AI, knowledge engineering workflows have been widely used given the importance of high-quality knowledge for reliable intelligent agents. Meanwhile, the scope of knowledge engineering, as apparent from its target tasks and use cases, has been shifting, togeth… ▽ More Knowledge engineering is the process of creating and maintaining knowledge-producing systems. Throughout the history of computer science and AI, knowledge engineering workflows have been widely used given the importance of high-quality knowledge for reliable intelligent agents. Meanwhile, the scope of knowledge engineering, as apparent from its target tasks and use cases, has been shifting, together with its paradigms such as expert systems, semantic web, and language modeling. The intended use cases and supported user requirements between these paradigms have not been analyzed globally, as new paradigms often satisfy prior pain points while possibly introducing new ones. The recent abstraction of systemic patterns into a boxology provides an opening for aligning the requirements and use cases of knowledge engineering with the systems, components, and software that can satisfy them best. This paper proposes a vision of harmonizing the best practices in the field of knowledge engineering by leveraging the software engineering methodology of creating reference architectures. We describe how a reference architecture can be iteratively designed and implemented to associate user needs with recurring systemic patterns, building on top of existing knowledge engineering workflows and boxologies. We provide a six-step roadmap that can enable the development of such an architecture, providing an initial design and outcome of the definition of architectural scope, selection of information sources, and analysis. We expect that following through on this vision will lead to well-grounded reference architectures for knowledge engineering, will advance the ongoing initiatives of organizing the neurosymbolic knowledge engineering space, and will build new links to the software architectures and data science communities. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 23 pages, 4 figures, 2 tables, camera-ready version, accepted for Transactions on Graph Data and Knowledge (TGDK)

Oberwolfach Report: Scalar Curvature Stability

Authors: Brian Allen

Abstract: Although scalar curvature is the simplest curvature invariant, our understanding of scalar curvature has not matured to the same level as Ricci or sectional curvature. Despite this fact, many rigidity phenomenon have been established which give some of the strongest insights into scalar curvature. Important examples include Geroch's conjecture, the positive mass theorem, and Llarull's theorem. In… ▽ More Although scalar curvature is the simplest curvature invariant, our understanding of scalar curvature has not matured to the same level as Ricci or sectional curvature. Despite this fact, many rigidity phenomenon have been established which give some of the strongest insights into scalar curvature. Important examples include Geroch's conjecture, the positive mass theorem, and Llarull's theorem. In order to further understand scalar curvature we ask corresponding geometric stability questions, where the hypotheses of the rigidity phenomenon are relaxed, and one would like to show that Riemannian manifolds which satisfy the relaxed conditions are close to the rigid objects in some topology. In this note we will survey what is known for scalar curvature stability, discuss what the questions are in this area, and introduce important tools which have been useful so far. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 7 pages, this report was written after the Oberwolfach Workshop entitled Analysis, Geometry and Topology of Positive Scalar Curvature Metrics which took place from February 19-23, 2024

Explainable Transformer Prototypes for Medical Diagnoses

Authors: Ugur Demir, Debesh Jha, Zheyuan Zhang, Elif Keles, Bradley Allen, Aggelos K. Katsaggelos, Ulas Bagci

Abstract: Deployments of artificial intelligence in medical diagnostics mandate not just accuracy and efficacy but also trust, emphasizing the need for explainability in machine decisions. The recent trend in automated medical image diagnostics leans towards the deployment of Transformer-based architectures, credited to their impressive capabilities. Since the self-attention feature of transformers contribu… ▽ More Deployments of artificial intelligence in medical diagnostics mandate not just accuracy and efficacy but also trust, emphasizing the need for explainability in machine decisions. The recent trend in automated medical image diagnostics leans towards the deployment of Transformer-based architectures, credited to their impressive capabilities. Since the self-attention feature of transformers contributes towards identifying crucial regions during the classification process, they enhance the trustability of the methods. However, the complex intricacies of these attention mechanisms may fall short of effectively pinpointing the regions of interest directly influencing AI decisions. Our research endeavors to innovate a unique attention block that underscores the correlation between 'regions' rather than 'pixels'. To address this challenge, we introduce an innovative system grounded in prototype learning, featuring an advanced self-attention mechanism that goes beyond conventional ad-hoc visual explanation techniques by offering comprehensible visual insights. A combined quantitative and qualitative methodological approach was used to demonstrate the effectiveness of the proposed method on the large-scale NIH chest X-ray dataset. Experimental results showed that our proposed method offers a promising direction for explainability, which can lead to the development of more trustable systems, which can facilitate easier and rapid adoption of such technology into routine clinics. The code is available at www.github.com/NUBagcilab/r2r_proto. △ Less

Submitted 11 March, 2024; originally announced March 2024.

89 New Ultracool Dwarf Co-Moving Companions Identified With The Backyard Worlds: Planet 9 Citizen Science Project

Authors: Austin Rothermich, Jacqueline K. Faherty, Daniella Bardalez-Gagliuffi, Adam C. Schneider, J. Davy Kirkpatrick, Aaron M. Meisner, Adam J. Burgasser, Marc Kuchner, Katelyn Allers, Jonathan Gagné, Dan Caselden, Emily Calamari, Mark Popinchalk, Genaro Suárez, Roman Gerasimov, Christian Aganze, Emma Softich, Chin-Chun Hsu, Preethi Karpoor, Christopher A. Theissen, Jon Rees, Rosario Cecilio-Flores-Elie, Michael C. Cushing, Federico Marocco, Sarah Casewell , et al. (21 additional authors not shown)

Abstract: We report the identification of 89 new systems containing ultracool dwarf companions to main sequence stars and white dwarfs, using the citizen science project Backyard Worlds: Planet 9 and cross-reference between Gaia and CatWISE2020. Thirty-two of these companions and thirty-three host stars were followed up with spectroscopic observations, with companion spectral types ranging from M7-T9 and ho… ▽ More We report the identification of 89 new systems containing ultracool dwarf companions to main sequence stars and white dwarfs, using the citizen science project Backyard Worlds: Planet 9 and cross-reference between Gaia and CatWISE2020. Thirty-two of these companions and thirty-three host stars were followed up with spectroscopic observations, with companion spectral types ranging from M7-T9 and host spectral types ranging from G2-M9. These systems exhibit diverse characteristics, from young to old ages, blue to very red spectral morphologies, potential membership to known young moving groups, and evidence of spectral binarity in 9 companions. Twenty of the host stars in our sample show evidence for higher order multiplicity, with an additional 11 host stars being resolved binaries themselves. We compare this sample's characteristics with those of the known stellar binary and exoplanet populations, and find our sample begins to fill in the gap between directly imaged exoplanets and stellary binaries on mass ratio-binding energy plots. With this study, we increase the population of ultracool dwarf companions to FGK stars by $\sim$42\%, and more than triple the known population of ultracool dwarf companions with separations larger than 1,000 au, providing excellent targets for future atmospheric retrievals. △ Less

Submitted 11 March, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 61 pages, 11 figures, 11 tables. Accepted for publication in AJ

Pulsar timing array source ensembles

Authors: Bruce Allen, Serena Valtolina

Abstract: The stochastic gravitational wave background for pulsar timing arrays is often modeled by a Gaussian ensemble which is isotropic and unpolarized. However, the Universe has a discrete set of polarized gravitational wave sources at specific sky locations. Can we trust that the Gaussian ensemble is an accurate description? To investigate this, we explicitly construct an ensemble containing $N$ indivi… ▽ More The stochastic gravitational wave background for pulsar timing arrays is often modeled by a Gaussian ensemble which is isotropic and unpolarized. However, the Universe has a discrete set of polarized gravitational wave sources at specific sky locations. Can we trust that the Gaussian ensemble is an accurate description? To investigate this, we explicitly construct an ensemble containing $N$ individual binary sources with circular orbits. The orbital inclination angles are randomly distributed, hence the individual sources are elliptically polarized. We then compute the first two moments of the Hellings and Downs correlation, as well as the pulsar-averaged correlation mean and (cosmic) variance. The first moments are the same as for a previously studied ensemble of circularly polarized sources. However, the second moments, and hence the variances, are different for the two ensembles. While neither discrete source model is exactly described by a Gaussian ensemble, we show that in the limit of large $N$, the differences are small. △ Less

Submitted 1 May, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: Final published version

Journal ref: Phys. Rev. D 109, 083038 (2024)

Optimization and Identification of Lattice Quantizers

Authors: Erik Agrell, Daniel Pook-Kolb, Bruce Allen

Abstract: Lattices with minimal normalized second moments are designed using a new numerical optimization algorithm. Starting from a random lower-triangular generator matrix and applying stochastic gradient descent, all elements are updated towards the negative gradient, which makes it the most efficient algorithm proposed so far for this purpose. A graphical illustration of the theta series, called theta i… ▽ More Lattices with minimal normalized second moments are designed using a new numerical optimization algorithm. Starting from a random lower-triangular generator matrix and applying stochastic gradient descent, all elements are updated towards the negative gradient, which makes it the most efficient algorithm proposed so far for this purpose. A graphical illustration of the theta series, called theta image, is introduced and shown to be a powerful tool for converting numerical lattice representations into their underlying exact forms. As a proof of concept, optimized lattices are designed in dimensions up to 16. In all dimensions, the algorithm converges to either the previously best known lattice or a better one. The dual of the 15-dimensional laminated lattice is conjectured to be optimal in its dimension and its exact normalized second moment is computed. △ Less

Submitted 23 June, 2024; v1 submitted 3 January, 2024; originally announced January 2024.

Conceptual Engineering Using Large Language Models

Authors: Bradley P. Allen

Abstract: We describe a method, based on Jennifer Nado's definition of classification procedures as targets of conceptual engineering, that implements such procedures using a large language model. We then apply this method using data from the Wikidata knowledge graph to evaluate concept definitions from two paradigmatic conceptual engineering projects: the International Astronomical Union's redefinition of… ▽ More We describe a method, based on Jennifer Nado's definition of classification procedures as targets of conceptual engineering, that implements such procedures using a large language model. We then apply this method using data from the Wikidata knowledge graph to evaluate concept definitions from two paradigmatic conceptual engineering projects: the International Astronomical Union's redefinition of PLANET and Haslanger's ameliorative analysis of WOMAN. We discuss implications of this work for the theory and practice of conceptual engineering. The code and data can be found on GitHub. △ Less

Submitted 30 November, 2023; originally announced December 2023.

Comments: 4 pages, 1 figure, Extended Abstract accepted for presentation at the 5th Conference on Philosophy of Artificial Intelligence (PhiAI 2023), for associated code and data see https://github.com/bradleypallen/zero-shot-classifiers-for-conceptual-engineering

ACM Class: I.2.7; I.2.4

Glued lattices are better quantizers than $K_{12}$

Authors: Erik Agrell, Daniel Pook-Kolb, Bruce Allen

Abstract: 40 years ago, Conway and Sloane proposed using the highly symmetrical Coxeter-Todd lattice $K_{12}$ for quantization, and estimated its second moment. Since then, all published lists identify $K_{12}$ as the best 12-dimensional lattice quantizer. Surprisingly, $K_{12}$ is not optimal: we construct two new 12-dimensional lattices with lower normalized second moments. The new lattices are obtained b… ▽ More 40 years ago, Conway and Sloane proposed using the highly symmetrical Coxeter-Todd lattice $K_{12}$ for quantization, and estimated its second moment. Since then, all published lists identify $K_{12}$ as the best 12-dimensional lattice quantizer. Surprisingly, $K_{12}$ is not optimal: we construct two new 12-dimensional lattices with lower normalized second moments. The new lattices are obtained by gluing together 6-dimensional lattices. △ Less

Submitted 23 June, 2024; v1 submitted 1 December, 2023; originally announced December 2023.

Exploring a Modification of $d_p$ Convergence

Authors: Brian Allen, Edward Bryden

Abstract: In the work by M. C. Lee, A. Naber, and R. Neumayer a beautiful $\varepsilon$-regularity theorem is proved under small negative scalar curvature and entropy bounds. In that paper, the $d_p$ distance for Riemannian manifolds is introduced and the quantitative stability results are given in terms of this notion of distance, with important examples showing why other existing notions of convergence ar… ▽ More In the work by M. C. Lee, A. Naber, and R. Neumayer a beautiful $\varepsilon$-regularity theorem is proved under small negative scalar curvature and entropy bounds. In that paper, the $d_p$ distance for Riemannian manifolds is introduced and the quantitative stability results are given in terms of this notion of distance, with important examples showing why other existing notions of convergence are not adequate in their setting. Due to the presence of an entropy bound, the possibility of long, thin splines forming along a sequence of Riemannian manifolds whose scalar curvature is becoming almost positive is ruled out. It is important to rule out such examples since the $d_p$ distance is not well behaved in the presences of splines that persist in the limit. Since there are many geometric stability conjectures where we want to allow for the presence of splines that persist in the limit, it is crucial to be able to modify the $d_p$ distance to retain its positive qualities and prevent it from being sensitive to splines. In this paper we explore one such modification of the $d_p$ distance and give a theorem which allows one to estimate the modified $d_p$ distance, which we expect to be useful in practice. △ Less

Submitted 24 June, 2024; v1 submitted 22 November, 2023; originally announced November 2023.

Comments: 22 pages, To appear in Nonlinear Analysis

Using Causal Threads to Explain Changes in a Dynamic System

Authors: Robert B. Allen

Abstract: We explore developing rich semantic models of systems. Specifically, we consider structured causal explanations about state changes in those systems. Essentially, we are developing process-based dynamic knowledge graphs. As an example, we construct a model of the causal threads for geological changes proposed by the Snowball Earth theory. Further, we describe an early prototype of a graphical inte… ▽ More We explore developing rich semantic models of systems. Specifically, we consider structured causal explanations about state changes in those systems. Essentially, we are developing process-based dynamic knowledge graphs. As an example, we construct a model of the causal threads for geological changes proposed by the Snowball Earth theory. Further, we describe an early prototype of a graphical interface to present the explanations. Unlike statistical approaches to summarization and explanation such as Large Language Models (LLMs), our approach of direct representation can be inspected and verified directly. △ Less

Submitted 19 November, 2023; originally announced November 2023.

Comments: 2023 ICADL

Efficient Algorithms for Monte Carlo Particle Transport on AI Accelerator Hardware

Authors: John Tramm, Bryce Allen, Kazutomo Yoshii, Andrew Siegel, Leighton Wilson

Abstract: The recent trend toward deep learning has led to the development of a variety of highly innovative AI accelerator architectures. One such architecture, the Cerebras Wafer-Scale Engine 2 (WSE-2), features 40 GB of on-chip SRAM, making it a potentially attractive platform for latency- or bandwidth-bound HPC simulation workloads. In this study, we examine the feasibility of performing continuous ener… ▽ More The recent trend toward deep learning has led to the development of a variety of highly innovative AI accelerator architectures. One such architecture, the Cerebras Wafer-Scale Engine 2 (WSE-2), features 40 GB of on-chip SRAM, making it a potentially attractive platform for latency- or bandwidth-bound HPC simulation workloads. In this study, we examine the feasibility of performing continuous energy Monte Carlo (MC) particle transport on the WSE-2 by porting a key kernel from the MC transport algorithm to Cerebras's CSL programming model. New algorithms for minimizing communication costs and for handling load balancing are developed and tested. The WSE-2 is found to run 130 times faster than a highly optimized CUDA version of the kernel run on an NVIDIA A100 GPU -- significantly outpacing the expected performance increase given the difference in transistor counts between the architectures. △ Less

Submitted 6 November, 2023; v1 submitted 3 November, 2023; originally announced November 2023.

ACM Class: D.1.3; J.2

Kernel-Elastic Autoencoder for Molecular Design

Authors: Haote Li, Yu Shee, Brandon Allen, Federica Maschietto, Victor Batista

Abstract: We introduce the Kernel-Elastic Autoencoder (KAE), a self-supervised generative model based on the transformer architecture with enhanced performance for molecular design. KAE is formulated based on two novel loss functions: modified maximum mean discrepancy and weighted reconstruction. KAE addresses the long-standing challenge of achieving valid generation and accurate reconstruction at the same… ▽ More We introduce the Kernel-Elastic Autoencoder (KAE), a self-supervised generative model based on the transformer architecture with enhanced performance for molecular design. KAE is formulated based on two novel loss functions: modified maximum mean discrepancy and weighted reconstruction. KAE addresses the long-standing challenge of achieving valid generation and accurate reconstruction at the same time. KAE achieves remarkable diversity in molecule generation while maintaining near-perfect reconstructions on the independent testing dataset, surpassing previous molecule-generating models. KAE enables conditional generation and allows for decoding based on beam search resulting in state-of-the-art performance in constrained optimizations. Furthermore, KAE can generate molecules conditional to favorable binding affinities in docking applications as confirmed by AutoDock Vina and Glide scores, outperforming all existing candidates from the training dataset. Beyond molecular design, we anticipate KAE could be applied to solve problems by generation in a wide range of applications. △ Less

Submitted 23 March, 2024; v1 submitted 12 October, 2023; originally announced October 2023.

Knowledge Engineering using Large Language Models

Authors: Bradley P. Allen, Lise Stork, Paul Groth

Abstract: Knowledge engineering is a discipline that focuses on the creation and maintenance of processes that generate and apply knowledge. Traditionally, knowledge engineering approaches have focused on knowledge expressed in formal languages. The emergence of large language models and their capabilities to effectively work with natural language, in its broadest sense, raises questions about the foundatio… ▽ More Knowledge engineering is a discipline that focuses on the creation and maintenance of processes that generate and apply knowledge. Traditionally, knowledge engineering approaches have focused on knowledge expressed in formal languages. The emergence of large language models and their capabilities to effectively work with natural language, in its broadest sense, raises questions about the foundations and practice of knowledge engineering. Here, we outline the potential role of LLMs in knowledge engineering, identifying two central directions: 1) creating hybrid neuro-symbolic knowledge systems; and 2) enabling knowledge engineering in natural language. Additionally, we formulate key open research questions to tackle these directions. △ Less

Submitted 1 October, 2023; originally announced October 2023.

Comments: 19 pages, 2 figures, accepted in Transactions on Graph Data and Knowledge

Data-Driven Approaches to Searches for the Technosignatures of Advanced Civilizations

Authors: T. Joseph W. Lazio, S. G. Djorgovski, Andrew Howard, Curt Cutler, Sofia Z. Sheikh, Stefano Cavuoti, Denise Herzing, Kiri Wagstaff, Jason T. Wright, Vishal Gajjar, Kevin Hand, Umaa Rebbapragada, Bruce Allen, Erica Cartmill, Jacob Foster, Dawn Gelino, Matthew J. Graham, Giuseppe Longo, Ashish A. Mahabal, Lior Pachter, Vikram Ravi, Gerald Sussman

Abstract: Humanity has wondered whether we are alone for millennia. The discovery of life elsewhere in the Universe, particularly intelligent life, would have profound effects, comparable to those of recognizing that the Earth is not the center of the Universe and that humans evolved from previous species. There has been rapid growth in the fields of extrasolar planets and data-driven astronomy. In a relati… ▽ More Humanity has wondered whether we are alone for millennia. The discovery of life elsewhere in the Universe, particularly intelligent life, would have profound effects, comparable to those of recognizing that the Earth is not the center of the Universe and that humans evolved from previous species. There has been rapid growth in the fields of extrasolar planets and data-driven astronomy. In a relatively short interval, we have seen a change from knowing of no extrasolar planets to now knowing more potentially habitable extrasolar planets than there are planets in the Solar System. In approximately the same interval, astronomy has transitioned to a field in which sky surveys can generate 1 PB or more of data. The Data-Driven Approaches to Searches for the Technosignatures of Advanced Civilizations_ study at the W. M. Keck Institute for Space Studies was intended to revisit searches for evidence of alien technologies in light of these developments. Data-driven searches, being able to process volumes of data much greater than a human could, and in a reproducible manner, can identify *anomalies* that could be clues to the presence of technosignatures. A key outcome of this workshop was that technosignature searches should be conducted in a manner consistent with Freeman Dyson's "First Law of SETI Investigations," namely "every search for alien civilizations should be planned to give interesting results even when no aliens are discovered." This approach to technosignatures is commensurate with NASA's approach to biosignatures in that no single observation or measurement can be taken as providing full certainty for the detection of life. Areas of particular promise identified during the workshop were (*) Data Mining of Large Sky Surveys, (*) All-Sky Survey at Far-Infrared Wavelengths, (*) Surveys with Radio Astronomical Interferometers, and (*) Artifacts in the Solar System. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: Final Report prepared for the W. M. Keck Institute for Space Studies (KISS), http://kiss.caltech.edu/workshops/technosignatures/technosignatures.html ; eds. Lazio, Djorgovski, Howard, & Cutler; The study leads gratefully acknowledge the outstanding support of Michele Judd, KISS Executive Director, and her dedicated staff, who made the study experience invigorating and enormously productive

Tuning of Nuclear Spectroscopic Telescope Array (NuSTAR) Application Specific Integrated Circuits (ASICs) to improve low energy threshold of future Hard X-ray Imaging Detectors

Authors: Daniel P. Violette, Branden Allen, Jaesub Hong, Hiromasa Miyasaka, Jonathan Grindlay

Abstract: Detector commanding, processing and readout of spaceborne instrumentation is often accomplished with Application Specific Integrated Circuits (ASICs). The ASIC designed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission (NuASIC) enables future tiled CdZnTe (CZT) detector array readout for x-ray detectors such as the High Resolution Energetic X-ray Imager (HREXI). Modified NuASIC gain s… ▽ More Detector commanding, processing and readout of spaceborne instrumentation is often accomplished with Application Specific Integrated Circuits (ASICs). The ASIC designed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission (NuASIC) enables future tiled CdZnTe (CZT) detector array readout for x-ray detectors such as the High Resolution Energetic X-ray Imager (HREXI). Modified NuASIC gain settings have been implemented for HREXI's broader targeted imaging energy range (3-300 keV) compared to NuSTAR (2-79 keV), which may require updated NuASIC internal parameters for optimal energy resolution. To reach HREXI's targeted low energy threshold, we have also enabled the NuASIC's "Charge Pump Mode" (CPM), which introduces an additional tuning parameter. In this paper, we describe the mechanics of the NuASIC's adjustable parameters and use our recently developed ASIC Test Stand (ATS) to probe a "bare" NuASIC using its internal test pulser. We record the effects of parameter tuning on the device's electronics noise and low energy threshold and report the optimal set of parameters for HREXI's updated gain setting. We detail a semi-automated procedure to derive the optimal parameters for each of HREXI's large area, closely tiled NuASIC/CZT detectors to expedite instrument integration. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: 20 pages, 11 figures, published in JATIS

Journal ref: JATIS Vol. 9, Issue 01, 016001 (January 2023)

Answers to frequently asked questions about the pulsar timing array Hellings and Downs curve

Authors: Joseph D. Romano, Bruce Allen

Abstract: We answer frequently asked questions (FAQs) about the Hellings and Downs correlation curve -- the "smoking-gun" signature that pulsar timing arrays (PTAs) have detected gravitational waves (GWs). Many of these questions arise from inadvertently applying intuition about the effects of GWs on LIGO-like detectors to the case of pulsar timing, where not all of that intuition applies. This is because E… ▽ More We answer frequently asked questions (FAQs) about the Hellings and Downs correlation curve -- the "smoking-gun" signature that pulsar timing arrays (PTAs) have detected gravitational waves (GWs). Many of these questions arise from inadvertently applying intuition about the effects of GWs on LIGO-like detectors to the case of pulsar timing, where not all of that intuition applies. This is because Earth-based detectors, like LIGO and Virgo, have arms that are short (km scale) compared to the wavelengths of the GWs that they detect (approximately 100-10,000 km). In contrast, PTAs respond to GWs whose wavelengths (tens of light-years) are much shorter than their arms (a typical PTA pulsar is hundreds to thousands of light-years from Earth). To demonstrate this, we calculate the time delay induced by a passing GW along an Earth-pulsar baseline (a "one-arm, one-way" detector) and compare it in the "short-arm" (LIGO-like) and "long-arm" (PTA) limits. This provides qualitative and quantitative answers to many questions about the Hellings and Downs curve. The resulting FAQ sheet should help in understanding the "evidence for GWs" recently announced by several PTA collaborations. △ Less

Submitted 29 January, 2024; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: 30 pages, 12 figures; substantially expanded compared to the previous version

Symmetry in models of natural selection

Authors: Benjamin Allen

Abstract: Symmetry arguments are frequently used -- often implicitly -- in mathematical modeling of natural selection. Symmetry simplifies the analysis of models and reduces the number of distinct population states to be considered. Here, I introduce a formal definition of symmetry in mathematical models of natural selection. This definition applies to a broad class of models that satisfy a minimal set of a… ▽ More Symmetry arguments are frequently used -- often implicitly -- in mathematical modeling of natural selection. Symmetry simplifies the analysis of models and reduces the number of distinct population states to be considered. Here, I introduce a formal definition of symmetry in mathematical models of natural selection. This definition applies to a broad class of models that satisfy a minimal set of assumptions, using a framework developed in previous works. In this framework, population structure is represented by a set of sites at which alleles can live, and transitions occur via replacement of some alleles by copies of others. A symmetry is defined as a permutation of sites that preserves probabilities of replacement and mutation. The symmetries of a given selection process form a group, which acts on population states in a way that preserves the Markov chain representing selection. Applying classical results on group actions, I formally characterize the use of symmetry to reduce the states of this Markov chain, and obtain bounds on the number of states in the reduced chain. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 21 pages, 4 figures

MSC Class: 92D15

Research Protocol for the Google Health Digital Well-being Study

Authors: Daniel McDuff, Andrew Barakat, Ari Winbush, Allen Jiang, Felicia Cordeiro, Ryann Crowley, Lauren E. Kahn, John Hernandez, Nicholas B. Allen

Abstract: The impact of digital device use on health and well-being is a pressing question to which individuals, families, schools, policy makers, legislators, and digital designers are all demanding answers. However, the scientific literature on this topic to date is marred by small and/or unrepresentative samples, poor measurement of core constructs (e.g., device use, smartphone addiction), and a limited… ▽ More The impact of digital device use on health and well-being is a pressing question to which individuals, families, schools, policy makers, legislators, and digital designers are all demanding answers. However, the scientific literature on this topic to date is marred by small and/or unrepresentative samples, poor measurement of core constructs (e.g., device use, smartphone addiction), and a limited ability to address the psychological and behavioral mechanisms that may underlie the relationships between device use and well-being. A number of recent authoritative reviews have made urgent calls for future research projects to address these limitations. The critical role of research is to identify which patterns of use are associated with benefits versus risks, and who is more vulnerable to harmful versus beneficial outcomes, so that we can pursue evidence-based product design, education, and regulation aimed at maximizing benefits and minimizing risks of smartphones and other digital devices. We describe a protocol for a Digital Well-Being (DWB) study to help answer these questions. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Identifying and Consolidating Knowledge Engineering Requirements

Authors: Bradley P. Allen, Filip Ilievski, Saurav Joshi

Abstract: Knowledge engineering is the process of creating and maintaining knowledge-producing systems. Throughout the history of computer science and AI, knowledge engineering workflows have been widely used because high-quality knowledge is assumed to be crucial for reliable intelligent agents. However, the landscape of knowledge engineering has changed, presenting four challenges: unaddressed stakeholder… ▽ More Knowledge engineering is the process of creating and maintaining knowledge-producing systems. Throughout the history of computer science and AI, knowledge engineering workflows have been widely used because high-quality knowledge is assumed to be crucial for reliable intelligent agents. However, the landscape of knowledge engineering has changed, presenting four challenges: unaddressed stakeholder requirements, mismatched technologies, adoption barriers for new organizations, and misalignment with software engineering practices. In this paper, we propose to address these challenges by developing a reference architecture using a mainstream software methodology. By studying the requirements of different stakeholders and eras, we identify 23 essential quality attributes for evaluating reference architectures. We assess three candidate architectures from recent literature based on these attributes. Finally, we discuss the next steps towards a comprehensive reference architecture, including prioritizing quality attributes, integrating components with complementary strengths, and supporting missing socio-technical requirements. As this endeavor requires a collaborative effort, we invite all knowledge engineering researchers and practitioners to join us. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Null Distance and Gromov-Hausdorff Convergence of Warped Product Spacetimes

Authors: Brian Allen

Abstract: What is the analogous notion of Gromov-Hausdorff convergence for sequences of spacetimes? Since a Lorentzian manifold is not inherently a metric space, one cannot simply use the traditional definition. One approach offered by Sormani and Vega \cite{SV} is to define a metric space structure on a spacetime by means of the null distance. Then one can define convergence of spacetimes using the usual d… ▽ More What is the analogous notion of Gromov-Hausdorff convergence for sequences of spacetimes? Since a Lorentzian manifold is not inherently a metric space, one cannot simply use the traditional definition. One approach offered by Sormani and Vega \cite{SV} is to define a metric space structure on a spacetime by means of the null distance. Then one can define convergence of spacetimes using the usual definition of Gromov-Hausdorff convergence. In this paper we explore this approach by giving many examples of sequences of warped product spacetimes with the null distance converging in the Gromov-Hausdorff sense. In addition, we give an optimal convergence theorem which shows that under natural geometric hypotheses a sequence of warped product spacetimes converge to a specific limiting warped product spacetime. The examples given further serve to show that the hypotheses of this convergence theorem are optimal. △ Less

Submitted 16 October, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: v2: 32 pages, 1 figure, To appear in General Relativity and Gravitation

On the Stability of Llarull's Theorem in Dimension Three

Authors: Brian Allen, Edward Bryden, Demetre Kazaras

Abstract: Llarull's Theorem states that any Riemannian metric on the $n$-sphere which has scalar curv{\-}ature greater than or equal to $n(n-1)$, and whose distance function is bounded below by the unit sphere's, is isometric to the unit sphere. Gromov later posed the {\emph{Spherical Stability Problem}}, which probes the flexibility of this fact. We give a resolution to this problem in dimension $3$. Infor… ▽ More Llarull's Theorem states that any Riemannian metric on the $n$-sphere which has scalar curv{\-}ature greater than or equal to $n(n-1)$, and whose distance function is bounded below by the unit sphere's, is isometric to the unit sphere. Gromov later posed the {\emph{Spherical Stability Problem}}, which probes the flexibility of this fact. We give a resolution to this problem in dimension $3$. Informally, the main result asserts that a sequence of Riemannian $3$-spheres whose distance functions are bounded below by the unit sphere's with uniformly bounded Cheeger isoperimetric constant and scalar curvatures tending to $6$ must approach the round $3$-sphere in the volume preserving Sormani-Wenger Intrinsic Flat sense. The argument is based on a proof of Llarull's Theorem due to Hirsch-Kazaras-Khuri-Zhang using spacetime harmonic functions. △ Less

Submitted 24 November, 2023; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: v2: Minor changes made to the exposition of the introduction, 22 pages, 3 figures, comments welcome

MSC Class: 53

The International Pulsar Timing Array checklist for the detection of nanohertz gravitational waves

Authors: Bruce Allen, Sanjeev Dhurandhar, Yashwant Gupta, Maura McLaughlin, Priyamvada Natarajan, Ryan M. Shannon, Eric Thrane, Alberto Vecchio

Abstract: Pulsar timing arrays (PTAs) provide a way to detect gravitational waves at nanohertz frequencies. In this band, the most likely signals are stochastic, with a power spectrum that rises steeply at lower frequencies. Indeed, the observation of a common red noise process in pulsar-timing data suggests that the first credible detection of nanohertz-frequency gravitational waves could take place within… ▽ More Pulsar timing arrays (PTAs) provide a way to detect gravitational waves at nanohertz frequencies. In this band, the most likely signals are stochastic, with a power spectrum that rises steeply at lower frequencies. Indeed, the observation of a common red noise process in pulsar-timing data suggests that the first credible detection of nanohertz-frequency gravitational waves could take place within the next few years. The detection process is complicated by the nature of the signals and the noise: the first observational claims will be statistical inferences drawn at the threshold of detectability. To demonstrate that gravitational waves are creating some of the noise in the pulsar-timing data sets, observations must exhibit the Hellings and Downs curve -- the angular correlation function associated with gravitational waves -- as well as demonstrating that there are no other reasonable explanations. To ensure that detection claims are credible, the International Pulsar Timing Array (IPTA) has a formal process to vet results prior to publication. This includes internal sharing of data and processing pipelines between different PTAs, enabling independent cross-checks and validation of results. To oversee and validate any detection claim, the IPTA has also created an eight-member Detection Committee (DC) which includes four independent external members. IPTA members will only publish their results after a formal review process has concluded. This document is the initial DC checklist, describing some of the conditions that should be fulfilled by a credible detection. At the present time none of the PTAs have a detection claim; therefore this document serves as a road map for the future. △ Less

Submitted 31 May, 2023; v1 submitted 9 April, 2023; originally announced April 2023.

Comments: 6 pages

Imaging effects due to pixel distortions in CdZnTe (CZT) detectors -- results from the HREXI Calibration Facility (HCF)

Authors: Arkadip Basak, Branden Allen, Jaesub Hong, Daniel P. Violette, Jonathan Grindlay

Abstract: ProtoEXIST2 (P2) was a prototype imaging X-ray detector plane developed for wide-field Time Domain Astrophysics (TDA) in the 5 - 200 keV energy band. It was composed of an 8 $\times$ 8 array of 5 mm thick, 2cm $\times$ 2cm pixelated (32 $\times$ 32) CdZnTe (CZT) detectors with a 0.6 mm pitch that utilize the NuSTAR ASIC(NuASIC) for readout. During the initial detector development process leading u… ▽ More ProtoEXIST2 (P2) was a prototype imaging X-ray detector plane developed for wide-field Time Domain Astrophysics (TDA) in the 5 - 200 keV energy band. It was composed of an 8 $\times$ 8 array of 5 mm thick, 2cm $\times$ 2cm pixelated (32 $\times$ 32) CdZnTe (CZT) detectors with a 0.6 mm pitch that utilize the NuSTAR ASIC(NuASIC) for readout. During the initial detector development process leading up to post-flight examination of the entire detector plane, distortions in expected pixel positions and shapes were observed in a significant fraction of the detectors. The HREXI (High Resolution Energetic X-ray Imager) Calibration Facility (HCF) was designed and commissioned to improve upon these early experiments and to rapidly map out and characterize pixel non-uniformities and defects within CZT detector planes at resolutions down to 50 $\rm μ$m. Using this facility, the sub-pixel level detector response of P2 was measured at 100 $\rm μ$m resolution and analyzed to extract and evaluate the area and profile of individual pixels, their morphology across the entire P2 detector plane for comparison with previous measurements and to provide additional characterization. In this article, we evaluate the imaging performance of a coded-aperture telescope using the observed pixel morphology for P2 detectors. This investigation will serve as an initial guide for detector selection in the development of HREXI detector planes, for the future implementation of the 4pi X-Ray Imaging Observatory (4piXIO) mission which aims to provide simultaneous and continuous imaging of the full sky ($\rm 4π$ sr) in the 3-200 keV energy band with $\rm \simeq$ 2 arcmin angular resolution and $\simeq$ 10 arcsec source localization, as well as other, future coded-aperture instruments. △ Less

Submitted 17 March, 2023; originally announced March 2023.

Journal ref: J. Astron. Telesc. Instrum. Syst. 9(1), 018005 (2023)

Deep Einstein@Home all-sky search for continuous gravitational waves in LIGO O3 public data

Authors: B. Steltner, M. A. Papa, H. -B. Eggenstein, R. Prix, M. Bensch, B. Allen, B. Machenschalk

Abstract: We present the results of an all-sky search for continuous gravitational waves in the public LIGO O3 data. The search covers signal frequencies $20$ Hz $\leq f \leq 800$ Hz and a spin-down range down to $-2.6\times 10^{-9}$ Hz s$^{-1}$1, motivated by detectability studies on synthetic populations of Galactic neutron stars. This search is the most sensitive all-sky search to date in this frequency/… ▽ More We present the results of an all-sky search for continuous gravitational waves in the public LIGO O3 data. The search covers signal frequencies $20$ Hz $\leq f \leq 800$ Hz and a spin-down range down to $-2.6\times 10^{-9}$ Hz s$^{-1}$1, motivated by detectability studies on synthetic populations of Galactic neutron stars. This search is the most sensitive all-sky search to date in this frequency/spin-down region. The initial search was performed using the first half of the public LIGO O3 data (O3a), utilizing Graphical Processing Units provided in equal parts by the volunteers of the Einstein@Home computing project and by the ATLAS cluster. After a hierarchical follow-up in seven stages, 12 candidates remain. Six are discarded at the eighth stage, by using the remaining O3 LIGO data (O3b). The surviving six can be ascribed to continuous-wave fake signals present in the LIGO data for validation purposes. We recover these fake signals with very high accuracy with our last stage search, which coherently combines all O3 data. Based on our results, we set upper limits on the gravitational wave amplitude $h_0$, and translate these in upper limits on the neutron star ellipticity and on the $r$-mode amplitude. The most stringent upper limits are at $203$ Hz, with $h_0=8.1 \times 10^{-26}$ at the 90% confidence level. Our results exclude neutron stars rotating faster than $5$ ms with ellipticities greater than $5\times 10^{-8} \left[{d\over{100~\textrm{pc}}}\right]$ within a distance $d$ from Earth and $r$-mode amplitudes $α\geq 10^{-5} \left[{d\over{100~\textrm{pc}}}\right]$ for neutron stars spinning faster than $150$ Hz. △ Less

Submitted 19 June, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

Comments: Accepted for publication in The Astrophysical Journal on 31 May 2023. 13 pages, 10 figures, 3 tables

Journal ref: ApJ 952 55 (2023)

Nonlinear social evolution and the emergence of collective action

Authors: Benjamin Allen, Abdur-Rahman Khwaja, James L. Donahue, Cassidy Lattanzio, Yulia A. Dementieva, Christine Sample

Abstract: Organisms from microbes to humans engage in a variety of social behaviors, which affect fitness in complex, often nonlinear ways. The question of how these behaviors evolve has consequences ranging from antibiotic resistance to human origins. However, evolution with nonlinear social interactions is challenging to model mathematically, especially in combination with spatial, group, and/or kin assor… ▽ More Organisms from microbes to humans engage in a variety of social behaviors, which affect fitness in complex, often nonlinear ways. The question of how these behaviors evolve has consequences ranging from antibiotic resistance to human origins. However, evolution with nonlinear social interactions is challenging to model mathematically, especially in combination with spatial, group, and/or kin assortment. We derive a mathematical condition for natural selection with synergistic interactions among any number of individuals. This result applies to populations with arbitrary (but fixed) spatial or network structure, group subdivision, and/or mating patterns. In this condition, nonlinear fitness effects are ascribed to collectives, and weighted by a new measure of collective relatedness. For weak selection, this condition can be systematically evaluated by computing branch lengths of ancestral trees. We apply this condition to pairwise games between diploid relatives, and to dilemmas of collective help or harm among siblings and on spatial networks. Our work provides a rigorous basis for extending the notion of ``actor", in the study of social evolution, from individuals to collectives. △ Less

Submitted 13 June, 2024; v1 submitted 28 February, 2023; originally announced February 2023.

Comments: 12 page main text, with 72 page Supplementary Information. Published in PNAS Nexus. Previously posted as a preprint entitled "Natural selection for collective action"

MSC Class: 92D15; 91A22

Journal ref: (2024). PNAS nexus, 3(4), pgae131

The Voronoi Region of the Barnes-Wall Lattice $Λ_{16}$

Authors: Daniel Pook-Kolb, Erik Agrell, Bruce Allen

Abstract: We give a detailed description of the Voronoi region of the Barnes-Wall lattice $Λ_{16}$, including its vertices, relevant vectors, and symmetry group. The exact value of its quantizer constant is calculated, which was previously only known approximately. To verify the result, we estimate the same constant numerically and propose a new very simple method to quantify the variance of such estimates,… ▽ More We give a detailed description of the Voronoi region of the Barnes-Wall lattice $Λ_{16}$, including its vertices, relevant vectors, and symmetry group. The exact value of its quantizer constant is calculated, which was previously only known approximately. To verify the result, we estimate the same constant numerically and propose a new very simple method to quantify the variance of such estimates, which is far more accurate than the commonly used jackknife estimator. △ Less

Submitted 15 January, 2023; originally announced January 2023.

Comments: 8 pages, 2 figures. This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

MSC Class: 11H06; 11H56; 52B20; 52C07

Journal ref: IEEE Journal on Selected Areas in Information Theory, vol. 4, pp. 16-23, 2023

A deep learning approach to using wearable seismocardiography (SCG) for diagnosing aortic valve stenosis and predicting aortic hemodynamics obtained by 4D flow MRI

Authors: Mahmoud E. Khani, Ethan M. I. Johnson, Aparna Sodhi, Joshua Robinson, Cynthia K. Rigsby, Bradly D. Allen, Michael Markl

Abstract: In this paper, we explored the use of deep learning for the prediction of aortic flow metrics obtained using 4D flow MRI using wearable seismocardiography (SCG) devices. 4D flow MRI provides a comprehensive assessment of cardiovascular hemodynamics, but it is costly and time-consuming. We hypothesized that deep learning could be used to identify pathological changes in blood flow, such as elevated… ▽ More In this paper, we explored the use of deep learning for the prediction of aortic flow metrics obtained using 4D flow MRI using wearable seismocardiography (SCG) devices. 4D flow MRI provides a comprehensive assessment of cardiovascular hemodynamics, but it is costly and time-consuming. We hypothesized that deep learning could be used to identify pathological changes in blood flow, such as elevated peak systolic velocity Vmax in patients with heart valve diseases, from SCG signals. We also investigated the ability of this deep learning technique to differentiate between patients diagnosed with aortic valve stenosis (AS), non-AS patients with a bicuspid aortic valve (BAV), non-AS patients with a mechanical aortic valve (MAV), and healthy subjects with a normal tricuspid aortic valve (TAV). In a study of 77 subjects who underwent same-day 4D flow MRI and SCG, we found that the Vmax values obtained using deep learning and SCGs were in good agreement with those obtained by 4D flow MRI. Additionally, subjects with TAV, BAV, MAV, and AS could be classified with ROC-AUC values of 92%, 95%, 81%, and 83%, respectively. This suggests that SCG obtained using low-cost wearable electronics may be used as a supplement to 4D flow MRI exams or as a screening tool for aortic valve disease. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 16 pages, 4 figures

Spinning up a Daze: TESS Uncovers a Hot Jupiter orbiting the Rapid-Rotator TOI-778

Authors: Jake Clark, Brett Addison, Jack Okumura, Sydney Vach, Alexis Heitzmann, Joseph Rodriguez, Duncan Wright, Mathieu Clerte, Carolyn Brown, Tara Fetherolf, Robert Wittenmyer, Peter Plavchan, Stephen Kane, Jonathan Horner, John Kielkopf, Avi Shporer, C. Tinney, Liu Hui-Gen, Sarah Ballard, Brendan Bowler, Matthew Mengel, George Zhou, Annette Lee, Avelyn David, Jessica Heim , et al. (46 additional authors not shown)

Abstract: NASA's Transiting Exoplanet Survey Satellite (TESS) mission, has been uncovering a growing number of exoplanets orbiting nearby, bright stars. Most exoplanets that have been discovered by TESS orbit narrow-line, slow-rotating stars, facilitating the confirmation and mass determination of these worlds. We present the discovery of a hot Jupiter orbiting a rapidly rotating ($v\sin{(i)}= 35.1\pm1.0$km… ▽ More NASA's Transiting Exoplanet Survey Satellite (TESS) mission, has been uncovering a growing number of exoplanets orbiting nearby, bright stars. Most exoplanets that have been discovered by TESS orbit narrow-line, slow-rotating stars, facilitating the confirmation and mass determination of these worlds. We present the discovery of a hot Jupiter orbiting a rapidly rotating ($v\sin{(i)}= 35.1\pm1.0$km/s) early F3V-dwarf, HD115447 (TOI-778). The transit signal taken from Sectors 10 and 37 of TESS's initial detection of the exoplanet is combined with follow-up ground-based photometry and velocity measurements taken from Minerva-Australis, TRES, CORALIE and CHIRON to confirm and characterise TOI-778b. A joint analysis of the light curves and the radial velocity measurements yield a mass, radius, and orbital period for TOI-778b of $2.76^{+0.24}_{-0.23}$Mjup, $1.370\pm0.043$Rjup and $\sim4.63$ days, respectively. The planet orbits a bright ($V = 9.1$mag) F3-dwarf with $M=1.40\pm0.05$Msun, $R=1.70\pm0.05$Rsun, and $\log g=4.05\pm0.17$. We observed a spectroscopic transit of TOI-778b, which allowed us to derive a sky-projected spin-orbit angle of $18^{\circ}\pm11^{\circ}$, consistent with an aligned planetary system. This discovery demonstrates the capability of smaller aperture telescopes such as Minerva-Australis to detect the radial velocity signals produced by planets orbiting broad-line, rapidly rotating stars. △ Less

Submitted 30 April, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 18 pages, 9 figures, and 4 tables. Submitted to the Astronomical Journal

Journal ref: AJ 165 207 (2023)

Exact calculation of quantizer constants for arbitrary lattices

Authors: Daniel Pook-Kolb, Bruce Allen, Erik Agrell

Abstract: We present an algorithm for the exact computer-aided construction of the Voronoi cells of lattices with known symmetry group. Our algorithm scales better than linearly with the total number of faces and is applicable to dimensions beyond 12, which previous methods could not achieve. The new algorithm is applied to the Coxeter-Todd lattice $K_{12}$ as well as to a family of lattices obtained from l… ▽ More We present an algorithm for the exact computer-aided construction of the Voronoi cells of lattices with known symmetry group. Our algorithm scales better than linearly with the total number of faces and is applicable to dimensions beyond 12, which previous methods could not achieve. The new algorithm is applied to the Coxeter-Todd lattice $K_{12}$ as well as to a family of lattices obtained from laminating $K_{12}$. By optimizing this family, we obtain a new best 13-dimensional lattice quantizer (among the lattices with published exact quantizer constants). △ Less

Submitted 18 February, 2024; v1 submitted 14 October, 2022; originally announced November 2022.

Comments: 24 pages, 7 figures

MSC Class: 11H06; 11H56; 52B20; 52C07

Stability of the positive mass theorem and torus rigidity theorems under integral curvature bounds

Authors: Brian Allen, Edward Bryden, Demetre Kazaras

Abstract: Work of D. Stern and Bray-Kazaras-Khuri-Stern provide differential-geometric identities which relate the scalar curvature of Riemannian 3-manifolds to global invariants in terms of harmonic functions. These quantitative formulas are useful for stability results and show promise for more applications of this type. In this paper, we analyze harmonic maps to flat model spaces in order to address conj… ▽ More Work of D. Stern and Bray-Kazaras-Khuri-Stern provide differential-geometric identities which relate the scalar curvature of Riemannian 3-manifolds to global invariants in terms of harmonic functions. These quantitative formulas are useful for stability results and show promise for more applications of this type. In this paper, we analyze harmonic maps to flat model spaces in order to address conjectures concerning the geometric stability of the positive mass theorem and the Geroch conjecture. By imposing integral Ricci curvature and isoperimetric bounds, we leverage the previously mentioned formulas to establish strong control on these harmonic maps. When the mass of an asymptotically flat manifold is sufficiently small or when a Riemannian torus has almost non-negative scalar curvature, we upgrade the maps to diffeomorphisms and give quantitative Hölder closeness to the model spaces. △ Less

Submitted 9 October, 2022; originally announced October 2022.

Variational quantum iterative power algorithms for global optimization

Authors: Thi Ha Kyaw, Micheline B. Soley, Brandon Allen, Paul Bergold, Chong Sun, Victor S. Batista, Alán Aspuru-Guzik

Abstract: We introduce a family of variational quantum algorithms called quantum iterative power algorithms (QIPA) that outperform existing hybrid near-term quantum algorithms of the same kind. We demonstrate the capabilities of QIPA as applied to three different global-optimization numerical experiments: the ground-state optimization of the $H_2$ molecular dissociation, search of the transmon qubit ground-… ▽ More We introduce a family of variational quantum algorithms called quantum iterative power algorithms (QIPA) that outperform existing hybrid near-term quantum algorithms of the same kind. We demonstrate the capabilities of QIPA as applied to three different global-optimization numerical experiments: the ground-state optimization of the $H_2$ molecular dissociation, search of the transmon qubit ground-state, and biprime factorization. Since our algorithm is hybrid, quantum/classical technologies such as error mitigation and adaptive variational ansatzes can easily be incorporated into the algorithm. Due to the shallow quantum circuit requirements, we anticipate large-scale implementation and adoption of the proposed algorithm across current major quantum hardware. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 17 pages, 7 figures

Hellings and Downs correlation of an arbitrary set of pulsars

Authors: Bruce Allen, Joseph D. Romano

Abstract: Pulsar timing arrays (PTAs) detect gravitational waves (GWs) via the correlations they induce in the arrival times of pulses from different pulsars. We assume that the GWs are described by a Gaussian ensemble. The mean correlation $h^2 μ_{\rm u}(γ)$ as a function of the angle $γ$ between the directions to two pulsars was predicted by Hellings and Downs (HD) in 1983. The variance… ▽ More Pulsar timing arrays (PTAs) detect gravitational waves (GWs) via the correlations they induce in the arrival times of pulses from different pulsars. We assume that the GWs are described by a Gaussian ensemble. The mean correlation $h^2 μ_{\rm u}(γ)$ as a function of the angle $γ$ between the directions to two pulsars was predicted by Hellings and Downs (HD) in 1983. The variance $σ^2_{\rm tot}(γ)$ in this correlation was recently calculated by Allen[11] for a single noise-free pulsar pair at angle $γ$, which shows that after averaging over many pairs, the variance reduces to an intrinsic cosmic variance $σ^2_{\rm cos}(γ)$. Here, we extend this to an $arbitrary$ set of pulsars at specific sky locations, with pulsar pairs binned by $γ$. We derive the linear combination of pulsar-pair correlations which is the optimal estimator of the HD correlation for each bin, illustrating our methods with plots of the expected range of variation away from the HD curve, for the sets of pulsars monitored by three active PTA collaborations. We compute the variance of and the covariance between these binned estimates, and show that these reduce to the cosmic variance and covariance $s(γ,γ')$ respectively, in the many-pulsar limit. The likely fluctuations away from the HD curve $μ_{\rm u}(γ)$ are strongly correlated/anticorrelated in the three angular regions where $μ_{\rm u}(γ)$ is successively positive, negative, and positive. We also construct the optimal estimator of the squared strain $h^2$. When there are very many pulsar pairs, this determines $h^2$ with arbitrary precision because PTAs probe an infinite set of GW modes. To assess observed deviations away from the HD curve, we characterize several $χ^2$ goodness-of-fit statistics. We also show how pulsar noise and measurement noise can be included. △ Less

Submitted 24 August, 2023; v1 submitted 15 August, 2022; originally announced August 2022.

Comments: Final published version. Note that the arXiv version of the Abstract has been shortened to fit arXiv requirements

Journal ref: Phys. Rev. D 108, 043026 (2023)

Deep Learning for Classification of Thyroid Nodules on Ultrasound: Validation on an Independent Dataset

Authors: Jingxi Weng, Benjamin Wildman-Tobriner, Mateusz Buda, Jichen Yang, Lisa M. Ho, Brian C. Allen, Wendy L. Ehieli, Chad M. Miller, Jikai Zhang, Maciej A. Mazurowski

Abstract: Objectives: The purpose is to apply a previously validated deep learning algorithm to a new thyroid nodule ultrasound image dataset and compare its performances with radiologists. Methods: Prior study presented an algorithm which is able to detect thyroid nodules and then make malignancy classifications with two ultrasound images. A multi-task deep convolutional neural network was trained from 127… ▽ More Objectives: The purpose is to apply a previously validated deep learning algorithm to a new thyroid nodule ultrasound image dataset and compare its performances with radiologists. Methods: Prior study presented an algorithm which is able to detect thyroid nodules and then make malignancy classifications with two ultrasound images. A multi-task deep convolutional neural network was trained from 1278 nodules and originally tested with 99 separate nodules. The results were comparable with that of radiologists. The algorithm was further tested with 378 nodules imaged with ultrasound machines from different manufacturers and product types than the training cases. Four experienced radiologists were requested to evaluate the nodules for comparison with deep learning. Results: The Area Under Curve (AUC) of the deep learning algorithm and four radiologists were calculated with parametric, binormal estimation. For the deep learning algorithm, the AUC was 0.69 (95% CI: 0.64 - 0.75). The AUC of radiologists were 0.63 (95% CI: 0.59 - 0.67), 0.66 (95% CI:0.61 - 0.71), 0.65 (95% CI: 0.60 - 0.70), and 0.63 (95%CI: 0.58 - 0.67). Conclusion: In the new testing dataset, the deep learning algorithm achieved similar performances with all four radiologists. The relative performance difference between the algorithm and the radiologists is not significantly affected by the difference of ultrasound scanner. △ Less

Submitted 4 May, 2023; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: Clinical Imaging (2023)

Simulating spin dynamics with quantum computers

Authors: Jarrett L. Lancaster, D. Brysen Allen

Abstract: IBM quantum computers are used to simulate the dynamics of small systems of interacting quantum spins. For time-independent systems with $N\leq 2$ spins, we construct circuits which compute the exact time evolution at arbitrary times and allow measurement of spin expectation values. Larger systems with spins require a the Lie-Trotter decomposition of the time evolution operator, and we investigate… ▽ More IBM quantum computers are used to simulate the dynamics of small systems of interacting quantum spins. For time-independent systems with $N\leq 2$ spins, we construct circuits which compute the exact time evolution at arbitrary times and allow measurement of spin expectation values. Larger systems with spins require a the Lie-Trotter decomposition of the time evolution operator, and we investigate the case of $N=3$ explicitly. Basic readout-error mitigation is introduced to improve the quality of results run on these noisy devices. The simulations provide an interesting experimental component to the standard treatment of quantum spin in an undergraduate quantum mechanics course. △ Less

Submitted 29 June, 2022; originally announced July 2022.

Comments: 9 pages, 7 figures

The coalescent in finite populations with arbitrary, fixed structure

Authors: Benjamin Allen, Alex McAvoy

Abstract: The coalescent is a stochastic process representing ancestral lineages in a population undergoing neutral genetic drift. Originally defined for a well-mixed population, the coalescent has been adapted in various ways to accommodate spatial, age, and class structure, along with other features of real-world populations. To further extend the range of population structures to which coalescent theory… ▽ More The coalescent is a stochastic process representing ancestral lineages in a population undergoing neutral genetic drift. Originally defined for a well-mixed population, the coalescent has been adapted in various ways to accommodate spatial, age, and class structure, along with other features of real-world populations. To further extend the range of population structures to which coalescent theory applies, we formulate a coalescent process for a broad class of neutral drift models with arbitrary -- but fixed -- spatial, age, sex, and class structure, haploid or diploid genetics, and any fixed mating pattern. Here, the coalescent is represented as a random sequence of mappings $\mathcal{C} = \left(C_t\right)_{t=0}^\infty$ from a finite set $G$ to itself. The set $G$ represents the ``sites'' (in individuals, in particular locations and/or classes) at which these alleles can live. The state of the coalescent, $C_t:G \to G$, maps each site $g \in G$ to the site containing $g$'s ancestor, $t$ time-steps into the past. Using this representation, we define and analyze coalescence time, coalescence branch length, mutations prior to coalescence, and stationary probabilities of identity-by-descent and identity-by-state. For low mutation, we provide a recipe for computing identity-by-descent and identity-by-state probabilities via the coalescent. Applying our results to a diploid population with arbitrary sex ratio $r$, we find that measures of genetic dissimilarity, among any set of sites, are scaled by $4r(1-r)$ relative to the even sex ratio case. △ Less

Submitted 28 June, 2024; v1 submitted 6 July, 2022; originally announced July 2022.

Comments: 71 pages, 2 figures

MSC Class: 92D15

Journal ref: Theoretical Population Biology 158: 150-169 (2024)

Artificial Intelligence Techniques for Next-Generation Mega Satellite Networks

Authors: Bassel Al Homssi, Kosta Dakic, Ke Wang, Tansu Alpcan, Ben Allen, Russell Boyce, Sithamparanathan Kandeepan, Akram Al-Hourani, Walid Saad

Abstract: Space communications, particularly massive satellite networks, re-emerged as an appealing candidate for next generation networks due to major advances in space launching, electronics, processing power, and miniaturization. However, massive satellite networks rely on numerous underlying and intertwined processes that cannot be truly captured using conventionally used models, due to their dynamic an… ▽ More Space communications, particularly massive satellite networks, re-emerged as an appealing candidate for next generation networks due to major advances in space launching, electronics, processing power, and miniaturization. However, massive satellite networks rely on numerous underlying and intertwined processes that cannot be truly captured using conventionally used models, due to their dynamic and unique features such as orbital speed, inter-satellite links, short pass time, and satellite footprint, among others. Hence, new approaches are needed to enable the network to proactively adjust to the rapidly varying conditions associated within the link. Artificial intelligence (AI) provides a pathway to capture these processes, analyze their behavior, and model their effect on the network. This article introduces the application of AI techniques for integrated terrestrial satellite networks, particularly massive satellite network communications. It details the unique features of massive satellite networks, and the overarching challenges concomitant with their integration into the current communication infrastructure. Moreover, this article provides insights into state-of-the-art AI techniques across various layers of the communication link. This entails applying AI for forecasting the highly dynamic radio channel, spectrum sensing and classification, signal detection and demodulation, inter-satellite and satellite access network optimization, and network security. Moreover, future paradigms and the mapping of these mechanisms onto practical networks are outlined. △ Less

Submitted 16 September, 2023; v1 submitted 2 June, 2022; originally announced July 2022.

Multistep Automated Data Labelling Procedure (MADLaP) for Thyroid Nodules on Ultrasound: An Artificial Intelligence Approach for Automating Image Annotation

Authors: Jikai Zhang, Maciej M. Mazurowski, Brian C. Allen, Benjamin Wildman-Torbiner

Abstract: Machine learning (ML) for diagnosis of thyroid nodules on ultrasound is an active area of research. However, ML tools require large, well-labelled datasets, the curation of which is time-consuming and labor-intensive. The purpose of our study was to develop and test a deep-learning-based tool to facilitate and automate the data annotation process for thyroid nodules; we named our tool Multistep Au… ▽ More Machine learning (ML) for diagnosis of thyroid nodules on ultrasound is an active area of research. However, ML tools require large, well-labelled datasets, the curation of which is time-consuming and labor-intensive. The purpose of our study was to develop and test a deep-learning-based tool to facilitate and automate the data annotation process for thyroid nodules; we named our tool Multistep Automated Data Labelling Procedure (MADLaP). MADLaP was designed to take multiple inputs included pathology reports, ultrasound images, and radiology reports. Using multiple step-wise modules including rule-based natural language processing, deep-learning-based imaging segmentation, and optical character recognition, MADLaP automatically identified images of a specific thyroid nodule and correctly assigned a pathology label. The model was developed using a training set of 378 patients across our health system and tested on a separate set of 93 patients. Ground truths for both sets were selected by an experienced radiologist. Performance metrics including yield (how many labeled images the model produced) and accuracy (percentage correct) were measured using the test set. MADLaP achieved a yield of 63% and an accuracy of 83%. The yield progressively increased as the input data moved through each module, while accuracy peaked part way through. Error analysis showed that inputs from certain examination sites had lower accuracy (40%) than the other sites (90%, 100%). MADLaP successfully created curated datasets of labeled ultrasound images of thyroid nodules. While accurate, the relatively suboptimal yield of MADLaP exposed some challenges when trying to automatically label radiology images from heterogeneous sources. The complex task of image curation and annotation could be automated, allowing for enrichment of larger datasets for use in machine learning development. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Variance of the Hellings-Downs Correlation

Authors: Bruce Allen

Abstract: Gravitational waves (GWs) create correlations in the arrival times of pulses from different pulsars. The expected correlation $μ(γ)$ as a function of the angle $γ$ between the directions to two pulsars was calculated by Hellings and Downs for an isotropic and unpolarized GW background, and several pulsar timing array (PTA) collaborations are working to observe these. We ask: given a set of noise-f… ▽ More Gravitational waves (GWs) create correlations in the arrival times of pulses from different pulsars. The expected correlation $μ(γ)$ as a function of the angle $γ$ between the directions to two pulsars was calculated by Hellings and Downs for an isotropic and unpolarized GW background, and several pulsar timing array (PTA) collaborations are working to observe these. We ask: given a set of noise-free observations, are they consistent with that expectation? To answer this, we calculate the expected variance $σ^2(γ)$ in the correlation for a single GW point source, as pulsar pairs with fixed separation angle $γ$ are swept around the sky. We then use this to derive simple analytic expressions for the variance produced by a set of discrete point sources uniformly scattered in space for two cases of interest: (1) point sources radiating GWs at the same frequency, generating confusion noise, and (2) point sources radiating GWs at distinct non-overlapping frequencies. By averaging over all pulsar sky positions at fixed separation angle $γ$, we show how this variance may be cleanly split into cosmic variance and pulsar variance, also demonstrating that measurements of the variance can provide information about the nature of GW sources. In a series of technical appendices, we calculate the mean and variance of the Hellings-Downs correlation for an arbitrary (polarized) point source, quantify the impact of neglecting pulsar terms, and calculate the pulsar and cosmic variance for a Gaussian ensemble. The mean and variance of the Gaussian ensemble may also be obtained from the previous discrete-source confusion-noise model in the limit of a high density of weak sources. △ Less

Submitted 24 February, 2023; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: Final published version

Journal ref: PRD 107, 043018 (2023)

Federated Learning Enables Big Data for Rare Cancer Boundary Detection

Authors: Sarthak Pati, Ujjwal Baid, Brandon Edwards, Micah Sheller, Shih-Han Wang, G Anthony Reina, Patrick Foley, Alexey Gruzdev, Deepthi Karkada, Christos Davatzikos, Chiharu Sako, Satyam Ghodasara, Michel Bilello, Suyash Mohan, Philipp Vollmuth, Gianluca Brugnara, Chandrakanth J Preetha, Felix Sahm, Klaus Maier-Hein, Maximilian Zenk, Martin Bendszus, Wolfgang Wick, Evan Calabrese, Jeffrey Rudie, Javier Villanueva-Meyer , et al. (254 additional authors not shown)

Abstract: Although machine learning (ML) has shown promise in numerous domains, there are concerns about generalizability to out-of-sample data. This is currently addressed by centrally sharing ample, and importantly diverse, data from multiple sites. However, such centralization is challenging to scale (or even not feasible) due to various limitations. Federated ML (FL) provides an alternative to train acc… ▽ More Although machine learning (ML) has shown promise in numerous domains, there are concerns about generalizability to out-of-sample data. This is currently addressed by centrally sharing ample, and importantly diverse, data from multiple sites. However, such centralization is challenging to scale (or even not feasible) due to various limitations. Federated ML (FL) provides an alternative to train accurate and generalizable ML models, by only sharing numerical model updates. Here we present findings from the largest FL study to-date, involving data from 71 healthcare institutions across 6 continents, to generate an automatic tumor boundary detector for the rare disease of glioblastoma, utilizing the largest dataset of such patients ever used in the literature (25,256 MRI scans from 6,314 patients). We demonstrate a 33% improvement over a publicly trained model to delineate the surgically targetable tumor, and 23% improvement over the tumor's entire extent. We anticipate our study to: 1) enable more studies in healthcare informed by large and diverse data, ensuring meaningful results for rare diseases and underrepresented populations, 2) facilitate further quantitative analyses for glioblastoma via performance optimization of our consensus model for eventual public release, and 3) demonstrate the effectiveness of FL at such scale and task complexity as a paradigm shift for multi-site collaborations, alleviating the need for data sharing. △ Less

Submitted 25 April, 2022; v1 submitted 22 April, 2022; originally announced April 2022.

Comments: federated learning, deep learning, convolutional neural network, segmentation, brain tumor, glioma, glioblastoma, FeTS, BraTS

Testing Station for Fast Screening of Through Silicon Via (TSV)-enabled Application Specific Integrated Circuits (ASICs) for Hard X-ray Imaging Detectors

Authors: Daniel P. Violette, Branden Allen, Jaesub Hong, Hiromasa Miyasaka, Jonathan Grindlay

Abstract: Application Specific Integrated Circuits (ASICs) are used in space-borne instruments for signal processing and detector readout. The electrical interface of these ASICs to frontend printed circuit boards (PCBs) is commonly accomplished with wire bonds. Through Silicon Via (TSV) technology has been proposed as an alternative interconnect technique that will reduce assembly complexity of ASIC packag… ▽ More Application Specific Integrated Circuits (ASICs) are used in space-borne instruments for signal processing and detector readout. The electrical interface of these ASICs to frontend printed circuit boards (PCBs) is commonly accomplished with wire bonds. Through Silicon Via (TSV) technology has been proposed as an alternative interconnect technique that will reduce assembly complexity of ASIC packaging by replacing wire bonding with flip-chip bonding. TSV technology is advantageous in large detector arrays where TSVs enable close detector tiling on all sides. Wafer-level probe card testing of TSV ASICs is frustrated by solder balls introduced onto the ASIC surface for flip-chip bonding that hamper alignment. Therefore, we developed the ASIC Test Stand (ATS) to enable rapid screening and characterization of individual ASIC die. We successfully demonstrated ATS operation on ASICs originally developed for CdZnTe detectors on the Nuclear Spectroscopic and Telescope Array (NuSTAR) mission that were later modified with TSVs in a via-last process. We tested both back-side blind-TSVs and front-side through-TSVs, with results from internal test pulser measurements that demonstrate performance equal to or exceeding the probe card wafer-level testing data. The ATS can easily be expanded or duplicated in order to parallelize ASIC screening for large area imaging detectors of future space programs. △ Less

Submitted 6 April, 2022; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 12 pages, 9 figures, submitted to SPIE Journal of Astronomical Telescopes, Instruments, and Systems (February 2022)

Performance of Random Template Banks

Authors: Bruce Allen

Abstract: When searching for new gravitational-wave or electromagnetic sources, the $n$ signal parameters (masses, sky location, frequencies,...) are unknown. In practice, one hunts for signals at a discrete set of points in parameter space, called a template bank. These may be constructed systematically as a lattice, or alternatively, by placing templates at randomly selected points in parameter space. Her… ▽ More When searching for new gravitational-wave or electromagnetic sources, the $n$ signal parameters (masses, sky location, frequencies,...) are unknown. In practice, one hunts for signals at a discrete set of points in parameter space, called a template bank. These may be constructed systematically as a lattice, or alternatively, by placing templates at randomly selected points in parameter space. Here, we calculate the fraction of signals lost by an $n$-dimensional random template bank (compared to a very finely spaced bank). This fraction is compared to the corresponding loss fraction for the best possible lattice-based template banks containing the same number of grid points. For dimensions $n<4$ the lattice-based template banks significantly outperform the random ones. However, remarkably, for dimensions $n>8$, the difference is negligible. In high dimensions, random template banks outperform the best known lattices. △ Less

Submitted 31 May, 2022; v1 submitted 5 March, 2022; originally announced March 2022.

Comments: Final published version. 9 pages, 3 figures, 1 table

Journal ref: Phys. Rev. D 105, 102003 (2022)

Scaling K2. V. Statistical Validation of 60 New Exoplanets From K2 Campaigns 2-18

Authors: Jessie L. Christiansen, Sakhee Bhure, Jon K. Zink, Kevin K. Hardegree-Ullman, Britt Duffy Adkins, Christina Hedges, Timothy D. Morton, Allyson Bieryla, David R. Ciardi, William D. Cochran, Courtney D. Dressing, Mark E. Everett, Howard Isaacson, John H. Livingston, Carl Ziegler, Perry Berlind, Michael L. Calkins, Gilbert A. Esquerdo, David W. Latham, Michael Endl, Phillip J. MacQueen, Benjamin J. Fulton, Lea A. Hirsch, Andrew W. Howard, Lauren M. Weiss , et al. (17 additional authors not shown)

Abstract: The NASA K2 mission, salvaged from the hardware failures of the Kepler telescope, has continued Kepler's planet-hunting success. It has revealed nearly 500 transiting planets around the ecliptic plane, many of which are the subject of further study, and over 1000 additional candidates. Here we present the results of an ongoing project to follow-up and statistically validate new K2 planets, in part… ▽ More The NASA K2 mission, salvaged from the hardware failures of the Kepler telescope, has continued Kepler's planet-hunting success. It has revealed nearly 500 transiting planets around the ecliptic plane, many of which are the subject of further study, and over 1000 additional candidates. Here we present the results of an ongoing project to follow-up and statistically validate new K2 planets, in particular to identify promising new targets for further characterization. By analyzing the reconnaissance spectra, high-resolution imaging, centroid variations, and statistical likelihood of the signals of 91 candidates, we validate 60 new planets in 46 systems. These include: a number of planets amenable to transmission spectroscopy (K2-384 f, K2-387 b, K2-390 b, K2-403 b, and K2-398 c), emission spectroscopy (K2-371 b, K2-370 b, and K2-399 b), and both (K2-405 b and K2-406 b); several systems with planets in or close to mean motion resonances (K2-381, K2-398) including a compact, TRAPPIST-1-like system of five small planets orbiting a mid-M dwarf (K2-384); an ultra-short period sub-Saturn in the hot Saturn desert (K2-399 b); and a super-Earth orbiting a moderately bright (V=11.93), metal-poor ([Fe/H]=-0.579+/-0.080) host star (K2-408 b). In total we validate planets around 4 F stars, 26 G stars, 13 K stars, and 3 M dwarfs. In addition, we provide a list of 37 vetted planet candidates that should be prioritized for future follow-up observation in order to be confirmed or validated. △ Less

Submitted 8 March, 2022; v1 submitted 3 March, 2022; originally announced March 2022.

Comments: 50 pages, 14 figures, 7 tables, accepted in AJ; replaced to fix error in planet name (thanks Travis!) and error in figure label