Bayesian Estimates of Astronomical Time Delays between Gravitationally Lensed Stochastic Light Curves

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.

Accepted for publication in the Annals of Applied Statistics

References in corpus (11)

• Modeling the Time Variability of SDSS Stripe 82 Quasars as a Damped Random Walk
• Is Quasar Optical Variability a Damped Random Walk?
• Time Delay Cosmography
• Strong Lens Time Delay Challenge: II. Results of TDC1
• Overall Objective Priors
• COSMOGRAIL XVI: Time delays for the quadruply imaged quasar DES J0408-5354 with high-cadence photometric monitoring
• Discovery of 5000 Active Galactic Nuclei behind the Magellanic Clouds
• $H_0$ from ten well-measured time delay lenses
• A Two-Year Time Delay for the Lensed Quasar SDSS J1029+2623
• Bayesian Estimates of Astronomical Time Delays between Gravitationally Lensed Stochastic Light Curves
• Rgbp: An R Package for Gaussian, Poisson, and Binomial Random Effects Models with Frequency Coverage Evaluations

Cited by in corpus (17)

• H0LiCOW - IX. Cosmographic analysis of the doubly imaged quasar SDSS 1206+4332 and a new measurement of the Hubble constant
• Time Delay Cosmography
• Strong lensing time-delay cosmography in the 2020s
• Time Delay Lens Modelling Challenge
• Astrometric requirements for strong lensing time-delay cosmography
• The Magnificent Five Images of Supernova Refsdal: Time Delay and Magnification Measurements
• A Repelling-Attracting Metropolis Algorithm for Multimodality
• Modeling Stochastic Variability in Multi-Band Time Series Data
• How proper are Bayesian models in the astronomical literature?
• Bayesian Estimates of Astronomical Time Delays between Gravitationally Lensed Stochastic Light Curves
• Robust and Accurate Inference via a Mixture of Gaussian and Student's t Errors
• Rapid Optical Variations Correlated with X-rays in the 2015 Second Outburst of V404 Cygni (GS 2023$+$338)
• nazgul: A statistical approach to gamma-ray burst localization. Triangulation via non-stationary time-series models
• TD-CARMA: Painless, accurate, and scalable estimates of gravitational-lens time delays with flexible CARMA processes
• GausSN: Bayesian Time-Delay Estimation for Strongly Lensed Supernovae
• Six Maxims of Statistical Acumen for Astronomical Data Analysis
• A Robust Bayesian Meta-Analysis for Estimating the Hubble Constant via Time Delay Cosmography