NANOGrav 15-year gravitational-wave background methods

2024-05-01
Johnson, Aaron D.
Meyers, Patrick M.
Baker, Paul T.
Cornish, Neil J.
Hazboun, Jeffrey S.
Littenberg, Tyson B.
Romano, Joseph D.
Taylor, Stephen R.
Vallisneri, Michele
Vigeland, Sarah J.
Olum, Ken D.
Siemens, Xavier
Ellis, Justin A.
van Haasteren, Rutger
Hourihane, Sophie
Agazie, Gabriella
Anumarlapudi, Akash
Archibald, Anne M.
Arzoumanian, Zaven
Blecha, Laura
Brazier, Adam
Brook, Paul R.
Burke-Spolaor, Sarah
Becsy, Bence
Casey-Clyde, J. Andrew
Charisi, Maria
Chatterjee, Shami
Chatziioannou, Katerina
Cohen, Tyler
Cordes, James M.
Crawford, Fronefield
Cromartie, H. Thankful
Crowter, Kathryn
DeCesar, Megan E.
Demorest, Paul B.
Dolch, Timothy
Drachler, Brendan
Ferrara, Elizabeth C.
Fiore, William
Fonseca, Emmanuel
Freedman, Gabriel E.
Garver-Daniels, Nate
Gentile, Peter A.
Glaser, Joseph
Good, Deborah C.
Gultekin, Kayhan
Jennings, Ross J.
Jones, Megan L.
Kaiser, Andrew R.
Kaplan, David L.
Kelley, Luke Zoltan
Kerr, Matthew
Key, Joey S.
Laal, Nima
Lam, Michael T.
Lamb, William G.
Lazio, T. Joseph W.
Lewandowska, Natalia
Liu, Tingting
Lorimer, Duncan R.
Luo, Jing
Lynch, Ryan S.
Ma, Chung-Pei
Madison, Dustin R.
McEwen, Alexander
McKee, James W.
McLaughlin, Maura A.
McMann, Natasha
Meyers, Bradley W.
Mingarelli, Chiara M. F.
Mitridate, Andrea
Ng, Cherry
Nice, David J.
Ocker, Stella Koch
Pennucci, Timothy T.
Perera, Benetge B. P.
Pol, Nihan S.
Radovan, Henri A.
Ransom, Scott M.
Ray, Paul S.
Sardesai, Shashwat C.
Schmiedekamp, Carl
Schmiedekamp, Ann
Schmitz, Kai
Shapiro-Albert, Brent J.
Simon, Joseph
Siwek, Magdalena S.
Stairs, Ingrid H.
Stinebring, Daniel R.
Stovall, Kevin
Susobhanan, Abhimanyu
Swiggum, Joseph K.
Turner, Jacob E.
Ünal, Caner
Wahl, Haley M.
Witt, Caitlin A.
Young, Olivia
Pulsar timing arrays (PTAs) use an array of millisecond pulsars to search for gravitational waves in the nanohertz regime in pulse time of arrival data. This paper presents rigorous tests of PTA methods, examining their consistency across the relevant parameter space. 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 Enterprise and 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 in gravitational wave searches, reducing the required resources facilitating the computation of Bayes factors via thermodynamic integration and sampling a large number of realizations for computing Bayesian false-alarm probabilities. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis methods. For the Bayesian analysis, we test prior recovery, simulation recovery, and Bayes factors. For the frequentist analysis, we test that the 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 gravitational wave background detection and parameter estimation perspectives. The tests presented here validate current analyses of PTA data.
PHYSICAL REVIEW D
Citation Formats
A. D. Johnson et al., “NANOGrav 15-year gravitational-wave background methods,” PHYSICAL REVIEW D, vol. 109, no. 10, pp. 0–0, 2024, Accessed: 00, 2025. [Online]. Available: https://hdl.handle.net/11511/116550.