Search

We present a new method to combine multimass equilibrium dynamical models and pulsar timing data to constrain the mass distribution and remnant populations of Milky Way globular clusters (GCs). We first apply this method to 47 Tuc, a cluster for which there exists an abundance of stellar kinematic data and which is also host to a large population of millisecond pulsars. We demonstrate that the pulsar timing data allow us to place strong constraints on the overall mass distribution and remnant populations even without fitting on stellar kinematics. Our models favor a small population of stellar-mass BHs in this cluster (with a total mass of $446^{+75}_{-72} \mathrm{M_\odot}$), arguing against the need for a large ($ > 2000 \ \mathrm{M_\odot}$) central intermediate-mass black hole. We then apply the method to Terzan 5, a heavily obscured bulge cluster which hosts the largest population of millisecond pulsars of any Milky Way GC and for which the collection of conventional stellar kinematic data is very limited. We improve existing constraints on the mass distribution and structural parameters of this cluster and place stringent constraints on its black hole content, finding an upper limit on the mass in BHs of $\sim 4000 \ \mathrm{M_\odot}$. This method allows us to probe the central dynamics of GCs even in the absence of stellar kinematic data and can be easily applied to other GCs with pulsar timing data, for which datasets will continue to grow with the next generation of radio telescopes., Comment: 33 pages, 19 figures. Accepted for publication in ApJ

The Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR) is a proposed Canadian-led 1m-class space telescope that will carry out ultraviolet and blue-optical wide-field imaging, spectroscopy, and photometry. CASTOR will provide an essential bridge in the post-Hubble era, preventing a protracted UV-optical gap in space astronomy and enabling an enormous range of discovery opportunities from the solar system to the nature of the Cosmos, in conjunction with the other great wide-field observatories of the next decade (e.g., Euclid, Roman, Vera Rubin). FORECASTOR (Finding Optics Requirements and Exposure times for CASTOR) will supply a coordinated suite of mission-planning tools that will serve as the one-stop shop for proposal preparation, data reduction, and analysis for the CASTOR mission. We present the first of these tools: a pixel-based, user-friendly, extensible, multi-mission exposure time calculator (ETC) built in Python, including a modern browser-based graphical user interface that updates in real time. We then provide several illustrative examples of FORECASTOR's use that advance the design of planned legacy surveys for the CASTOR mission: a search for the most massive white dwarfs in the Magellanic Clouds; a study of the frequency of flaring activity in M stars, their distribution and impacts on habitability of exoplanets; mapping the proper motions of faint stars in the Milky Way; wide and deep galaxy surveys; and time-domain studies of active galactic nuclei., Comment: Updated references and acknowledgements to match published version. 24 pages, 16 figures, 3 tables, published in AJ

Populations of stellar-mass black holes (BHs) in globular clusters (GCs) influence their dynamical evolution and have important implications on one of the main formation channels for gravitational wave sources. Inferring the size of these populations remains difficult, however. In this work, multimass models of 34 Milky Way GCs, first presented in Dickson et al., are used to explore the present-day BH populations. Direct constraints on both the total and visible mass components provided by several observables allow these models to accurately determine the distribution of the dark mass (including BHs) within clusters, as we demonstrate in a proof-of-concept fitting of the models to mock observations extracted from Monte Carlo cluster models. New constraints on the BH population retained to the present-day in each cluster are inferred from our models. We find that BH mass fractions ranging from 0 to 1 per cent of the total mass are typically required to explain the observations, except for Omega Cen, for which we infer a mass fraction above 5 per cent, in agreement with previous works. Relationships between the dark remnant populations and other cluster parameters are examined, demonstrating a clear anti-correlation between the amount of BHs and mass segregation between visible stars, as well as a correlation between remnant mass fractions and the dynamical age of clusters. Our inferred BH populations are in good agreement overall with other recent studies using different methodologies, but with notable discrepancies for individual clusters., Comment: 18 pages, 9 figures, MNRAS accepted

The distribution of stars and stellar remnants (white dwarfs, neutron stars, black holes) within globular clusters holds clues about their formation and long-term evolution, with important implications for their initial mass function (IMF) and the formation of black hole mergers. In this work, we present best-fitting multimass models for 37 Milky Way globular clusters, which were inferred from various datasets, including proper motions from Gaia EDR3 and HST, line-of-sight velocities from ground-based spectroscopy and deep stellar mass functions from HST. We use metallicity dependent stellar evolution recipes to obtain present-day mass functions of stars and remnants from the IMF. By dynamically probing the present-day mass function of all objects in a cluster, including the mass distribution of remnants, these models allow us to explore in detail the stellar (initial) mass functions of a large sample of Milky Way GCs. We show that, while the low-mass mass function slopes are strongly dependent on the dynamical age of the clusters, the high-mass slope ($\alpha_3; m > 1 M_\odot$) is not, indicating that the mass function in this regime has generally been less affected by dynamical mass loss. Examination of this high-mass mass function slope suggests an IMF in this mass regime consistent with a Salpeter IMF is required to reproduce the observations. This high-mass IMF is incompatible with a top-heavy IMF, as has been proposed recently. Finally, based on multimass model fits to our sample of Milky Way GCs, no significant correlation is found between the high-mass IMF slope and cluster metallicity., Comment: 21 pages, 10 figures, MNRAS accepted

We have determined stellar mass functions of 120 Milky Way globular clusters and massive LMC/SMC star clusters based on a comparison of archival Hubble Space Telescope photometry with a large grid of direct N-body simulations. We find a strong correlation of the global mass function slopes of star clusters with both their internal relaxation times as well as their lifetimes. Once dynamical effects are being accounted for, the mass functions of most star clusters are compatible with an initial mass function described by a broken power-law distribution $N(m) \sim m^\alpha$ with break masses at 0.4 M$_\odot$ and 1.0 M$_\odot$ and mass function slopes of $\alpha_{Low}=-0.3$ for stars with masses $m<0.4$ M$_\odot$, $\alpha_{High}=-2.30$ for stars with $m>1.0$ M$_\odot$ and $\alpha_{Med}=-1.65$ for intermediate-mass stars. Alternatively, a log-normal mass function with a characteristic mass $\log M_C = -0.36$ and width $\sigma_C=0.28$ for low-mass stars and a power-law mass function for stars with $m>1$ M$_\odot$ also fits our data. We do not find a significant environmental dependency of the initial mass function with either cluster mass, density, global velocity dispersion or metallicity. Our results lead to a larger fraction of high-mass stars in globular clusters compared to canonical Kroupa/Chabrier mass functions, increasing the efficiency of self-enrichment in clusters and helping to alleviate the mass budget problem of multiple stellar populations in globular clusters. By comparing our results with direct N-body simulations we finally find that only simulations in which most black holes are ejected by natal birth kicks correctly reproduce the observed correlations., Comment: 20 pages, 9 figures, accepted for publication in MNRAS, data is available from http://people.smp.uq.edu.au/HolgerBaumgardt/globular/index.html