Estimating microlensing parameters from observables and stellar isochrones with pyLIMASS

E. Bachelet,M. Hundertmark,S. Calchi Novati

2024-05-04

Abstract:We present pyLIMASS, a novel algorithm for estimating the physical properties of the lensing system in microlensing events. The main idea of pyLIMASS is to combine all available information regarding the microlensing event, defined as observables, and to estimate the parameter distributions of the system, such as the lens mass and distance. The algorithm is based on isochrones for the stars model and combine the observables using a Gaussian Mixtures approach. After describing the mathematical formalism and its implementation, we discuss the algorithm's performance on simulated and published events. Generally, the pyLIMASS estimations are in good agreement (i.e., within 1-$\sigma$) with the results of the selected published events, making it an effective tool to estimate the lens properties and their distribution. The applicability of the method was tested by using a catalog of realistically simulated events that could be observed by the future Galactic Bulge Time Domain Survey of the Nancy Grace Roman Space Telescope. By solely using constraints from the Roman lightcurves and images, pyLIMASS estimates the masses of the lens of the Roman catalog with a median precision of 20% with almost no bias.

Instrumentation and Methods for Astrophysics,Earth and Planetary Astrophysics,Astrophysics of Galaxies,Solar and Stellar Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **How to estimate the physical parameters of gravitational lensing events from observational data, especially the mass and distance of the lensing celestial body**. Gravitational microlensing is a powerful astronomical phenomenon that can be used to detect the mass distribution of various celestial bodies in the Milky Way. However, it is usually very difficult to directly measure the physical properties (such as mass, distance) of these celestial bodies. Existing methods rely on light - curve modeling and other auxiliary observations, but these methods often need to combine multiple observational data and have large uncertainties.

To solve this problem, the author proposes a new algorithm **pyLIMASS**, which aims to estimate the parameter distribution of the lensing system, such as lensing mass, distance, etc., by combining all available observational information (defined as "observables") and using the Gaussian Mixture Model (GM). Specifically, the main contributions of pyLIMASS include:

1. **Combining multiple observational information**: The algorithm can integrate multiple observational data from light - curves, parallax, the brightness of the source star and the lensing star, etc.
2. **Using the Gaussian Mixture Model**: Fit the distribution of observational data through the Gaussian Mixture Model, so as to estimate the physical parameters of the lensing system more accurately.
3. **Introducing stellar isochrones as prior information**: Use the stellar isochrone model as prior information to help constrain the physical parameters of the lensing star.
4. **Efficiently handling complex distributions**: The Gaussian Mixture Model can effectively handle multimodal distributions, which are common in gravitational microlensing events because there may be degeneracy or competing models between parameters.

Through these methods, pyLIMASS can estimate the mass of the lensing celestial body with high precision only using the light - curve and image data of the Roman Space Telescope, with a median precision of 20% and almost no bias. This makes pyLIMASS an effective tool for estimating the characteristics and distribution of lensing celestial bodies in gravitational microlensing events.

### Formula Summary

- **Einstein ring crossing time \(t_E\)**:
\[
t_E=\frac{\theta_E}{\mu_{\text{rel}}}
\]
where \(\mu_{\text{rel}}=\mu_L - \mu_S\) is the relative proper motion vector between the lens and the source star.

- **Relative parallax \(\pi_{\text{rel}}\)**:
\[
\pi_{\text{rel}}=\pi_L-\pi_S
\]

- **Einstein angular radius \(\theta_E\)**:
\[
\theta_E = \sqrt{\kappa M_L\pi_{\text{rel}}}
\]
where \(\kappa\) is a constant and \(M_L\) is the lensing mass (in solar masses).

- **Probability density function of the Gaussian Mixture Model**:
\[
p(y|\theta)=\sum_{k = 1}^G\phi_kN(\mu_k,\Sigma_k,y)
\]
where \(\theta=(\phi_k,\mu_k,\Sigma_k)\) is a set of hyperparameters, \(\mu_k\) and \(\Sigma_k\) are the mean vector and covariance matrix of the \(k\)-th multivariate Gaussian component respectively, and \(\phi_k\) is the weight of the \(k\)-th Gaussian component.

Through these formulas and methods, pyLIMASS can effectively estimate the characteristics of the lensing celestial body in gravitational microlensing events, thus providing a new tool for understanding the celestial body distribution in the Milky Way.