An Accurate Modeling of Nano-hertz Gravitational Wave Signal from Eccentric Supermassive Binary Black Holes: An Essential Step Toward a Robust Discovery

TL;DR

This work develops a robust framework to model the nanohertz SGWB from eccentric SMBHB populations by integrating a multi-scale galaxy– SMBH framework with an accurate time-domain eccentric waveform model ESIGMAHM. The approach produces SGWB maps and introduces the spectral covariance matrix $C_{\rm N}(f_1,f_2)$ to capture inter-frequency correlations induced by multi-harmonic emission from eccentric binaries, offering a novel observable beyond the isotropic spectrum $\Omega_{\rm GW}(f)$. Results show that eccentricity can leave measurable imprints on both the SGWB density and its spectral covariance, with stronger off-diagonal structure for higher eccentricity, and that next-generation PTAs (e.g., SKA-scale with ~2000 pulsars) could distinguish different eccentricity distributions, especially under low red-noise conditions. This framework thus provides a path to jointly infer SMBHB population properties and break degeneracies with environmental effects, guiding future PTA analyses toward more accurate SMBHB astrophysics.

Abstract

The stochastic gravitational wave background (SGWB) in the nanohertz (nHz) regime, detectable by pulsar timing arrays (PTAs), provides a promising probe of the cosmic population of supermassive black hole binaries (SMBHBs). These binaries are expected to retain significant eccentricity throughout their evolution. We present a new technique to model the nHz SGWB by incorporating eccentricity into a multi-scale adaptive simulation framework. Using the time-domain eccentric waveform model ESIGMAHM, we generate realistic GW signals from astrophysical populations of SMBHBs. Unlike circular binaries, eccentric systems emit across multiple frequencies, introducing spectral correlations between frequency bins. These correlations provide a novel observational signature of the eccentricity distribution of the SMBHB population. In this work, we adopt simplified power-law models for the eccentricity distribution. While this does not capture the full complexity of galactic environments, it effectively highlights the key features of GW emission from eccentric binaries and their imprint on the SGWB. Our approach advances nHz GW signal modeling by incorporating eccentricity at small scales, enabling more realistic predictions and offering a new avenue for probing SMBHB astrophysics with future PTA observations.

Figures (12)

• Figure 1: Schematic diagram summarizing our simulation technique, showing the multi-scale framework linking Gpc-scale structure (MICECAT), Mpc-scale SMBH–galaxy co-evolution (ROMULUS25), and sub-pc binary dynamics including environment and eccentricity, to model the nHz SGWB. The SMBHB population is generated via a Monte Carlo sampling procedure based on the MICECAT simulated galaxy catalog. For each binary, a time-domain gravitational waveform is computed using eccentric waveform modeling. These waveforms are then truncated to match the chosen observation time and are subsequently Fourier transformed. Finally, the SGWB map is constructed by summing the GW energy density contributions from all sources within each sky pixel.
• Figure 2: Time-domain waveforms generated by InspiralESIGMAHM along with their Fourier transforms for a circular binary ($e = 0$). $h^{\mathbb{R}}_{\ell m}$ and $h^{\mathbb{I}}_{\ell m}$ represent the real and imaginary parts of the $(\ell,m)$ GW mode, while $|\tilde{h}^{\mathbb{R}}_{\ell m}|$ and $|\tilde{h}^{\mathbb{I}}_{\ell m}|$ denote the absolute values of their respective Fourier transforms. The $(2,2)$ mode is shown along with the next three dominant modes: $(2,1)$, $(3,3)$, and $(4,4)$.
• Figure 3: Time-domain waveforms generated by InspiralESIGMAHM along with their Fourier transforms for an eccentric binary ($e$ = 0.5). $h^{\mathbb{R}}_{\ell m}$ and $h^{\mathbb{I}}_{\ell m}$ represent the real and imaginary parts of the $(\ell,m)$ GW mode, while $|\tilde{h}^{\mathbb{R}}_{\ell m}|$ and $|\tilde{h}^{\mathbb{I}}_{\ell m}|$ denote the absolute values of their respective Fourier transforms. The $(2,2)$ mode is shown along with the next three dominant modes: $(2,1)$, $(3,3)$, and $(4,4)$.
• Figure 4: Eccentricity distributions of SMBHBs modeled using a power-law probability density function, $P(e) \propto e^{\zeta}$. Distributions are shown for three values of the power-law index: $\zeta = -1$ (favoring nearly circular binaries), $\zeta = 0$ (uniform distribution), and $\zeta = 1$ (favoring highly eccentric binaries).
• Figure 5: SGWB spectrum, $\Omega_{\rm GW}(f)$, for three different eccentricity distributions characterized by $\zeta = -1$, $\zeta = 0$, and $\zeta = 1$. The solid lines represent the median $\Omega_{\rm GW}(f)$ obtained from 300 Monte Carlo realizations, while the shaded regions denote the $68\%$ credible intervals.