-
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.
-
LISA Galactic Binaries in the Roman Galactic Bulge Time-Domain Survey
Authors:
Matthew C. Digman,
Chris M. Hirata
Abstract:
Short-period Galactic white dwarf binaries detectable by LISA are the only guaranteed persistent sources for multi-messenger gravitational-wave astronomy. Large-scale surveys in the 2020s present an opportunity to conduct preparatory science campaigns to maximize the science yield from future multi-messenger targets. The Nancy Grace Roman Space Telescope Galactic Bulge Time Domain Survey will (in…
▽ More
Short-period Galactic white dwarf binaries detectable by LISA are the only guaranteed persistent sources for multi-messenger gravitational-wave astronomy. Large-scale surveys in the 2020s present an opportunity to conduct preparatory science campaigns to maximize the science yield from future multi-messenger targets. The Nancy Grace Roman Space Telescope Galactic Bulge Time Domain Survey will (in its Reference Survey design) image seven fields in the Galactic Bulge approximately 40,000 times each. Although the Reference Survey cadence is optimized for detecting exoplanets via microlensing, it is also capable of detecting short-period white dwarf binaries. In this paper, we present forecasts for the number of detached short-period binaries the Roman Galactic Bulge Time Domain Survey will discover and the implications for the design of electromagnetic surveys. Although population models are highly uncertain, we find a high probability that the baseline survey will detect of order ~5 detached white dwarf binaries. The Reference Survey would also have a $\gtrsim20\%$ chance of detecting several known benchmark white dwarf binaries at the distance of the Galactic Bulge.
△ Less
Submitted 30 December, 2022;
originally announced December 2022.
-
Parameter Estimation for Stellar-Origin Black Hole Mergers In LISA
Authors:
Matthew C. Digman,
Neil J. Cornish
Abstract:
The population of stellar origin black hole binaries (SOBHBs) detected by existing ground-based gravitational wave detectors is an exciting target for the future space-based Laser Interferometer Space Antenna (LISA). LISA is sensitive to signals at significantly lower frequencies than ground-based detectors. SOBHB signals will thus be detected much earlier in their evolution, years to decades befo…
▽ More
The population of stellar origin black hole binaries (SOBHBs) detected by existing ground-based gravitational wave detectors is an exciting target for the future space-based Laser Interferometer Space Antenna (LISA). LISA is sensitive to signals at significantly lower frequencies than ground-based detectors. SOBHB signals will thus be detected much earlier in their evolution, years to decades before they merge. The mergers will then occur in the frequency band covered by ground-based detectors. Observing SOBHBs years before merger can help distinguish between progenitor models for these systems. We present a new Bayesian parameter estimation algorithm for LISA observations of SOBHBs that uses a time-frequency (wavelet) based likelihood function. Our technique accelerates the analysis by several orders of magnitude compared to the standard frequency domain approach and allows for an efficient treatment of non-stationary noise.
△ Less
Submitted 8 December, 2022;
originally announced December 2022.
-
LISA Gravitational Wave Sources in A Time-Varying Galactic Stochastic Background
Authors:
Matthew C. Digman,
Neil J. Cornish
Abstract:
A unique challenge for data analysis with the Laser Interferometer Space Antenna (LISA) is that the noise backgrounds from instrumental noise and astrophysical sources will change significantly over both the year and the entire mission. Variations in the noise levels will be on time scales comparable to, or shorter than, the time most signals spend in the detector's sensitive band. The variation i…
▽ More
A unique challenge for data analysis with the Laser Interferometer Space Antenna (LISA) is that the noise backgrounds from instrumental noise and astrophysical sources will change significantly over both the year and the entire mission. Variations in the noise levels will be on time scales comparable to, or shorter than, the time most signals spend in the detector's sensitive band. The variation in the amplitude of the galactic stochastic GW background from galactic binaries as the antenna pattern rotates relative to the galactic center is a particularly significant component of the noise variation. LISA's sensitivity to different source classes will therefore vary as a function of sky location and time. The variation will impact both overall signal-to-noise and the efficiency of alerts to EM observers to search for multi-messenger counterparts.
△ Less
Submitted 29 June, 2022;
originally announced June 2022.
-
Fast Bayesian analysis of individual binaries in pulsar timing array data
Authors:
Bence Bécsy,
Neil J. Cornish,
Matthew C. Digman
Abstract:
Searching for gravitational waves in pulsar timing array data is computationally intensive. The data is unevenly sampled, and the noise is heteroscedastic, necessitating the use of a time-domain likelihood function with attendant expensive matrix operations. The computational cost is exacerbated when searching for individual supermassive black hole binaries, which have a large parameter space due…
▽ More
Searching for gravitational waves in pulsar timing array data is computationally intensive. The data is unevenly sampled, and the noise is heteroscedastic, necessitating the use of a time-domain likelihood function with attendant expensive matrix operations. The computational cost is exacerbated when searching for individual supermassive black hole binaries, which have a large parameter space due to the additional pulsar distance, phase offset and noise model parameters needed for each pulsar. We introduce a new formulation of the likelihood function which can be used to make the Bayesian analysis significantly faster. We divide the parameters into projection and shape parameters. We then accelerate the exploration of the projection parameters by more than four orders of magnitude by precomputing the expensive inner products for each set of shape parameters. The projection parameters include nuisance parameters such as the gravitational wave phase offset at each pulsar. In the new scheme, these troublesome nuisance parameters are efficiently marginalized over using multiple-try Markov chain Monte Carlo sampling as part of a Metropolis-within-Gibbs scheme. The acceleration provided by our method will become increasingly important as pulsar timing datasets rapidly grow. Our method also makes sophisticated analyses more tractable, such as searches for multiple binaries, or binaries with non-negligible eccentricities.
△ Less
Submitted 20 June, 2022; v1 submitted 14 April, 2022;
originally announced April 2022.
-
COSMIC Variance in Binary Population Synthesis
Authors:
Katelyn Breivik,
Scott Coughlin,
Michael Zevin,
Carl L. Rodriguez,
Kyle Kremer,
Claire S. Ye,
Jeff J. Andrews,
Michael Kurkowski,
Matthew C. Digman,
Shane L. Larson,
Frederic A. Rasio
Abstract:
The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observat…
▽ More
The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite: COSMIC which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signal observable by the Laser Interferometer Space Antenna (LISA). We find that $\sim10^8$ compact binaries reside in the Milky Way today, while $\sim10^4$ of them may be resolvable by LISA.
△ Less
Submitted 10 September, 2020; v1 submitted 3 November, 2019;
originally announced November 2019.
-
(Not as) Big as a Barn: Upper Bounds on Dark Matter-Nucleus Cross Sections
Authors:
Matthew C. Digman,
Christopher V. Cappiello,
John F. Beacom,
Christopher M. Hirata,
Annika H. G. Peter
Abstract:
Critical probes of dark matter come from tests of its elastic scattering with nuclei. The results are typically assumed to be model-independent, meaning that the form of the potential need not be specified and that the cross sections on different nuclear targets can be simply related to the cross section on nucleons. For point-like spin-independent scattering, the assumed scaling relation is…
▽ More
Critical probes of dark matter come from tests of its elastic scattering with nuclei. The results are typically assumed to be model-independent, meaning that the form of the potential need not be specified and that the cross sections on different nuclear targets can be simply related to the cross section on nucleons. For point-like spin-independent scattering, the assumed scaling relation is $σ_{χA} \propto A^2 μ_A^2 σ_{χN}\propto A^4 σ_{χN}$, where the $A^2$ comes from coherence and the $μ_A^2\simeq A^2 m_N^2$ from kinematics for $m_χ\gg m_A$. Here we calculate where model independence ends, i.e., where the cross section becomes so large that it violates its defining assumptions. We show that the assumed scaling relations generically fail for dark matter-nucleus cross sections $σ_{χA} \sim 10^{-32}-10^{-27}\;\text{cm}^2$, significantly below the geometric sizes of nuclei, and well within the regime probed by underground detectors. Last, we show on theoretical grounds, and in light of existing limits on light mediators, that point-like dark matter cannot have $σ_{χN}\gtrsim10^{-25}\;\text{cm}^2$, above which many claimed constraints originate from cosmology and astrophysics. The most viable way to have such large cross sections is composite dark matter, which introduces significant additional model dependence through the choice of form factor. All prior limits on dark matter with cross sections $σ_{χN}>10^{-32}\;\text{cm}^2$ with $m_χ\gtrsim 1\;\text{GeV}$ must therefore be re-evaluated and reinterpreted.
△ Less
Submitted 30 December, 2022; v1 submitted 24 July, 2019;
originally announced July 2019.
-
Forecasting Super-Sample Covariance in Future Weak Lensing Surveys with SuperSCRAM
Authors:
Matthew C. Digman,
Joseph E. McEwen,
Christopher M. Hirata
Abstract:
The observable universe contains density perturbations on scales larger than any finite volume survey. Perturbations on scales larger than a survey can measure degrade its power to constrain cosmological parameters. The dependence of survey observables such as the weak lensing power spectrum on these long-wavelength modes results in super-sample covariance. Accurately forecasting parameter constra…
▽ More
The observable universe contains density perturbations on scales larger than any finite volume survey. Perturbations on scales larger than a survey can measure degrade its power to constrain cosmological parameters. The dependence of survey observables such as the weak lensing power spectrum on these long-wavelength modes results in super-sample covariance. Accurately forecasting parameter constraints for future surveys requires accurately accounting for the super-sample effects. If super-sample covariance is in fact a major component of the survey error budget, it may be necessary to investigate mitigation strategies that constrain the specific realization of the long-wavelength modes. We present a Fisher matrix based formalism for approximating the magnitude of super-sample covariance and the effectiveness of mitigation strategies for realistic survey geometries. We implement our formalism in the public code SuperSCRAM: Super-Sample Covariance Reduction and Mitigation. We illustrate SuperSCRAM with an example application, where the modes contributing to super-sample covariance in the WFIRST weak lensing survey are constrained by the low-redshift galaxy number counts in the wider LSST footprint. We find that super-sample covariance increases the volume of the error ellipsoid in 7D cosmological parameter space by a factor of 4.5 relative to Gaussian statistical errors only, but our simple mitigation strategy more than halves the contamination, to a factor of 2.0.
△ Less
Submitted 26 April, 2019;
originally announced April 2019.