Astronomy &Astrophysics A&A 674, A89 (2023) https://doi.org/10.1051/0004-6361/202346166 © The Authors 2023 MOA-2022-BLG-249Lb: Nearby microlensing super-Earth planet detected from high-cadence surveys Cheongho Han1, Andrew Gould2,3, Youn Kil Jung4,5, Ian A. Bond6, Weicheng Zang7,8, Sun-Ju Chung4,7, Michael D. Albrow9, Kyu-Ha Hwang4, Yoon-Hyun Ryu4, In-Gu Shin7, Yossi Shvartzvald10, Hongjing Yang8, Jennifer C. Yee7, Sang-Mok Cha4,11, Doeon Kim1, Dong-Jin Kim4, Seung-Lee Kim4, Chung-Uk Lee4, Dong-Joo Lee4, Yongseok Lee4,11, Byeong-Gon Park4, Richard W. Pogge3, (The KMTNet collaboration) Shude Mao8, Wei Zhu8, (Microlensing Astronomy Probe Collaboration) Fumio Abe12, Richard Barry13, David P. Bennett13,14, Aparna Bhattacharya13,14, Hirosame Fujii12, Akihiko Fukui15,16, Ryusei Hamada17, Yuki Hirao17, Stela Ishitani Silva14,18, Yoshitaka Itow12, Rintaro Kirikawa17, Iona Kondo17, Naoki Koshimoto19, Yutaka Matsubara12, Sho Matsumoto17, Shota Miyazaki17, Yasushi Muraki12, Arisa Okamura17, Greg Olmschenk13, Clément Ranc20, Nicholas J. Rattenbury21, Yuki Satoh17, Takahiro Sumi17, Daisuke Suzuki17, Taiga Toda17, Mio Tomoyoshi17, Paul J. Tristram22, Aikaterini Vandorou13,14, Hibiki Yama17, and Kansuke Yamashita17, (The MOA Collaboration) (Affiliations can be found after the references) Received 16 February 2023 / Accepted 5 April 2023 ABSTRACT Aims. We investigate the data collected by the high-cadence microlensing surveys during the 2022 season in search of planetary signals appearing in the light curves of microlensing events. From this search, we find that the lensing event MOA-2022-BLG-249 exhibits a brief positive anomaly that lasted for about one day, with a maximum deviation of ∼0.2 mag from a single-source, single-lens model. Methods. We analyzed the light curve under the two interpretations of the anomaly: one originated by a low-mass companion to the lens (planetary model) and the other originated by a faint companion to the source (binary-source model). Results. We find that the anomaly is better explained by the planetary model than the binary-source model. We identified two solutions rooted in the inner–outer degeneracy and for both of them, the estimated planet-to-host mass ratio, q ∼ 8 × 10−5, is very small. With the constraints provided by the microlens parallax and the lower limit on the Einstein radius, as well as the blend-flux constraint, we find that the lens is a planetary system, in which a super-Earth planet, with a mass of (4.83 ± 1.44) M⊕, orbits a low-mass host star, with a mass of (0.18 ± 0.05) M�, lying in the Galactic disk at a distance of (2.00 ± 0.42) kpc. The planet detection demonstrates the elevated microlensing sensitivity of the current high-cadence lensing surveys to low-mass planets. Key words. planets and satellites: general 1. Introduction The microlensing method of finding planets has various advan- tages that can serve as a complement to other planet detection methods. In particular, it provides a unique tool for detecting planets belonging to faint stars because the lensing character- istics do not depend on the object’s light. This method is also useful for detecting outer planets because of the high microlens- ing sensitivity to planets lying at around the Einstein radius, which approximately corresponds to the snow line of a planetary system. We refer to the review paper of Gaudi (2012) for a dis- cussion of the various advantages of the microlensing method. Another important advantage of the microlensing method is its high sensitivity to low-mass planets (Bennett & Rhie 1996). In general, the microlensing signal of a planet appears as a discontinuous perturbation to the smooth and symmetric lensing light curve produced by the host of the planet (Mao & Paczyński 1991; Gould & Loeb 1992). The amplitude of the planetary microlensing signal weakly depends on the planet-to- host mass ratio, q, although the duration of the signal becomes shorter in proportion to q−1/2 (Han 2006). This implies that the microlensing sensitivity can extend to lower-mass planets as the observational cadence becomes higher. The observational cadence of microlensing surveys has been greatly enhanced over the 2010s, thanks to the replacement of the cameras installed on the telescopes of previously established surveys of the Microlensing Observations in Astrophysics survey (MOA: Bond et al. 2011) and the Optical Gravitational Lensing Experiment (OGLE: Udalski et al. 2015). Its new cameras now have very wide fields of view (FOVs). In addition, a new sur- vey of the Korea Microlensing Telescope Network (KMTNet: Kim et al. 2016) was initiated. With these high-cadence surveys, A89, page 1 of 10 Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication. https://www.aanda.org https://doi.org/10.1051/0004-6361/202346166 https://www.edpsciences.org/en/ https://creativecommons.org/licenses/by/4.0 https://www.aanda.org/subscribe-to-open-faqs mailto:subscribers@edpsciences.org A&A 674, A89 (2023) Table 1. Low-mass microlensing planets. Planet Type Reference OGLE-2005-BLG-390Lb super-Earth Beaulieu et al. (2006) MOA-2007-BLG-192Lb super-Earth Bennett et al. (2008) MOA-2009-BLG-266Lb super-Earth Muraki et al. (2011) MOA-2011-BLG-262Lb super-Earth Bennett et al. (2014) MOA-2013-BLG-605Lb super-Earth Sumi et al. (2016) OGLE-2013-BLG-0341Lb terrestrial planet Gould et al. (2014) OGLE-2016-BLG-1195Lb Earth-mass planet Shvartzvald et al. (2017), Bond et al. (2017), Vandorou et al. (2023) OGLE-2016-BLG-1928L terrestrial-mass rogue planet Mróz et al. (2020) OGLE-2017-BLG-0482Lb super-Earth Han et al. (2018) OGLE-2017-BLG-1806Lb super-Earth Zang et al. (2023) KMT-2017-BLG-0428Lb super-Earth Zang et al. (2023) KMT-2017-BLG-1003Lb super-Earth Zang et al. (2023) KMT-2017-BLG-1194Lb super-Earth Zang et al. (2023) OGLE-2018-BLG-0532Lb super-Earth Ryu et al. (2020) OGLE-2018-BLG-0677Lb super-Earth Herrera-Martín et al. (2020) OGLE-2018-BLG-0977Lb super-Earth Hwang et al. (2022) OGLE-2018-BLG-1185Lb super-Earth Kondo et al. (2021) KMT-2018-BLG-1025Lb super-Earth Han et al. (2021) KMT-2018-BLG-1988Lb super-Earth Han et al. (2022a) KMT-2018-BLG-0029Lb super-Earth Gould et al. (2020) OGLE-2019-BLG-0960Lb super-Earth Yee et al. (2021) OGLE-2019-BLG-1053Lb terrestrial planet Zang et al. (2021b) KMT-2019-BLG-0253Lb super-Earth Hwang et al. (2022) KMT-2019-BLG-1367Lb super-Earth Zang et al. (2023) KMT-2019-BLG-1806Lb super-Earth Zang et al. (2023) KMT-2020-BLG-0414Lb Earth-mass planet Zang et al. (2021a) KMT-2021-BLG-0912Lb super-Earth Han et al. (2022b) KMT-2021-BLG-1391Lb super-Earth Ryu et al. (2022) lensing events can be observed with a cadence down to 0.25 h, compared to the one-day cadence provided by earlier surveys. The detection rate of very low-mass planets has greatly increased with the enhanced sensitivity to very short anoma- lies in microlensing light curves from this higher observational cadence. In Table 1, we list the discovered microlensing plan- ets with masses below that of a super-Earth planet, together with brief comments of the planet types and related references. Among these 28 planets, 5 are terrestrial planets with masses similar to that of Earth and the other 23 are super-Earth planets with masses higher than Earth’s but substantially below those of ice giants in the Solar System, (i.e., Uranus and Neptune). To be noted is that 22 planets (79%) have been detected since the full operation of the current high-cadence lensing surveys. Very low-mass planets detected before the era of high-cadence sur- vey were found using a specially designed observational strategy, in which survey groups focused on detecting lensing events and followup groups densely observed the events found by the survey groups with the employment of multiple narrow-FOV telescopes. However, the detection rate of low-mass planets based on this strategy was low because of the limited number of events that could be observed by followup groups. In contrast, high-cadence surveys can densely monitor all lensing events without the need of extra followup observations. In this paper, we report the discovery of a super-Earth planet found from inspection of the 2022 season microlensing data collected by the KMTNet and MOA surveys. The planet was dis- covered by analyzing the light curve of the microlensing event MOA-2022-BLG-249, for which a very short-term anomaly was covered by the survey data despite its weak deviation. We check various interpretations of the signal and confirm its planetary origin. We present our analysis according to the following organi- zation. In Sect. 2, we describe the observations of the planetary lensing event and the data obtained from these observations. In Sect. 3, we depict the characteristics of the event and the anomaly appearing in the lensing light curve. We present the analyses of the light curve conducted under various interpretations of the anomaly and investigate higher order effects that affect the lensing-magnification pattern. We identify the source star of the event and check the feasibility of measuring the angular Einstein radius in Sect. 4. We estimate the physical parameters in Sect. 5. We summarize the results of the analysis and present our conclusions in Sect. 6. 2. Observations and data The microlensing event MOA-2022-BLG-249 occurred on a source lying toward the Galactic bulge field at (RA,Dec)J2000 = (17:55:27.73–28:18:21.82), (l, b) = (+1.◦65,−1.◦53). The mag- nification of the source flux induced by lensing was first found by the MOA group on 2022 May 22, which corresponds to the abridged heliocentric Julian date of HJD′ ≡ HJD − 2450000 = 9721.48. The KMTNet group identified the event at HJD′ = 9721.63 (4 h after the MOA discovery) and designated the event as KMT-2022-BLG-0874. Hereafter, we refer to the event as MOA-2022-BLG-249 in accordance with the convention of the A89, page 2 of 10 Han, C., et al.: A&A proofs, manuscript no. aa46166-23 Fig. 1. Light curve of the microlensing event MOA-2022-BLG-249. The arrow marked by tanom in the second panel indicates the location of the anomaly. Top panel shows the enlarged view around the anomaly region. The solid and dashed curves drawn over the data points are 1L1S models obtained with (parallax model) and without (standard model) the consideration of microlens-parallax effects. Two lower panels shows the residuals from the two models. microlensing community using the event ID reference of the first discovery group. The event lies approximately 100′′ outside the footprint of the OGLE survey. In any case, there are no data from the survey because the OGLE telescope was shut down during most of the 2022 season due to the Covid-19 pandemic. The source location corresponds to a sub-prime field of the MOA survey and, thus, the coverage of the event is relatively sparse. In contrast, the source was in the KMTNet prime fields of BLG02 and BLG42, where observations were conducted with a high combined cadence of 0.25 h. Thus, the light curve of the event was densely covered by the KMTNet data. The event was addi- tionally observed by a survey of the Microlensing Astronomy Probe (MAP) collaboration, with a cadence of one or two points per night. The source flux gradually increased until the lensing light curve reached its peak on 2022 May 27 (HJD′ ∼ 9727) and then returned to the baseline. The duration of the event is very long, and the lensing magnification lasted throughout the whole 2022 bulge season. The event was observed with the use of multiple telescopes operated by the individual survey and followup groups. The MOA group utilized the 1.8 m telescope of the Mt. John Obser- vatory in New Zealand, the KMTNet group made the use of the three identical 1.6 m telescopes lying at the Siding Spring Obser- vatory in Australia (KMTA), the Cerro Tololo Interamerican Observatory in Chile (KMTC), and the South African Astro- nomical Observatory in South Africa (KMTS). The MAP group used the 3.6 m Canada-France-Hawaii Telescope (CFHT) in Hawaii. Data reduction and photometry of the event were carried out using the photometry pipelines of the individual groups, and the error bars of the individual data sets were readjusted using the routine described in Yee et al. (2012). Figure 1 shows the lensing light curve of MOA-2022-BLG- 249. The solid and dashed curves drawn over the data points are single-source single-lens (1L1S) models obtained from mod- eling with (parallax model) and without (standard model) the consideration of microlens-parallax effects (Gould 1992). A detailed discussion on the parallax effects is presented in Sect. 3. Although the observed light curve appears to be well described by the 1L1S model, we find that there exists a brief anomaly appearing at tanom ∼ 9733, which corresponds to about six days after the peak of the light curve. The upper panel of Fig. 1 shows the enlarged view of the region around the anomaly. The anomaly exhibited a positive deviation from the 1L1S model, and it lasted for about one day with a maximum deviation of ∆I ∼ 0.2 mag. The anomaly was mostly covered by the combi- nation of the KMTS and KMTA data sets and the region just before the major deviation was additionally covered by the two CFHT data points. The anomaly during the time gaps among the KMTS and KMTA coverage could have been covered by the KMTC data set, but the Chilean site was clouded out during the five-day period around the time of the anomaly. Similarly, the MOA group did not observe this field during a 12 day interval that included the anomaly. 3. Light curve analysis It is known that a brief positive anomaly in a lensing light curve can arise via two channels: one by a low-mass companion to the lens (Mao & Paczyński 1991; Gould & Loeb 1992) and the other by a faint companion to the source (Gaudi 1998). In this section, we present the analysis of the lensing light curve conducted to reveal the nature of the anomaly. Details of the analysis based on the single-lens binary-source (1L2S) and the binary-lens and single-source (2L1S) interpretations are presented in Sects. 3.1 and 3.2, respectively. The analysis under each interpretation was carried out in search for a lensing solution, which represents a set of lens- ing parameters describing the observed lensing light curve. The common lensing parameters for both 2L1S and 1L2S models are (t0, u0, tE, ρ), which represent the time of the closest lens- source approach, the projected lens-source separation scaled to the angular Einstein radius (θE) at t0 (impact parameter), the event time scale, and the source radius scaled to θE (normalized source radius), respectively. Besides these basic parameters, a 2L1S modeling requires extra parameters of (s, q, α), where the first two parameters represent the projected separation (scaled to θE) and mass ratio between the lens components, respec- tively, and the last parameter denotes the angle of the source trajectory as measured from the binary-lens axis. A 1L2S model also requires additional parameters, including (t0,2, u0,2, ρ2, qF), which represent the closest approach time, impact parameter and the normalized radius of the secondary source, and the flux ratio between the secondary and primary source stars, respectively. In the 1L2S model, we designate the time of closest approach, the impact parameter and normalized radius of the primary source as t0,1, u0,1 and ρ1, respectively, to distinguish them from those describing the secondary source. In the modeling, we take the microlens-parallax effects into consideration because the event lasted for a significant fraction of a year. For a long time-scale event like MOA-2022-BLG-249, the deviation of the source motion from rectilinear caused by the orbital motion of Earth around the sun can be substantial (Gould 1992). In order to consider these effects in the modeling, we add two extra lensing parameters (πE,N , πE,E), which denote the north and east components of the microlens-lens parallax vector πE = (πrel/θE)(µ/µ), respectively. Here, µ represents the vector of the A89, page 3 of 10 A&A 674, A89 (2023) Table 2. 1L2S model parameters. Parameter u0,1 > 0 u0,1 < 0 χ2/d.o.f. 9865.2/9800 9863.6/9800 t0,1 (HJD′) 9727.102 ± 0.007 9727.093 ± 0.006 t0,2 (HJD′) 9733.058 ± 0.009 9733.050 ± 0.008 u0,1 (10−3) 23.5 ± 0.9 −22.1 ± 0.8 u0,2 (10−3) −0.7 ± 0.2 −0.9 ± 0.2 tE (days) 143.36 ± 4.11 147.62 ± 3.95 ρ1 (10−3) 20.35 ± 1.75 18.72 ± 1.82 ρ2 (10−3) 1.54 ± 0.09 1.51 ± 0.09 πE,N −0.397 ± 0.040 −0.488 ± 0.044 πE,E 0.262 ± 0.015 0.251 ± 0.015 qF 0.0048 ± 0.0002 0.0052 ± 0.00021 Notes. HJD′ = HJD − 2450000. Fig. 2. Comparison of the 2L1S and 1L2S models. Lower two panels show the residuals from the individual models. relative lens-source proper motion and πrel denotes the relative lens-source parallax, which is related to the distance to the lens, DL, and source, DS, by πrel = AU(1/DL−1/DS). In each parallax modeling, we checked a pair of solutions with u0 > 0 and u0 < 0. 3.1. 1L2S model The 1L2S modeling was carried via a downhill approach using the Markov chain Monte Carlo (MCMC) method because the lensing magnification smoothly changes with the variation of the 1L2S lensing parameters. The initial parameters of (t0, u0, tE) were given by adopting the values obtained from the 1L1S modeling, while those related to the source companion, that is, (t0,2, u0,2, ρ2, qF), were given by considering the location and magnitude of the anomaly. We refer to Hwang et al. (2013) for details of the 1L2S modeling. The lensing parameters of the u0 > 0 and u0 < 0 solutions and their χ2 values of the fit together with the degree of freedom (d.o.f.) are listed in Table 2. It is found that the solution with u0,1 < 0 results in a slightly better fit than the solution with u0,1 > 0, by ∆χ2 = 1.6. The model curve of the u0,1 < 0 solution and its residual in the region of the anomaly are shown in Fig. 2. It is found that the 1L2S Fig. 3. Lens system configurations. Two upper panels show the configu- rations of the inner and outer 2L1S solutions with u0 > 0 and the bottom panel shows the configuration of the 1L2S solution with u0,1 < 0. In each of the panels showing the 2L1S configurations, the red cuspy fig- ures represent caustics, the line with an arrow is the source trajectory, and grey curves encompassing the caustic are equi-magnification con- tours. In the panel of the 1L2S solution, the black filled dot represent the lens, and the blue and red curves denote the trajectories of the primary (marked by S 1) and secondary (marked by S 1) source stars, respectively. models approximately delineate the observed anomaly, but they leave slight residual both in the rising and falling parts of the anomaly. In particular, the negative residuals in the rising part of the anomaly appears both in the KMTS and CFHT data sets, suggesting that these residuals are likely to be real. The lens-system configuration of the u0,1 < 0 model is shown in the bottom panel of Fig. 3, in which the arrowed curves marked in blue and red represent the trajectories of the primary (labeled as “S 1”) and secondary source (labeled as “S 2”) stars, respectively. According to the 1L2S interpretation, the anomaly was produced by the close approach of a secondary source to the lens. The secondary source is very faint, and its flux is ∼0.5% of the flux from the primary source. 3.2. 2L1S model The 2L1S modeling was conducted in two steps. In the first step, we searched for the binary-lens parameters s and q via a grid approach with multiple starting values of the source tra- jectory angle α, while we found the other lensing parameters via a downhill approach. We then constructed a ∆χ2 map on the (s, q) parameter plane and identified a pair of degenerate solutions resulting from the “inner-outer” degeneracy (Gaudi & Gould 1997). In the second step, we refined the lensing parame- ters of the individual local solutions by allowing all parameters to vary. In Table 3, we list the lensing parameters of the inner and outer 2L1S solutions, for each of which there are a pair of solu- tions with u0 > 0 and u0 < 0. Among the solutions, it was found that the inner solution with u0 > 0 yields the best fit to the data. From the comparison of the 2L1S fit with that of 1L2S fit, it is found that the anomaly is better explained by the 2L1S interpretation than the 1L2S interpretation. In Fig. 2, we draw the model curve of the inner 2L1S solution (with u0 > 0) and its residual, showing that the residual of the 1L2S model around the anomaly does not appear in the residual of the 2L1S model. From a comparison of the fits, it is found that the 2L1S model provides a better fit to the data than the 1L2S model by A89, page 4 of 10 Han, C., et al.: A&A proofs, manuscript no. aa46166-23 Table 3. Parameters of 2L1S models (parallax only). Parameter Inner Outer u0 > 0 u0 < 0 u0 > 0 u0 < 0 χ2/d.o.f. 9682.4/9799 9690.5/9799 9695.9/9799 9703.6/9799 t0 (HJD′) 9727.192 ± 0.006 9727.196 ± 0.006 9727.192 ± 0.007 9727.195 ± 0.006 u0 0.023 ± 0.001 −0.023 ± 0.001 0.023 ± 0.001 −0.022 ± 0.001 tE (days) 133.81 ± 2.97 131.90 ± 2.95 134.11 ± 2.89 133.82 ± 2.75 s 1.086 ± 0.002 1.091 ± 0.002 0.967 ± 0.003 0.961 ± 0.002 q (10−5) 7.55 ± 0.44 8.94 ± 0.52 8.01 ± 0.47 9.0 ± 0.50 α (rad) 3.638 ± 0.002 −3.599 ± 0.002 3.638 ± 0.002 −3.599 ± 0.002 ρ (10−3) <1.2 <1.2 <1.2 <1.2 πE,N −0.491 ± 0.038 −0.561 ± 0.043 −0.465 ± 0.039 −0.548 ± 0.0425 πE,E 0.260 ± 0.016 0.276 ± 0.014 0.267 ± 0.014 0.276 ± 0.0146 ∆χ2 = 181.2, indicating that the origin of the perturbation is a low-mass companion to the lens rather than a faint companion to the source. The lens-system configurations of the inner and outer 2L1S solutions, with u0 > 0 values, are shown in the two upper panels of Fig. 3. According to the inner and outer solutions, the anomaly was produced by the source passages through the regions lying on the side close to and farther from the plane- tary caustic, respectively. The inner and outer solutions can be viewed as “wide” and “close” solutions, respectively, arising due to the similarity between the central caustics induced by a wide planet and a close planet: a “close-wide” degeneracy (Griest & Safizadeh 1998). Yee et al. (2021) pointed out that the transition between the outer-inner and close-wide degeneracies is continu- ous, and Hwang et al. (2022) introduced an analytic expression for the relation between the binary separations of the inner (sin) and outer (sout) solutions: s† = (sin × sout)1/2 = (u2 anom + 4)1/2 + uanom 2 . (1) Here uanom = (τ2 anom + u2 0)1/2 represents the lens-source separa- tion at the time of the anomaly, tanom, and τanom = (tanom − t0)/tE. It is found that the value of s† estimated from the planet sepa- rations (sin, sout) = (1.086, 0.967), that is, s† = (sin × sout)1/2 = 1.024, matches the value estimated from the lensing parame- ters (t0, u0, tE, tanom) very well, that is, s† = [(u2 anom + 4)1/2 + uanom]/2 = 1.024. The estimated companion-to-primary mass ratio, q ∼ 8 × 10−5, is very low for both the inner and outer solu- tions, while the event time scale, tE ∼ 134 days, is substantially longer than the several weeks of typical Galactic lensing events (Han & Gould 2003). The normalized source radius cannot be accurately measured because the source did not cross the caus- tic and only the upper limit, ρmax = 1.2 × 10−3, can be placed. See the scatter plot of the MCMC points on the (u0, ρ) parameter plane presented in Fig. 4. We note that the degeneracy between the 2L1S and 1L2S models was able to be securely resolved thanks to the nature of the event with a high peak magnification and an acute source tra- jectory angle. In this case, the duration of the anomaly increases by a factor |1/ sinα| (Yee et al. 2021), which corresponds to a factor of 2.1 in the case of MOA-2022-BLG-249. That is to say, if this anomaly occurred with a right angle, that is, α ∼ 90◦ or 270◦, then the anomaly would have been half as short and the data might not have been good enough to distinguish the 2L1S model from the 1L2S model. For events with acute trajectory Fig. 4. Scatter plot of points in the MCMC chain on the (u0, ρ) param- eter plane obtained from the 2L1S modeling. The color coding is set to designate points with ≤1σ (red), ≤2σ (yellow), ≤3σ (green), ≤4σ (cyan), and ≤5σ (blue). angles, the magnification is lower than the peak magnification by a factor | sinα|. This factor is 0.47 (0.8 magnitudes) in the case of MOA-2022-BLG-249. 3.3. Microlens-parallax effects It has been found that considering parallax effects is important for the precise description of the observed light curve. This is somewhat expected from the long time scale of the event. The improvement of the fit with the parallax effect is huge, that is, by ∆χ2 = 5300 with respect to the model obtained under the assumption that the relative lens-source motion is rectilinear. The inner and outer 2L1S solutions result in similar values of the parallax parameters of (πE,N , πE,E) ∼ (−0.48, 0.27). We checked the solidness of the parallax measurement by inspecting the consistency of the parallax parameters measured from the 1L1S and 2L1S modeling. Figure 5 shows the scat- ter plots of MCMC points of the 1L1S solution and the inner and outer 2L1S solutions with positive u0 values. The parallax A89, page 5 of 10 A&A 674, A89 (2023) Table 4. Parameters of 2L1S models (orbit+parallax). Parameter Inner Outer u0 > 0 u0 < 0 u0 > 0 u0 < 0 χ2/d.o.f. 9682.1/9801 9690.0/9801 9682.3/9801 9686.1/9801 t0 (HJD′) 9727.191 ± 0.006 9727.198 ± 0.007 9727.188 ± 0.006 9727.190 ± 0.006 u0 0.023 ± 0.001 −0.022 ± 0.001 0.023 ± 0.001 −0.022 ± 0.001 tE (days) 134.50 ± 2.78 133.99 ± 3.12 135.56 ± 3.14 135.94 ± 2.35 s 1.073 ± 0.009 1.077 ± 0.011 0.942 ± 0.007 0.932 ± 0.009 q (10−5) 7.22 ± 0.72 8.89 ± 0.93 7.17 ± 0.69 8.07 ± 0.60 α (rad) 3.646 ± 0.015 −3.614 ± 0.016 3.642 ± 0.015 −3.615 ± 0.010 ρ (10−3) <1.2 <1.2 <1.2 <1.2 πE,N −0.464 ± 0.037 −0.557 ± 0.045 −0.474 ± 0.038 −0.529 ± 0.040 πE,E 0.270 ± 0.015 0.276 ± 0.015 0.259 ± 0.014 0.276 ± 0.014 ds/dt (yr−1) 0.705 ± 0.54 0.832 ± 0.621 1.70 ± 0.42 2.021 ± 0.554 dα/dt (yr−1) −0.58 ± 0.93 0.964 ± 1.009 −0.29 ± 0.94 0.975 ± 0.630 Fig. 5. Scatter plots of the MCMC points in the chains of the 1L1S, and the inner and outer 2L1S solutions on the (πE,E , πE,N) parameter plane. In all cases, we present plots of the solutions with u0 > 0, while the solutions with u0 < 0 result in similar plots. The dotted circles are drawn at every 0.2πE interval. The color coding is same as that used in Fig. 4. modeling was conducted by excluding the data around the per- turbation (9730 < HJD′ < 9736) because the parallactic Earth motion has a long-term effect on the lensing light curve. From this check, it was found that all the tested models result in con- sistent parallax parameters and this indicates that the parallax parameters are securely measured. We also checked the effect of the planetary orbital motion on the πE measurement because the planet might have moved during the 6-day period between the peak and the planetary per- turbation and this could affect the lens system configuration. For this check, we tested additional models considering the planetary motion by including two orbital parameters of (ds/dt, dα/dt), which represent the annual change rates of the planetary sep- aration and source trajectory angle, respectively. The lensing parameters of the solutions considering the lens-orbital motion are listed in Table 4. It is found that the lens-orbital motion does not have a significant effect on the microlens-parallax parame- ters. This can be seen in Fig. 6, where we present the scatter plots Fig. 6. Scatter plot of the MCMC chain on the parameter planes of higher-order parameters of (πE,N , πE,E , ds/dt, dα/dt) for the inner 2L1S solution with u0 > 0 obtained considering both microlens-parallax and lens-orbital effects. The plot on the (πE,E , πE,N) plane in the upper right inset is presented for the direct comparison with the plots presented in Fig. 5. of the MCMC points on the (πE,N , πE,E , ds/dt, dα/dt) parame- ter planes for the inner u0 > 0 solution, considering both the microlens-parallax and lens-orbital effects. The plots show that the uncertainties of the orbital parameters, that is, ds/dt and dα/dt, are very large. Although there are some variations of the plots in the orbital-parameter space for the other solutions, that is, the inner solution with u0 < 0 and outer solutions with u0 > 0 and u0 < 0, the variation of the parallax parameters is minor. Furthermore, the parallax parameters are similar to the values determined without considering the lens-orbital effect. These results indicate that the effect of the lens-orbital motion on the light curve is minor. We additionally checked the possibility that the parallax effect is imitated by the orbital effect induced by a source com- panion for which its luminosity contribution to the lensing light A89, page 6 of 10 Han, C., et al.: A&A proofs, manuscript no. aa46166-23 Fig. 7. Results of the xallarap modeling. Left panel shows the χ2 values of the xallarap fits as a function of the source orbital period and the right panel shows the lower limit of R = Q3/(1 + Q)2 as a function of the period. The dashed horizontal line in the left panel indicates the χ2 value of the parallax fit. curve is negligible: xallarap effects (Griest & Hu 1992; Han & Gould 1997; Smith et al. 2002). For this check, we conducted an additional modeling with the consideration of xallarap effects. Following the parameterization of Dong et al. (2009), the xal- larap modeling was done by including five extra parameters of (ξE,N , ξE,E , P, ψ, i). Here the first two parameters (ξE,N , ξE,E) denote the north and east components of the xallarap vector, ξE, respectively, and the other parameters represent the period, phase angle, and inclination of the binary-source orbit, respectively. The magnitude of the xallarap vector, ξE = (ξ2 E,N + ξ2 E,E)1/2, is related to the semi-major axis, a, of the source orbit by ξE = aS/r̂E, where r̂E = DSθE denotes the physical Einstein radius projected onto the source plane, aS = aMS,2/(MS,1 + MS,2), and (MS,1,MS,2) are the masses of source components. Combined with the Kepler’s law, the mass ratio between the source stars, Q = MS,2/MS,1, follows the relation (Dong et al. 2009) R = Q3 (1 + Q)2 = (aS/AU)3 (P/yr)2(MS,1/M�) . (2) The result of the xallarap modeling is presented in Fig. 7, in which the left and right panels show χ2 value of the xallarap fit and the lower limit of R = Q3/(1 + Q)2 value with respect to the orbital period of the source, respectively. For the computation of R, we adopted the mass of the primary source of MS ,1 ∼ 1 M� and distance to the source of DS = 8 kpc. For the angular Ein- stein radius, we adopted the lower limit of θE,min ∼ 0.46 mas because R ∝ aS ∝ θE, and thus the lower limit of the R value results from the lower limit of θE. The procedure of θE,min deter- mination is discussed in Sect. 4. From the comparison of the χ2 values between the xallarap, χ2 xallarap, and parallax, χ2 parallax, solutions, it is found that χ2 xallarap is higher than χ2 parallax for solu- tions with P < 1 yr, almost same as χ2 parallax for solutions with P ∼ 1 yr, and slightly lower than χ2 parallax for solutions with P > 1 yr. For the solutions with P > 1 yr, the χ2 difference ∆χ2 = χ2 parallax−χ 2 xallarap . 6.8 is very minor with three additional dof, and this corresponds to about 11% probability even assum- ing Gaussian statistics. Furthermore, the ratio R & 40 for these solutions and, thus, the mass ratio ratio is Q & 40, implying that mass of the source companion is MS,2 ∼ 40 M�, which corre- sponds to that of a black hole and thus unphysical. The results of the xallarap models indicate that there is no evidence that the light curve is affected by xallarap effects and, more impor- tantly, no evidence that the parallax signal is actually due to systematics. Therefore, we conclude that the parallax signal is real. As discussed in Sect. 4, the source may be a disk star lying in front of the bulge and the source distance may be smaller than the adopted value of 8 kpc, but this has little impact on the conclusions resulting from the xallarap modeling. 4. Source star and Einstein radius In this section, we define the source star of the lensing event not only for the purpose of fully characterizing the event, but also for constraining the lensing observable of the angular Einstein radius. The value of θE is estimated from the normalized source radius and angular source radius θ∗ by θE = θ∗/ρ. Although the ρ value cannot be measured for MOA-2022-BLG-249 due to the non-caustic-crossing nature of the anomaly, it is possible to con- strain its upper limit, which yields the lower limit of the Einstein radius, that is, θE,min = θ∗/ρmax. We specified the type of the source star by measuring its color and magnitude. For this specification, we first placed the source in the instrumental color-magnitude diagram (CMD) of stars around the source by measuring the V- and I-band magni- tudes of the source by regressing the light curve data measured in the individual passbands with respect to the lensing magni- fication estimated by the model. We then calibrated the source color and magnitude using the centroid of the red giant clump (RGC), for which its extinction-corrected (de-reddened) color and magnitude are known, as a reference (Yoo et al. 2004), that is, (V − I, I)0,S = (V − I, I)0,RGC + [(V − I, I)S − (V − I, I)RGC]. (3) Here, (V − I, I)S and (V − I, I)RGC denote the instrumental colors and magnitudes of the source and RGC centroid, respectively, and (V − I, I)0,S and (V − I, I)0,RGC indicate their corresponding de-reddened values. Figure 8 shows the locations of the source and RGC cen- troid in the instrumental CMD. We also marked the location of the blend. The instrumental color and magnitude are (V − I, I)S = (2.785 ± 0.017, 20.700 ± 0.001) for the source and (V − I, I)RGC = (3.285, 16.800) for the RGC centroid. With the known de-reddened color and magnitude of the RGC centroid, (V − I, I)0,RGC = (1.060, 14.530) (Bensby et al. 2013; Nataf et al. 2013), we estimated the de-reddened color and magnitude of the source of: (V − I, I)0,S = (0.560 ± 0.017, 18.430 ± 0.001). (4) According to the estimated color and magnitude, the source is a mid to late F-type main-sequence star and it probably lies in the disk in front of the bulge, although it could be a rare star in the bulge. We also estimated the de-reddened color and brightness of the blend as (V − I, I)0,b = (0.76, 19.31), assuming that the blend lies behind most of the dust, that is, in or near the bulge. We checked the possibility that the lens is the major source of the A89, page 7 of 10 A&A 674, A89 (2023) Fig. 8. Source location with respect to the red giant clump (RGC) in the instrumental CMD. The location of the blend is marked. blended flux. For this check, we measured the astrometric offset between the centroids of the source measured before and at the time of the lensing magnification. Considering that this offset is measured in the same season, it is expected that the offset would be smaller than the measurement uncertainty if the lens is the blend. The measured offset in the KMTC image is ∆θ = (169.2± 44.7) mas. This 3.8σ offset is confirmed by the offset ∆θ = (80± 10) mas measured in the CFHT data taken with seeing of 0.45′′– 0.55′′. This indicates that the blend is caused by a nearby star lying close to the source rather than the lens. With the specification of the source, we then estimate the angular radius of the source. For this, we first converted the mea- sured V − I color into V − K color with the use of the Bessell & Brett (1988) relation and then estimated the angular radius of the source using the (V − K,V)–θ∗ relation of Kervella et al. (2004). We estimated that the source has an angular radius of: θ∗ = 0.551 ± 0.040 µas, (5) and this yields the minimum values of the angular Einstein radius: θE,min = θ∗ ρmax = 0.46 mas, (6) and the relative lens-source proper motion: µmin = θE,min tE = 1.25 mas yr−1. (7) 5. Physical lens parameters The physical parameters of a lens are constrained by measuring the lensing observables of an event. These observables include the event time scale, tE, Einstein radius, θE, and microlens parallax vector, πE = (πE,N , πE,E), and the mass and distance to the lens are determined from the combination of these observables as M = θE κπE ; DL = AU πEθE + πS , (8) respectively (Gould 2000). Here κ = 4G/(c2AU) and πS = AU/DS is the parallax of the source. For MOA-2022-BLG- 249, the values of tE and πE are securely measured, but the value of θE cannot be measured and only its lower limit is constrained, making it difficult to analytically estimate M and DL using the relations in Eq. (8). We, therefore, estimate the physical lens parameters by conducting a Bayesian analysis based on the measured lensing observables and other available constraints. In the first step of the Bayesian analysis, we generated a large number (2 × 108) of artificial lensing events, for which the loca- tions the lens and source and their relative proper motion were assigned on the basis of a Galactic model and the lens masses were allocated on the basis of a mass-function model by conduct- ing a Monte Carlo simulation. In the simulation, we adopted the models of the Galaxy and lens mass function described in Jung et al. (2021) and Jung et al. (2018), respectively. For each simu- lated event, we computed the lensing observables corresponding to the values of M, DL, DS, and µ by tE = θE µ ; θE = (κMπrel)1/2; πE = πrel θE . (9) In the second step, we imposed a weight wi = exp(−χ2/2) to each artificial event and constructed posteriors of M and DL. In this procedure, the χ2 value is calculated as χ2 = ( tE,i − tE σtE )2 + 2∑ j,k=1 b jk(πE, j,i − πE,i)(πE,k,i − πE,i), (10) where (πE,1, πE,2)i = (πE,N , πE,E)i is expressed in two compo- nent form, (tE,i,πE,i) are the observables of each simulated event, (tE,πE) represent the measured observables, σtE, is the uncer- tainty in the tE measurement, and b jk is the inverse covariance matrix of πE. We refer to Eqs. (10) and (11) of Gould et al. (2022). Finally, we imposed the constraint of the Einstein radius by setting wi = 0 for events with θE ≤ θE,min. We imposed an additional constraint provided by the fact that the flux from the lens cannot be greater than the blend flux. Imposing this blend-flux constraint may be important because the distance to the lens expected from the large value of the measured microlens parallax, πE = (π2 E,N + π2 E,E)1/2 ∼ 0.55, is small. In order to impose this constraint, we computed the lens magnitude as: IL = MI,L + 5 log ( DL pc ) − 5 + AI,L, (11) where MI,L is the absolute I-band magnitude of a star cor- responding to the lens mass, and AI,L represents the I-band extinction to the lens. For the computation of AI,L, we modeled the extinction as AI,L = AI,tot [ 1 − exp ( − |z| hz,dust )] , (12) where AI,tot = 2.49 is the total I-band extinction toward the field, hz,dust = 100 pc is the vertical scale height of dust, z = DL sin b + z0, b is the Galactic latitude, and z0 = 15 pc is vertical position of the Sun above the Galactic plane. Figure 9 shows the Bayesian posteriors of the host mass (top panel), Mhost, distances to the planetary system (middle panel), and source (bottom panel). We present two sets of pos- terior: one with (shaded distribution) and the other without A89, page 8 of 10 Han, C., et al.: A&A proofs, manuscript no. aa46166-23 Fig. 9. Bayesian posteriors of the lens mass, distance to the lens and source. In each panel, the shaded and and unshaded distributions are obtained with (shaded) and without (unshaded) imposing the blend-flux constraint, respectively. Table 5. Physical lens parameters. Parameter Inner Outer Mhost (M�) 0.18 ± 0.05 ← Mplanet (M⊕) 4.83 ± 1.44 ← DL (kpc) 2.00 ± 0.42 ← a⊥ (AU) 1.63 ± 0.35 1.45 ± 0.31 (unshaded distribution) the blend-flux constraint. The poste- rior distributions show that the physical parameters are tightly defined despite the limited information on the angular Einstein radius. In Table 5, we summarize the estimated physical param- eters, in which the median values are presented as representative values and the uncertainties are estimated as 16% and 84% of the Bayesian posterior distributions. Here, the planet mass is estimated as Mplanet = qMhost and the projected planet-host sep- aration is computed by a⊥ = sθEDL. From the fact that the lower mass limit and the upper distance limit estimated using the analytic relations in Eq. (8) based on the lower limit of the Einstein radius, that is, Mmin = θE,min/κπE ∼ 0.12 M� and DL,max = AU/(πEθE,min + πS) ∼ 2.4 kpc, match well the corre- sponding limits of the Bayesian posteriors indicates that these limits of the physical parameters are set by the combination of the constraints provided by θE,min and πE. On the other hand, the upper limit of the mass and lower limit of the distance are set by the blend-flux constraint. This can be seen from the comparison of the Bayesian posteriors obtained, with and without imposing the lens flux constraint. It turns out that the lens is a planetary system, in which a low-mass planet orbits a low-mass host star lying in the Galactic disk. The estimated mass of the planet, Mplanet ∼ 4.8 M⊕, indi- cates that the planet is a super-Earth, and the detection of the system aptly demonstrates the elevated microlensing sensitiv- ity to low-mass planets with the increase of the observational cadence. The estimated mass of the host, Mhost ∼ 0.18 M�, and distance, DL ∼ 2.0 kpc, indicate that the host of the planet is a very low-mass M dwarf lying in the Galactic disk. Finding planets belonging to such low-mass stars using other methods is difficult because of the faintness of host stars, and thus the dis- covered planetary system well demonstrates the usefulness of the microlensing method in finding planets with low-mass host stars. The planetary system lies at a substantially closer distance than those of typical microlensing planets, which usually lie either in the bulge or in the portion of the disk that is closer to the bulge than to the Sun. We checked the hypothesis that the source is in the disk by additionally conducting a Bayesian analysis, in which we assumed that DS = 7 kpc and the dispersion of the source motion is negligible, as in the case of disk stars. We found that this analysis results in similar posteriors to those presented in Fig. 9, indicating that the uncertain source location has little effect on the result. 6. Summary and conclusion We analyzed the microlensing event MOA-2022-BLG-249, for which the light curve exhibits a brief positive anomaly with a duration of ∼1 day and a maximum deviation of ∼0.2 mag from a single-source single-lens model. We tested both the planetary and binary-source origins, which are the two channels producing a short-term positive anomaly in a lensing light curve. We found that the anomaly was produced by a planetary companion to the lens, rather than a binary companion to the source. We identified two solutions rooted in the inner-outer degeneracy, for both of which the estimated planet-to-host mass ratio, q ∼ 8 × 10−5, is very small. With the constraints pro- vided by the microlens parallax, along with the lower limit of the Einstein radius taken together with the blend-flux constraint, we find that the lens is a planetary system, in which a super- Earth planet, with a mass of (4.83 ± 1.44) M⊕, orbits a low-mass host star, with a mass of (0.18 ± 0.05) M�, lying in the Galactic disk at a distance (2.00 ± 0.42) kpc. The planet detection aptly demonstrates the elevated microlensing sensitivity of the current high-cadence lensing surveys to low-mass planets. Acknowledgements. Work by C.H. was supported by the grants of National Research Foundation of Korea (2020R1A4A2002885 and 2019R1A2C2085965). This research has made use of the KMTNet system operated by the Korea Astron- omy and Space Science Institute (KASI) and the data were obtained at three host sites of CTIO in Chile, SAAO in South Africa, and SSO in Australia. This research was supported by the Korea Astronomy and Space Science Institute under the R&D program (Project No. 2023-1-832-03) supervised by the Min- istry of Science and ICT. The MOA project is supported by JSPS KAKENHI Grant Number JSPS24253004, JSPS26247023, JSPS23340064, JSPS15H00781, JP16H06287, and JP17H02871. J.C.Y., I.G.S., and S.J.C. acknowledge support from NSF Grant No. AST-2108414. Y.S. acknowledges support from NSF Grant No. 2020740. This research uses data obtained through the Telescope Access Program (TAP), which has been funded by the TAP member institutes. W. Zang, H.Y., S.M., and W. Zhu acknowledge support by the National Science Founda- tion of China (Grant No. 12133005). W. Zang acknowledges the support from the Harvard-Smithsonian Center for Astrophysics through the CfA Fellowship. W. Zhu acknowledges the science research grants from the China Manned Space Project with No. CMS-CSST-2021-A11. C.R. was supported by the Research fellowship of the Alexander von Humboldt Foundation. References Beaulieu, J. -P., Bennett, D. P., Fouqué, P., et al. 2006, Nature, 439, 437 Bennett, D. P., & Rhie, S. H. 1996, ApJ, 472, 660 Bennett, D. P., Bond, I. A., Udalski, A., et al. 2008, ApJ, 684, 663 A89, page 9 of 10 http://linker.aanda.org/10.1051/0004-6361/202346166/1 http://linker.aanda.org/10.1051/0004-6361/202346166/2 http://linker.aanda.org/10.1051/0004-6361/202346166/3 A&A 674, A89 (2023) Bennett, D. P., Batista, V., Bond, I. A., et al. 2014, ApJ, 785, 155 Bensby, T. Yee, J.C., Feltzing, S., et al. 2013, A&A, 549, A147 Bessell, M. S., & Brett, J. M. 1988, PASP, 100, 1134 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 Dong, S., Gould, A., Udalski, A., et al. 2009, ApJ, 695, 970 Gaudi, B. S. 1998, ApJ, 506, 533 Gaudi, B. S. 2012, ARA&A, 50, 411 Gaudi, B. S., & Gould, A. 1997, ApJ, 486, 85 Gaudi, B. S., & Han, C. 2004, ApJ, 611, 528 Gould, A. 1992, ApJ, 392, 442 Gould, A. 2000, ApJ, 542, 785 Gould, A., & Loeb, A. 1992, ApJ, 396, 104 Gould, A., Udalski, A., Shin, I.-G., et al. 2014, Science, 345, 46 Gould, A., Ryu, Y.-H., Calchi Novati, S., et al. 2020, JKAS, 53, 9 Gould, A., Han, C., Weicheng, Z., et al. 2022, A&A, 664, A13 Griest, K., & Hu, W. 1992, ApJ, 397, 362 Griest, K., & Safizadeh, N. 1998, ApJ, 500, 37 Han, C. 2006, ApJ, 638, 1080 Han, C., & Gould, A. 1997, ApJ, 480, 196 Han, C., & Gould, A. 2003, ApJ, 502, 172 Han, C., Hirao, Y., Udalski, A., et al. 2018, AJ, 155, 211 Han, C., Udalski, A., Lee, C.-U., et al. 2021, A&A, 649, A90 Han, C., Bond, I. A., Yee, J. C., et al. 2022a, A&A, 658, A94 Han, C., Gould, A., Albrow, M. D., et al. 2022b, A&A, 658, A62 Herrera-Martín, A., Albrow, M. D., Udalski, A., et al. 2020, AJ, 159, 256 Hwang, K.-H., Choi, J.-Y., Bond, I. A., et al. 2013, ApJ, 778, 55 Hwang, K.-H., Zang, W., Gould, A., et al. 2022, AJ, 163, 43 Jung, Y. K., Udalski, A., Gould, A., et al. 2018, AJ, 155, 219 Jung, Y. K., Han, C., Udalski, A., et al. 2021, AJ, 161, 293 Kervella, P., Thévenin, F., Di Folco, E., & Ségransan, D. 2004, A&A, 426, 29 Kim, S.-L., Lee, C.-U., Park, B.-G., et al. 2016, JKAS, 49, 37 Kondo, I., Yee, J. C., Bennett, D. P., et al. 2021, AJ, 162, 77 Mao, S., & Paczyński, B. 1991, ApJ, 374, 37 Mróz, P., Poleski, R., Gould, A., et al. 2020, ApJ, 903, L11 Muraki, Y., Han, C., Bennett, D. P., et al. 2011, ApJ, 741, 22 Nataf, D. M., Gould, A., Fouqué, P., et al. 2013, ApJ, 769, 88 Ryu, Y.-H., Udalski, A., Yee, J. C., et al. 2020, AJ, 160, 183 Ryu, Y.-H., Jung, Y. K., Yang, H., et al. 2022, AJ, 164, 180 Shvartzvald, Y., Yee, J. C., Calchi Novati, S., et al. 2017, ApJ, 840, L3 Smith, M. C., Mao, S., Woźniak, P., et al. 2002, MNRAS, 336, 670 Sumi, T., Udalski, A., Bennett, D. P., et al. 2016, ApJ, 825, 112 Udalski, A., Szymański, M. K., & Szymański, G. 2015, Acta Astron., 65, 1 Yee, J. C., Shvartzvald, Y., Gal-Yam, A., et al. 2012, ApJ, 755, 102 Yee, J. C., Zang, W., Udalski, A., et al. 2021, AJ, 162, 180 Yoo, J., DePoy, D. L., Gal-Yam, A., et al. 2004, ApJ, 603, 139 Vandorou, A., Dang, L., Bennett, D. P., et al. 2023, AJ, submitted [arXiv:2302.01168] Zang, W., Han, C., Kondo, I., et al. 2021a, RAA, 21, 239 Zang, W., Hwang, K.-H., Udalski, A., et al. 2021b, AJ, 162, 163 Zang, W., Jung, Y. K., Yang, H., et al. 2023, AJ, 165, 103 1 Department of Physics, Chungbuk National University, Cheongju 28644, Republic of Korea e-mail: cheongho@astroph.chungbuk.ac.kr 2 Max-Planck-Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany 3 Department of Astronomy, Ohio State University, 140 W. 18th Ave., Columbus, OH 43210, USA 4 Korea Astronomy and Space Science Institute, Daejon 34055, Republic of Korea 5 Korea University of Science and Technology, Korea, (UST), 217 Gajeong-ro, Yuseong-gu, Daejeon, 34113, Republic of Korea 6 Institute of Natural and Mathematical Science, Massey University, Auckland 0745, New Zealand 7 Center for Astrophysics, Harvard & Smithsonian, 60 Garden St., Cambridge, MA 02138, USA 8 Department of Astronomy, Tsinghua University, Beijing 100084, PR China 9 University of Canterbury, Department of Physics and Astronomy, Private Bag 4800, Christchurch 8020, New Zealand 10 Department of Particle Physics and Astrophysics, Weizmann Insti- tute of Science, Rehovot 76100, Israel 11 School of Space Research, Kyung Hee University, Yongin, Kyeonggi 17104, Republic of Korea 12 Institute for Space-Earth Environmental Research, Nagoya Univer- sity, Nagoya 464-8601, Japan 13 Code 667, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 14 Department of Astronomy, University of Maryland, College Park, MD 20742, USA 15 Department of Earth and Planetary Science, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 16 Instituto de Astrofísica de Canarias, Vía Láctea s/n, 38205 La Laguna, Tenerife, Spain 17 Department of Earth and Space Science, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan 18 Department of Physics, The Catholic University of America, Washington, DC 20064, USA 19 Department of Astronomy, Graduate School of Science, The Uni- versity of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 20 Sorbonne Université, CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France 21 Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand 22 University of Canterbury Mt. John Observatory, PO Box 56, Lake Tekapo 8770, New Zealand A89, page 10 of 10 http://linker.aanda.org/10.1051/0004-6361/202346166/4 http://linker.aanda.org/10.1051/0004-6361/202346166/5 http://linker.aanda.org/10.1051/0004-6361/202346166/6 http://linker.aanda.org/10.1051/0004-6361/202346166/7 http://linker.aanda.org/10.1051/0004-6361/202346166/8 http://linker.aanda.org/10.1051/0004-6361/202346166/8 http://linker.aanda.org/10.1051/0004-6361/202346166/9 http://linker.aanda.org/10.1051/0004-6361/202346166/10 http://linker.aanda.org/10.1051/0004-6361/202346166/11 http://linker.aanda.org/10.1051/0004-6361/202346166/12 http://linker.aanda.org/10.1051/0004-6361/202346166/13 http://linker.aanda.org/10.1051/0004-6361/202346166/14 http://linker.aanda.org/10.1051/0004-6361/202346166/15 http://linker.aanda.org/10.1051/0004-6361/202346166/16 http://linker.aanda.org/10.1051/0004-6361/202346166/17 http://linker.aanda.org/10.1051/0004-6361/202346166/18 http://linker.aanda.org/10.1051/0004-6361/202346166/19 http://linker.aanda.org/10.1051/0004-6361/202346166/20 http://linker.aanda.org/10.1051/0004-6361/202346166/21 http://linker.aanda.org/10.1051/0004-6361/202346166/22 http://linker.aanda.org/10.1051/0004-6361/202346166/23 http://linker.aanda.org/10.1051/0004-6361/202346166/24 http://linker.aanda.org/10.1051/0004-6361/202346166/25 http://linker.aanda.org/10.1051/0004-6361/202346166/26 http://linker.aanda.org/10.1051/0004-6361/202346166/27 http://linker.aanda.org/10.1051/0004-6361/202346166/28 http://linker.aanda.org/10.1051/0004-6361/202346166/29 http://linker.aanda.org/10.1051/0004-6361/202346166/30 http://linker.aanda.org/10.1051/0004-6361/202346166/31 http://linker.aanda.org/10.1051/0004-6361/202346166/32 http://linker.aanda.org/10.1051/0004-6361/202346166/33 http://linker.aanda.org/10.1051/0004-6361/202346166/34 http://linker.aanda.org/10.1051/0004-6361/202346166/35 http://linker.aanda.org/10.1051/0004-6361/202346166/36 http://linker.aanda.org/10.1051/0004-6361/202346166/37 http://linker.aanda.org/10.1051/0004-6361/202346166/38 http://linker.aanda.org/10.1051/0004-6361/202346166/39 http://linker.aanda.org/10.1051/0004-6361/202346166/40 http://linker.aanda.org/10.1051/0004-6361/202346166/41 http://linker.aanda.org/10.1051/0004-6361/202346166/42 http://linker.aanda.org/10.1051/0004-6361/202346166/43 http://linker.aanda.org/10.1051/0004-6361/202346166/44 http://linker.aanda.org/10.1051/0004-6361/202346166/45 http://linker.aanda.org/10.1051/0004-6361/202346166/46 http://linker.aanda.org/10.1051/0004-6361/202346166/47 http://linker.aanda.org/10.1051/0004-6361/202346166/48 http://linker.aanda.org/10.1051/0004-6361/202346166/49 https://arxiv.org/abs/2302.01168 http://linker.aanda.org/10.1051/0004-6361/202346166/51 http://linker.aanda.org/10.1051/0004-6361/202346166/52 http://linker.aanda.org/10.1051/0004-6361/202346166/53 mailto:cheongho@astroph.chungbuk.ac.kr MOA-2022-BLG-249Lb: Nearby microlensing super-Earth planet detected from high-cadence surveys 1 Introduction 2 Observations and data 3 Light curve analysis 3.1 1L2S model 3.2 2L1S model 3.3 Microlens-parallax effects 4 Source star and Einstein radius 5 Physical lens parameters 6 Summary and conclusion Acknowledgements References