Data Analysis problem set
OlympiadIOAA 2023 | Chorzów, Poland | Data Analysis 2 | 75 pt

2023 DA2 | Weighing an isolated black hole with astrometric microlensing

In 2022, two groups reported an isolated black hole from the microlensing event OGLE-2011-BLG-0462. Here we analyse Hubble Space Telescope data to reproduce their result.

In gravitational microlensing, light from a distant source star is bent and magnified by an intervening lens. The angular Einstein radius depends on the lens mass MM and distance DℓD_\ell:

θE=4GMc2Ds−DℓDsDℓ,\theta_E = \sqrt{\frac{4GM}{c^2}\frac{D_s - D_\ell}{D_sD_\ell}},

where Ds≈8D_s\approx8 kpc for a source in the Galactic bulge.

(a) Calculate θE\theta_E in mas for a 1 M⊙1\,M_\odot lens at 1 kpc. (2)

When the lens and source are separated by θ=u θE\theta = u\,\theta_E, two images form on the line through them, at angular distances θ±=12(u±u2+4)θE\theta_\pm = \tfrac12\left(u\pm\sqrt{u^2+4}\right)\theta_E from the lens, with magnifications

A±=12(u2+2uu2+4±1).A_\pm = \frac12\left(\frac{u^2+2}{u\sqrt{u^2+4}}\pm1\right).

(b) Telescopes cannot resolve the two images, only their brightness-weighted centroid. Derive the angular separation θc\theta_c of the centroid from the lens as a function of uu and θE\theta_E. (8)

(c) Derive the deflection Δθ=θc−θ\Delta\theta = \theta_c - \theta of the centroid from the unlensed source position. What is it when u≈0u\approx0? (4)

The source and lens move relative to each other, so both magnification and centroid position change with time. The OGLE light curve gives an Einstein timescale tE=247t_E = 247 days (time to move one θE\theta_E), a peak on 21 July 2011 (HJD 2455763), and a minimum separation u0≈0u_0\approx0. Assume rectilinear relative motion. The table gives HST positions of the source relative to background objects.

(d) Plot the East and North positions against time. (10)

(e) The observed motion is the source's rectilinear proper motion plus the microlensing deflection. Calculate the proper motion (mas/yr) and its uncertainty in East and North. (8)

(f) Subtract the proper motion, then calculate and plot the total deflection against uu. Neglect the uncertainty of the proper motion. (20)

(g) Determine θE\theta_E and its uncertainty. (Hint: linearise Δθ\Delta\theta.) (16)

(h) Including the Earth's orbital motion gives the microlensing parallax πE=(πℓ−πs)/θE=0.095±0.009\pi_E = (\pi_\ell - \pi_s)/\theta_E = 0.095\pm0.009. Find the lens mass and its uncertainty. (7)

HST astrometry of the source star
HJDEast (mas)North (mas)
2455765.22.58 ± 0.137.29 ± 0.16
2455865.72.32 ± 0.125.44 ± 0.24
2456179.70.46 ± 0.141.62 ± 0.08
2456195.80.88 ± 0.361.56 ± 0.77
2456426.2−1.02 ± 0.21−0.94 ± 0.12
2456587.7−2.04 ± 0.07−1.88 ± 0.40
2456956.6−4.54 ± 0.25−5.16 ± 0.29
2457995.2−11.14 ± 0.12−15.14 ± 0.17

Saved in this browser only, no account needed.

Reading resources

Material that covers the physics behind this problem.