Showing 1–50 of 260 results for author: Baker, T

The Anomalous Acceleration of PSR J2043+1711: Long-Period Orbital Companion or Stellar Flyby?

Authors: Thomas Donlon II, Sukanya Chakrabarti, Michael T. Lam, Daniel Huber, Daniel Hey, Enrico Ramirez-Ruiz, Benjamin Shappee, David L. Kaplan, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Paul R. Brook, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels, Peter A. Gentile , et al. (31 additional authors not shown)

Abstract: Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4σ$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr… ▽ More Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4σ$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in \textit{Gaia} DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6\%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 80 kyr, a companion mass of 0.3 $M_\odot$ (indicative of a white dwarf or low-mass star), and a semi-major axis of 2000 AU. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios. △ Less

Submitted 8 July, 2024; originally announced July 2024.

K-mouflage at high k: extending the reach of $\texttt{Hi-COLA}$

Authors: Ashim Sen Gupta, Bartolomeo Fiorini, Tessa Baker

Abstract: The $\texttt{Hi-COLA}$ code is an efficient dark matter simulation suite that flexibly handles the Horndeski family of modified gravity models. In this work we extend the scope of $\texttt{Hi-COLA}$ to accommodate Horndeski theories with K-mouflage screening, allowing for the computation of matter power spectra in the non-linear regime in these models. We explore the boost of the dark matter power… ▽ More The $\texttt{Hi-COLA}$ code is an efficient dark matter simulation suite that flexibly handles the Horndeski family of modified gravity models. In this work we extend the scope of $\texttt{Hi-COLA}$ to accommodate Horndeski theories with K-mouflage screening, allowing for the computation of matter power spectra in the non-linear regime in these models. We explore the boost of the dark matter power spectrum relative to GR-$Λ$CDM in K-mouflage gravity, and also discuss how large-scale structure computations change between the Einstein and Jordan frames. A dissection of the relative contributions of the modified background, linear growth, fifth force, and the conformal factor (a new inclusion to $\texttt{Hi-COLA}$) to the boost factor is presented. The ability of $\texttt{Hi-COLA}$ to run with general Horndeski models and multiple screening mechanisms makes it an ideal tool for testing gravity with upcoming galaxy survey data. △ Less

Submitted 30 June, 2024; originally announced July 2024.

Comments: 31 pages, 10 figures, 1 table. Comments welcome. The Hi-COLA code can be found at https://github.com/Hi-COLACode/Hi-COLA

Matter Power Spectra in Modified Gravity: A Comparative Study of Approximations and $N$-Body Simulations

Authors: Benjamin Bose, Ashim Sen Gupta, Bartolomeo Fiorini, Guilherme Brando, Farbod Hassani, Tessa Baker, Lucas Lombriser, Baojiu Li, Cheng-Zong Ruan, Cesar Hernandez-Aguayo, Luis Atayde, Noemi Frusciante

Abstract: Testing gravity and the concordance model of cosmology, $Λ$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and… ▽ More Testing gravity and the concordance model of cosmology, $Λ$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and K-mouflage. In particular, we compare the predictions from full $N$-body simulations, two $N$-body codes with approximate time integration schemes, a parametrised modified $N$-body implementation and the analytic halo model reaction approach. We find the modification to the $Λ$CDM spectrum is in $2\%$ agreement for $z\leq1$ and $k\leq 1~h/{\rm Mpc}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 20 pages, 4 figures, 4 tables

Spectral Introspection Identifies Group Training Dynamics in Deep Neural Networks for Neuroimaging

Authors: Bradley T. Baker, Vince D. Calhoun, Sergey M. Plis

Abstract: Neural networks, whice have had a profound effect on how researchers study complex phenomena, do so through a complex, nonlinear mathematical structure which can be difficult for human researchers to interpret. This obstacle can be especially salient when researchers want to better understand the emergence of particular model behaviors such as bias, overfitting, overparametrization, and more. In N… ▽ More Neural networks, whice have had a profound effect on how researchers study complex phenomena, do so through a complex, nonlinear mathematical structure which can be difficult for human researchers to interpret. This obstacle can be especially salient when researchers want to better understand the emergence of particular model behaviors such as bias, overfitting, overparametrization, and more. In Neuroimaging, the understanding of how such phenomena emerge is fundamental to preventing and informing users of the potential risks involved in practice. In this work, we present a novel introspection framework for Deep Learning on Neuroimaging data, which exploits the natural structure of gradient computations via the singular value decomposition of gradient components during reverse-mode auto-differentiation. Unlike post-hoc introspection techniques, which require fully-trained models for evaluation, our method allows for the study of training dynamics on the fly, and even more interestingly, allow for the decomposition of gradients based on which samples belong to particular groups of interest. We demonstrate how the gradient spectra for several common deep learning models differ between schizophrenia and control participants from the COBRE study, and illustrate how these trajectories may reveal specific training dynamics helpful for further analysis. △ Less

Submitted 17 June, 2024; originally announced June 2024.

The NANOGrav 15 yr Data Set: Chromatic Gaussian Process Noise Models for Six Pulsars

Authors: Bjorn Larsen, Chiara M. F. Mingarelli, Jeffrey S. Hazboun, Aurelien Chalumeau, Deborah C. Good, Joseph Simon, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Paul R. Brook, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels, Peter A. Gentile, Joseph Glaser, Ross J. Jennings , et al. (39 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characteriz… ▽ More Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect achromatic noise properties in each pulsar. The models we compare include existing DM models used by NANOGrav and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar dataset. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from $\log_{10}A_{\text{RN}} = -14.1^{+0.1}_{-0.1}$ to $-14.7^{+0.3}_{-0.5}$, and the spectral index broadens from $γ_{\text{RN}} = 2.6^{+0.5}_{-0.4}$ to $γ_{\text{RN}} = 3.5^{+1.2}_{-0.9}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs. △ Less

Submitted 23 May, 2024; originally announced May 2024.

The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background

Authors: Gabriella Agazie, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Lucas Brown, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels , et al. (58 additional authors not shown)

Abstract: The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analyt… ▽ More The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 σ- 1.9 σ$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 σ- 2.1 σ$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 10 pages, 8 figures, 1 appendix, submitted to ApJ

Shareability of steering in 2-producible states

Authors: Qiu-Cheng Song, Travis J. Baker, Howard M. Wiseman

Abstract: Quantum steering is the phenomenon whereby one party (Alice) proves entanglement by "steering'' the system of another party (Bob) into distinct ensembles of states, by performing different measurements on her subsystem. Here, we investigate steering in a network scenario involving $n$ parties, who each perform local measurements on part of a global quantum state, that is produced using only two-pa… ▽ More Quantum steering is the phenomenon whereby one party (Alice) proves entanglement by "steering'' the system of another party (Bob) into distinct ensembles of states, by performing different measurements on her subsystem. Here, we investigate steering in a network scenario involving $n$ parties, who each perform local measurements on part of a global quantum state, that is produced using only two-party entangled states, and mixing with ancillary separable states. We introduce three scenarios which can be straightforwardly implemented in standard quantum optics architecture, which we call random $\frac{n}{2}$-pair entanglement, random pair entanglement and semi-random pair entanglement. We study steerability of the states across two-party marginals which arise in the three scenarios, and derive analytically the necessary and sufficient steering criteria for different sets of measurement settings. Strikingly, using the semi-random pair entanglement construction, one party can steer every one of the $n-1$ other parties, for arbitrarily large $n$, using only two measurements. Finally, exploiting symmetry, we study various small network configurations (three or four parties) in the three scenarios, under different measurements and produced by different two-party entangled states. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 11 pages, 9 figures, Close to published version

Journal ref: Physical Review A 108, 012216 (2023)

Multiscale Neuroimaging Features for the Identification of Medication Class and Non-Responders in Mood Disorder Treatment

Authors: Bradley T. Baker, Mustafa S. Salman, Zening Fu, Armin Iraji, Elizabeth Osuch, Jeremy Bockholt, Vince D. Calhoun

Abstract: In the clinical treatment of mood disorders, the complex behavioral symptoms presented by patients and variability of patient response to particular medication classes can create difficulties in providing fast and reliable treatment when standard diagnostic and prescription methods are used. Increasingly, the incorporation of physiological information such as neuroimaging scans and derivatives int… ▽ More In the clinical treatment of mood disorders, the complex behavioral symptoms presented by patients and variability of patient response to particular medication classes can create difficulties in providing fast and reliable treatment when standard diagnostic and prescription methods are used. Increasingly, the incorporation of physiological information such as neuroimaging scans and derivatives into the clinical process promises to alleviate some of the uncertainty surrounding this process. Particularly, if neural features can help to identify patients who may not respond to standard courses of anti-depressants or mood stabilizers, clinicians may elect to avoid lengthy and side-effect-laden treatments and seek out a different, more effective course that might otherwise not have been under consideration. Previously, approaches for the derivation of relevant neuroimaging features work at only one scale in the data, potentially limiting the depth of information available for clinical decision support. In this work, we show that the utilization of multi spatial scale neuroimaging features - particularly resting state functional networks and functional network connectivity measures - provide a rich and robust basis for the identification of relevant medication class and non-responders in the treatment of mood disorders. We demonstrate that the generated features, along with a novel approach for fast and automated feature selection, can support high accuracy rates in the identification of medication class and non-responders as well as the identification of novel, multi-scale biomarkers. △ Less

Submitted 12 February, 2024; originally announced February 2024.

LISA Definition Study Report

Authors: Monica Colpi, Karsten Danzmann, Martin Hewitson, Kelly Holley-Bockelmann, Philippe Jetzer, Gijs Nelemans, Antoine Petiteau, David Shoemaker, Carlos Sopuerta, Robin Stebbins, Nial Tanvir, Henry Ward, William Joseph Weber, Ira Thorpe, Anna Daurskikh, Atul Deep, Ignacio Fernández Núñez, César García Marirrodriga, Martin Gehler, Jean-Philippe Halain, Oliver Jennrich, Uwe Lammers, Jonan Larrañaga, Maike Lieser, Nora Lützgendorf , et al. (86 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the e… ▽ More The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 155 pages, with executive summary and table of contents

Low-Rank Learning by Design: the Role of Network Architecture and Activation Linearity in Gradient Rank Collapse

Authors: Bradley T. Baker, Barak A. Pearlmutter, Robyn Miller, Vince D. Calhoun, Sergey M. Plis

Abstract: Our understanding of learning dynamics of deep neural networks (DNNs) remains incomplete. Recent research has begun to uncover the mathematical principles underlying these networks, including the phenomenon of "Neural Collapse", where linear classifiers within DNNs converge to specific geometrical structures during late-stage training. However, the role of geometric constraints in learning extends… ▽ More Our understanding of learning dynamics of deep neural networks (DNNs) remains incomplete. Recent research has begun to uncover the mathematical principles underlying these networks, including the phenomenon of "Neural Collapse", where linear classifiers within DNNs converge to specific geometrical structures during late-stage training. However, the role of geometric constraints in learning extends beyond this terminal phase. For instance, gradients in fully-connected layers naturally develop a low-rank structure due to the accumulation of rank-one outer products over a training batch. Despite the attention given to methods that exploit this structure for memory saving or regularization, the emergence of low-rank learning as an inherent aspect of certain DNN architectures has been under-explored. In this paper, we conduct a comprehensive study of gradient rank in DNNs, examining how architectural choices and structure of the data effect gradient rank bounds. Our theoretical analysis provides these bounds for training fully-connected, recurrent, and convolutional neural networks. We also demonstrate, both theoretically and empirically, how design choices like activation function linearity, bottleneck layer introduction, convolutional stride, and sequence truncation influence these bounds. Our findings not only contribute to the understanding of learning dynamics in DNNs, but also provide practical guidance for deep learning engineers to make informed design decisions. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Efficient prescription to search for linear gravitational wave memory from hyperbolic black hole encounters and its application to the NANOGrav 12.5-year dataset

Authors: Subhajit Dandapat, Abhimanyu Susobhanan, Lankeswar Dey, A. Gopakumar, Paul T. Baker, Philippe Jetzer

Abstract: Burst with memory events are potential transient gravitational wave sources for the maturing pulsar timing array (PTA) efforts. We provide a computationally efficient prescription to model pulsar timing residuals induced by supermassive black hole pairs in general relativistic hyperbolic trajectories employing a Keplerian-type parametric solution. Injection studies have been pursued on the resulti… ▽ More Burst with memory events are potential transient gravitational wave sources for the maturing pulsar timing array (PTA) efforts. We provide a computationally efficient prescription to model pulsar timing residuals induced by supermassive black hole pairs in general relativistic hyperbolic trajectories employing a Keplerian-type parametric solution. Injection studies have been pursued on the resulting bursts with linear GW memory (LGWM) events with simulated datasets to test the performance of our pipeline, followed by its application to the publicly available NANOGrav 12.5-year (NG12.5) dataset. Given the absence of any evidence of LGWM events within the real NG12.5 dataset, we impose $95\%$ upper limits on the PTA signal amplitude as a function of the sky location of the source and certain characteristic frequency ($n$) of the signal. The upper limits are computed using a signal model that takes into account the presence of intrinsic timing noise specific to each pulsar, as well as a common, spatially uncorrelated red noise, alongside the LGWM signal. Our investigations reveal that the $95\%$ upper limits on LGWM amplitude, marginalized over all other parameters, is 3.48 $\pm 0.51 \ μ$s for $n>3.16$ nHz. This effort should be relevant for constraining both burst and memory events in the upcoming International Pulsar Timing Array data releases. △ Less

Submitted 16 May, 2024; v1 submitted 5 February, 2024; originally announced February 2024.

Comments: 20 pages, 11 figures

Journal ref: Phys. Rev. D 109, 103018, 2024

Energy Efficiency Maximization in IRS-enabled Phase Cooperative PS-SWIPT based Self-sustainable IoT Network

Authors: Haleema Sadia, Ahmad Kamal Hassan, Ziaul Haq Abbas, Ghulam Abbas, Thar Baker

Abstract: Power splitting based simultaneous wireless information and power transfer (PS-SWIPT) appears to be a promising solution to support future self-sustainable Internet of Things (SS-IoT) networks. However, the performance of these networks is constrained by radio frequency signal strength and channel impairments. To address this challenge, intelligent reflecting surfaces (IRSs) are introduced in PS-S… ▽ More Power splitting based simultaneous wireless information and power transfer (PS-SWIPT) appears to be a promising solution to support future self-sustainable Internet of Things (SS-IoT) networks. However, the performance of these networks is constrained by radio frequency signal strength and channel impairments. To address this challenge, intelligent reflecting surfaces (IRSs) are introduced in PS-SWIPT based SS-IoT networks to improve network efficiency by controlling signal reflections. In this article, an IRS-enabled phase cooperative framework is proposed to improve energy efficiency (EE) of the IoT network $({\mathtt {I}}^{net})$ using phase shifts of the user network $({\mathtt {U}^{net})}$, without constraining hardware resources at ${\mathtt {U}^{net}}$. We exploit transmit beamforming (BF) at access points (APs) and phase shifts optimization at the IRS end with phase effective cooperation between APs to enhance ${\mathtt {I}}^{net}$ EE performance. The maximization problem turns out to be NP-hard, so first, an alternating optimization (AO) is solved for the ${\mathtt {U}^{net}}$ using low computational complexity heuristic BF approaches, namely, transmit minimum-mean-square-error and zero-forcing BF, and phase optimization is performed using semidefinite relaxation (SDR) approach. To combat the computational complexity of AO, we also propose an alternative solution by exploiting heuristic BF schemes and an iterative algorithm, i.e., the element-wise block-coordinate descent method for phase shifts optimization. Next, EE maximization is solved for the ${\mathtt {I}^{net}}$ by optimizing the PS ratio and active BF vectors by exploiting optimal phase shifts of the ${\mathtt {U}}^{net}$. Simulation results confirm that employing IRS phase cooperation in PS-SWIPT based SS-IoT networks can significantly improve EE performance of ${\mathtt {I}^{net}}$ without constraining resources. △ Less

Submitted 10 January, 2024; originally announced January 2024.

The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest , et al. (74 additional authors not shown)

Abstract: Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correl… ▽ More Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 11 pages, 5 figures

The clustering of dark sirens' invisible host galaxies

Authors: Charles Dalang, Tessa Baker

Abstract: Dark sirens are a powerful way to infer cosmological and astrophysical parameters from the combination of gravitational wave sirens and galaxy catalogues. Importantly, the method relies on the completeness of the galaxy catalogues being well modelled. A magnitude-limited catalogue will always be incomplete to some extent, requiring a completion scheme to avoid biasing the parameter inference. Stan… ▽ More Dark sirens are a powerful way to infer cosmological and astrophysical parameters from the combination of gravitational wave sirens and galaxy catalogues. Importantly, the method relies on the completeness of the galaxy catalogues being well modelled. A magnitude-limited catalogue will always be incomplete to some extent, requiring a completion scheme to avoid biasing the parameter inference. Standard methods include homogeneous and multiplicative completion, which have the advantage of simplicity but underestimate or overestimate the amplitude of structure at low completeness, respectively. In this work, we propose a new method to complete galaxy catalogues which uses clustering information to incorporate knowledge of the large scale structure into the dark sirens method. We find that if the structure of the true number of galaxies is sufficiently well preserved in the catalogue, our estimator can perform drastically better than both homogeneous and multiplicative completion. We lay the foundations for a maximally informative dark sirens analysis and discuss its limitations. △ Less

Submitted 9 July, 2024; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 29 pages, 11 figures, adapted to match published version in JCAP

Fast production of cosmological emulators in modified gravity: the matter power spectrum

Authors: Bartolomeo Fiorini, Kazuya Koyama, Tessa Baker

Abstract: We test the convergence of fast simulations based on the COmoving Lagrangian Acceleration (COLA) method for predictions of the matter power spectrum, specialising our analysis in the redshift range $1 \le z \le 1.65$, relevant to high-redshift spectroscopic galaxy surveys. We then focus on the enhancement of the matter power spectrum in modified gravity (MG), the boost factor, using the Dvali-Gaba… ▽ More We test the convergence of fast simulations based on the COmoving Lagrangian Acceleration (COLA) method for predictions of the matter power spectrum, specialising our analysis in the redshift range $1 \le z \le 1.65$, relevant to high-redshift spectroscopic galaxy surveys. We then focus on the enhancement of the matter power spectrum in modified gravity (MG), the boost factor, using the Dvali-Gabadadze-Porrati (DGP) theory as a test case but developing a general approach that can be applied to other MG theories. After identifying the minimal simulation requirements for accurate DGP boost factors, we design and produce a COLA simulation suite that we use to train a neural network emulator for the DGP boost factor. Using MG-AREPO simulations as a reference, we estimate the emulator accuracy to be of $\sim 3\%$ up to $k=5 \, h {\rm Mpc}^{-1}$ at $0 \leq z \leq 2$. We make the emulator publicly available at: https://github.com/BartolomeoF/nDGPemu. △ Less

Submitted 20 November, 2023; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 24 pages, 17 figures, Comments are welcome - Accepted for publication in JCAP

The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Belinda D. Cheeseboro, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Lankeswar Dey, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara , et al. (63 additional authors not shown)

Abstract: The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for contin… ▽ More The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity $e_0<0.5$ and the symmetric mass ratio $η>0.1$. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of $88.1\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $η>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations. △ Less

Submitted 15 January, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

How to Detect an Astrophysical Nanohertz Gravitational-Wave Background

Authors: Bence Bécsy, Neil J. Cornish, Patrick M. Meyers, Luke Zoltan Kelley, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Tyler Cohen, James M. Cordes, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch , et al. (71 additional authors not shown)

Abstract: Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, nam… ▽ More Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them. △ Less

Submitted 1 December, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 14 pages, 8 figures, version matching published paper

Journal ref: ApJ 959 9 (2023)

Magnification and evolution bias of transient sources: GWs and SNIa

Authors: Stefano Zazzera, José Fonseca, Tessa Baker, Chris Clarkson

Abstract: Third-generation gravitational wave (GW) observatories such as the Einstein Telescope and Cosmic Explorer, together with the LSST survey at the Vera Rubin Observatory, will yield an abundance of extra-galactic transient objects. This opens the exciting possibility of using GW sources and Supernovae Type Ia (SNIa) as luminosity distance tracers of large-scale structure for the first time. The large… ▽ More Third-generation gravitational wave (GW) observatories such as the Einstein Telescope and Cosmic Explorer, together with the LSST survey at the Vera Rubin Observatory, will yield an abundance of extra-galactic transient objects. This opens the exciting possibility of using GW sources and Supernovae Type Ia (SNIa) as luminosity distance tracers of large-scale structure for the first time. The large volumes accessible to these surveys imply that we may need to include relativistic corrections, such as lensing and Doppler magnification. However, the amplitude of these effects depends on the magnification and evolution biases of the transient sources, which are not yet understood. In this paper we develop comprehensive frameworks to address and model these biases for both populations of transient objects; in particular, we define how to compute these biases for GW sources. We then analyse the impact of magnification and evolution biases on the relativistic corrections and on the angular power spectrum of these sources. We show that correct modelling and implementation of these biases is crucial for measuring the cross-correlations of transient sources at higher redshifts. △ Less

Submitted 5 July, 2024; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 28 pages, 14 figures, 4 tables, comments welcome

Journal ref: JCAP05(2024)095

Testing the nature of gravitational wave propagation using dark sirens and galaxy catalogues

Authors: Anson Chen, Rachel Gray, Tessa Baker

Abstract: The dark sirens method enables us to use gravitational wave events without electromagnetic counterparts as tools for cosmology and tests of gravity. Furthermore, the dark sirens analysis code gwcosmo can now robustly account for information coming from both galaxy catalogues and the compact object mass distribution. We present here an extension of the gwcosmo code and methodology to constrain para… ▽ More The dark sirens method enables us to use gravitational wave events without electromagnetic counterparts as tools for cosmology and tests of gravity. Furthermore, the dark sirens analysis code gwcosmo can now robustly account for information coming from both galaxy catalogues and the compact object mass distribution. We present here an extension of the gwcosmo code and methodology to constrain parameterized deviations from General Relativity that affect the propagation of gravitational waves. We show results of our analysis using data from the GWTC-3 gravitational wave catalogues, in preparation for application to the O4 observing run. After testing our pipelines using the First Two Years mock data set, we reanalyse 46 events from GWTC-3, and combine the posterior for BBH and NSBH sampling results for the first time. We obtain joint constraints on H0 and parameterized deviations from General Relativity in the Power Law + Peak BBH population model. With increased galaxy catalogue support in the future, our work sets the stage for dark sirens to become a powerful tool for testing gravity. △ Less

Submitted 25 September, 2023; v1 submitted 7 September, 2023; originally announced September 2023.

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

Scalable multiparty steering based on a single pair of entangled qubits

Authors: Alex Pepper, Travis. J. Baker, Yuanlong Wang, Qiu-Cheng Song, Lynden. K. Shalm, Varun. B. Varma, Sae Woo Nam, Nora Tischler, Sergei Slussarenko, Howard. M. Wiseman, Geoff. J. Pryde

Abstract: The distribution and verification of quantum nonlocality across a network of users is essential for future quantum information science and technology applications. However, beyond simple point-to-point protocols, existing methods struggle with increasingly complex state preparation for a growing number of parties. Here, we show that, surprisingly, multiparty loophole-free quantum steering, where o… ▽ More The distribution and verification of quantum nonlocality across a network of users is essential for future quantum information science and technology applications. However, beyond simple point-to-point protocols, existing methods struggle with increasingly complex state preparation for a growing number of parties. Here, we show that, surprisingly, multiparty loophole-free quantum steering, where one party simultaneously steers arbitrarily many spatially separate parties, is achievable by constructing a quantum network from a set of qubits of which only one pair is entangled. Using these insights, we experimentally demonstrate this type of steering between three parties with the detection loophole closed. With its modest and fixed entanglement requirements, this work introduces a scalable approach to rigorously verify quantum nonlocality across multiple parties, thus providing a practical tool towards developing the future quantum internet. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Joint cosmological and gravitational-wave population inference using dark sirens and galaxy catalogues

Authors: Rachel Gray, Freija Beirnaert, Christos Karathanasis, Benoît Revenu, Cezary Turski, Anson Chen, Tessa Baker, Sergio Vallejo, Antonio Enea Romano, Tathagata Ghosh, Archisman Ghosh, Konstantin Leyde, Simone Mastrogiovanni, Surhud More

Abstract: In the absence of numerous gravitational-wave detections with confirmed electromagnetic counterparts, the "dark siren" method has emerged as a leading technique of gravitational-wave cosmology. The method allows redshift information of such events to be inferred statistically from a catalogue of potential host galaxies. Due to selection effects, dark siren analyses necessarily depend on the mass d… ▽ More In the absence of numerous gravitational-wave detections with confirmed electromagnetic counterparts, the "dark siren" method has emerged as a leading technique of gravitational-wave cosmology. The method allows redshift information of such events to be inferred statistically from a catalogue of potential host galaxies. Due to selection effects, dark siren analyses necessarily depend on the mass distribution of compact objects and the evolution of their merger rate with redshift. Informative priors on these quantities will impact the inferred posterior constraints on the Hubble constant ($H_0$). It is thus crucial to vary these unknown distributions during an $H_0$ inference. This was not possible in earlier analyses due to the high computational cost, restricting them to either excluding galaxy catalogue information, or fixing the gravitational-wave population mass distribution and risking introducing bias to the $H_0$ measurement. This paper introduces a significantly enhanced version of the Python package GWCOSMO, which allows joint estimation of cosmological and compact binary population parameters. This thereby ensures the analysis is now robust to a major source of potential bias. The gravitational-wave events from the Third Gravitational-Wave Transient Catalogue are reanalysed with the GLADE+ galaxy catalogue, and an updated, more reliable measurement of $H_0=69^{+12}_{-7}$ km s$^{-1}$ Mpc$^{-1}$ is found (maximum a posteriori probability and 68% highest density interval). This improved method will enable cosmological analyses with future gravitational-wave detections to make full use of the information available (both from galaxy catalogues and the compact binary population itself), leading to promising new independent bounds on the Hubble constant. △ Less

Submitted 20 November, 2023; v1 submitted 4 August, 2023; originally announced August 2023.

Comments: 31 pages, 11 figures, accepted for publication in JCAP

Synthetic pre-training for neural-network interatomic potentials

Authors: John L. A. Gardner, Kathryn T. Baker, Volker L. Deringer

Abstract: Machine learning (ML) based interatomic potentials have transformed the field of atomistic materials modelling. However, ML potentials depend critically on the quality and quantity of quantum-mechanical reference data with which they are trained, and therefore developing datasets and training pipelines is becoming an increasingly central challenge. Leveraging the idea of "synthetic" (artificial) d… ▽ More Machine learning (ML) based interatomic potentials have transformed the field of atomistic materials modelling. However, ML potentials depend critically on the quality and quantity of quantum-mechanical reference data with which they are trained, and therefore developing datasets and training pipelines is becoming an increasingly central challenge. Leveraging the idea of "synthetic" (artificial) data that is common in other areas of ML research, we here show that synthetic atomistic data, themselves obtained at scale with an existing ML potential, constitute a useful pre-training task for neural-network interatomic potential models. Once pre-trained with a large synthetic dataset, these models can be fine-tuned on a much smaller, quantum-mechanical one, improving numerical accuracy and stability in computational practice. We demonstrate feasibility for a series of equivariant graph-neural-network potentials for carbon, and we carry out initial experiments to test the limits of the approach. △ Less

Submitted 24 July, 2023; originally announced July 2023.

The NANOGrav 12.5-year Data Set: Search for Gravitational Wave Memory

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis , et al. (65 additional authors not shown)

Abstract: We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a sign… ▽ More We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately $3.3 \times 10^{-14}$. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613$-$0200. This suggests that this pulsar has some excess noise which can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 29 pages, 5 figures

The NANOGrav 15-year Gravitational-Wave Background Analysis Pipeline

Authors: Aaron D. Johnson, Patrick M. Meyers, Paul T. Baker, Neil J. Cornish, Jeffrey S. Hazboun, Tyson B. Littenberg, Joseph D. Romano, Stephen R. Taylor, Michele Vallisneri, Sarah J. Vigeland, Ken D. Olum, Xavier Siemens, Justin A. Ellis, Rutger van Haasteren, Sophie Hourihane, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde , et al. (71 additional authors not shown)

Abstract: This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{P… ▽ More This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data. △ Less

Submitted 7 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 10 figures, 1 table; Companion paper to "The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background"; For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan DeCesar, Paul B. Demorest, Matthew C. Digman, Timothy Dolch, Brendan Drachler , et al. (74 additional authors not shown)

Abstract: Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav… ▽ More Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of $\sim$4 nHz, and another at $\sim$170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings-Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 23 pages, 13 figures, 2 tables. Accepted for publication in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore , et al. (68 additional authors not shown)

Abstract: The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these… ▽ More The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Paul T. Baker, Bence Bécsy, Laura Blecha, Alexander Bonilla, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler , et al. (89 additional authors not shown)

Abstract: The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary popul… ▽ More The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations. △ Less

Submitted 18 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

The NANOGrav 15-year Data Set: Search for Signals from New Physics

Authors: Adeela Afzal, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Jose Juan Blanco-Pillado, Laura Blecha, Kimberly K. Boddy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie , et al. (98 additional authors not shown)

Abstract: The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic string… ▽ More The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 74 pages, 31 figures, 4 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. Decesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca , et al. (66 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated… ▽ More Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white noise parameters and two red noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of $7\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars

Authors: Gabriella Agazie, Md Faisal Alam, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Victoria Bonidie, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, Christopher Chapman, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler , et al. (75 additional authors not shown)

Abstract: We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival… ▽ More We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by three years, now spanning nearly 16 years for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic gravitational-wave background. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 90 pages, 74 figures, 6 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Becsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler, Megan E. DeCesar , et al. (89 additional authors not shown)

Abstract: We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectr… ▽ More We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3σ$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4σ$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

Lightweight Attribute Localizing Models for Pedestrian Attribute Recognition

Authors: Ashish Jha, Dimitrii Ermilov, Konstantin Sobolev, Anh Huy Phan, Salman Ahmadi-Asl, Naveed Ahmed, Imran Junejo, Zaher AL Aghbari, Thar Baker, Ahmed Mohamed Khedr, Andrzej Cichocki

Abstract: Pedestrian Attribute Recognition (PAR) deals with the problem of identifying features in a pedestrian image. It has found interesting applications in person retrieval, suspect re-identification and soft biometrics. In the past few years, several Deep Neural Networks (DNNs) have been designed to solve the task; however, the developed DNNs predominantly suffer from over-parameterization and high com… ▽ More Pedestrian Attribute Recognition (PAR) deals with the problem of identifying features in a pedestrian image. It has found interesting applications in person retrieval, suspect re-identification and soft biometrics. In the past few years, several Deep Neural Networks (DNNs) have been designed to solve the task; however, the developed DNNs predominantly suffer from over-parameterization and high computational complexity. These problems hinder them from being exploited in resource-constrained embedded devices with limited memory and computational capacity. By reducing a network's layers using effective compression techniques, such as tensor decomposition, neural network compression is an effective method to tackle these problems. We propose novel Lightweight Attribute Localizing Models (LWALM) for Pedestrian Attribute Recognition (PAR). LWALM is a compressed neural network obtained after effective layer-wise compression of the Attribute Localization Model (ALM) using the Canonical Polyadic Decomposition with Error Preserving Correction (CPD-EPC) algorithm. △ Less

Submitted 16 June, 2023; originally announced June 2023.

The observed number counts in luminosity distance space

Authors: José Fonseca, Stefano Zazzera, Tessa Baker, Chris Clarkson

Abstract: Next generation surveys will provide us with an unprecedented number of detections of supernovae Type Ia and gravitational wave merger events. Cross-correlations of such objects offer novel and powerful insights into the large-scale distribution of matter in the universe. Both of these sources carry information on their luminosity distance, but remain uninformative about their redshifts; hence the… ▽ More Next generation surveys will provide us with an unprecedented number of detections of supernovae Type Ia and gravitational wave merger events. Cross-correlations of such objects offer novel and powerful insights into the large-scale distribution of matter in the universe. Both of these sources carry information on their luminosity distance, but remain uninformative about their redshifts; hence their clustering analyses and cross-correlations need to be carried out in luminosity distance space, as opposed to redshift space. In this paper, we calculate the full expression for the number count fluctuation in terms of a perturbation to the observed luminosity distance. We find the expression to differ significantly from the one commonly used in redshift space. Furthermore, we present a comparison of the number count angular power spectra between luminosity distance and redshift spaces. We see a wide divergence between the two at large scales, and we note that lensing is the main contribution to such differences. On such scales and at higher redshifts the difference between the angular power spectra in luminosity distance and redshift spaces can be roughly 50$\%$. We also investigate cross-correlating different redshift bins using different tracers, i.e. one in luminosity distance space and one in redshift, simulating the cross-correlation angular power spectrum between background gravitational waves/supernovae and foreground galaxies. Finally, we show that in a cosmic variance limited survey, the relativistic corrections to the density-only term ought to be included. △ Less

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

Comments: 34 pages, 9 figures, 1 table. Agrees with published version

Journal ref: JCAP08(2023)050

Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

Authors: M. Falxa, S. Babak, P. T. Baker, B. Bécsy, A. Chalumeau, S. Chen, Z. Chen, N. J. Cornish, L. Guillemot, J. S. Hazboun, C. M. F. Mingarelli, A. Parthasarathy, A. Petiteau, N. S. Pol, A. Sesana, S. B. Spolaor, S. R. Taylor, G. Theureau, M. Vallisneri, S. J. Vigeland, C. A. Witt, X. Zhu, J. Antoniadis, Z. Arzoumanian, M. Bailes , et al. (102 additional authors not shown)

Abstract: The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi… ▽ More The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit. △ Less

Submitted 19 March, 2023; originally announced March 2023.

Peculiar velocity effects on the Hubble constant from time-delay cosmography

Authors: Charles Dalang, Martin Millon, Tessa Baker

Abstract: Two major challenges of contemporary cosmology are the Hubble tension and the cosmic dipole tension. At the crossroad of these, we investigate the impact of peculiar velocities on estimations of the Hubble constant from time-delay cosmography. We quantify the bias on the inference of the Hubble constant due to peculiar velocities of the lens, the source and of the observer. The former two, which m… ▽ More Two major challenges of contemporary cosmology are the Hubble tension and the cosmic dipole tension. At the crossroad of these, we investigate the impact of peculiar velocities on estimations of the Hubble constant from time-delay cosmography. We quantify the bias on the inference of the Hubble constant due to peculiar velocities of the lens, the source and of the observer. The former two, which may cancel from one system to another, affect the determination of the angular diameter distances in the time-delay formula, and reconstructed quantities like the angle to the source, via a lens model. On the other hand, the peculiar velocity of the observer, which is a debated quantity in the context of the cosmic dipole tension, systematically affects observed angles through aberration, redshifts, angular diameter distance and reconstructed quantities. We compute in detail the effect of these peculiar velocities on the inference of the Hubble constant to linear order in the peculiar velocities for the seven lenses of the H0LiCOW/TDCOSMO collaboration. The bias generated by the observer's peculiar velocity alone can reach $1.15\%$ for the lenses which are well aligned with it. This results in a $0.25 \%$ bias for the seven combined lenses. Assuming a typical peculiar velocity of $300$ km s$^{-1}$ for the lens and the source galaxies, these add an additional random uncertainty, which can reach $1\%$ for an individual lens, but reduces to $0.24\%$ for the full TDCOSMO sample. The picture may change if peculiar velocities turn out to be larger than expected. Any time-delay cosmography program which aims for percent precision on the Hubble constant may need to take this source of systematic bias into account. This is especially so for future ground-based surveys which cover a fraction of the celestial sphere that is well aligned with the observer's peculiar velocity. △ Less

Submitted 2 October, 2023; v1 submitted 13 January, 2023; originally announced January 2023.

Comments: 19 pages, 8 figures, 2 tables, online version adapted to make it consistent with published version

The NANOGrav 12.5-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis, E. C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman , et al. (53 additional authors not shown)

Abstract: Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data s… ▽ More Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data set. We created new methods to accurately model the uncertainties on pulsar distances in our analysis, and we implemented new techniques to account for a common red noise process in pulsar timing array data sets while searching for deterministic gravitational wave signals, including continuous waves. As we found no evidence for continuous waves in our data, we placed 95\% upper limits on the strain amplitude of continuous waves emitted by these sources. At our most sensitive frequency of 7.65 nanohertz, we placed a sky-averaged limit of $h_0 < $ $(6.82 \pm 0.35) \times 10^{-15}$, and $h_0 <$ $(2.66 \pm 0.15) \times 10^{-15}$ in our most sensitive sky location. Finally, we placed a multi-messenger limit of $\mathcal{M} <$ $(1.41 \pm 0.02) \times 10^9 M_\odot$ on the chirp mass of the supermassive black hole binary candidate 3C~66B. △ Less

Submitted 6 June, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 22 pages, 12 figures. Accepted by ApJL

On the power of one pure steered state for EPR-steering with a pair of qubits

Authors: Qiu-Cheng Song, Travis J. Baker, Howard M. Wiseman

Abstract: As originally introduced, the EPR phenomenon was the ability of one party (Alice) to steer, by her choice between two measurement settings, the quantum system of another party (Bob) into two distinct ensembles of pure states. As later formalized as a quantum information task, EPR-steering can be shown even when the distinct ensembles comprise mixed states. Consider the scenario where Alice and Bob… ▽ More As originally introduced, the EPR phenomenon was the ability of one party (Alice) to steer, by her choice between two measurement settings, the quantum system of another party (Bob) into two distinct ensembles of pure states. As later formalized as a quantum information task, EPR-steering can be shown even when the distinct ensembles comprise mixed states. Consider the scenario where Alice and Bob each have a qubit and Alice performs dichotomic projective measurements. In this case, the states in the ensembles to which she can steer form the surface of an ellipsoid ${\cal E}$ in Bob's Bloch ball. Further, let the steering ellipsoid ${\cal E}$ have nonzero volume. It has previously been shown that if Alice's first measurement setting yields an ensemble comprising two pure states, then this, plus any one other measurement setting, will demonstrate EPR-steering. Here we consider what one can say if the ensemble from Alice's first setting contains only one pure state $\mathsf{p}\in{\cal E}$, occurring with probability $p_\mathsf{p}$. Using projective geometry, we derive the necessary and sufficient condition analytically for Alice to be able to demonstrate EPR-steering of Bob's state using this and some second setting, when the two ensembles from these lie in a given plane. Based on this, we show that, for a given ${\cal E}$, if $p_\mathsf{p}$ is high enough [$p_{\sf p} > p_{\rm max}^{\cal E} \in [0,1)$] then any distinct second setting by Alice is sufficient to demonstrate EPR-steering. Similarly we derive a $p_{\rm min}^{\cal E}$ such that $p_\mathsf{p}>p_{\rm min}^{\cal E}$ is necessary for Alice to demonstrate EPR-steering using only the first setting and some other setting. Moreover, the expressions we derive are tight; for spherical steering ellipsoids, the bounds coincide: $p_{\rm max}^{\cal E} = p_{\rm min}^{\cal E}$. △ Less

Submitted 3 April, 2024; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: 30 pages, 9 figures, v2: Close to published version

Journal ref: New J. Phys. 25 (2023) 053005

The Hitchhiker's guide to the galaxy catalog approach for gravitational wave cosmology

Authors: Jonathan R. Gair, Archisman Ghosh, Rachel Gray, Daniel E. Holz, Simone Mastrogiovanni, Suvodip Mukherjee, Antonella Palmese, Nicola Tamanini, Tessa Baker, Freija Beirnaert, Maciej Bilicki, Hsin-Yu Chen, Gergely Dálya, Jose Maria Ezquiaga, Will M. Farr, Maya Fishbach, Juan Garcia-Bellido, Tathagata Ghosh, Hsiang-Yu Huang, Christos Karathanasis, Konstantin Leyde, Ignacio Magaña Hernandez, Johannes Noller, Gregoire Pierra, Peter Raffai , et al. (6 additional authors not shown)

Abstract: We outline the ``dark siren'' galaxy catalog method for cosmological inference using gravitational wave (GW) standard sirens, clarifying some common misconceptions in the implementation of this method. When a confident transient electromagnetic counterpart to a GW event is unavailable, the identification of a unique host galaxy is in general challenging. Instead, as originally proposed by Schutz (… ▽ More We outline the ``dark siren'' galaxy catalog method for cosmological inference using gravitational wave (GW) standard sirens, clarifying some common misconceptions in the implementation of this method. When a confident transient electromagnetic counterpart to a GW event is unavailable, the identification of a unique host galaxy is in general challenging. Instead, as originally proposed by Schutz (1986), one can consult a galaxy catalog and implement a dark siren statistical approach incorporating all potential host galaxies within the localization volume. Trott & Hunterer 2021 recently claimed that this approach results in a biased estimate of the Hubble constant, $H_0$, when implemented on mock data, even if optimistic assumptions are made. We demonstrate explicitly that, as previously shown by multiple independent groups, the dark siren statistical method leads to an unbiased posterior when the method is applied to the data correctly. We highlight common sources of error possible to make in the generation of mock data and implementation of the statistical framework, including the mismodeling of selection effects and inconsistent implementations of the Bayesian framework, which can lead to a spurious bias. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 19 pages, 11 figures. Submitted to ApJ

An unusual pulse shape change event in PSR J1713+0747 observed with the Green Bank Telescope and CHIME

Authors: Ross J. Jennings, James M. Cordes, Shami Chatterjee, Maura A. McLaughlin, Paul B. Demorest, Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Paul R. Brook, Tyler Cohen, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Timothy Dolch, Elizabeth C. Ferrara, Emmanuel Fonseca, Deborah C. Good, Jeffrey S. Hazboun, Megan L. Jones, David L. Kaplan, Michael T. Lam, T. Joseph W. Lazio, Duncan R. Lorimer, Jing Luo, Ryan S. Lynch , et al. (19 additional authors not shown)

Abstract: The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter G… ▽ More The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter Green Bank Telescope (GBT) in a three-year period encompassing the shape change event, between February 2020 and February 2023. As of February 2023, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying TOA residuals display a strong non-monotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency, $ν$) nor a change in dispersion measure (DM) alone (which would produce a delay proportional to $ν^{-2}$). However, it does bear some resemblance to the two previous "chromatic timing events" observed in J1713+0747 (Demorest et al. 2013; Lam et al. 2016), as well as to a similar event observed in PSR J1643-1224 in 2015 (Shannon et al. 2016). △ Less

Submitted 31 January, 2024; v1 submitted 21 October, 2022; originally announced October 2022.

Comments: 18 pages, 8 figures. Submitted to ApJ. Data available at https://doi.org/10.5281/zenodo.7236459

Testing unit root non-stationarity in the presence of missing data in univariate time series of mobile health studies

Authors: Charlotte Fowler, Xiaoxuan Cai, Justin T. Baker, Jukka-Pekka Onnela, Linda Valeri

Abstract: The use of digital devices to collect data in mobile health (mHealth) studies introduces a novel application of time series methods, with the constraint of potential data missing at random (MAR) or missing not at random (MNAR). In time series analysis, testing for stationarity is an important preliminary step to inform appropriate later analyses. The augmented Dickey-Fuller (ADF) test was develope… ▽ More The use of digital devices to collect data in mobile health (mHealth) studies introduces a novel application of time series methods, with the constraint of potential data missing at random (MAR) or missing not at random (MNAR). In time series analysis, testing for stationarity is an important preliminary step to inform appropriate later analyses. The augmented Dickey-Fuller (ADF) test was developed to test the null hypothesis of unit root non-stationarity, under no missing data. Beyond recommendations under data missing completely at random (MCAR) for complete case analysis or last observation carry forward imputation, researchers have not extended unit root non-stationarity testing to a context with more complex missing data mechanisms. Multiple imputation with chained equations, Kalman smoothing imputation, and linear interpolation have also been proposed for time series data, however such methods impose constraints on the autocorrelation structure, and thus impact unit root testing. We propose maximum likelihood estimation and multiple imputation using state space model approaches to adapt the ADF test to a context with missing data. We further develop sensitivity analysis techniques to examine the impact of MNAR data. We evaluate the performance of existing and proposed methods across different missing mechanisms in extensive simulations and in their application to a multi-year smartphone study of bipolar patients. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Testing gravitational wave propagation with multiband detections

Authors: Tessa Baker, Enrico Barausse, Anson Chen, Claudia de Rham, Mauro Pieroni, Gianmassimo Tasinato

Abstract: Effective field theories (EFT) of dark energy (DE) -- built to parameterise the properties of DE in an agnostic manner -- are severely constrained by measurements of the propagation speed of gravitational waves (GW). However, GW frequencies probed by ground-based interferometers lie around the typical strong coupling scale of the EFT, and it is likely that the effective description breaks down bef… ▽ More Effective field theories (EFT) of dark energy (DE) -- built to parameterise the properties of DE in an agnostic manner -- are severely constrained by measurements of the propagation speed of gravitational waves (GW). However, GW frequencies probed by ground-based interferometers lie around the typical strong coupling scale of the EFT, and it is likely that the effective description breaks down before even reaching that scale. We discuss how this leaves the possibility that an appropriate ultraviolet completion of DE scenarios, valid at scales beyond an EFT description, can avoid present constraints on the GW speed. Instead, additional constraints in the lower frequency LISA band would be harder to escape, since the energies involved are orders of magnitude lower. By implementing a method based on GW multiband detections, we show indeed that a single joint observation of a GW150914-like event by LISA and a terrestrial interferometer would allow one to constrain the speed of light and gravitons to match to within $10^{-15}$. Multiband GW observations can therefore firmly constrain scenarios based on the EFT of DE, in a robust and unambiguous way. △ Less

Submitted 27 March, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 11 pages, 4 figures. Updated to match the version published in JCAP

Journal ref: JCAP03(2023)044

Implementation of an efficient Bayesian search for gravitational wave bursts with memory in pulsar timing array data

Authors: Jerry Sun, Paul T. Baker, Aaron D. Johnson, Dustin R. Madison, Xavier Siemens

Abstract: The standard Bayesian technique for searching pulsar timing data for gravitational wave (GW) bursts with memory (BWMs) using Markov Chain Monte Carlo (MCMC) sampling is very computationally expensive to perform. In this paper, we explain the implementation of an efficient Bayesian technique for searching for BWMs. This technique makes use of the fact that the signal model for Earth-term BWMs (BWMs… ▽ More The standard Bayesian technique for searching pulsar timing data for gravitational wave (GW) bursts with memory (BWMs) using Markov Chain Monte Carlo (MCMC) sampling is very computationally expensive to perform. In this paper, we explain the implementation of an efficient Bayesian technique for searching for BWMs. This technique makes use of the fact that the signal model for Earth-term BWMs (BWMs passing over the Earth) is fully factorizable. We estimate that this implementation reduces the computational complexity by a factor of 100. We also demonstrate that this technique gives upper limits consistent with published results using the standard Bayesian technique, and may be used to perform all of the same analyses that standard MCMC techniques can perform. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 19 pages, 3 figures, 1 table. Submitted to Astrophysical Journal

$\texttt{Hi-COLA}$: Fast, approximate simulations of structure formation in Horndeski gravity

Authors: Bill S. Wright, Ashim Sen Gupta, Tessa Baker, Georgios Valogiannis, Bartolomeo Fiorini

Abstract: We introduce $\texttt{Hi-COLA}$, a code designed to run fast, approximate $\textit{N}$-body simulations of non-linear structure formation in reduced Horndeski gravity. Given an input Lagrangian, $\texttt{Hi-COLA}$ dynamically constructs the appropriate field equations and consistently solves for the cosmological background, linear growth, and screened fifth force of that theory. Hence… ▽ More We introduce $\texttt{Hi-COLA}$, a code designed to run fast, approximate $\textit{N}$-body simulations of non-linear structure formation in reduced Horndeski gravity. Given an input Lagrangian, $\texttt{Hi-COLA}$ dynamically constructs the appropriate field equations and consistently solves for the cosmological background, linear growth, and screened fifth force of that theory. Hence $\texttt{Hi-COLA}$ is a general, adaptable, and useful tool that allows the mildly non-linear regime of many Horndeski theories to be investigated for the first time, at low computational cost. In this work, we first describe the screening approximations and simulation setup of $\texttt{Hi-COLA}$ for theories with Vainshtein screening. We validate the code against traditional $\textit{N}$-body simulations for cubic Galileon gravity, finding $2.5\%$ agreement up to $k_{\rm max}=1.2~h/{\rm Mpc}$. To demonstrate the flexibility of $\texttt{Hi-COLA}$, we additionally run the first simulations of an extended shift-symmetric gravity theory. We use the consistency and modularity of $\texttt{Hi-COLA}$ to dissect how the modified background, linear growth, and screened fifth force all contribute to departures from $Λ$CDM in the non-linear matter power spectrum. $\texttt{Hi-COLA}$ can be found at https://github.com/Hi-COLACode/Hi-COLA . △ Less

Submitted 22 March, 2023; v1 submitted 4 September, 2022; originally announced September 2022.

Comments: 51 pages, 20 figures, 5 tables, comments welcome. Updated to match the version published in JCAP

Journal ref: JCAP03(2023)040

Optimized Laser Models with Heisenberg-Limited Coherence and Sub-Poissonian Beam Photon Statistics

Authors: L. A. Ostrowski, T. J. Baker, S. N. Saadatmand, H. M. Wiseman

Abstract: Recently it has been shown that it is possible for a laser to produce a stationary beam with a coherence (quantified as the mean photon number at spectral peak) which scales as the fourth power of the mean number of excitations stored within the laser, this being quadratically larger than the standard or Schawlow-Townes limit [1]. Moreover, this was analytically proven to be the ultimate quantum l… ▽ More Recently it has been shown that it is possible for a laser to produce a stationary beam with a coherence (quantified as the mean photon number at spectral peak) which scales as the fourth power of the mean number of excitations stored within the laser, this being quadratically larger than the standard or Schawlow-Townes limit [1]. Moreover, this was analytically proven to be the ultimate quantum limit (Heisenberg limit) scaling under defining conditions for CW lasers, plus a strong assumption about the properties of the output beam. In Ref. [2], we show that the latter can be replaced by a weaker assumption, which allows for highly sub-Poissonian output beams, without changing the upper bound scaling or its achievability. In this Paper, we provide details of the calculations in Ref. [2], and introduce three new families of laser models which may be considered as generalizations of those presented in that work. Each of these families of laser models is parameterized by a real number, $p$, with $p=4$ corresponding to the original models. The parameter space of these laser families is numerically investigated in detail, where we explore the influence of these parameters on both the coherence and photon statistics of the laser beams. Two distinct regimes for the coherence may be identified based on the choice of $p$, where for $p>3$, each family of models exhibits Heisenberg-limited beam coherence, while for $p<3$, the Heisenberg limit is no longer attained. Moreover, in the former regime, we derive formulae for the beam coherence of each of these three laser families which agree with the numerics. We find that the optimal parameter is in fact $p\approx4.15$, not $p=4$. △ Less

Submitted 1 May, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: This is a companion manuscript to the letter entitled "No Tradeoff between Coherence and Sub-Poissonianity for Heisenberg-Limited Lasers", arxiv:2208.14081. 22 pages, 11 figures

Journal ref: Physical Review A 107, 053702 (2023)

No Tradeoff between Coherence and Sub-Poissonianity for Heisenberg-Limited Lasers

Authors: L. A. Ostrowski, T. J. Baker, S. N. Saadatmand, H. M. Wiseman

Abstract: The Heisenberg limit to laser coherence $\mathfrak{C}$ -- the number of photons in the maximally populated mode of the laser beam -- is the fourth power of the number of excitations inside the laser. We generalize the previous proof of this upper bound scaling by dropping the requirement that the beam photon statistics be Poissonian (i.e., Mandel's $Q=0$). We then show that the relation between… ▽ More The Heisenberg limit to laser coherence $\mathfrak{C}$ -- the number of photons in the maximally populated mode of the laser beam -- is the fourth power of the number of excitations inside the laser. We generalize the previous proof of this upper bound scaling by dropping the requirement that the beam photon statistics be Poissonian (i.e., Mandel's $Q=0$). We then show that the relation between $\mathfrak{C}$ and sub-Poissonianity ($Q<0$) is win-win, not a tradeoff. For both regular (non-Markovian) pumping with semi-unitary gain (which allows $Q\xrightarrow{}-1$), and random (Markovian) pumping with optimized gain, $\mathfrak{C}$ is maximized when $Q$ is minimized. △ Less

Submitted 1 May, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: This is a companion letter to the manuscript entitled "Optimized Laser Models with Heisenberg-Limited Coherence and Sub-Poissonian Beam Photon Statistics", arxiv:2208.14082. 6 pages, 2 figures

Journal ref: Physical Review Letters 130, 183602 (2023)

State space model multiple imputation for missing data in non-stationary multivariate time series with application in digital Psychiatry

Authors: Xiaoxuan Cai, Xinru Wang, Li Zeng, Habiballah Rahimi Eichi, Dost Ongur, Lisa Dixon, Justin T. Baker, Jukka-Pekka Onnela, Linda Valeri

Abstract: Mobile technology enables unprecedented continuous monitoring of an individual's behavior, social interactions, symptoms, and other health conditions, presenting an enormous opportunity for therapeutic advancements and scientific discoveries regarding the etiology of psychiatric illness. Continuous collection of mobile data results in the generation of a new type of data: entangled multivariate ti… ▽ More Mobile technology enables unprecedented continuous monitoring of an individual's behavior, social interactions, symptoms, and other health conditions, presenting an enormous opportunity for therapeutic advancements and scientific discoveries regarding the etiology of psychiatric illness. Continuous collection of mobile data results in the generation of a new type of data: entangled multivariate time series of outcome, exposure, and covariates. Missing data is a pervasive problem in biomedical and social science research, and the Ecological Momentary Assessment (EMA) using mobile devices in psychiatric research is no exception. However, the complex structure of multivariate time series introduces new challenges in handling missing data for proper causal inference. Data imputation is commonly recommended to enhance data utility and estimation efficiency. The majority of available imputation methods are either designed for longitudinal data with limited follow-up times or for stationary time series, which are incompatible with potentially non-stationary time series. In the field of psychiatry, non-stationary data are frequently encountered as symptoms and treatment regimens may experience dramatic changes over time. To address missing data in possibly non-stationary multivariate time series, we propose a novel multiple imputation strategy based on the state space model (SSMmp) and a more computationally efficient variant (SSMimpute). We demonstrate their advantages over other widely used missing data strategies by evaluating their theoretical properties and empirical performance in simulations of both stationary and non-stationary time series, subject to various missing mechanisms. We apply the SSMimpute to investigate the association between social network size and negative mood using a multi-year observational smartphone study of bipolar patients, controlling for confounding variables. △ Less

Submitted 12 April, 2023; v1 submitted 28 June, 2022; originally announced June 2022.

Cosmology with the Laser Interferometer Space Antenna

Authors: Pierre Auclair, David Bacon, Tessa Baker, Tiago Barreiro, Nicola Bartolo, Enis Belgacem, Nicola Bellomo, Ido Ben-Dayan, Daniele Bertacca, Marc Besancon, Jose J. Blanco-Pillado, Diego Blas, Guillaume Boileau, Gianluca Calcagni, Robert Caldwell, Chiara Caprini, Carmelita Carbone, Chia-Feng Chang, Hsin-Yu Chen, Nelson Christensen, Sebastien Clesse, Denis Comelli, Giuseppe Congedo, Carlo Contaldi, Marco Crisostomi , et al. (155 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations exten… ▽ More The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Report number: LISA CosWG-22-03

AI for Next Generation Computing: Emerging Trends and Future Directions

Authors: Sukhpal Singh Gill, Minxian Xu, Carlo Ottaviani, Panos Patros, Rami Bahsoon, Arash Shaghaghi, Muhammed Golec, Vlado Stankovski, Huaming Wu, Ajith Abraham, Manmeet Singh, Harshit Mehta, Soumya K. Ghosh, Thar Baker, Ajith Kumar Parlikad, Hanan Lutfiyya, Salil S. Kanhere, Rizos Sakellariou, Schahram Dustdar, Omer Rana, Ivona Brandic, Steve Uhlig

Abstract: Autonomic computing investigates how systems can achieve (user) specified control outcomes on their own, without the intervention of a human operator. Autonomic computing fundamentals have been substantially influenced by those of control theory for closed and open-loop systems. In practice, complex systems may exhibit a number of concurrent and inter-dependent control loops. Despite research into… ▽ More Autonomic computing investigates how systems can achieve (user) specified control outcomes on their own, without the intervention of a human operator. Autonomic computing fundamentals have been substantially influenced by those of control theory for closed and open-loop systems. In practice, complex systems may exhibit a number of concurrent and inter-dependent control loops. Despite research into autonomic models for managing computer resources, ranging from individual resources (e.g., web servers) to a resource ensemble (e.g., multiple resources within a data center), research into integrating Artificial Intelligence (AI) and Machine Learning (ML) to improve resource autonomy and performance at scale continues to be a fundamental challenge. The integration of AI/ML to achieve such autonomic and self-management of systems can be achieved at different levels of granularity, from full to human-in-the-loop automation. In this article, leading academics, researchers, practitioners, engineers, and scientists in the fields of cloud computing, AI/ML, and quantum computing join to discuss current research and potential future directions for these fields. Further, we discuss challenges and opportunities for leveraging AI and ML in next generation computing for emerging computing paradigms, including cloud, fog, edge, serverless and quantum computing environments. △ Less

Submitted 5 March, 2022; originally announced March 2022.

Comments: Accepted for Publication in Elsevier IoT Journal, 2022

Measuring the propagation speed of gravitational waves with LISA

Authors: Tessa Baker, Gianluca Calcagni, Anson Chen, Matteo Fasiello, Lucas Lombriser, Katarina Martinovic, Mauro Pieroni, Mairi Sakellariadou, Gianmassimo Tasinato, Daniele Bertacca, Ippocratis D. Saltas

Abstract: The propagation speed of gravitational waves, $c_T$, has been tightly constrained by the binary neutron star merger GW170817 and its electromagnetic counterpart, under the assumption of a frequency-independent $c_T$. Drawing upon arguments from Effective Field Theory and quantum gravity, we discuss the possibility that modifications of General Relativity allow for transient deviations of $c_T$ fro… ▽ More The propagation speed of gravitational waves, $c_T$, has been tightly constrained by the binary neutron star merger GW170817 and its electromagnetic counterpart, under the assumption of a frequency-independent $c_T$. Drawing upon arguments from Effective Field Theory and quantum gravity, we discuss the possibility that modifications of General Relativity allow for transient deviations of $c_T$ from the speed of light at frequencies well below the band of current ground-based detectors. We motivate two representative Ansätze for $c_T(f)$, and study their impact upon the gravitational waveforms of massive black hole binary mergers detectable by the LISA mission. We forecast the constraints on $c_T(f)$ obtainable from individual systems and a population of sources, from both inspiral and a full inspiral-merger-ringdown waveform. We show that LISA will enable us to place stringent independent bounds on departures from General Relativity in unexplored low-frequency regimes, even in the absence of an electromagnetic counterpart. △ Less

Submitted 13 July, 2022; v1 submitted 1 March, 2022; originally announced March 2022.

Comments: 38+20 pages, 22 figures. V2: Analysis improved, to be published in JCAP

Report number: LISA CosWG-22-02