Cosmological parameters encoding our understanding of the expansion history of the Universe can be constrained by the accurate estimation of time delays arising in gravitationally lensed systems. We propose TD-CARMA, a Bayesian method to estimate cosmological time delays by modelling the observed and irregularly sampled light curves as realizations of a Continuous Auto-Regressive Moving Average (CARMA) process. Our model accounts for heteroskedastic measurement errors and microlensing, an additional source of independent extrinsic long-term variability in the source brightness. The semi-separable structure of the CARMA covariance matrix allows for fast and scalable likelihood computation using Gaussian Process modeling. We obtain a sample from the joint posterior distribution of the model parameters using a nested sampling approach. This allows for ``painless'' Bayesian Computation, dealing with the expected multi-modality of the posterior distribution in a straightforward manner and not requiring the specification of starting values or an initial guess for the time delay, unlike existing methods. In addition, the proposed sampling procedure automatically evaluates the Bayesian evidence, allowing us to perform principled Bayesian model selection. TD-CARMA is parsimonious, and typically includes no more than a dozen unknown parameters. We apply TD-CARMA to six doubly lensed quasars HS 2209+1914, SDSS J1001+5027, SDSS J1206+4332, SDSS J1515+1511, SDSS J1455+1447, SDSS J1349+1227, estimating their time delays as $-21.96 \pm 1.448$, $120.93 \pm 1.015$, $111.51 \pm 1.452$, $210.80 \pm 2.18$, $45.36 \pm 1.93$ and $432.05 \pm 1.950$ respectively. These estimates are consistent with those derived in the relevant literature, but are typically two to four times more precise.

翻译：宇宙膨胀历史的理解依赖于对宇宙学参数的约束，而精确估计引力透镜系统中的时间延迟是其中的关键。本文提出TD-CARMA贝叶斯方法，通过将观测到的非均匀采样光变曲线建模为连续自回归滑动平均（CARMA）过程的实现，来估计宇宙学时间延迟。该模型同时考虑了异方差测量误差和微引力透镜效应（一种源亮度中额外的独立外部长周期变异性）。利用CARMA协方差矩阵的半可分离结构，我们通过高斯过程建模实现了快速且可扩展的似然计算。采用嵌套采样方法获取模型参数联合后验分布的样本，实现了"无痛"贝叶斯计算：既能直接处理后验分布预期存在的多峰性，又无需像现有方法那样指定起始值或时间延迟的初始猜测。此外，该采样过程可自动计算贝叶斯证据，从而支持基于原则的贝叶斯模型选择。TD-CARMA模型简洁紧凑，通常仅包含不超过12个未知参数。我们将TD-CARMA应用于六个双透镜类星体HS 2209+1914、SDSS J1001+5027、SDSS J1206+4332、SDSS J1515+1511、SDSS J1455+1447和SDSS J1349+1227，估计其时间延迟分别为-21.96±1.448天、120.93±1.015天、111.51±1.452天、210.80±2.18天、45.36±1.93天和432.05±1.950天。这些估计值与相关文献结果一致，但精度通常提高了2到4倍。