pith. sign in

arxiv: 2605.17654 · v1 · pith:JTMJJHHBnew · submitted 2026-05-17 · 🌌 astro-ph.SR

Four-Body Gravitational Microlensing Events Involving Both a Binary Lens and a Binary Source

Cheongho Han , Chung-Uk Lee , Andrzej Udalski , Michael D. Albrow , Sun-Ju Chung , Andrew Gould , Youn Kil Jung , Kyu-Ha~Hwang
show 26 more authors
Yoon-Hyun Ryu Yossi Shvartzvald In-Gu Shin Jennifer C. Yee Weicheng Zang Hongjing Yang Doeon Kim Dong-Jin Kim Seung-Lee Kim Dong-Joo Lee Sang-Mok Cha Yongseok Lee Byeong-Gon Park Richard W. Pogge Przemek Mr\'oz Micha{\l} K. Szyma\'nskim Jan Skowron Rados{\l}aw Poleski Igor Soszy\'nski Pawe{\l} Pietrukowicz Szymon Koz{\l}owski Krzysztof A. Rybicki Patryk Iwanek Krzysztof Ulaczyk Marcin Wrona Mariusz Gromadzki Mateusz J. Mr\'oz
This is my paper

Pith reviewed 2026-05-19 22:06 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords gravitational microlensingbinary lensbinary sourcefour-body eventscaustic crossingsKMTNetOGLElight curve analysis
0
0 comments X

The pith

Three microlensing events are best fit by four-body models with both binary lenses and binary sources.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper re-examines three anomalous microlensing events from KMTNet and OGLE data that previous models could not fully explain. It demonstrates that configurations with a binary lens and a binary source, known as 2L2S, account for the full light-curve structures including weak residuals, late re-brightenings, and isolated anomalies. These models outperform simpler 2L1S or 3L1S alternatives. The sources are identified as pairs of main-sequence stars with specific spectral types, while the lenses appear to be low-mass binaries, one possibly including a brown dwarf. This work emphasizes the need for complex modeling in preparation for precise future observations.

Core claim

The authors show that the events KMT-2021-BLG-0209, KMT-2021-BLG-0901, and OGLE-2025-BLG-0356 can be described by four-body gravitational microlensing involving a binary lens and a binary source. In each case, the secondary source interacts with the lens caustics to produce the observed features that simpler models left as residuals.

What carries the argument

The 2L2S model, which incorporates caustic crossings by both components of the binary source in addition to the binary lens.

If this is right

  • The lenses are low-mass binaries, with one system having a brown-dwarf companion.
  • Source stars are binary pairs of types (G8V, M3V), (G8V, K2V), and (G6V, G8V).
  • These add to the sample of characterized 2L2S microlensing events.
  • Systematic testing of multi-body models will be essential for Roman Space Telescope data.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Many other microlensing anomalies might be reinterpreted as binary sources rather than more complex lens structures.
  • High-precision photometry from space telescopes could reveal even more such four-body events.
  • Binary sources may be more common in microlensing samples than currently assumed.

Load-bearing premise

The light-curve anomalies are due to caustic interactions with a faint secondary source rather than unmodeled systematics or alternative configurations.

What would settle it

Spectroscopic or high-resolution imaging confirmation that the source stars are not binaries with the predicted colors and magnitudes would disprove the 2L2S models for these events.

Figures

Figures reproduced from arXiv: 2605.17654 by Andrew Gould, Andrzej Udalski, Byeong-Gon Park, Cheongho Han, Chung-Uk Lee, Doeon Kim, Dong-Jin Kim, Dong-Joo Lee, Hongjing Yang, Igor Soszy\'nski, In-Gu Shin, Jennifer C. Yee, Krzysztof A. Rybicki, Krzysztof Ulaczyk, Kyu-Ha~Hwang, Marcin Wrona, Mariusz Gromadzki, Mateusz J. Mr\'oz, Michael D. Albrow, Micha{\l} K. Szyma\'nskim Jan Skowron, Patryk Iwanek, Pawe{\l} Pietrukowicz, Przemek Mr\'oz, Rados{\l}aw Poleski, Richard W. Pogge, Sang-Mok Cha, Seung-Lee Kim, Sun-Ju Chung, Szymon Koz{\l}owski, Weicheng Zang, Yongseok Lee, Yoon-Hyun Ryu, Yossi Shvartzvald, Youn Kil Jung.

Figure 1
Figure 1. Figure 1: Light curve of the microlensing event KMT-2021-BLG-0209. The lower panel shows the full light curve, while the upper panel presents a zoomed-in view of the caustic-crossing region. The dotted and solid curves represent the best-fit models based on the 2L1S and 2L2S configurations, respectively. The two insets in the lower panel display the source trajectories relative to the caustic structure. The left ins… view at source ↗
Figure 2
Figure 2. Figure 2: Lensing light curve, system configurations, and model curves of the lensing event KMT-2021-BLG-0901. Notations are similar to those of [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Light curve of the lensing event OGLE-2025-BLG-0356. The lower panel shows two sets of lens-system configurations: the 2L2S configuration on the left and the 3L1S configuration on the right. In the right 2L2S inset, the red and blue open circles indicate the positions of the binary source stars at t1. The circles are not drawn to scale with respect to the source size. well reproduced by a 2L1S model, where… view at source ↗
Figure 4
Figure 4. Figure 4: Locations of the binary-source stars in the color–magnitude diagram (CMD). The grey and brown CMDs are constructed from KMTC and HST observations, respectively. plotted on the data as a dotted curve in [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Bayesian posterior distribution of the primary lens mass. The blue and red curves show the contributions from the disk and bulge lens populations, respectively. The solid vertical line indicates the median of the posterior, and the dotted lines indicate the 1σ range. observables with prior knowledge of the Galactic pop￾ulations that can act as lenses and sources, including their spatial distributions and k… view at source ↗
read the original abstract

We present detailed analyses of three anomalous microlensing events--KMT-2021-BLG-0209, KMT-2021-BLG-0901, and OGLE-2025-BLG-0356--identified from a systematic re-examination of KMTNet light curves for which previous modeling attempts failed or left persistent residuals. Although all three events show caustic-related features consistent with binary-lens microlensing, we find that their full light-curve structures can be described by four-body configurations that required four-body configurations involving a binary lens and a binary source. In KMT-2021-BLG-0209, weak caustic-exit residuals arise from a faint companion source undergoing an additional caustic interaction. In KMT-2021-BLG-0901, a late-time re-brightening is produced when the secondary source encounters the resonant caustic long after the primary. For OGLE-2025-BLG-0356, we test the degeneracy between 3L1S and 2L2S interpretations of a short isolated anomaly and find that the 2L2S model provides a significantly better fit. Source colors and magnitudes indicate binary sources composed of (G8V, M3V), (G8V, K2V), and (G6V, G8V) stars for the three events, respectively. Bayesian inference suggests that the lenses are predominantly low-mass binaries, including one system (KMT-2021-BLG-0901) with a companion consistent with a brown dwarf. These events add to the growing sample of well-characterized 2L2S systems and underscore the importance of systematically testing complex models, particularly in anticipation of the high-precision microlensing data expected from the Roman Space Telescope survey.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The paper presents detailed modeling of three anomalous microlensing events (KMT-2021-BLG-0209, KMT-2021-BLG-0901, OGLE-2025-BLG-0356) previously showing persistent residuals under simpler models. It claims that four-body 2L2S configurations (binary lens plus binary source) fully describe the light curves, with the secondary source explaining specific features such as weak caustic-exit residuals, late-time re-brightening, and a short isolated anomaly. The 2L2S models are reported to yield significantly better fits than 2L1S or 3L1S alternatives; source colors/magnitudes constrain the binaries to (G8V, M3V), (G8V, K2V), and (G6V, G8V) pairs; and Bayesian inference indicates predominantly low-mass lenses, including one possible brown-dwarf companion.

Significance. If the statistical preference for the 2L2S solutions holds after proper accounting for extra parameters, the work adds three well-characterized 2L2S events to the literature and illustrates the value of systematically testing higher-order models ahead of Roman Space Telescope data. The use of source color/magnitude constraints and Bayesian lens inference on standard microlensing equations is a strength, but the overall impact is tempered by the absence of quantified fit improvements in the provided material.

major comments (1)
  1. [Abstract and modeling sections] Abstract and modeling sections: The central claim that 2L2S models 'provide significantly better fits' than 3L1S or 2L1S alternatives is load-bearing for the interpretation of the residuals as secondary-source caustic crossings. However, no Δχ² values, degrees-of-freedom changes, or reduced-χ² comparisons are reported despite the addition of at least four extra parameters (flux ratio, source separation vector, and possible independent proper motion). This prevents assessment of whether the improvement exceeds the Occam penalty and leaves open the possibility that residuals arise from unmodeled systematics or higher-order effects rather than a physical binary source.
minor comments (1)
  1. [Abstract] The abstract refers to 'four-body configurations that required four-body configurations'—a minor phrasing redundancy that should be cleaned up for clarity.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful and constructive review. We agree that the absence of explicit quantitative comparisons of fit quality is a limitation in the submitted manuscript and will address this directly in revision. Below we respond to the major comment.

read point-by-point responses

  1. Referee: [Abstract and modeling sections] Abstract and modeling sections: The central claim that 2L2S models 'provide significantly better fits' than 3L1S or 2L1S alternatives is load-bearing for the interpretation of the residuals as secondary-source caustic crossings. However, no Δχ² values, degrees-of-freedom changes, or reduced-χ² comparisons are reported despite the addition of at least four extra parameters (flux ratio, source separation vector, and possible independent proper motion). This prevents assessment of whether the improvement exceeds the Occam penalty and leaves open the possibility that residuals arise from unmodeled systematics or higher-order effects rather than a physical binary source.

    Authors: We accept this criticism. The submitted manuscript does not report the requested Δχ², Δdof, or reduced-χ² values, which prevents a transparent evaluation of whether the improvement justifies the added parameters. In the revised version we will add a table (or subsection) in the modeling section for each event that lists the minimum χ² for the best 2L1S, 3L1S, and 2L2S solutions, the number of additional free parameters in the 2L2S model, and the resulting Δχ² per degree of freedom. We will also include a brief discussion of whether the improvement exceeds a simple Occam penalty. The physical case for the 2L2S interpretation rests on both the statistical improvement and the ability of the secondary source to account for specific, localized residuals (weak caustic-exit feature, late-time re-brightening, isolated anomaly) while remaining consistent with the measured source colors and magnitudes. Nevertheless, we agree that the quantitative statistics must be shown explicitly before the claim of a 'significantly better fit' can be fully assessed. revision: yes

Circularity Check

0 steps flagged

No significant circularity: standard microlensing model fitting with external priors

full rationale

The paper performs direct light-curve modeling of three events using established 2L2S, 3L1S, and 2L1S configurations drawn from standard microlensing equations. Model preference is assessed via fit quality to observed data, with source colors/magnitudes and Bayesian inference relying on external stellar population priors rather than quantities defined solely by the paper's own fitted parameters. No step reduces by construction to a self-definition, fitted input renamed as prediction, or load-bearing self-citation chain; the central claim that 2L2S models describe the full curves better is a data-driven comparison, not an internal tautology. The derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The modeling rests on standard point-mass lens and rectilinear-motion assumptions plus many free parameters fitted directly to the light curves; no new entities are postulated.

free parameters (2)
  • lens masses, separations, and velocities
    Multiple parameters fitted to reproduce the observed light-curve features for each event.
  • source magnitudes and colors
    Derived from photometry to classify the binary sources.
axioms (1)
  • standard math Standard gravitational microlensing assumptions of point-mass lenses and straight-line relative motion
    Invoked throughout the light-curve modeling described in the abstract.

pith-pipeline@v0.9.0 · 6051 in / 1409 out tokens · 43679 ms · 2026-05-19T22:06:27.320086+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

29 extracted references · 29 canonical work pages · 4 internal anchors

  1. [1]

    Alard, C., & Lupton, R. H. 1998, ApJ, 503, 325, doi: 10.1086/305984

    work page doi:10.1086/305984 1998

  2. [2]

    R., 2009, @doi [ ] 10.1111/j.1365-2966.2009.15167.x , http://adsabs.harvard.edu/abs/2009MNRAS.398..607V 398, 607

    Albrow, M. D., Horne, K., Bramich, D. M., et al. 2009, MNRAS, 397, 2099, doi: 10.1111/j.1365-2966.2009.15098.x 12

    work page doi:10.1111/j.1365-2966.2009.15098.x 2009

  3. [3]

    W., Allsman, R

    Alcock, C., Akerlof, C. W., Allsman, R. A., et al. 1993, Nature, 365, 621, doi: 10.1038/365621a0

    work page doi:10.1038/365621a0 1993

  4. [4]

    1993, Nature, 365, 623, doi: 10.1038/365623a0

    Aubourg, E., Bareyre, P., Br´ ehin, S., et al. 1993, Nature, 365, 623, doi: 10.1038/365623a0

    work page doi:10.1038/365623a0 1993

  5. [5]

    C., Feltzing, S., et al

    Bensby, T., Yee, J. C., Feltzing, S., et al. 2013, A&A, 549, A247, doi: 10.1051/0004-6361/201220678

    work page doi:10.1051/0004-6361/201220678 2013

  6. [6]

    1991, A&A, 248, 485

    Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485

    work page 1991

  7. [7]

    2003, ApJ, 592, 172, doi: 10.1086/375706

    Han, C., & Gould, A. 2003, ApJ, 592, 172, doi: 10.1086/375706

    work page doi:10.1086/375706 2003

  8. [8]

    2025a, PASP, 137, 104501, doi: 10.1088/1538-3873/ae06c8

    Han, C., Lee, C.-U., Udalski, A., et al. 2025a, PASP, 137, 104501, doi: 10.1088/1538-3873/ae06c8

    work page doi:10.1088/1538-3873/ae06c8

  9. [9]

    A., et al

    Han, C., Udalski, A., Bond, I. A., et al. 2024, A&A, 692, A221, doi: 10.1051/0004-6361/202451806

    work page doi:10.1051/0004-6361/202451806 2024

  10. [10]

    K., et al

    Han, C., Udalski, A., Jung, Y. K., et al. 2023, A&A, 670, A172, doi: 10.1051/0004-6361/202245525

    work page doi:10.1051/0004-6361/202245525 2023

  11. [11]

    2021, AJ, 161, 270, doi: 10.3847/1538-3881/abf4d0

    Han, C., Udalski, A., Kim, D., et al. 2021, AJ, 161, 270, doi: 10.3847/1538-3881/abf4d0

    work page doi:10.3847/1538-3881/abf4d0 2021

  12. [12]

    A., Watson, A

    Holtzman, J. A., Watson, A. M., Baum, W. A., et al. 1998, AJ, 115, 1946, doi: 10.1086/300336

    work page doi:10.1086/300336 1998

  13. [13]

    K., Han, C., Udalski, A., et al

    Jung, Y. K., Han, C., Udalski, A., et al. 2021, AJ, 161, 293, doi: 10.3847/1538-3881/abf8bd

    work page doi:10.3847/1538-3881/abf8bd 2021

  14. [14]

    K., Udalski, A., Bond, I

    Jung, Y. K., Udalski, A., Bond, I. A., et al. 2017, ApJ, 841, 75, doi: 10.3847/1538-4357/aa7057

    work page doi:10.3847/1538-4357/aa7057 2017

  15. [15]

    K., Zang, W., Han, C., et al

    Jung, Y. K., Zang, W., Han, C., et al. 2022, AJ, 164, 262, doi: 10.3847/1538-3881/ac9c5c

    work page doi:10.3847/1538-3881/ac9c5c 2022

  16. [16]

    2004, A&A, 426, 29, doi: 10.1051/0004-6361:20035930

    Kervella, P., Th´ evenin, F., Di Folco, E., & S´ egransan, D. 2004, A&A, 426, 29, doi: 10.1051/0004-6361:20035930

    work page doi:10.1051/0004-6361:20035930 2004

  17. [17]

    2016, Journal of Korean Astronomical Society, 49, 37, doi: 10.5303/JKAS.2016.49.1.37

    Kim, S.-L., Lee, C.-U., Park, B.-G., et al. 2016, JKAS, 49, 37, doi: 10.5303/JKAS.2016.49.1.37

    work page doi:10.5303/jkas.2016.49.1.37 2016

  18. [18]

    M., Gould, A., Fouqu´ e, P., et al

    Nataf, D. M., Gould, A., Fouqu´ e, P., et al. 2013, ApJ, 769, 88, doi: 10.1088/0004-637X/769/2/88

    work page doi:10.1088/0004-637x/769/2/88 2013

  19. [19]

    T., Gaudi, B

    Penny, M. T., Gaudi, B. S., Kerins, E., et al. 2019, ApJS, 241, 34, doi: 10.3847/1538-4365/aafb69

    work page doi:10.3847/1538-4365/aafb69 2019

  20. [20]

    A., Henry, T

    Raghavan, D., McAlister, H. A., Henry, T. J., et al. 2010, ApJS, 190, 1, doi: 10.1088/0067-0049/190/1/1

    work page internal anchor Pith review doi:10.1088/0067-0049/190/1/1 2010

  21. [21]

    Wide-Field InfrarRed Survey Telescope-Astrophysics Focused Telescope Assets WFIRST-AFTA 2015 Report

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

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1503.03757 2015

  22. [22]

    P., Bond, I

    Sumi, T., Bennett, D. P., Bond, I. A., et al. 2013, ApJ, 778, 150, doi: 10.1088/0004-637X/778/2/150

    work page doi:10.1088/0004-637x/778/2/150 2013

  23. [23]

    B., & Crotts, A

    Tomaney, A. B., & Crotts, A. P. S. 1996, AJ, 112, 2872, doi: 10.1086/118228

    work page doi:10.1086/118228 1996

  24. [24]

    The Optical Gravitational Lensing Experiment. Real Time Data Analysis Systems in the OGLE-III Survey

    Udalski, A. 2003, Acta Astronomica, 53, 291, doi: 10.48550/arXiv.astro-ph/0401123

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/0401123 2003

  25. [25]

    1992, Acta Astronomica, 42, 253

    Mateo, M. 1992, Acta Astronomica, 42, 253

    work page 1992

  26. [26]

    OGLE-IV: Fourth Phase of the Optical Gravitational Lensing Experiment

    Udalski, A., Szyma´ nski, M. K., & Szyma´ nski, G. 2015, Acta Astronomica, 65, 1, doi: 10.48550/arXiv.1504.05966 Wo´ zniak, P. 2000, Acta Astronomica, 50, 421, doi: 10.48550/arXiv.astro-ph/0012143

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1504.05966 2015

  27. [27]

    C., Hwang, K.-H., et al

    Yang, H., Yee, J. C., Hwang, K.-H., et al. 2024, MNRAS, 528, 11, doi: 10.1093/mnras/stad3672

    work page doi:10.1093/mnras/stad3672 2024

  28. [28]

    C., Shvartzvald, Y., Gal-Yam, A., et al

    Yee, J. C., Shvartzvald, Y., Gal-Yam, A., et al. 2012, ApJ, 755, 102, doi: 10.1088/0004-637X/755/2/102

    work page doi:10.1088/0004-637x/755/2/102 2012

  29. [29]

    2021, AJ, 162,

    Yoo, J., DePoy, D. L., A., G.-Y., et al. 2004, ApJ, 603, 13, doi: 10.1086/381241

    work page doi:10.1086/381241 2004