A differentiable binary microlensing model using adaptive contour integration method
arXiv:2501.07268 · doi:10.3847/1538-3881/adb1b2
Abstract
We present microlux, which is a Jax-based code that can compute the binary microlensing light curve and its derivatives both efficiently and accurately. The key feature of microlux is the implementation of a modified version of the adaptive sampling algorithm that was originally proposed by V. Bozza to account for the finite-source effect most efficiently. The efficiency and accuracy of microlux have been verified across the relevant parameter space for binary microlensing. As a differentiable code, microlux makes it possible to apply gradient-based algorithms to the search and posterior estimation of the microlensing modeling. As an example, we use microlux to model a real microlensing event and infer the model posterior via both Fisher information matrix and Hamiltonian Monte Carlo, neither of which would have been possible without the access to accurate model gradients.
arXiv comment14 pages, 7 figures. Accepted by AJ. GitHub repo at https://github.com/CoastEgo/microlux
References in corpus (30)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Array Programming with NumPy
- emcee: The MCMC Hammer
- The NumPy array: a structure for efficient numerical computation
- The NASA Exoplanet Archive: Data and Tools for Exoplanet Research
- Real-Time Difference Imaging Analysis of MOA Galactic Bulge Observations During 2000
- Discovery of a Jupiter/Saturn Analog with Gravitational Microlensing
- A Review of automatic differentiation and its efficient implementation
- Detecting Earth-Mass Planets with Gravitational Microlensing
- Predictions of the WFIRST Microlensing Survey I: Bound Planet Detection Rates
- Planetary Detection Efficiency of the Magnification 3000 Microlensing Event OGLE-2004-BLG-343
- An isolated mass gap black hole or neutron star detected with astrometric microlensing
- Extended-Source Effect and Chromaticity in Two-Point-Mass Microlensing
- Microlensing with advanced contour integration algorithm: Green's theorem to third order, error control, optimal sampling and limb darkening
- Hexadecapole Approximation in Planetary Microlensing
- A Terrestrial Planet in a ~1 AU Orbit Around One Member of a ~15 AU Binary
- VBBinaryLensing: a public package for microlensing light curve computation
- Systematic KMTNet Planetary Anomaly Search, Paper I: OGLE-2019-BLG-1053Lb, A Buried Terrestrial Planet
- Modeling microlensing events with MulensModel
- Systematic errors as a source of mass discrepancy in black hole microlensing event OGLE-2011-BLG-0462
- Spitzer as Microlens Parallax Satellite: Mass and Distance Measurements of Binary Lens System OGLE-2014-BLG-1050L
- Adaptive Contouring -- an efficient way to calculate microlensing light curves of extended sources
- Real-Time Likelihood-Free Inference of Roman Binary Microlensing Events with Amortized Neural Posterior Estimation
- Measuring microlensing parallax via simultaneous observations from Chinese Space Station Telescope and Roman Telescope
- KMT-2019-BLG-0371 and the Limits of Bayesian Analysis
- Twinkle: A GPU-based binary-lens microlensing code with contour integration method
- Improved Aberth-Ehrlich root-finding algorithm and its further application for Binary Microlensing
- Fast computation of quadrupole and hexadecapole approximations in microlensing with a single point-source evaluation
- On the Perturbative Picture and the Chang-Refsdal Lens Approximation for Planetary Microlensing
- MAGIC: Microlensing Analysis Guided by Intelligent Computation