PopSyCLE: A New Population Synthesis Code for Compact Object Microlensing Events

Abstract

We present a new Milky Way microlensing simulation code, dubbed PopSyCLE (Population Synthesis for Compact object Lensing Events). PopSyCLE is the first resolved microlensing simulation to include a compact object distribution derived from numerical supernovae explosion models and both astrometric and photometric microlensing effects. We demonstrate the capabilities of PopSyCLE by investigating the optimal way to find black holes (BHs) with microlensing. Candidate BHs have typically been selected from wide-field photometric microlensing surveys, such as OGLE, by selecting events with long Einstein crossing times ($t_E>120$ days). These events can be selected at closest approach and monitored astrometrically in order to constrain the mass of each lens; PopSyCLE predicts a BH detection rate of ~40% for such a program. We find that the detection rate can be enhanced to ~85% by selecting events with both $t_E>120$ days and a microlensing parallax of ${\pi }_E<0.08$. Unfortunately, such a selection criterion cannot be applied during the event as ${\pi }_E$ requires both pre- and post-peak photometry. However, historical microlensing events from photometric surveys can be revisited using this new selection criteria in order to statistically constrain the abundance of BHs in the Milky Way. The future WFIRST microlensing survey provides both precise photometry and astrometry and will yield individual masses of $\mathcal{O}(100-1000) black holes, which is at least an order of magnitude more than is possible with individual candidate follow-up with current facilities. The resulting sample of BH masses from WFIRST will begin to constrain the shape of the black hole present-day mass function, BH multiplicity, and BH kick velocity distributions.

33 pages, 22 figures, accepted to ApJ

References in corpus (24)

• The NumPy array: a structure for efficient numerical computation
• The Gaia mission
• Modules for Experiments in Stellar Astrophysics (MESA)
• X-ray Properties of Black-Hole Binaries
• New evolutionary models for pre-main sequence and main sequence low-mass stars down to the hydrogen-burning limit
• Primordial Black Holes as Dark Matter
• Limits on the Macho Content of the Galactic Halo from the EROS-2 Survey of the Magellanic Clouds
• Constraining the Maximum Mass of Neutron Stars From Multi-Messenger Observations of GW170817
• Dynamical modelling of the Galactic bulge and bar: the Milky Way's bar pattern speed, stellar, and dark matter mass distribution
• The Structure of the Milky Way's Bar Outside the Bulge
• No large population of unbound or wide-orbit Jupiter-mass planets
• The Initial-Final Mass Relation: Direct Constraints at the Low Mass End
• Investigating stellar-mass black hole kicks
• Stellar Proper Motions in the Galactic Bulge from deep HST ACS/WFC Photometry
• Kinematics at the Edge of the Galactic Bulge: Evidence for Cylindrical Rotation
• Relativistic deflection of background starlight measures the mass of a nearby white dwarf star
• The Bulge Radial Velocity Assay (BRAVA): I. Sample Selection and a Rotation Curve
• The Galactic distribution of X-ray binaries and its implications for compact object formation and natal kicks
• OGLE-III Microlensing Events and the Structure of the Galactic Bulge
• Hunting Black Holes with Gaia
• A Search for Stellar-Mass Black Holes via Astrometric Microlensing
• Extinction law in the range 0.4-4.8 $\mu$m and the 8620 $\mathring{A}$ DIB towards the stellar cluster Westerlund 1
• Parameter degeneracies and (un)predictability of gravitational microlensing events
• Gravitational Microlensing by Neutron Stars and Radio Pulsars: Event Rates, Timescale Distributions, and Mass Measurements