skip to main content


Title: Gravitational Microlensing Rates in Milky Way Globular Clusters
Abstract

Many recent observational and theoretical studies suggest that globular clusters (GCs) host compact object populations large enough to play dominant roles in their overall dynamical evolution. Yet direct detection, particularly of black holes and neutron stars, remains rare and limited to special cases, such as when these objects reside in close binaries with bright companions. Here we examine the potential of microlensing detections to further constrain these dark populations. Based on state-of-the-art GC models from theCMC Cluster Catalog, we estimate the microlensing event rates for black holes, neutron stars, white dwarfs (WDs), and, for comparison, also for M dwarfs in Milky Way GCs, as well as the effects of different initial conditions on these rates. Among compact objects, we find that WDs dominate the microlensing rates, simply because they largely dominate by numbers. We show that microlensing detections are in general more likely in GCs with higher initial densities, especially in clusters that undergo core collapse. We also estimate microlensing rates in the specific cases of M22 and 47 Tuc using our best-fitting models for these GCs. Because their positions on the sky lie near the rich stellar backgrounds of the Galactic bulge and the Small Magellanic Cloud, respectively, these clusters are among the Galactic GCs best suited for dedicated microlensing surveys. The upcoming 10 yr survey with the Rubin Observatory may be ideal for detecting lensing events in GCs.

 
more » « less
Award ID(s):
2001751 2108624
NSF-PAR ID:
10365711
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
928
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 181
Size(s):
["Article No. 181"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The globular cluster 47 Tucanae (47 Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars (N≳ 106) and high density. Here we present detailed and self-consistent 47 Tuc models computed with theCluster Monte Carlocode (CMC). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47 Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions (as in, e.g., Kroupa 2001), and an initial Elson profile (Elson et al. 1987) that is overfilling the cluster’s tidal radius. We include new prescriptions inCMCfor the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47 Tuc; our best-fit model contains ∼50 millisecond pulsars, 70% of which are formed through giant collisions and tidal captures. Our models also suggest that 47 Tuc presently contains up to ∼200 stellar-mass black holes, ∼5 binary black holes, ∼15 low-mass X-ray binaries, and ∼300 cataclysmic variables.

     
    more » « less
  2. Abstract

    The dense environments in the cores of globular clusters (GCs) facilitate many strong dynamical encounters among stellar objects. These encounters have been shown to be capable of ejecting stars from the host GC, whereupon they become runaway stars, or hypervelocity stars (HVSs) if unbound to the galactic potential. We study high-speed stellar ejecta originating from GCs by using Monte CarloN-body models, in particular focusing on binary–single encounters involving compact objects. We pair our model-discriminated populations with observational catalogs of Milky Way GCs (MWGCs) to compose a present-day Galactic population of stellar ejecta. We find that these kinds of encounters can accelerate stars to velocities in excess of 2000 km s−1, to speeds beyond the previously predicted limits for ejecta from star-only encounters and in the same regime of Galactic center ejections. However, the same ejections can only account for 1.5%–20% of the total population of stellar runaways, and only 0.0001%–1% of HVS, with similar relative rates found for runaway white dwarfs. We also provide credible regions for ejecta from 149 MWGCs, which we hope will be useful as supplementary evidence when pairing runaway stars with origin GCs.

     
    more » « less
  3. Abstract

    Uncertainty in the initial–final mass relation (IFMR) has long been a problem in understanding the final stages of massive star evolution. One of the major challenges of constraining the IFMR is the difficulty of measuring the mass of nonluminous remnant objects (i.e., neutron stars and black holes). Gravitational-wave detectors have opened the possibility of finding large numbers of compact objects in other galaxies, but all in merging binary systems. Gravitational lensing experiments using astrometry and photometry are capable of finding compact objects, both isolated and in binaries, in the Milky Way. In this work we improve the Population Synthesis for Compact object Lensing Events (PopSyCLE)microlensing simulation code in order to explore the possibility of constraining the IFMR using the Milky Way microlensing population. We predict that the Roman Space Telescope’s microlensing survey will likely be able to distinguish different IFMRs based on the differences at the long end of the Einstein crossing time distribution and the small end of the microlensing parallax distribution, assuming the small (πE≲ 0.02) microlensing parallaxes characteristic of black hole lenses are able to be measured accurately. We emphasize that future microlensing surveys need to be capable of characterizing events with small microlensing parallaxes in order to place the most meaningful constraints on the IFMR.

     
    more » « less
  4. Abstract

    Gravitational-wave detections are enabling measurements of the rate of coalescences of binaries composed of two compact objects—neutron stars and/or black holes. The coalescence rate of binaries containing neutron stars is further constrained by electromagnetic observations, including Galactic radio binary pulsars and short gamma-ray bursts. Meanwhile, increasingly sophisticated models of compact objects merging through a variety of evolutionary channels produce a range of theoretically predicted rates. Rapid improvements in instrument sensitivity, along with plans for new and improved surveys, make this an opportune time to summarise the existing observational and theoretical knowledge of compact-binary coalescence rates.

     
    more » « less
  5. ABSTRACT

    Stellar-mass black holes (BHs) can be retained in globular clusters (GCs) until the present. Simulations of GC evolution find that the relaxation driven mass-loss rate is elevated if BHs are present, especially near dissolution. We capture this behaviour in a parametrized mass-loss rate, bench marked by results from N-body simulations, and use it to evolve an initial GC mass function (GCMF), similar to that of young massive clusters in the Local Universe, to an age of 12 Gyr. Low-metallicity GCs ([Fe/H] ≲ −1.5) have the highest mass-loss rates, because of their relatively high BH masses, which combined with their more radial orbits and stronger tidal field in the past explains the high turnover mass of the GCMF ($\sim 10^5\, {\rm M}_\odot$ ) at large Galactic radii ($\gtrsim 10\, {\rm kpc}$ ). The turnover mass at smaller Galactic radii is similar because of the upper mass truncation of the initial GCMF and the lower mass-loss rate due to the higher metallicities. The density profile in the Galaxy of mass lost from massive GCs ($\gtrsim 10^{5}\, {\rm M}_\odot$ ) resembles that of nitrogen-rich stars in the halo, confirming that these stars originated from GCs. We conclude that two-body relaxation is the dominant effect in shaping the GCMF from a universal initial GCMF, because including the effect of BHs reduces the need for additional disruption mechanisms.

     
    more » « less