Mohit Raj Sah, Akash Maurya, Suvodip Mukherjee, Prayush Kumar, Vida Saeedzadeh, Arif Babul, Chandra Kant Mishra, Kaushik Paul, Thomas R. Quinn, Michael Tremmel

The stochastic gravitational wave background (SGWB) in the nanohertz (nHz) regime, detectable by pulsar timing arrays (PTAs), offers a promising avenue to probe the cosmic population of supermassive black hole binaries (SMBHBs). These SMBHBs are expected to retain substantial eccentricity throughout their evolution due to their formation history. In this study, we propose a new modeling scenario of the nHz SGWB by incorporating the eccentricity of SMBHBs into a multi-scale adaptive simulation-based framework. We employ a time-domain eccentric waveform model, \esigmahm{}, to generate realistic gravitational wave (GW) signals from an astrophysical population of SMBHB, including physical effects from sub-parsec scales to Gpc scales. The eccentric inspiraling binary, unlike circular binaries, emits GW signal in multiple frequencies. As a consequence, the SGWB energy density in each frequency bin is not independent; instead, the presence of eccentricity introduces a spectral correlation between different frequencies. We show that these spectral correlations are absent for circular binaries but become increasingly significant for populations with higher eccentricities. Our novel approach can capture this effect and opens up the window towards measuring this with a high signal-to-noise ratio with future observations. This work develops the frontier of nHz signal modeling using eccentricity at small scales and can model realistic nHz signal, which will be essential for robust inference from future observations to shed light on the astrophysical properties of SMBHBs.