Adaptive Optics Imaging Can Break the Central Caustic Cusp Approach Degeneracy in High-magnification Microlensing Events Sean K. Terry1 , David P. Bennett2,3 , Aparna Bhattacharya2,3, Naoki Koshimoto2,3 , Jean-Philippe Beaulieu4 , Joshua W. Blackman4 , Ian A. Bond5 , Andrew A. Cole4 , Jessica R. Lu1 , Jean Baptiste Marquette6,7 , Clément Ranc8 , Natalia Rektsini4, and Aikaterini Vandorou2,3 1 Department of Astronomy, University of California Berkeley, Berkeley, CA 94701, USA; sean.terry@berkeley.edu 2 Code 667, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 3 Department of Astronomy, University of Maryland, College Park, MD 20742, USA 4 School of Natural Sciences, University of Tasmania, Private Bag 37, Hobart, TAS, 70001, Australia 5 School of Mathematical and Computational Sciences, Massey University, Auckland 0745, New Zealand 6 Laboratoire d’Astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, allée Geoffroy Saint-Hilaire, F-33615 Pessac, France 7 (Associated with) Sorbonne Université, UPMC et CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis Bd Arago, F-75014 Paris, France 8 Zentrum für Astronomie der Universität Heidelberg, Astronomiches Rechen-Institut, Mönchhofstr. 12–14, D-69120 Heidelberg, Germany Received 2022 June 5; revised 2022 September 13; accepted 2022 September 26; published 2022 October 26 Abstract We report new results for the gravitational microlensing target OGLE-2011-BLG-0950 from adaptive optics images using the Keck Observatory. The original analysis by Choi et al. and reanalysis by Suzuki et al. report degenerate solutions between planetary and stellar binary lens systems. This particular case is the most important type of degeneracy for exoplanet demographics because the distinction between a planetary mass or stellar binary companion has direct consequences for microlensing exoplanet statistics. The 8 and 10 yr baselines allow us to directly measure a relative proper motion of 4.20± 0.21 mas yr−1, confirming the detection of the lens star system and ruling out the planetary companion models that predict a ∼4× smaller relative proper motion. The Keck data also rule out the wide stellar binary solution unless one of the components is a stellar remnant. The combination of the lens brightness and close stellar binary light-curve parameters yields primary and secondary star masses of M 1.12A 0.09 0.11= - + and M M0.47B 0.10 0.13 ☉= - + at a distance of D 6.70L 0.30 0.55= - + kpc and a projected separation of 0.39 0.04 0.05 - + au. Assuming that the predicted proper motions are measurably different, the high-resolution imaging method described here can be used to disentangle this degeneracy for events observed by the Roman exoplanet microlensing survey using Roman images taken near the beginning or end of the survey. Unified Astronomy Thesaurus concepts: Gravitational microlensing (672); Binary lens microlensing (2136); High- resolution microlensing event imaging (2138); Binary stars (154); Galactic bulge (2041) 1. Introduction Gravitational microlensing enables the detection of stars and exoplanets in a wide range of environments and distances along the line between Earth and the central region of the Galaxy. In addition to main-sequence stars and exoplanets, more exotic systems like a Jupiter analog orbiting a white dwarf (WD; Blackman et al. 2021) and an isolated black hole or neutron star (NS; Lam et al. 2022; Sahu et al. 2022) have recently been published. Koshimoto et al. (2021b) recently used exoplanet microlensing detections to show that there is no significant dependence on planet frequency with galactocentric distance. This implies that planets residing in the Galactic bulge are likely similar to planets near the solar system. While the number of detected microlensing events per year has been steadily rising, a unique circumstance of binary lens microlensing events is that they can possess different types of degeneracies. The well-known “close–wide” degeneracy occurs when the central caustic shape between a closely separated binary lens and either of the two caustic shapes in a widely separated binary lens are essentially identical (Albrow et al. 2001). There are degeneracies involving the binary mass ratio and finite size of the source star for low- magnification events, as well as degeneracies between planetary caustic perturbations and extreme flux ratio binary events (Gaudi 1998). A recent paper (Zhang et al. 2022) describes the “offset” degeneracy, which combines the close– wide degeneracy with other degeneracies relating models with different lens separations. Fortunately, many of these degeneracies can be mitigated by obtaining high-accuracy and well-sampled photometry during the microlensing light curve, particularly during caustic crossings or close approaches. Many of these degeneracies involve minor differences in the lens separation in the plane of the sky that are smaller than the uncertainty due to the unmeasured uncertainty along the line of sight, but in some cases, the degeneracy can involve very different separations. In the 30- event exoplanet demographics study of Suzuki et al. (2016, hereafter S16), there were eight planetary systems with lens separation degeneracies that were too small to have an important effect on the demographics (Gould et al. 2006; Bennett 2008; Janczak et al. 2010; Bachelet et al. 2012; Suzuki et al. 2013; Bennett et al. 2014, 2018; Nagakane et al. 2017), two planetary events with strong close–wide degen- eracies (Dong et al. 2009a; Fukui et al. 2015), and the single event presented in this paper, where the primary degeneracy is between planetary and stellar binary models (Han & Gaudi 2008). The Astronomical Journal, 164:217 (14pp), 2022 November https://doi.org/10.3847/1538-3881/ac9518 © 2022. The Author(s). Published by the American Astronomical Society. Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. 1 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0002-9881-4760 https://orcid.org/0000-0002-9881-4760 https://orcid.org/0000-0002-9881-4760 mailto:sean.terry@berkeley.edu http://astrothesaurus.org/uat/672 http://astrothesaurus.org/uat/2136 http://astrothesaurus.org/uat/2138 http://astrothesaurus.org/uat/2138 http://astrothesaurus.org/uat/154 http://astrothesaurus.org/uat/2041 https://doi.org/10.3847/1538-3881/ac9518 https://crossmark.crossref.org/dialog/?doi=10.3847/1538-3881/ac9518&domain=pdf&date_stamp=2022-10-26 https://crossmark.crossref.org/dialog/?doi=10.3847/1538-3881/ac9518&domain=pdf&date_stamp=2022-10-26 http://creativecommons.org/licenses/by/4.0/ Choi et al. (2012, hereafter C12) gave further descriptions of this degeneracy, which is related to the source star trajectory approaching a central caustic cusp due to either a planetary or stellar companion to the host. In the planetary case, the source star with a large source radius crossing time (t*) passes by the two strong cusps that bracket a negative perturbation region at a trajectory angle α∼ 90° with respect to the planet–star axis. The effect of the weaker cusp between these two cusps is not obvious in the light curve. For the stellar binary case, the source star, with a much smaller t*, passes two adjacent cusps that bracket a weaker negative perturbation region caused by a diamond-shaped caustic at a trajectory angle α∼ 45° with respect to the binary axis. This degeneracy is severe because regardless of how well the perturbation is sampled by the data, the interpretation of the planetary or binary solution is generally limited by the systematics of the photometry. In contrast to the degeneracies mentioned previously, this degeneracy is likely the most important of its kind when it comes to exoplanet demographics because the ambiguity of the lens system parameters from the light-curve modeling have a direct consequence on the exoplanet statistics that are drawn from the data (i.e., planetary systems versus nonplanetary systems). Finally, we refer to this degeneracy in a more specific manner than C12; we denote it the “central caustic cusp approach” degeneracy. Since the source trajectory in both models approaches a caustic, finite source effects can be observed that allow the source radius crossing time, t*, to be measured. With knowledge of t*, an estimate of the lens– source relative proper motion, μrel, can be made for each model. Further, high-resolution follow-up observations can be made years after the event to directly measure μrel (Batista et al. 2015; Bennett et al.2015; Bhattacharya et al. 2018; Terry et al. 2021) and compare the direct measurement to the light-curve models to determine which interpretation is correct. As we detail in Section 2.1.1, C12 recognized the existence of systematic errors in the photometry from several data sets, which led them to conclude that a Δχ2∼ 105 between the planetary and stellar binary solutions was not significant. Object OGLE-2011-BLG-0950 is included in S16 and Suzuki et al. (2018), which is one of the largest statistical studies of the microlensing exoplanet population. The event was reanalyzed by S16 using optimized photometry from several data sets, which largely resolved the systematic photometry error problem that was present in the C12 analysis. This reanalysis resulted in a much larger uncertainty in the measurement of t* for the stellar binary models. The new modeling work of S16 and the current study show that the models with smaller and uncertain t* values are favored mainly by three data sets (with corresponding Δχ2): Microlensing Observations in Astrophysics (MOA; Δχ2∼ 57), CTIO-I (Δχ2∼ 49), and Danish (Δχ2∼ 63). Further details about the differences between C12, S16, and our new analysis are given in Section 2.1. As mentioned previously, an ambiguous event like OGLE- 2011-BLG-0950 is particularly important for population statistics because ignoring or accepting a target like this could bias results in cases like S16 that aimed to measure the cold exoplanet mass ratio function. We note that another microlen- sing event (OGLE-2011-BLG-0526) exhibiting the same degeneracy was found in the same observing season as OGLE-2011-BLG-0950. It was claimed by C12 that this implies that this degeneracy may be common. Furthermore, a retrospective search through a 9 yr sample (2006–2014)9 of microlensing events from MOA (Bond et al. 2001; Sumi et al. 2003) yields at least three events (MOA-2012-BLG-201, MOA-bin-65, and MOA-2014-BLG-051) that show some evidence of this central caustic cusp approach degeneracy. We note that all of the events included in this 9 yr sample are vetted and classified by eye. Therefore, it may be the case that additional events exhibiting the central caustic cusp approach degeneracy are not identified in this 9 yr sample. Lastly, it is expected that some fraction of the microlensing events that the upcoming Roman Galactic Exoplanet Survey discovers will exhibit this degeneracy. This paper is organized as follows. In Section 2, we perform improved photometry of the light curve and compare our updated best-fit solutions with previous studies of the target. In Section 3, we describe the Keck adaptive optics (AO) follow- up analysis that confirms the stellar binary solution. Section 4 details our lens–source relative proper motion and flux ratio measurements from the 2019 and 2021 epochs. In Sections 5 and 6, we discuss the identification of the lens star and a subsequent search for a luminous lens star companion. We report the lens system physical parameters in Section 7. Lastly, we discuss the overall results and conclude the paper in Section 8. 2. Updated Light-curve Modeling The high-magnification event OGLE-2011-BLG-0950/ MOA-2011-BLG-336, located at R.A.= 17:57:16.63, decl.= −32:39:57.0 and Galactic coordinates (l, b= (−1.93, −4.05)), was alerted by the Optical Gravitational Lensing Experiment (OGLE; Udalski et al. 1993, 2015) on 2011 July 11 and MOA on 2011 July 31. The perturbation was well sampled near the peak of the light curve, and a total of 15 telescopes performed observations at various times throughout the event. Telescopes in New Zealand (Auckland 0.4 m, FCO 0.4 m, Kumeu 0.4 m), Chile (CTIO 1.3 m, Danish 1.54 m), Israel (WISE 1.0 m), Australia (FTS 2.0 m, PST 0.3 m), and Hawaii (FTN 2.0 m) performed follow-up observations around the high-magnifica- tion peak on 2011 August 13. The measurements from each observatory can be seen in Figure 1, with the colored list of telescopes corresponding to the colors of each data point. There have been several improvements to the photometric reduction process since the C12 analysis; therefore, we have rereduced the photometry for several data sets. We have used the updated photometry methods described in Bond et al. (2001, 2017) to reduce the data from the MOA 1.8 m telescope and the SMARTS telescope at CTIO. The SMARTS-CTIO data were reduced using difference imaging photometry (Bond et al. 2001, 2017), and the MOA data were corrected for errors due to chromatic differential refraction (Bennett et al. 2012; Bond et al. 2017). The OGLE data have also been rereduced and included in our new data sets. This rereduction procedure is similar to the reanalysis performed by Suzuki et al. (2016), which we describe further in Section 2.1.2. The updated light-curve modeling follows the image- centered ray shooting method of Bennett & Rhie (1996) and Bennett (2010). The three fundamental microlensing para- meters that are modeled for a single lens are the Einstein radius crossing time, tE, and the time and distance of closest approach between the source and lens center of mass, t0 and u0, 9 https://exoplanetarchive.ipac.caltech.edu/docs/MOAMission.html 2 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. https://exoplanetarchive.ipac.caltech.edu/docs/MOAMission.html https://exoplanetarchive.ipac.caltech.edu/docs/MOAMission.html https://exoplanetarchive.ipac.caltech.edu/docs/MOAMission.html respectively. For binary lenses, there are three additional parameters to model: the binary lens mass ratio, q; their separation, s, in units of the Einstein radius; and the angle between the source star trajectory and the binary lens axis, α. As mentioned earlier, an additional parameter can be modeled if finite source effects are observed; this parameter is the source radius crossing time, t*. The resulting best-fit models show the same fourfold degeneracy that C12 and S16 described. Our best-fit planetary and stellar binary solutions differ by Δχ2∼ 27, with nearly identical χ2 values for s< 1 and s> 1 within both solutions. The degeneracy resulting from our updated modeling is more severe than what C12 found (Δχ2∼ 105) and is in agreement with what S16 found (Δχ2∼ 20). We discuss the differences between our results, the C12 results, and the S16 results in Section 2.1. Figure 1 shows the best-fit stellar binary model (s< 1), and Table 1 shows the parameters of our best-fit close and wide models for both the planetary and stellar binary solutions. Figure 1. The OGLE-2011-BLG-0950 photometry with the updated best-fit light-curve model (second column of Table 1, the stellar binary with s < 1). The middle panel shows an enlarged view of the peak, and the bottom panel shows the residuals to the best-fit stellar binary solution. Table 1 Best-fit Model Parameters Stellar Binary Planetary Parameter s < 1 s > 1 s < 1 s > 1 tE (days) 65.586 ± 0.721 101.870 ± 8.195 68.345 ± 0.817 68.153 ± 0.752 t0 (HJD′) 5786.3965 ± 0.0005 5786.3925 ± 0.0005 5786.3969 ± 0.0005 5786.3959 ± 0.0005 u0 (10 −3) 8.460 ± 0.106 5.443 ± 0.431 8.244 ± 0.110 8.612 ± 0.101 s 0.0768 ± 0017 21.9678 ± 1.2473 0.7257 ± 0.0104 1.3668 ± 0.0191 α (rad) −0.719 ± 0.005 −0.718 ± 0.004 −1.531 ± 0.003 −1.532 ± 0.002 q 0.417 ± 0.115 1.446 ± 0.231 (5.395 ± 0.271) × 10−4 (5.371 ± 0.269) × 10−4 t* (days) 0.0856 ± 0.0882 0.0635 ± 0.0416 0.3136 ± 0.0041 0.3129 ± 0.0037 IS 19.405 ± 0.062 19.397 ± 0.061 19.461 ± 0.069 19.457 ± 0.067 VS 21.144 ± 0.072 21.136 ± 0.072 21.199 ± 0.077 21.195 ± 0.075 Fit χ2 7046.21 7047.68 7020.62 7019.52 3 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. An estimate of the lens–source relative proper motion, μrel, can be made if the angular size of the source can be determined. In order to measure the source radius, we need to determine the extinction-corrected source magnitude and color. To achieve this, the SMARTS-CTIO V- and I-band data were calibrated to the OGLE-III catalog (Szymański et al. 2011), and then we measured the red clump centroid at Vrc− Irc= 2.11, Irc= 15.85, following the method of Bennett et al. (2010). Using the bulge red clump giant magnitude, color, and distance from Nataf et al. (2013), we find I- and V-band extinction of AI= 1.33 and AV= 2.16. This gives an extinction-corrected magnitude of Is0= 18.12± 0.06 and Vs0= 19.00± 0.07. We then use a modified version of the surface brightness relations from Boyajian et al. (2014), using stars spanning the range in colors that are relevant for microlensing targets: V I Ilog 2 0.5014 0.4197 0.2 . 1s s s0 0 0( ) ( ) ( )q = + - -* This yields an angular source size of θ* = 0.93± 0.11 μas for the stellar binary solution and θ* = 0.90± 0.10 μas for the planetary solution. Further, through the relation μrel= θ*/t*, we can determine the lens–source relative proper motion for both stellar binary and planetary interpretations. For the close stellar binary solution, μrel,G= 3.95± 4.10 mas yr−1, and for the wide planetary solution, μrel,G= 1.05± 0.20 mas yr−1, where the subscript “G” refers to the calculation being made in the inertial geocentric reference frame that moves with the Earthʼs velocity at the time of the microlensing event. Figure 2 shows the central caustic for the close planetary and binary models, along with the source trajectory. The source size, ρ* = t*/tE (in θE units), is ∼2.7× larger for the planetary model than for the stellar binary model because the magnifica- tion induced by the planetary model cusps is weaker than in the stellar binary case. Therefore, the source must pass closer to the cusps to get the same signal. However, if the source passes closer to the cusps, this produces sharper light-curve features, unless the t* value is increased to smooth them out. This is a generic feature of the degeneracy. Finally, our new light-curve modeling results are consistent with the results of S16, and our results show smaller best-fit values for the mass ratios and larger tE values than what C12 reported. These differences are carefully examined in Section 2.1. Details of the inferred physical parameters for the lens system are given in Section 7. 2.1. Comparison to Previous Studies 2.1.1. C12 As mentioned previously, C12 performed the original light- curve modeling for this event. This first modeling effort used many data sets derived from earlier iterations of the photo- metric pipelines (OGLE, MOA, CTIO, and other μFUN data sets). Clearly, C12 was aware of systematics in some of those reductions, as evidenced by their conclusion that the Δχ2∼ 105 favoring the planetary solution was not significant, when this large a difference would typically be considered significant enough to rule out less likely solutions. The authors stated that the systematic residuals of the data from the planetary model are larger than the difference between the planetary and binary models. Our analysis using rereduced photometry for several data sets has largely removed these systematic errors (following S16; Section 2.1.2), which results in a significantly smaller Δχ2∼ 27, only slightly favoring the planetary companion solutions. Further investigation of the C12 analysis shows evidence that several of the earlier photometric data sets contributed to their spurious measurement of ρ*. This also led to smaller error bars on ρ* (see Table 1 of C12). We find that the finite source effect is largely unconstrained for the stellar binary models (see t* in Table 1) from our analysis of the rereduced photometry. We also note that although C12 reported measurements of both components of the microlensing parallax, πEE and πEN, the error bars for their estimates are of order ∼100%, which we do not consider a significant detection of πE. We include parallax in our modeling and also do not find a significant measurement Figure 2. Central caustic for the best planetary (left) and binary (right) models. The source size and trajectory are denoted by the green circle and red solid line, respectively. 4 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. of πE; the best-fit parallax values for the stellar binary case give πEE= −0.040± 0.034 and πEN= −0.046± 0.233. Lastly, the best-fit tE that we find for the wide binary model, 101.9 days, is larger than the corresponding tE that C12 reported. For the wide binary models, most of the light curve sees only the effect of one lens. So, the effective tE for the event is reduced by q1 + or 1 q 1+ , depending on which star has a close approach with the source. While we use the same coordinate system for all of the models presented in Table 1, it appears that C12 made a change in the coordinate system for their wide binary model. 2.1.2. S16 The light-curve photometry was reanalyzed by S16 for their statistical study of the cold exoplanet population from MOA events detected between 2007 and 2012. Their analysis included the rereduced data from the observatories listed in the previous section. Using this optimized photometry, S16 was able to remove many of the systematic photometry issues that were present in the C12 analysis, which resulted in many S16 results contradicting the best-fit parameters reported by C12, particularly for the stellar binary solutions. They reported a significantly smaller Δχ2 between the stellar binary and planetary solutions, a Δχ2∼ 20 that favors the planetary solutions. Our reanalysis of the light-curve photometry follows that of the S16 analysis and gives Δχ2∼ 27 between the stellar binary and planetary solutions. While this result is consistent with S16, we note that the target would have been formally classified as a planetary event in the S16 statistical analysis, given our Δχ2> 25. However, an event very close to the S16 selection criteria warranted a careful investigation that included both possibilities in a Bayesian analysis (as S16 conducted). The use of optimized photometry led the S16 analysis to properly conclude a lack of constraint on t* in the stellar binary solutions. As can be seen in Table 3 of S16, the authors reported uncertainties of ∼100% on the t* measurement for the stellar binary solutions. We find similarly large t* uncertainties in our best-fit values for the stellar binary solutions (Table 1). Finally, although the less certain t* values lead to larger errors on the μrel estimates for the binary models, the relative difference between these values and the planetary models t*, as well as the inferred μrel, is large enough (i.e., ∼6σ) such that the direct measurement of μrel with Keck (Section 4) remains unambiguous (Figure 6). 3. Multi-epoch High-resolution Imaging with Keck The target OGLE-2011-BLG-0950 was observed with the NIRC2 instrument on Keck II in the Kshort band (λc= 2.146 μm, hereafter Ks) on 2019 May 27. The target was also observed with the OSIRIS imager on Keck I in the Kprime band (λc= 2.115 μm, hereafter Kp) on 2021 July 14. The 2019 Ks-band data have an average point-spread function (PSF) FWHM of 66.2 mas. The 2021 Kp-band data have an average PSF FWHM of 66.8 mas, very similar to the 2019 Ks- band data. Both epochs used the same tip/tilt guide star of R magnitude ∼15 at a separation of ∼5 5 to the north of the target. Although the 2019 NIRC2 data appear to be of equal or slightly better quality than 2021 OSIRIS, there are minor systematic artifacts on the PSF shape due to imperfect AO correcting on the NIRC2 system. These types of PSF systematics have been successfully modeled in the past on highly blended targets (Terry et al. 2021). We regard the astrometry and photometry results from the 2019 data as reliable because they are consistent with the 2021 OSIRIS data (Section 3.3). Ultimately, both data sets are consistent with the stellar binary interpretation for the lens system. The 2019 data suggested that we had detected the lens star at a proper motion only consistent with the stellar binary models; therefore, we reobserved the target in 2021 to confirm that this star had the appropriate proper motion to be the lens system. For the 2019 Ks-band observations, both the NIRC2 wide and narrow cameras were used. The pixel scales for the wide and narrow cameras are 39.69 and 9.94 mas pixel−1, respec- tively. Generally, the wide camera is used for photometric calibration to Vista Variables in the Via Lactea (VVV; as described below), and the narrow camera is used to make the precise measurement of the lens and source star separations. All of the images were taken using the Keck II laser guide star AO system. For the narrow data, we combined 15 flat-field frames, six dark frames, and 15 sky frames for calibrating the science frames. A total of 21 Ks-band science frames with an integration time of 60 s frame–1 were reduced using the Keck AO Imaging (KAI) data reduction pipeline (Lu 2022) to correct instrumental aberrations and geometric distortion (Ghez et al. 2008; Lu et al. 2008; Yelda et al. 2010; Service et al. 2016). Further, a coadd of four wide camera images were used for photometric calibration to images from the VVV survey (Minniti et al. 2010) following the procedures outlined in Beaulieu et al. (2018). The wide camera images were flat- fielded, dark current–corrected, and stacked using the SWarp software (Bertin 2010). We performed astrometry and photo- metry on the coadded wide camera image using SExtractor (Bertin & Arnouts 1996) and subsequently calibrated the narrow camera images to the wide camera image by matching 80 bright isolated stars in the frames. The uncertainty resulting from this procedure is 0.07 mag. For the 2021 Kp-band data, we combined 40 flat-field frames, 10 dark frames, and 15 sky frames for calibrating our science images. The 24 Kp-band science frames with an integration time of 60 s frame–1 were reduced with the same KAI pipeline. The combined frame can be seen in the upper left panel of Figure 3, which has a PSF FWHM of ∼67 mas. It is worth noting that the astrometric distortion solution for the OSIRIS imager has not yet been made publicly available. In development of the distortion solution (M. Freeman et al. 2022, in preparation), it has been shown that the absolute distortion at the pixels located closest to the centroids of the source and lens are [dx, dy]= (−0.219, 0.097) pixels for the source and [dx, dy]= (−0.236, 0.102) pixels for the lens (M. Freeman 2022, private communication). The pixel scale for the OSIRIS imager is 9.95 mas pixel−1, so the difference in the measured distortion at these two locations translates to 0.17 mas on-sky. This is significantly smaller than the astrometric errors calculated from the DAOPHOT_MCMC+Jackknife analysis (Section 3.1.1, Table 3); thus, we conclude that the 2021 OSIRIS astrometry is not significantly affected by unaccounted-for geometric distortions. 3.1. PSF Fitting In the binary solutions, since the lens and source stars have a predicted separation of 0.49× FWHM and 0.65× FWHM in 2019 and 2021, respectively, we expect the stars to be partially resolved, so it is necessary to use a PSF fitting routine to 5 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. measure both targets separately. Following the methods of Bhattacharya et al. (2018) and Terry et al. (2021), we use a modified version of the DAOPHOT-II package (Stetson 1987), which we call DAOPHOT_MCMC, to run Markov Chain Monte Carlo (MCMC) sampling on the pixel grid encompass- ing the blended targets. Further details of the MCMC routine are given in Terry et al. (2021). The stellar profile is clearly extended in both the NIRC2 and OSIRIS data, as can be seen in Figures 3 and 4. Using DAOPHOT_MCMC to fit a single-star PSF to the blend produces the residual seen in the lower left panel of Figure 3, which shows a strong signal due to the extended flux from the blended star. Rerunning the routine in the two-star fitting mode produces a significantly better fit, as expected, with a χ2 improvement of Δχ2= 2462. The two-star residual is nearly featureless, as can be seen in the lower right panel of Figure 3. Table 2 shows the calibrated magnitudes for the two stars of KSSE= 17.02± 0.05 and KNNW= 16.83± 0.05, where the subscript SSE represents the south–southeast star, and NNW represents the north–northwest star. The uncertainties are derived from the “MCMC+Jackknife method” described in Section 3.1.1. Using the VVV extinction calculator (Gonzalez et al. 2011) and the Nishiyama et al. (2009) extinction law, we find a K-band extinction of AK= 0.20± 0.06. From our reanalysis of the light-curve modeling (Section 2), we measure a source color of (V− I)S= 1.74± 0.09, which leads to an extinction-corrected color of (V− I)S,0= 0.88± 0.09. Further, we use the color–color relations of Kenyon & Hartmann (1995) and the I-band magnitude, IS= 19.405, to predict a source K- band magnitude of KS= 17.16± 0.08. This predicted magni- tude is fainter than both stars detected in the Keck epochs, ∼1σ fainter than the SSE star and ∼5σ fainter than the NNW star. While the predicted K-band source magnitude is roughly consistent with the SSE star being the source, it is not as definitive as the typical result for this procedure because the lens and source usually have measurably different brightnesses (Bennett et al. 2020; Bhattacharya et al. 2021; Terry et al. 2021). Because of this potential ambiguity, we perform an additional analysis using a Galactic model to confirm our tentative identification of the source star (Section 5). 3.1.1. Astrometric and Photometric Errors with MCMC+Jackknife The full details of DAOPHOT_MCMC are given in Terry et al. (2021). Here we outline the new modifications we make that include the collation and propagation of astrometry/ photometry errors through iterative MCMC runs on all jackknife frames. First, we generate the individual jackknife frames with the reduce.jackknife() function inside the KAI pipeline. Then we run the standard PSF function inside DAOPHOT on each jackknife frame to generate an empirical PSF associated with each jackknife frame. The same five reference stars within Figure 3. Upper left panel: coadded sum of 24 60 s OSIRIS Kp-band images from 2021. Upper right panel: close-up of blended source and lens stars. Lower left panel: single-star PSF residual clearly showing a signal from the blend. Lower right panel: two-star PSF residual showing a smooth subtraction. The color bar represents the intensity (counts) in the bottom panel residual images. 6 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. 4″ and 1 mag of the target were used in each jackknife frame to build the PSF model. It is necessary to generate a different empirical PSF for each jackknife frame because the shape of the PSF varies (sometimes significantly) between Keck AO images, and this PSF variation is precisely what we want to capture in the astrometric and photometric errors. We then employ an iterative scheme that runs DAOPHOT_MCMC on each jackknife frame. Finally, the output best-fit values and errors from each MCMC are combined in the jackknife error calculation (i.e., Equation (3) from Bhattacharya et al. 2021). These errors are reported for the astrometry and flux ratios in Tables 3 and 4. 3.2. 2019 NIRC2 Analysis As mentioned previously, the 2019 Ks-band images have a PSF FWHM of ∼66 mas, and the average Strehl ratio (SR) across the 21 science frames is SR ∼ 0.33. In an attempt to minimize the effect of PSF systematics, a careful selection of PSF reference stars was made to build the empirical PSF models for this epoch. In testing the different PSF models, we selected between four and nine reference stars with magnitudes −1.03× more likely to be the lens than the SSE object. For the wide scenario, the NNW object is ∼4× more likely to be the lens than the SSE object. However, we note that for the wide solution, the positions of the stars would be slightly different than the positions displayed in Figure 5. The results of this analysis are consistent with what we find in Section 3; the predicted K-band magnitudes imply that the source star is the SSE object. Finally, it is indeed true that we cannot completely rule out the possibility that the SSE star is the lens. Nevertheless, our final results are little affected by this identification because both candidate stars are of similar brightness; thus, the resulting lens properties will ultimately be similar. Figure 6 shows the posterior probability distribution for the lens–source relative proper motion as directly measured by our Keck high-resolution data shown in red. Included in the figure is the prior probability distribution of relative proper motion derived using a Galactic model as described in Koshimoto et al. (2021a) and shown in blue. The μrel,G estimate from the planetary model is given by the orange region, and the 2σ upper limit for the stellar binary model is given by the black solid line and arrow. Although the best-fit planetary solution gives an unusually small μrel,G= 1.05± 0.20 mas yr−1, it could not be completely ruled out until our direct measurement 8 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. Figure 5. Probability distribution for the north and east components of lens–source relative proper motion (μrel) using the Galactic model from Koshimoto et al. (2021a). The lens positions (NNW and SSE) are plotted in black and given by the relative motion of the two stars detected in the 2021 OSIRIS data. This implies that the NNW star is >3× more likely to be the lens than the SSE star. Figure 6. Probability distribution for μrel derived using a Galactic model from Koshimoto et al. (2021a; blue, where dark shades are the central 68% of the distribution, and light shades are the central 95% of the distribution). The posterior distribution from the direct measurement with Keck is shown in red. The orange hatched region indicates the predicted μrel from the best-fit planetary solution, and the black solid line with arrow shows the 2σ upper limit on μrel from the best-fit stellar binary solution. 9 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. of μrel,G= 4.06± 0.22 mas yr−1 with Keck confirmed the stellar binary solution. 6. A Search for a Wide Lens Companion Since we have confirmed the stellar binary solution (Section 3) and identified the luminous lens star (Section 4), we further investigate the possibility of resolving the companion to the primary lens. Given the best-fit close and wide solutions from the light-curve modeling (Section 2 and Table 1), the 2D projected separation between both lens components is ∼0.1× θE for the close scenario and ∼21.9× θE for the wide scenario. This translates to on-sky separations of 0.06 and 26.05 mas, respectively. Clearly, the projected separation for the close stellar binary scenario is too small to be detected in either epoch, at a separation of ∼0.001× FWHM. However, the projected separation for the wide binary scenario, of ∼0.42× FWHM, should allow the companion to be detectable if it is luminous. We conduct two independent searches for this companion in the following subsections. 6.1. Residual Pixel Grid We first analyze the residual pixel grid after the initial two- star PSF fitting is performed. The 2019 residual shows two oversubtracted regions and one undersubtracted region. The oversubtracted regions are both approximately 4× 4 pixels in size and ∼3.5σ below the mean pixel intensity inside the total PSF radius of 20 pixels. The undersubtracted region is approximately 5× 6 pixels in size and 2.5σ above the mean pixel intensity inside the same PSF radius. The undersubtracted region may indicate the presence of an additional luminous source at this location; however, its separation from the lens is inconsistent with the best wide binary solution from the light- curve modeling. The location is measured at ∼78 mas from the luminous lens component and ∼69 mas from the source, which is 3–4× larger than the expected separation from the wide binary solution. As noted previously, the 2019 data have PSF systematics due to imperfect AO correcting. Similar over- and undersubtracted regions are observed in a majority of similar brightness stars in the central 5″ of the frame, particularly a similar undersubtracted region of similar size to that observed on the blended targets. There are no similarly blended stars of approximately equal magnitude near the target to perform a direct comparison of, so we focus this residual grid analysis on the target stars themselves. This analysis shows that the likely cause of the noise regions we detect in the 2019 residual is due to systematic errors on the PSF. In contrast, the 2021 residual pixel grid is significantly smoother, with at most ∼0.8σ deviation above the mean pixel intensity within the same PSF radius as the 2019 data. Additionally, the empirical PSF models we derive for both the 2019 and 2021 Jackknife frames are uniformly circular. The largest nonuniformity in PSF shape used for any of the Jackknife frames is (FWHMY− FWHMX)= 0.023 pixels. For the close binary interpretation, two stars separated by ∼0.001× FWHM are effectively located at the same pixel location in our data, which means the PSF fitting and extraction accounts for the combined flux of both lens stars. Lastly, we conclude that our search through the residual pixel grids shows little to no definitive evidence for an additional source of flux at the expected distance from the luminous lens component for the wide binary scenario. 6.2. Three-star MCMC Search We conduct a three-star DAOPHOT_MCMC search on the 2021 OSIRIS stacked frame with two constraints. The first constraint is the magnitude of the separation between star 2 (the luminous lens component) and star 3 (the nonluminous or less luminous lens component). The second constraint we impose is the magnitude of the separation between star 1 (the source star) and star 3 (the non-/less luminous lens component). The latter constraint is imposed in order to prevent the MCMC from searching in locations that are disallowed by the wide binary best-fit parameters (Table 1). As mentioned in the previous section, there is evidence of PSF systematics in the 2019 NIRC2 data; therefore, we chose to omit the 2019 epoch from the three- star MCMC analysis. For the 2021 OSIRIS data, we used a separation constraint of 2.65± 0.25 pixels or 26.35± 2.49mas for the star 2/star 3 positions and a separation constraint of 4.45± 1.05 pixels or 44.24± 10.44 mas for the star 1/star 3 positions. The three-star MCMC results give positions of the source and luminous lens component that are in agreement (within 1σ) with our two-star MCMC analysis. The distribution for the position of the non-/less luminous lens component is shown in the left panel of Figure 7 as the white and blue contours. The possible positions for the undetected component are approxi- mately perpendicular to the source-luminous lens component separation vector, as is required to be consistent with the best- fit wide binary solution from the light-curve modeling. The right panel of Figure 7 shows the posterior distribution for the calibrated K-band magnitude of the dark lens component, K 22.2 0.8 1.4= - + . This corresponds to a star 3/star 2 flux ratio of f 6.88 10R 4.92 7.12 3= ´- + - . Lastly, we compare the DAOPHOT_MCMC best-fit χ2 values between this three-star analysis and the two-star analysis (Section 3.1). The two-star solution gives a marginally smaller best-fit χ2, with a difference of Δχ2 ∼ 2. Since the flux ratio between the lens compo- nents is very small in the three-star DAOPHOT_MCMC run, the result is a best-fit χ2 that is nearly identical to the two-star DAOPHOT_MCMC result. Additionally, we can compare the flux ratio distribution for the lens compo- nents as given by the DAOPHOT_MCMC analysis with the mass ratio that is given by the best-fit wide binary solution from the light-curve modeling (Section 2). Using empirical mass–luminosity relations (Henry & McCarthy 1993; Delfosse et al. 2000) with an assumed 0.1 mag uncertainty, we find that nearly all flux ratios given by the DAOPHOT_MCMC posterior distribution are inconsistent with the expected flux ratio of fR ∼ 0.44 that corresponds to a mass ratio of q ∼ 1.45 from the wide binary solution (via light-curve modeling). This allows us to rule out the wide binary scenario if we assume that both lens components are luminous. 6.3. A Nonluminous Primary Lens Although the analyses in Sections 6.1 and 6.2 suggest that the close stellar binary scenario is preferred, we cannot completely rule out the possibi- lity that a primary lens object or widely separated lens companion is below the detection threshold in the Keck data. If the less massive component is “dark” ( i .e., a WD or NS) , then it would remain undetectable 10 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. at any separation from its luminous counterpart. For clarification, in the current section and Section 6.4, we use the term “primary” to describe the (less massive) lens component that the source trajectory comes nearest to (see q > 1 wide solution in Table 1) . Given the best-fit mass ratio from the light-curve modeling, q∼ 1.45, and mass–luminosity relations (Henry & McCarthy 1993) assuming all detected flux in Keck is from the companion star, the companion mass is calculated as M M0.96B 0.09 0.11 ☉= - + . Further, this gives a primary lens mass of M M0.66A 0.08 0.10 ☉= - + . The lens system distance is calculated to be D 4.80L 0.30 0.50= - + kpc for this case. This may be consistent with a WD primary orbited by a main-sequence companion. 6.4. A Nonluminous Companion The alternative scenario to a dark primary lens is a dark companion orbiting a luminous primary lens. Again using the best-fit wide binary mass ratio and mass–luminosity relations assuming all of the measured flux is coming from the primary star, the primary lens mass is estimated to be MA= M1.08 0.09 0.11 ☉- + , orbited by a dark companion of mass M M1.57B 0.11 0.15 ☉= - + . For this scenario, the inferred lens system distance is farther, at D 5.79L 0.35 0.61= - + kpc. The companion in this system may be consistent with a WD or NS; however, it is very unlikely that an NS could form (through Type II supernovae, for example) and maintain a companion (Burrows 1987). We conduct a search of all public X-ray survey catalogs and find no objects in the vicinity of OGLE- 2011-BLG-0950. The nearest unrelated transient object, XTE J1755–324, classified as an X-ray nova (Revnivtsev et al. 1998), is located 42′ from the microlensing target. We note that it is unlikely that an NS in a wide orbit would emit X-rays. Further, the WD companion scenario would give an object that is at or above the Chandrasekhar limit. Such WDs are quite rare and do not remain stable for very long (Hillman et al. 2016), further reducing the likelihood for the wide orbiting WD companion interpretation. To conclude, we find no strong evidence of an additional widely separated lens object in either epoch. For the 2021 data, a three-star MCMC analysis gives a possible wide orbiting dark object with a brightness of K 22.2 0.8 1.4= - + , which is below the detection limit. Given the best-fit mass ratio for the wide binary model, this implies either a WD primary lens with a main- sequence companion or a main-sequence primary lens with a WD/NS companion. Nevertheless, both wide stellar binary scenarios are less preferred than the close stellar binary solution by the DAOPHOT_MCMC best-fit χ2, expected flux ratio, and WD/NS formation scenarios. 7. Lens System Properties As a result of the new direct measurement of the lens–source relative proper motion, we have successfully broken the central caustic cusp approach degeneracy for this event. The original fourfold degeneracy has now become a single degeneracy between the close and wide stellar binary solutions. Further, the results of our search for a luminous lens companion (Section 6) give strong evidence that the close stellar binary interpretation is the correct solution. Working from this point, we use the Keck lens flux and mass–luminosity relations (Henry & McCarthy 1993; Henry et al. 1999; Delfosse et al. 2000) in order to constrain the stellar binary lens distance. Given a Galactic latitude of b= −4°.05 and lens system distance of ∼7 kpc, the lens is likely to be behind most of the interstellar dust that is in the foreground of the source. We can describe the foreground extinction as follows: A e e A 1 1 , 5i L D b h D b h i S, sin sin , L S dust dust ( ) ∣ ( ) ∣ ∣ ( ) ∣= - - - - where i represents the passbands I, V, and K. We assume a dust scale height of hdust= 0.10± 0.02 kpc. Additionally, we can use the θE value inferred from the direct measurement of μrel, along with a mass–distance relation assuming we know the Figure 7. Left: three-star MCMC distributions for the positions of the source star and luminous lens component (black contours), as well as the dark/faint lens component (blue/white contours). Right: posterior probability distribution for the K-band magnitude of the dark/faint lens component. The central 68.3% of the distribution is shown in dark blue, and the remaining central 95.4% of the distribution is in light blue. 11 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. distance to the source (Bennett 2008; Gaudi 2012), M c G D D D D4 , 6L S L S L 2 E 2 ( )q= - where ML is the lens mass, G and c are the gravitational constant and speed of light, and DL and DS are the distance to the lens and source, respectively. Figure 8 shows the measured mass and distance of the binary lens. The red curve represents the constraint from the mass–luminosity relation, with dashed lines representing the error from the Keck lens flux measure- ment. For this close stellar binary case, the empirical mass– luminosity relation was numerically calculated for each star in the lens system, considering that the measured lens flux with Keck is the combination of two luminous stars with a nonnegligible mass ratio q. Additionally, the θE constraint from the direct measurement of μrel from Keck is shown in green, with dashed lines representing the error on the θE measurement. Lastly, we include in Figure 8 the estimated mass and distance values from the largely unconstrained microlensing parallax measurement as given by the stellar binary light-curve modeling. The OGLE-2011-BLG-0950 lens system is located at a distance of ∼6.7 kpc and has a mass ratio of q= 0.42± 0.12 and total mass of M M1.59TOT 0.05 0.08 ☉= - + . The combination of the relatively large error on the mass ratio from the light-curve modeling and the relatively small error on the total mass from the Keck imaging gives a primary mass of M M1.12A 0.09 0.11 ☉= - + and secondary mass of M M0.47B 0.10 0.13 ☉= - + . These masses are consistent with a K dwarf orbiting a star near the top of the main sequence. The 2D projected separation for the close binary is measured to be a 0.39 0.04 0.05=^ - + au for the binary stars. Table 5 gives all of the lens system parameters along with their 2σ ranges for the close stellar binary solution. 8. Discussion and Conclusions Our follow-up high-resolution observations of the microlen- sing target OGLE-2011-BLG-0950 have allowed us to make a Table 5 Close Stellar Binary Lens System Properties Parameter Units Values and Rms 2σ Range Angular Einstein radius (θE) mas 0.76 ± 0.08 0.61–0.92 Geocentric lens–source relative proper motion (μrel,G) mas yr−1 4.06 ± 0.22 3.62–4.50 Primary mass (MA) Me 1.12 0.09 0.11 - + 0.94–1.34 Secondary mass (MB) Me 0.47 0.10 0.13 - + 0.28–0.74 2D separation (a⊥) au 0.39 0.04 0.05 - + 0.31–0.49 Lens distance (DL) kpc 6.70 0.30 0.55 - + 6.10–7.80 Source distance (DS) kpc 9.17 0.45 1.07 - + 8.27–11.31 Figure 8. Close stellar binary mass–distance relation for OGLE-2011-BLG-0950L with constraints from the Keck K-band lens flux measurement in red and the angular Einstein radius from the direct measurement of μrel in green. The individual masses are measured to be M 1.12A 0.09 0.11= - + and M M0.47B 0.10 0.13 ☉= - + . 12 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. direct measurement of flux from the lens star(s), as well as a precise determination of the direction and amplitude of the lens–source relative proper motion. We are able to successfully break the central caustic cusp approach degeneracy by showing that the lens–source relative proper motion directly measured with Keck is only compatible with the stellar binary solutions for the lens system. Further, the probability distribution estimates for the lens–source relative proper motion derived using a Galactic model show very low μrel values for the planetary companion models. This low probability for the planetary μrel values could have been noticed a priori. Ultimately, the subsequent Keck observations have now completely ruled out any planetary companion models. Additionally, an analysis of the 2021 OSIRIS data and their residuals favors the close stellar binary solution; however, we cannot fully rule out the wide binary scenario if the binary includes a stellar remnant. We modified the PSF fitting routine DAOPHOT_MCMC (Terry et al. 2021) to calculate Jackknife (i.e., drop-one frame) errors for astrometry and photometry. We used these DAO- PHOT_MCMC+Jackknife error bars for our final analysis. These pipelines, or something similar, can be used in future analyses of highly blended microlensing targets and will likely form the foundation for the Roman mass measurement method. Roman will collect the precursor or follow-up data that are needed to enable direct lens detections. The maximum time baseline between two Roman epochs will be ∼5 yr. We note that Bennett et al. (2006) and Dong et al. (2009a) measured the lens–source separations (with HST) for the first two planetary microlensing events less than 2 yr after the events; these measurements were recently confirmed by Bennett et al. (2020) and A. Bhattacharya et al. (2022, in preparation). With at least 100× more images during its microlensing survey, Roman should reliably measure lens–source separations for many detected events. The study of S16 includes a reanalysis of OB 110950 using optimized photometry, in which they reported the same fourfold degeneracy found in C12. Importantly, the S16 modeling finds the correct best-fit parameters for the stellar binary models, which we confirm in this work. To avoid biasing the statistical results, S16 included both stellar binary and planetary companion possibilities in a Bayesian analysis; they formally concluded that the planetary companion solution is slightly (Δχ2∼ 20) preferred over the stellar binary solution. Our new results that firmly place OB 110950 in the stellar binary regime would at most have only a minor effect (<4%) on the broken power-law function that was reported in S16. Finally, as described in C12, this degeneracy is severe in the sense that two significantly different-shaped central caustics can lead to indistinguishable light-curve perturbations (C12 reported that their Δχ2∼ 105 is not significant because of systematics in the photometry). The reanalysis of the light curve by S16 and this work shows that the degeneracy becomes more severe with the use of newly rereduced photometry. The fact that C12 identified two events from a single observing season with this degeneracy, in addition to ∼three known events in a 9 yr MOA sample that also show evidence of this degeneracy, implies that this type of degeneracy may be common and will likely be encountered in a nonnegligible fraction of Roman microlensing detections. The techniques described in this work should allow for many of the events exhibiting this degeneracy to be reconciled, given measurably different μrel values for the degenerate solutions. The Keck Telescope observations and data analysis were supported by NASA Keck PI Data Award 80NSSC18K0793, administered by the NASA Exoplanet Science Institute. Data presented herein were obtained at the W. M. Keck Observatory from telescope time allocated to the National Aeronautics and Space Administration through the agency’s scientific partner- ship with the California Institute of Technology and the University of California. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. D.P.B. and A.B. were supported by NASA through grant NASA- 80NSSC18K0274. This work was supported by the University of Tasmania through the UTAS Foundation and the endowed Warren Chair in Astronomy and the ANR COLD-WORLDS (ANR-18-CE31-0002). N.K. was supported by the JSPS overseas research fellowship. Work by C.R. was supported by a research fellowship of the Alexander von Humboldt Foundation. I.A.B. is supported by a Marsden Grant from the Royal Society of New Zealand (MAU-1901). This research was also supported in part by the Australian Government through the Australian Research Council Discovery Program (project No. 200101909) grant awarded to A.A.C. and J.-P.B. Some of this research has made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under the Exoplanet Exploration Program. This work made use of data from the Astro Data Lab at NSF’s OIR Lab, which is operated by the Association of Universities for Research in Astronomy (AURA), Inc., under a cooperative agreement with the National Science Foundation. The authors thank Dr. Matthew Freeman for providing OSIRIS distortion details. S.K.T. thanks Gary Bennett for helpful software discussions and comments on this manuscript. The authors would like to thank the anonymous referee for helpful comments that improved the structure and led to a stronger manuscript. Software: daophot_mcmc (Terry et al. 2021), genulens (Koshimoto & Ranc 2021), KAI (Lu 2022), Matplotlib (Hunter 2007), Numpy (Oliphant 2006). ORCID iDs Sean K. Terry https://orcid.org/0000-0002-5029-3257 David P. Bennett https://orcid.org/0000-0001-8043-8413 Naoki Koshimoto https://orcid.org/0000-0003-2302-9562 Jean-Philippe Beaulieu https://orcid.org/0000-0003- 0014-3354 Joshua W. Blackman https://orcid.org/0000-0001-5860-1157 Ian A. Bond https://orcid.org/0000-0002-8131-8891 Andrew A. Cole https://orcid.org/0000-0003-0303-3855 Jessica R. Lu https://orcid.org/0000-0001-9611-0009 Jean Baptiste Marquette https://orcid.org/0000-0002-7901-7213 Clément Ranc https://orcid.org/0000-0003-2388-4534 Aikaterini Vandorou https://orcid.org/0000-0002-9881-4760 References Albrow, M. D., An, J., Beaulieu, J.-P., et al. 2001, ApJ, 549, 759 Bachelet, E., Shin, I.-G., Han, C., et al. 2012, ApJ, 754, 73 Batista, V., Beaulieu, J. P., Bennett, D. P., et al. 2015, ApJ, 808, 170 Beaulieu, J. P., Batista, V., Bennett, D. P., et al. 2018, AJ, 155, 78 13 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0002-5029-3257 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0001-8043-8413 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-2302-9562 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0003-0014-3354 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0001-5860-1157 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0002-8131-8891 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0003-0303-3855 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0001-9611-0009 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0002-7901-7213 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0003-2388-4534 https://orcid.org/0000-0002-9881-4760 https://orcid.org/0000-0002-9881-4760 https://orcid.org/0000-0002-9881-4760 https://orcid.org/0000-0002-9881-4760 https://orcid.org/0000-0002-9881-4760 https://orcid.org/0000-0002-9881-4760 https://orcid.org/0000-0002-9881-4760 https://orcid.org/0000-0002-9881-4760 https://doi.org/10.1086/319437 https://ui.adsabs.harvard.edu/abs/2001ApJ...549..759A/abstract https://doi.org/10.1088/0004-637X/754/1/73 https://ui.adsabs.harvard.edu/abs/2012ApJ...754...73B/abstract https://doi.org/10.1088/0004-637X/808/2/170 https://ui.adsabs.harvard.edu/abs/2015ApJ...808..170B/abstract https://doi.org/10.3847/1538-3881/aaa293 https://ui.adsabs.harvard.edu/abs/2018AJ....155...78B/abstract Bennett, D. P. 2008, Exoplanets (Berlin: Springer), 47 Bennett, D. P. 2010, ApJ, 716, 1408 Bennett, D. P., Anderson, J., Bond, I. A., Udalski, A., & Gould, A. 2006, ApJL, 647, L171 Bennett, D. P., Batista, V., Bond, I. A., et al. 2014, ApJ, 785, 155 Bennett, D. P., Bhattacharya, A., Anderson, J., et al. 2015, ApJ, 808, 169 Bennett, D. P., Bhattacharya, A., Beaulieu, J.-P., et al. 2020, AJ, 159, 68 Bennett, D. P., & Rhie, S. H. 1996, ApJ, 472, 660 Bennett, D. P., Rhie, S. H., Nikolaev, S., et al. 2010, ApJ, 713, 837 Bennett, D. P., Sumi, T., Bond, I. A., et al. 2012, ApJ, 757, 119 Bennett, D. P., Udalski, A., Bond, I. A., et al. 2018, AJ, 156, 113 Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 Bertin, E. 2010, SWarp: Resampling and Co-adding FITS Images Together, Astrophysics Source Code Library, ascl:1010.068 Bhattacharya, A., Beaulieu, J. P., Bennett, D. P., et al. 2018, AJ, 156, 289 Bhattacharya, A., Bennett, D. P., Beaulieu, J. P., et al. 2021, AJ, 162, 60 Blackman, J. W., Beaulieu, J. P., Bennett, D. P., et al. 2021, Natur, 598, 272 Bond, I. A., Abe, F., Dodd, R. J., et al. 2001, MNRAS, 327, 868 Bond, I. A., Bennett, D. P., Sumi, T., et al. 2017, MNRAS, 469, 2434 Boyajian, T. S., van Belle, G., & von Braun, K. 2014, AJ, 147, 47 Burrows, A. 1987, ApJL, 318, L57 Choi, J.-Y., Shin, I.-G., Han, C., et al. 2012, ApJ, 756, 48 Delfosse, X., Forveille, T., Ségransan, D., et al. 2000, A&A, 364, 217 Dong, S., Bond, I. A., Gould, A., et al. 2009a, ApJ, 698, 1826 Dong, S., Gould, A., Udalski, A., et al. 2009b, ApJ, 695, 970 Fukui, A., Gould, A., Sumi, T., et al. 2015, ApJ, 809, 74 Gaudi, B. S. 1998, ApJ, 506, 533 Gaudi, B. S. 2012, ARA&A, 50, 411 Ghez, A. M., Salim, S., Weinberg, N. N., et al. 2008, ApJ, 689, 1044 Gonzalez, O. A., Rejkuba, M., Minniti, D., et al. 2011, A&A, 534, L14 Gould, A., Udalski, A., An, D., et al. 2006, ApJL, 644, L37 Han, C., & Gaudi, B. S. 2008, ApJ, 689, 53 Henry, T. J., Franz, O. G., Wasserman, L. H., et al. 1999, ApJ, 512, 864 Henry, T. J., & McCarthy, D. W. J. 1993, AJ, 106, 773 Hillman, Y., Prialnik, D., Kovetz, A., & Shara, M. M. 2016, ApJ, 819, 168 Hunter, J. D. 2007, CSE, 9, 90 Janczak, J., Fukui, A., Dong, S., et al. 2010, ApJ, 711, 731 Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 Koshimoto, N., Baba, J., & Bennett, D. P. 2021a, ApJ, 917, 78 Koshimoto, N., Bennett, D. P., Suzuki, D., & Bond, I. A. 2021b, ApJL, 918, L8 Koshimoto, N., & Ranc, C. 2021, koshimoto/genulens v1.1: A Tool for Gravitational Microlensing Events Simulation, Zenodo, doi:10.5281/ zenodo.4898012 Lam, C. Y., Lu, J. R., Udalski, A., et al. 2022, ApJL, 933, L23 Lu, J. 2022, Keck-DataReductionPipelines/KAI: v1.0.0 Release of KAI, Zenodo, doi:10.5281/zenodo.6522913 Lu, J., Ghez, A., Hornstein, S. D., et al. 2008, ApJ, 690, 1463 Minniti, D., Lucas, P. W., Emerson, J. P., et al. 2010, NewA, 15, 433 Nagakane, M., Sumi, T., Koshimoto, N., et al. 2017, AJ, 154, 35 Nataf, D. M., Gould, A., Fouqué, P., et al. 2013, ApJ, 769, 88 Nishiyama, S., Tamura, M., Hatano, H., et al. 2009, ApJ, 696, 1407 Oliphant, T. E. 2006, A Guide to NumPy, Vol. 1 (USA: Trelgol Publishing) Revnivtsev, M., Gilfanov, M., & Churazov, E. 1998, A&A, 339, 483 Sahu, K. C., Anderson, J., Casertano, S., et al. 2022, ApJ, 933, 83 Service, M., Lu, J. R., Campbell, R., et al. 2016, PASP, 128, 095004 Stetson, P. B. 1987, PASP, 99, 191 Sumi, T., Abe, F., Bond, I. A., et al. 2003, ApJ, 591, 204 Suzuki, D., Bennett, D. P., Ida, S., et al. 2018, ApJL, 869, L34 Suzuki, D., Bennett, D. P., Sumi, T., et al. 2016, ApJ, 833, 145 Suzuki, D., Udalski, A., Sumi, T., et al. 2013, ApJ, 780, 123 Szymański, M. K., Udalski, A., Soszyński, I., et al. 2011, AcA, 61, 83 Terry, S. K., Bhattacharya, A., Bennett, D. P., et al. 2021, AJ, 161, 54 Udalski, A., Szymanski, M., Kaluzny, J., et al. 1993, AcA, 43, 289 Udalski, A., Szymański, M., & Szymański, G. 2015, AcA, 65, 1 Yelda, S., Lu, J. R., Ghez, A. M., et al. 2010, ApJ, 725, 331 Zhang, K., Gaudi, B. S., & Bloom, J. S. 2022, NatAs, 6, 782 14 The Astronomical Journal, 164:217 (14pp), 2022 November Terry et al. https://ui.adsabs.harvard.edu/abs/2008exop.book...47B/abstract https://doi.org/10.1088/0004-637X/716/2/1408 https://ui.adsabs.harvard.edu/abs/2010ApJ...716.1408B/abstract https://doi.org/10.1086/507585 https://ui.adsabs.harvard.edu/abs/2006ApJ...647L.171B/abstract https://doi.org/10.1088/0004-637X/785/2/155 https://ui.adsabs.harvard.edu/abs/2014ApJ...785..155B/abstract https://doi.org/10.1088/0004-637X/808/2/169 https://ui.adsabs.harvard.edu/abs/2015ApJ...808..169B/abstract https://doi.org/10.3847/1538-3881/ab6212 https://ui.adsabs.harvard.edu/abs/2020AJ....159...68B/abstract https://doi.org/10.1086/178096 https://ui.adsabs.harvard.edu/abs/1996ApJ...472..660B/abstract https://doi.org/10.1088/0004-637X/713/2/837 https://ui.adsabs.harvard.edu/abs/2010ApJ...713..837B/abstract https://doi.org/10.1088/0004-637X/757/2/119 https://ui.adsabs.harvard.edu/abs/2012ApJ...757..119B/abstract https://doi.org/10.3847/1538-3881/aad59c https://ui.adsabs.harvard.edu/abs/2018AJ....156..113B/abstract https://doi.org/10.1051/aas:1996164 https://ui.adsabs.harvard.edu/abs/1996A&AS..117..393B/abstract http://arxiv.org/abs/1010.068 https://doi.org/10.3847/1538-3881/aaed46 https://ui.adsabs.harvard.edu/abs/2018AJ....156..289B/abstract https://doi.org/10.3847/1538-3881/abfec5 https://ui.adsabs.harvard.edu/abs/2021AJ....162...60B/abstract https://doi.org/10.1038/s41586-021-03869-6 https://ui.adsabs.harvard.edu/abs/2021Natur.598..272B/abstract https://doi.org/10.1046/j.1365-8711.2001.04776.x https://ui.adsabs.harvard.edu/abs/2001MNRAS.327..868B/abstract https://doi.org/10.1093/mnras/stx1049 https://ui.adsabs.harvard.edu/abs/2017MNRAS.469.2434B/abstract https://doi.org/10.1088/0004-6256/147/3/47 https://ui.adsabs.harvard.edu/abs/2014AJ....147...47B/abstract https://doi.org/10.1086/184937 https://ui.adsabs.harvard.edu/abs/1987ApJ...318L..57B/abstract https://doi.org/10.1088/0004-637X/756/1/48 https://ui.adsabs.harvard.edu/abs/2012ApJ...756...48C/abstract https://ui.adsabs.harvard.edu/abs/2000A&A...364..217D/abstract https://doi.org/10.1088/0004-637X/698/2/1826 https://ui.adsabs.harvard.edu/abs/2009ApJ...698.1826D/abstract https://doi.org/10.1088/0004-637X/695/2/970 https://ui.adsabs.harvard.edu/abs/2009ApJ...695..970D/abstract https://doi.org/10.1088/0004-637X/809/1/74 https://ui.adsabs.harvard.edu/abs/2015ApJ...809...74F/abstract https://doi.org/10.1086/306256 https://ui.adsabs.harvard.edu/abs/1998ApJ...506..533G/abstract https://doi.org/10.1146/annurev-astro-081811-125518 https://ui.adsabs.harvard.edu/abs/2012ARA&A..50..411G/abstract https://doi.org/10.1086/592738 https://ui.adsabs.harvard.edu/abs/2008ApJ...689.1044G/abstract https://doi.org/10.1051/0004-6361/201117959 https://ui.adsabs.harvard.edu/abs/2011A&A...534L..14G/abstract https://doi.org/10.1086/505421 https://ui.adsabs.harvard.edu/abs/2006ApJ...644L..37G/abstract https://doi.org/10.1086/592723 https://ui.adsabs.harvard.edu/abs/2008ApJ...689...53H/abstract https://doi.org/10.1086/306793 https://ui.adsabs.harvard.edu/abs/1999ApJ...512..864H/abstract https://doi.org/10.1086/116685 https://ui.adsabs.harvard.edu/abs/1993AJ....106..773H/abstract https://doi.org/10.3847/0004-637X/819/2/168 https://ui.adsabs.harvard.edu/abs/2016ApJ...819..168H/abstract https://ui.adsabs.harvard.edu/abs/2016ApJ...819..168H/abstract https://doi.org/10.1109/MCSE.2007.55 https://ui.adsabs.harvard.edu/abs/2007CSE.....9...90H/abstract https://doi.org/10.1088/0004-637X/711/2/731 https://ui.adsabs.harvard.edu/abs/2010ApJ...711..731J/abstract https://doi.org/10.1086/192235 https://ui.adsabs.harvard.edu/abs/1995ApJS..101..117K/abstract https://doi.org/10.3847/1538-4357/ac07a8 https://ui.adsabs.harvard.edu/abs/2021ApJ...917...78K/abstract https://doi.org/10.3847/2041-8213/ac17ec https://ui.adsabs.harvard.edu/abs/2021ApJ...918L...8K/abstract https://ui.adsabs.harvard.edu/abs/2021ApJ...918L...8K/abstract https://doi.org/10.5281/zenodo.4898012 https://doi.org/10.5281/zenodo.4898012 https://doi.org/10.3847/2041-8213/ac7442 https://ui.adsabs.harvard.edu/abs/2022ApJ...933L..23L/abstract https://doi.org/10.5281/zenodo.6522913 https://doi.org/10.1088/0004-637X/690/2/1463 https://ui.adsabs.harvard.edu/abs/2009ApJ...690.1463L/abstract https://doi.org/10.1016/j.newast.2009.12.002 https://ui.adsabs.harvard.edu/abs/2010NewA...15..433M/abstract https://doi.org/10.3847/1538-3881/aa74b2 https://ui.adsabs.harvard.edu/abs/2017AJ....154...35N/abstract https://doi.org/10.1088/0004-637X/769/2/88 https://ui.adsabs.harvard.edu/abs/2013ApJ...769...88N/abstract https://doi.org/10.1088/0004-637X/696/2/1407 https://ui.adsabs.harvard.edu/abs/2009ApJ...696.1407N/abstract https://ui.adsabs.harvard.edu/abs/1998A&A...339..483R/abstract https://doi.org/10.3847/1538-4357/ac739e https://ui.adsabs.harvard.edu/abs/2022ApJ...933...83S/abstract https://doi.org/10.1088/1538-3873/128/967/095004 https://ui.adsabs.harvard.edu/abs/2016PASP..128i5004S/abstract https://doi.org/10.1086/131977 https://ui.adsabs.harvard.edu/abs/1987PASP...99..191S/abstract https://doi.org/10.1086/375212 https://ui.adsabs.harvard.edu/abs/2003ApJ...591..204S/abstract https://doi.org/10.3847/2041-8213/aaf577 https://ui.adsabs.harvard.edu/abs/2018ApJ...869L..34S/abstract https://doi.org/10.3847/1538-4357/833/2/145 https://ui.adsabs.harvard.edu/abs/2016ApJ...833..145S/abstract https://doi.org/10.1088/0004-637X/780/2/123 https://ui.adsabs.harvard.edu/abs/2014ApJ...780..123S/abstract https://ui.adsabs.harvard.edu/abs/2011AcA....61...83S/abstract https://doi.org/10.3847/1538-3881/abcc60 https://ui.adsabs.harvard.edu/abs/2021AJ....161...54T/abstract https://ui.adsabs.harvard.edu/abs/1993AcA....43..289U/abstract https://ui.adsabs.harvard.edu/abs/2015AcA....65....1U/abstract https://doi.org/10.1088/0004-637X/725/1/331 https://ui.adsabs.harvard.edu/abs/2010ApJ...725..331Y/abstract https://doi.org/10.1038/s41550-022-01671-6 https://ui.adsabs.harvard.edu/abs/2022NatAs...6..782Z/abstract 1. Introduction 2. Updated Light-curve Modeling 2.1. Comparison to Previous Studies 2.1.1. C12 2.1.2. S16 3. Multi-epoch High-resolution Imaging with Keck 3.1. PSF Fitting 3.1.1. Astrometric and Photometric Errors with MCMC+Jackknife 3.2.2019 NIRC2 Analysis 3.3.2021 OSIRIS Analysis 4. Lens–Source Relative Proper Motion 5. Source and Lens Star Identification 6. A Search for a Wide Lens Companion 6.1. Residual Pixel Grid 6.2. Three-star MCMC Search 6.3. A Nonluminous Primary Lens 6.4. A Nonluminous Companion 7. Lens System Properties 8. Discussion and Conclusions References