Parameter estimation of eccentric massive black hole binaries with LISA and its cosmological implications

Future space-based gravitational wave (GW) observatories such as LISA will detect massive black hole binaries (MBHBs), which are expected to be accompanied by electromagnetic counterparts, thereby providing bright standard sirens for cosmology. The orbital eccentricity of MBHBs can significantly improve the parameter estimation of GWs because the multiple harmonics induced by eccentricity provide additional information and help break down the degeneracies among waveform parameters. In this paper, we use the EccentricFD waveform and construct 5-year GW event catalogs for LISA under three population models (popIII, Q3d and Q3nod). For the three models, we find that an initial eccentricity of $e_0=0.4$ at $10^{-4}$ Hz yields improvements in sky localization and distance inference by a factor of $\mathcal{O}(10)$ in the best cases. As a consequence, the average number of bright sirens increases substantially: from 8 to 11 (PopIII), 6 to 12 (Q3d) and 13 to 24 (Q3nod). This increase in event number, together with enhanced localization and distance inference, leads to tighter cosmological constraints. In the $\Lambda$CDM model, for instance, the relative uncertainty on $H_0$ is reduced from $8.17\%$ to $4.35\%$ for the Q3d model, corresponding to an improvement of approximately $47\%$. We also investigate the improvement in constraints on the dark energy equation of state and modified GW propagation when combining bright sirens with the latest cosmic microwave background data. These results demonstrate that eccentricity is a remarkably significant feature in GW detection and parameter estimation, enabling more accurate measurements of the Universe with future space-based observatories.