A Close Binary Lens Revealed by the Microlensing Event Gaia20bof

During the last 25 yr, hundreds of binary stars and planets have been discovered toward the Galactic bulge by microlensing surveys. Thanks to a new generation of large-sky surveys, it is now possible to regularly detect microlensing events across the entire sky. The OMEGA Key Projet at the Las Cumbres Observatory carries out automated follow-up observations of microlensing events alerted by these surveys with the aim of identifying and characterizing exoplanets as well as stellar remnants. In this study, we present the analysis of the binary lens event Gaia20bof. By automatically requesting additional observations, the OMEGA Key Project obtained dense time coverage of an anomaly near the peak of the event, allowing characterization of the lensing system. The observed anomaly in the lightcurve is due to a binary lens. However, several models can explain the observations. Spectroscopic observations indicate that the source is located at ≤2.0 kpc, in agreement with the parallax measurements from Gaia. While the models are currently degenerate, future observations, especially the Gaia astrometric time series as well as high-resolution imaging, will provide extra constraints to distinguish between them.


INTRODUCTION
The gravitational microlensing effect (Einstein 1936) has been used for more than twenty years to detect faint objects in the Milky Way.Originally used to probe the nature of the dark matter in the galactic halo (Paczyński 1986) by observing toward the Magellanic Clouds, microlensing surveys are now focused on the Galactic Bulge, where the event rate is the highest.Thanks to the Optical Gravitational Lensing Experiment (Udalski et al. 2015), the Microlensing Observation in Astrophysics (Sumi et al. 2003), the Korea Microlensing Telescope Network (Kim et al. 2016) and follow-up teams1 , more than hundred exoplanets have been detected2 .More recently, the first isolated stellar-mass black hole has been detected by microlensing using precise astrometry from the Hubble Space Telescope (Sahu et al. 2022;Lam et al. 2022;Mróz et al. 2022).
A new generation of large-sky surveys opens the possibility to detect microlensing across the entire sky.The combination of large field-of-view, high spatial resolution, low limiting magnitudes and rapid data processing offer the possibility to discover microlensing event in the Galactic Disk.This opens the opportunity to study the galactic population of lenses that are difficult to observe otherwise.In particular, the Gaia mission of the European Space agency has detected hundreds of microlensing events (Wyrzykowski et al. 2023).Events observed by Gaia are often of special interest for several reasons.The duration of events in the Galactic Disk is usually longer, allowing the detection of additional effects such as the microlensing parallax (Gould 1994) or the orbital motion of the lens (Wyrzykowski et al. 2020), which provide unique constraints on the mass and distances of the lenses.Moreover, the future astrometric time series delivered by the Gaia mission offers the possibility to measure the angular Einstein ring radius θ E via the astrometric microlensing signal (Rybicki et al. 2018), providing a mass/distance relation for a large fraction of lenses.By design, these surveys generally deliver a weekly cadence of observation, which is not dense enough to accurately catch anomalous features in microlensing events.Therefore, the follow-up of these events is of paramount importance to ensure their characterisation (Tsapras 2018).
As with most microlensing events detected by the Gaia mission, Gaia20bof (equatorial: (α = 184.61816• , δ = −63.49726• , J2000), galactic (l = 299.26406• ,b = −0.86052)) is located in the Galactic Disk.The event was announced by the Gaia Science Alerts (GSA3 hereafter) (Hodgkin et al. 2021) on 30 March 2020 and the microlensing nature of the event was confirmed via spectroscopic classification from Southern African Large Telescope, SALT (Ihanec et al. 2020), on the 10 June 2020, because the data displays a stellar spectrum with only absoprtion lines.In response to this alert, several groups started to observe this event and identified an anomalous feature at the event peak.We present the different observations collected for this event in Section 2. The modelling of the photometric observations presented in Section 3.2 reveals that the lightcurves can be equally well explained by different competing models.Section 3.3 presents the analysis of the different spectra collected and confirms the measurement from Gaia revealing that the source is a red subgiant located at 2 kpc.As discussed in Section 4 and Section 5, this implies that the lens is relatively close and that additional follow-up observations will place additional constraints on the models in the near future.

Gaia
Gaia's photometric measurements consist of a wide G-band (Jordi et al. 2010) obtained with roughly monthly cadence while Gaia scans the sky.The photometry is publicly available on the GSA web page (Hodgkin et al. 2021).We used the procedure described in Kruszyńska et al. (2022) to obtain the photometric errors of the 163 measurements used in this study.

The OMEGA Key Project
The OMEGA Key Project is an international collaboration that performs automatic follow-up of microlensing events detected by large-sky surveys.The primary goal is the characterization of cold planets and stellar remnants in the entire sky.Microlensing candidates are first collected from various channels, including from GSA, in the Microlensing Observing Platform (MOP)4 .The MOP system is a Target and Observation Manager, or TOM system, built with the TOM Toolkit package (Street et al. 2018).As soon as data are available, this system automatically fits a single lens model (including microlensing parallax, see for example Gould (2004)) to all of the data available for each event, using the pyLIMA modelling software (Bachelet et al. 2017).Typically, there are more potential targets at any given time than there are telescope resources to observe them and it is therefore necessary to prioritize targets.MOP incorporates an algorithm that prioritize targets based on their position in the sky and their current status (Hundertmark et al. 2018).For instance, events located towards the galactic Bulge (255 • ≤ α ≤ 275 • and −36 • ≤ δ ≤ −22 • ) are observed only if they are sufficiently sensitive to planets according to Hundertmark et al. (2018).This is motivated by the fact that this region is regularly monitored by OGLE, MOA and KMTNet surveys.Otherwise, higher cadence imaging and spectroscopic observations are requested automatically by the MOP system, tailored to current parameters and phase of the evolution of each event(∼ daily).Most of the observations are collected via the Las Cumbres Observatory (LCO) automatic robotic telescopes Total 746 a.These observations were taken solely at the baseline and hence were not used for the modelling.
network (Brown et al. 2013).The observing strategy consists of regular monitoring of the event in two bands, namely SDSSg' and SDSS-i' (with the exception of events located towards the Galactic Bulge, where the SDSS-g' band is replaced by the SDSS-r' band, due to the higher extinction in these fields), with a cadence depending of the event priority and Einstein ring crossing time t E .Observing in two bands serves two primary purposes.First, it measures the chromaticity of ongoing events.Indeed, microlensing events are achromatic phenomena (to first order) and therefore the color evolution contributes to exclude astrophysical false-positive detections (such as Be stars) but also helps to distinguish between microlensing models (namely, the double source/double lens scenario, see below).Second, as emphasized by Yoo et al. (2004), observations in (at least) two bands allows the estimation of the angular source radius, a key component in estimating the mass of the lenses.We note that the cadence in the SDSS-i' band is twice that of the other bands in order to increase the sampling of the lightcurve.If at any given time, models predict a high planet sensitivity (during a high magnification event for example) (Hundertmark et al. 2018), a 15-minute cadence mode is triggered for the next 48 hours to ensure that any potential anomalies are well sampled.If a target is predicted to exceed a brightness threshold of 17 mag, several low-resolution spectra (R∼1000) using the FLOYDS instruments (Brown et al. 2013) are also requested, to help with the source characterization, see for example Fukui et al. (2019) as well as to classify contaminants in the alerts stream, such Young-Stellar Objects.If the event gets bright enough (i.e.V ≤ 11 mag), high-resolution NRES spectra (R∼55000) (Siverd et al. 2018) can also be triggered to obtain a precise estimation of the source spectral parameters, such as in the event Gaia19bld (Bachelet et al. 2022).
As listed in Table 1, hundreds of images have been collected in the SDSS-g' and SDSS-i' bands, from La Silla in Chile (LSC), Siding Spring in Australia (COJ) and the South African Astronomical Observatory in South Africa (CPT) LCO sites.Some of these images come from precursor observing programs prior to OMEGA.We note that, with the exception of the LSC site that remained closed for a long time, the COVID-19 pandemic only mildly impacted this observing program at LCO.

Other follow-up data
The target has been followed up photometrically with a network of small telescopes under the umbrella of the OPTICON-RadioNET Pilot program of the European Commission, e.g.Wyrzykowski et al. (2020).For Gaia20bof we used LCO, TRT and Skynet's PROMPT-MO and PROMPT5 telescopes, the summary of these observations is gathered in Table 1.
TRT stands for Thai Robotic Telescopes, a 0.7m telescope equipped with Andor iKon-L DW936 BV CCD camera with resolution 0.8 arcseconds/pixel.

Spectroscopic data
We collected two spectra during the course of the event: one is from 10-m Southern African Large Telescope (SALT; Buckley et al. 2006) equiped with the Robert Stobie Spectrograph (Burgh et al. 2003;Kobulnicky et al. 2003), and the second is from the X-Shooter intrument (Vernet et al. 2011) mounted on the ESO 8-m Very Large Telescope (VLT).
The SALT/RSS low-resolution spectrum was obtained on 2020 June 6 (JD∼2459006) at event magnification ∼1.85 , with the exposure time 600 s.The longslit mode was used with the slit width 1.5 arcsec, grating pg0300.The wavelength range of the obtained spectrum covers 370-930 nm giving the resolving power R ∼ 350 and average signal-to-noise ratio S/N=160.It has been reduced in a standard way using PySALT6 software (bias subtraction and flat-field correction; Crawford et al. 2010) and then wavelength and flux calibration was applied using standard IRAF routines thanks to having calibrating data for Ar comparison lamp and spectrophotometric standard star.
The VLT/X-Shooter spectrum was obtained at the baseline of the microlensing event on January 7, 20227 for three wavelength channels: UVB (300 − 559.5 nm), VIS (559.5 − 1024 nm) and NIR (1024 − 2480 nm).We integrated with 191 s, 220 s and 300 s of exposure time for UVB, VIS and NIR channel, respectively.The resolution of this spectrum is R ∼ 10000 in UVB part (at slit width 1.0 arcsec) and R ∼ 15000 in VIS (at slit width 0.7 arcsec) and NIR (at slit width 0.6 arcsec) parts, on average.It has been reduced with the dedicated EsoReflex8 pipeline (v.2.9.1).For the calibration of UVB, VIS wavelengths, ThAr lamp was used, while for NIR -a set of Ar, Hg, Ne and Xe lamps.Due to the poor quality of some spectral ranges and low signal-to-noise ratio, only the parts between 350 − 552 nm, 560 − 745 nm, 780 − 914 nm and 1133 − 1350 nm, 1450 − 1800 nm, 1950 − 2356 nm were used in further analysis.

Photometric reductions
The follow-up observations were reduced using the Black Hole TOM (BHTOM) infrastructure9 , which utilised CCDPhot suite of image processing and CPCS photometry calibration tools (Zieliński et al. 2019(Zieliński et al. , 2020)).The photometry was calibrated to the Gaia Synthetic Photometry catalog Gaia synthetic SDSS magnitudes (Montegriffo et al. 2022).We note that the magnitudes are in the AB system (Oke & Gunn 1983).

Modelling of microlensing event
As can be seen in Figure 1, the OMEGA follow-up data clearly reveal an anomaly around JD ∼ 2458970.First, we explored the binary lens versus the binary source interpretation (Dominik et al. 2019) and found that the latter can be safely ruled out because of the high χ 2 value (about two times higher than the binary lens solutions).However, the anomaly is mild, with no clear sign of a caustic crossing and looks similar to the approach of a Chang-Refsdal lens (Chang & Refsdal 1979).Such lightcurves are notoriously known for presenting strong degeneracies in terms of lens geometries, see for instance Dominik (1999); Han & Gaudi (2008); Shvartzvald et al. (2016); Bozza et al. (2016), implying that several models can reproduce the observations accurately, which is exactly the case here.A binary model is parametrised by the time t 0 of the minimum impact parameters u 0 relative to the centre of mass of the binary lens.The angular Einstein ring radius crossing time is defined as t E = θ E µ rel , where µ rel is the (geocentric) relative proper motion.We also considered ρ = θ ⋆ /θ E , the normalised source radius where θ * is the angular source radius, in the models, but this parameter is not constrained in the present case, because the source trajectories are far from the caustic.For this reason, we did not consider limb-darkening effects.The models also include the normalised angular projected separation between the two components of the lens s, and their mass ratio q.Finally, the angle between the lens trajectory and the binary axes is defined by α (counter-clockwise).Because of the event duration (t E ≥ 50 d), we also consider the microlensing parallax vector π E = (π EN , π EE ) (Gould 2004) and set the time of reference t 0,par = t 0 for all models.The norm of the parallax vector is directly related to the lens and source distances via: with π LS = 1/D l −1/D s , while the parallax vector is co-linear to the lens-source relative proper motion projected in the North-East plan of the sky.The first steps of modelling have been done using the RTModel infrastructure10 , which located eight minima.We refine these solutions with a Monte-Carlo Markov Chain (MCMC) exploration on each of them as reported in the Table 3.2.The posterior exploration has been performed with the updated version 2.0 of pyLIMA (Bachelet et al. 2017), which maximizes the total log-likelihood L: where N telescopes is the total number of telescopes, N data is the total number of photometric observations of a given telescope, σ t are the uncertainties associated with the measured flux f t at the time t and m t is the corresponding microlensing model flux.
The version 2.0 of pyLIMA introduces the option of rescaling the flux uncertainties σ i,t of each telescope i during the MCMC exploration using A detailed presentation of the new features available with the latest release of pyLIMA will be the subject of a separate paper (Bachelet et al., in prep.).The MCMC exploration is performed with the emcee package (Foreman-Mackey et al. 2013) with a fixed number of 10000 chains and 32 walkers.The convergence of the chains have been analyzed separately for each model.A summary of the best fitting models found is displayed in Table 3.2, and the best model is presented in Figure 1.For each parameter set, the 16, 50 and 84 percentiles of the chains are shown.We note that the reported χ 2 values for each model are obtained from a gradient fit started at the best MCMC chain-position (and includes the rescaled uncertainties obtained from the chains).The Wide+ model (i.e.s ≥ 1) is currently slightly favoured and points towards a stellar binary lens interpretation, but the planetary companion scenario remains possible at this time.

Spectroscopic properties of the source
Based on the low-resolution SALT/RSS spectrum we were able to classified the source as a reddened GK-type star (Ihanec et al. 2020).Data shows the stellar spectrum with only absorption lines, therefore, the changes of the brightness of Gaia20bof could be explained by the gravitational microlensing phenomenon.Therefore, we decided to continue with extensive photometric follow-up monitoring of Gaia20bof event.
The high-resolution spectrum from X-Shooter has been used for determining of the atmospheric parameters, i.e., effective temperature T eff , surface gravity log g and metallicity [M/H], of the source star.It was possible thanks to a spectral line fitting method by using iSpec11 package (Blanco-Cuaresma et al. 2014;Blanco-Cuaresma 2019).We used the SPECTRUM12 radiative transfer code to generate a set of synthetic spectra based on a grid of MARCS models (Gustafsson et al. 2008), solar abundances taken from Grevesse et al. (2007) and line list with atomic data from Gaia-ESO Survey (GESv6; Heiter et al. 2021).The GESv6 line list covers the wavelength range from 420 to 920 nm so this method uses only UVB and VIS part of the X-Shooter spectrum.The best-matching fit was found for the following parameters: T eff = (5533 ± 89) K, log g = (3.54± 0.19), [M/H] = (−0.51± 0.07) dex, and is presented together with the observational X-Shooter data in Fig. 2.
Moreover, both spectra have been modeled with a template matching method using Spyctres, similarly to Bachelet et al. (2022); Bachelet (2024).The latest version of Spyctres includes an update of the extinction law from Wang & Chen (2019).Briefly, we modeled the two spectra using the stellar template library from Kurucz (1993) as well as the SED at the time of spectra acquisition, including the source magnification A(t).This allows an accurate flux calibration and ultimately the estimation of A V and the stellar parameters.With this method, the final solution for the source star parameters was found as T eff = (5297 ± 30) K, log g = (3.50+0.30  −0.25 ), [M/H] = (−0.7 +0.3 −0.1 ) dex, together with the line-of-sight extinction in V band: A V = 1.55 +0.03 −0.04 mag.Both spectra and the results of template matching are visible in Figure 3 while the modelling results are displayed in Table 3.The results obtained in spectroscopic analysis are in good agreement with the measurements from the GaiaDR3 release, where T eff = (5434±16) K, log g = (3.53±0.02),[M/H] = (−0.520±0.001)dex (GaiaDR3 6054150372473485696; Gaia Collaboration et al. 2016Collaboration et al. , 2022)).

Distance to the source
Using the derived parameters and the PARSEC stellar isochrones13 (Bressan et al. 2012), we found that the source is most likely an old subgiant of G2 spectral type, as can be seen in Figure 4. Combining the source luminosity, the source apparent magnitude in the G band and the extinction law derived from spectral analysis, we found that the source angular radius θ * = 5.4 ± 0.2 µas is independent of the source age and distance.We note that the angular radius value is also confirmed by the color-radius relation from Boyajian et al. (2014) with θ * = 5.2 ± 0.2 µas (see Table 3.2).
Based on the stellar isochrones, the extinction obtained from the spectra fits and assuming the source magnitude G = (16.00± 0.05) mag, we estimated the source distance to be D S = 2.1 +0.5 −0.2 kpc, as can be seen in the right panel of Fig. 4. We note that using different source magnitude bands (like SDSS-i') and models leads us to similar conclusion.
Moreover, assuming the typical absolute magnitude M V = (3.0 ± 0.3) mag (Straižys 1992) for metal-poor G2-type subgiant, we also calculated the distance to the source analytically.This way, we obtained the value of D S = (1.95 ± 0.35) kpc.
These estimates are in excellent agreement with the parallax measurements from Gaia: D S = 2.33 +0.03 −0.02 kpc (distance_gspphot).In addition, the value of the distance inferred for the Gaia20bof parallax published in GaiaDR2 (Bailer-Jones et al. 2018) is D S = 4.32 +2.94  −1.81 kpc, while the updated value for the parallax published in GaiaEDR3 assuming geometric and photogeometric  observed blend.Therefore, the lens system is a low mass and relatively close binary, but its exact nature will be revealed with additional observations in the near future.

Gaia astrometry
The next data release from the Gaia Mission (DR4, ∼ 2025), will include astrometric time series that will help constrain the models.While the photometric lightcurves are almost identical for all of the models, the astrometric microlensing signals (Walker 1995;Dominik & Sahu 2000) can differ significantly.In Figure 5 we show the astrometric microlensing signals as seen by Gaia for the eight models presented, assuming θ E = 1 mas.It is clear that the exquisite astrometric precision of Gaia will allow the selection of the most plausible model (or, at least, eliminate most of them) as well as measuring θ E .This has already been done for the event Gaia16aye (Wyrzykowski et al. 2020) 14 , which is similar in brightness to Gaia20bof, but presents more caustic crossing features.In the case of Gaia16aye, the early access to the astrometric time series confirmed both the microlensing model and the angular Einstein ring radius measurement of ∼ 3 mas.

High resolution imaging
By directly measuring the lens fluxes and/or the relative proper motions, high resolution imagers have been able to place strong constraints on microlensing systems, see for example Beaulieu (2018); Bhattacharya et al. (2018); Vandorou et al. (2020).In the case of Gaia20bof, high resolution imaging will also bring new constraints.First, the models predict different blend magnitudes in different bands, as listed in the Table 3.2.A direct measurement of the lens flux, as well as a better analysis of the blend, will limit the number of possible scenarios.Moreover, the lens/source relative proper motions are significantly different between the models.Indeed, given that ⃗ π E = π E ⃗ µ LS the (geocentric) relative proper motion can be written as (Gould 2004): All models predict different timescales t E and parallax vectors π E (especially the Wide models).Therefore, observations of the lens and the source at different epochs will discriminate between models and ultimately provide a measurement of θ E (Vandorou et al. 2020).Assuming a typical value of θ E ∼ 1 mas, high-resolution observations could start in about 5 years, with a lens and source separation δ ≥ 10 mas for all models except Wide-.

Radial velocity measurements
As soon as the lens and source separate, it will be possible to conduct radial-velocity follow-up of the host that should lead to the full characterization of the system in combination with the two methods described previously.The radial-velocity monitoring of microlensing lenses is challenging, due to the faintness of the host.However, it has been done in the past on at least two occasions.Yee et al. (2016) used the Keck High-resolution Echelle Spectrometer (Vogt et al. 1994) and the Magellan Inamori Kyocera Echelle spectrometer (Bernstein et al. 2003) to measure the radial velocity signal of the host of the microlensing event OGLE-2009-BLG-020 (Skowron et al. 2011).Their data confirms and refine the stellar binary parameters from the original publication.A second radial velocity test has been made by Boisse et al. (2015) on the microlensing predictions of OGLE-2011-BLG-0417 (Shin et al. 2012).In this case, Boisse et al. (2015) did not measure any modulations in the radial velocity of the host, indicating strong tensions with the microlensing models later confirmed by high-resolution imaging (Santerne et al. 2016).Ultimately, Bachelet et al. (2018) identified a competitive microlensing model that explains the lack of modulations.In the case of Gaia20bof, radial velocity observations will be challenging, as the lens cannot be brighter than the measured blend (ib ≤ 18 mag).However, it could be done with the most recent spectrographs, such as ESO ESPRESSO (Pepe et al. 2021).

CONCLUSIONS
We report the detection of a close (≤ 2 kpc) binary lens by combining space-based and ground-based time series photometry.Spectroscopic data indicates that the source is a relatively close (D S ≤ 2 kpc) red sub-giant, in good agreement with parallax measurements from Gaia.The photometric lightcurves can be explained by a binary lens model.However, several degenerate models can reproduce the observations.We discuss several methods to distinguish between the models in the future.First, high-resolution imaging can confirm or reject the model prediction on the relative proper-motion direction, as well as give a measurement of θ E .Indeed, with a lens distance D L ≤ 2 kpc, it is almost certain that the lens will be observable with current facilities such as Keck.High-resolution images would also provide a measurement of the host flux that will need to be compared with the blend values reported for each model.Secondly, the Gaia astrometric time series, expected with the Gaia DR 4, would be extremely useful for constraining the lens properties.Indeed, the eight models predict significantly different astrometric shifts that should be measurable if θ E is not too small (i.e.θ E ≥ 0.5 mas).In this case, the astrometric time series will allow the estimation of the lens mass via the direct measurement of θ E .Finally, we discussed the possibility of measuring the radial velocity of the host.This would be challenging, as we expect the host to be faint, but could be done with the most precise instruments.

Figure 1 .
Figure 1.Lightcurves of the event observed from the ground and Gaia.The Wide+ model is presented (blue solid line, the dashed line represents the model seen by Gaia) but all models produce almost identical lightcurves.The parallax signal between Gaia and Earth is almost not visible.All datasets have been artificially shifted to Gaia magnitudes for presentation purposes.The presented uncertainties take into account the rescaling parameters derived during the posterior exploration, see the main text for more details.

Figure 3 .
Figure3.The two spectra from SALT (black) and X-Shooter (red), as well as the best models, are visible.The gray vertical lines indicates telluric bands, where the data were not used for the modelling.

Figure 4 .
Figure 4. (Left) PARSEC stellar isochrones for 9, 9.5 and 10 Gyr with a fixed metallicity of -0.5.The source is most likely an old subgiant.(Middle) Source angular radius based on the spectra models as a function of the source distance and age (color-coded).The points are derived from the template matching MCMC chains and the stellar isochrones.(Right) Source distance based on the Gaia source measurements from the model (G s ∼ 16 mag) and the extinction estimated from the template matching modeling.Black dash lines indicates the 1-σ confidence region from the parallax measurement of Gaia.

Figure 5 .
Figure 5. (Left) The source trajectories (blue), caustics (red) and critical curves (black) for the eight models.Note that the centre of mass of the lenses is kept fixed at (0,0).The microlensing astrometric deflections in right ascension (Middle) and in declination (Right) versus time are also displayed.The color indicates the observed Gaia G magnitude.
EB gratefully acknowledge support from NASA grant 80NSSC19K0291.EB's work was carried out within the framework of the ANR project COLD-WORLDS supported by the French National Agency for Research with the reference ANR-18-CE31-0002.This work was authored by employees of Caltech/IPAC under Contract No. 80GSFC21R0032 with the National Aeronautics and Space Administration.This research has made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program.This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia),processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium).Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France.YT acknowledges the support of DFG priority program SPP 1992 "Exploring the Diversity of Extrasolar Planets" (TS 356/3-1).This work is supported by Polish NCN grants: Daina No. 2017/27/L/ST9/03221, Harmonia No. 2018/30/M/ST9/00311 and MNiSW grant DIR/WK/2018/12.The BHTOM project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No. 101004719 (OPTICON RadioNet Pilot, ORP).BHTOM acknowledges the following people who helped with its development: Patrik Sivak, Kacper Raciborski, Piotr Trzcionkowski and AKOND company.This paper made use of the Whole Sky Database (wsdb) created by Sergey Koposov and maintained at the Institute of Astronomy, Cambridge by Sergey Koposov, Vasily Belokurov and Wyn Evans with financial support from the Science & Technology Facilities Council (STFC) and the European Research Council (ERC), with the use of the Q3C software ( http://adsabs.harvard.edu/abs/2006ASPC..351..735K ).RFJ acknowledges support for this project provided by ANID's Millennium Science Initiative through grant ICN12_009, awarded to the Millennium Institute of Astrophysics (MAS), and by ANID's Basal project FB210003.DAHB acknowledges support from the South African National Research Foundation.PAW acknowledges financial support from the NRF (grant no.129359) and UCT.The SALT observations were obtained under the SALT Large Science Programme on transients (2018-2-LSP-001; PI: DAHB).Polish participation in SALT is funded by grant number MEiN nr 2021/WK/01.

Table 1 .
Summary of the observations.
The telescope is located at 259.3045 deg West and 26.6955 deg North.
PROMPT-MO is a 0.4 m RC telescope equipped with Apogee USB CCD with resolution 0.598 arcseconds/pixel.The telescope is located at 243.011 deg West and 31.638deg South in the Meckering Observatory, Australia.PROMPT5 is another Skynet Network telescope and is a Ritchey-Chretien 0.41 m telescope operating using Apogee CCD camera at 0.8 arcseconds/pixel resolution.The telescope is located in CTIO Chile at 70.8053889 deg West and 30.1676389 deg South.

Table 2 .
Best models from the modelling of the lightcurves.The microlensing parameters and errors corresponds to the 16, 50 and 84 percentiles of the MCMC chains.The χ 2 corresponds to the maximum-likelihood and are based on rescaled uncertainties rescaling.

Table 3 .
Summary of the source properties from the spectral analysis and GaiaDR3.