Prospects for constraining interacting dark energy cosmology with gravitational-wave bright sirens detected by future SKA-era pulsar timing arrays (2210.04000v4)

Abstract: We explore the constraints on cosmological parameters in interacting dark energy (IDE) models described by energy transfer rates $Q = \beta H \rho_{\rm de}$ and $Q = \beta H \rho_{\rm c}$, using simulated gravitational-wave (GW) bright siren data from pulsar timing arrays (PTAs) and the Planck 2018 cosmic microwave background (CMB) data. In particular, we simulate a future PTA observation in the FAST/SKA era with 20 millisecond pulsars (MSPs), each having 20\,ns white noise over a 10-year observation span, and demonstrate that this mock dataset significantly improves the constraint precision of key cosmological parameters such as the Hubble constant $H_0$, matter density $\Omega_m$, and the coupling parameter $\beta$. For the IDE model $Q = \beta H \rho_{\rm de}$, PTA data alone provides tighter constraints on these parameters than the CMB data alone, primarily due to the high sensitivity of GW standard sirens in probing the late universe. Combining PTA and CMB data further enhances the constraints by 43.6\% for $H_0$, 43.2\% for $\Omega_m$, and 44.7\% for $\beta$, relative to using CMB data alone. In contrast, for $Q = \beta H \rho_{\rm c}$, the CMB data alone constrains $\beta$ more tightly than the PTA data, due to the stronger impact of this interaction in the early universe. Nevertheless, the PTA+\,CMB combination still yields improvements of 13.3\% for $H_0$, 22.7\% for $\Omega_m$, and 18.2\% for $\beta$. Increasing the number of MSPs in the PTA further tightens all parameter constraints in both IDE models. Our results highlight the great potential of future PTA observations for significantly improving cosmological parameter estimation in IDE models, offering critical insights into the nature of dark energy and its interaction with dark matter.