Open Access
September 2017 Bayesian estimates of astronomical time delays between gravitationally lensed stochastic light curves
Hyungsuk Tak, Kaisey Mandel, David A. van Dyk, Vinay L. Kashyap, Xiao-Li Meng, Aneta Siemiginowska
Ann. Appl. Stat. 11(3): 1309-1348 (September 2017). DOI: 10.1214/17-AOAS1027

Abstract

The gravitational field of a galaxy can act as a lens and deflect the light emitted by a more distant object such as a quasar. Strong gravitational lensing causes multiple images of the same quasar to appear in the sky. Since the light in each gravitationally lensed image traverses a different path length from the quasar to the Earth, fluctuations in the source brightness are observed in the several images at different times. The time delay between these fluctuations can be used to constrain cosmological parameters and can be inferred from the time series of brightness data or light curves of each image. To estimate the time delay, we construct a model based on a state-space representation for irregularly observed time series generated by a latent continuous-time Ornstein–Uhlenbeck process. We account for microlensing, an additional source of independent long-term extrinsic variability, via a polynomial regression. Our Bayesian strategy adopts a Metropolis–Hastings within Gibbs sampler. We improve the sampler by using an ancillarity-sufficiency interweaving strategy and adaptive Markov chain Monte Carlo. We introduce a profile likelihood of the time delay as an approximation of its marginal posterior distribution. The Bayesian and profile likelihood approaches complement each other, producing almost identical results; the Bayesian method is more principled but the profile likelihood is simpler to implement. We demonstrate our estimation strategy using simulated data of doubly- and quadruply-lensed quasars, and observed data from quasars Q0957$+$561 and J1029$+$2623.

Citation

Download Citation

Hyungsuk Tak. Kaisey Mandel. David A. van Dyk. Vinay L. Kashyap. Xiao-Li Meng. Aneta Siemiginowska. "Bayesian estimates of astronomical time delays between gravitationally lensed stochastic light curves." Ann. Appl. Stat. 11 (3) 1309 - 1348, September 2017. https://doi.org/10.1214/17-AOAS1027

Information

Received: 1 February 2016; Revised: 1 January 2017; Published: September 2017
First available in Project Euclid: 5 October 2017

zbMATH: 1380.62275
MathSciNet: MR3709561
Digital Object Identifier: 10.1214/17-AOAS1027

Keywords: adaptive MCMC , ancillarity-sufficiency interweaving strategy , Gibbs sampler , Gravitational lensing , J1029$+$2623 , LSST , microlensing , Ornstein–Uhlenbeck process , profile likelihood , Q0957$+$561 , quasar

Rights: Copyright © 2017 Institute of Mathematical Statistics

Vol.11 • No. 3 • September 2017
Back to Top