Showing 1–28 of 28 results for author: Siwek, M

Eccentric Signatures of Stellar-Mass Binary Black Holes with Circumbinary Disks in LISA

Authors: Isobel M. Romero-Shaw, Samir Goorachurn, Magdalena Siwek, Christopher J. Moore

Abstract: Binary black holes may have circumbinary disks if formed through common-envelope evolution or within gaseous environments. Disks can drive binaries into wider and more eccentric orbits, while gravitational waves harden and circularise them. We combine cutting-edge evolution prescriptions for disk-driven binaries with well-known equations for gravitational-wave-driven evolution, and study the evolu… ▽ More Binary black holes may have circumbinary disks if formed through common-envelope evolution or within gaseous environments. Disks can drive binaries into wider and more eccentric orbits, while gravitational waves harden and circularise them. We combine cutting-edge evolution prescriptions for disk-driven binaries with well-known equations for gravitational-wave-driven evolution, and study the evolution of stellar-mass binary black holes. We find that binaries are driven by their disk to an equilibrium eccentricity, $0.2\lesssim e_\mathrm{eq} \lesssim0.5$, that dominates their evolution. Once they transition to the GW-dominated regime their eccentricity decreases rapidly; we find that binaries with long-lived disks will likely be observed in LISA with detectable eccentricities $\sim 10^{-2}$ at $0.01$ Hz, with the precise value closely correlating with the binary's initial mass ratio. This may lead binary black holes with CBDs observed in LISA to be confused with dynamically-formed binary black holes. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 7 pages, 3 figures. Submitted to MNRAS Letters. Comments welcome!

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

Signatures of Circumbinary Disk Dynamics in Multi-Messenger Population Studies of Massive Black Hole Binaries

Authors: Magdalena Siwek, Luke Zoltan Kelley, Lars Hernquist

Abstract: We investigate the effect of cutting-edge circumbinary disk (CBD) evolution models on massive black hole binary (MBHB) populations and the gravitational wave background (GWB). We show that CBD-driven evolution leaves a tell-tale signature in MBHB populations, by driving binaries towards an equilibrium eccentricity that depends on binary mass ratio. We find high orbital eccentricities (… ▽ More We investigate the effect of cutting-edge circumbinary disk (CBD) evolution models on massive black hole binary (MBHB) populations and the gravitational wave background (GWB). We show that CBD-driven evolution leaves a tell-tale signature in MBHB populations, by driving binaries towards an equilibrium eccentricity that depends on binary mass ratio. We find high orbital eccentricities ($e_{\rm b} \sim 0.5$) as MBHBs enter multi-messenger observable frequency bands. The CBD-induced eccentricity distribution of MBHB populations in observable bands is independent of the initial eccentricity distribution at binary formation, erasing any memory of eccentricities induced in the large-scale dynamics of merging galaxies. Our results suggest that eccentric MBHBs are the rule rather than the exception in upcoming transient surveys, provided that CBDs regularly form in MBHB systems. We show that the GWB amplitude is sensitive to CBD-driven preferential accretion onto the secondary, resulting in an increase in GWB amplitude $A_{\rm yr^{-1}}$ by over 100\% with just 10\% Eddington accretion. As we self consistently allow for binary hardening and softening, we show that CBD-driven orbital expansion does not diminish the GWB amplitude, and instead increases the amplitude by a small amount. We further present detection rates and population statistics of MBHBs with $M_{\rm b} \gtrsim 10^6 \, M_{\odot}$ in LISA, showing that most binaries have equal mass ratios and can retain residual eccentricities up to $e_{\rm b} \sim 10^{-3}$ due to CBD-driven evolution. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 13 pages, 11 figures. Submitted to MNRAS

The Santa Barbara Binary-Disk Code Comparison

Authors: Paul C. Duffell, Alexander J. Dittmann, Daniel J. D'Orazio, Alessia Franchini, Kaitlin M. Kratter, Anna B. T. Penzlin, Enrico Ragusa, Magdalena Siwek, Christopher Tiede, Haiyang Wang, Jonathan Zrake, Adam M. Dempsey, Zoltan Haiman, Alessandro Lupi, Michal Pirog, Geoffrey Ryan

Abstract: We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed… ▽ More We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, $0.2 (a^2 Ωおーむ_B /νにゅー)$ binary orbits, where $νにゅー$ is the kinematic viscosity. We did not investigate dependence on binary mass ratio, eccentricity, disk temperature, or disk viscosity; therefore, these benchmarks may act as guides towards expanding converged solutions to the wider parameter space but might need to be updated in a future study that investigates dependence on system parameters. We find the most major discrepancies between codes resulted from the dimensionality of the setup (3D vs 2D disks). Beyond this, we find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D. △ Less

Submitted 20 June, 2024; v1 submitted 20 February, 2024; originally announced February 2024.

Comments: ApJ Accepted

Long-term Evolution of Binary Orbits Induced by Circumbinary Disks

Authors: Ruggero Valli, Christopher Tiede, Alejandro Vigna-Gómez, Jorge Cuadra, Magdalena Siwek, Jing-Ze Ma, Daniel J. D'Orazio, Jonathan Zrake, Selma E. de Mink

Abstract: Circumbinary disks are found in a variety of astrophysical scenarios, spanning binary star formation to accreting supermassive black hole binaries. The interaction with a circumbinary disk can yield opposite effects on the binary orbit leading to circularization, or exciting the eccentricity, widening the orbit or shrinking it and facilitating mergers. We present a new formalism for the long-term… ▽ More Circumbinary disks are found in a variety of astrophysical scenarios, spanning binary star formation to accreting supermassive black hole binaries. The interaction with a circumbinary disk can yield opposite effects on the binary orbit leading to circularization, or exciting the eccentricity, widening the orbit or shrinking it and facilitating mergers. We present a new formalism for the long-term evolution of the disk-binary interaction based on the results of recent suites of hydrodynamic simulations, which resolve the complex geometry of the gas in the vicinity of the binary and fully account for the gravitational and accretion forces. We release a python package, \texttt{spindler}, that implements our model. We show that, unless the mass reservoir feeding the disk is comparable to the mass of the binary, accretion onto the binary depletes the disk mass before inducing a significant change in orbital separation or mass ratio. This finding implies that, in most scenarios, interaction with a circumbinary disk is not an efficient mechanism to shrink the orbit of the binary. However, as long as the mass of the disk is at least a few percent of the mass of the binary, the interaction can excite the eccentricity up to an equilibrium value, and induce a statistical correlation between mass ratio and eccentricity. We consider the applicability of our model to a variety of astrophysical scenarios: during star formation, in evolved stellar binaries, triples and in supermassive black hole binaries. We discuss the theoretical and observational implications of our predictions. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 18 pages, 12 figures. Submitted. Comments welcome. Spindler available as a Python package on Zenodo: https://doi.org/10.5281/zenodo.10529060

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 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)

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

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

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.

Orbital Evolution of Binaries in Circumbinary Disks

Authors: Magdalena Siwek, Rainer Weinberger, Lars Hernquist

Abstract: We present the to-date largest parameter space exploration of binaries in circumbinary disks (CBDs), deriving orbital evolution prescriptions for eccentric, unequal mass binaries from our suite of hydrodynamic simulations. In all cases, binary eccentricities evolve towards steady state values that increase with mass ratio, and saturate at an equilibrium eccentricity $e_{\rm b, eq} \sim 0.5$ in the… ▽ More We present the to-date largest parameter space exploration of binaries in circumbinary disks (CBDs), deriving orbital evolution prescriptions for eccentric, unequal mass binaries from our suite of hydrodynamic simulations. In all cases, binary eccentricities evolve towards steady state values that increase with mass ratio, and saturate at an equilibrium eccentricity $e_{\rm b, eq} \sim 0.5$ in the large mass ratio regime, in line with resonant theory. For binaries accreting at their combined Eddington limit, a steady state eccentricity can be achieved within a few Megayears. Once at their steady state eccentricities, binaries with $q_{\rm b} \gtrsim 0.3$ evolve towards coalescence, while lower mass ratio systems expand due to CBD torques. We discuss implications for population studies of massive black hole binaries, protostars in binary systems, and post-common envelope binaries observed by ground-based gravitational wave detectors. △ Less

Submitted 3 February, 2023; originally announced February 2023.

Comments: 12 pages, 8 figures. Submitted to MNRAS

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

Orbital and Radiative Properties of Wandering Intermediate-Mass Black Holes in the ASTRID Simulation

Authors: Emma Jane Weller, Fabio Pacucci, Yueying Ni, Nianyi Chen, Tiziana Di Matteo, Magdalena Siwek, Lars Hernquist

Abstract: Intermediate-Mass Black Holes (IMBHs) of $10^3-10^6 \, M_\odot$ are commonly found at the center of dwarf galaxies. Simulations and observations convincingly show that a sizable population of IMBHs could wander off-center in galaxies. We use the cosmological simulation ASTRID to study the orbital and radiative properties of wandering IMBHs in massive galaxies at $z\sim3$. We find that this populat… ▽ More Intermediate-Mass Black Holes (IMBHs) of $10^3-10^6 \, M_\odot$ are commonly found at the center of dwarf galaxies. Simulations and observations convincingly show that a sizable population of IMBHs could wander off-center in galaxies. We use the cosmological simulation ASTRID to study the orbital and radiative properties of wandering IMBHs in massive galaxies at $z\sim3$. We find that this population of black holes has large orbital inclinations ($60^\circ\pm22^\circ$) with respect to the principal plane of the host. The eccentricity of their orbits is also significant ($0.6\pm0.2$) and decreases with time. Wandering IMBHs undergo spikes of accretion activity around the pericenter of their orbits, with rates $10^{-3}-10^{-5}$ times the Eddington rate and a median accretion duty cycle of $\sim 12\%$. Their typical spectral energy distribution peaks in the infrared at $\sim 11 \, μみゅー\rm m$ rest-frame. Assuming a standard value of $10\%$ for the matter-to-energy radiative efficiency, IMBHs reach $2-10$ keV X-ray luminosities $>10^{37} \, \mathrm{erg\,s^{-1}}$ for $\sim10\%$ of the time. This luminosity corresponds to fluxes $>10^{-15} \, \mathrm{erg \, s^{-1} \, cm^{-2}}$ within $10$ Mpc. They could be challenging to detect because of competing emissions from X-ray binaries and the interstellar medium. X-ray luminosities $> 10^{41} \, \mathrm{erg \, s^{-1}}$, in the hyper-luminous X-ray sources (HLXs) regime, are reached by $\sim 7\%$ of the IMBHs. These findings suggest that HLXs are a small subset of the wandering IMBH population, which is characterized by luminosities $10^3-10^4$ times fainter. Dedicated surveys are needed to assess the demographics of this missing population of black holes. △ Less

Submitted 2 April, 2023; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS. This is the final version of the manuscript. 9 pages, 7 figures

Preferential Accretion and Circumbinary Disk Precession in Eccentric Binary Systems

Authors: Magdalena Siwek, Rainer Weinberger, Diego Munoz, Lars Hernquist

Abstract: We present a suite of high resolution hydrodynamic simulations of binaries immersed in circumbinary accretion disks (CBDs). For the first time, we investigate the preferential accretion rate as a function of both eccentricity $e_{\rm b}$ and mass ratio $q_{\rm b}$ in a densely sampled parameter space, finding that when compared with circular binaries, mass ratios grow more efficiently in binaries… ▽ More We present a suite of high resolution hydrodynamic simulations of binaries immersed in circumbinary accretion disks (CBDs). For the first time, we investigate the preferential accretion rate as a function of both eccentricity $e_{\rm b}$ and mass ratio $q_{\rm b}$ in a densely sampled parameter space, finding that when compared with circular binaries, mass ratios grow more efficiently in binaries on moderately eccentric orbits ($0.0 \lesssim e_{\rm b} \lesssim 0.4$), and high eccentricities ($e_{\rm b} \gtrsim 0.6 $) suppress mass ratio growth. We suggest that this non-monotonic preferential accretion behaviour may produce an observable shift in the mass ratio distributions of stellar binaries and massive black hole binaries. We further find that the response of a CBD can be divided into three regimes, depending on eccentricity and mass ratio: (i) CBDs around circular binaries always precess freely, whereas CBDs around eccentric binaries either (ii) undergo forced precession or (iii) remain locked at an angle with respect to the binary periapsis. Forced precession in eccentric binaries is associated with strong modulation of individual accretion rates on the precession timescale, a potentially observable signature in accreting binaries with short orbital periods. We provide CBD locking angles and precession rates as a function of $e_{\rm b}$ and $q_{\rm b}$ for our simulation suite. △ Less

Submitted 21 November, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 15 pages, 14 figures. Published in Monthly Notices of the Royal Astronomical Society

The NANOGrav 12.5-year data set: Search for Non-Einsteinian Polarization Modes in theGravitational-Wave Background

Authors: Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Bence Becsy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Dallas M. DeGan, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Nathan Garver-Daniels, Peter A. Gentile , et al. (46 additional authors not shown)

Abstract: We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (ten… ▽ More We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (tensor transverse, TT) correlations. Specifically, we find ST correlations with a signal-to-noise ratio of 2.8 that are preferred over TT correlations (Hellings and Downs correlations) with Bayesian odds of about 20:1. However, the significance of ST correlations is reduced dramatically when we include modeling of the Solar System ephemeris systematics and/or remove pulsar J0030$+$0451 entirely from consideration. Even taking the nominal signal-to-noise ratios at face value, analyses of simulated data sets show that such values are not extremely unlikely to be observed in cases where only the usual TT modes are present in the GWB. In the absence of a detection of any polarization mode of gravity, we place upper limits on their amplitudes for a spectral index of $γがんま= 5$ and a reference frequency of $f_\text{yr} = 1 \text{yr}^{-1}$. Among the upper limits for eight general families of metric theories of gravity, we find the values of $A^{95\%}_{TT} = (9.7 \pm 0.4)\times 10^{-16}$ and $A^{95\%}_{ST} = (1.4 \pm 0.03)\times 10^{-15}$ for the family of metric spacetime theories that contain both TT and ST modes. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 24 pages, 18 figures, 3 appendices. Please send any comments/questions to Nima Laal (laaln@oregonstate.edu)

Journal ref: The Astrophysical Journal Letters, vol. 923, no. 2, p. L22, Dec. 2021

The effect of differential accretion on the Gravitational Wave Background and the present day MBH Binary population

Authors: Magdalena S. Siwek, Luke Zoltan Kelley, Lars Hernquist

Abstract: Massive black hole binaries (MBHBs) form as a consequence of galaxy mergers. However, it is still unclear whether they typically merge within a Hubble time, and how accretion may affect their evolution. These questions will be addressed by pulsar timing arrays (PTAs), which aim to detect the GW background (GWB) emitted by MBHBs during the last Myrs of inspiral. Here we investigate the influence of… ▽ More Massive black hole binaries (MBHBs) form as a consequence of galaxy mergers. However, it is still unclear whether they typically merge within a Hubble time, and how accretion may affect their evolution. These questions will be addressed by pulsar timing arrays (PTAs), which aim to detect the GW background (GWB) emitted by MBHBs during the last Myrs of inspiral. Here we investigate the influence of differential accretion on MBHB merger rates, chirp masses and the resulting GWB spectrum. We evolve a MBHB sample from the Illustris hydrodynamic cosmological simulation using semi-analytic models and for the first time self-consistently evolve their masses with binary accretion models. In all models, MBHBs coalesce with median total masses up to $1.5 \times 10^8 M_{\odot}$, up to $3-4$ times larger than in models neglecting accretion. In our model with the largest plausible impact, the median mass ratio of coalescing MBHBs increases by a factor $3.6$, the coalescence rate by $52.3\%$, and the GWB amplitude by a factor $4.0$, yielding a dimensionless GWB strain $A_{yr^{-1}} = 1 \times 10^{-15}$. Our model that favours accretion onto the primary MBH reduces the median mass ratio of coalescing MBHBs by a factor of $2.9$, and yields a GWB amplitude $A_{yr^{-1}} = 3.1 \times 10^{-16}$. This is nearly indistinguishable from our model neglecting accretion, despite higher MBHB masses at coalescence. \textbf{We further predict binary separation and mass ratio distributions of stalled MBHBs in the low-redshift universe, and find that these depend sensitively on binary accretion models. This presents the potential for combined EM and GW observational constraints on merger rates and accretion models of MBHB populations.} △ Less

Submitted 2 June, 2020; v1 submitted 18 May, 2020; originally announced May 2020.

Comments: Updated to appropriately reference previous work by Sesana & Barausse et al

Observational limits on the X-ray emission from the bubble nebula surrounding Ho IX X-1

Authors: Rajath Sathyaprakash, Timothy P. Roberts, Magdalena M. Siwek

Abstract: Optical and radio observations of shock-ionised bubble nebulae surrounding ultraluminous X-ray sources (ULXs) suggest that they are powered by jets or super-critical outflows presumably launched from the ULX accretion disc. Recent simulations of these systems have shown that the shocked wind can emit thermal X-rays with estimated luminosities $\sim 10^{36}$ erg s$^{-1}$. In this work, we investiga… ▽ More Optical and radio observations of shock-ionised bubble nebulae surrounding ultraluminous X-ray sources (ULXs) suggest that they are powered by jets or super-critical outflows presumably launched from the ULX accretion disc. Recent simulations of these systems have shown that the shocked wind can emit thermal X-rays with estimated luminosities $\sim 10^{36}$ erg s$^{-1}$. In this work, we investigated whether it is possible to detect and spatially resolve the X-ray emission from these systems using archival Chandra observations of the ULX Holmberg IX X-1. This source is an ideal target to study for two reasons: it is surrounded by an optical bubble nebula with a large spatial extent ($\sim$ 400 pc) that can easily be resolved with Chandra. Further, it has a hard X-ray continuum that is easily distinguishable from the expected soft thermal emission from the nebula. However, a spectral and photometric analysis on stacked Chandra observations of the source reveals that there is no strong evidence for an X-ray bubble associated with it, to a limiting luminosity of $\sim 2 \times 10^{36}$ erg s$^{-1}$. The detection of such X-ray nebulae may be possible with future X-ray missions such as Athena, which would provide useful constraints on the kinematics of the outflow. Finally, our observations also emphasise that the nebular emission does not contribute significantly to the residuals in the X-ray spectrum of the source, which are more likely to be linked to processes localised to the ULX. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: 9 pages, 6 figures; accepted for publication in MNRAS

Host Galaxies of Type Ic and Broad-lined Type Ic Supernovae from the Palomar Transient Factory: Implication for Jet Production

Authors: Maryam Modjaz, Federica B. Bianco, Magdalena Siwek, Shan Huang, Daniel A. Perley, David Fierroz, Yu-Qian Liu, Iair Arcavi, Avishay Gal-Yam, Nadia Blagorodnova, Bradley S. Cenko, Alexei V. Filippenko, Mansi M. Kasliwal, S. R. Kulkarni, Steve Schulze, Kirsty Taggart, Weikang Zhen

Abstract: Unlike the ordinary supernovae (SNe) some of which are hydrogen and helium deficient (called Type Ic SNe), broad-lined Type Ic SNe (SNe Ic-bl) are very energetic events, and all SNe coincident with bona fide long duration gamma-ray bursts (LGRBs) are of Type Ic-bl. Understanding the progenitors and the mechanism driving SN Ic-bl explosions vs those of their SNe Ic cousins is key to understanding t… ▽ More Unlike the ordinary supernovae (SNe) some of which are hydrogen and helium deficient (called Type Ic SNe), broad-lined Type Ic SNe (SNe Ic-bl) are very energetic events, and all SNe coincident with bona fide long duration gamma-ray bursts (LGRBs) are of Type Ic-bl. Understanding the progenitors and the mechanism driving SN Ic-bl explosions vs those of their SNe Ic cousins is key to understanding the SN-GRB relationship and jet production in massive stars. Here we present the largest set of host-galaxy spectra of 28 SNe Ic and 14 SN Ic-bl, all discovered before 2013 by the same untargeted survey, namely the Palomar Transient Factory (PTF). We carefully measure their gas-phase metallicities, stellar masses (M*s) and star-formation rates (SFRs) by taking into account recent progress in the metallicity field and propagating uncertainties correctly. We further re-analyze the hosts of 10 literature SN-GRBs using the same methods and compare them to our PTF SN hosts with the goal of constraining their progenitors from their local environments by conducting a thorough statistical comparison, including upper limits. We find that the metallicities, SFRs and M*s of our PTF SN Ic-bl hosts are statistically comparable to those of SN-GRBs, but significantly lower than those of the PTF SNe Ic. The mass-metallicity relations as defined by the SNe Ic-bl and SN-GRBs are not significantly different from the same relations as defined by the SDSS galaxies, in contrast to claims by earlier works. Our findings point towards low metallicity as a crucial ingredient for SN Ic-bl and SN-GRB production since we are able to break the degeneracy between high SFR and low metallicity. We suggest that the PTF SNe Ic-bl may have produced jets that were choked inside the star or were able break out of the star as unseen low-luminosity or off-axis GRBs. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: 47 pages, 8 figures, submitted to ApJ in Nov 2018

Optical and X-ray luminosity of expanding nebulae around ultraluminous X-ray sources

Authors: Magdalena Siwek, Aleksander Sadowski, Ramesh Narayan, Timothy P. Roberts, Roberto Soria

Abstract: We have performed a set of simulations of expanding, spherically symmetric nebulae inflated by winds from accreting black holes in ultraluminous X-ray sources (ULXs). We implemented a realistic cooling function to account for free-free and bound-free cooling. For all model parameters we considered, the forward shock in the interstellar medium becomes radiative at a radius $\sim $ 100 pc. The emiss… ▽ More We have performed a set of simulations of expanding, spherically symmetric nebulae inflated by winds from accreting black holes in ultraluminous X-ray sources (ULXs). We implemented a realistic cooling function to account for free-free and bound-free cooling. For all model parameters we considered, the forward shock in the interstellar medium becomes radiative at a radius $\sim $ 100 pc. The emission is primarily in the optical and UV, and the radiative luminosity is about 50% of the total kinetic luminosity of the wind. In contrast, the reverse shock in the wind is adiabatic so long as the terminal outflow velocity of the wind $v_{\rm w} \sim 0.003c$. The shocked wind in these models radiates in X-rays, but with a luminosity of only $\sim 10^{35} \rm\,erg\,s^{-1}$. For wind velocities $v_{\rm w} \sim 0.001c$, the shocked wind becomes radiative, but it is no longer hot enough to produce X-rays. Instead it emits in optical and UV, and the radiative luminosity is comparable to 100% of the wind kinetic luminosity. We suggest that measuring the optical luminosities and putting limits on the X-ray and radio emission from shock-ionized ULX bubbles may help in estimating the mass outflow rate of the central accretion disk and the velocity of the outflow. △ Less

Submitted 11 May, 2017; originally announced May 2017.

Comments: 11 pages, 7 figures