Assessing the Impact of Binary Systems on Microlensing Using SPISEA and PopSyCLE Population Simulations

TL;DR

This study enhances microlensing population simulations by incorporating binary systems using SPISEA and PopSyCLE, revealing that binaries significantly influence observed event characteristics and should be more thoroughly modeled.

Contribution

The paper introduces a method to include binary systems in microlensing simulations, improving the match between simulated and observed event distributions.

Findings

55% of microlensing events involve binaries.

Including binaries shifts the Einstein crossing time distribution closer to observations.

Most binary events appear as single-lens, single-source lightcurves.

Abstract

Gravitational microlensing provides a unique opportunity to probe the mass distribution of stars, black holes, and other objects in the Milky Way. Population simulations are necessary to interpret results from microlensing surveys. The contribution from binary objects is often neglected or minimized in analysis of observations and simulations despite the high percentage of binary systems and microlensing's ability to probe binaries. To simulate the population effects we added multiple systems to Stellar Population Interface for Stellar Evolution and Atmospheres (SPISEA), which simulates stellar clusters. We then inject these multiples into Population Synthesis for Compact-object Lensing Events (PopSyCLE), which simulates Milky Way microlensing surveys. When making OGLE observational selection criteria, we find that 55% of observed microlensing events involve a binary system. Specifically, 14.5% of events have a multiple-lens and a single source, 31.7% have a single lens and a multiple-source, and 8.8% have a multiple-lens and a multiple-source. The majority of these events have photometric lightcurves that appear single and are fit well by a single-lens, single-source model. This suggests that binary source and binary lens-binary source models should be included more frequently in event analysis. The mean Einstein crossing time shifts from 19.1 days for single events only to 21.3 days for singles and multiple events, after cutting binary events with multiple peaks. The Einstein crossing time distribution of singles and single-peaked multiple events is better aligned with observed distributions from OGLE (arXiv:1707.07634) than singles alone, indicating that multiple systems are a significant missing piece between simulations and reality.