An Isolated Stellar-mass Black Hole Detected through Astrometric Microlensing

We report the ﬁ rst 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 ; 270 days ) , high-magni ﬁ cation microlensing event MOA-2011-BLG-191 / OGLE-2011-BLG-0462 ( 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 de ﬂ ection 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 ± 1.3 M e and a distance of 1.58 ± 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, con ﬁ rms 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 ∼ 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 ﬁ rst for an isolated stellar-mass BH using any technique.


Black Holes in Binary Systems
Stars with initial masses greater than ∼20 M 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 8 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 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 andCasares &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 "elec- † PLANET Collaboration, MiNDSTEp Consortium ‡ OGLE Collaboration § MiNDSTEp Consortium, RoboNet Collaboration ¶ RoboNet Collaboration * * PLANET Collaboration tromagnetically measured" BH masses show a distribution peaking near 7-8 M , with few if any below ∼5 M (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.3M (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 .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. 2021, and references therein) and measuring their masses.At the high-mass end, the BH mass distribution falls off above ∼10 M , and very few electromagnetic BH masses are known above ∼15 M , the only exception in the Milky Way being an updated mass determination of 21.1±2.2M 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 and 17 ± 4 M 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 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,b), the inferred masses of the BHs among the pre-merger systems range from ≈6 to 95 M , with two low-mass outliers among the secondary components at ∼2.6 and 2.8 M , which could be either BHs or NSs.

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 astromet-ric 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.

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 θ E , given by where M lens is the mass of the lens, and 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 θ 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 µ LS between lens and source, we can straightforwardly define t E , the time for the source to traverse an angular distance of θ E in the barycentric reference frame, as t E = θ E /µ LS (more details in §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 non-negligible 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 θ E is proportional to the square root of the lens mass.So, all else being equal, events due to massive compact objects would preferentially tend to be characterized by longer event durations (t E 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, non-luminous BH lens with a large Einstein radius-but it could alternatively be due simply to an unusually slow-moving, faint, lowmass 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 θ 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.

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 θ θ θ S and θ θ θ L denoting the angular positions of the source and lens as seen by the observer, we can define a dimensionless source-lens separation 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 which depends solely on u = |u u 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 §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 where F base ≡ F S + F B is the baseline flux and g ≡ F B /F S is the specific blend ratio.

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 (Gould 1992;Paczyński 1998).In contrast to the photometric microlensing signature, the astrometric signature is explicitly proportional to the angular Einstein radius θ E .Moreover, while the magnification diverges (for a point-like source) as u → 0, the astrometric shift becomes maximal for u = √ 2. The light magnification falls rapidly with increasing u: for large separations, u 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), andBramich (2018).
The photometric and astrometric signatures of a microlensing event are connected because they arise from the same source-lens trajectory, u u 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 θ E , as well as the orientation angle of the trajectory.

Parallax Effect and Proper Motion
For constructing the source-lens trajectory u u u(t), we have to consider the proper motions µ µ µ S and µ µ µ L of the source and lens objects, as well as their parallaxes, π S and π L .
Let γ γ γ(t) (1 AU) denote the projection of the 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 (cf.An et al. 2002;Gould 2004) so that for θ θ θ(t) ≡ θ θ θ S (t) − θ θ θ L (t), one finds where µ µ µ LS ≡ µ µ µ L −µ µ µ S and π LS ≡ π L − π S are the relative proper motion and relative parallax between lens and source, while (θ Consequently, with the microlensing parallax parameter π E = π LS /θ E , u u u(t) ≡ θ θ θ(t)/θ E takes the form (cf. Dominik et al. 2019) where as well as By construction, δγ γ γ(t 0 ) = 0 and δ γ γ γ(t 0 ) = 0.By choosing t 0 so that u u u 0 ⊥ u u u 0 , we can write u u u(t) in its components toward northern and eastern directions as where , and ψ denotes the direction angle of u u u 0 measured from north toward east.Alternatively, the trajectory can be parameterized as where , and ψ denotes the direction angle of −µ µ µ LS = µ µ µ S − µ µ µ L measured from north toward east.Note that the starred . Source-lens trajectory u u u(t) as seen by the observer, showing the effect of annual parallax.At epoch t0, the tangent to the source-lens trajectory u u u0 is not (anti-)parallel to the direction of the lens-source proper motion μ µ µ LS , but differs by πE γ γ γ(t0) [Equation ( 11)], related to the orbital velocity of the Earth at t0.Consequently, we distinguish the direction angles ψ, ψ , and ϕLS, referring to u u u0, and ∓μ µ µ LS , respectively.Furthermore, tE = 1/| u u u0| and t E = θE/|µ µ µLS|.
quantities, ψ , 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., ψ eq , ψ ecl , ψ gal (and ψ eq , ψ ecl , ψ gal ) for equatorial, ecliptic, and galactic 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, gives the position angle of the proper motion of the lens with respect to the source µ µ µ LS , measured from equatorial north toward east.We illustrate the geometry of the source-lens trajectory u u u(t) in Figure 1.

Measuring the Lens Mass
As stated above, the only useful physical parameter in a typical microlensing event is the timescale t E = θ E /µ 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, π 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 π E ≡ π LS /θ E can be inferred.Coincidently, 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 θ E , so that by combining photometric and astrometric observations, M lens , π LS , and µ µ µ LS become fully decoupled.Specifically, with π E from the photometry and θ E from the astrometry, the definition of θ E , Equation (1), immediately gives us where κ = 4G/[c 2 (1 AU)] 8.144 mas M −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, π LS = π E θ E , can then be used to estimate the distance to the lens, using Equation (2).As a bonus, the event timescale gives a direct 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.

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 , 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 θ E .In this case, the size of the angular Einstein ring is θ E 4 mas, so that the maximum astrometric shift is ∼1.4 mas (occurring at a separation of u = √ 2), and the maximum light magnification is a factor of ∼20 (at closest angular approach).
As Figure 2 illustrates, and as discussed in §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 θ 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 fig- ure 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.
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-SNR HST imaging, as discussed below.Therefore the deflection is detectable at 1σ per epoch for lens masses down to ∼ 0.5 M , 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.

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 capabil-Figure 2. Astrometric shift (top panel) and light magnification (bottom panel) for a microlensing event produced by a 5 M black hole at a distance of 2 kpc passing in front of a background star at 8 kpc.The assumed minimum impact parameter is u0 = 0.05.The maximum astrometric shift of the source is ∼1.4 mas, at u = √ 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.
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 = √ 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.
ity to carry out sub-milliarcsecond astrometry through a variety of techniques.For example, high-S/N HST observations of isolated sources were used to achieve sub-milliarcsecond 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 microarcseconds, in order to measure the transverse velocity of M31 (Sohn et al. 2012).Spatialscan techniques have been used to achieve an astrometric accuracy of ∼30 microarcseconds (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 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 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 to 0.68 mas.The timescales of those events were 60 to 160 days, and there were no detections of astrometric deflections.In 2009, we began a multi-cycle 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.
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 neighbor- ing 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.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 groundbased survey observations, and more typical seeing, re-spectively.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.

HST OBSERVATIONS
The  event satisfies all the selection criteria for our HST follow-up program described in §3.1, and thus we triggered our observing sequence.Our project had a "non-disruptive" target-ofopportunity status, requiring a lead time of about two to three 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 0.08), so that the expected astrometric deflection was δ 0.04 θ 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 §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 six years, 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 −1 .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 timedependent 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 as-trometric 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.

Image Processing
We used the flat-fielded and CTE-corrected ( flc) images produced by the Space Telescope Science Institute pipeline reductions (Sahu et al. 2021;Dressel et al. 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).

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 website2 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-year course of our observations due 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  at the center of the reference frame at (x, y) = (1000, 1000) at the 2016.0epoch, 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 that we derive represent absolute proper motions.After this initial set-up 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 ≥ 18.The agreement is imperfect, of course, since the HST observations have a 6-year baseline and have higher S/N in individual measurements, resulting in higher accuracy in proper-motion measurements, particularly at fainter magnitudes. Te agreement is better for the brighter sample, for which Gaia proper-motion errors are typically 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 years compared to the 3 years 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 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  .Proper motions of the reference stars used in our analysis derived from our HST observations with a 6year 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.
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 referencestar 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 the histograms of the residuals of each measurement from that star's proper motion solution along the RA and Dec 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 pixel (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 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 dot.The nearby bright star 0. 4 away and about 20 times brighter than the unmagnified source is shown as a red dot.Magenta points 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 dot).Bottom right panel: proper motions of the bulge-like stars (black crosses), the reference stars (magenta), and the source (green dot).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.Proper-motion errors of the reference stars.The y-axis shows the proper-motion error, defined as σµ = [σ 2 µx + σ 2 µy ] 0.5 , where σµx and σµy are the proper-motion errors along the x and y axes, which are parallel to right ascension and declination, 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.have similar parallactic motion as the source star.We note that the source parallax is small to begin with ( 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 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 wellsampled, 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 pixel (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 crucial in improving the analysis/results.The resulting astrometric positions of the source were used for further analysis as described in the next section.

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  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.
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.So, 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 §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 years of observations with HST.There is also no detectable blending as described in more detail in §8.2.The source at baseline brightness has apparent magnitudes of m F606W = 21.946± 0.014 and m F814W = 19.581± 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 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.

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 highmagnification event, it attracted intensive monitoring by a number of additional ground-based telescopesespecially 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 re-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 seven days around peak magnification (bottom panel).Superposed is our model fit to the light curve, from the analysis described in the next two sections.

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.

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 which 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.
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 to 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.
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 keep it as a variable in our analysis.The final values ( §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 difficult to include such effects deterministically because of the imperfect knowl-egde of the blending (at the sub-percent level) for individual images.

Trajectory
As discussed in §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 π E and ϕ 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, ϕ LS corresponds to the position angle (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 §2.2).
In principle, a sufficiently accurate light curve can provide good constraints on both π E and ϕ 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 π E and ϕ LS are strongly correlated.In addition, the derived value of π 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 π E ranges from 0.07 to 0.12, with larger values corresponding to larger values of ϕ LS , ranging from 330 • to 358 • .The reason is that increasing ϕ 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 π E has to increase with ϕ 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 π E and ϕ LS can be derived.Fortunately, astrometry provides a robust independent estimate of ϕ LS , allowing us to break this degeneracy and determine the two quantities separately.

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 θ E , with u RA increasing to the east.We have used our actual final model described below in §8.2 for this illustration.The top panel shows the motion of the lens, with π E = 0.0894 and ϕ LS = 342.• 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, com- puted using the JPL ephemeris of the Earth. 3Red dots 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 §8.3.)The plot shows that near the closest angular approach-and thus the peak magnificationthe relative path is substantially affected by parallax; however, the astrometric deflection is very small at this time (see §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 3 https://ssd.jpl.nasa.gov/horizons/app.html#/ the late-time deflection directions robustly constrain the orientation, ϕ LS , of the lens trajectory.

Constraining the Lens Trajectory Orientation
Since the parallactic effect is unimportant for constraining ϕ LS , we first fitted the photometry using a light-curve model that neglects parallax.The resultant model (with t 0 − 2450000 = 5763.33,t E = 231.56days, A max = 372.62,and g = 19.3)predicts the total deflection in units of θ E as a function of time, through Equation ( 6).
As described in §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 both by the proper motion of the source, and its deflections, which are a function of θ E and ϕ LS .We fitted a model to the positions, whose parameters are the x (RA) and y (Dec) components of the proper motion, θ E , and ϕ LS .This fit resulted in values of θ E = 5.2 ± 0.5 mas and ϕ LS = 337.• 9 ± 5. • 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 RA and Dec 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 RA versus Dec 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.The solid black line passes through the origin, at the ϕ 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 ϕ LS = 354.• 8 (shown by the solid magenta line), the allowed range of ϕ 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.

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.In addition to OGLE and MOA, photometric timeseries data were obtained at several more observatories (see Table 2).These data typically cover a relatively narrow range of ∼20 days around the peak, beíng 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 signif-icant 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 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.uRA and uDec are the RA and Dec components of the lens-source separation, in units of θE.Bottom panel: enlarged view of the lens trajectory respect to the source around closest angular approach.Black lines correspond to the u0,− solution, which is the preferred solution, and the red lines show the less-preferred u0,+ 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 u0,+ solution (shown by the solid red line) passes on the north side of the source, but the trajectory of the u0,− solution (shown by the black solid line) passes to the south.The u0,− solution is used in our analysis presented here.It is worth noting, however, that the derived θE, and hence the mass of the lens, is nearly identical for both solutions.ϕLS is the position angle of the barycentric proper motion of the lens with respect to the source (i.e., no parallax), as shown here.
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 SIGNALMEN microlensing anomaly detector (Dominik et al. 2007(Dominik et al. , 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 where σ denotes the reported error bar, κ is a scaling factor, and σ 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 κ and σ 0 in Table 5.If the size of the error bars does not vary substantially, there is a degeneracy between κ and σ 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,− .The microlensing parallax parameter π E is constrained by the wing of the light curve and much less sensitive to the peak region, and therefore is mostly constrained by the microlensing survey data.From various combinations of data sets, we consistently find π E = 0.10 ± 0.02, but some variation in the trajectory angle ψ, correlated with t E and the blend fraction, yielding visually indistinguishable model light curves.
However, the angle of lens-source proper motion, ϕ LS , follows robustly from the astrometric data (see §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 semi-major axis is parallel to µ µ µ LS .If we restrict this angle to the range 333 • ≤ ϕ LS ≤ 343 • , as suggested by the astrometric data, the photometric light curve does not change substantially, and we find π E 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 ϕ LS and π E in the photometric solution.However, restricting the trajectory angle ϕ LS as robustly derived in the last section from the orientation of the centroid shifts to the range 333  The combined constraint from photometry and astrometry is used in our subsequent analysis.
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 θ E and π 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-the-sky approach described in §7.2, which makes it easier to work with, and also show the actual paths of the lens and the source, and 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 §9.2, the best constraints on the source distance come instead from photometry and high-resolution spectroscopy.
From these parameters, the undeflected paths of the source and lens can 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 §8.1).Also, in order to avoid undue impact from any secular variations in photometric responses, we only include OGLE and MOA photometric measurements within ±2 years from the peak of the event.We adopt the approximate values of ϕ 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 dispersion in the measurements, and validated each solution based on a reasonable match to the astrometric data.In our final model, most photometric points are still at an uncertainty below 10 mmag.This solution has a total astrometric χ 2 of 136 with 106 points.(The solution with nominal weights has a higher astrometric χ 2 of 149.)The parameters of the final model are given in Table 4.
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 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 as well as 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 fit the HST photometry.The resulting fit is excellent (see Figure 13), and yields baseline source magnitudes of m F606W = 21.946± 0.012 and m F814W = 19.581± 0.012, with corresponding blending factors for the HST photometry of g = −0.012± 0.015 and −0.006 ± 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 a Notes: (1) Undeflected proper motion of the source in Gaia EDR3 absolute frame; (2) Angular Einstein radius; (3) Angular Einstein radius θE divided by absolute value of lens-source proper motion µ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 θE; (7) Impact parameter in units of θE without parallax motion; (8) Impact parameter derived using the model-fit parameters above and after including parallax motion, in units of θE; (9) Time of closest angular approach derived as above and after including parallax motion.
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.

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 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 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 4 give the proper motions of the source.We use the values of θ E , t E , and ϕ LS given in the next three rows to determine the proper motion of the lens with respect to the source as −2.10 ± 0.22 and 6.66 ± 0.67 mas yr −1 in RA and Dec, respectively.The resulting absolute proper motions of the lens are given below in Table 6.

PROPERTIES OF THE SOURCE
As described in § §2.2-2.3, the mass determination for the lens does not depend upon the individual distances to the lens and source, but only on the relative lenssource parallax, π LS , and the Einstein ring radius, θ Equantities 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 and lens on sky 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 0.00271.Thus it is desirable to estimate the angular diameter of the source, to verify that it is consistent with the pointsource light-curve modeling adopted in the previous section.
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.

High-Resolution Spectroscopy
In a target-of-opportunity program focusing on highdispersion 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 modelatmosphere analysis of the spectra, obtaining an effective temperature and surface gravity of T eff = 5382 ± 92 K and log g = 3.80 ± 0.13, consistent with a late Gtype subgiant.The metallicity was found to be slightly above solar, at [Fe/H] = +0.26± 0.14, and the radial velocity is +134.0km 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 Figure 18.Predicted and measured positions of the source for our adopted joint astrometric and photometric fit.Small dots 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.
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.

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 tool4 (Bressan et al. 2012;Marigo et al. 2017); and "BaSTI" isochrones5 (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 Note-The adopted error bar (in magnitude) becomes σ = (κ σ) 2 + σ 2 0 , where σ is the reported error bar of the photometric data.We have applied range constraints κ > 0.1 and σ0 ≥ 10 −5 .An asterisk ( * ) indicates that the value is at the range boundary.
(red filled circle with error bars).Superposed are PAR-SEC 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 +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 λ /A V given for the WFC3 filters at the PARSEC website (derived from the interstellar extinction curves of Cardelli et al. 1989 andO'Donnell 1994, with R V = 3.1), this color excess cor-responds to an extinction of A F814W = 3.33 mag.Thus, using an apparent magnitude at baseline of m F814W = 19.581from Table 1, we find a true distance modulus of (m − M ) 0 = 13.86,giving a linear distance of D S = 5.9 ± 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 2.5, as deduced by Nataf et al. (2013) for the Galactic bulge region, were used, the distance would increase by ∼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 ≈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 ∼0.0038 mas.
Location of the MOA-11-191/OGLE-11-462 source star in the log g versus 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 to 10 Gyr.The orange curves are BaSTI isochrones for [M/H] = +0.3 and the same range of ages.

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 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 point-like source.
Averaging the latter over an extent δ explicitly gives For u 1, one finds A(u) 1/u and thereby A (u) 2/u 3 .Consequently, we find the largest differ-ences between finite and point-like source for the smallest u.For a source of angular radius θ , we find a source size parameter ρ ≡ θ /θ E , and if u 0 ρ , we can approximate the magnification difference by With the angular diameter of the source 2 θ 0.0038 mas and the angular Einstein radius θ E = 5.2 mas, we find ρ ≡ θ /θ E = 0.00037 while u 0 = 0.00271.Inserting into Equation ( 20) gives ∆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 20, comparing 2.5 log[(A + g)/(1 + g)] for A and A + ∆A gives a difference of 4 mmag between finite source and point-like 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 π 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.

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.
An object with a mass this large cannot be a single (or double) NS or white dwarf.It can only be a BHor 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.

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 π LS = 0.463 ± 0.051 mas (from §8), and the distance to the source (from §9.2), D S = 5.9 ± 1.3 kpc.
These values yield a lens distance of D L = 1.58 ± 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.

From the Baseline Magnitude
Main-sequence stars with a mass in the range 7.1 ± 1.3 M have a range of I-band absolute magnitudes of M I = −1.5 ± 0.3, based on the tabulation 6 of stellar properties compiled by Pecaut & Mamajek (2013).
The lens distance of 1.58 ± 0.18 kpc, from the previous subsection, corresponds to a true distance modulus of (m − M ) 0 = 11.0 ± 0.2.Thus a main-sequence lens of the measured mass would have an unreddened apparent magnitude of I 9.5.The I-band extinction of the source (from §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 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.

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 6 Online version of 2021 March 2, at http://www.pas.rochester.edu/ ∼ emamajek/EEM dwarf UBVIJHK colors Teff.txt θ E = 5.18 ± 0.51 mas and t E = 270.7 ± 11.2 days imply that the proper motion of the lens with respect to the source is 6.99 ± 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 years 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 7 from the finalepoch F814W images (in the individual un-resampled 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 PSF to optimize the fit (again, in the un-resampled 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% 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 abso-7 From https://www.stsci.edu/hst/instrumentation/wfc3/data-analysis/psf lute magnitude is M I 10.2 (corresponding to a mainsequence star of about 0.2 M or less, or a white dwarf of ≈0.6 M with a cooling age older than ∼10 7 years).
These considerations leave a single BH-or a BH with a mass greater than ∼3.55 M in a binary paired with another BH, or with a NS or faint star-as the only viable possibilities for the lens.

Constraints on Binarity
Having established that the MOA-11-191/OGLE-11-462 lens is an extremely faint, or non-luminous, object with a mass of 7.1 ± 1.3 M , 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, θ 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 θ E , the caustic extends to (Schneider & Weiss 1986;Erdl & Schneider 1993;Dominik 1999) Comparing this to u 0 0.00271 yields the constraint d < 0.07 for the equal-mass case, with less than 11 % variation over the range q ∈ [0.4,1].With D L θ 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 0.00271 and blend ratio g 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 ± 0.55 M , the merging timescale for such a system is 10 7 years (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 θ E referring to the mass of the object identified, in the wide-binary case one finds the caustic extending to (e.g., Dominik 1999) 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 Observatory, Kumeu Observatory, Vintage Lane Observatory, the Wise Observatory 0.46m, and the Weizmann Institute.Setting the requirement that the χ 2 due to the planet is >50, we find that the sensitivity limits can be approximated with a piecewise function for log q > −0.4 0.3413 log q + 1.677 for − 4.8 < log q ≤ −0.4 0.04 for log q ≤ −4.8 . .

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 ± 1.3 M .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.

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.

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.

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 to 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 BP − G RP 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 −1 (Table 6), it would take only 0.2 Myr for the lens to traverse 20 , so it's possible that the lens originated much farther away.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 Galactic north 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 or +100 km s −1 changes V and W only by about ±0.2 and ±2.8 km s −1 , 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 the older population with red points.The older population, which has proper-motion dispersions of σ µ,l = 2.59 mas yr −1 and σ µ,b = 1.59 mas yr −1 , is dynamically hotter than the younger population, which has dispersions of σ µ,l = 1.08 mas yr −1 and σ µ,b = 1.05 mas yr −1 .
The BH itself is a prominent outlier, predominantly in b, relative to the older and younger populations (∼4 and 6σ, respectively).This suggests 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 −1 , corresponding to a tangentialvelocity offset of ∼45 km s −1 .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  12. DISCUSSION AND FUTURE WORK

The Formation of Isolated Black Holes
As reviewed in §1.1, it is generally believed that only stars above ∼20 M collapse to form BHs. To form an isolated ∼7 M 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. 1999), 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 a SN, and its mass loss and SN kick detached the binary.for a sample of stars at similar distances to that of the black hole, selected as described in the text.Typical uncertainties are about ±0.05 mas yr −1 in both coordinates.Blue points represent a younger main-sequence population with Gaia colors of GBP −GRP < 1.0, and red points 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 black hole, which is a significant outlier compared to the surrounding population.
The MOA-11-191/OGLE-11-462 BH's inferred mass of ∼7.1 M lies remarkably close to the peak of the Galactic X-ray-binary BH mass distribution (see §1.1).This may suggest a similar evolutionary path with mass loss through close-binary interactions, but a larger sample of mass measurements for isolated BHs is needed to distinguish these different formation scenarios.

Future Work
There is no known X-ray or radio source at the position of  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 §1.2, and these would be worth attempting.
Finally, we note the potential of future facilities, such as the Very Large Telescope Interferometer (VLTI), 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 square-degree region in the 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 wideangle 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. 2018a).In the bulge, it will complement the Roman data by filling in gaps in the Roman survey cadence (Street et al. 2018b;Strader et al. 2018).Rubin discoveries may need high-resolution timeseries imaging follow-up for full characterization of the events.
Note: 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.
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,8 under proposal IDs 12322, 12670, 12986, 13458, 14783 and 15318.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 Astrometry of Long-Duration Microlensing Events

Figure 4 .
Figure 4. HST image in the F814W (I-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 WNW, 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.

Figure 5 .
Figure 5. 2. 1 × 2. 0 cutouts around MOA-11-191/OGLE-11-462 as observed by HST at all 8 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.

Figure 6
Figure6.Proper motions of the reference stars used in our analysis derived from our HST observations with a 6year 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.

Figure 7 .
Figure 7. Left panel: color-magnitude diagram [mF606W versus (mF606W − mF814W)] for all stars in the HST/WFC3 field.The main-sequence turnoff occurs at mF606W22, 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 dot.The nearby bright star 0. 4 away and about 20 times brighter than the unmagnified source is shown as a red dot.Magenta points 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 dot).Bottom right panel: proper motions of the bulge-like stars (black crosses), the reference stars (magenta), and the source (green dot).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.

Figure 8 .
Figure 8. Motions of nine representative astrometric reference stars in right ascension and declination.The red lines are linear fits to the proper motions of the stars.
Figure 9.Proper-motion errors of the reference stars.The y-axis shows the proper-motion error, defined as σµ = [σ 2 µx + σ 2 µy ] 0.5 , where σµx and σµy are the proper-motion errors along the x and y axes, which are parallel to right ascension and declination, 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.

Figure 10 .
Figure10.Histogram of the position residuals for all the reference stars in RA and Dec 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 Figure9, reference stars cover a range of magnitudes around the unamplified source magnitude.

Figure 11 .
Figure11.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 Figure12, 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.

Figure 13 .
Figure 13.Top panel: Photometry of MOA-11-191/OGLE-11-462, obtained with HST over 6 years at 8 epochs in the F606W (V ) and F814W (I) filters (black filled circles), along with our final model fits from §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-year 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 §8.2 for more details and discussion).

Figure 14 .
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 upper panel and over a zoomed-in region covering seven days around peak magnification in bottom panel.All data have been transformed to OGLE I magnitudes according to the inferred baseline magnitudes and blend ratio from the common model.

8. 1 .
First Approach: All Photometric Data Sets, and Robustness of Parallax Measurement

Figure 16 .
Figure16.The top two panels show the average values of the measured deflections in RA and Dec at each HST epoch.We use the values of t0, tE, and u0 as derived from a light-curve fit without parallax, and fit for the proper motion of the source, θE, and ϕLS.The solid black line is the best fit with θE = 5.2 mas and ϕLS = 337.• 9.The bottom left panel shows the total deflections at each epoch.The bottom right panel shows the RA versus Dec 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 ϕLS angle derived above.The dotted blue lines indicate the allowed range from astrometry based on the uncertainties.The allowed range of ϕLS from photometry alone is shown by the dotted magenta lines, the best fit value being ϕLS = 354.• 8 shown by the solid magenta line.The combined constraint from photometry and astrometry is used in our subsequent analysis.

Figure 17 .
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 t0.The small dots 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.

Figure 20 .
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, then the pixels are stacked.(D) The original minus the optimized PSF.(E) The previous panel with a neighbor added to the previous frame 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.

Figure 21 .
Figure21.Color-magnitude diagram for a sample of stars in the Galactic disk at similar distances to that of the black hole, 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.

Figure 22 .
Figure22.Proper-motion components from Gaia EDR3 for a sample of stars at similar distances to that of the black hole, selected as described in the text.Typical uncertainties are about ±0.05 mas yr −1 in both coordinates.Blue points represent a younger main-sequence population with Gaia colors of GBP −GRP < 1.0, and red points 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 black hole, which is a significant outlier compared to the surrounding population.

Table 2 .
Journal of HST Wide Field Camera 3 Observations a Individual exposure times ranged from a minimum of 60 s at Epoch 1, to a maximum of 285 s at later epochs.

Table 4 .
Parameters of the Full Fit to Astrometry and Photometry

Table 5 .
Revised Error Bars of Photometric Data

Table 6 .
Properties of the(1) This paper, §10.1; (2) This paper, §10.2; (3) This paper, Table4;(4) This paper, derived from model in Table4; 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 due to unknown actual radial velocity.11.3.Space Velocity, Population Constraints, Natal Kick