-
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.
-
PINT: Maximum-likelihood estimation of pulsar timing noise parameters
Authors:
Abhimanyu Susobhanan,
David Kaplan,
Anne Archibald,
Jing Luo,
Paul Ray,
Timothy Pennucci,
Scott Ransom,
Gabriella Agazie,
William Fiore,
Bjorn Larsen,
Patrick O'Neill,
Rutger van Haasteren,
Akash Anumarlapudi,
Matteo Bachetti,
Deven Bhakta,
Chloe Champagne,
H. Thankful Cromartie,
Paul Demorest,
Ross Jennings,
Matthew Kerr,
Sasha Levina,
Alexander McEwen,
Brent Shapiro-Albert,
Joseph Swiggum
Abstract:
PINT is a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework within PINT to characterize the single-pulsar noise processes present in pulsar timing datasets. This framework enables the parameter estimation for both uncorr…
▽ More
PINT is a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework within PINT to characterize the single-pulsar noise processes present in pulsar timing datasets. This framework enables the parameter estimation for both uncorrelated and correlated noise processes as well as the model comparison between different timing and noise models in a computationally inexpensive way. We demonstrate the efficacy of the new framework by applying it to simulated datasets as well as a real dataset of PSR B1855+09. We also describe the new features implemented in PINT since it was first described in the literature.
△ Less
Submitted 19 June, 2024; v1 submitted 3 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.
-
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.
-
The SpacePy space science package at 12 years
Authors:
Jonathan T. Niehof,
Steven K. Morley,
Daniel T. Welling,
Brian A. Larsen
Abstract:
For over a decade, the SpacePy project has contributed open-source solutions for the production and analysis of heliophysics data and simulation results. Here we introduce SpacePy's functionality for the scientific user and present relevant design principles. We examine recent advances and the future of SpacePy in the broader scientific Python ecosystem, concluding with some of the work that has u…
▽ More
For over a decade, the SpacePy project has contributed open-source solutions for the production and analysis of heliophysics data and simulation results. Here we introduce SpacePy's functionality for the scientific user and present relevant design principles. We examine recent advances and the future of SpacePy in the broader scientific Python ecosystem, concluding with some of the work that has used SpacePy.
△ Less
Submitted 9 October, 2022; v1 submitted 22 August, 2022;
originally announced August 2022.
-
Bibliometric Indicators of Young Authors in Astrophysics: Can Later Stars be Predicted?
Authors:
Frank Havemann,
Birger Larsen
Abstract:
We test 16 bibliometric indicators with respect to their validity at the level of the individual researcher by estimating their power to predict later successful researchers. We compare the indicators of a sample of astrophysics researchers who later co-authored highly cited papers before their first landmark paper with the distributions of these indicators over a random control group of young aut…
▽ More
We test 16 bibliometric indicators with respect to their validity at the level of the individual researcher by estimating their power to predict later successful researchers. We compare the indicators of a sample of astrophysics researchers who later co-authored highly cited papers before their first landmark paper with the distributions of these indicators over a random control group of young authors in astronomy and astrophysics. We find that field and citation-window normalisation substantially improves the predicting power of citation indicators. The two indicators of total influence based on citation numbers normalised with expected citation numbers are the only indicators which show differences between later stars and random authors significant on a 1% level. Indicators of paper output are not very useful to predict later stars. The famous $h$-index makes no difference at all between later stars and the random control group.
△ Less
Submitted 11 April, 2014;
originally announced April 2014.
-
Mixed dark matter with low-mass bosons
Authors:
Georg B. Larsen,
Jes Madsen
Abstract:
We calculate the linear power spectrum for a range of mixed dark matter (MDM) models assuming a massive (few eV) boson, $φ$, instead of a neutrino as the hot component. We consider both the case where the hot dark matter (HDM) particle is a boson and the cold component is some other unknown particle, and the case where there is only one dark matter particle, a boson, with the cold dark matter (C…
▽ More
We calculate the linear power spectrum for a range of mixed dark matter (MDM) models assuming a massive (few eV) boson, $φ$, instead of a neutrino as the hot component. We consider both the case where the hot dark matter (HDM) particle is a boson and the cold component is some other unknown particle, and the case where there is only one dark matter particle, a boson, with the cold dark matter (CDM) component in a Bose condensate. Models resembling the latter type could arise from neutrino decays - we discuss some variants of this idea. The power spectra for MDM models with massive bosons are almost identical to neutrino MDM models for a given mass fraction of HDM if the bosons are distinct from their antiparticles ($φ\neq\barφ$) and have a temperature like that of neutrinos, whereas models with $φ=\barφ$ tend to overproduce small-scale structure.
△ Less
Submitted 24 January, 1996;
originally announced January 1996.
-
Mixed dark matter with neutrino chemical potentials
Authors:
Georg B. Larsen,
Jes Madsen
Abstract:
Mixed dark matter models with one low-mass (e.g. 2.4eV) neutrino flavor are shown to be in good agreement with observations if the neutrinos have non-zero chemical potentials. This agreement holds (except for the problem with a low-age Universe) even for high values of the Hubble-parameter. Massless neutrinos with non-zero chemical potentials may reconcile cold dark matter with observations. Som…
▽ More
Mixed dark matter models with one low-mass (e.g. 2.4eV) neutrino flavor are shown to be in good agreement with observations if the neutrinos have non-zero chemical potentials. This agreement holds (except for the problem with a low-age Universe) even for high values of the Hubble-parameter. Massless neutrinos with non-zero chemical potentials may reconcile cold dark matter with observations. Some fine-tuning is required to avoid problems with Big Bang nucleosynthesis.
△ Less
Submitted 24 July, 1995; v1 submitted 9 February, 1995;
originally announced February 1995.