Abstract The cosmic distance duality relation (CDDR), expressed as d L ( z ) = (1 + z ) 2 D A ( z ), is a fundamental relation in modern cosmology. In this work, we apply a method to test the CDDR using simulated strongly lensed gravitational-wave (SLGW) signals from massive binary black holes as observed by proposed space-based detector networks. Our analysis is conducted under the point-mass lens model, considering the strong lensing scenario that produces two images. We generate 90 days of simulated SLGW data for 10 events based on the Population III stellar formation model, with source redshifts in the range z s ∈ 2, 6 and lens redshifts in z L ∈ 0.2, 1. The deviation of CDDR is parameterized by η 1 ( z ) = 1 + η 0 z and η 2 ( z ) = 1 + η 0 z /(1 + z ), and we incorporate the deviation parameter η 0 directly into the waveform model. Parameter estimation is performed within a Bayesian statistical framework, combining simulated data from both Taiji and LISA. For a single lensed event, the joint Taiji+LISA analysis improves the measurement precision of η 0 by roughly a factor of 2 compared with Taiji-only observations. By combining 10 simulated events, the population-level constraints on η 0 , quantified by the half width of the 95% credible interval, reach approximately 2.61 × 10 −4 (1.72 × 10 −4 ) for the η 1 ( z ) parameterization and 1.22 × 10 −3 (6.86 × 10 −4 ) for η 2 ( z ) in the Taiji-only (Taiji+LISA) scenario, respectively. The inferred values of η 0 remain consistent with η 0 = 0 within the estimated uncertainties, with no statistically significant evidence for deviations from the CDDR at the achieved precision. These results demonstrate the significant advantage of joint space-based observations for high-precision tests of the CDDR.
Yuan et al. (Wed,) studied this question.