# An Isolated Stellar-mass Black Hole Detected through Astrometric Microlensing

Kailash C. Sahu, Jay Anderson, Stefano Casertano, Howard E. Bond, Andrzej Udalski, Martin Dominik, Annalisa Calamida, Andrea Bellini, Thomas M. Brown, Marina Rejkuba, et al.

► To cite this version:

Kailash C. Sahu, Jay Anderson, Stefano Casertano, Howard E. Bond, Andrzej Udalski, et al. An Isolated Stellar-mass Black Hole Detected through Astrometric Microlensing. The Astrophysical Journal, 2022, 933, 10.3847/1538-4357/ac739e. insu-03748250

HAL Id: insu-03748250

https://insu.hal.science/insu-03748250

Submitted on 9 Aug 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

An Isolated Stellar-mass Black Hole Detected through Astrometric Microlensing*

Kailash C. Sahu<fnref n="1" />,⁷⁰ ID, Jay Anderson<fnref n="1" /> ID, Stefano Casertano<fnref n="1" />, Howard E. Bond<fnref n="1" />,<fnref n="2" /> ID, Andrzej Udalski<fnref n="3" />,⁷¹ ID, Martin Dominik<fnref n="4" />,⁷² ID, Annalisa Calamida<fnref n="1" /> ID, Andrea Bellini<fnref n="1" /> ID, Thomas M. Brown<fnref n="1" /> ID, Marina Rejkuba<fnref n="5" /> ID, Varun Bajaj<fnref n="1" />, Noé Kains<fnref n="6" />,⁷³ ID, Henry C. Ferguson<fnref n="1" /> ID, Chris L. Fryer<fnref n="7" /> ID, Philip Yock<fnref n="8" /> ID, Przemek Mróz<fnref n="3" />, Szymon Kozłowski<fnref n="3" /> ID, Paweł Pietrukowicz<fnref n="3" /> ID, Radek Poleski<fnref n="3" /> ID, Jan Skowron<fnref n="3" /> ID, Igor Soszyński<fnref n="3" /> ID, Michał K. Szymański<fnref n="3" /> ID, Krzysztof Ulaczyk<fnref n="3" />,<fnref n="9" /> ID, Łukasz Wyrzykowski<fnref n="3" /> ID

(OGLE Collaboration),

Richard K. Barry<fnref n="10" /> ID, David P. Bennett<fnref n="10" />,<fnref n="11" /> ID, Ian A. Bond<fnref n="12" />, Yuki Hirao<fnref n="13" /> ID, Stela Ishitani Silva<fnref n="10" />,<fnref n="14" /> ID, Iona Kondo<fnref n="13" /> ID, Naoki Koshimoto<fnref n="10" /> ID, Clément Ranc<fnref n="15" /> ID, Nicholas J. Rattenbury<fnref n="16" /> ID, Takahiro Sumi<fnref n="13" /> ID, Daisuke Suzuki<fnref n="13" /> ID, Paul J. Tristram<fnref n="17" />, Aikaterini Vandorou<fnref n="10" />,<fnref n="11" /> ID

(MOA Collaboration),

Jean-Philippe Beaulieu<fnref n="18" />,<fnref n="19" />, Jean-Baptiste Marquette<fnref n="20" />, Andrew Cole<fnref n="18" /> ID, Pascal Fouqué<fnref n="21" /> ID, Kym Hill<fnref n="18" />, Stefan Dieters<fnref n="18" />, Christian Coutures<fnref n="19" />, Dijana Dominis-Prester<fnref n="22" /> ID, Clara Bennett<fnref n="23" />, Etienne Bachelet<fnref n="24" /> ID, John Menzies<fnref n="25" />, Michael Albrow<fnref n="26" /> ID, Karen Pollard<fnref n="26" />

(PLANET Collaboration),

Andrew Gould<fnref n="27" />,<fnref n="28" />, Jennifer C. Yee<fnref n="29" /> ID, William Allen<fnref n="30" />, Leonardo A. Almeida<fnref n="31" />,³² ID, Grant Christie³³, John Drummond³⁴,³⁵, Avishay Gal-Yam³⁶ ID, Evgeny Gorbikov³⁷, Francisco Jablonski³⁸, Chung-Uk Lee³⁹ ID, Dan Maoz⁴⁰, Ilan Manulis⁴¹, Jennie McCormick⁴², Tim Natusch³³,⁴³, Richard W. Pogge<fnref n="28" /> ID, Yossi Shvartzvald⁴¹ ID

(μFUN Collaboration),

Uffe G. Jørgensen⁴⁴ ID, Khalid A. Alsubai⁴⁵, Michael I. Andersen⁴⁶, Valerio Bozza⁴⁷,⁴⁸ ID, Sebastiano Calchi Novati⁴⁹, Martin Burgdorf⁵⁰ ID, Tobias C. Hinse⁵¹,⁵² ID, Markus Hundertmark<fnref n="15" />,⁷³ ID, Tim-Oliver Husser⁵³, Eamonn Kerins⁵⁴ ID, Penelope Longa-Peña⁵⁵, Luigi Mancini<fnref n="27" />,⁵⁶,⁵⁷ ID, Matthew Penny⁵⁸, Sohrab Rahvar⁵⁹ ID, Davide Ricci⁶⁰ ID, Sedighe Sajadian⁶¹ ID, Jesper Skottfelt⁶² ID, Colin Snodgrass⁶³,⁷³ ID, John Southworth⁶⁴ ID, Jeremy Tregloan-Reed⁶⁵, Joachim Wambsganss<fnref n="15" /> ID, Olivier Wertz⁶⁶

(MiNDSTEp Consortium),

and

Yiannis Tsapras<fnref n="15" /> ID, Rachel A. Street<fnref n="24" /> ID, D. M. Bramich⁶⁷,⁶⁸, Keith Horne<fnref n="4" />,⁷⁴ ID, and Iain A. Steele⁶⁹

(RoboNet Collaboration)

<footnote n="1">¹ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA; ksahu@stsci.edu</footnote>

<footnote n="2">² Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA</footnote>

<footnote n="3">³ Astronomical Observatory, University of Warsaw, Al. Ujazdowskie 4, 00-478 Warszawa, Poland</footnote>

<footnote n="4">⁴ University of St Andrews, Centre for Exoplanet Science, SUPA School of Physics & Astronomy, North Haugh, St Andrews, KY16 9SS, UK</footnote>

<footnote n="5">⁵ European Southern Observatory, Karl-Schwarzschild-Straße 2, D-85748 Garching bei München, Germany</footnote>

<footnote n="6">⁶ Department of Physics & Astronomy, Barnard College, Columbia University, 3009 Broadway, New York, NY 10027, USA</footnote>

<footnote n="7">⁷ Center for Theoretical Astrophysics, Los Alamos National Laboratory, Los Alamos, NM 87545, USA</footnote>

<footnote n="8">⁸ Department of Physics, University of Auckland, Auckland, New Zealand</footnote>

<footnote n="9">⁹ Department of Physics, University of Warwick, Gibbet Hill Road, Coventry, CV4 7AL, UK</footnote>

<footnote n="10">¹⁰ Code 667, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA</footnote>

<footnote n="11">¹¹ Department of Astronomy, University of Maryland, College Park, MD 20742, USA</footnote>

<footnote n="12">¹² Institute of Natural and Mathematical Sciences, Massey University, Auckland 0745, New Zealand</footnote>

<footnote n="13">¹³ Department of Earth and Space Science, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan</footnote>

<footnote n="14">¹⁴ Department of Physics, The Catholic University of America, Washington, DC 20064, USA</footnote>

<footnote n="15">¹⁵ Zentrum für Astronomie der Universität Heidelberg, Astronomisches Rechen-Institut, Mönchhofstr. 12-14, D-69120 Heidelberg, Germany</footnote>

<footnote n="16">¹⁶ Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand</footnote>

<footnote n="17">¹⁷ University of Canterbury Mt. John Observatory, P.O. Box 56, Lake Tekapo 8770, New Zealand</footnote>

<footnote n="18">¹⁸ School of Natural Sciences, University of Tasmania, Private Bag 37 Hobart, Tasmania 7001, Australia</footnote>

<footnote n="19">¹⁹ Sorbonne Université, UPMC Univ Paris 6 et CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98 bis bd Arago, F-75014 Paris, France</footnote>

<footnote n="20">²⁰ Laboratoire d'astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, allée Geoffroy SaintHilaire, F-33615 Pessac, France</footnote>

<footnote n="21">²¹ Université de Toulouse, UPS-OMP, IRAP, Toulouse, France</footnote>

<footnote n="22">²² Faculty of Physics, University of Rijeka, HR-51000 Rijeka, Croatia</footnote>

<footnote n="23">²³ Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA</footnote>

<footnote n="24">²⁴ Las Cumbres Observatory Global Telescope Network, 6740 Cortona Drive, Suite 102, Goleta, CA 93117, USA</footnote>

<footnote n="25">²⁵ South African Astronomical Observatory, PO Box 9, Observatory 7935, South Africa</footnote>

<footnote n="26">²⁶ University of Canterbury, Department of Physics & Astronomy, Private Bag 4800, Christchurch 8020, New Zealand</footnote>

<footnote n="27">²⁷ Max-Planck-Institute for Astronomy, Königstuhl 17, D-69117 Heidelberg, Germany</footnote>

<footnote n="28">²⁸ Department of Astronomy, Ohio State University, 140 W. 18th Avenue, Columbus, OH 43210, USA</footnote>

<footnote n="29">²⁹ Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA</footnote>

<footnote n="30">³⁰ Vintage Lane Observatory, Blenheim, New Zealand</footnote>

<footnote n="31">³¹ Escola de Ciências e Tecnologia, Universidade Federal do Rio Grande do Norte, Natal—RN, 59072-970, Brazil</footnote>

<footnote>* This research is based in part on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555.</footnote>
³² Programa de Pós-Graduação em Física, Universidade do Estado do Rio Grande do Norte, Mossoró—RN, 59610-210, Brazil

³³ Auckland Observatory, Auckland, New Zealand

³⁴ Possum Observatory, Patutahi, New Zealand

³⁵ Centre for Astrophysics, University of Southern Queensland, Toowoomba, Queensland 4350, Australia

³⁶ Benoziyo Center for Astrophysics, Weizmann Institute of Science, 76100 Rehovot, Israel

³⁷ School of Physics and Astronomy, Raymond and Beverley Sackler Faculty of Exact Sciences, Tel-Aviv University, Tel Aviv 69978, Israel

³⁸ Instituto Nacional de Pesquisas Espaciais, Astrophysics Division, Sao Jose dos Campos, Brazil

³⁹ Korea Astronomy and Space Science Institute, Daejon 34055, Republic of Korea

⁴⁰ School of Physics and Astronomy, Tel-Aviv University, Tel-Aviv 6997801, Israel

⁴¹ Department of Particle Physics and Astrophysics, Weizmann Institute of Science, 76100 Rehovot, Israel

⁴² Farm Cove Observatory, Centre for Backyard Astrophysics, Pakuranga, Auckland, New Zealand

⁴³ Institute for Radio Astronomy and Space Research (IRASR), AUT University, Auckland, New Zealand

⁴⁴ Centre for ExoLife Sciences, Niels Bohr Institute, University of Copenhagen, Øster Voldgade 5, 1350 Copenhagen, Denmark

⁴⁵ Qatar Environment and Energy Research Institute (QEERI), HBKU, Qatar Foundation, Doha, Qatar

⁴⁶ Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, DK-2100 Copenhagen, Denmark

⁴⁷ Dipartimento di Fisica “E.R. Caianiello,” Università di Salerno, Via Giovanni Paolo II 132, I-84084 Fisciano, Italy

⁴⁸ Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, Napoli, Italy

⁴⁹ IPAC, Mail Code 100-22, Caltech, 1200 E. California Blvd., Pasadena, CA 91125, USA

⁵⁰ Universität Hamburg, Faculty of Mathematics, Informatics and Natural Sciences, Department of Earth Sciences, Meteorological Institute, Bundesstraße 55, D-20146 Hamburg, Germany

⁵¹ Institute of Astronomy, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudziadzka 5, 87-100 Toruń, Poland

⁵² Chungnam National University, Department of Astronomy, Space Science and Geology, Daejeon, Republic of Korea

⁵³ Institut fur Astrophysik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, D-37077 Göttingen, Germany

⁵⁴ Jodrell Bank Centre for Astrophysics, Alan Turing Building, University of Manchester, Manchester, M13 9PL, UK

⁵⁵ Centro de Astronomia, Universidad de Antofagasta, Avenida Angamos 601, Antofagasta 1270300, Chile

⁵⁶ Department of Physics, University of Rome “Tor Vergata,” Via della Ricerca Scientifica 1, I-00133 Roma, Italy

⁵⁷ INAF—Astrophysical Observatory of Turin, Via Osservatorio 20, I-10025 Pino Torinese, Italy

⁵⁸ Louisiana State University, 261-B Nicholson Hall, Tower Drive, Baton Rouge, LA 70803-4001, USA

⁵⁹ Department of Physics, Sharif University of Technology, PO Box 11155-9161, Tehran, Iran

⁶⁰ INAF—Padova Astronomical Observatory, Vicolo dell'Osservatorio 5, I-35122 Padova, Italy

⁶¹ Department of Physics, Isfahan University of Technology, Isfahan 84156-83111, Iran

⁶² Centre for Electronic Imaging, Department of Physical Sciences, The Open University, Milton Keynes, MK7 6AA, UK

⁶³ Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh, EH9 3HJ, UK

⁶⁴ Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK

⁶⁵ Instituto de Investigación en Astronomía y Ciencias Planetarias, Universidad de Atacama, Copayapu 485, Copiapó, Atacama, Chile

⁶⁶ Space Sciences, Technologies, and Astrophysics Research (STAR) Institute, University of Liège, Liège, Belgium

⁶⁷ Center for Astro, Particle, and Planetary Physics, New York University Abu Dhabi, P.O. Box 129188, Saadiyat Island, Abu Dhabi, UAE
⁶⁸ Division of Engineering, New York University Abu Dhabi, P.O. Box 129188, Saadiyat Island, Abu Dhabi, UAE

⁶⁹ Astrophysics Research Institute, Liverpool John Moores University, Liverpool, L3 5RF, UK

Received 2022 January 28; revised 2022 May 19; accepted 2022 May 24; published 2022 July 6

Abstract

We report the first unambiguous detection and mass measurement of an isolated stellar-mass black hole (BH). We used the Hubble Space Telescope (HST) to carry out precise astrometry of the source star of the long-duration ($t_E \simeq 270$ days), high-magnification microlensing event MOA-2011-BLG-191/OGLE-11-462 (hereafter designated as MOA-11-191/OGLE-11-462), in the direction of the Galactic bulge. HST imaging, conducted at eight epochs over an interval of 6 yr, reveals a clear relativistic astrometric deflection of the background star's apparent position. Ground-based photometry of MOA-11-191/OGLE-11-462 shows a parallactic signature of the effect of Earth's motion on the microlensing light curve. Combining the HST astrometry with the ground-based light curve and the derived parallax, we obtain a lens mass of $7.1 \pm 1.3 M_\odot$ and a distance of $1.58 \pm 0.18$ kpc. We show that the lens emits no detectable light, which, along with having a mass higher than is possible for a white dwarf or neutron star, confirms its BH nature. Our analysis also provides an absolute proper motion for the BH. The proper motion is offset from the mean motion of Galactic disk stars at similar distances by an amount corresponding to a transverse space velocity of $\sim 45$ km s$^{-1}$, suggesting that the BH received a “natal kick” from its supernova explosion. Previous mass determinations for stellar-mass BHs have come from radial velocity measurements of Galactic X-ray binaries and from gravitational radiation emitted by merging BHs in binary systems in external galaxies. Our mass measurement is the first for an isolated stellar-mass BH using any technique.

<figure><img src="image_3.png" /></figure>

<footnote n="70">⁷⁰ PLANET Collaboration, MiNDSTEP Consortium.</footnote>

<footnote n="71">⁷¹ OGLE Collaboration.</footnote>

<footnote n="72">⁷² MiNDSTEP Consortium, RoboNet Collaboration.</footnote>

<footnote n="73">⁷³ RoboNet Collaboration.</footnote>

<footnote n="74">⁷⁴ PLANET Collaboration.</footnote>
*Unified Astronomy Thesaurus concepts: Black holes (162); Gravitational microlensing (672)*

*Supporting material: data behind figures*

# 1. Measuring the Masses of Black Holes

## 1.1. Black Holes in Binary Systems

Stars with initial masses greater than ~20 $M_\odot$ are expected to end their lives as black holes (BHs; e.g., Fryer & Kalogera 2001; Woosley et al. 2002; Heger et al. 2003; Spera et al. 2015; Sukhbold et al. 2016). Objects of these masses constitute roughly 0.1% of all stars, leading to the expectation that the Galaxy should now contain of the order of ≈10⁸ BHs (Shapiro & Teukolsky 1983; van den Heuvel 1992; Brown & Bethe 1994; Samland 1998).

However, the actual detection of stellar-mass BHs is observationally challenging, and determining their masses is even more so. BHs have been identified in the Galaxy and Local Group through X-ray emission due to accretion in short-period binary systems, most of them soft X-ray transients. In such cases, dynamical masses of the BHs can be measured or estimated through radial velocity measurements and light-curve modeling for the optical companion stars (the techniques are reviewed by Remillard & McClintock 2006; Casares & Jonker 2014). Masses of nearly two dozen BHs in X-ray binary systems have been determined using these methods, with varying degrees of precision. These “electromagnetically measured” BH masses show a distribution peaking near 7–8 $M_\odot$, with few if any below ~5 $M_\odot$ (e.g., Özel et al. 2010; Farr et al. 2011; Kreidberg et al. 2012; Corral-Santana et al. 2016). This suggests that a “mass gap” exists between the lowest-mass BHs and the highest measured masses of neutron stars (NSs) in binary radio pulsars of ~2.1–2.3 $M_\odot$ (Linares et al. 2018; Cromartie et al. 2020). Recently, however, a few nonaccreting or weakly accreting BHs have been discovered in longer-period spectroscopic binaries in the field (e.g., Thompson et al. 2019; Jayasinghe et al. 2021) and in globular clusters (Giesers et al. 2019), lying in the NS–BH gap with dynamical masses of ~3–4.5 $M_\odot$. Precision astrometry of nearby stars by Gaia shows the promise of detecting additional wide-binary systems containing quiescent BHs (e.g., Chawla et al. 2021; Janssens et al. 2022, and references therein) and measuring their masses. At the high-mass end, the BH mass distribution falls off above ~10 $M_\odot$, and very few electromagnetic BH masses are known above ~15 $M_\odot$, the only exception in the Milky Way being an updated mass determination of 21.1 ± 2.2 $M_\odot$ for the BH in Cygnus X-1 (Miller-Jones et al. 2021). Among extragalactic X-ray binaries, BH masses as high as 15.65 ± 1.45 $M_\odot$ and 17 ± 4 $M_\odot$ have been reported for M33 X-7 (Orosz et al. 2007) and NGC 300 X-1 (Binder et al. 2021), respectively. An even higher mass of at least 23.1 $M_\odot$ was reported for the compact object in IC 10 X-1 (Silverman & Filippenko 2008), but this has been questioned (Laycock et al. 2015).

The first detections of gravitational waves (GWs) by the Laser Interferometry Gravitational Wave Observatory (LIGO) and Virgo Collaboration (Abbott et al. 2016) revealed a population of massive merging binary BHs, BH–NS pairs, and binary NSs at extragalactic distances. In the source catalogs from the third observing run of the Advanced LIGO and Advanced Virgo Collaboration (Abbott et al. 2021a, 2021b), the inferred masses of the BHs among the premerger systems range from ≈6 to 95 $M_\odot$, with two low-mass outliers among

the secondary components at ~2.6 and 2.8 $M_\odot$, which could be either BHs or NSs.

## 1.2. Isolated Black Holes

The electromagnetic and GW mass measurements described above are all for BHs in binary systems, including those undergoing mass accretion or mergers. However, there are reasons to believe that a substantial fraction of stellar-mass BHs are single, rather than belonging to binaries. First, about 30% of massive stars are born single (see Sana et al. 2012; de Mink et al. 2014). Second, in a close-binary system, the pair may enter into a common envelope and merge before the supernova (SN) explosion (e.g., Fryer et al. 1999; Zhang & Fryer 2001; Tutukov et al. 2011; Dominik et al. 2012). Lastly, in a wide binary, the “natal kick” imparted to the companion by the SN event may be large enough to detach the two components, producing an isolated BH (e.g., Tauris & van den Heuvel 2006; Belczynski et al. 2016). Being less altered by interactions with companions, single BHs potentially provide a more direct probe of BH formation than those in binaries.

Isolated BHs are extremely difficult to detect directly. They emit no light of their own, and the accretion rate from the interstellar medium is generally likely to be too low to produce detectable X-ray or radio emission (see, however, Agol & Kamionkowski 2002; Fender et al. 2013; Tsuna & Kawanaka 2019; Scarcella et al. 2021, for the case of isolated BHs in dense environments). In fact, until now, no isolated stellar-mass BH has ever been unambiguously found within our Galaxy or elsewhere.

Microlensing is the only available method for measuring the masses of isolated BHs. *Astrometric microlensing*—the relativistic deflection of the apparent position of a background star when a compact object passes in front of it—provides a direct method for measuring the masses of BH lenses. High spatial resolution interferometric observations of microlensing events (Dong et al. 2019; Zang et al. 2020) and observations of rare events where the lens passes over the surface of the source (Yoo et al. 2004) can also yield the masses of BH lenses. In this paper, we describe how the technique of astrometric microlensing is used to determine masses. We discuss our ongoing program of astrometric measurements of microlensing events with the Hubble Space Telescope (HST). Then, we report the first detection of an isolated BH and our measurement of its mass.

# 2. Measuring the Masses of Isolated Black Holes with Astrometric Microlensing

## 2.1. Microlensing Events and Black Hole Candidates

A microlensing event occurs when a star or compact object (the lens) passes almost exactly in front of a background star (the source). As predicted by general relativity (Einstein 1936), the lens magnifies the image of the source, producing an apparent amplification of its brightness. The lens also slightly shifts the apparent position of the source (Miyamoto & Yoshii 1995; Høg et al. 1995; Walker 1995)—an analog of the deflection of stellar images during the 1919 solar eclipse
(Dyson et al. 1920), which provided support for the general theory of relativity.

Microlensing survey programs, including OGLE (Udalski et al. 2015), MOA (Bond et al. 2001), and KMTNet (Kim et al. 2016), carry out photometric monitoring of rich stellar fields in the Galactic bulge. These surveys typically detect >2000 events toward the Galactic bulge annually. To date, more than 30,000 microlensing events have been discovered and monitored by these survey programs.

The characteristic scale of gravitational microlensing is provided by the angular Einstein radius $\theta_E$, given by

$$ \theta_E \equiv \sqrt{\frac{4GM_{\text{lens}}}{c^2} \frac{\pi_{\text{LS}}}{1 \text{ au}}}, \qquad (1) $$

where $M_{\text{lens}}$ is the mass of the lens and

$$ \pi_{\text{LS}} \equiv \pi_{\text{L}} - \pi_{\text{S}} = (1 \text{ au}) \left( \frac{1}{D_{\text{L}}} - \frac{1}{D_{\text{S}}} \right) \qquad (2) $$

is the relative lens-source parallax, with $D_L$ and $D_S$ being the distances from the observer to the lens and to the source, respectively.

The Einstein radius $\theta_E$, however, cannot be obtained directly from the magnification light curve, whose only characteristic that carries a physical dimension is the timescale. Given the relative proper motion $\mu_{LS}$ between lens and source, we can straightforwardly define $t_E^*$, the time for the source to traverse an angular distance of $\theta_E$ in the barycentric reference frame, as $t_E^* = \theta_E/\mu_{LS}$ (more details in Section 2.4). The distribution of the timescale $t_E^*$ of the observed events peaks around 25 days, with $t_E^*$ ranging from a fraction of a day to several hundred days (e.g., Wyrzykowski et al. 2015).

If BHs constitute a small but nonnegligible fraction of the total stellar mass of the Galaxy, as described above, then a few of the observed microlensing events are expected to be due to BHs. Equation (1) shows that the angular Einstein radius $\theta_E$ is proportional to the square root of the lens mass. Hence, all else being equal, events due to massive compact objects would preferentially tend to be characterized by longer event durations ($t_E \gtrsim 150$ days), combined with an apparent lack of light contribution from the lens. However, a degeneracy between lens mass and proper motion remains. Thus, a long-duration event with no light contribution from the lens could arise from a high-mass, nonluminous BH lens with a large Einstein radius —but it could alternatively be due simply to an unusually slow-moving, faint, low-mass ordinary star. If the lens were a luminous massive star, this would generally be recognizable through the contribution of its light, particularly when observations are available in two different bandpasses.

Indeed, several long-duration OGLE and MOA microlensing events have been suggested as being due to BHs (e.g., Bennett et al. 2002; Mao et al. 2002; Minniti et al. 2015; Wyrzykowski & Mandel 2020). However, these claims remain statistical in nature, being derived from assumptions about the transverse velocity distributions.

The degeneracy between mass and relative velocity can be lifted if precise astrometry is added to the photometry of the microlensing event. The size of the expected astrometric shift is small—of the order of milliarcseconds—but it is proportional to the angular Einstein radius $\theta_E$, and therefore if this small shift can be measured, then the mass of the lens can be determined unambiguously, as described in detail below.

## 2.2. Photometric Microlensing

Photometric microlensing is the apparent transient brightening that results as a background source passes almost directly behind a foreground lens (see reviews by Paczyński 1996; Gaudi 2012, and Tsapras 2018).

With $\theta_S$ and $\theta_L$ denoting the angular positions of the source and lens as seen by the observer, we can define a dimensionless source–lens separation

$$ \mathbf{u} \equiv \frac{\theta_S - \theta_L}{\theta_E}. \qquad (3) $$

As the lens intervenes near the line of sight from the observer to the source, the gravitational bending of light leads to a time-varying magnification

$$ A(u) = \frac{u^2 + 2}{u\sqrt{u^2 + 4}}, \qquad (4) $$

which depends solely on $u = |\mathbf{u}|$. This expression holds as long as the finite angular size of the source star can be neglected, which we will assume in the following discussion. We will demonstrate in Section 9.3 that this is a valid approximation for the case analyzed in this paper.

If $F_S$ is the intrinsic source flux and $F_B$ the background flux contributed by any other objects not resolved from the observed source star, the observed flux of the target for a specific telescope and filter is given by

$$ F(t) = F_S A[u(t)] + F_B = F_{\text{base}} \frac{A[u(t)] + g}{1 + g}, \qquad (5) $$

where $F_{\text{base}} \equiv F_S + F_B$ is the baseline flux and $g \equiv F_B/F_S$ is the specific blend ratio.

## 2.3. Astrometric Microlensing

Microlensing also produces an astrometric shift of the apparent position of the source. If we assume that we can observe the centroid of light formed by the images of the source without any contribution from other bodies such as the lens or other neighboring stars, its shift is described by the vector

$$ \delta(\mathbf{u}) = \frac{\mathbf{u}}{u^2 + 2} \theta_E \qquad (6) $$

(Gould 1992; Paczyński 1998). In contrast to the photometric microlensing signature, the astrometric signature is explicitly proportional to the angular Einstein radius $\theta_E$.

Moreover, while the magnification diverges (for a pointlike source) as $u \to 0$, the astrometric shift becomes maximal for $u = \sqrt{2}$. The light magnification falls rapidly with increasing $u$: for large separations, $u \gg 1$, the brightness enhancement, $A(u) - 1$, falls as $1/u^4$. On the other hand, the centroid shift, given by Equation (6), decreases more slowly with $u$, and for large separations it falls only as $1/u$. The astrometric perturbation thus has a considerably longer duration than the photometric signal. For more details, see Dominik & Sahu (2000), Sahu et al. (2014), and Bramich (2018).

The photometric and astrometric signatures of a microlensing event are connected because they arise from the same source–lens trajectory, $\mathbf{u}(t)$. Specifically, if a fit to the light curve of a microlensing event already yields $u$ as a function of time, the astrometric data then provide a direct measurement of
the angular Einstein radius $\theta_E$, as well as the orientation angle of the trajectory.

## 2.4. Parallax Effect and Proper Motion

For constructing the source–lens trajectory $\mathbf{u}(t)$, we have to consider the proper motions $\mu_S$ and $\mu_L$ of the source and lens objects, as well as their parallaxes, $\pi_S$ and $\pi_L$.

Let $\gamma(t)$ (1 au) denote the projection of Earth's orbit onto a plane perpendicular to the direction toward the source star. The apparent geocentric positions of source and lens star are then given by (see An et al. 2002; Gould 2004)

$$ 
\begin{aligned}
\theta_S(t) &= \theta_{S,0} + (t - t_0)\mu_S - \pi_S \gamma(t), \\
\theta_L(t) &= \theta_{L,0} + (t - t_0)\mu_L - \pi_L \gamma(t),
\end{aligned}
\quad (7) $$

so that for $\theta(t) \equiv \theta_S(t) - \theta_L(t)$ one finds

$$ \theta(t) = (\theta_S - \theta_L)_0 - (t - t_0)\mu_{LS} + \pi_{LS} \gamma(t), \quad (8) $$

where $\mu_{LS} \equiv \mu_L - \mu_S$ and $\pi_{LS} \equiv \pi_L - \pi_S$ are the relative proper motion and relative parallax between lens and source, while $(\theta_S - \theta_L)_0 \equiv \theta_{S,0} - \theta_{L,0}$.

Consequently, with the microlensing parallax parameter $\pi_E = \pi_{LS}/\theta_E$, $\mathbf{u}(t) \equiv \theta(t)/\theta_E$ takes the form (see Dominik et al. 2019)

$$ \mathbf{u}(t) = \mathbf{u}_0 + (t - t_0)\dot{\mathbf{u}}_0 + \pi_E \delta\gamma(t), \quad (9) $$

where

$$ \mathbf{u}_0 \equiv \mathbf{u}(t_0) = \frac{(\theta_S - \theta_L)_0}{\theta_E} + \pi_E \gamma(t_0), \quad (10) $$

$$ \dot{\mathbf{u}}_0 \equiv \dot{\mathbf{u}}(t_0) = -\frac{\mu_{LS}}{\theta_E} + \pi_E \dot{\gamma}(t_0), \quad (11) $$

as well as

$$ \delta\gamma(t) = \gamma(t) - \gamma(t_0) - (t - t_0)\dot{\gamma}(t_0). \quad (12) $$

By construction, $\delta\gamma(t_0) = 0$ and $\delta\dot{\gamma}(t_0) = 0$.

By choosing $t_0$ so that $\mathbf{u}_0 \perp \dot{\mathbf{u}}_0$, we can write $\mathbf{u}(t)$ in its components toward northern and eastern directions as

$$ 
\begin{aligned}
u_n(t) &= \frac{t - t_0}{t_E} \cos \psi - u_0 \sin \psi + \pi_E \delta\gamma_n(t), \\
u_e(t) &= \frac{t - t_0}{t_E} \sin \psi + u_0 \cos \psi + \pi_E \delta\gamma_e(t),
\end{aligned}
\quad (13) $$

where $u_0 \equiv |\mathbf{u}_0|$, $t_E = 1/|\dot{\mathbf{u}}_0|$, and $\psi$ denotes the direction angle of $\dot{\mathbf{u}}_0$ measured from north toward east. Alternatively, the trajectory can be parameterized as

$$ 
\begin{aligned}
u_n(t) &= \frac{t - t_0^*}{t_E^*} \cos \psi^* - u_0^* \sin \psi^* + \pi_E \gamma_n(t), \\
u_e(t) &= \frac{t - t_0^*}{t_E^*} \sin \psi^* + u_0^* \cos \psi^* + \pi_E \gamma_e(t),
\end{aligned}
\quad (14) $$

where $u_0^* \equiv |\theta_S(t_0^*) - \theta_L(t_0^*)|/\theta_E$, $t_E^* = \theta_E/|\mu_{LS}|$, and $\psi^*$ denotes the direction angle of $-\mu_{LS} = \mu_S - \mu_L$ measured from north toward east. Note that the starred quantities, $\psi^*$, $t_E^*$, $t_0^*$, and $u_0^*$ refer to parameters in the barycentric reference frame, whereas the corresponding unstarred quantities refer to parameters as seen by an observer on Earth.

This applies to any orthonormal reference frame; the only difference is in the specific angle, e.g., $\psi_{eq}$, $\psi_{ecl}$, and $\psi_{gal}$ (and $\psi_{eq}^*$, $\psi_{ecl}^*$, and $\psi_{gal}^*$) for equatorial, ecliptic, and galactic

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Source–lens trajectory $\mathbf{u}(t)$ as seen by the observer, showing the effect of annual parallax. At epoch $t_0$, the tangent to the source–lens trajectory $\dot{\mathbf{u}}_0$ is not (anti)parallel to the direction of the lens–source proper motion $\hat{\mu}_{LS}$ but differs by $\pi_E \dot{\gamma}(t_0)$ (Equation (11)), related to the orbital velocity of Earth at $t_0$. Consequently, we distinguish the direction angles $\psi$, $\psi^*$, and $\varphi_{LS}$, referring to $\dot{\mathbf{u}}_0$ and $\mp \hat{\mu}_{LS}$, respectively. Furthermore, $t_E = 1/|\dot{\mathbf{u}}_0|$ and $t_E^* = \theta_E/|\mu_{LS}|$.</figcaption></figure>

coordinates, respectively, which are related by a rotation of the coordinate axes at the target position.

While most discussions of photometric microlensing events choose an ecliptic coordinate frame, the observed astrometric data are more easily described in an equatorial coordinate frame, and therefore we adopt the latter in the following analysis. In this frame,

$$ \varphi_{LS} = \psi_{eq}^* + 180^\circ \quad (15) $$

gives the position angle (PA) of the proper motion of the lens with respect to the source $\mu_{LS}$, measured from equatorial north toward east. We illustrate the geometry of the source–lens trajectory $\mathbf{u}(t)$ in Figure 1.

## 2.5. Measuring the Lens Mass

As stated above, the only useful physical parameter in a typical microlensing event is the timescale $t_E^* = \theta_E/\mu_{LS}$, which is the time it takes the source to traverse the radius of the Einstein ring, which itself depends on the lens mass, $M_{lens}$, and the lens–source parallax, $\pi_{LS}$.

However, for long-duration events, the annual parallax tends to lead to prominent departures in the photometric signature (Gould 1992; Alcock et al. 1995), so that a microlensing parallax parameter $\pi_E \equiv \pi_{LS}/\theta_E$ can be inferred. Coincidentally, the BH mass lenses tend to imply such long-duration events. On the other hand, the astrometric signature is proportional to the angular Einstein radius $\theta_E$, so that by combining photometric and astrometric observations, $M_{lens}$, $\pi_{LS}$, and $\mu_{LS}$ become fully decoupled. Specifically, with $\pi_E$ from the photometry and $\theta_E$ from the astrometry, the definition of $\theta_E$, Equation (1), immediately gives us

$$ M_{\text{lens}} = \frac{\theta_E c^2 (1 \text{ au})}{\pi_E 4G} = \frac{\theta_E}{\kappa \pi_E}, \quad (16) $$

where $\kappa = 4G/[c^2(1 \text{ au})] \simeq 8.144 \text{ mas } M_{\odot}^{-1}$.

In the case of microlensing toward the Galactic bulge, the source often lies at the distance of the bulge itself, which can be verified from its baseline position in a color–magnitude diagram (CMD). If spectroscopic observations are available in addition to baseline photometry—as is the case for the event discussed in this paper—a more accurate source distance can be determined. The lens–source relative parallax, $\pi_{LS} = \pi_E \theta_E$, can then be used to estimate the distance to the lens, using Equation (2). As a bonus, the event timescale gives a direct
<figure><img src="image_5.png" /></figure>

**Figure 2.** Astrometric shift (top panel) and light magnification (bottom panel) for a microlensing event produced by a 5 $M_\odot$ BH at a distance of 2 kpc passing in front of a background star at 8 kpc. The assumed minimum impact parameter is $u_0 = 0.05$. The maximum astrometric shift of the source is ~1.4 mas, at $u = \sqrt{2}$, and the maximum magnification is ~20, at the time of closest angular approach. Note the much longer duration for the astrometric shift, compared to that of the light magnification.

measure of the relative transverse velocity of lens and source (which, for stellar remnants, might include “kicks” received in SN explosions). This method thus provides independent measurements of three separate physical parameters of the lens: its mass, distance, and transverse velocity.

## 2.6. Characteristics of Astrometric Deflections

Some features of astrometric deflections under various scenarios are described by Dominik & Sahu (2000). To illustrate a typical case, we show in Figure 2 the calculated astrometric shifts and light magnification for a nominal event of a BH lens of mass 5 $M_\odot$, at a distance of 2 kpc from the Sun, passing in front of a background source situated in the Galactic bulge at a distance of 8 kpc. The closest angular approach is assumed to be at a separation of 0.05$\theta_E$. In this case, the size of the angular Einstein ring is $\theta_E \simeq 4$ mas, so that the maximum astrometric shift is ~1.4 mas (occurring at a separation of $u = \sqrt{2}$) and the maximum light magnification is a factor of ~20 (at closest angular approach).

As Figure 2 illustrates, and as discussed in Section 2.3, the duration of the astrometric deflection is considerably longer than that of the photometric magnification. This makes it necessary to carry out the astrometric measurements over a longer time interval than the photometry. Although the deflection measured at any given epoch provides in principle an estimate of $\theta_E$, it is necessary to observe at multiple epochs in order to separate the shifts caused by microlensing from those caused by the proper motion of the source; observations at a late epoch are particularly useful for this purpose. The figure also shows that the astrometric shift is close to zero at the time of highest magnification; therefore, observations near the photometric peak are also very useful to constrain the source proper motion.

<figure><img src="image_6.png" /></figure>

**Figure 3.** Maximum astrometric shift of a source in the Galactic bulge at 8 kpc as a function of lens mass, for lens distances of 2, 4, and 6 kpc. Since the maximum astrometric shift occurs at a lens-source separation of $u = \sqrt{2}$, while the maximum light magnification occurs at the minimum value of $u$, all high-magnification microlensing events will pass through the point of maximum astrometric deflection. High-S/N imaging with HST allows measurements of astrometric shifts with a precision of ~0.2 mas per observation epoch, shown by the dotted line at the bottom.

In Figure 3, we plot the maximum astrometric shifts for a source in the Galactic bulge at 8 kpc, as functions of lens mass. The lenses are assumed to be located at distances of 2 and 4 kpc (“disk” lenses), and 6 kpc (“bulge” lenses). The dotted line at the bottom shows the nominal astrometric precision of 0.2 mas achievable with high signal-to-noise ratio (S/N) HST imaging, as discussed below. Therefore, the deflection is detectable at $1\sigma$ per epoch for lens masses down to ~0.5 $M_\odot$, except at lens distances larger than 6 kpc. The most favorable situation for a precise mass measurement, of course, would be for a nearby, high-mass lens.

## 2.7. High-precision Astrometry

Although an unambiguous determination of lens mass is possible from a combination of photometry and astrometry, as we have just discussed, the expected astrometric shifts are extremely small, of the order of milliarcseconds or less. HST has demonstrated its capability to carry out submilliarcsecond astrometry through a variety of techniques. For example, high-S/N HST observations of isolated sources were used to achieve submilliarcsecond accuracy, leading to the measurement of proper motions for several distant hypervelocity stars (Brown et al. 2015). A collection of sources was used as probes to achieve an astrometric accuracy of ~12 µas, in order to measure the transverse velocity of M31 (Sohn et al. 2012). Spatial-scan techniques have been used to achieve an astrometric accuracy of ~30 µas (Casertano et al. 2016; Riess et al. 2018) in the trigonometric parallax of Cepheids used for accurate determination of $H_0$ and to measure the distance to the globular cluster NGC 6397 (Brown et al. 2018). Recently, our group used the astrometric microlensing technique to measure the mass of the nearby white dwarf Stein 2051 B, achieving an
<table><caption>**Table 1**
Basic Data for MOA-11-191/OGLE-11-462 Microlensing Event</caption><thead><tr><th>Parameter</th><th>Value</th><th>Sources and Notes<sup>a</sup></th></tr></thead><tbody><tr><td>Event designation (MOA)</td><td>MOA-2011-BLG-191</td><td>(1)</td></tr><tr><td>Event designation (OGLE)</td><td>OGLE-2011-BLG-0462</td><td>(1)</td></tr><tr><td>J2000 R.A., α</td><td>17:51:40.2082</td><td>(2)</td></tr><tr><td>J2000 decl., δ</td><td>-29:53:26.502</td><td>(2)</td></tr><tr><td>Galactic coordinates, (l, b)</td><td>359°86, -1°62</td><td>(2)</td></tr><tr><td>Baseline F606W magnitude</td><td>21.943 ± 0.014</td><td>(3)</td></tr><tr><td>Baseline F814W magnitude</td><td>19.578 ± 0.012</td><td>(3)</td></tr><tr><td>Baseline (F606W - F814W) color</td><td>2.365 ± 0.026</td><td>(3)</td></tr><tr><td>Peak magnification, A<sub>max</sub></td><td>369</td><td>(4)</td></tr><tr><td>Date of peak magnification, t<sub>0</sub></td><td>2011 July 20.825</td><td>(4)</td></tr><tr><td>Timescale, t<sub>E</sub></td><td>270.7 ± 11.2 days</td><td>(4)</td></tr></tbody></table>

**Note.**

<sup>a</sup> Sources and notes: (1) MOA and OGLE websites; the event was first alerted by MOA. (2) This paper, from astrometric analysis in Section 5.2 in Gaia EDR3 frame at average epoch 2013.5. (3) This paper, Vegamag scale, from photometric analysis in Section 5.3. (4) This paper, from Table 6.

astrometric precision of ~0.2 mas per epoch (Sahu et al. 2017). Kains et al. (2017) looked for astrometric deflections in HST observations of 10 microlensing events with timescales of <50 days. They achieved an astrometric precision of 0.2 mas per epoch (but did not detect any deflections). From the ground, Zurlo et al. (2018) used the Very Large Telescope (VLT) to measure the mass of Proxima Centauri through astrometric microlensing. Lu et al. (2016) employed the Keck telescope to look specifically for isolated BHs by monitoring three microlensing events, where they achieved a final positional error of 0.26–0.68 mas. The timescales of those events were 60–160 days, and there were no detections of astrometric deflections.

### 3. In Search of Isolated Black Holes with HST

#### 3.1. Astrometry of Long-duration Microlensing Events

In 2009, we began a multicycle HST program of astrometry of long-duration microlensing events in the direction of the Galactic bulge in order to detect isolated BHs and measure their masses. Our aim is to select events having timescales ≳200 days, light curves showing no evidence for a light contribution by a luminous lens, and preferably a high magnification factor. We then obtain high-resolution HST imaging as the events proceed, in order to measure the astrometric deflections of the background sources. To date we have monitored eight long-duration events. For some of them, there is no clear detection of an astrometric signal, but our data analysis is still in progress, and the results will be discussed in separate publications. In the present paper we analyze and discuss our findings for an event that clearly shows a large astrometric deflection, consistent with a high-mass lens.

#### 3.2. MOA-2011-BLG-191/OGLE-2011-BLG-0462

MOA-2011-BLG-191/OGLE-2011-BLG-0462 (hereafter designated MOA-11-191/OGLE-11-462) was a long-duration and high-magnification microlensing event in the direction of the Galactic bulge. It was discovered independently by both

<figure><img src="image_7.png" /><figcaption>**Figure 4.** HST image in the F814W (l-band) filter of an 8" × 8" region centered on MOA-11-191/OGLE-11-462, obtained at our final epoch in 2017 August. North is at the top, east on the left. Encircled in green is the source star, now returned to baseline luminosity. The site is resolved into the source, a much brighter neighboring star 0"′4 to the W–NW, and several nearby fainter stars. The inner cyan circle has a diameter of 1", corresponding to the typical best seeing in ground-based microlensing survey images; the outer cyan circle's diameter is 2", which is not unusual seeing. The source, bright neighbor, and several fainter stars are generally blended in ground-based frames, and the blending increases with seeing.</figcaption></figure>

MOA and OGLE ground-based microlensing survey programs and announced by both teams nearly simultaneously on 2011 June 2, through their public-alert websites.<fnref n="75" /> The target was also covered by the Wise Microlensing Survey. Table 1 gives details of this remarkable event.

MOA-11-191/OGLE-11-462 occurred in an extremely crowded Galactic bulge field, less than 2° from the Galactic center. The observed peak magnification factor of this event was only about 20 in the ground-based data, but this was strongly diluted by blending with neighboring stars. It soon became apparent, based on findings disseminated through internal communications in the microlensing groups, that the undiluted event actually had an extremely high magnification factor, approaching 400. Blending also made the apparent timescale of the event appear shorter than the actual value, which was inferred to be longer than 200 days. It was clear from the ground-based observations that there was blending, for two reasons. First, the light curve for a typical event has a characteristic shape that is completely determined by the timescale and the maximum magnification, except for distortions due, e.g., to the lens–source relative parallax. The shape of the observed light curve was inconsistent with the expected shape unless the light at baseline was highly diluted by a blend, thus implying that the real magnification was much larger than the observed value. Second, as the source brightened, its centroid position in the ground-based images was seen to change, again consistent with blending with a neighboring star.

<footnote n="75">⁷⁵ MOA alerts: https://www.massey.ac.nz/~iabond/moa/alerts. OGLE alerts: http://ogle.astrouw.edu.pl/ogle4/ews/2011/ews.html.</footnote>
<table><caption>**Table 2**
Journal of HST Wide Field Camera 3 Observations</caption><thead><tr><th>Epoch</th><th>Date</th><th>MJD</th><th>Year</th><th>Proposal ID</th><th>No. Frames in F606W<sup>a</sup></th><th>No. Frames in F814W<sup>a</sup></th></tr></thead><tbody><tr><td>1</td><td>2011 Aug 8</td><td>55,781.7</td><td>2011.600</td><td>GO-12322</td><td>4</td><td>5</td></tr><tr><td>2</td><td>2011 Oct 31</td><td>55,865.2</td><td>2011.829</td><td>GO-12670</td><td>3</td><td>4</td></tr><tr><td>3</td><td>2012 Sep 9</td><td>56,179.2</td><td>2012.689</td><td>GO-12670</td><td>3</td><td>4</td></tr><tr><td>4</td><td>2012 Sep 25</td><td>56,195.3</td><td>2012.733</td><td>GO-12986</td><td>3</td><td>4</td></tr><tr><td>5</td><td>2013 May 13</td><td>56,425.8</td><td>2013.364</td><td>GO-12986</td><td>3</td><td>4</td></tr><tr><td>6</td><td>2013 Oct 22</td><td>56,587.2</td><td>2013.806</td><td>GO-13458</td><td>3</td><td>4</td></tr><tr><td>7</td><td>2014 Oct 26</td><td>56,956.1</td><td>2014.816</td><td>GO-13458</td><td>3</td><td>4</td></tr><tr><td>8</td><td>2017 Aug 29</td><td>57,994.7</td><td>2017.660</td><td>GO-14783</td><td>3</td><td>4</td></tr></tbody></table>

**Note.**

Note that this shift is due simply to blending and scales with the separation of the two stars; it is unrelated to the much smaller relativistic deflection of the source itself, which is discussed below.

Figure 4 shows an 8'' × 8'' region centered on the source, as imaged by us in the F814W (I-band) filter by HST with its Wide Field Camera 3 (WFC3). The source star is encircled in green. A conspicuous neighbor, nearly 20 times brighter than the unmagnified source, lies at a separation of only 0.''4. The cyan circles in the figure have diameters of 1'' and 2'', corresponding to the generally best seeing in the ground-based survey observations and more typical seeing, respectively. Thus, in ground-based images, the source is indeed blended with the bright neighbor and a number of fainter stars, depending on the seeing.

High-magnification events are generally very sensitive to perturbations due to planets around the lensing objects (Mao & Paczyński 1991; Griest & Safizadeh 1998). Thus, considerable interest was aroused by MOA-11-191/OGLE-11-462 among groups engaged in searches for such planets. As a result, intensive photometric monitoring of this event was carried out by multiple groups, providing valuable data for our analysis.

## 4. HST Observations

The MOA-11-191/OGLE-11-462 event satisfies all the selection criteria for our HST follow-up program described in Section 3.1, and thus we triggered our observing sequence. Our project had a “nondisruptive” target-of-opportunity status, requiring a lead time of about 2–3 weeks from activation to the first observations. The first-epoch HST data were obtained on 2011 August 8, some 19 days after the peak light magnification on 2011 July 20. The magnification was still reasonably high (~12, corresponding to $u \simeq 0.08$), so that the expected astrometric deflection was $\delta \simeq 0.04 \theta_E$ (see Equation (6)), i.e., close to zero at this epoch, but its correct value is taken into account in the model described in Section 8. Subsequent HST observations indicated departure from a linear proper motion for the source. Thus, we continued the imaging, ultimately over an interval of over 6 yr, long enough for robust separation of the relativistic deflection from proper motion. Table 2 gives the HST observing log.

All our HST observations were obtained with the UVIS channel of WFC3, whose CCD detectors provide a plate scale of 39.6 mas pixel⁻¹. To avoid buffer dumps during the orbital visibility period and thus maximize observing efficiency, we used the UVIS2-2K2C-SUB subarray, giving a field of view (FOV) of 80'' × 80''. This FOV is large enough to provide

dozens of nearby astrometric reference stars surrounding the primary target.

The WFC3 detectors are subject to an increasing amount of degradation of their charge transfer efficiency (CTE) as they are exposed to the space environment. The chosen subarray aperture places the target in the middle of the left half of the UVIS2 CCD, which lessens the impact of imperfect CTE relative to a placement closer to the center of the FOV. Nevertheless, a time-dependent correction for CTE must still be applied in the astrometric analysis of the images.

Our HST observations were taken at a total of eight epochs, strategically scheduled for measurement and characterization of the astrometric deflections. At each epoch, we obtained images in two filters (to verify the achromatic nature of the event, and to test for blending by very close companions): “V” (F606W) and “I” (F814W). At the initial epoch, when the source was bright, we obtained nine exposures, four in F606W and five in F814W. At each subsequent epoch, using longer integration times because of the fading of the source, we obtained seven exposures, three in F606W and four in F814W. Individual exposure times were adjusted to take into account the brightness of the source and the orbital visibility of HST and ranged from a minimum of 60 s at the first epoch to a maximum of 285 s at the later epochs. The telescope pointing was dithered by ~200 pixels (~8'') between individual exposures; this allowed retention of a common set of reference stars in all the exposures, in order to mitigate errors in the distortion solution. To maximize the S/N for the most crucial astrometric measurements, we separated Epochs 3 and 4 by only 16 days in 2012 September, around the time when the deflection was expected to be near maximum.

Figure 5 zooms in on the field around the source in Figure 4, showing a 2''1 × 2''0 region as observed at all eight epochs. The bright source is marked with an arrow in the Epoch 1 (top left) image, and it can be seen to fade in the subsequent frames. The astrometric deflection was highest at Epochs 3 and 4, even though the photometric magnification was only about 10% at this epoch. There was very little photometric change in the subsequent epochs, but the astrometric deflections remained detectable until Epoch 7, demonstrating the need for astrometric monitoring over a much longer duration than the photometric variability period.

## 5. HST Data Analysis

### 5.1. Image Processing

We used the flat-fielded and CTE-corrected (f1c) images produced by the Space Telescope Science Institute pipeline

<footnote><sup>a</sup> Individual exposure times ranged from a minimum of 60 s at Epoch 1 to a maximum of 285 s at later epochs.</footnote>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** $2''1 \times 2''0$ cutouts around MOA-11-191/OGLE-11-462 as observed by HST at all eight epochs. Exposures were taken in F606W and F814W from 2011 to 2017; see Table 2 for details. The source star is marked by an arrow in the first-epoch (top left) image. At this epoch, on 2011 August 8, the magnification was a factor of ~12. The maximum astrometric deflection occurred at Epochs 3 and 4, when the photometric magnification was only about 10%. In 2017, the source had returned very close to its unmagnified brightness and undeflected position.</figcaption></figure>

reductions (Sahu et al. 2021; Dressel 2021) for the analysis. As noted above, WFC3 suffered from increasingly poor CTE during this period, so it was essential to take it into account. The _flc_ products were produced using the v2.0 pixel-based CTE model described by Anderson (2021).

## 5.2. Astrometric Analysis

To measure stellar positions in individual frames, we used an updated version of the star-measuring algorithm described in Anderson & King (2006). The routine goes through each exposure pixel by pixel and identifies as a potential star any local maximum that is sufficiently bright and isolated. The routine uses the spatially variable effective point-spread functions (PSFs) provided at the WFC3/UVIS website<fnref n="76" /> to fit the PSF to the star images in the individual _flc_ exposures, in order to determine a position and flux for each star in the raw pixel frame of that exposure. Finally, the positions are corrected for geometric distortion using the distortion solutions provided by Bellini et al. (2011).

As the positions of individual stars are expected to change during the ~6 yr course of our observations owing to their proper motions, we needed to determine their proper motions to properly specify the reference frame. For this, we began with the Gaia Early Data Release 3 (EDR3; Gaia Collaboration et al. 2021) positions and motions for the bright but unsaturated HST stars in the field. The reference frame was constructed to place the bright star close to MOA-11-191/OGLE-11-462 at the center of the reference frame at $(x, y) = (1000, 1000)$ at the 2016.0 epoch, with a plate scale of 40 mas pixel$^{-1}$ and north up. (Note that the Gaia catalog could be incomplete in this region because of the high source density.) Using the Gaia positions and motions, we determined the position for each Gaia star in this frame at each epoch in order to properly transform the distortion-corrected observations at that epoch into the reference frame. This ensures that the proper motions

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Proper motions of the reference stars used in our analysis derived from our HST observations with a 6 yr baseline, versus the Gaia proper motions. Red points represent stars brighter than $G = 18$, for which the Gaia errors are typically <0.2 mas yr$^{-1}$, and blue points represent fainter stars for which the errors are larger. Our analysis uses the Gaia reference frame, so it is natural that there is good agreement between them. But the individual HST measurements have higher precision, particularly at fainter magnitudes.</figcaption></figure>

that we derive represent absolute proper motions. After this initial setup of the reference frame based on the brighter stars, we incorporated high-precision HST stars that were too faint to be found with high precision in the Gaia catalog and solved for their accurate positions and motions. We then used their time-dependent positions to improve the reference frame. Even after allowing individual solutions to improve on the basis of HST observations, there remains very good agreement between our proper motions and those of Gaia. Figure 6 plots the proper motions of our reference stars derived from our HST observations against the Gaia proper motions, where the red points are for brighter stars with $G < 18$ and blue points are for fainter stars with $G \ge 18$. The agreement is imperfect, of course, since the HST observations have a 6 yr baseline and have higher S/N in individual measurements, resulting in higher accuracy in proper-motion measurements, particularly at fainter magnitudes. The agreement is better for the brighter sample, for which Gaia proper-motion errors are typically

<footnote n="76">⁷⁶ https://www.stsci.edu/hst/instrumentation/wfc3/data-analysis/psf</footnote>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Left panel: CMD ($m_{F606W}$ vs. $(m_{F606W} - m_{F814W})$) for all stars in the HST/WFC3 field. The main-sequence turnoff occurs at $m_{F606W} \simeq 22$, above which the disk and the bulge split into two sequences: the redder stars are mainly bulge objects, while the bluer ones are mainly disk objects. Stars marked by blue squares are selected as “disk-like” stars, and the stars marked by black crosses are selected as “bulge-like.” The position of the (unmagnified) microlensed source is shown as a green circle. The nearby bright star 0''4 away and about 20 times brighter than the unmagnified source is shown as a red circle. Magenta circles are the astrometric reference stars used in our analysis. Top right panel: proper motions of the disk-like stars (blue), the reference stars (magenta), and the source (green circle). Bottom right panel: proper motions of the bulge-like stars (black crosses), the reference stars (magenta), and the source (green circle). Only bulge-like stars (shown in magenta) with brightness similar to the observed brightness of the star were used in the final astrometric transformations. This ensures that the reference stars are very similar to the target source star, and uncorrected CTE and small parallax effects should cancel in the differential astrometric measurements.</figcaption></figure>

smaller ($<0.2$ mas yr$^{-1}$). Note that the MOA-11-191/OGLE-11-462 source itself is too faint at baseline for inclusion in the Gaia catalog.

In the next step of our analysis, we used only the HST observations because the Gaia measurements have much higher uncertainties for the fainter stars and the HST observations have a longer baseline of 6 yr compared to the 3 yr of Gaia. In this step of the transformation, we used stars (1) with brightness

similar to the average brightness of the target, (2) with color similar to the source's color (see Figure 7), and (3) lying within 350 pixels of the source. The first criterion minimizes any residual shift caused by CTE effects. We note that we already used the most recent CTE correction software for our analysis. Since the CTE effects on the position measurements are differential, using linear transformations based on stars of similar brightness should remove any residual CTE effects. (It
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Motions of nine representative astrometric reference stars in R.A. and decl. The red lines are linear fits to the proper motions of the stars.</figcaption></figure>

is worth noting here that the images with the highest astrometric deflection were taken when WFC3/UIVS was young and when CTE losses were small.) The second criterion ensures that the stars used in the transformation belong to the bulge, which helps in minimizing errors due to parallax, as described in more detail below. The third criterion minimizes residuals in the distortion solution.

We employed an iterative procedure to measure the positions and proper motions of the stars, starting from the revised values in each iteration. We rejected the highest-sigma point after each complete iteration. We repeated this procedure until the highest-sigma point was no more than a preset tolerance, for which we adopted 6σ. Only a small number of points were

rejected by this procedure, mostly affected by cosmic-ray hits on the detector. Then, at each epoch the reference-star positions were corrected for proper motion, and the positions of the source were determined relative to this adjusted frame. The estimated uncertainty in the position of the source star relative to the adjusted frame is ~0.4 mas in each individual exposure.

As an illustration, Figure 8 shows the proper motions as measured for nine representative stars. Figure 9 shows the errors in the proper-motion measurements of the reference stars. We note that all the reference stars are within about 1.5 mag of each other (See Figure 7). Figure 10 shows the histograms of the residuals of each measurement from that star's proper-motion solution along the R.A. and decl.
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Proper-motion errors of the reference stars. The y-axis shows the proper-motion error, defined as $\sigma_{\mu} = [\sigma_{\mu x}^2 + \sigma_{\mu y}^2]^{0.5}$, where $\sigma_{\mu x}$ and $\sigma_{\mu y}$ are the proper-motion errors along the x- and y-axes, which are parallel to R.A. and decl., respectively. Reference stars cover a range of 1.5 mag around the (unamplified) magnitude of the source; their proper-motion uncertainties vary significantly from object to object, with a modest systematic increase with magnitude.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Histogram of the position residuals for all the reference stars in R.A. and decl. from their proper-motion solution, as measured in each image separately, and scaled with the measured dispersion for that star. The scaled residuals closely follow a standard normal distribution (shown by the red curves). As noted in the caption to Figure 9, reference stars cover a range of magnitudes around the unamplified source magnitude.</figcaption></figure>

directions. Both distributions are consistent with a Gaussian distribution (the red curve). As shown by previous similar studies, the final reference-frame positions are expected to be internally accurate to better than 0.01 pixels (Anderson et al. 2008; Bellini et al. 2015).

We specifically solved for the proper motions of the reference stars, but we ignored their parallaxes. The reason for adopting this approach is the following.

Our choice of reference stars ensures that a large fraction of them belong to the Galactic bulge and hence have similar parallactic motion to the source star. We note that the source parallax is small to begin with ($\lesssim 0.2$ mas). Then, the source position is referenced to stars chosen to be at comparable distance, so any remaining impact of the source parallax on the

<figure><img src="image_14.png" /><figcaption>**Figure 11.** The stacked F606W (left) and F814W (right) PSFs, with the 10-real-pixel radius shown in green. Note that the F606W PSF has a diffraction spike bump very close to the location of the source star. The bump for the other diffraction spikes does not show up in the residual images shown in the second row of Figure 12, which implies that it is subtracted well beneath the source, making the position measurements of the source more robust. The F814W PSF has a very strong radial gradient (along with azimuthal structure) at the location of the source. Without subtracting a high-fidelity model PSF, there would be some impact on the measurement of the source positions. And as the source moves relative to the neighbor, then the source would move across the PSF halo, which could introduce artificial shifts if the PSF of the bright neighbor is not subtracted well.</figcaption></figure>

astrometry or photometry of the event is expected to be negligible.

As described earlier, there is a bright star ~10 pixels away from our source. The target star is close to the brightness of this neighbor in the first two epochs and slowly fades to its nominal brightness, at which point it is about 3 mag fainter than the bright neighbor. For accurate astrometry of the source, we wanted to make sure that the position measurements of the source are not affected by the presence of this bright neighbor. So we wanted to subtract the PSF of the bright star before measuring the positions of the source at every epoch. However, subtracting the bright star is not just a matter of subtracting a standard PSF. The separation is ~10 pixels, and the available library PSFs go out to 12 pixels and are tapered and not very accurate in the wings. Thus, we needed to make a more extended PSF model.

To make an extended PSF, we carefully selected stars that (1) are within ~350 pixels of the bright neighbor, (2) have brightness and color similar to the neighbor, and (3) are fairly isolated. We found 18 such stars (excluding the bright star itself), which provided a good sample to make the required extended PSFs. We used the images of these 18 stars to produce a separate well-sampled, extended PSF for each individual exposure. For illustration, Figure 11 shows the stacked PSFs in F606W and F814W. We have taken particular care to make sure that the PSF is well characterized in the wings since the source lies in the wings of the bright star, and subtracting the wings correctly is crucial for accurate astrometry.

We then took this PSF model for each exposure and subtracted it from the neighbor star in each exposure. Figure 12 shows the original images (first and third rows) and the subtracted images (second and fourth rows) in the F606W and F814W filters. The residuals are very small, particularly in the wings of the PSF. We found that the astrometric position of the source changes by ~0.03 pixels (1.2 mas) after this subtraction, which could have a significant effect on the mass determination of the lens, so this extra step of neighbor subtraction was
<figure><img src="image_15.png" /><figcaption>**Figure 12.** The source star lies in the PSF wings of a bright neighbor, marked in these frames with green circles with a radius of 10 pixels. Since the available WFC3 library PSFs do not extend to this large a radius, a special PSF extending to 20 pixels was constructed using the same HST images. This PSF was used to subtract the neighbor before measuring the position of the source in each exposure. The top row shows the original stacked F606W images of Epochs 1–8 (from left to right), and the second row shows the PSF-subtracted images. The third row shows the original stacked F814W images, and the bottom row shows the subtracted images. The stacks are overbinned by a factor of 2 to show details. Care was taken to assure that the PSF is well characterized in the wings of the neighbor star where the source lies.</figcaption></figure>

crucial in improving the analysis/results. The resulting astrometric positions of the source were used for further analysis as described in the next section.

### 5.3. Photometric Analysis

In addition to the measured positions, the analysis algorithm provides PSF-based photometry of all the stars in the field. To set a calibrated zero-point, we used standard aperture photometry to determine the fluxes of a few isolated stars in the field within an aperture with a 10-pixel radius. These fluxes were then corrected to an infinite aperture, using encircled-energy measurements from Calamida et al. (2021), and the photometric zero-point in the image headers (`PHOTFLAM`) was used to convert these fluxes to the Vegamag scale. The mean difference between these values and the values obtained by the PSF fitting was then applied to all of the PSF magnitudes to convert them to Vegamag.

As described above, the source star lies on the wings of the PSF of the neighboring bright star. Hence, for accurate photometry of the source, it was critical to correctly subtract the light contribution from the bright neighbor. The photometry for the source was carried out after subtracting the superposed flux from the neighbor star using a high-fidelity PSF as described above. The resultant time-series HST photometry of the source is shown in Figure 13, along with the model light curve described below in Section 8. There is no detectable color change as the event progresses and the star fades: the color of the source has remained constant to within 0.01 mag during the 6 yr of observations with HST. There is also no detectable blending as described in more detail in Section 8.2.

<figure><img src="image_16.png" /><figcaption>**Figure 13.** Top panel: photometry of MOA-11-191/OGLE-11-462, obtained with HST over 6 yr at eight epochs in the F606W (V) and F814W (I) filters (black filled circles), along with our final model fits from Section 8.2 (red and blue curves). Errors on the photometry are smaller than the plotting points. The second panel plots the observed values of $V - I$, showing that the source color remained constant to within ~0.01 mag during the entire 6 yr duration. The third and fourth panels show that the residuals, $V_{obs} - V_{model}$ and $I_{obs} - I_{model}$, at the epochs of the HST observations are consistent with zero within the measurement uncertainties (see Section 8.2 for more details and discussion).</figcaption></figure>

(The data used to create this figure are available.)
<table><caption>**Table 3**
Journal of Ground-based Photometry of MOA-11-191/OGLE-11-462</caption><thead><tr><th>Data Set</th><th>Telescope Location</th><th>Aperture (m)</th><th>Filter</th><th>No. of Observations</th><th>Date Range (HJD – 2,450,000)</th></tr></thead><tbody><tr><td>MOA</td><td>New Zealand</td><td>1.8</td><td>R</td><td>46040</td><td>3824.093...4441.117</td></tr><tr><td>OGLE</td><td>Chile</td><td>1.3</td><td>I</td><td>15546</td><td>5260.855...8787.509</td></tr><tr><td>Wise Survey</td><td>Israel</td><td>1.0</td><td>I</td><td>953</td><td>5658.529...5722.543</td></tr><tr><td>Danish DFOSC</td><td>Chile</td><td>1.54</td><td>I</td><td>921</td><td>5744.804...5782.554</td></tr><tr><td>Danish LuckyCam</td><td>Chile</td><td>1.54</td><td>broad</td><td>10</td><td>5738.624...5742.706</td></tr><tr><td>MONET North</td><td>Texas, USA</td><td>1.2</td><td>I</td><td>214</td><td>5762.722...5764.824</td></tr><tr><td>Faulkes North</td><td>Hawaii, USA</td><td>2.0</td><td>SDSS i'</td><td>99</td><td>5763.777...5768.955</td></tr><tr><td>Liverpool</td><td>Canary Islands, Spain</td><td>2.0</td><td>SDSS i'</td><td>254</td><td>5739.521...5768.438</td></tr><tr><td>SAAO 1.0 m</td><td>South Africa</td><td>1.0</td><td>I</td><td>611</td><td>5751.264...5777.315</td></tr><tr><td>SAAO 1.0 m</td><td>South Africa</td><td>1.0</td><td>V</td><td>44</td><td>5758.454...5765.385</td></tr><tr><td>U. Tasmania</td><td>Australia</td><td>1.0</td><td>I</td><td>60</td><td>5761.055...5769.110</td></tr><tr><td>CTIO</td><td>Chile</td><td>1.3</td><td>I</td><td>226</td><td>5757.511...5772.785</td></tr><tr><td>CTIO</td><td>Chile</td><td>1.3</td><td>V</td><td>27</td><td>5757.515...5763.807</td></tr><tr><td>Auckland</td><td>New Zealand</td><td>0.4</td><td>R</td><td>160</td><td>5759.873...5778.893</td></tr><tr><td>Farm Cove</td><td>New Zealand</td><td>0.35</td><td>unfiltered</td><td>37</td><td>5741.856...5761.856</td></tr><tr><td>Kumeu Obs.</td><td>New Zealand</td><td>0.35</td><td>R</td><td>63</td><td>5759.803...5762.134</td></tr><tr><td>Vintage Lane</td><td>New Zealand</td><td>0.4</td><td>unfiltered</td><td>60</td><td>5762.792...5767.941</td></tr><tr><td>Weizmann</td><td>Israel</td><td>0.4</td><td>I</td><td>167</td><td>5762.281...5764.423</td></tr><tr><td>Wise</td><td>Israel</td><td>0.46</td><td>I</td><td>142</td><td>5762.266...5763.456</td></tr></tbody></table>

Note. The data were acquired by the MOA (Bond et al. 2001), OGLE (Udalski et al. 2015), Wise Microlensing Survey (Shvartzvald et al. 2016), MiNDSTEp (Dominik et al. 2010), RoboNet (Tsapras et al. 2009), PLANET (Sackett et al. 2004), and μFUN (Gould et al. 2006) teams.

The source at baseline brightness has apparent magnitudes of $m_{F606W} = 21.943 \pm 0.014$ and $m_{F814W} = 19.578 \pm 0.012$.

We made a stack of all the images for each filter for every epoch. The HST images allow us to detect and measure magnitudes of stars as faint as $V \sim 25$. Figure 7 shows the CMD based on this photometry, where we also show the position of the source and the bright neighbor $0''4$ away.

## 6. Ground-based Light Curve

MOA-11-191/OGLE-11-462 was monitored photometrically by several ground-based observatories. The coverage by MOA, OGLE, and Wise Microlensing Survey extended over several years. Moreover, as a high-magnification event, it attracted intensive monitoring by a number of additional ground-based telescopes—especially around the time of peak brightness, where the microlensing light curve is sensitive to planet detection. Table 3 gives a journal of the photometric observations and data used in our analysis. We use the data reduced by the surveys and other groups (Udalski et al. 1992; Bond et al. 2001; Sackett et al. 2004; Gould et al. 2006; Tsapras et al. 2009; Dominik et al. 2010).

Figure 14 shows the light curve, both over a 300-day interval (top panel) and zooming in on the 7 days around peak magnification (bottom panel). Superposed is our model fit to the light curve, from the analysis described in the next two sections.

## 7. Blending, Relative Parallax, and Lens Trajectory

The combination of ground-based photometric monitoring with long-term astrometric and photometric measurements from HST affords the ability to constrain all aspects of this event and obtain high-quality measurements of its parameters. In this section, we describe some characteristics of the data, illustrating the key information that can be obtained from photometry and astrometry separately through heuristic considerations.

### 7.1. Photometric Blending

Ground-based photometry of MOA-11-191/OGLE-11-462 suffers from significant amounts of blending with neighboring stars, as shown in our HST images (Figure 4). Moreover, the amount of blending changes markedly depending on image quality. There is a bright neighbor star only $0''4$ away, along with two fainter stars within $0''5$ of the source. At larger separations, there are three more stars within $1''$, whose combined brightness is greater than that of the baseline source, and there are several more stars within $1''5$ that are also brighter than the source. Since in the available ground-based imaging data the measured image quality is seldom better than $1''$, the ground-based photometry will always include the light from at least the three closest stars within $0''5$. In a fraction of the data, taken under poorer seeing conditions, the source is blended with an increasing number of neighbors.

Our HST images allow us to place constraints on the expected blending parameter, *g*, defined as the ratio of the flux from neighbors included in the photometry to the flux from the unmagnified source itself (Equation (5)), in the ground-based observations. The bright neighbor is 18.88 times brighter than the source at baseline in F814W, and the contribution from the two fainter stars is an additional 0.19 times that of the source. Thus, the expected blending factor due to these three stars is $g = 19.07$ in the F814W filter. Since the bright neighbor is similar in color to the source (see Figure 7), we adopt this value of *g* for both OGLE (I band) and MOA (R band) data as an initial estimate, but we keep it as a variable in our analysis. The final values (Section 5.2) differ significantly between OGLE and MOA, possibly due in part to differences in processing between the two data sets.

The effect of blending can be reduced substantially by basing the photometry on difference images. It can be further reduced by restricting the analysis to images taken under good seeing. It is obvious, however, that variable blending will affect the noise characteristics of ground-based photometry; it is
<figure><img src="image_17.png" /><figcaption>**Figure 14.** Ground-based photometric observations of MOA-11-191/OGLE-11-462, along with a best-fitting model light curve, shown over a 300-day interval in the top panel and over a zoomed-in region covering 7 days around peak magnification in the bottom panel. All data have been transformed to OGLE *I* magnitudes according to the inferred baseline magnitudes and blend ratio from the common model.</figcaption></figure>

difficult to include such effects deterministically because of the imperfect knowledge of the blending (at the subpercent level) for individual images.

## 7.2. Heuristic Considerations

### 7.2.1. Photometric Constraints on Parallax and Lens Trajectory

As discussed in Section 2.3, the light curve of a long-duration microlensing event such as MOA-11-191/OGLE-11-462 can show distortion by the relative parallactic motions of the source and lens (e.g., Gould 1992; Alcock et al. 1995). Specifically, the light curve is sensitive to $\pi_E$ and $\phi_{LS}$ because their combination modifies the relative path of source and lens and thus the shape of the light curve. We note that, in our formalism, $\phi_{LS}$ corresponds to the PA of the path of the lens relative to the source without parallax in equatorial coordinates (not to be confused with the instantaneous path of the lens at the time of closest angular approach in ecliptic coordinates; see Section 2.2).

In principle, a sufficiently accurate light curve can provide good constraints on both $\pi_E$ and $\phi_{LS}$. However, as discussed in the previous subsection, the photometry is significantly affected by blending. We attempted to model the light curve alone but

found that it can be fitted with a range of parameter combinations, in which the values of $\pi_E$ and $\phi_{LS}$ are strongly correlated. In addition, the derived value of $\pi_E$ varies with the specific subset of photometric data chosen for analysis, as well as with the assumed blending factor for those data. The derived value of $\pi_E$ ranges from 0.07 to 0.12, with larger values corresponding to larger values of $\phi_{LS}$, ranging from 330° to 358°. The reason is that increasing $\phi_{LS}$ makes the lens move in a more northerly direction as seen in the bottom panel of Figure 15. Since parallax is predominantly in the east-west direction, in order to produce a fixed change in $u$, the value of $\pi_E$ has to increase with $\phi_{LS}$, so that the change in position due to parallax can compensate for a more northerly motion of the lens. Several different combinations of these quantities can reproduce the observed light curve, with differences between solutions of the order of 1 mmag at early and late times and ~5 mmag near the peak. Systematic differences in the data at this level could be caused by small variations in blending associated with changes in the ground-based seeing, or other minor secular variations in the photometry. Therefore, we conclude that, when photometry alone is used to constrain the parameters of the event, only a reliable joint constraint on $\pi_E$ and $\phi_{LS}$ can be derived. Fortunately, astrometry provides a
<figure><img src="image_18.png" /></figure>

**Figure 15.** Top panel: path of the lens with respect to the source, the position of the source being fixed at (0, 0). The direction of motion is shown by an arrow. The black straight line is the lens path without parallax, and the wavy line is the path including parallax. Red points mark the epochs of HST observations. $u_{RA}$ and $u_{Dec}$ are the R.A. and decl. components of the lens-source separation, in units of $\theta_E$. Bottom panel: enlarged view of the lens trajectory with respect to the source around closest angular approach. Black lines correspond to the $u_{0,+}$ solution, which is the preferred solution, and the red lines show the less preferred $u_{0,+}$ solution. In both cases, the barycentric trajectories (no parallax), shown by dotted lines, pass on the north side of the source. The geocentric trajectory (i.e., the trajectory as seen by an observer on Earth, which includes parallax) of the $u_{0,+}$ solution (shown by the solid red line) passes on the north side of the source, but the trajectory of the $u_{0,-}$ solution (shown by the black solid line) passes to the south. The $u_{0,-}$ solution is used in our analysis presented here. It is worth noting, however, that the derived $\theta_E$, and hence the mass of the lens, is nearly identical for both solutions. $\varphi_{LS}$ is the PA of the barycentric proper motion of the lens with respect to the source (i.e., no parallax), as shown here.

robust independent estimate of $\varphi_{LS}$, allowing us to break this degeneracy and determine the two quantities separately.

### 7.2.2. Astrometric Deflection and Orientation of the Relative Motion

In order to understand how astrometry can constrain the direction of motion of the lens, it is useful to consider an illustrative plot of the motion of the lens relative to the source, as shown in Figure 15. North is at the top, east on the left; the ($u_{RA}$, $u_{Dec}$) coordinates give the position of the lens relative to the source in units of $\theta_E$, with $u_{RA}$ increasing to the east. We have used our actual final model described below in Section 8.2 for this illustration. The top panel shows the motion of the lens, with $\pi_E = 0.0894$ and $\varphi_{LS} = 342^\circ.5$, and an impact parameter of $u_0 = 0.00271$. The straight line represents the proper motion of the lens with respect to the source, while the wavy line adds the parallactic motion, computed using the JPL ephemeris of Earth.<fnref n="77" /> Red circles show the position of the lens at the eight epochs of our HST observations.

The bottom panel in Figure 15 shows an enlarged view of the lens trajectory near the source position. The dotted black line represents the proper motion of the lens with respect to the source without parallax, while the solid black line includes the parallactic effect. (The red lines correspond to a less preferred $u_{0,+}$ solution described in Section 8.3.) The plot shows that near the closest angular approach—and thus the peak magnification—the relative path is substantially affected by parallax; however, the astrometric deflection is very small at this time (see Section 2.3). Since the source deflection is always in the direction of the line joining the instantaneous position of the lens to the undeflected position of the source, the directions of the source deflections at late times will remain nearly constant with little parallax effect; thus, the late-time deflection directions robustly constrain the orientation, $\varphi_{LS}$, of the lens trajectory.

### 7.2.3. Constraining the Lens Trajectory Orientation

Since the parallactic effect is unimportant for constraining $\varphi_{LS}$, we first fitted the photometry using a light-curve model that neglects parallax. The resultant model (with $t_0 = 2,450,000 = 5763.33$, $t_E = 231.56$ days, $A_{max} = 372.62$, and $g = 19.3$) predicts the total deflection in units of $\theta_E$ as a function of time, through Equation (6).

As described in Section 5.2, we have accurate measurements of the (x, y) positions of the source at the eight epochs of HST observations. These positions are affected by both the proper motion of the source and its deflections, which are a function of $\theta_E$ and $\varphi_{LS}$. We fitted a model to the positions, whose parameters are the x (R.A.) and y (decl.) components of the proper motion, $\theta_E$, and $\varphi_{LS}$. This fit resulted in values of $\theta_E = 5.2 \pm 0.5$ mas and $\varphi_{LS} = 337^\circ.9 \pm 5^\circ.0$.

Components of the resultant deflections as a function of time, after subtracting the best-fitting proper motion, are shown in the top two panels of Figure 16. These panels plot the deflections in the R.A. and decl. directions. The bottom left panel shows the total amount of deflection, again as a function of time. Note that the total deflection reaches a maximum of ~2 mas in late 2012.

The bottom right panel of Figure 16 plots the R.A. versus decl. deflections. These deflections are always along the line

<footnote n="77">⁷⁷ https://ssd.jpl.nasa.gov/horizons/app.html#/</footnote>
<figure><img src="image_19.png" /></figure>

**Figure 16.** The top two panels show the average values of the measured deflections in R.A. and decl. at each HST epoch. We use the values of $t_0$, $t_E$, and $u_0$ as derived from a light-curve fit without parallax, and fit for the proper motion of the source, $\theta_E$, and $\phi_{LS}$. The solid black line is the best fit with $\theta_E = 5.2$ mas and $\phi_{LS} = 337^{\circ}9$. The bottom left panel shows the total deflections at each epoch. The bottom right panel shows the R.A. vs. decl. deflections. These deflections are always along the line joining the lens to the source, and thus at large separations their direction is opposite to the direction of the relative motion of the lens with respect to the source. The solid black line passes through the origin, at the $\phi_{LS}$ angle derived above. The dotted blue lines indicate the allowed range from astrometry based on the uncertainties. The allowed range of $\phi_{LS}$ from photometry alone is shown by the dotted magenta lines, the best-fit value being $\phi_{LS} = 354^{\circ}8$ shown by the solid magenta line. The combined constraint from photometry and astrometry is used in our subsequent analysis.

joining the lens to the source, and thus at large separations their direction is opposite to the direction of the relative motion of the lens. The solid black line passes through the origin, at the $\phi_{LS}$ angle derived above. The dotted blue lines indicate the allowed range based on the uncertainties. Although the best-fit value from photometry alone is $\phi_{LS} = 354^{\circ}8$ (shown by the solid magenta line), the allowed range of $\phi_{LS}$ (shown by the dashed magenta lines) has a flat probability distribution. The combined constraint from photometry and astrometry is used in our subsequent analysis.

Note that for clarity in Figure 16 we have shown the mean deflections at each epoch, and not the individual measurements. In particular, this avoids a confusing overlap of points in the bottom right panel, where the deflections are not a monotonic function of time. The individual measurements are shown in the next section.

# 8. Full Modeling of the Photometric and Astrometric Data

In this section we give full details of our analysis, carried out independently by several coauthors using different parameterizations, all leading to a consistent final model of the event.

## 8.1. First Approach: All Photometric Data Sets, and Robustness of Parallax Measurement

In addition to OGLE and MOA, photometric time-series data were obtained at several more observatories (see Table 2). These data typically cover a relatively narrow range of ~20 days around the peak, being primarily aimed at searching for planet-related distortions of the light curve. These photometric series, unlike those obtained by the survey programs, do not provide significant constraints on the lens-source model—especially since each set of observations can have a different baseline magnitude and blending parameter. Nevertheless, in our initial approach, we included all the data in our analysis but also carried out analyses separately for “OGLE-only,” “MOA-only,” and “OGLE+MOA-only” data sets.

Ground-based time-series photometry is susceptible to systematic noise, and this must not be mistaken for real features of the light curve. To improve the robustness of our solution, we model the photometric uncertainties, and moreover force the model to follow the bulk of the data by explicitly down-weighting outliers. As implemented for the SIGNAL-MEN microlensing anomaly detector (Dominik et al. 2007, 2019), we specifically adopt a bi-square weight function
with regard to the median residual and a Gaussian distribution for the uncertainties with revised standard deviation in magnitude of

$$ \tilde{\sigma} = \sqrt{(\kappa \sigma)^2 + \sigma_0^2}, \quad (17) $$

where $\sigma$ denotes the reported error bar, $\kappa$ is a scaling factor, and $\sigma_0$ corresponds to a systematic error added in quadrature. For the plot of the various data sets as shown in Figure 14, we give the respective estimated values of $\kappa$ and $\sigma_0$ in Table 4. If the size of the error bars does not vary substantially, there is a degeneracy between $\kappa$ and $\sigma_0$, and either of the parameters provides modified constant error bars (while it does not matter which).

We find two viable models, significantly only distinguished by the sign of $u_0$, and will refer in the following to the model with $u_0, \dots$. The microlensing parallax parameter $\pi_E$ is constrained by the wing of the light curve and much less sensitive to the peak region, and therefore it is mostly constrained by the microlensing survey data. From various combinations of data sets, we consistently find $\pi_E = 0.10 \pm 0.02$ but some variation in the trajectory angle $\psi$, correlated with $t_E$ and the blend fraction, yielding visually indistinguishable model light curves.

However, the angle of lens-source proper motion, $\phi_{LS}$, follows robustly from the astrometric data (see Section 7.2.2), given that the centroid shift to first order (i.e., neglecting the small distortion caused by parallax) traces an ellipse (a highly flattened ellipse resembling a line in our case) whose semimajor axis is parallel to $\mu_{LS}$. If we restrict this angle to the range $333^\circ \le \phi_{LS} \le 343^\circ$, as suggested by the astrometric data, the photometric light curve does not change substantially, and we find $\pi_E \simeq 0.086$.

We emphasize here that it is incorrect to say that there is a discrepancy between the paths determined from photometry and astrometry, since there is a correlation between $\phi_{LS}$ and $\pi_E$ in the photometric solution. However, restricting the trajectory angle $\phi_{LS}$ as robustly derived in the last section from the orientation of the centroid shifts to the range $333^\circ \le \phi_{LS} \le 343^\circ$ suggests $\pi_E \simeq 0.086$.

## 8.2. Second Approach: Simultaneous Fit of Photometric and Astrometric Data

We now turn to a full analysis in which we fit the astrometric and photometric data simultaneously in order to obtain all of the parameters. Such a solution is important, since the crucial parameters of $\theta_E$ and $\pi_E$ are derived from two different types of data. A simultaneous solution is also essential for a correct estimate of the uncertainties in the model parameters.

We follow the same plane-of-sky approach described in Section 7.2, which makes it easier to work with, and also show the actual paths of the lens and the source, as well as the deflections. We follow a different parameterization procedure where the model parameters we optimize contain all terms needed to characterize the positions of the lens and the source on the sky as a function of time; these include the reference positions and proper motions of both lens and source, their relative parallax, and the angular Einstein radius of the lens. In principle, the source parallax is also needed; however, its parallax in the reference system we use is close to zero, and it is not meaningfully constrained by the observations. As discussed below in Section 9.2, the best constraints on the source distance come instead from photometry and high-resolution spectroscopy.

<table><caption>**Table 4**
Revised Error Bars of Photometric Data</caption><thead><tr><th>Data Set</th><th>κ</th><th>σ<sub>0</sub></th></tr></thead><tbody><tr><td>OGLE</td><td>1.25</td><td>0.005</td></tr><tr><td>MOA</td><td>1.18</td><td>0.0025</td></tr><tr><td>Wise Survey</td><td>0.64</td><td>0.009</td></tr><tr><td>Danish 1.54 m DFOSC</td><td>2.75</td><td>10<sup>-5</sup> (*)</td></tr><tr><td>Danish 1.54 m LuckyCam</td><td>0.1 (*)</td><td>0.050</td></tr><tr><td>MONET North 1.2 m</td><td>0.57</td><td>0.006</td></tr><tr><td>Faulkes North 2.0 m</td><td>0.1 (*)</td><td>0.014</td></tr><tr><td>Liverpool 2.0 m</td><td>1.17</td><td>0.008</td></tr><tr><td>SAAO 1.0 m I</td><td>1.42</td><td>0.0014</td></tr><tr><td>SAAO 1.0 m V</td><td>0.59</td><td>0.009</td></tr><tr><td>Tasmania 1.0 m</td><td>0.28</td><td>0.017</td></tr><tr><td>CTIO 1.3 m I</td><td>0.62</td><td>10<sup>-5</sup> (*)</td></tr><tr><td>CTIO 1.3 m V</td><td>1.01</td><td>10<sup>-5</sup> (*)</td></tr><tr><td>Auckland 0.4 m</td><td>0.91</td><td>0.007</td></tr><tr><td>Farm Cove 0.35 m</td><td>1.00</td><td>0.003</td></tr><tr><td>Kumeu 0.35 m</td><td>1.87</td><td>10<sup>-5</sup> (*)</td></tr><tr><td>Vintage Lane 0.4 m</td><td>0.1 (*)</td><td>0.011</td></tr><tr><td>Weizmann 0.4 m</td><td>0.80</td><td>10<sup>-5</sup> (*)</td></tr><tr><td>Wise 0.46 m</td><td>0.22</td><td>0.009</td></tr></tbody></table>

**Note.** The adopted error bar (in magnitude) becomes $\tilde{\sigma} = \sqrt{(\kappa \sigma)^2 + \sigma_0^2}$, where $\sigma$ is the reported error bar of the photometric data. We have applied range constraints $\kappa > 0.1$ and $\sigma_0 \ge 10^{-5}$. An asterisk (*) indicates that the value is at the range boundary.

From these parameters, the undeflected paths of the source and lens can be determined. The deflection of the lensed image of the source is then computed, and the resulting deflected source positions are matched to the observed positions. The same calculation also yields the source magnification; in order to match the observed photometry, the model must include a baseline magnitude of the source and a blending parameter for each photometric data set. Consistent with the previous approach, we found that most of the photometric data sets cover too short a time interval to yield meaningful constraints on the event parameters in the presence of significant blending; therefore, we limit the model optimization to the MOA and OGLE photometric data sets and validate the resulting model for the other data sets separately (see Section 8.1). In addition, in order to avoid undue impact from any secular variations in photometric responses, we only include OGLE and MOA photometric measurements within $\pm 2$ yr from the peak of the event. We adopt the approximate values of $\phi_{LS}$, $t_E$, $t_0$, and $u_0$ from the analysis of the previous section as our initial estimates, but we leave all parameters free in the optimization; the baseline magnitude and blending parameters for MOA and OGLE are also separately optimized.

As discussed in the previous subsection, the results of the optimization depend to some extent on the relative weighting of astrometry and photometry. Because the number of photometric measurements greatly exceeds that of astrometric measurements and the nominal photometric uncertainties are very small, an optimization using nominal errors disproportionately weights photometry, resulting in a poor match to the astrometry. In order to obtain a more balanced weighting of astrometry and photometry, we scaled the photometric errors by different amounts in different temporal bins making sure that the scaled errors are compatible with the statistical dispersions in the measurements, and validated the solution based on a reasonable match to the astrometric data. In our final
<figure><img src="image_20.png" /><figcaption>**Figure 17.** Representation of the reconstructed motions of the lens and of the source on the plane of the sky. The reference point is the undeflected position of the source at time $t_0$. The small circles are the individual measurements from HST images (green for F606W, orange for F814W); the larger red triangles represent the average positions for each epoch. The cyan squares are the fitted positions at each epoch. The gray lines connect the undeflected source and lens positions at each HST epoch.</figcaption></figure>

model, most photometric points are still at an uncertainty below 10 mmag. This solution has a total astrometric $\chi^2$ of 136 with 106 points. (The solution with nominal weights has a higher astrometric $\chi^2$ of 149.) The parameters of the final model are given in Table 5.

Figure 17 shows the reconstructed motion of the lens (magenta) and of the source (black) in the plane of the sky, based on our adopted model. The predicted apparent source trajectory is shown by the green solid line. The astrometric measurements are shown individually by the small filled circles, and their epoch averages are shown as red triangles. Cyan squares show the model position at each epoch; gray lines connect the model lens position to the undeflected and deflected source positions at each epoch.

Figure 18 presents the measured and predicted source positions separately for the x and y coordinates. To improve the legibility of the plot, the fitted proper motion of the source has been subtracted from both model and data; therefore, the points shown represent the deflection of the source. The black line is our final adopted model, which takes photometry and astrometry into account.

We can now check the consistency of the final model with our measured HST photometry in the F606W and F814W filters and also constrain the amount of blending. We used the final model to calculate the magnifications at the HST observation epochs and varied the baseline magnitudes and the blending factor to

<table><caption>**Table 5**
Parameters of the Full Fit to Astrometry and Photometry</caption><thead><tr><th>Parameter</th><th>Units</th><th>Value</th><th>Uncertainty (1σ)</th><th>Notes<sup>a</sup></th></tr></thead><tbody><tr><td>μ<sub>S</sub> (R.A.)</td><td>mas yr<sup>-1</sup></td><td>-2.263</td><td>0.029</td><td>(1)</td></tr><tr><td>μ<sub>S</sub> (decl.)</td><td>mas yr<sup>-1</sup></td><td>-3.597</td><td>0.030</td><td>(1)</td></tr><tr><td>θ<sub>E</sub></td><td>mas</td><td>5.18</td><td>0.51</td><td>(2)</td></tr><tr><td>t<sub>E</sub><sup>*</sup></td><td>days</td><td>270.7</td><td>11.2</td><td>(3)</td></tr><tr><td>φ<sub>LS</sub></td><td>deg</td><td>342.5</td><td>4.9</td><td>(4)</td></tr><tr><td>t<sub>0</sub><sup>*</sup> (HJD - 2,450,000.0)</td><td>days</td><td>5765.00</td><td>0.87</td><td>(5)</td></tr><tr><td>π<sub>E</sub></td><td></td><td>0.0894</td><td>0.0135</td><td>(6)</td></tr><tr><td>u<sub>0</sub><sup>*</sup></td><td></td><td>0.0422</td><td>0.0072</td><td>(7)</td></tr><tr><td>MOA baseline <i>R</i><br>magnitude</td><td>mag</td><td>16.5147</td><td>0.0016</td><td></td></tr><tr><td>MOA blending parameter</td><td></td><td>16.07</td><td>0.66</td><td></td></tr><tr><td>OGLE baseline <i>I</i><br>magnitude</td><td>mag</td><td>16.4063</td><td>0.0015</td><td></td></tr><tr><td>OGLE blending parameter</td><td></td><td>18.80</td><td>0.79</td><td></td></tr><tr><td>Derived parameters:</td><td></td><td></td><td></td><td></td></tr><tr><td>u<sub>0</sub></td><td></td><td>0.00271</td><td></td><td>(8)</td></tr><tr><td>t<sub>0</sub> (HJD - 2,450,000.0)</td><td>days</td><td>5763.32</td><td></td><td>(9)</td></tr></tbody></table>

**Note.**

<footnote><sup>a</sup> Notes: (1) Undeflected proper motion of the source in Gaia EDR3 absolute frame. (2) Angular Einstein radius. (3) Angular Einstein radius $\theta_E$ divided by absolute value of lens-source proper motion $\mu_{LS}$. (4) Orientation angle of the lens proper motion relative to the source (N through E). (5) Time of closest angular approach without parallax motion. (6) Relative parallax of lens and source in units of $\theta_E$. (7) Impact parameter in units of $\theta_E$ without parallax motion. (8) Impact parameter derived using the model-fit parameters above and after including parallax motion, in units of $\theta_E$. (9) Time of closest angular approach derived as above and after including parallax motion.</footnote>
<figure><img src="image_21.png" /><figcaption>**Figure 18.** Predicted and measured positions of the source for our adopted joint astrometric and photometric fit. Small circles show the individual measurements from HST images (green for F606W, orange for F814W), while red triangles with error bars show the average and uncertainty at each HST epoch. The black line is our final adopted model. The fitted proper motion for the source has been subtracted from both model and measurements, in order to allow a better scaling of the plot.
(The data used to create this figure are available.)</figcaption></figure>

fit the HST photometry. The resulting fit is excellent (see Figure 13) and yields baseline source magnitudes of $m_{F606W} = 21.943 \pm 0.012$ and $m_{F814W} = 19.578 \pm 0.012$, with corresponding blending factors for the HST photometry of $g = -0.012 \pm 0.015$ and $-0.006 \pm 0.012$, respectively. The bottom two panels of Figure 13 show the residuals of the observations relative to the model in F606W and F814W at the epochs of HST observations, where we have assumed the minimum physically allowed value of $g = 0$. (It makes little difference if we instead use the slightly negative values of $g$ from the formal fit.) The stringent constraints on blending and lack of color variation make it unlikely that our deflection measurements could be affected by blending with a binary companion or a field star lying within the HST PSF.

## 8.3. Lens Motion

The path of the lens in the sky plane, as derived from the above analysis, is shown in Figure 15. We note that there are two solutions corresponding to $u_{0,+}$ and $u_{0,-}$ (Gould 2004). However, at $t_0$, the angular separation between the source and the lens is dominated by the parallactic motion of the lens (the separation caused by the parallactic motion is ~0.04, compared to the actual value of $u_0 \simeq 0.00271$). Thus, the path without parallax in both solutions lies on the $u_{0,+}$ side.

The paths of the lens for the $u_{0,+}$ and $u_{0,-}$ solutions are shown in Figure 15, where the dashed lines show the path of the lens without parallax, and the solid lines show the path of the lens after the parallactic motion is taken into account. The respective paths are separated only by 0.02 mas at the time of

maximum magnification, and they quickly merge. The $u_{0,-}$ solution is the preferred solution and is the one used in our analysis. However, we verified that both solutions provide practically identical deflections (since the deflection measurements are at much higher $u$, where the two paths nearly merge), and the results are the same for all practical purposes in both solutions.

The first two rows of Table 5 give the proper motions of the source. We use the values of $\theta_E, v_E^*$, and $\phi_{LS}$ given in the next three rows to determine the proper motion of the lens with respect to the source as $-2.10 \pm 0.22$ mas yr$^{-1}$ and $6.66 \pm 0.67$ mas yr$^{-1}$ in R.A. and decl., respectively. The resulting absolute proper motions of the lens are given below in Table 6.

## 9. Properties of the Source

As described in Sections 2.2–2.3, the mass determination for the lens does not depend on the individual distances to the lens and source, but only on the relative lens–source parallax, $\pi_{LS}$, and the Einstein ring radius, $\theta_E$—quantities that are directly determined from the light curve and the measured astrometric deflections. However, we still need an estimate of the distance to the source, $D_S$, in order to determine the distance to the lens, $D_L$, which is discussed below. Moreover, MOA-11-191/ OGLE-11-462 was a very high magnification event, with an impact parameter of only $u_0 \simeq 0.00271$. Thus, it is desirable to estimate the angular diameter of the source, to verify that it is consistent with the point-source light-curve modeling adopted in the previous section.
<table><caption>**Table 6**
Properties of the MOA-11-191/OGLE-11-462 Black Hole Lens</caption><thead><tr><th>Property</th><th>Value</th><th>Sources and Notes<sup>a</sup></th></tr></thead><tbody><tr><td>Mass, M<sub>lens</sub></td><td>7.1 ± 1.3 M<sub>c</sub></td><td>(1)</td></tr><tr><td>Distance, D<sub>L</sub></td><td>1.58 ± 0.18 kpc</td><td>(2)</td></tr><tr><td>Einstein ring radius, θ<sub>E</sub></td><td>5.18 ± 0.51 mas</td><td>(3)</td></tr><tr><td>Proper motion, (μ<sub>α</sub>, μ<sub>δ</sub>)</td><td>(−4.36 ± 0.22, +3.06 ± 0.66) mas yr<sup>−1</sup></td><td>(4)</td></tr><tr><td>Galactic position, (X, Y, Z)</td><td>(−4, −1580, −45) pc</td><td>(5)</td></tr><tr><td>Space velocities, (V, W)</td><td>(+3, +40) km s<sup>−1</sup></td><td>(6)</td></tr></tbody></table>

**Note.**

a Sources and notes: (1) This paper, Section 10.1. (2) This paper, Section 10.2. (3) This paper, Table 5. (4) This paper, derived from model in Table 5; absolute proper-motion components in Gaia EDR3 J2000 frame. (5) Galactic position relative to Sun: X in direction of Galactic rotation, Y in direction away from Galactic center, Z perpendicular to Galactic plane toward north Galactic pole. (6) Space velocity components relative to Sun, assuming zero radial velocity: V in direction of Galactic rotation, W toward north Galactic pole; U component toward Galactic center is undetermined owing to unknown actual radial velocity.

In this section we use results of ground-based spectroscopy of the magnified source and HST photometry at its baseline brightness, to estimate its distance and angular diameter.

## 9.1. High-resolution Spectroscopy

In a target-of-opportunity program focusing on high-dispersion spectroscopy of high-amplitude microlensing events in the Galactic bulge, Bensby et al. (2013) obtained four spectra of MOA-11-191/OGLE-11-462 on 2011 July 20 and 21. These observations were made almost precisely around the dates of maximum magnification. The authors employed spectrographs on three different large telescopes, VLT, Magellan, and Keck I, at spectral resolutions ranging from 46,000 to 90,000. The aim of their program was to determine chemical compositions of dwarfs and subgiant stars in the Galactic bulge.

Bensby et al. (2013) carried out a non-LTE model atmosphere analysis of the spectra, obtaining an effective temperature and surface gravity of $T_{eff} = 5382 \pm 92$ K and $\log g = 3.80 \pm 0.13$, consistent with a late G-type subgiant. The metallicity was found to be slightly above solar, at $[Fe/H] = +0.26 \pm 0.14$, and the radial velocity is $+134.0$ km s$^{-1}$.

As discussed above (e.g., Figures 4 and 5), the HST images show that MOA-11-191/OGLE-11-462 is accompanied by a neighboring star only 0.4 away, which would have been included in the spectrograph apertures. However, at the dates of the observations, MOA-11-191/OGLE-11-462 was ~3.2 mag brighter than the neighbor in the F606W bandpass; thus, the contamination was about 5%, which should have a minimal effect on the spectroscopic investigation.

As shown in Figure 7, the neighbor has a color similar to that of the baseline source but is more luminous. We verified its small effect on the spectroscopy by artificially contaminating the source spectrum with a 5% flux from a giant star with the same effective temperature and metallicity, but with $\log g = 0$. This produced a smaller change in the line strengths of the combined light than would a change of $\log g$ for the source star by its 0.13 dex uncertainty.

## 9.2. Distance and Angular Diameter

We estimate the distance to the source using two independent sets of theoretical stellar isochrones: “Padova” isochrones, which we obtained using the PARSEC web tool<sup>78</sup> (Bressan et al. 2012; Marigo et al. 2017), and “BaSTI” isochrones<sup>79</sup>

<figure><img src="image_22.png" /><figcaption>**Figure 19.** Location of the MOA-11-191/OGLE-11-462 source star in the $\log g$ vs. effective temperature diagram (red filled circle with error bars), using stellar parameters from Bensby et al. (2013). Superposed are theoretical isochrones from two sources. The blue curves show PARSEC isochrones, calculated for the source’s measured metallicity of $[M/H] = +0.26$, for ages of 2–10 Gyr. The orange curves are BaSTI isochrones for $[M/H] = +0.3$ and the same range of ages.</figcaption></figure>

(Hidalgo et al. 2018; Pietrinferni et al. 2021). In Figure 19 we plot the position of the source in the distance-independent $\log g$ versus $\log T_{eff}$ plane (red filled circle with error bars). Superposed are PARSEC isochrones, calculated for the measured metallicity of $[M/H] = +0.26$ (blue curves), and BaSTI isochrones, for a metallicity of $[M/H] = +0.3$ (orange curves). Both sets of isochrones agree well around the location of the source star, yielding an age of about 4 Gyr, but with considerable uncertainty.

At the position of the source in the diagram, a star on the PARSEC 4 Gyr isochrone has absolute magnitudes (Vegamag scale) in the native HST/WFC3 bandpasses of $M_{F606W} = +3.04$ and $M_{F814W} = +2.40$ and an intrinsic color of $(m_{F606W} - m_{F814W})_0 = 0.64$. The BaSTI models—which have been updated recently, using new filter throughput curves released in 2020 October by the WFC3 team (Calamida et al. 2021)—yield similar absolute magnitudes of $+3.01$ and

<footnote n="78">78 Version 3.6 at http://stev.oapd.inaf.it/cmd.xs</footnote>

<footnote n="79">79 http://basti-iac.oa-abruzzo.inaf.it/index.html</footnote>
+2.38 and an intrinsic color of 0.63. Taking the averages and using the observed baseline color of $m_{F606W} - m_{F814W} = 2.365$ (Table 1), we find a color excess of $E(m_{F606W} - m_{F814W}) = 1.73$. Based on the values of $A_\lambda/A_V$ given for the WFC3 filters at the PARSEC website (derived from the interstellar extinction curves of Cardelli et al. 1989 and O'Donnell 1994, with $R_V = 3.1$), this color excess corresponds to an extinction of $A_{F814W} = 3.33$ mag. Thus, using an apparent magnitude at baseline of $m_{F814W} = 19.581$ from Table 1, we find a true distance modulus of $(m - M)_0 = 13.86$, giving a linear distance of $D_S = 5.9 \pm 1.3$ kpc. The quoted error is the formal uncertainty (dominated by the uncertainty in the spectroscopic log g), but there likely are additional systematic errors, given the large amount of extinction, the assumption of a standard value of $R_V$, the high sensitivity to log g, and other sources. For example, if a low value of $R_V \simeq 2.5$, as deduced by Nataf et al. (2013) for the Galactic bulge region, were used, the distance would increase by $\sim 2.8$ kpc.

The nominal 5.9 kpc distance places the source on the near edge of the Galactic bulge. However, it cannot be ruled out that the source is actually located closer to the center of the Galactic bulge, at a distance of $\simeq 8$ kpc. The star's high radial velocity is consistent with stars in the bulge but not large enough to rule out disk membership.

Taking the spectroscopy, photometry, and derived distance at face value, these parameters correspond to an angular diameter of the source of $\sim 0.0038$ mas.

### 9.3. Finite-source Effects Are Negligible

We can estimate the effect of the finite size of a small source on the light curve by expanding $A(u)$ around $u = u_0$, giving

$$A(u) \simeq A(u_0) + A'(u_0)(u - u_0) + \frac{1}{2} A''(u_0)(u - u_0)^2, \quad (18)$$

where the prime symbol denotes the derivative. If we average over the face of a spherically symmetric source star, the term with the first derivative cancels out, leaving the quadratic term as the first one that contributes to the difference between finite and pointlike source. Averaging the latter over an extent $\delta$ explicitly gives

$$(\Delta A)_\delta(\delta) = \frac{1}{4\delta} A''(u_0) \int_{u_0-\delta}^{u_0+\delta} (u - u_0)^2 du = \frac{1}{6} A''(u_0) \delta^2. \quad (19)$$

For $u \ll 1$, one finds $A(u) \simeq 1/u$ and thereby $A''(u) \simeq 2/u^3$. Consequently, we find the largest differences between finite and pointlike sources for the smallest $u$. For a source of angular radius $\theta_*$, we find a source size parameter $\rho_* \equiv \theta_*/\theta_E$, and if $u_0 \gg \rho_*$, we can approximate the magnification difference by

$$\begin{aligned} \Delta A &\simeq \frac{1}{\rho_*} \int_0^{\rho_*} (\Delta A)_\delta (\sqrt{\rho_*^2 - \eta^2}) d\eta \\ &= \frac{1}{6} A''(u_0) \rho_*^2 \int_0^1 (1 - \xi^2) d\xi \\ &\simeq \frac{2}{9} \frac{\rho_*^2}{u_0^3}. \end{aligned} \quad (20)$$

With the angular diameter of the source $2\theta_* \simeq 0.0038$ mas and the angular Einstein radius $\theta_E = 5.2$ mas, we find $\rho_* \equiv \theta_*/\theta_E = 0.00037$ while $u_0 = 0.00271$. Inserting into Equation (20) gives $\Delta A = 1.5$ as compared to $A(u_0) = 369$ at closest angular approach between lens and source, in close agreement with a numerical evaluation of the average of the exact magnification function $A(u)$. With a blend ratio $g \simeq 20$, comparing $2.5 \log[(A + g)/(1 + g)]$ for A and $A + \Delta A$ gives a difference of 4 mmag between finite source and pointlike source, which is of the order of the systematic errors of the photometry.

Accounting for the finite source size will thus result in a very small change in $u_0$, which has very little effect on other parameters, evidenced by their practically identical values for our $u_{0,+}$ and $u_{0,-}$ solutions. Similar to the case of the long-duration microlensing event OGLE-2014-BLG-1186 discussed by Dominik et al. (2019), the measurement of the parallax parameter $\pi_E$ comes from the wings of the photometric light curve, not from the peak. Additionally, none of the astrometric deflection measurements are close to the peak. With the finite size of the source making little difference, we can safely neglect any limb-darkening effects.

## 10. Nature of the Lens

In this section we discuss the nature of the lensing object of the MOA-11-191/OGLE-11-462 microlensing event. We consider the mass of the lens, discuss constraints on its optical luminosity, and consider whether it is a single object or could be a binary system.

### 10.1. Mass

As shown by Equation (16), the mass of the MOA-11-191/ OGLE-11-462 lens can be determined from the values of its angular Einstein radius, $\theta_E$, and the relative lens-source parallax, $\pi_E$. In Section 8, we derived $\theta_E = 5.18 \pm 0.51$ mas and $\pi_E = 0.089 \pm 0.014$. These yield a mass of $M_{\text{lens}} = 7.1 \pm 1.3 M_\odot$.

An object with a mass this large cannot be a single (or double) NS or white dwarf. It can only be a BH—or an ordinary star, or conceivably a binary (or higher multiple) containing stars, BHs, and/or other compact companions. To distinguish between these possibilities, we consider the observational constraints on the luminosity and binarity of the lens.

### 10.2. Lens Distance

The mass determination for the lens is independent of assumptions about the distances to the lens and source. However, to constrain the optical luminosity of the lens, we first do need to determine its distance. This can be obtained from Equation (2), using the value of $\pi_{LS} = 0.463 \pm 0.051$ mas (from Section 8) and the distance to the source (from Section 9.2) $D_S = 5.9 \pm 1.3$ kpc.
<figure><img src="image_23.png" /><figcaption>**Figure 20.** Panels showing 10 × 10-pixel stacks in the reference frame of 4 × 250 s F814W images from the 2017 epoch, centered on the source. (a) The original `_flc` data with only the bright neighbor to the upper right subtracted. (b) Stack of the best-fit library PSF. (c) The original minus the library PSF. The subtraction is done in the `_flc` frame, and then the pixels are stacked. (d) The original minus the optimized PSF. (e) Same as the previous panel, but with a neighbor added that has 10% of the source brightness and offset by 42 mas. (f) The source+neighbor pixels fitted with a single-star PSF; the centroid moves up and to the right to account for the neighbor.</figcaption></figure>

These values yield a lens distance of $D_L = 1.58 \pm 0.18$ kpc. Note that the derived distance of the lens is only a weak function of the adopted source distance. If the source were at the distance of the Galactic center, the lens distance would only increase to ~1.7 kpc.

## 10.3. Luminosity Constraints

### 10.3.1. From the Baseline Magnitude

Main-sequence stars with a mass in the range $7.1 \pm 1.3 M_\odot$ have a range of *I*-band absolute magnitudes of $M_I = -1.5 \pm 0.3$, based on the tabulation<fnref n="80" /> of stellar properties compiled by Pecaut & Mamajek (2013).

The lens distance of $1.58 \pm 0.18$ kpc, from the previous subsection, corresponds to a true distance modulus of $(m - M)_0 = 11.0 \pm 0.2$. Thus, a main-sequence lens of the measured mass would have an unreddened apparent magnitude of $I \simeq 9.5$. The *I*-band extinction of the source (from Section 9.2) is ~3.3 mag, which is likely an overestimate for the lens, since it is considerably nearer than the source and may suffer less extinction. Nevertheless, this amount of extinction gives an expected apparent magnitude of a main-sequence lens of about $I \simeq 12.8$, or brighter if the extinction is lower. The star would be yet brighter if it has evolved off the main sequence to a higher luminosity. Since the baseline F814W magnitude of the source (plus lens) is ~19.6 (Table 1), a main-sequence or brighter stellar lens is conclusively ruled out.

### 10.3.2. Direct Limits from Final-epoch HST Imaging

A much stronger constraint on its optical luminosity comes from the fact that the lens was not detected in our final-epoch HST frames. The parameters $\theta_E = 5.18 \pm 0.51$ mas and $t_E = 270.7 \pm 11.2$ days imply that the proper motion of the lens with respect to the source is $6.99 \pm 0.75$ mas yr$^{-1}$. At peak magnification in 2011 July, the separation between lens and source was ~0.01 mas. Therefore, at the epoch of our last HST observation in 2017 (6.1 yr after the peak), the lens-source separation had increased to ~42.6 mas. At this large a separation, a luminous lens would cause the PSF of the source in our images to show clear signs of elongation. In order to search for such a distortion, we subtracted a “library” F814W PSF<fnref n="81" /> from the final-epoch F814W images (in the individual

unresampled `_flc` frames). We saw no indication of elongation.

In order to place the most stringent limit possible on the lens brightness, we fit the above-referenced PSF to stars of similar brightness and color to that of the source and adjusted the shape of the PSF to optimize the fit (again, in the unresampled frame). We saw a ~1% variation of the PSF, which is typical of minor, “breathing”-related changes in telescope focus.

Figure 20 shows stacks of the source in our 2017 data at successive stages of analysis and illustrates the process of searching for photometric evidence of the lens. In particular, panel (d) shows the subtraction with the improved PSF. For illustration purposes, we show in panel (f) the residual pattern expected for a lens that is 10% of the brightness of the source. There is no such signature visible in panel (d). We quantify this as follows.

Using the improved PSF, we modeled the source image as a superposition of two stars separated by 42 mas (about 1 detector pixel) along the implied direction of proper motion. We did this simultaneously in the 4 × 250 s F814W images from 2017. We found a best-fit flux for the lens of −1% of the source flux. We performed the same fits to similar-brightness neighboring stars and found a distribution between −2% and +1%. Finally, we added to each neighbor star a companion of 2.5% at the presumed location and were able to recover the added flux at a level of 2.5% ± 1%. We would clearly detect a lens if it were there.

We therefore conclude that the brightest the lens could be at the presumed location in 2017 is ~1% of the source flux. We also explored how bright a lens could be at an offset of 42 mas in *any* direction from the source. We were able to rule out any lens brighter than 2.5% of the source flux.

A luminosity of ~1% of that of the source corresponds to an apparent F814W or *I* magnitude of ~24.6. At the distance and extinction of the lens, the limit on its absolute magnitude is $M_I \gtrsim 10.2$ (corresponding to a main-sequence star of about $0.2 M_\odot$ or less, or a white dwarf of $\approx 0.6 M_\odot$ with a cooling age older than ~$10^7$ yr).

These considerations leave a single BH—or a BH with a mass greater than ~3.55 $M_\odot$ in a binary paired with another BH, or with an NS or faint star—as the only viable possibilities for the lens.

## 10.4. Constraints on Binarity

Having established that the MOA-11-191/OGLE-11-462 lens is an extremely faint, or nonluminous, object with a mass

<footnote n="80">⁸⁰ Online version of 2021 March 2, at http://www.pas.rochester.edu/~emamajek/EEM_dwarf_UBVIRHK_colors_Teff.txt.</footnote>

<footnote n="81">⁸¹ From https://www.stsci.edu/hst/instrumentation/wfc3/data-analysis/psf.</footnote>
of $7.1 \pm 1.3 M_\odot$, we now consider in more detail whether the lens could be a binary system rather than a single, isolated BH.

MOA-11-191/OGLE-11-462 was a high-magnification event, making its light curve especially sensitive to lens binarity. As described above, the event was intensely monitored photometrically by several microlensing groups in search of planetary-companion signatures, and no manifestation of a binary lens was seen.

In fact, lens binarity will strongly affect both the photometric and the astrometric signature if the angular separation between the components is around the angular Einstein radius, $\theta_E$, and it therefore has to be either much smaller or much larger. In either case, such a binary lens would produce a diamond-shaped caustic.

For a close binary, this caustic would be located around the center of mass of the lens, and both the photometric and astrometric microlensing signatures would be strongly altered if the source–lens trajectory passed over it or came close. Specifically, for a mass ratio $q$ and a separation $d\theta_E$, the caustic extends to (Schneider & Weiß 1986; Erdl & Schneider 1993; Dominik 1999)

$$s_{\pm}^{c} = \pm \frac{2q}{(1+q)^2} d^2. \qquad (21)$$

Comparing this to $u_0 \simeq 0.00271$ yields the constraint $d < 0.07$ for the equal-mass case, with less than 11% variation over the range $q \in [0.4, 1]$. With $D_L \theta_E = 8.2$ au, this implies a separation smaller than 0.6 au. Considering a systematic noise floor of ~5 mmag, deviations from a point-lens light curve exceed this for the given $u_0 \simeq 0.00271$ and blend ratio $g \simeq 20$ if $d > 0.022$, so that our data imply an upper limit to the separation of ~0.18 au.

In such a case, since one of the components has to have a mass of at least $3.55 \pm 0.55 M_\odot$, the merging timescale for such a system is $\lesssim 10^7$ yr (Carroll & Ostlie 2006). Given that there is no evidence of active star formation in this location, this seems unlikely.

In the case of widely separated binary components, the observed properties of the lens refer to one of them, leading to the same values as by assuming that the lens is a single isolated body, so that the mass determination of the putative BH remains valid.

With $\theta_E$ referring to the mass of the object identified, in the wide-binary case one finds the caustic extending to (e.g., Dominik 1999)

$$s_{\pm}^{w} = \pm \frac{q}{d^2}. \qquad (22)$$

Considering that any companion should be dark, as well as the noise floor of our photometric data, we exclude any companion more massive than 10% of the identified lens within 230 au.

To further quantify the limits on a binary companion to the lens, we performed a Rhie et al. (2000)-type analysis in which we simulate hypothetical binary-lens light curves with the same epochs and properties as the real data. We then fit those light curves with a point-lens model to determine the potential signal from a given lens companion. For this test, we use a subset of the data that provides good coverage of the overall light curve and that is sufficient to provide dense coverage of the peak, where a potential companion would have the greatest effect (Griest & Safizadeh 1998). Specifically, we use the OGLE and MOA survey data sets and the microFUN data from Auckland

<figure><img src="image_24.png" /><figcaption>Figure 21. CMD for a sample of stars in the Galactic disk at similar distances to that of the BH, selected from Gaia EDR3 as described in the text. The population at this location contains main-sequence stars of a younger population and older stars with redder colors.</figcaption></figure>

Observatory, Kumeu Observatory, Vintage Lane Observatory, the Wise Observatory 0.46 m, and the Weizmann Institute. Setting the requirement that the $\chi^2$ due to the planet is >50, we find that the sensitivity limits can be approximated with a piecewise function

$$|\log(d_{\text{lim}})| = \begin{cases} 1.54 & \text{for } \log q > -0.4 \\ 0.3413 \log q + 1.677 & \text{for } -4.8 < \log q \le -0.4 \\ 0.04 & \text{for } \log q \le -4.8 \end{cases}. \qquad (23)$$

Thus, for $q = 0.1$, this analysis also excludes companions with separations $0.4$ au $\lesssim a_\perp \lesssim 180$ au.

## 10.5. The Lens Is an Isolated Black Hole

To summarize this section, we have found that the lens that produced the MOA-11-191/OGLE-11-462 microlensing event has a mass of $7.1 \pm 1.3 M_\odot$. There are stringent limits on its optical luminosity, ruling out the possibility that the lens could be an ordinary star of that mass. Its mass is greater than possible for a white dwarf or NS, or for a binary pair of white dwarfs and/or NSs. Even if the lens is a binary, at least one of its components must still be a BH—but only under unlikely circumstances could it actually be a binary at all. In the rest of our discussion, we will assume that the lens is an isolated, single BH.

## 11. Properties of the Black Hole

Table 6 summarizes the inferred properties of the BH lens that produced the MOA-11-191/OGLE-11-462 photometric and astrometric microlensing event. In this section we discuss some properties of the lens and its Galactic environment.
<figure><img src="image_25.png" /><figcaption>**Figure 22.** Proper-motion components from Gaia EDR3 for a sample of stars at similar distances to that of the BH, selected as described in the text. Typical uncertainties are about ±0.05 mas yr⁻¹ in both coordinates. Blue points represent a younger main-sequence population with Gaia colors of G<sub>BP</sub> − G<sub>RP</sub> < 1.0, and red points represent an older and dynamically hotter population. The directions of Galactic longitude and latitude are indicated at the lower left. The black point with error bars shows the proper motion of the BH, which is a significant outlier compared to the surrounding population.</figcaption></figure>

## 11.1. Galactic Location

The BH lies in the direction of, but is considerably closer than, the Galactic bulge. The fifth row in Table 6 gives the rectangular components of its Galactic position relative to the Sun. The BH is located in the Galactic disk at a distance from us of ~1580 pc toward the Galactic center and about 45 pc below the Galactic plane. This places it approximately between the Scutum–Centaurus and Sagittarius–Carina spiral arms of the Milky Way, in the terminology of Reid et al. (2019). There are no conspicuous star-forming regions, young stellar objects, or SN remnants at this location.

## 11.2. The Surrounding Stellar Population

Most of the stars contained in the small FOV of our HST frames lie at considerably larger distances than the BH itself. In order to study a larger sample of the relatively nearby stellar population surrounding the spatial location of the BH, we selected stars from Gaia EDR3 that satisfy the following criteria: they lie within 20' of the BH on the sky, have trigonometric parallaxes with fractional errors of less than 10% that fall in the range 0.55–0.70 mas, and are brighter than G = 17.

Figure 21 shows the CMD of this population in the Gaia magnitude and color system. The CMD is afflicted by severe differential extinction, which blurs the stellar sequences. Nevertheless, it shows two prominent components: a moderately young and sparse population of main-sequence stars bluer than G<sub>BP</sub> − G<sub>RP</sub> ≃ 1.0, and a denser older component of redder stars. There is little or no evidence for an extremely young population of massive stars, at the BH's location. However, with a transverse proper motion of ~45 km s⁻¹ (Table 6), it would take only 0.2 Myr for the lens to traverse 20', so it is possible that the lens originated much farther away.

## 11.3. Space Velocity, Population Constraints, Natal Kick

The sixth row in Table 6 gives the components of the BH's space velocity relative to the Sun in the V and W directions (in the directions of Galactic rotation and toward the north Galactic pole, respectively). The third component of the space velocity, U, is undetermined since the radial velocity of the BH is unknown; changing the radial velocity to −100 and +100 km s⁻¹ changes V and W only by about ±0.2 and ±2.8 km s⁻¹, respectively. The relatively low value of V indicates that the BH is not a member of an extreme Galactic halo population and may be moving in a roughly circular Galactic orbit (unless the radial velocity is large). However, it does have a moderately high W velocity perpendicular to the Galactic plane, so the orbit has a modest inclination to the Galactic plane. Apart from these considerations, there are few constraints on the age of the BH, except that it has likely had time to move far from its birthplace.

Figure 22 plots the absolute proper motion of the BH (black filled circle with error bars, from the fourth row of Table 6) and compares it with the proper motions of the neighboring Gaia stars from the sample described above. The directions of Galactic longitude and latitude are indicated at the lower left, showing that the dispersions are higher in the longitude direction than perpendicular to it. Members of the “young” population from Figure 21 are plotted with blue points, and members of the older population are plotted with red points. The older population, which has proper-motion dispersions of σ<sub>μ,l</sub> = 2.59 mas yr⁻¹ and σ<sub>μ,b</sub> = 1.59 mas yr⁻¹, is dynamically hotter than the younger population, which has dispersions of σ<sub>μ,l</sub> = 1.08 mas yr⁻¹ and σ<sub>μ,b</sub> = 1.05 mas yr⁻¹.

The BH itself is a prominent outlier, predominantly in b, relative to the older and younger populations (~4σ and 6σ, respectively). This suggests that it may have received a “natal kick” from the SN explosion associated with its birth—assuming that the BH arose from the surrounding population. Relative to the mean proper motion of the younger population, the BH’s motion is offset by about 6.0 mas yr⁻¹, corresponding to a tangential velocity offset of ~45 km s⁻¹. However, considering that its radial velocity is unknown, we cannot exclude the possibility that the BH might simply be passing through the surrounding population and could conceivably be significantly older. In a recent paper, Andrews & Kalogera (2022) have given an extensive discussion of the possible natal kick of MOA-11-191/OGLE-11-462.

# 12. Discussion and Future Work

## 12.1. The Formation of Isolated Black Holes

As reviewed in Section 1.1, it is generally believed that only stars above ~20 M<sub>☉</sub> collapse to form BHs. To form an isolated ~7 M<sub>☉</sub> BH, such as the one discussed in this paper, either its progenitor was a single star that experienced considerable mass loss through stellar winds, it was produced through fallback in a weak SN explosion (e.g., Fryer et al. 2012), or it was born in a close binary that lost mass through a common-envelope mass ejection. In the binary scenario, the binary evidently became unbound, presumably when its lower-mass companion itself became an SN, and its mass loss and SN kick detached the binary.

The MOA-11-191/OGLE-11-462 BH's inferred mass of ~1 M<sub>☉</sub> lies remarkably close to the peak of the Galactic X-ray-binary BH mass distribution (see Section 1.1). This may suggest a similar evolutionary path with mass loss through
close-binary interactions, but a larger sample of mass measure-
ments for isolated BHs is needed to distinguish these different
formation scenarios.

12.2. *Future Work*

There is no known X-ray or radio source at the position of MOA-11-191/OGLE-11-462. Considering the high extinction toward the BH, it is possible that it lies in a region of high density of interstellar matter. If so, the accretion rate of material from the interstellar matter by the BH could conceivably be large enough for it to be detected in deep X-ray and/or radio observations, as discussed in Section 1.2, and these would be worth attempting.

Finally, we note the potential of future facilities, such as the Very Large Telescope Interferometer (VLTI), the Nancy Grace Roman Space Telescope, and the Rubin Observatory, in this area. The GRAVITY instrument with VLTI is currently being upgraded to further improve its sensitivity (Eisenhauer 2019), which will make it possible to measure θ_E from spatially resolved microlensed images (Dong et al. 2019) of a large number of microlensed sources in the Milky Way. Roman will detect and characterize several thousand microlensing events (Spergel et al. 2015; Penny et al. 2019); with the same aperture size as HST, its astrometric capabilities are expected to be similar (WFIRST Astrometry Working Group et al. 2019). According to current plans, Roman will observe microlensing events over a ~2 deg² region in the Galactic bulge with a ~15-minute cadence for several months each year; thus, it will make thousands of photometric and astrometric observations of every long-duration event, leading to very high accuracy in both kinds of measurements. With such data, Roman should detect numerous BHs, yielding exquisite determinations of their masses, distances, and velocities, and providing important insights to our understanding of the formation and evolution of BHs. Rubin's wide-angle survey will provide deep, long-baseline photometry and consistent data reduction, making it ideal for the discovery and photometric characterization of microlensing events (Sajadian & Poleski 2019), probing many different lines of sight (and hence stellar populations) across the Milky Way and Magellanic Clouds (Street et al. 2018b). In the bulge, it will complement the Roman data by filling in gaps in the Roman survey cadence (Street et al. 2018a; Strader et al. 2018). Rubin discoveries may need high-resolution time-series imaging follow-up for full characterization of the events.

Based in part on observations made with the NASA/ESA Hubble Space Telescope, obtained at STScI, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. Support for this research was provided by NASA through grants from STScI. HST data used in this paper are available from the Mikulski Archive for Space Telescopes at STScI,<fnref n="82" /> under proposal IDs 12322, 12670, 12986, 13458, and 14783.

This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/ gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/ dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participat- ing in the Gaia Multilateral Agreement.

The MOA project is supported by JSPS KAK-ENHI grant Nos. JSPS24253004, JSPS26247023, JSPS23340064, JSPS15H00781, JP16H06287, 17H02871, and 19KK0082.

We acknowledge the help and dedication of the late Dr. John Greenhill, who, as part of the PLANET collaboration, played a key role in the efforts of collaboration in general, and in particular running the Tasmania observatory operations.

K.C.S. spent some time at the European Southern Observa-
tory, Institute for Advanced Study, and Harvard-Smithsonian
Center for Astrophysics, where some of this work was done,
and thanks them for their hospitality.

Ł.W. acknowledges support from the Polish NCN grant Daina No. 2017/27/L/ST9/03221.

J.-P.B. acknowledges support by the University of Tasmania through the UTAS Foundation and the endowed Warren Chair in Astronomy, and support by ANR COLD-WORLDS (ANR-18-CE31-0002) at Le Centre National de la Recherche Scientifique (CNRS) in Paris and the Laboratoire d'astrophysique de Bordeaux.

U.G.J. acknowledges funding from the European Union H2020-MSCA-ITN-2019 under grant No. 860470 (CHAMELEON) and from the Novo Nordisk Foundation Interdisciplinary Synergy Programme grant no. NNF19OC0057374.

L.M. acknowledges support from the ``Fondi di Ricerca Scientificad'Ateneo 2021'' of the University of Rome ``Tor Vergata.''

D.D.P. acknowledges support by the University of Rijeka through the grant uniri-prirod-18-48.

T.C.H. acknowledges financial support from the National Research Foundation (NRF; No. 2019R111A1A 01059609).

Y.T. acknowledges the support of DFG priority program SPP 1992 "Exploring the Diversity of Extrasolar Planets" (TS 356/3-1).

*Note added:* After our paper was submitted, a study by Lam et al. (2022) that includes an independent investigation of MOA-11-191/OGLE-11-462 was posted on arXiv. We have not used any of the measurements or results from that paper in our analysis.

Facilities: HST (WFC3).

ORCID iDs

Kailash C. Sahu https://orcid.org/0000-0001-6008-1955
Jay Anderson https://orcid.org/0000-0003-2861-3995
Howard E. Bond https://orcid.org/0000-0003-1377-7145
Andrzej Udalski https://orcid.org/0000-0001-5207-5619
Martin Dominik https://orcid.org/0000-0002-3202-0343
Annalisa Calamida https://orcid.org/0000-0002-0882-7702
Andrea Bellini https://orcid.org/0000-0003-3858-637X
Thomas M. Brown https://orcid.org/0000-0002-1793-9968
Marina Rejkuba https://orcid.org/0000-0002-6577-2787
Noé Kains https://orcid.org/0000-0001-8803-6769
Henry C. Ferguson https://orcid.org/0000-0001-7113-2738
Chris L. Fryer https://orcid.org/0000-0003-2624-0056
Philip Yock https://orcid.org/0000-0001-9716-7752
Szymon Kozłowski https://orcid.org/0000-0003-4084-880X
Paweł Pietrukowicz https://orcid.org/0000-0002-2339-5899
Radek Poleski https://orcid.org/0000-0002-9245-6368
Jan Skowron https://orcid.org/0000-0002-2335-1730
Igor Soszyński https://orcid.org/0000-0002-7777-0842
Michał K. Szymański https://orcid.org/0000-0002-0548-8995
Krzysztof Ulaczyk https://orcid.org/0000-0001-6364-408X
Łukasz Wyrzykowski https://orcid.org/0000-0002-9658-6151

<footnote n="82">⁸² https://archive.stsci.edu/hst/search.php</footnote>
Richard K. Barry https://orcid.org/0000-0003-4916-0892

David P. Bennett https://orcid.org/0000-0001-8043-8413

Yuki Hirao https://orcid.org/0000-0003-4776-8618

Stela Ishitani Silva https://orcid.org/0000-0003-2267-1246

Iona Kondo https://orcid.org/0000-0002-3401-1029

Naoki Koshimoto https://orcid.org/0000-0003-2302-9562

Clément Ranc https://orcid.org/0000-0003-2388-4534

Nicholas J. Rattenbury https://orcid.org/0000-0001-5069-319X

Takahiro Sumi https://orcid.org/0000-0002-4035-5021

Daisuke Suzuki https://orcid.org/0000-0002-5843-9433

Aikaterini Vandorou https://orcid.org/0000-0002-9881-4760

Andrew Cole https://orcid.org/0000-0003-0303-3855

Pascal Fouqué https://orcid.org/0000-0002-1436-7351

Dijana Dominis-Prester https://orcid.org/0000-0002-9880-5039

Étienne Bachelet https://orcid.org/0000-0002-6578-5078

Michael Albrow https://orcid.org/0000-0003-3316-4012

Jennifer C. Yee https://orcid.org/0000-0001-9481-7123

Leonardo A. Almeida https://orcid.org/0000-0002-3817-6402

Avishay Gal-Yam https://orcid.org/0000-0002-3653-5598

Chung-Uk Lee https://orcid.org/0000-0003-0043-3025

Richard W. Pogge https://orcid.org/0000-0003-1435-3053

Yossi Shvartzvald https://orcid.org/0000-0003-1525-5041

Uffe G. Jørgensen https://orcid.org/0000-0001-7303-914X

Valerio Bozza https://orcid.org/0000-0003-4590-0136

Martin Burgdorf https://orcid.org/0000-0002-5854-4217

Tobias C. Hinse https://orcid.org/0000-0001-8870-3146

Markus Hundertmark https://orcid.org/0000-0003-0961-5231

Eamonn Kerins https://orcid.org/0000-0002-1743-4468

Luigi Mancini https://orcid.org/0000-0002-9428-8732

Sohrab Rahvar https://orcid.org/0000-0002-7084-5725

Davide Ricci https://orcid.org/0000-0002-9790-0552

Sedighe Sajadian https://orcid.org/0000-0002-2859-1071

Jesper Skottfelt https://orcid.org/0000-0003-1310-8283

Colin Snodgrass https://orcid.org/0000-0001-9328-2905

John Southworth https://orcid.org/0000-0002-3807-3198

Joachim Wambsganss https://orcid.org/0000-0002-8365-7619

Yiannis Tsapras https://orcid.org/0000-0001-8411-351X

Rachel A. Street https://orcid.org/0000-0001-6279-0552

Keith Horne https://orcid.org/0000-0003-1728-0304

## References

Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016, *PhRvL*, 116, 061102

Abbott, R., Abbott, T. D., Abraham, S., et al. 2021a, *PhRvX*, 11, 021053

Abbott, R., Abbott, T. D., Acernese, F., et al. 2021b, arXiv:2111.03606v1

Agol, E., & Kamionkowski, M. 2002, *MNRAS*, 334, 553

Alcock, C., Allisman, R. A., Alves, D., et al. 1995, *ApJ*, 454, L125

An, J. H., Albrow, M. D., Beaulieu, J.-P., et al. 2002, *ApJ*, 572, 521

Anderson, J. 2021, WFC3 Instrument Science Report, 2021-13 (https://www.stsci.edu/hst/instrumentation/wfc3/documentation/instrument-science-reports-isrs)

Anderson, J., & King, I. R. 2006, ACS Instrument Science Report, 2006-01 (https://www.stsci.edu/hst/instrumentation/acs/documentation/instrument-science-reports-isrs)

Anderson, J., King, I. R., Richer, H. B., et al. 2008, *AJ*, 135, 2114

Andrews, J. J., & Kalogera, V. 2022, *ApJ*, 930, 159

Belczynski, K., Repetto, S., Holz, D. E., et al. 2016, *ApJ*, 819, 108

Bellini, A., Anderson, J., & Bedin, L. R. 2011, *PASP*, 123, 622

Bellini, A., Vesperini, E., Piotto, G., et al. 2015, *ApJL*, 810, L13

Bennett, D. P., Becker, A. C., Quinn, J. L., et al. 2002, *ApJ*, 579, 639

Bensby, T., Yee, J. C., Feltzing, S., et al. 2013, *A&A*, 549, A147

Binder, B. A., Sy, J. M., Eracleous, M., et al. 2021, *ApJ*, 910, 74

Bond, I. A., Abe, F., Dodd, R. J., et al. 2001, *MNRAS*, 327, 868

Bramich, D. M. 2018, *A&A*, 618, A44

Bressan, A., Marigo, P., Girardi, L., et al. 2012, *MNRAS*, 427, 127

Brown, G. E., & Bethe, H. A. 1994, *ApJ*, 423, 659

Brown, T. M., Casertano, S., Strader, J., et al. 2018, *ApJ*, 856, L6

Brown, W. R., Anderson, J., Gnedin, O. Y., et al. 2015, *ApJ*, 804, 49

Calamida, A., Mack, J., Medina, J., et al. 2021, WFC3 Instrument Science Report, 2021-04 (https://www.stsci.edu/hst/instrumentation/wfc3/documentation/instrument-science-reports-isrs)

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, *ApJ*, 345, 245

Carroll, B. W., & Ostlie, D. A. 2006, An Introduction to Modern (II ed.; San Francisco, CA: Addison-Wesley)

Casares, J., & Jonker, P. G. 2014, *SSRv*, 183, 223

Casertano, C., Riess, A. G., Anderson, J., et al. 2016, *ApJ*, 825, 11

Chawla, S., Chatterjee, S., Breivik, K., et al. 2021, arXiv:2110.05979

Corral-Santana, J. M., Casares, J., Muñoz-Darias, T., et al. 2016, *A&A*, 587, A61

Cromartie, H. T., Fonseca, E., Ransom, S. M., et al. 2020, *NatAs*, 4, 72

de Mink, S. E., Sana, H., Langer, N., et al. 2014, *ApJ*, 782, 7

Dominik, M. 1999, *A&A*, 349, 108

Dominik, M., Bachelet, E., Bozza, V., et al. 2019, *MNRAS*, 484, 5608

Dominik, M., Belczynski, K., Fryer, C., et al. 2012, *ApJ*, 759, 52

Dominik, M., Jørgensen, U. G., Rattenbury, N. J., et al. 2010, *AN*, 331, 671

Dominik, M., Rattenbury, N. J., Allan, A., et al. 2007, *MNRAS*, 380, 792

Dominik, M., & Sahu, K. C. 2000, *ApJ*, 534, 213

Dong, S., Mérand, A., Delplancke-Ströbele, F., et al. 2019, *ApJ*, 871, 70

Dressel, L. 2021, Wide Field Camera 3 Instrument Handbook for Cycle 29, v13.0 (Baltimore MD: STScI), https://hst-docs.stsci.edu/wfc3ibh

Dyson, F. W., Eddington, A. S., & Davidson, C. 1920, *RSPTA*, 220, 291

Einstein, A. 1936, *Sci*, 84, 506

Eisenhauer, F. 2019, The Very Large Telescope in 2030 (Munich: ESO Garching)

Erdl, H., & Schneider, P. 1993, *A&A*, 268, 453

Farr, W. M., Sravan, N., Cantrell, A., et al. 2011, *ApJ*, 741, 103

Fender, R. P., Maccarone, T. J., & Heywood, I. 2013, *MNRAS*, 430, 1538

Fryer, C. L., Belczynski, K., Wiktorowicz, G., et al. 2012, *ApJ*, 749, 91

Fryer, C. L., & Kalogera, V. 2001, *ApJ*, 554, 548

Fryer, C. L., Woosley, S. E., & Hartmann, D. H. 1999, *ApJ*, 526, 152

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2021, *A&A*, 649, A1

Gaudí, B. S. 2012, *ARA&A*, 50, 411

Giesers, B., Kamann, S., Dreizler, S., et al. 2019, *A&A*, 632, A3

Gould, A. 1992, *ApJ*, 392, 442

Gould, A. 2004, *ApJ*, 606, 319

Gould, A., Udalski, A., An, D., et al. 2006, *ApJL*, 644, L37

Griest, K., & Safizadeh, N. 1998, *ApJ*, 500, 37

Heger, A., Fryer, C. L., Woosley, S. E., et al. 2003, *ApJ*, 591, 288

Hidalgo, S. L., Pietrinferni, A., Cassisi, S., et al. 2018, *ApJ*, 856, 125

Hög, E., Novikov, I. D., & Polanarev, A. G. 1995, *A&A*, 294, 287

Janssens, S., Shenar, T., Sana, H., et al. 2022, *A&A*, 658, A129

Jayasinghe, T., Stanek, K. Z., Thompson, T. A., et al. 2021, *MNRAS*, 504, 2577

Kains, N., Calamida, A., Sahu, K. C., et al. 2017, *ApJ*, 843, 145

Kim, S.-L., Lee, C.-U., Park, B.-G., et al. 2016, *JKAS*, 49, 37

Kreidberg, L., Bailyn, C. D., Far, W. M., et al. 2012, *ApJ*, 757, 36

Lam, C. Y., Lu, J. R., Udalski, A., et al. 2022, *ApJL*, 933, L23

Laycock, S. G. T., Maccarone, T. J., & Christodoulou, D. M. 2015, *MNRAS*, 452, L31

Linares, M., Shahbaz, T., & Casares, J. 2018, *ApJ*, 859, 54

Lu, J. R., Sinukoff, E., Ofek, E. O., et al. 2016, *ApJ*, 830, 41

Mao, S., & Paczyński, B. 1991, *ApJL*, 374, L37

Mao, S., Smith, M. C., Woźniak, P., et al. 2002, *MNRAS*, 329, 349

Marigo, P., Girardi, L., Bressan, A., et al. 2017, *ApJ*, 835, 77

Miller-Jones, J. C. A., Bahramian, A., Orosz, J. A., et al. 2021, *Sci*, 371, 1046

Minniti, D., Contreras Ramos, R., Alonso-García, J., et al. 2015, *ApJL*, 810, L20

Miyamoto, M., & Yoshii, Y. 1995, *AJ*, 110, 1427

Nataf, D. M., Gould, A., Fouqé, P., et al. 2013, *ApJ*, 769, 88

O'Donnell, J. E. 1994, *ApJ*, 422, 158

Orosz, J. A., McClintock, J. E., Narayan, R., et al. 2007, *Natur*, 449, 872

Özel, F., Psaltis, D., Narayan, R., et al. 2010, *ApJ*, 725, 1918

Paczyński, B. 1996, *ARA&A*, 34, 419

Paczyński, B. 1998, *ApJL*, 494, L23

Pecaut, M. J., & Mamajek, E. E. 2013, *ApJS*, 208, 9

Penny, M. T., Gaudí, B. S., Kerins, E., et al. 2019, *ApJS*, 241, 3

Pietrinferni, A., Hidalgo, S., Cassisi, S., et al. 2021, *ApJ*, 908, 102

Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2019, *ApJ*, 885, 131

Remillard, R. A., & McClintock, J. E. 2006, *ARA&A*, 44, 49

Rhie, S. H., Bennett, D. P., Becker, A. C., et al. 2000, *ApJ*, 533, 378

Riess, A. G., Casertano, S., Yuan, W., et al. 2018, *ApJ*, 855, 136

Sackett, P. D., Albrow, M. D., Beaulieu, J.-P., et al. 2004, *IAU Symp.* 213, *Bioastronomy 2002: Life Among the Stars* (San Francisco, CA: ASP), 35
Sahu, K. C., Anderson, J., Casertano, S., et al. 2017, *Sci*, 356, 1046

Sahu, K. C., Bond, H. E., Anderson, J., et al. 2014, *ApJ*, 782, 89

Sahu, K. C., Anderson, J., Baggett, S., et al. 2021, WFC3 Data Handbook, v5.0 (Baltimore MD: STScI), https://hst-docs.stsci.edu/wfc3dhb

Sajadian, S., & Poleski, R. 2019, *ApJ*, 871, 205

Samland, M. 1998, *ApJ*, 496, 155

Sana, H., de Mink, S. E., de Koter, A., et al. 2012, *Sci*, 337, 444

Scarcella, F., Gaggero, D., Connors, R., et al. 2021, *MNRAS*, 505, 4036

Schneider, P., & Weißenbach, A. 1986, A&A, 164, 237

Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White Dwarfs, and Neutron Stars: The Physics of Compact Objects (New York: Wiley-Interscience), 663

Shvartzvald, Y., Maoz, D., Udalski, A., et al. 2016, *MNRAS*, 457, 4089

Silverman, J. M., & Filippenko, A. V. 2008, *ApJ*, 678, L17

Sohn, S. T., Anderson, J., & van der Marel, R. P. 2012, *ApJ*, 753, 7

Spera, M., Mapelli, M., & Bressan, A. 2015, *MNRAS*, 451, 4086

Spergel, D., Gehrels, N., Baltay, C., et al. 2015, arXiv:1503.03757

Strader, J., Aydi, E., Britt, C., et al. 2018, arXiv:1811.12433

Street, R. A., Lund, M. B., Donachie, M., et al. 2018a, arXiv:1812.04445

Street, R. A., Lund, M. B., Khakpash, S., et al. 2018b, arXiv:1812.03137

Sukhbold, T., Ertl, T., Woosley, S. E., et al. 2016, *ApJ*, 821, 38

Tauris, T. M., & van den Heuvel, E. P. J. 2006, Compact Stellar X-ray Sources (Cambridge: Cambridge University Press), 623

Thompson, T. A., Kochanek, C. S., Stanek, K. Z., et al. 2019, *Sci*, 366, 637

Tsapras, Y., Street, R., Horne, K., et al. 2009, *AN*, 330, 4

Tsuna, D., & Kawanaka, N. 2019, *MNRAS*, 488, 2099

Tutukov, A. V., Fedorova, A. V., & Cherepashchuk, A. M. 2011, *ARep*, 55, 247

Udalski, A., Szymanski, M., Kaluzny, J., et al. 1992, *AcA*, 42, 253

Udalski, A., Szymański, M. K., & Szymański, G. 2015, *AcA*, 65, 1

van den Heuvel, E. P. J. 1992, in Environment Observation and Climate Modelling through International Space Projects, ed. T. Guyenne & J. Hunt (Noordwijk: ESA/ESTEC), 77

Walker, M. A. 1995, *ApJ*, 453, 37

WFIRST Astrometry Working Group, Sanderson, R. E., Bellini, R. E., et al. 2019, *JATIS*, 5, 044005

Woooley, S. E., Heger, A., & Weaver, T. A. 2002, *RvMP*, 74, 1015

Wyrzykowski, Ł., & Mandel, I. 2020, *A&A*, 636, A20

Wyrzykowski, Ł., Rynkiewicz, A. E., Skowron, J., et al. 2015, *ApJS*, 216, 12

Tsapras, Y. 2018, *Geosc*, 8, 365

Yoo, J., DePoy, D. L., Gal-Yam, A., et al. 2004, *ApJ*, 603, 139

Zang, W., Dong, S., Gould, A., et al. 2020, *ApJ*, 897, 180

Zhang, W., & Fryer, C. L. 2001, *ApJ*, 550, 357

Zurlo, A., Gratton, R., Mesa, D., et al. 2018, *MNRAS*, 480, 236