Search for: All records

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 $\sim 2Z_{\odot }$as opposed to the total population of galaxies which peaks at $\sim 0.6Z_{\odot }$. 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. 

Abstract A population of compact object binaries emitting gravitational waves that are not individually resolvable will form a stochastic gravitational-wave signal. While the expected spectrum over population realizations is well known from Phinney, its higher-order moments have not been fully studied before or computed in the case of arbitrary binary evolution. We calculate analytic scaling relationships as a function of gravitational-wave frequency for the statistical variance, skewness, and kurtosis of a stochastic gravitational-wave signal over population realizations due to finite source effects. If the time derivative of the binary orbital frequency can be expressed as a power law in frequency, we find that these moment quantities also take the form of power-law relationships. We also develop a numerical population synthesis framework against which we compare our analytic results, finding excellent agreement. These new scaling relationships provide physical context to understanding spectral fluctuations in a gravitational-wave background signal and may provide additional information that can aid in explaining the origin of the nanohertz-frequency signal observed by pulsar timing array campaigns. 

ABSTRACT Time-domain data sets of many varieties can be prone to statistical outliers that result from instrumental or astrophysical anomalies. These can impair searches for signals within the time series and lead to biased parameter estimation. Versatile outlier mitigation methods tuned toward multimessenger time-domain searches for supermassive binary black holes have yet to be fully explored. In an effort to perform robust outlier isolation with low computational costs, we propose a Gibbs sampling scheme. This provides structural simplicity to outlier modelling and isolation, as it requires minimal modifications to adapt to time-domain modelling scenarios with pulsar-timing array or photometric data. We robustly diagnose outliers present in simulated pulsar-timing data sets, and then further apply our methods to pulsar J1909−3744 from the NANOGrav 9-year Data set. We also explore the periodic binary-AGN candidate PG1302−102 using data sets from the Catalina Real-time Transient Survey, All-Sky Automated Survey for Supernovae, and the Lincoln Near-Earth Asteroid Research. We present our findings and outline future work that could improve outlier modelling and isolation for multimessenger time-domain searches. 

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 data set. 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. By modeling the angular power distribution as a sum over spherical harmonics (where the coefficients are not bound to always generate positive power everywhere), we find that the Bayesian 95% upper limit on the level of dipole anisotropy is (C_l=1/C_l=0) < 27%. This is similar to the upper limit derived under the constraint of positive power everywhere, indicating that the dipole may be close to the data-informed regime. By contrast, the constraints on anisotropy at higher spherical-harmonic multipoles are strongly prior dominated. 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 data set and show that this data set 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. 

ABSTRACT The search for gravitational waves using Pulsar Timing Arrays (PTAs) is a computationally expensive complex analysis that involves source-specific noise studies. As more pulsars are added to the arrays, this stage of PTA analysis will become increasingly challenging. Therefore, optimizing the number of included pulsars is crucial to reduce the computational burden of data analysis. Here, we present a suite of methods to rank pulsars for use within the scope of PTA analysis. First, we use the maximization of the signal-to-noise ratio as a proxy to select pulsars. With this method, we target the detection of stochastic and continuous gravitational wave signals. Next, we present a ranking that minimizes the coupling between spatial correlation signatures, namely monopolar, dipolar, and Hellings & Downs correlations. Finally, we also explore how to combine these two methods. We test these approaches against mock data using frequentist and Bayesian hypothesis testing. For equal-noise pulsars, we find that an optimal selection leads to an increase in the log-Bayes factor two times steeper than a random selection for the hypothesis test of a gravitational wave background versus a common uncorrelated red noise process. For the same test but for a realistic European PTA (EPTA) data set, a subset of 25 pulsars selected out of 40 can provide a log-likelihood ratio that is 89 % of the total, implying that an optimally selected subset of pulsars can yield results comparable to those obtained from the whole array. We expect these selection methods to play a crucial role in future PTA data combinations. 

ABSTRACT We present a study of optically selected dual Active Galactic Nuclei (AGN) with projected separations of 3–97 kpc. Using multiwavelength (MWL) information (optical, X-ray, mid-IR), we characterized the intrinsic nuclear properties of this sample and compared them with those of isolated systems. Among the 124 X-ray-detected AGN candidates, 52 appear in pairs and 72 as single X-ray sources. Through MWL analysis, we confirmed the presence of the AGN in >80 per cent of the detected targets in pairs (42 out of 52). X-ray spectral analysis confirms the trend of increasing AGN luminosity with decreasing separation, suggesting that mergers may have contributed to triggering more luminous AGN. Through X-ray/mid-IR ratio versus X-ray colours, we estimated a fraction of Compton-thin AGN (with 1022 cm−2 < NH < 1024 cm−2) of about 80 per cent, while about 16 per cent are Compton-thick sources (with NH > 1024 cm−2). These fractions of obscured sources are larger than those found in samples of isolated AGN, confirming that pairs of AGN show higher obscuration. This trend is further confirmed by comparing the de-reddened [O iii] emission with the observed X-ray luminosity. However, the derived fraction of Compton-thick sources in this sample at the early stages of merging is lower than that reported for late-merging dual-AGN samples. Comparing NH from X-rays with that derived from E(B − V) from narrow-line regions, we found that the absorbing material is likely to be associated with the torus or broad-line regions. We also explored the X-ray detection efficiency of dual-AGN candidates, finding that, when observed properly (at on-axis positions and with long exposures), X-ray data represent a powerful way to confirm and investigate dual-AGN systems. 

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 NANOGrav’s recent 12.5 yr 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 nHz, we placed a sky-averaged limit ofh₀< (6.82 ± 0.35) × 10⁻¹⁵, andh₀< (2.66 ± 0.15) × 10⁻¹⁵in our most sensitive sky location. Finally, we placed a multimessenger limit of $\mathit{}<(1.41\pm 0.02)\times {10}^9{M}_{\odot }$on the chirp mass of the supermassive black hole binary candidate 3C 66B. 

ABSTRACT Supermassive black hole binaries (SMBHBs) are a natural outcome of galaxy mergers and should form frequently in galactic nuclei. Sub-parsec binaries can be identified from their bright electromagnetic emission, e.g. Active Galactic Nuclei (AGNs) with Doppler shifted broad emission lines or AGN with periodic variability, as well as from the emission of strong gravitational radiation. The most massive binaries (with total mass >108M⊙) emit in the nanohertz band and are targeted by Pulsar Timing Arrays (PTAs). Here we examine the synergy between electromagnetic and gravitational wave signatures of SMBHBs. We connect both signals to the orbital dynamics of the binary and examine the common link between them, laying the foundation for joint multimessenger observations. We find that periodic variability arising from relativistic Doppler boost is the most promising electromagnetic signature to connect with GWs. We delineate the parameter space (binary total mass/chirp mass versus binary period/GW frequency) for which joint observations are feasible. Currently multimessenger detections are possible only for the most massive and nearby galaxies, limited by the sensitivity of PTAs. However, we demonstrate that as PTAs collect more data in the upcoming years, the overlapping parameter space is expected to expand significantly. 

ABSTRACT The International Pulsar Timing Array 2nd data release is the combination of data sets 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 per cent upper limits on their amplitude h95. The most sensitive frequency is 10 nHz with h95 = 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.