Eccentricity effects on the SMBH GW background

We study how eccentricity affects the gravitational wave (GW) spectrum from supermassive black hole (SMBH) binaries. We develop a fast and accurate semi-analytic method for computing the GW spectra, the distribution for the spectral fluctuations and the correlations between different frequencies. As GW emission circularizes binaries, the suppression of the signal strength due to eccentricity will be relevant for signals from wider binaries emitting at lower frequencies. Such a feature is present in the signal observed at pulsar timing arrays. We find that, when orbital decay of the SMBH binaries is driven by GWs only, the shape of the observed signal prefers highly eccentric binaries e _2\, nHz = 0.83+0.04-0.05. However, when environmental effects are included, the initial eccentricity can be significantly lowered yet the scenario with purely circular binaries is mildly disfavored.