skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Friday, May 16 until 2:00 AM ET on Saturday, May 17 due to maintenance. We apologize for the inconvenience.


Title: Host galaxy demographics of individually detectable supermassive black-hole binaries with pulsar timing arrays
Abstract Massive black hole binaries (MBHBs) produce gravitational waves (GWs) that are detectable with pulsar timing arrays. We determine the properties of the host galaxies of simulated MBHBs at the time they are producing detectable GW signals. The population of MBHB systems we evaluate is from theIllustriscosmological simulations taken in tandem with post processing semi-analytic models of environmental factors in the evolution of binaries. Upon evolving to the GW frequency regime accessible by pulsar timing arrays, we calculate the detection probability of each system using a variety of different values for pulsar noise characteristics in a plausible near-future International Pulsar Timing Array dataset. We find that detectable systems have host galaxies that are clearly distinct from the overall binary population and from most galaxies in general. With conservative noise factors, we find that host stellar metallicity, for example, peaks at 2 Z as opposed to the total population of galaxies which peaks at 0.6 Z . Additionally, the most detectable systems are much brighter in magnitude and more red in color than the overall population, indicating their likely identity as large ellipticals with diminished star formation. These results can be used to develop effective search strategies for identifying host galaxies and electromagnetic counterparts following GW detection by pulsar timing arrays.  more » « less
Award ID(s):
2307719 2007993 2146016 2125764 2020265
PAR ID:
10563209
Author(s) / Creator(s):
; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Classical and Quantum Gravity
Volume:
42
Issue:
2
ISSN:
0264-9381
Format(s):
Medium: X Size: Article No. 025021
Size(s):
Article No. 025021
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We use ALMA observations of CO(2–1) in 13 massive (M*≳ 1011M) poststarburst galaxies atz∼ 0.6 to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The poststarburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (Data Release 14) based on their spectral shapes, as part of the Studying QUenching at Intermediate-z Galaxies: Gas, angu L ar momentum, and Evolution ( SQuIGG L E ) program. Early results showed that two poststarburst galaxies host large H2reservoirs despite their low inferred star formation rates (SFRs). Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, with M H 2 10 9 M. Given their high stellar masses, this mass limit corresponds to an average gas fraction of f H 2 M H 2 / M * 7 % or ∼14% using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with theDn4000 spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched ≲150 Myr prior to observation host detectable CO(2–1) emission, while older poststarburst galaxies are undetected. The large H2reservoirs and low SFRs in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100–200 Myr, the SQuIGG L E galaxies require the additional and efficient heating or removal of cold gas to bring their low SFRs in line with standard H2scaling relations. 
    more » « less
  2. 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 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 the achromatic noise properties in each pulsar. The models we compare include existing DM models used by the North American Nanohertz Observatory for Gravitational waves (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 data set. 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 RN = 14.1 0.1 + 0.1 to 14.7 0.5 + 0.3 , and the spectral index broadens from γ RN = 2.6 0.4 + 0.5 to γ RN = 3.5 0.9 + 1.2 . 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. 
    more » « less
  3. Abstract We investigate the stellar mass–black hole mass ( * BH ) relation with type 1 active galactic nuclei (AGNs) down to BH = 10 7 M , corresponding to a ≃ −21 absolute magnitude in rest-frame ultraviolet, atz= 2–2.5. Exploiting the deep and large-area spectroscopic survey of the Hobby–Eberly Telescope Dark Energy Experiment (HETDEX), we identify 66 type 1 AGNs with BH ranging from 107–1010Mthat are measured with single-epoch virial method using Civemission lines detected in the HETDEX spectra. * of the host galaxies are estimated from optical to near-infrared photometric data taken with Spitzer, the Wide-field Infrared Survey Explorer, and ground-based 4–8 m class telescopes byCIGALEspectral energy distribution (SED) fitting. We further assess the validity of SED fitting in two cases by host-nuclear decomposition performed through surface brightness profile fitting on spatially resolved host galaxies with the James Webb Space Telescope/NIRCam CEERS data. We obtain the * BH relation covering the unexplored low-mass ranges of BH 10 7 10 8 M , and conduct forward modeling to fully account for the selection biases and observational uncertainties. The intrinsic * BH relation atz∼ 2 has a moderate positive offset of 0.52 ± 0.14 dex from the local relation, suggestive of more efficient black hole growth at higher redshift even in the low-mass regime of BH 10 7 10 8 M . Our * BH relation is inconsistent with the BH suppression at the low- * regime predicted by recent hydrodynamic simulations at a 98% confidence level, suggesting that feedback in the low-mass systems may be weaker than those produced in hydrodynamic simulations. 
    more » « less
  4. Abstract We present the stellar population properties of 69 short gamma-ray burst (GRB) host galaxies, representing the largest uniformly modeled sample to date. Using theProspectorstellar population inference code, we jointly fit photometry and/or spectroscopy of each host galaxy. We find a population median redshift of z = 0.64 0.32 + 0.83 (68% confidence), including nine photometric redshifts atz≳ 1. We further find a median mass-weighted age oftm= 0.8 0.53 + 2.71 Gyr, stellar mass of log(M*/M) = 9.69 0.65 + 0.75 , star formation rate of SFR = 1.44 1.35 + 9.37 Myr−1, stellar metallicity of log(Z*/Z) = 0.38 0.42 + 0.44 , and dust attenuation of A V = 0.43 0.36 + 0.85 mag (68% confidence). Overall, the majority of short GRB hosts are star-forming (≈84%), with small fractions that are either transitioning (≈6%) or quiescent (≈10%); however, we observe a much larger fraction (≈40%) of quiescent and transitioning hosts atz≲ 0.25, commensurate with galaxy evolution. We find that short GRB hosts populate the star-forming main sequence of normal field galaxies, but do not include as many high-mass galaxies as the general galaxy population, implying that their binary neutron star (BNS) merger progenitors are dependent on a combination of host star formation and stellar mass. The distribution of ages and redshifts implies a broad delay-time distribution, with a fast-merging channel atz> 1 and a decreased neutron star binary formation efficiency from high to low redshifts. If short GRB hosts are representative of BNS merger hosts within the horizon of current gravitational wave detectors, these results can inform future searches for electromagnetic counterparts. All of the data and modeling products are available on the Broadband Repository for Investigating Gamma-ray burst Host Traits website. 
    more » « less
  5. Abstract Pulsar timing array (PTA) searches for gravitational waves (GWs) aim to detect a characteristic correlation pattern in the timing residuals of galactic millisecond pulsars. This pattern is described by the PTA overlap reduction function (ORF) Γ a b ( ξ a b ) , which is known as the Hellings–Downs (HD) curve in general relativity (GR). In theories of modified gravity, the HD curve often receives corrections. Assuming, e.g. a subluminal GW phase velocity, one finds a drastically enhanced ORF in the limit of small angular separations between pulsaraand pulsarbin the sky, ξ a b 0 . In particular, working in harmonic space and performing an approximate resummation of all multipole contributions, the auto correlation coefficientΓaaseems to diverge. In this paper, we confirm that this divergence is unphysical and provide an exact and analytical expression forΓaain dependence of the pulsar distanceLaand the GW phase velocity v ph . In the GR limit and assuming a large pulsar distance, our expression reduces to Γ a a = 1 . In the case of subluminal phase velocity, we show that the regularization of the naive divergent result is a finite-distance effect, meaning thatΓaascales linearly withfLa, wherefis the GW frequency. For superluminal phase velocity (subluminal group velocity), which is relevant in the case of massive gravity, we correct an earlier analytical result forΓab. Our results pave the way for fitting modified-gravity theories with nonstandard phase velocity to PTA data, which requires a proper understanding of the auto correlation coefficientΓaa
    more » « less