Binary Black Hole Mergers Research Articles

Dropping Anchor: Understanding the Populations of Binary Black Holes with Random and Aligned-spin Orientations

The relative spin orientations of black holes (BHs) in binaries encode their evolutionary history: BHs assembled dynamically should have isotropically distributed spins, while spins of BHs originating in the field should be aligned with the orbital angular momentum. In this article, we introduce a simple population model for these dynamical and field binaries that uses spin orientations as an anchor to disentangle these two evolutionary channels. We then analyze binary BH mergers in the Third Gravitational-Wave Transient Catalog (GWTC-3) and ask whether BHs from the isotropic-spin population possess different distributions of mass ratios, spin magnitudes, or redshifts from the preferentially aligned-spin population. We find no compelling evidence that binary BHs in GWTC-3 have different source-property distributions depending on their spin alignment, but we do find that the dynamical and field channels cannot both have mass-ratio distributions that strongly favor equal masses. We give an example of how this can be used to provide insights into the various processes that drive these BHs to merge. We also find that the current detections are insufficient in extracting differences in spin magnitude or redshift distributions of isotropic and aligned-spin populations.

Merging binary black holes formed through double-core evolution

Context. To date, various formation channels of merging events have been heavily explored with the detection of nearly 100 double black hole (BH) merger events reported by the LIGO-Virgo-KAGRA (LVK) Collaboration. In this paper, we systematically investigate an alternative formation scenario: binary BHs (BBHs) formed through double helium stars (hereafter, “double-core evolution channel”). In this scenario, two helium stars (He-rich stars) could be the outcome of the classical isolated binary evolution scenario with and without the common envelope (CE) phase (i.e., CE channel and stable mass transfer channel) or, alternatively, of massive close binaries evolving chemically homogeneously (i.e., CHE channel). Aims. We study the properties (i.e., the chirp masses and the effective spins) of BBHs formed through the double-core evolution and investigate the impact of different efficiencies of angular momentum transport within massive He-rich stars on double-core evolution. Methods. We performed detailed stellar structure and binary evolution calculations that take into account internal rotation and mass loss of He-rich stars as well as tidal interactions in binaries. We systematically studied the parameter space of initial binary He-rich stars, including the initial mass and metallicity of He-rich stars as well as initial orbital periods. Apart from direct core collapse with mass and angular momentum conserved, we also follow the framework in Batta & Ramirez-Ruiz (2019, ArXiv e-prints [arXiv:1904.04835]) to estimate the mass and spin of the resulting BHs. Results. We show that the radii of massive He-rich stars decrease as a function of time, which comes mainly from mass loss and mixing in high metallicity and from mixing in low metallicity. For double He-rich stars with equal masses in binaries, we find that tides start to be at work on the zero age helium main sequence (i.e., the time when a He-rich star starts to burn helium in the core, which is analogous to zero age main sequence for core hydrogen burning) for initial orbital periods not longer than 1.0 day, depending on the initial metallicities. In addition to the stellar mass-loss rate and tidal interactions in binaries, we find that the role of the angular momentum transport efficiency in determining the resulting BH spins becomes stronger when considering BH progenitors originated from a higher metal-metallicity environment. We highlight that the double-core evolution scenario does not always produce fast-spinning BBHs and compare the properties of the BBHs reported from the LVK with our modeling. Conclusions. After detailed binary calculations of double-core evolution, we have confirmed that the spin of the BH is not only determined by the interplay of the binary’s different initial conditions (metallicity, mass, and orbital period) but is also dependent on the angular momentum transport efficiency within its progenitor. We predict that with the sensitivity improvements to the LVK’s next observing run (O4), the sample of merging BBHs will contain more sources with positive but moderate (even high) χeff and part of the events will likely show to have been formed through the double-core evolution channel.

Dictionary Learning: A Novel Approach to Detecting Binary Black Holes in the Presence of Galactic Noise with LISA.

The noise produced by the inspiral of millions of white dwarf binaries in the MilkyWay may pose a threat to one of the main goals of the space-based LISA mission: the detection of massive black hole binary mergers. We present a novel study for reconstruction of merger waveforms in the presence of Galactic confusion noise using dictionary learning. We discuss the limitations of untangling signals from binaries with total mass from 10^{2} M_{⊙} to 10^{4} M_{⊙}. Our method proves extremely successful for binaries with total mass greater than ∼3×10^{3} M_{⊙} up to redshift 3 in conservative scenarios, and up to redshift 7.5 in optimistic scenarios. In addition, consistently good waveform reconstruction of merger events is found if the signal-to-noise ratio is approximately 5 or greater.

Biorthogonal harmonics for the decomposition of gravitational radiation. I. Angular modes, completeness, and the introduction of adjoint-spheroidal harmonics

The estimation of radiative modes is a central problem in gravitational wave theory, with essential applications in signal modeling and data analysis. This problem is complicated by most astrophysically relevant systems not having modes that are analytically tractable. A ubiquitous workaround is to use not modes, but multipole moments defined by spin-weighted spherical harmonics. However, spherical multipole moments are only related to the modes of systems without angular momentum. As a result, they can obscure the underlying physics of astrophysically relevant systems, such as binary black hole merger and ringdown. In such cases, spacetime angular momentum means that radiative modes are not spherical, but spheroidal in nature. Here, we work through various problems related to spheroidal harmonics. We show for the first time that spheroidal harmonics are not only capable of representing arbitrary gravitational wave signals, but that they also possess a kind of orthogonality not used before in general relativity theory. Along the way we present a new class of spin-weighted harmonic functions dubbed ``adjoint-spheroidal'' harmonics. These new functions may be used for the general estimation of spheroidal multipole moments via complete biorthogonal decomposition (in the angular domain). By construction, adjoint-spheroidal harmonics suppress mode-mixing effects known to plague spherical harmonic decomposition; as a result, they better approximate a system's true radiative modes. We discuss potential applications of these results. Lastly, we briefly comment on the challenges posed by the analogous problem with Teukolsky's radial functions.

Gravitational recoil from binary black hole mergers in scalar field clouds

In vacuum, the gravitational recoil of the final black hole from the merger of two black holes depends exclusively on the mass ratio and spins of the coalescing black holes, and on the eccentricity of the binary. If matter is present, accretion by the merging black holes may modify significantly their masses and spins, altering both the dynamics of the binary and the gravitational recoil of the remnant black hole. This paper considers such a scenario. We investigate the effects on the kick of the final black hole from immersing the binary in a scalar field cloud. We consider two types of configurations: one with nonspinning and unequal-mass black holes, and a second with equal mass and spinning holes. For both types, we investigate how the gravitational recoil of the final black hole changes as we vary the energy density of the scalar field. We find that the accretion of the scalar field by the merging black holes could have a profound effect. For the nonspinning, unequal-mass binary black holes, the kicks are in general larger than in the vacuum case, with speeds of $\ensuremath{\sim}1,200\text{ }\text{ }\mathrm{km}/\mathrm{s}$ for binaries with mass ratio $2\ensuremath{\mathbin:}1$, 1 order of magnitude larger than in vacuum. For equal mass, binaries with black holes with spins aligned with the orbital angular momentum, kicks larger than in vacuum are also found. For systems with spins in the superkick configuration, the scalar field triggers a similar dependence of the kicks with the entrance angle at merger as in the vacuum case but in this case depending on the strength of the scalar field.

Observational Imprints of Enhanced Scalar Power on Small Scales in Ultra Slow Roll Inflation and Associated Non-Gaussianities

The discovery of gravitational waves from merging binary black holes has generated considerable interest in examining whether these black holes could have a primordial origin. If a significant number of black holes have to be produced in the early universe, the primordial scalar power spectrum should have an enhanced amplitude on small scales, when compared to the COBE normalized values on the large scales that is strongly constrained by the anisotropies in the cosmic microwave background. In the inflationary scenario driven by a single, canonical scalar field, such power spectra can be achieved in models that permit a brief period of ultra slow roll inflation during which the first slow roll parameter decreases exponentially. In this review, we shall consider a handful of such inflationary models as well as a reconstructed scenario and examine the extent of formation of primordial black holes and the generation of secondary gravitational waves in these cases. We shall also discuss the strength and shape of the scalar bispectrum and the associated non-Gaussianity parameter that arise in such situations. We shall conclude with an outlook wherein we discuss the wider implications of the increased strengths of the non-Gaussianities on smaller scales.

Understanding the high-mass binary black hole population from stable mass transfer and super-Eddington accretion in bpass

ABSTRACT With the remarkable success of the Ligo-Virgo-Kagra consortium in detecting binary black hole mergers, it has become possible to use the population properties to constrain our understanding of the progenitor stars’ evolution. The most striking features of the observed primary black hole mass distributions are the extended tail up to 100 M⊙ and an excess of masses at 35 M⊙. Currently, isolated binary population synthesis have difficulty explaining these features. Using the well-tested bpass detailed stellar binary evolution models to determine mass transfer stability, accretion rates, and remnant masses, we postulate that stable mass transfer with super-Eddington accretion is responsible for the extended tail. These systems are able to merge within the Hubble time due to more stable mass transfer at higher donor masses with higher mass ratios and spin-orbit coupling, allowing the orbits to shrink sufficiently. Furthermore, we find that in bpass the 35 M⊙ excess is not due to pulsational pair-instability, as previously thought, but a feature caused by stable mass transfer, whose regime is limited by the mass transfer stability, quasi-homogeneous evolution, and stellar winds. These findings are at odds with those from other population synthesis codes but in agreement with other recent studies using detailed binary evolution models.

Gaussian processes for glitch-robust gravitational-wave astronomy

ABSTRACT Interferometric gravitational-wave observatories have opened a new era in astronomy. The rich data produced by an international network enable detailed analysis of the curved space-time around black holes. With nearly 100 signals observed so far and thousands expected in the next decade, their population properties enable insights into stellar evolution and the expansion of our Universe. However, the detectors are afflicted by transient noise artefacts known as ‘glitches’ which contaminate the signals and bias inferences. Of the 90 signals detected to date, 18 were contaminated by glitches. This feasibility study explores a new approach to transient gravitational-wave data analysis using Gaussian processes, which model the underlying physics of the glitch-generating mechanism rather than the explicit realization of the glitch itself. We demonstrate that if the Gaussian process kernel function can adequately model the glitch morphology, we can recover the parameters of simulated signals. Moreover, we find that the Gaussian processes kernels used in this work are well suited to modelling long-duration glitches which are most challenging for existing glitch-mitigation approaches. Finally, we show how the time-domain nature of our approach enables a new class of time-domain tests of General Relativity, performing a re-analysis of the inspiral-merger-ringdown test on the first observed binary black hole merger. Our investigation demonstrates the feasibility of the Gaussian processes as an alternative to the traditional framework but does not yet establish them as a replacement. Therefore, we conclude with an outlook on the steps needed to realize the full potential of the Gaussian process approach.

Prospects for LISA to detect a gravitational-wave background from first order phase transitions

First order phase transitions in the early universe could produce a gravitational-wave background that might be detectable by the Laser Interferometer Space Antenna (LISA). Such an observation would provide evidence for physics beyond the Standard Model. We study the ability of LISA to observe a gravitational-wave background from phase transitions in the presence of an extragalactic foreground from binary black hole mergers throughout the universe, a galactic foreground from white dwarf binaries, and LISA noise. Modelling the phase transition gravitational wave background as a double broken power law, we use the deviance information criterion as a detection statistic, and Fisher matrix and Markov Chain Monte Carlo methods to assess the measurement accuracy of the parameters of the power spectrum. While estimating all the parameters associated with the gravitational-wave backgrounds, foregrounds, and LISA noise, we find that LISA could detect a gravitational-wave background from phase transitions with a peak frequency of 1 mHz and normalized energy density amplitude of Ωp ≃ 3 × 10-11. With Ωp ≃ 10-10, the signal is detectable if the peak frequency is in the range 4 × 10-4 to 9 × 10-3 Hz, and the peak amplitude and frequency can be estimated to an accuracy of 10% to 1%.

Gravitational wave inference on a numerical-relativity simulation of a black hole merger beyond general relativity

We apply common gravitational wave inference procedures on binary black hole merger waveforms beyond general relativity. We consider dynamical Chern-Simons gravity, a modified theory of gravity with origins in string theory and loop quantum gravity. This theory introduces an additional parameter $\ell$, corresponding to the length-scale below which beyond-general-relativity effects become important. We simulate data based on numerical relativity waveforms produced under an approximation to this theory, which differ from those of general relativity in the strongly nonlinear merger regime. We consider a system with parameters similar to GW150914 with different values of $\ell$ and signal-to-noise ratios. We perform two analyses of the simulated data. The first is a template-based analysis that uses waveforms derived under general relativity and allows us to identify degeneracies between the two waveform morphologies. The second is a morphology-independent analysis based on BayesWave that does not assume that the signal is consistent with general relativity. The BayesWave analysis faithfully reconstructs the simulated signals. However, waveform models derived under general relativity are unable to fully mimic the simulated modified-gravity signals and such a deviation would be identifiable with existing inference tools. Depending on the magnitude of the deviation, we find that the templated analysis can under perform the morphology-independent analysis in fully recovering simulated beyond-GR waveforms even for achievable signal-to-noise ratios $\gtrsim 20{-}30$.

Isolated and dynamical black hole mergers withB-POP: the role of star formation and dynamics, star cluster evolution, natal kicks, mass and spins, and hierarchical mergers

ABSTRACTThe current interpretation of LIGO–Virgo–KAGRA data suggests that the primary mass function of merging binary black holes (BBHs) at redshift z ≲ 1 contains multiple structures, while spins are relatively low. Theoretical models of BBH formation in different environments can provide a key to interpreting the population of observed mergers, but they require the simultaneous treatment of stellar evolution and dynamics, galaxy evolution, and general relativity. We present B-POP, a population synthesis tool to model BBH mergers formed in the field or via dynamical interactions in young, globular, and nuclear clusters. Using B-POP, we explore how black hole (BH) formation channels, star cluster evolution, hierarchical mergers, and natal BH properties affect the population of BBH mergers. We find that the primary mass distribution of BBH mergers extends beyond $M_1 \simeq 200\, {}$ M⊙, and the effective spin parameter distribution hints at different natal spins for single and binary BHs. Observed BBHs can be interpreted as members of a mixed population comprised of $\sim 34{{\ \rm per\ cent}} \, {}(66{{\ \rm per\ cent}})$ isolated (dynamical) BBHs, with the latter likely dominating at redshift z > 1. Hierarchical mergers constitute the $4.6\!-\!7.9{{\ \rm per\ cent}}$ of all mergers in the reference model, dominating the primary mass distribution beyond $M_1 \gt 65\, {}$ M⊙. The inclusion of cluster mass-loss and expansion causes an abrupt decrease in the probability for mergers beyond the third generation to occur. Considering observational biases, we find that $2.7\!-\!7.5{{\ \rm per\ cent}}$ of mock mergers involve intermediate-mass black hole (IMBH) seeds formed via stellar collisions. Comparing this percentage to observed values will possibly help us to constrain IMBH formation mechanisms.

Measuring properties of primordial black hole mergers at cosmological distances: Effect of higher order modes in gravitational waves

Primordial black holes (PBHs) may form from the collapse of matter overdensities shortly after the Big Bang. One may identify their existence by observing gravitational wave (GW) emissions from merging PBH binaries at high redshifts $z\gtrsim 30$, where astrophysical binary black holes (BBHs) are unlikely to merge. The next-generation ground-based GW detectors, Cosmic Explorer and Einstein Telescope, will be able to observe BBHs with total masses of $\mathcal{O}(10-100)~M_{\odot}$ at such redshifts. This paper serves as a companion paper of arXiv:2108.07276, focusing on the effect of higher-order modes (HoMs) in the waveform modeling, which may be detectable for these high redshift BBHs, on the estimation of source parameters. We perform Bayesian parameter estimation to obtain the measurement uncertainties with and without HoM modeling in the waveform for sources with different total masses, mass ratios, orbital inclinations and redshifts observed by a network of next-generation GW detectors. We show that including HoMs in the waveform model reduces the uncertainties of redshifts and masses by up to a factor of two, depending on the exact source parameters. We then discuss the implications for identifying PBHs with the improved single-event measurements, and expand the investigation of the model dependence of the relative abundance between the BBH mergers originating from the first stars and the primordial BBH mergers as shown in arXiv:2108.07276.

First machine learning gravitational-wave search mock data challenge

We present the results of the first Machine Learning Gravitational-Wave Search Mock Data Challenge. For this challenge, participating groups had to identify gravitational-wave signals from binary black hole mergers of increasing complexity and duration embedded in progressively more realistic noise. The final of the 4 provided datasets contained real noise from the O3a observing run and signals up to a duration of 20 s with the inclusion of precession effects and higher order modes. We present the average sensitivity distance and run-time for the 6 entered algorithms derived from 1 month of test data unknown to the participants prior to submission. Of these, 4 are machine learning algorithms. We find that the best machine learning based algorithms are able to achieve up to 95% of the sensitive distance of matched-filtering based production analyses for simulated Gaussian noise at a false-alarm rate (FAR) of one per month. In contrast, for real noise, the leading machine learning search achieved 70%. For higher FARs the differences in sensitive distance shrink to the point where select machine learning submissions outperform traditional search algorithms at FARs $\ensuremath{\ge}200$ per month on some datasets. Our results show that current machine learning search algorithms may already be sensitive enough in limited parameter regions to be useful for some production settings. To improve the state-of-the-art, machine learning algorithms need to reduce the false-alarm rates at which they are capable of detecting signals and extend their validity to regions of parameter space where modeled searches are computationally expensive to run. Based on our findings we compile a list of research areas that we believe are the most important to elevate machine learning searches to an invaluable tool in gravitational-wave signal detection.

Black hole merger simulations in wave dark matter environments

The interaction of binary black hole mergers with their environments can be studied using numerical relativity simulations. These start only a short finite time before merger, at which point appropriate initial conditions must be imposed. A key task is therefore to identify the configuration that is appropriate for the binary and its environment at this stage of the evolution. In this work we study the behavior of wave dark matter around equal mass black hole binaries, finding that there is a preferred, quasistationary profile that persists and grows over multiple orbits, in contrast to heavier mass dark matter where any overdensity tends to be dispersed by the binary motion. While different initial configurations converge to the preferred quasistationary one after several orbits, unwanted transient oscillations are generated in the process, which may have an impact on the signal in short simulation runs. We also point out that naively superimposing the matter onto a circular binary results in artificially eccentric orbits due to the matter backreaction, which is an effect of the initial conditions and not a signature of dark matter. We discuss the further work required so that comparison of waveforms obtained with environments to vacuum cases can be done in a meaningful way.

Eccentric pairs: Analytic gravitational waves from binary black holes in elliptic orbits

Gravitational waves (GW) from eccentric binaries have intricate signals encoding important features about the location, creation and evolution of the sources. Eccentricity shortens the merger time, making the emitted GW statistically predominant in the observed data once detectors will reach the required sensitivity. We present a novel implementation of fully analytical GW templates from eccentric binary black hole (BBH) mergers within the Wolfram Mathematica software. We increase the accuracy by identifying and minimizing the possible source of errors. We start with an overview of the physics involved in eccentric mergers, then assemble the strain for the inspiral by employing up to six post-Newtonian (PN) corrections. We complete the eccentric inspiral with the quasi-circular Backwards one Body (BOB) merger model in frequency, amplitude and phase, then we build the hybrid GW strain for the whole evolution of the binary. For low eccentricity, we reach coincidence in the overlap, with no ambiguity in the time interval, a remarkable improvement from the usual matching techniques. For high-eccentricity, we compensate for the implicit quasi-circular assumption of the BOB approach, by introducing a small rescaling in amplitude. Our streamlined implementation is relevant for the new field of GW astronomy and is straightforward to understand, use and extend, offering researchers in the field a valuable open resource tool.

Systematic bias on the inspiral-merger-ringdown consistency test due to neglect of orbital eccentricity

The inspiral-merger-ringdown (IMR) consistency test checks the consistency of the final mass and final spin of a binary black hole merger remnant, independently inferred via the inspiral and merger-ringdown parts of the waveform. As binaries are expected to be nearly circularized when entering the frequency band of ground-based detectors, tests of general relativity (GR) currently employ quasicircular waveforms. We quantify the effect of residual orbital eccentricity on the IMR consistency test. We find that eccentricity causes a significant systematic bias in the inferred final mass and spin of the remnant black hole at an orbital eccentricity (defined at $10$ Hz) of $e_0 \gtrsim 0.1$ in the LIGO band (for a total binary mass in the range $65$-$200 \,M_{\odot}$). For binary black holes observed by Cosmic Explorer (CE), the systematic bias becomes significant for $e_0 \gtrsim 0.015$ (for $200$-$600 \,M_{\odot}$ systems). This eccentricity-induced bias on the final mass and spin leads to an apparent inconsistency in the IMR consistency test, manifesting as a false violation of GR. Hence, eccentric corrections to waveform models are important for constructing a robust test of GR, especially for third-generation detectors. We also estimate the eccentric corrections to the relationship between the inspiral parameters and the final mass and final spin; they are shown to be quite small.

Novel Ringdown Amplitude-Phase Consistency Test.

The ringdown signal emitted during a binary black hole coalescence can be modeled as a linear superposition of the characteristic damped modes of the remnant black hole that get excited during the merger phase. While checking the consistency of the measured frequencies and damping times against theKerr BH spectrum predicted by general relativity (GR) is a cornerstone of strong-field tests of gravity, the consistency of measured excitation amplitudes and phases have been largely left unexplored. For a nonprecessing, quasicircular binary black hole merger, we find that GR predicts a narrow region in the space of mode amplitude ratio and phase difference, independently of the spin of the binary components. Using this unexpected result, we develop a new null test of strong-field gravity which demands that the measured amplitudes and phases of different ringdown modes should lie within this narrow region predicted by GR. We call this the amplitude-phase consistency test and introduce a procedure for performing it using information from the ringdown signal. Lastly, we apply this test to the GW190521 event, using the multimodal ringdown parameters inferred by Capano etal. [arXiv:2105.05238]. While ringdown measurements errors for this event are large, we show that GW190521 is consistent with the amplitude-phase consistency test. Our test is particularly well suited for accommodating multiple loud ringdown detections as those expected in the near future, and can be used complementarily to standard black-hole spectroscopy as a proxy for modified gravity, compact objects other than black holes, binary precession and eccentricity.

A Light in the Dark: Searching for Electromagnetic Counterparts to Black Hole–Black Hole Mergers in LIGO/Virgo O3 with the Zwicky Transient Facility

The accretion disks of active galactic nuclei (AGNs) are promising locations for the merger of compact objects detected by gravitational wave (GW) observatories. Embedded within a baryon-rich, high-density environment, mergers within AGNs are the only GW channel where an electromagnetic (EM) counterpart must occur (whether detectable or not). Considering AGNs with unusual flaring activity observed by the Zwicky Transient Facility (ZTF), we describe a search for candidate EM counterparts to binary black hole (BBH) mergers detected by LIGO/Virgo in O3. After removing probable false positives, we find nine candidate counterparts to BBH mergers during O3 (seven in O3a, two in O3b) with a p-value of 0.0019. Based on ZTF sky coverage, AGN geometry, and merger geometry, we expect ≈3(N BBH/83)(f AGN/0.5) potentially detectable EM counterparts from O3, where N BBH is the total number of observed BBH mergers and f AGN is the fraction originating in AGNs. Further modeling of breakout and flaring phenomena in AGN disks is required to reduce our false-positive rate. Two of the events are also associated with mergers with total masses >100 M ⊙, which is the expected rate for O3 if hierarchical (large-mass) mergers occur in the AGN channel. Candidate EM counterparts in future GW observing runs can be better constrained by coverage of the Southern sky as well as spectral monitoring of unusual AGN flaring events in LIGO/Virgo alert volumes. A future set of reliable AGN EM counterparts to BBH mergers will yield an independent means of measuring cosmic expansion (H 0) as a function of redshift.

Multi-messenger studies with the Pierre Auger Observatory

The combination of data from observatories measuring ultra-high energy cosmic rays, photons, neutrinos and gravitational waves has provided new insights into the most extreme phenomena in the Universe. Sharing information within a broad community is the foundation of the multi-messenger approach. The Pierre Auger Observatory, the world’s largest cosmic ray detector, provides sensitivity to photons and neutrinos above 1017eV, thus contributing significantly to this joint effort. The latest results from diffuse and targeted searches will be reviewed here, along with results from follow-up analyses and future perspectives. In particular, preliminary limits on photon fluence from a selection of gravitational wave sources detected by LIGO/Virgo and results of the search for ultra-high energy neutrinos from binary black hole mergers will be presented.

Splitting of the atomic hydrogen 21-cm line induced by gravitational waves from a supermassive black hole binary

We investigate the characteristic of the hyperfine splitting of a ground-state hydrogen atom subjected to a field of a gravitational wave. In our consideration, hydrogen atoms are in the source’s local wave zone. The interaction between the atoms and the gravitational wave induces the shifts of some energy levels in the hyperfine structure, computed by the first-order perturbation theory. The splitting of energy levels leads to more than one hyperfine transition. Our analysis suggests that a gravitational wave emitted by supermassive black hole binaries could be strong enough to induce a significant effect on the hyperfine splitting. This radio-wave spectrum provides the gravitational-wave signature and thus offers a way of identifying the merger of supermassive black hole binaries.