No Arabic abstract
Milli-second pulsars with highly stable periods can be considered as very precise clocks and can be used for pulsar timing array (PTA) which attempts to detect nanoheltz gravitational waves (GWs) directly. Main sources of nanoheltz GWs are supermassive black hole (SMBH) binaries which have sub-pc-scale orbits. On the other hand, a SMBH binary which is in an earlier phase and has pc-scale orbit emits ultra-low-frequency ($lesssim 10^{-9},mathrm{Hz}$) GWs cannot be detected with the conventional methodology of PTA. Such binaries tend to obtain high eccentricity, possibly $sim 0.9$. In this paper, we develop a formalism for extending constraints on GW amplitudes from single sources obtained by PTA toward ultra-low frequencies considering the waveform expected from an eccentric SMBH binary. GWs from an eccentric binaries are contributed from higher harmonics and, therefore, have a different waveform those from a circular binary. Furthermore, we apply our formalism to several hypothetical SMBH binaries at the center of nearby galaxies, including M87, using the constraints from NANOGravs 11-year data set. For a hypothetical SMBH binary at the center of M87, the typical upper limit on the mass ratio is $0.16$ for eccentricity of $0.9$ and semi-major axis of $a=1~mathrm{pc}$, assuming the binary phase to be the pericenter.
Gravitational waves (GWs) in the nano-hertz band are great tools for understanding the cosmological evolution of supermassive black holes (SMBHs) in galactic nuclei. We consider SMBH binaries in high-$z$ ultra-luminous infrared galaxies (ULIRGs) as sources of a stochastic GW background (GWB). ULIRGs are likely associated with gas-rich galaxy mergers containing SMBHs that possibly occur at most once in the life of galaxies, unlike multiple dry mergers at low redshift. Adopting a well-established sample of ULIRGs, we study the properties of the GWB due to coalescing binary SMBHs in these galaxies. Since the ULIRG population peaks at $z>1.5$, the amplitude of the GWB is not affected even if BH mergers are delayed by as long as $sim $ 10 Gyrs. Despite the rarity of the high-$z$ ULIRGs, we find a tension with the upper limits from Pulsar Timing Array (PTA) experiments. This result suggests that if a fraction $f_{rm m,gal}$ of ULIRGs are associated with SMBH binaries, then no more than $20 f_{rm m,gal}(lambda_{rm Edd}/0.3)^{5/3}(t_{rm life}/30~{rm Myr})~%$ of the binary SMBHs in ULIRGs can merge within a Hubble time, for plausible values of the Eddington ratio of ULIRGs ($lambda_{rm Edd}$) and their lifetime ($t_{rm life}$).
We probe ultra-low-frequency gravitational waves (GWs) with statistics of spin-down rates of milli-second pulsars (MSPs) by a method proposed in our prevous work (Yonemaru et al. 2016). The considered frequency range is $10^{-12}{rm Hz} lesssim f_{rm GW} lesssim 10^{-10}$Hz, which cannot be accessed by the conventional pulsar timing array. The effect of such low-frequency GWs appears as a bias to spin-down rates which has a quadrupole pattern in the sky. We use the skewness of the spin-down rate distribution and the number of MSPs with negative spin-down rates to search for the bias induced by GWs. Applying this method to 149 MSPs selected from the ATNF pulsar catalog, we derive upper bounds on the time derivative of the GW amplitudes of $dot{h} < 6.2 times 10^{-18}~{rm sec}^{-1}$ and $dot{h} < 8.1 times 10^{-18}~{rm sec}^{-1}$ in the directions of the Galactic Center and M87, respectively. Approximating the GW amplitude as $dot{h} sim 2 pi f_{rm GW} h$, the bounds translate into $h < 3 times 10^{-9}$ and $h < 4 times 10^{-9}$, respectively, for $f_{rm GW} = 1/(100~{rm yr})$. Finally, we give the implications to possible super-massive black hole binaries at these sites.
The advent of time domain astronomy is revolutionizing our understanding of the Universe. Programs such as the Catalina Real-time Transient Survey (CRTS) or the Palomar Transient Factory (PTF) surveyed millions of objects for several years, allowing variability studies on large statistical samples. The inspection of $approx$250k quasars in CRTS resulted in a catalogue of 111 potentially periodic sources, put forward as supermassive black hole binary (SMBHB) candidates. A similar investigation on PTF data yielded 33 candidates from a sample of $approx$35k quasars. Working under the SMBHB hypothesis, we compute the implied SMBHB merger rate and we use it to construct the expected gravitational wave background (GWB) at nano-Hz frequencies, probed by pulsar timing arrays (PTAs). After correcting for incompleteness and assuming virial mass estimates, we find that the GWB implied by the CRTS sample exceeds the current most stringent PTA upper limits by almost an order of magnitude. After further correcting for the implicit bias in virial mass measurements, the implied GWB drops significantly but is still in tension with the most stringent PTA upper limits. Similar results hold for the PTF sample. Bayesian model selection shows that the null hypothesis (whereby the candidates are false positives) is preferred over the binary hypothesis at about $2.3sigma$ and $3.6sigma$ for the CRTS and PTF samples respectively. Although not decisive, our analysis highlights the potential of PTAs as astrophysical probes of individual SMBHB candidates and indicates that the CRTS and PTF samples are likely contaminated by several false positives.
We present results from a controlled numerical experiment investigating the effect of stellar density gas on the coalescence of binary black holes (BBHs) and the resulting gravitational waves (GWs). This investigation is motivated by the proposed stellar core fragmentation scenario for BBH formation and the associated possibility of an electromagnetic counterpart to a BBH GW event. We employ full numerical relativity coupled with general-relativistic hydrodynamics and set up a $30 + 30 M_odot$ BBH (motivated by GW150914) inside gas with realistic stellar densities. Our results show that at densities $rho gtrsim 10^6 - 10^7 , mathrm{g , cm}^{-3}$ dynamical friction between the BHs and gas changes the coalescence dynamics and the GW signal in an unmistakable way. We show that for GW150914, LIGO observations conclusively rule out BBH coalescence inside stellar gas of $rho gtrsim 10^7 , mathrm{g,cm}^{-3}$. Typical densities in the collapsing cores of massive stars are in excess of this density. This excludes the fragmentation scenario for the formation of GW150914.
Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, $A_{rm c,yr}$, to be < $1.0times10^{-15}$ with 95% confidence. This limit excludes predicted ranges for $A_{rm c,yr}$ from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves.