Showing 1–50 of 167 results for author: Gieles, M

Probing populations of dark stellar remnants in the globular clusters 47 Tuc and Terzan 5 using pulsar timing

Authors: Peter J. Smith, Vincent Hénault-Brunet, Nolan Dickson, Mark Gieles, Holger Baumgardt

Abstract: 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 pul… ▽ More 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 $412^{+77}_{-146} \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 4500 \ \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. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 32 pages, 19 figures. Submitted to ApJ, comments welcome

Binary-single interactions with different mass ratios

Authors: Bruno Rando Forastier, Daniel Marín Pina, Mark Gieles, Simon Portegies Zwart, Fabio Antonini

Abstract: Dynamical interactions in star clusters are an efficient mechanism to produce the coalescing binary black holes (BBHs) that have been detected with gravitational waves (GWs). We want to understand how BBH coalescence can occur during - or after - binary-single interactions with different mass ratios. We perform gravitational scattering experiments of binary-single interactions using different mass… ▽ More Dynamical interactions in star clusters are an efficient mechanism to produce the coalescing binary black holes (BBHs) that have been detected with gravitational waves (GWs). We want to understand how BBH coalescence can occur during - or after - binary-single interactions with different mass ratios. We perform gravitational scattering experiments of binary-single interactions using different mass ratios of the binary components ($q_2\equiv m_2/m_1\le1$) and the incoming single ($q_3\equiv m_3/m_1$). We extract cross sections and rates for (i) GW capture during resonant interactions; (ii) GW inspiral in between resonant interactions and apply the results to different globular cluster conditions. We find that GW capture during resonant interactions is most efficient if $q_2\simeq q_3$ and that the mass-ratio distribution of BBH coalescence due to inspirals is $\propto m_1^{-1}q^{2.9+α}$, where $α$ is the exponent of the BH mass function. The total rate of GW captures and inspirals depends mostly on $m_1$ and is relatively insensitive to $q_2$ and $q_3$. We show that eccentricity increase by non-resonant encounters approximately doubles the rate of BBH inspiral in between resonant encounters. For a given GC mass and radius, the BBH merger rate in metal-rich GCs is approximately double that of metal-poor GCs, because of their (on average) lower BH masses ($m_1$) and steeper BH mass function, yielding binaries with lower $q$. Our results enable the translating from the mass-ratio distribution of dynamically formed BBH mergers to the underlying BH mass function. The additional mechanism that leads to a doubling of the inspirals provides an explanation for the reported high fraction of in-cluster inspirals in $N$-body models of clusters. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 14 pages, 7 figures

Capture of field stars by dark substructures

Authors: Jorge Peñarrubia, Raphaël Errani, Matthew G. Walker, Mark Gieles, Tjarda C. N. Boekholt

Abstract: We use analytical and $N$-body methods to study the capture of field stars by gravitating substructures moving across a galactic environment. We find that the majority of captured stars move on temporarily-bound orbits that escape from the substructure potential after a few orbital revolutions. In numerical experiments where a substructure model is immersed in a galaxy on a circular orbit, we also… ▽ More We use analytical and $N$-body methods to study the capture of field stars by gravitating substructures moving across a galactic environment. We find that the majority of captured stars move on temporarily-bound orbits that escape from the substructure potential after a few orbital revolutions. In numerical experiments where a substructure model is immersed in a galaxy on a circular orbit, we also find particles that remain bound to the substructure potential for indefinitely-long times. This population is absent from substructure models initially placed outside the galaxy on an eccentric orbit. We show that gravitational capture is most efficient in dwarf spheroidal galaxies (dSphs) on account of their low velocity dispersions and high stellar phase-space densities. In these galaxies `dark' sub-subhaloes which do not experience in-situ star formation may capture field stars and become visible as stellar overdensities with unusual properties: (i) they would have a large size for their luminosity, (ii) contain stellar populations indistinguishable from the host galaxy, and (iii) exhibit dark matter (DM)-dominated mass-to-light ratios. We discuss the nature of several `anomalous' stellar systems reported as star clusters in the Fornax and Eridanus II dSphs which exhibit some of these characteristics. A large population of DM sub-subhaloes with a mass function $d N/d M_\bullet\sim M_\bullet^{-α}$ generates stellar substructures with a luminosity function, $d N/d M_\star\sim M_\star^{-β}$, where $β=(2α+1)/3=1.6$ for $α=1.9$. Detecting and characterizing these objects in dSphs would provide unprecedented constraints on the particle mass and cross section of a large range of DM particle candidates. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 21 pages. Submitted

Dynamical formation of Gaia BH3 in the progenitor globular cluster of the ED-2 stream

Authors: Daniel Marín Pina, Sara Rastello, Mark Gieles, Kyle Kremer, Laura Fitzgerald, Bruno Rando Forastier

Abstract: Context. The star-black hole (S-BH) binary known as Gaia BH3, discovered by the Gaia Collaboration is chemically and kinematically associated with the metal-poor ED-2 stream in the Milky Way halo. Aims. We explore the possibility that Gaia BH3 was assembled dynamically in the progenitor globular cluster (GC) of the ED-2 stream. Methods. We used a public suite of star-by-star dynamical Monte Ca… ▽ More Context. The star-black hole (S-BH) binary known as Gaia BH3, discovered by the Gaia Collaboration is chemically and kinematically associated with the metal-poor ED-2 stream in the Milky Way halo. Aims. We explore the possibility that Gaia BH3 was assembled dynamically in the progenitor globular cluster (GC) of the ED-2 stream. Methods. We used a public suite of star-by-star dynamical Monte Carlo models to identify S-BH binaries in GCs with different initial masses and (half-mass) radii. Results. We show that a likely progenitor of the ED-2 stream was a relatively low-mass ($\lesssim10^5M_\odot$) GC with an initial half-mass radius of ~4 pc. Such a GC can dynamically retain a large fraction of its BH population and dissolve on the orbit of ED-2. From the suite of models we find that GCs produce ~ 3 - 30 S-BH binaries, approximately independent of initial GC mass and inversely correlated with initial cluster radius. Scaling the results to the Milky Way GC population, we find that ~75% of the S-BH binaries formed in GCs are ejected from their host GC, all in the early phases of evolution ($\lesssim1$ Gyr); these are expected to no longer be close to streams. The ~25% of S-BH binaries retained until dissolution are expected to form part of streams, such that for an initial mass of the progenitor of ED-2 of a few $10^4M_\odot$, we expect ~2-3 S-BH to end up in the stream. GC models with metallicities similar to Gaia BH3 ($\lesssim1\%$ solar) include S-BH binaries with similar BH masses ($\gtrsim30M_\odot$), orbital periods, and eccentricities. Conclusions. We predict the Galactic halo contains of order $10^5$ S-BH binaries that formed dynamically in GCs, a fraction of which may readily be detected in Gaia DR4. The detection of these sources provides valuable tests of BH dynamics in clusters and the possible role in formation of gravitational wave sources. △ Less

Submitted 9 July, 2024; v1 submitted 19 April, 2024; originally announced April 2024.

Comments: 5 pages, 3 figures, accepted in A&A

Stellar streams from black hole-rich star clusters

Authors: Daniel Roberts, Mark Gieles, Denis Erkal, Jason L. Sanders

Abstract: Nearly a hundred progenitor-less, thin stellar streams have been discovered in the Milky Way, thanks to Gaia and related surveys. Most streams are believed to have formed from star clusters and it was recently proposed that extended star clusters -- rich in stellar-mass black holes (BHs) -- are efficient in creating streams. To understand the nature of stream progenitors better, we quantify the di… ▽ More Nearly a hundred progenitor-less, thin stellar streams have been discovered in the Milky Way, thanks to Gaia and related surveys. Most streams are believed to have formed from star clusters and it was recently proposed that extended star clusters -- rich in stellar-mass black holes (BHs) -- are efficient in creating streams. To understand the nature of stream progenitors better, we quantify the differences between streams originating from star clusters with and without BHs using direct $N$-body models and a new model for the density profiles of streams based on time-dependent escape rates from clusters. The QSG (Quantifying Stream Growth) model facilitates the rapid exploration of parameter space and provides an analytic framework to understand the impact of different star cluster properties and escape conditions on the structure of streams. Using these models it is found that, compared to streams from BH-free clusters on the same orbit, streams of BH-rich clusters: (1) are approximately five times more massive; (2) have a peak density three times closer to the cluster 1 Gyr post-dissolution (for orbits of Galactocentric radius > 10 kpc), and (3) have narrower peaks and more extended wings in their density profile. We discuss other observable stream properties that are affected by the presence of BHs in their progenitor cluster, namely the width of the stream, its radial offset from the orbit, and the properties of the gap at the progenitor's location. Our results provide a step towards using stellar streams to constrain the BH content of dissolved (globular) star clusters. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 13 pages, 8 figures

HRMOS White Paper: Science Motivation

Authors: Laura Magrini, Thomas Bensby, Anna Brucalassi, Sofia Randich, Robin Jeffries, Gayandhi de Silva, Asa Skuladottir, Rodolfo Smiljanic, Oscar Gonzalez, Vanessa Hill, Nadege Lagarde, Eline Tolstoy, Jose' Maria Arroyo-Polonio, Martina Baratella, John R. Barnes, Giuseppina Battaglia, Holger Baumgardt, Michele Bellazzini, Katia Biazzo, Angela Bragaglia, Bradley Carter, Giada Casali, Gabriele Cescutti, Camilla Danielski, Elisa Delgado Mena , et al. (30 additional authors not shown)

Abstract: The High-Resolution Multi-Object Spectrograph (HRMOS) is a facility instrument that we plan to propose for the Very Large Telescope (VLT) of the European Southern Observatory (ESO), following the initial presentation at the VLT 2030 workshop held at ESO in June 2019. HRMOS provides a combination of capabilities that are essential to carry out breakthrough science across a broad range of active res… ▽ More The High-Resolution Multi-Object Spectrograph (HRMOS) is a facility instrument that we plan to propose for the Very Large Telescope (VLT) of the European Southern Observatory (ESO), following the initial presentation at the VLT 2030 workshop held at ESO in June 2019. HRMOS provides a combination of capabilities that are essential to carry out breakthrough science across a broad range of active research areas from stellar astrophysics and exoplanet studies to Galactic and Local Group archaeology. HRMOS fills a gap in capabilities amongst the landscape of future instrumentation planned for the next decade. The key characteristics of HRMOS will be high spectral resolution (R = 60000 - 80000) combined with multi-object (20-100) capabilities and long term stability that will provide excellent radial velocity precision and accuracy (10m/s). Initial designs predict that a SNR~100 will be achievable in about one hour for a star with mag(AB) = 15, while with the same exposure time a SNR~ 30 will be reached for a star with mag(AB) = 17. The combination of high resolution and multiplexing with wavelength coverage extending to relatively blue wavelengths (down to 380\,nm), makes HRMOS a spectrograph that will push the boundaries of our knowledge and that is envisioned as a workhorse instrument in the future. The science cases presented in this White Paper include topics and ideas developed by the Core Science Team with the contributions from the astronomical community, also through the wide participation in the first HRMOS Workshop (https://indico.ict.inaf.it/event/1547/) that took place in Firenze (Italy) in October 2021. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 88 pages, 39 figures. Comments and expressions of interest are welcome by contacting members of the Core Science Team

The influence of black holes on the binary population of the globular cluster Palomar 5

Authors: Long Wang, Mark Gieles, Holger Baumgardt, Chengyuan Li, Xiaoying Pang, Baitian Tang

Abstract: The discovery of stellar-mass black holes (BHs) in globular clusters (GCs) raises the possibility of long-term retention of BHs within GCs. These BHs influence various astrophysical processes, including merger-driven gravitational waves and the formation of X-ray binaries. They also impact cluster dynamics by heating and creating low-density cores. Previous N-body models suggested that Palomar 5,… ▽ More The discovery of stellar-mass black holes (BHs) in globular clusters (GCs) raises the possibility of long-term retention of BHs within GCs. These BHs influence various astrophysical processes, including merger-driven gravitational waves and the formation of X-ray binaries. They also impact cluster dynamics by heating and creating low-density cores. Previous N-body models suggested that Palomar 5, a low-density GC with long tidal tails, may contain more than 100 BHs. To test this scenario, we conduct N-body simulations of Palomar 5 with primordial binaries to explore the influence of BHs on binary populations and the stellar mass function. Our results show that primordial binaries have minimal effect on the long-term evolution. In dense clusters with BHs, the fraction of wide binaries with periods >$10^5$ days decreases, and the disruption rate is independent of the initial period distribution. Multi-epoch spectroscopic observations of line-of-sight velocity changes can detect most bright binaries with periods below $10^4$ days, significantly improving velocity dispersion measurements. Four BH-MS binaries in the model with BHs suggests their possible detection through the same observation method. Including primordial binaries leads to a flatter inferred mass function because of spatially unresolved binaries, leading to a better match of the observations than models without binaries, particularly in Palomar 5's inner region. Future observations should focus on the cluster velocity dispersion and binaries with periods of $10^4-10^5$ days in Palomar 5's inner and tail regions to constrain BH existence. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: 20 pages, 19 figures, submitted to MNRAS

A closer look at the binary content of NGC 1850

Authors: S. Saracino, S. Kamann, N. Bastian, M. Gieles, T. Shenar, N. Reindl, J. Müller-Horn, C. Usher, S. Dreizler, V. Hénault-Brunet

Abstract: Studies of young clusters have shown that a large fraction of O-/early B-type stars are in binary systems, where the binary fraction increases with mass. These massive stars are present in clusters of a few Myrs, but gradually disappear for older clusters. The lack of detailed studies of intermediate-age clusters has meant that almost no information is available on the multiplicity properties of s… ▽ More Studies of young clusters have shown that a large fraction of O-/early B-type stars are in binary systems, where the binary fraction increases with mass. These massive stars are present in clusters of a few Myrs, but gradually disappear for older clusters. The lack of detailed studies of intermediate-age clusters has meant that almost no information is available on the multiplicity properties of stars with M $<$ 4 $M_{\odot}$. In this study we present the first characterization of the binary content of NGC 1850, a 100 Myr-old massive star cluster in the Large Magellanic Cloud, relying on a VLT/MUSE multi-epoch spectroscopic campaign. By sampling stars down to M = 2.5 $M_{\odot}$, we derive a close binary fraction of 24 $\pm$ 5 \% in NGC 1850, in good agreement with the multiplicity frequency predicted for stars of this mass range. We also find a trend with stellar mass (magnitude), with higher mass (brighter) stars having higher binary fractions. We modeled the radial velocity curves of individual binaries using The Joker and constrained the orbital properties of 27 systems, $\sim$17\% of all binaries with reliable radial velocities in NGC 1850. This study has brought to light a number of interesting objects, such as four binaries showing mass functions f(M) $>$ 1.25 $M_{\odot}$. One of these, star #47, has a peculiar spectrum, explainable with the presence of two disks in the system, around the visible star and the dark companion, which is a black hole candidate. These results confirm the importance and urgency of studying the binary content of clusters of any age. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: 24 pages, 22 Figures, 1 Table; Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

Multimass modelling of Milky Way globular clusters -- II. present-day black hole populations

Authors: Nolan Dickson, Peter J. Smith, Vincent Hénault-Brunet, Mark Gieles, Holger Baumgardt

Abstract: 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-da… ▽ More 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. △ Less

Submitted 12 February, 2024; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 18 pages, 9 figures, MNRAS accepted

Demographics of three-body binary black holes in star clusters: implications for gravitational waves

Authors: Daniel Marín Pina, Mark Gieles

Abstract: To explain both the dynamics of a globular cluster and its production of gravitational waves from coalescing binary black holes, it is necessary to understand its population of dynamically-formed (or, `three-body') binaries. We provide a theoretical understanding of this population, benchmarked by direct $N$-body models. We find that $N$-body models of clusters on average have only one three-body… ▽ More To explain both the dynamics of a globular cluster and its production of gravitational waves from coalescing binary black holes, it is necessary to understand its population of dynamically-formed (or, `three-body') binaries. We provide a theoretical understanding of this population, benchmarked by direct $N$-body models. We find that $N$-body models of clusters on average have only one three-body binary at any given time. This is different from theoretical expectations and models of binary populations, which predict a larger number of binaries ($\sim 5$), especially for low-$N$ clusters ($\sim 100$), or in the case of two-mass models, low number of black holes. We argue that the presence of multiple binaries is suppressed by a high rate of binary-binary interactions, which efficiently ionise one of the binaries involved. These also lead to triple formation and potentially gravitational wave (GW) captures, which may provide an explanation for the recently reported high efficiency of in-cluster mergers in models of low-mass clusters ($\lesssim 10^5\,{\rm M}_\odot)$. △ Less

Submitted 5 December, 2023; v1 submitted 20 August, 2023; originally announced August 2023.

Comments: 12 pages, 16 figures, accepted for publication in MNRAS

Binary vision: The merging black hole binary mass distribution via iterative density estimation

Authors: Jam Sadiq, Thomas Dent, Mark Gieles

Abstract: Binary black hole (BBH) systems detected via gravitational-wave (GW) emission are a recently opened astrophysical frontier with many unknowns and uncertainties. Accurate reconstruction of the binary distribution with as few assumptions as possible is desirable for inference on formation channels and environments. Most population analyses have, though, assumed a power law in binary mass ratio $q$,… ▽ More Binary black hole (BBH) systems detected via gravitational-wave (GW) emission are a recently opened astrophysical frontier with many unknowns and uncertainties. Accurate reconstruction of the binary distribution with as few assumptions as possible is desirable for inference on formation channels and environments. Most population analyses have, though, assumed a power law in binary mass ratio $q$, and/or assumed a universal $q$ distribution regardless of primary mass. Kernel density estimation (KDE)-based methods allow us to dispense with such assumptions and directly estimate the joint binary mass distribution. We deploy a self-consistent iterative method to estimate this full BBH mass distribution, finding local maxima in primary mass consistent with previous investigations and a secondary mass distribution with a partly independent structure, inconsistent with both power laws and with a constant function of $q$. We find a weaker preference for near-equal mass binaries than in most previous investigations; instead, the secondary mass has its own "spectral lines" at slightly lower values than the primary, and we observe an anti-correlation between primary and secondary masses around the ~$10M_\odot$ peak. △ Less

Submitted 22 July, 2023; originally announced July 2023.

Comments: 13 pages, 11 figures

Stellar-mass black holes in the Hyades star cluster?

Authors: Stefano Torniamenti, Mark Gieles, Zephyr Penoyre, Tereza Jerabkova, Long Wang, Friedrich Anders

Abstract: Astrophysical models of binary-black hole mergers in the Universe require a significant fraction of stellar-mass black holes (BHs) to receive negligible natal kicks to explain the gravitational wave detections. This implies that BHs should be retained even in open clusters with low escape velocities ($\lesssim1~\mathrm{km \, s^{-1}}$). We search for signatures of the presence of BHs in the nearest… ▽ More Astrophysical models of binary-black hole mergers in the Universe require a significant fraction of stellar-mass black holes (BHs) to receive negligible natal kicks to explain the gravitational wave detections. This implies that BHs should be retained even in open clusters with low escape velocities ($\lesssim1~\mathrm{km \, s^{-1}}$). We search for signatures of the presence of BHs in the nearest open cluster to the Sun - the Hyades - by comparing density profiles of direct $N$-body models to data from $Gaia$. The observations are best reproduced by models with $2-3$ BHs at present. Models that never possessed BHs have an half-mass radius $\sim30\%$ smaller than the observed value, while those where the last BHs were ejected recently ($\lesssim150~$Myr ago) can still reproduce the density profile. In 50% of the models hosting BHs, we find BHs with stellar companion(s). Their period distribution peaks at $\sim10^3$ yr, making them unlikely to be found through velocity variations. We look for potential BH companions through large $Gaia$ astrometric and spectroscopic errors, identifying 56 binary candidates - none of which consistent with a massive compact companion. Models with $2-3$ BHs have an elevated central velocity dispersion, but observations can not yet discriminate. We conclude that the present-day structure of the Hyades requires a significant fraction of BHs to receive natal kicks smaller than the escape velocity of $\sim 3\, \mathrm{km \, s^{-1}}$ at the time of BH formation and that the nearest BHs to the Sun are in, or near, Hyades. △ Less

Submitted 7 August, 2023; v1 submitted 17 March, 2023; originally announced March 2023.

Comments: 20 pages, 15 figures, 5 tables. Accepted for publication in MNRAS. Comments welcome

N-enhancement in GN-z11: First evidence for supermassive stars nucleosynthesis in proto-globular clusters-like conditions at high redshift ?

Authors: C. Charbonnel, D. Schaerer, N. Prantzos, L. Ramírez-Galeano, T. Fragos, A. Kuruvandothi, R. Marques-Chaves, M. Gieles

Abstract: Unusually high N/O abundance ratios were recently reported for a very compact, intensively star-forming object GN-z11 at $z=10.6$ from JWST/NIRSpec observations. We present an empirical comparison with the C, N, and O abundance ratios in Galactic globular clusters (GCs) over a large metallicity range. We show that hot hydrogen-burning nucleosynthesis within supermassive stars (SMS) formed through… ▽ More Unusually high N/O abundance ratios were recently reported for a very compact, intensively star-forming object GN-z11 at $z=10.6$ from JWST/NIRSpec observations. We present an empirical comparison with the C, N, and O abundance ratios in Galactic globular clusters (GCs) over a large metallicity range. We show that hot hydrogen-burning nucleosynthesis within supermassive stars (SMS) formed through runaway collisions can consistently explain the observed abundances ratio in GN-z11 and in GCs. This suggests that a proto-globular cluster hosting a SMS could be at the origin of the strong N-enrichment in GN-z11. Our model predicts the behavior of N/O, C/O, and Ne/O ratios as a function of metallicity, which can be tested if high-$z$ objects similar to GN-z11 are detected with JWST in the future. Further studies and statistics will help differentiate the proto-GC scenario from the Wolf-Rayet scenario that we quantify with a population synthesis model, and shed more light on this peculiar object. △ Less

Submitted 17 April, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: Accepted by A&A Letters (revised version; in press)

Journal ref: A&A 673, L7 (2023)

Updated radial velocities and new constraints on the nature of the unseen source in NGC1850 BH1

Authors: Sara Saracino, Tomer Shenar, Sebastian Kamann, Nate Bastian, Mark Gieles, Christopher Usher, Julia Bodensteiner, Angela Kochoska, Jerome A. Orosz, Hugues Sana

Abstract: A black hole candidate orbiting a luminous star in the Large Magellanic Cloud young cluster NGC 1850 ($\sim100$Myr) has recently been reported based on radial velocity and light curve modelling. Subsequently, an alternative explanation has been suggested for the system: a bloated post-mass transfer secondary star (M$_{\rm initial} \sim 4-5M_{\odot}$, M$_{\rm current} \sim 1-2M_{\odot}$) with a mor… ▽ More A black hole candidate orbiting a luminous star in the Large Magellanic Cloud young cluster NGC 1850 ($\sim100$Myr) has recently been reported based on radial velocity and light curve modelling. Subsequently, an alternative explanation has been suggested for the system: a bloated post-mass transfer secondary star (M$_{\rm initial} \sim 4-5M_{\odot}$, M$_{\rm current} \sim 1-2M_{\odot}$) with a more massive, yet luminous companion (the primary). Upon reanalysis of the MUSE spectra, we found that the radial velocity variations originally reported were underestimated ($K_{\rm 2,revised} = 176\pm3$km/s vs $K_{\rm 2,original} = 140\pm3$km/s) because of the weighting scheme adopted in the full-spectrum fitting analysis. The increased radial velocity semi-amplitude translates into a system mass function larger than previously deduced ($f_{\rm revised}$=2.83$M_{\odot}$ vs $f_{\rm original}$=1.42$M_{\odot}$). By exploiting the spectral disentangling technique, we place an upper limit of 10\% of a luminous primary source to the observed optical light in NGC1850 BH1, assuming that the primary and secondary are the only components contributing to the system. Furthermore, by analysing archival near-infrared data, we find clues to the presence of an accretion disk in the system. These constraints support a low-mass post-mass transfer star but do not provide a definitive answer whether the unseen component in NGC1850 BH1 is indeed a black hole. These results predict a scenario where, if a primary luminous source of mass M $\ge 4.7M_{\odot}$, is present in the system (given the inclination and secondary mass constraints), it must be hidden in a optically thick disk to be undetected in the MUSE spectra. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 10 pages, 8 Figures and 2 Tables. Accepted for publication by MNRAS

The mass-loss rates of star clusters with stellar-mass black holes: implications for the globular cluster mass function

Authors: Mark Gieles, Oleg Y. Gnedin

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 parameterised mass-loss rate, benchmarked by results from $N$-body simulations, and use it to evolve an initial GC mass function (GCMF), sim… ▽ More 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 parameterised mass-loss rate, benchmarked 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]$\lesssim-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 ($\sim10^5\,{\rm M}_\odot$) at large Galactic radii ($\gtrsim 10\,$kpc). The turnover mass at smaller Galactic radii is similar as the result of the upper mass truncation of the initial GCMF and the lower mass-loss rate because of the higher metallicities. The density profile in the Galaxy of mass lost from massive GCs ($\gtrsim10^{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. △ Less

Submitted 26 April, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

Comments: 18 pages, 16 figures, accepted for MNRAS

Multimass modelling of Milky Way globular clusters -- I. Implications on their stellar initial mass function above 1 M$_{\odot}$

Authors: Nolan Dickson, Vincent Hénault-Brunet, Holger Baumgardt, Mark Gieles, Peter Smith

Abstract: 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 var… ▽ More 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 ($α_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. △ Less

Submitted 24 April, 2023; v1 submitted 2 March, 2023; originally announced March 2023.

Comments: 21 pages, 10 figures, MNRAS accepted

The multiplicity fraction in 202 open clusters from Gaia

Authors: J. Donada, F. Anders, C. Jordi, E. Masana, M. Gieles, G. I. Perren, L. Balaguer-Núñez, A. Castro-Ginard, T. Cantat-Gaudin, L. Casamiquela

Abstract: In this study, we estimate the fraction of binaries with high mass ratios for 202 open clusters in the extended solar neighbourhood (closer than 1.5 kpc from the Sun). This is one of the largest homogeneous catalogues of multiplicity fractions in open clusters to date, including the unresolved and total (close-binary) multiplicity fractions of main-sequence systems with mass ratio larger than… ▽ More In this study, we estimate the fraction of binaries with high mass ratios for 202 open clusters in the extended solar neighbourhood (closer than 1.5 kpc from the Sun). This is one of the largest homogeneous catalogues of multiplicity fractions in open clusters to date, including the unresolved and total (close-binary) multiplicity fractions of main-sequence systems with mass ratio larger than $0.6_{-0.15}^{+0.05}$. The unresolved multiplicity fractions are estimated applying a flexible mixture model to the observed Gaia colour-magnitude diagrams of the open clusters. Then we use custom Gaia simulations to account for the resolved systems and derive the total multiplicity fractions. The studied open clusters have ages between 6.6 Myr and 3.0 Gyr and total high-mass-ratio multiplicity fractions between 6% and 80%, with a median of 18%. The multiplicity fractions increase with the mass of the primary star, as expected. The average multiplicity fraction per cluster displays an overall decreasing trend with the open cluster age up to ages about 100 Myr, above which the trend increases. Our simulations show that most of this trend is caused by complex selection effects (introduced by the mass dependence of the multiplicity fraction and the magnitude limit of our sample). Furthermore, the multiplicity fraction is not significantly correlated with the clusters' position in the Galaxy. The spread in multiplicity fraction decreases significantly with the number of cluster members (used as a proxy for cluster mass). We also find that the multiplicity fraction decreases with metallicity, in line with recent studies using field stars. △ Less

Submitted 1 May, 2023; v1 submitted 26 January, 2023; originally announced January 2023.

Comments: 17 pages, 13 figures, accepted for publication in A&A after minor modifications

Journal ref: A&A 675, A89 (2023)

Dynamics in the outskirts of four Milky Way globular clusters: it's the tides that dominate

Authors: Zhen Wan, Anthony D. Arnold, William H. Oliver, Geraint F. Lewis, Holger Baumgardt, Mark Gieles, Vincent Hénault-Brunet, Thomas de Boer, Eduardo Balbinot, Gary Da Costa, Dougal Mackey, Denis Erkal, Annette Ferguson, Pete Kuzma, Elena Pancino, Jorge Penarrubia, Nicoletta Sanna, Antonio Sollima, Roeland P. van der Marel, Laura L. Watkins

Abstract: We present the results of a spectroscopic survey of the outskirts of 4 globular clusters -- NGC 1261, NGC 4590, NGC 1904, and NGC 1851 -- covering targets within 1 degree from the cluster centres, with 2dF/AAOmega on the Anglo-Australian Telescope (AAT) and FLAMES on the Very Large Telescope (VLT). We extracted chemo-dynamical information for individual stars, from which we estimated the veloc… ▽ More We present the results of a spectroscopic survey of the outskirts of 4 globular clusters -- NGC 1261, NGC 4590, NGC 1904, and NGC 1851 -- covering targets within 1 degree from the cluster centres, with 2dF/AAOmega on the Anglo-Australian Telescope (AAT) and FLAMES on the Very Large Telescope (VLT). We extracted chemo-dynamical information for individual stars, from which we estimated the velocity dispersion profile and the rotation of each cluster. The observations are compared to direct $N$-body simulations and appropriate {\sc limepy}/{\sc spes} models for each cluster to interpret the results. In NGC 1851, the detected internal rotation agrees with existing literature, and NGC 1261 shows some rotation signal beyond the truncation radius, likely coming from the escaped stars. We find that the dispersion profiles for both the observations and the simulations for NGC 1261, NGC 1851, and NGC 1904 do not decrease as the {\sc limepy}/{\sc spes} models predict beyond the truncation radius, where the $N$-body simulations show that escaped stars dominate; the dispersion profile of NGC 4590 follows the predictions of the {\sc limepy}/{\sc spes} models, though the data do not effectively extend beyond the truncation radius. The increasing/flat dispersion profiles in the outskirts of NGC 1261, NGC 1851 and NGC 1904, are reproduced by the simulations. Hence, the increasing/flat dispersion profiles of the clusters in question can be explained by the tidal interaction with the Galaxy without introducing dark matter. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: 17 pages, 10 figures. Accepted by MNRAS. Comments are welcome

The VLT-FLAMES Tarantula Survey: Observational evidence for two distinct populations of massive runaway stars in 30 Doradus

Authors: H. Sana, O. H. Ramírez-Agudelo, V. Hénault-Brunet, L. Mahy, L. A. Almeida, A. de Koter, J. M. Bestenlehner, C. J. Evans, N. Langer, F. R. N. Schneider, P. A. Crowther, S. E. de Mink, A. Herrero, D. J. Lennon, M. Gieles, J. Maíz Apellániz, M. Renzo, E. Sabbi, J. Th. van Loon, J. S. Vink

Abstract: Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active… ▽ More Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active star-forming region in the Local group. We use RV measurements of the O-type star populations in 30 Doradus obtained by the VLT-FLAMES Tarantula Survey and the Tarantula Massive Binary Monitoring to identify single and binary O-type runaways. We discuss their rotational properties and qualitatively compare observations with expectations of ejection scenarios. We identify 23 single and one binary O-type runaway objects, most of them outside the main star-forming regions in 30 Doradus. We find an overabundance of rapid rotators (vsini > 200km/s) among the runaway population, providing an explanation of the overabundance of rapidly rotating stars in the 30 Doradus field. Considerations of the projected rotation rates and runaway line-of-sight (los) velocities reveal a conspicuous absence of rapidly rotating (vsini > 210k/ms), fast moving (v_{los} > 60km/s) runaways, and suggest the presence of two different populations of runaway stars: a population of rapidly-spinning but slowly moving runaways and a population of fast moving but slowly rotating ones. These are detected with a ratio close to 2:1 in our sample. We argue that slowly moving but rapidly spinning runaways result from binary ejections, while rapidly moving but slowly spinning runaways could result from dynamical ejections. Given that detection biases will more strongly impact the slow-moving population, our results suggest that the binary evolution scenario dominates the current massive runaway population in 30 Doradus. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A Letters; 9 pages, 5 figures

Journal ref: A&A 668, L5 (2022)

Coalescing black hole binaries from globular clusters: mass distributions and comparison to gravitational wave data from GWTC-3

Authors: Fabio Antonini, Mark Gieles, Fani Dosopoulou, Debatri Chattopadhyay

Abstract: We use our cluster population model, cBHBd, to explore the mass distribution of merging black hole binaries formed dynamically in globular clusters. We include in our models the effect of mass growth through hierarchical mergers and compare the resulting distributions to those inferred from the third gravitational wave transient catalogue. We find that none of our models can reproduce the peak at… ▽ More We use our cluster population model, cBHBd, to explore the mass distribution of merging black hole binaries formed dynamically in globular clusters. We include in our models the effect of mass growth through hierarchical mergers and compare the resulting distributions to those inferred from the third gravitational wave transient catalogue. We find that none of our models can reproduce the peak at $m_1\simeq 10M_\odot$ in the primary black hole mass distribution that is inferred from the data. This disfavours a scenario where most of the sources are formed in globular clusters. On the other hand, a globular cluster origin can account for the inferred secondary peak at $m_1\simeq 35M_\odot$, which requires that the most massive clusters form with half-mass densities $ρ_{\rm h,0} \gtrsim 10^4 M_\odot \rm pc^{-3}$. Finally, we find that the lack of a high mass cut--off in the inferred mass distribution can be also explained by the repopulation of an initial mass gap through hierarchical mergers. Matching the inferred merger rate above $\simeq 50M_\odot$ requires both initial cluster densities $ρ_{\rm h,0} \gtrsim 10^4 M_\odot \rm pc^{-3}$, and that black holes form with nearly zero spin. A hierarchical merger scenario makes specific predictions for the appearance and position of multiple peaks in the black hole mass distribution, which can be tested against future data. △ Less

Submitted 29 March, 2023; v1 submitted 1 August, 2022; originally announced August 2022.

Comments: accepted by MNRAS; 11 pages, 5 figures

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Authors: Tomer Shenar, Hugues Sana, Laurent Mahy, Kareem El-Badry, Pablo Marchant, Norbert Langer, Calum Hawcroft, Matthias Fabry, Koushik Sen, Leonardo A. Almeida, Michael Abdul-Masih, Julia Bodensteiner, Paul A. Crowther, Mark Gieles, Mariusz Gromadzki, Vincent Henault-Brunet, Artemio Herrero, Alex de Koter, Patryk Iwanek, Szymon Kozłowski, Daniel J. Lennon, Jesus Maız Apellaniz, Przemysław Mroz, Anthony F. J. Moffat, Annachiara Picco , et al. (13 additional authors not shown)

Abstract: Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-ty… ▽ More Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: Accepted to Nature Astronomy, 64 pages, 15 figures, 2 tables; ESO press release: https://www.eso.org/public/news/eso2210/; Nat Asr paper URL: https://www.nature.com/articles/s41550-022-01730-y

A black hole detected in the young massive LMC cluster NGC 1850

Authors: S. Saracino, S. Kamann, M. G. Guarcello, C. Usher, N. Bastian, I. Cabrera-Ziri, M. Gieles, S. Dreizler, G. S. Da Costa, T. -O. Husser, V. Hénault-Brunet

Abstract: We report the detection of a black hole (NGC 1850 BH1) in the $\sim$100 Myr-old massive cluster NGC~1850 in the Large Magellanic Cloud. It is in a binary system with a main-sequence turn-off star (4.9 $\pm$ 0.4 M${_\odot}$), which is starting to fill its Roche Lobe and becoming distorted. Using 17 epochs of VLT/MUSE observations we detected radial velocity variations exceeding 300 km/s associated… ▽ More We report the detection of a black hole (NGC 1850 BH1) in the $\sim$100 Myr-old massive cluster NGC~1850 in the Large Magellanic Cloud. It is in a binary system with a main-sequence turn-off star (4.9 $\pm$ 0.4 M${_\odot}$), which is starting to fill its Roche Lobe and becoming distorted. Using 17 epochs of VLT/MUSE observations we detected radial velocity variations exceeding 300 km/s associated to the target star, linked to the ellipsoidal variations measured by OGLE-IV in the optical bands. Under the assumption of a semi-detached system, the simultaneous modelling of radial velocity and light curves constraints the orbital inclination of the binary to ($38 \pm 2$)$^{\circ}$, resulting in a true mass of the unseen companion of $11.1_{-2.4}^{+2.1}$ $M_{\odot}$. This represents the first direct dynamical detection of a black hole in a young massive cluster, opening up the possibility of studying the initial mass function and the early dynamical evolution of such compact objects in high-density environments. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 11 pages, 8 figures, 4 Tables. Accepted for publication by Monthly Notices of the Royal Astronomical Society

MUSE narrow field mode observations of the central kinematics of M15

Authors: Christopher Usher, Sebastian Kamann, Mark Gieles, Vincent Hénault-Brunet, Emanuele Dalessandro, Eduardo Balbinot, Antonio Sollima

Abstract: We present observations of the stellar kinematics of the centre of the core collapsed globular cluster M15 obtained with the MUSE integral field spectrograph on the VLT operating in narrow field mode. Thanks to the use of adaptive optics, we obtain a spatial resolution of 0.1arcsec and are able to reliably measure the radial velocities of 864 stars within 8 arcsec of the centre of M15 thus providi… ▽ More We present observations of the stellar kinematics of the centre of the core collapsed globular cluster M15 obtained with the MUSE integral field spectrograph on the VLT operating in narrow field mode. Thanks to the use of adaptive optics, we obtain a spatial resolution of 0.1arcsec and are able to reliably measure the radial velocities of 864 stars within 8 arcsec of the centre of M15 thus providing the largest sample of radial velocities ever obtained for the innermost regions of this system. Combined with previous observations of M15 using MUSE in wide field mode and literature data, we find that the central kinematics of M15 are complex with the rotation axis of the core of M15 offset from the rotation axis of the bulk of the cluster. While this complexity has been suggested by previous work, we confirm it at higher significance and in more detail. △ Less

Submitted 24 February, 2021; v1 submitted 23 February, 2021; originally announced February 2021.

Comments: 8 pages, 6 figures. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

A supra-massive population of stellar-mass black holes in the globular cluster Palomar 5

Authors: Mark Gieles, Denis Erkal, Fabio Antonini, Eduardo Balbinot, Jorge Peñarrubia

Abstract: Palomar 5 is one of the sparsest star clusters in the Galactic halo and is best-known for its spectacular tidal tails, spanning over 20 degrees across the sky. With N-body simulations we show that both distinguishing features can result from a stellar-mass black hole population, comprising ~20% of the present-day cluster mass. In this scenario, Palomar 5 formed with a `normal' black hole mass frac… ▽ More Palomar 5 is one of the sparsest star clusters in the Galactic halo and is best-known for its spectacular tidal tails, spanning over 20 degrees across the sky. With N-body simulations we show that both distinguishing features can result from a stellar-mass black hole population, comprising ~20% of the present-day cluster mass. In this scenario, Palomar 5 formed with a `normal' black hole mass fraction of a few per cent, but stars were lost at a higher rate than black holes, such that the black hole fraction gradually increased. This inflated the cluster, enhancing tidal stripping and tail formation. A gigayear from now, the cluster will dissolve as a 100% black hole cluster. Initially denser clusters end up with lower black hole fractions, smaller sizes, and no observable tails. Black hole-dominated, extended star clusters are therefore the likely progenitors of the recently discovered thin stellar streams in the Galactic halo. △ Less

Submitted 4 July, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: 33 pages, 7 figures, 3 tables, accepted for Nature Astronomy (update to match final published version to appear on 5 July 2021)

The dynamics of the globular cluster NGC3201 out to the Jacobi radius

Authors: Zhen Wan, William Oliver, Holger Baumgardt, Geraint Lewis, Mark Gieles, Vincent Hénault-Brunet, Thomas de Boer, Eduardo Balbinot, Gary Da Costa, Dougal Mackey

Abstract: As part of a chemo-dynamical survey of five nearby globular clusters with 2dF/AAOmega on the AAT, we have obtained kinematic information for the globular cluster NGC3201. Our new observations confirm the presence of a significant velocity gradient across the cluster which can almost entirely be explained by the high proper motion of the cluster. After subtracting the contribution of this perspecti… ▽ More As part of a chemo-dynamical survey of five nearby globular clusters with 2dF/AAOmega on the AAT, we have obtained kinematic information for the globular cluster NGC3201. Our new observations confirm the presence of a significant velocity gradient across the cluster which can almost entirely be explained by the high proper motion of the cluster. After subtracting the contribution of this perspective rotation, we found a remaining rotation signal with an amplitude of $\sim1\ km/s$ around a different axis to what we expect from the tidal tails and the potential escapers, suggesting that this rotation is internal and can be a remnant of its formation process. At the outer part, we found a rotational signal that is likely a result from potential escapers. The proper motion dispersion at large radii reported by Bianchini et al. has previously been attributed to dark matter. Here we show that the LOS dispersion between 0.5-1 Jacobi radius is lower, yet above the predictions from an N-body model of NGC3201 that we ran for this study. Based on the simulation, we find that potential escapers cannot fully explain the observed velocity dispersion. We also estimate the effect on the velocity dispersion of different amounts of stellar-mass black holes and unbound stars from the tidal tails with varying escape rates and find that these effects cannot explain the difference between the LOS dispersion and the N-body model. Given the recent discovery of tidal tail stars at large distances from the cluster, a dark matter halo is an unlikely explanation. We show that the effect of binary stars, which is not included in the N-body model, is important and can explain part of the difference in dispersion. We speculate that the remaining difference must be the result of effects not included in the N-body model, such as initial cluster rotation, velocity anisotropy and Galactic substructure. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: 14 pages, 13 figures, accepted by MNRAS

On the response of a star cluster to a tidal perturbation

Authors: Luis A. Martinez-Medina, Mark Gieles, Oleg Y. Gnedin, Hui Li

Abstract: We study the response of star clusters to individual tidal perturbations using controlled $N$-body simulations. We consider perturbations by a moving point mass and by a disc, and vary the duration of the perturbation as well as the cluster density profile. For fast perturbations (i.e. `shocks'), the cluster gains energy in agreement with theoretical predictions in the impulsive limit. For slow di… ▽ More We study the response of star clusters to individual tidal perturbations using controlled $N$-body simulations. We consider perturbations by a moving point mass and by a disc, and vary the duration of the perturbation as well as the cluster density profile. For fast perturbations (i.e. `shocks'), the cluster gains energy in agreement with theoretical predictions in the impulsive limit. For slow disc perturbations, the energy gain is lower, and this has previously been attributed to adiabatic damping. However, the energy gain due to slow perturbations by a point-mass is similar to, or larger than that due to fast shocks, which is not expected because adiabatic damping should be almost independent of the nature of the tides. We show that the geometric distortion of the cluster during slow perturbations is of comparable importance for the energy gain as adiabatic damping, and that the combined effect can qualitatively explain the results. The half-mass radius of the bound stars after a shock increases up to $\sim$7\% for low-concentration clusters, and decreases $\sim$3\% for the most concentrated ones. The fractional mass loss is a non-linear function of the energy gain, and depends on the nature of the tides and most strongly on the cluster density profile, making semi-analytic model predictions for cluster lifetimes extremely sensitive to the adopted density profile. △ Less

Submitted 10 August, 2022; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 13 pages, 15 figures. Accepted for publication in MNRAS

Merger rate of black hole binaries from globular clusters: theoretical error bars and comparison to gravitational wave data from GWTC-2

Authors: Fabio Antonini, Mark Gieles

Abstract: Black hole binaries formed dynamically in globular clusters are believed to be one of the main sources of gravitational waves in the Universe. Here, we use our new population synthesis code, cBHBd, to determine the redshift evolution of the merger rate density and masses of black hole binaries formed in globular clusters. We simulate $\sim 2$ million models to explore the parameter space that is r… ▽ More Black hole binaries formed dynamically in globular clusters are believed to be one of the main sources of gravitational waves in the Universe. Here, we use our new population synthesis code, cBHBd, to determine the redshift evolution of the merger rate density and masses of black hole binaries formed in globular clusters. We simulate $\sim 2$ million models to explore the parameter space that is relevant to real clusters and over all mass scales. We show that when uncertainties on the initial cluster mass function and density are properly taken into account, they become the two dominant factors in setting the theoretical error bars on merger rates. Other model parameters (e.g., natal kicks, black hole masses, metallicity) have virtually no effect on the local merger rate density, although they affect the masses of the merging black holes. Modelling the merger rate density as a function of redshift as $R(z)=R_0(1+z)^κ$ at $z<2$, and marginalizing over uncertainties, we find: $R_0=7.2^{+21.5}_{-5.5}{\rm Gpc^{-3}yr^{-1}}$ and $κ=1.6^{+0.4}_{-0.6}$ ($90\%$ credibility). The rate parameters for binaries that merge inside the clusters are ${R}_{\rm 0,in}=1.6^{+1.9}_{-1.0}{\rm Gpc^{-3}yr^{-1}}$ and $κ_{\rm in}=2.3^{+1.3}_{-1.0}$; $\sim 20\%$ of these form as the result of a gravitational-wave capture, implying that eccentric mergers from globular clusters contribute $\lesssim 0.4 \rm Gpc^{-3}yr^{-1}$ to the local rate. A comparison to the merger rate reported by LIGO-Virgo shows that a scenario in which most of the detected black hole mergers are formed in globular clusters is consistent with current constraints, and requires initial cluster half-mass densities $\gtrsim 10^4 M_\odot \rm pc^{-3}$. Such models also reproduce the inferred primary black hole mass distribution for masses $13-30 M_\odot$, but under-predict the data outside this range. △ Less

Submitted 7 November, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: Revised to include a comparison to the data from GWTC-2. Accepted in PRD; 23 pages, 13 figures

Journal ref: Phys. Rev. D 102, 123016 (2020)

Weighing stars from birth to death: mass determination methods across the HRD

Authors: Aldo Serenelli, Achim Weiss, Conny Aerts, George C. Angelou, David Baroch, Nate Bastian, Paul G. Beck, Maria Bergemann, Joachim M. Bestenlehner, Ian Czekala, Nancy Elias-Rosa, Ana Escorza, Vincent Van Eylen, Diane K. Feuillet, Davide Gandolfi, Mark Gieles, Leo Girardi, Yveline Lebreton, Nicolas Lodieu, Marie Martig, Marcelo M. Miller Bertolami, Joey S. G. Mombarg, Juan Carlos Morales, Andres Moya, Benard Nsamba , et al. (9 additional authors not shown)

Abstract: The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approac… ▽ More The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,\msun$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars. △ Less

Submitted 9 April, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: Invited review article for The Astronomy and Astrophysics Review. 146 pages, 16 figures, 11 tables. Accepted version by the Journal. It includes summary figure of accuracy/precision of methods for mass ranges and summary table for individual methods

The star cluster age function in the Galactic disc with Gaia DR2: Fewer old clusters and a low cluster formation efficiency

Authors: Friedrich Anders, Tristan Cantat-Gaudin, Irene Quadrino-Lodoso, Mark Gieles, Carme Jordi, Alfred Castro-Ginard, Lola Balaguer-Núñez

Abstract: We perform a systematic reanalysis of the age distribution of Galactic open star clusters. Using a catalogue of homogeneously determined ages for 834 open clusters contained in a 2 kpc cylinder around the Sun and characterised with astrometric and photometric data from the Gaia satellite, we find that it is necessary to revise earlier works that relied on data from the Milky Way Star Cluster surve… ▽ More We perform a systematic reanalysis of the age distribution of Galactic open star clusters. Using a catalogue of homogeneously determined ages for 834 open clusters contained in a 2 kpc cylinder around the Sun and characterised with astrometric and photometric data from the Gaia satellite, we find that it is necessary to revise earlier works that relied on data from the Milky Way Star Cluster survey. After establishing age-dependent completeness limits for our sample, we find that the cluster age function in the range $6.5 < \log t<10$ is compatible with a Schechter-type or broken power-law function, whose parameters we determine by MCMC fitting. Our best-fit values indicate an earlier drop of the age function (by a factor of $2-3$) with respect to the results obtained in the last five years, and are instead more compatible with results obtained in the early 2000s and radio observations of inner-disc clusters. Furthermore, we find a typical destruction time-scale of $\sim1.5$ Gyr for a $10^4\, {\rm M}_{\odot}$ cluster and a present-day cluster-formation rate of $0.55_{-0.15}^{+0.19}$ Myr$^{-1}$kpc$^{-2}$, suggesting that only $16_{-8}^{+11}$ \% of all stars born in the solar neighbourhood form in bound clusters. Accurate cluster-mass measurements are now needed to place more precise constraints on open-cluster formation and evolution models. △ Less

Submitted 4 December, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: 4 pages + 2 pages appendix, 4 figures, A&A Letters, in press. Notebooks including supplementary analysis: https://github.com/fjaellet/gaiadr2-caf

Journal ref: A&A 645, L2 (2021)

Star clusters near and far; tracing star formation across cosmic time

Authors: Angela Adamo, Peter Zeidler, J. M. Diederik Kruijssen, Mélanie Chevance, Mark Gieles, Daniela Calzetti, Corinne Charbonnel, Hans Zinnecker, Martin G. H. Krause

Abstract: Star clusters are fundamental units of stellar feedback and unique tracers of their host galactic properties. In this review, we will first focus on their constituents, i.e.\ detailed insight into their stellar populations and their surrounding ionised, warm, neutral, and molecular gas. We, then, move beyond the Local Group to review star cluster populations at various evolutionary stages, and in… ▽ More Star clusters are fundamental units of stellar feedback and unique tracers of their host galactic properties. In this review, we will first focus on their constituents, i.e.\ detailed insight into their stellar populations and their surrounding ionised, warm, neutral, and molecular gas. We, then, move beyond the Local Group to review star cluster populations at various evolutionary stages, and in diverse galactic environmental conditions accessible in the local Universe. At high redshift, where conditions for cluster formation and evolution are more extreme, we are only able to observe the integrated light of a handful of objects that we believe will become globular clusters. We therefore discuss how numerical and analytical methods, informed by the observed properties of cluster populations in the local Universe, are used to develop sophisticated simulations potentially capable of disentangling the genetic map of galaxy formation and assembly that is carried by globular cluster populations. △ Less

Submitted 3 June, 2020; v1 submitted 13 May, 2020; originally announced May 2020.

Comments: Final published version. 61 pages, 10 figures. To appear in Space Science Reviews (submitted January 31, 2020; accepted May 8, 2020), topical collection Star formation

The Physics of Star Cluster Formation and Evolution

Authors: Martin G. H. Krause, Stella S. R. Offner, Corinne Charbonnel, Mark Gieles, Ralf S. Klessen, Enrique Vazquez-Semadeni, Javier Ballesteros-Paredes, Philipp Girichidis, J. M. Diederik Kruijssen, Jacob L. Ward, Hans Zinnecker

Abstract: Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively… ▽ More Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work. △ Less

Submitted 15 May, 2020; v1 submitted 2 May, 2020; originally announced May 2020.

Comments: 53 pages, 7 figures, preprint of a review to appear in Space Science Reviews Topical collection "Star formation". Changes in this version: updated comments line to to acceptance of the paper by the journal

Time-domain study of the young massive cluster Westerlund 2 with the Hubble Space Telescope. I

Authors: E. Sabbi, M. Gennaro, J. Anderson, V. Bajaj, N. Bastian, J. S. Gallagher, III, M. Gieles, D. J. Lennon, A. Nota, K. C. Sahu, P. Zeidler

Abstract: Time-domain studies of pre-main sequence stars have long been used to investigate star properties during their early evolutionary phases and to trace the evolution of circumstellar environments. Historically these studies have been confined to the nearest, low-density, star forming regions. We used the Wide Field Camera 3 on board of the Hubble Space Telescope to extend, for the first time, the st… ▽ More Time-domain studies of pre-main sequence stars have long been used to investigate star properties during their early evolutionary phases and to trace the evolution of circumstellar environments. Historically these studies have been confined to the nearest, low-density, star forming regions. We used the Wide Field Camera 3 on board of the Hubble Space Telescope to extend, for the first time, the study of pre-main sequence variability to one of the few young massive clusters in the Milky Way, Westerlund 2. Our analysis reveals that at least 1/3 of the intermediate and low-mass pre-main sequence stars in Westerlund 2 are variable. Based on the characteristics of their light curves, we classified ~11% of the variable stars as weak-line T-Tauri candidates, ~ 52% as classical T-Tauri candidates, ~ 5% as dippers and ~26% as bursters. In addition, we found that 2% of the stars below 6Mo (~6% of the variables) are eclipsing binaries, with orbital periods shorter than 80 days. The spatial distribution of the different populations of variable pre-main sequence stars suggests that stellar feedback and UV-radiation from massive stars play an important role on the evolution of circumstellar and planetary disks. △ Less

Submitted 11 February, 2020; originally announced February 2020.

Comments: 19 pages, 15 figures, accepted for publication on ApJ

A closer look at the spur, blob, wiggle, and gaps in GD-1

Authors: T. J. L. de Boer, D. Erkal, M. Gieles

Abstract: The GD-1 stream is one of the longest and coldest stellar streams discovered to date, and one of the best objects for constraining the dark matter properties of the Milky Way. Using data from {\it Gaia} DR2 we study the proper motions, distance, morphology and density of the stream to uncover small scale perturbations. The proper motion cleaned data shows a clear distance gradient across the strea… ▽ More The GD-1 stream is one of the longest and coldest stellar streams discovered to date, and one of the best objects for constraining the dark matter properties of the Milky Way. Using data from {\it Gaia} DR2 we study the proper motions, distance, morphology and density of the stream to uncover small scale perturbations. The proper motion cleaned data shows a clear distance gradient across the stream, ranging from 7 to 12 kpc. However, unlike earlier studies that found a continuous gradient, we uncover a distance minimum at $\varphi_{1}\approx$-50 deg, after which the distance increases again. We can reliably trace the stream between -85$<\varphi_{1}<$15 deg, showing an even further extent to GD-1 beyond the earlier extension of \citet{Price-Whelan18a}. We constrain the stream track and density using a Boolean matched filter approach and find three large under densities and find significant residuals in the stream track lining up with these gaps. In particular, a gap is visible at $\varphi_{1}$=-3 deg, surrounded by a clear sinusoidal wiggle. We argue that this wiggle is due to a perturbation since it has the wrong orientation to come from a progenitor. We compute a total initial stellar mass of the stream segment of 1.58$\pm$0.07$\times$10$^{4}$ M$_{\odot}$. With the extended view of the spur in this work, we argue that the spur may be unrelated to the adjacent gap in the stream. Finally, we show that an interaction with the Sagittarius dwarf can create features similar to the spur. △ Less

Submitted 13 November, 2019; originally announced November 2019.

Comments: Submitted to MNRAS. comments are very welcome

On the black hole content and initial mass function of 47 Tuc

Authors: Vincent Hénault-Brunet, Mark Gieles, Jay Strader, Miklos Peuten, Eduardo Balbinot, Kaela E. K. Douglas

Abstract: The globular cluster (GC) 47 Tuc has recently been proposed to host an intermediate-mass black hole (IMBH) or a population of stellar-mass black holes (BHs). To shed light on its dark content, we present an application of self-consistent multimass models with a varying mass function and content of stellar remnants, which we fit to various observational constraints. Our best-fitting model successfu… ▽ More The globular cluster (GC) 47 Tuc has recently been proposed to host an intermediate-mass black hole (IMBH) or a population of stellar-mass black holes (BHs). To shed light on its dark content, we present an application of self-consistent multimass models with a varying mass function and content of stellar remnants, which we fit to various observational constraints. Our best-fitting model successfully matches the observables and correctly predicts the radial distribution of millisecond pulsars and their gravitational accelerations inferred from long-term timing observations. The data favours a population of BHs with a total mass of $430^{+386}_{-301}$ $M_{\odot}$, but the most likely model has very few BHs. Since our models do not include a central IMBH and accurately reproduce the observations, we conclude that there is currently no need to invoke the presence of an IMBH in 47 Tuc. The global present-day mass function inferred is significantly depleted in low-mass stars (power-law slope $α=-0.52^{+0.17}_{-0.16}$). Given the orbit and predicted mass-loss history of this massive GC, the dearth of low-mass stars is difficult to explain with a standard initial mass function (IMF) followed by long-term preferential escape of low-mass stars driven by two-body relaxation, and instead suggests that 47 Tuc may have formed with a bottom-light IMF. We discuss alternative evolutionary origins for the flat mass function and ways to reconcile this with the low BH retention fraction. Finally, by capturing the effect of dark remnants, our method offers a new way to probe the IMF in a GC above the current main-sequence turn-off mass, for which we find a slope of $-2.49\pm0.08$. △ Less

Submitted 25 October, 2019; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: 16 pages, 7 figures, accepted to MNRAS after minor revision

HAYDN -- High-precision AsteroseismologY of DeNse stellar fields (ESA Voyage 2050 White Paper)

Authors: Andrea Miglio, Leo Girardi, Frank Grundahl, Benoit Mosser, Nate Bastian, Angela Bragaglia, Karsten Brogaard, Gael Buldgen, William Chantereau, Bill Chaplin, Cristina Chiappini, Marc-Antoine Dupret, Patrick Eggenberger, Mark Gieles, Rob Izzard, Daisuke Kawata, Christoffer Karoff, Nadege Lagarde, Ted Mackereth, Demetrio Magrin, Georges Meynet, Eric Michel, Josefina Montalban, Valerio Nascimbeni, Arlette Noels , et al. (7 additional authors not shown)

Abstract: In the last decade, the Kepler and CoRoT space-photometry missions have demonstrated the potential of asteroseismology as a novel, versatile and powerful tool to perform exquisite tests of stellar physics, and to enable precise and accurate characterisations of stellar properties, with impact on both exoplanetary and Galactic astrophysics. Based on our improved understanding of the strengths and l… ▽ More In the last decade, the Kepler and CoRoT space-photometry missions have demonstrated the potential of asteroseismology as a novel, versatile and powerful tool to perform exquisite tests of stellar physics, and to enable precise and accurate characterisations of stellar properties, with impact on both exoplanetary and Galactic astrophysics. Based on our improved understanding of the strengths and limitations of such a tool, we argue for a new small/medium space mission dedicated to gathering high-precision, high-cadence, long photometric series in dense stellar fields. Such a mission will lead to breakthroughs in stellar astrophysics, especially in the metal poor regime, will elucidate the evolution and formation of open and globular clusters, and aid our understanding of the assembly history and chemodynamics of the Milky Way's bulge and few nearby dwarf galaxies. △ Less

Submitted 7 April, 2021; v1 submitted 14 August, 2019; originally announced August 2019.

Comments: White paper submitted in response to Voyage 2050 long-term plan in the ESA Science Programme, now published in Experimental Astronomy, see https://link.springer.com/article/10.1007/s10686-021-09711-1

Journal ref: Exp Astron (2021)

Massive stars in extremely metal-poor galaxies: A window into the past

Authors: M. Garcia, C. J. Evans, J. M. Bestenlehner, J. C. Bouret, N. Castro, M. Cerviño, A. W. Fullerton, M. Gieles, A. Herrero, A. de Koter, D. J. Lennon, J. Th. van Loon, F. Martins, S. E. de Mink, F. Najarro, I. Negueruela, H. Sana, S. Simón-Díaz, D. Szécsi, F. Tramper, J. Vink, A. Wofford

Abstract: Cosmic History has witnessed the lives and deaths of multiple generations of massive stars, all of them invigorating their host galaxies with ionizing photons, kinetic energy, fresh material and stellar-mass black holes. Ubiquitous engines as they are, Astrophysics needs a good understanding of their formation, evolution, properties and yields throughout the history of the Universe, and with decre… ▽ More Cosmic History has witnessed the lives and deaths of multiple generations of massive stars, all of them invigorating their host galaxies with ionizing photons, kinetic energy, fresh material and stellar-mass black holes. Ubiquitous engines as they are, Astrophysics needs a good understanding of their formation, evolution, properties and yields throughout the history of the Universe, and with decreasing metal content mimicking the environment at the earliest epochs. Ultimately, a physical model that could be extrapolated to zero metallicity would enable tackling long-standing questions such as "What did the First, very massive stars of the Universe look like?" or "What was their role in the re-ionization of the Universe?". Yet, most our knowledge of metal-poor massive stars is drawn from one single point in metallicity. Massive stars in the Small Magellanic Cloud (SMC, $\sim 1/5 Z_{\odot}$) currently serve as templates for low-metallicity objects in the early Universe, even though significant differences with respect to massive stars with poorer metal content have been reported. This White Paper summarizes the current knowledge on extremely (sub-SMC) metal poor massive stars, highlighting the most outstanding open questions and the need to supersede the SMC as standard. A new paradigm can be built from nearby extremely metal-poor galaxies that make a new metallicity ladder, but massive stars in these galaxies are out of reach to current observational facilities. Such task would require an L-size mission, consisting of a 10m-class space telescope operating in the optical and the ultraviolet ranges. Alternatively, we propose that ESA unites efforts with NASA to make the LUVOIR mission concept a reality, thus continuing the successful partnership that made Hubble Space Telescope one of the greatest observatories of all time. △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: A white paper submitted for the Voyage 2050 long-term plan in the ESA Science Programme. 21 pages, 1 table, 6 figures. arXiv admin note: substantial text overlap with arXiv:1903.05235

Supermassive stars as the origin of the multiple populations in globular clusters

Authors: Mark Gieles, Corinne Charbonnel

Abstract: Globular clusters (GCs) display anomalous light element abundances (HeCNONaMgAl), resembling the yields of hot-hydrogen burning, but there is no consensus yet on the origin of these ubiquitous multiple populations. We present a model in which a super-massive star (SMS, >10^3 Msun) forms via stellar collisions during GC formation and pollutes the intra-cluster medium. The growth of the SMS finds a… ▽ More Globular clusters (GCs) display anomalous light element abundances (HeCNONaMgAl), resembling the yields of hot-hydrogen burning, but there is no consensus yet on the origin of these ubiquitous multiple populations. We present a model in which a super-massive star (SMS, >10^3 Msun) forms via stellar collisions during GC formation and pollutes the intra-cluster medium. The growth of the SMS finds a balance with the wind mass loss rate, such that the SMS can produce a significant fraction of the total GC mass in processed material, thereby overcoming the so-called mass-budget problem that plagues other models. Because of continuous rejuvenation, the SMS acts as a `conveyer-belt' of hot-hydrogen burning yields with (relatively) low He abundances, in agreement with empirical constraints. Additionally, the amount of processed material per unit of GC mass correlates with GC mass, addressing the specific mass budget problem. We discuss uncertainties and tests of this new self-enrichment scenario. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 4 pages, to appear in proceedings of IAUS351 "Star Clusters: from the Milky Way to the Early Universe"

Modelling the Effects of Dark Matter Substructure on Globular Cluster Evolution with the Tidal Approximation

Authors: Jeremy J. Webb, Jo Bovy, Raymond G. Carlberg, Mark Gieles

Abstract: We present direct $N$-body simulations of tidally filling 30,000 ${\rm M}_\odot$ star clusters orbiting between 10 kpc and 100 kpc in galaxies with a range of dark matter substructure properties. The time-dependent tidal force is determined based on the combined tidal tensor of the galaxy's smooth and clumpy dark matter components, the latter of which causes fluctuations in the tidal field that ca… ▽ More We present direct $N$-body simulations of tidally filling 30,000 ${\rm M}_\odot$ star clusters orbiting between 10 kpc and 100 kpc in galaxies with a range of dark matter substructure properties. The time-dependent tidal force is determined based on the combined tidal tensor of the galaxy's smooth and clumpy dark matter components, the latter of which causes fluctuations in the tidal field that can heat clusters. The strength and duration of these fluctuations are sensitive to the local dark matter density, substructure fraction, sub-halo mass function, and the sub-halo mass-size relation. Based on the cold dark matter framework, we initially assume sub-halos are Hernquist spheres following a power-law mass function between $10^5$ and $10^{11} {\rm M}_\odot$ and find that tidal fluctuations are too weak and too short to affect star cluster evolution. Treating sub-halos as point masses, to explore how denser sub-halos affect clusters, we find that only sub-halos with masses greater than $10^{6} {\rm M}_\odot$ will cause cluster dissolution times to decrease. These interactions can also decrease the size of a cluster while increasing the velocity dispersion and tangential anisotropy in the outer regions via tidal heating. Hence increased fluctuations in the tidal tensor, especially fluctuations that are due to low-mass halos, do not necessarily translate into mass loss. We further conclude that the tidal approximation can be used to model cluster evolution in the tidal fields of cosmological simulations with a minimum cold dark matter sub-halo mass of $10^{6} {\rm M}_\odot$, as the effect of lower-mass sub-halos on star clusters is negligible. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 17 pages, 10 figures, Accepted for publication in MNRAS

Population synthesis of black hole binary mergers from star clusters

Authors: Fabio Antonini, Mark Gieles

Abstract: Black hole (BH) binary mergers formed through dynamical interactions in dense star clusters are believed to be one of the main sources of gravitational waves for Advanced LIGO and Virgo. Here we present a fast numerical method for simulating the evolution of star clusters with BHs, including a model for the dynamical formation and merger of BH binaries. Our method is based on Hénon's principle of… ▽ More Black hole (BH) binary mergers formed through dynamical interactions in dense star clusters are believed to be one of the main sources of gravitational waves for Advanced LIGO and Virgo. Here we present a fast numerical method for simulating the evolution of star clusters with BHs, including a model for the dynamical formation and merger of BH binaries. Our method is based on Hénon's principle of balanced evolution, according to which the flow of energy within a cluster must be balanced by the energy production inside its core. Because the heat production in the core is powered by the BHs, one can then link the evolution of the cluster to the evolution of its BH population. This allows us to construct evolutionary tracks of the cluster properties including its BH population and its effect on the cluster and, at the same time, determine the merger rate of BH binaries as well as their eccentricity distributions. The model is publicly available and includes the effects of a BH mass spectrum, mass-loss due to stellar evolution, the ejection of BHs due to natal and dynamical kicks, and relativistic corrections during binary-single encounters. We validate our method using direct $N$-body simulations, and find it to be in excellent agreement with results from recent Monte Carlo models of globular clusters. This establishes our new method as a robust tool for the study of BH dynamics in star clusters and the modelling of gravitational wave sources produced in these systems. Finally, we compute the rate and eccentricity distributions of merging BH binaries for a wide range of cluster initial conditions, spanning more than two orders of magnitude in mass and radius. △ Less

Submitted 1 April, 2020; v1 submitted 27 June, 2019; originally announced June 2019.

Comments: Accepted for publication in MNRAS

Globular clusters as probes of dark matter cusp-core transformations

Authors: M. D. A. Orkney, J. I. Read, James A. Petts, Mark Gieles

Abstract: Bursty star formation in dwarf galaxies can slowly transform a steep dark matter cusp into a constant density core. We explore the possibility that globular clusters (GCs) retain a dynamical memory of this transformation. To test this, we use the nbody6df code to simulate the dynamical evolution of GCs, including stellar evolution, orbiting in static and time-varying potentials for a Hubble time.… ▽ More Bursty star formation in dwarf galaxies can slowly transform a steep dark matter cusp into a constant density core. We explore the possibility that globular clusters (GCs) retain a dynamical memory of this transformation. To test this, we use the nbody6df code to simulate the dynamical evolution of GCs, including stellar evolution, orbiting in static and time-varying potentials for a Hubble time. We find that GCs orbiting within a cored dark matter halo, or within a halo that has undergone a cusp-core transformation, grow to a size that is substantially larger ($R_{\rm eff} > 10$ pc) than those in a static cusped dark matter halo. They also produce much less tidal debris. We find that the cleanest signal of an historic cusp-core transformation is the presence of large GCs with tidal debris. However, the effect is small and will be challenging to observe in real galaxies. Finally, we qualitatively compare our simulated GCs with the observed GC populations in the Fornax, NGC 6822, IKN and Sagittarius dwarf galaxies. We find that the GCs in these dwarf galaxies are systematically larger ($\langle R_{\rm eff}\rangle \simeq 7.8$ pc), and have substantially more scatter in their sizes, than in-situ metal rich GCs in the Milky Way and young massive star clusters forming in M83 ($\langle R_{\rm eff} \rangle \simeq 2.5$ pc). We show that the size, scatter and survival of GCs in dwarf galaxies are all consistent with them having evolved in a constant density core, or a potential that has undergone a cusp-core transformation, but not in a dark matter cusp. △ Less

Submitted 10 September, 2019; v1 submitted 11 June, 2019; originally announced June 2019.

Comments: 14 pages, 10 figures

Journal ref: MNRAS 488 (2019) 2977-2988

The Detailed Science Case for the Maunakea Spectroscopic Explorer, 2019 edition

Authors: The MSE Science Team, Carine Babusiaux, Maria Bergemann, Adam Burgasser, Sara Ellison, Daryl Haggard, Daniel Huber, Manoj Kaplinghat, Ting Li, Jennifer Marshall, Sarah Martell, Alan McConnachie, Will Percival, Aaron Robotham, Yue Shen, Sivarani Thirupathi, Kim-Vy Tran, Christophe Yeche, David Yong, Vardan Adibekyan, Victor Silva Aguirre, George Angelou, Martin Asplund, Michael Balogh, Projjwal Banerjee , et al. (239 additional authors not shown)

Abstract: (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the sc… ▽ More (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more. △ Less

Submitted 9 April, 2019; originally announced April 2019.

Comments: 9 chapters, 301 pages, 100 figures. This version of the DSC is a comprehensive update of the original version, released in 2016, which can be downloaded at arXiv:1606.00043. A detailed summary of the design of MSE is available in the MSE Book 2018, available at arXiv:1810.08695

Spherical models of star clusters with potential escapers

Authors: Ian Claydon, Mark Gieles, Anna Lisa Varri, Douglas C. Heggie, Alice Zocchi

Abstract: An increasing number of observations of the outer regions of globular clusters (GCs) have shown a flattening of the velocity dispersion profile and an extended surface density profile. Formation scenarios of GCs can lead to different explanations of these peculiarities, therefore the dynamics of stars in the outskirts of GCs are an important tool in tracing back the evolutionary history and format… ▽ More An increasing number of observations of the outer regions of globular clusters (GCs) have shown a flattening of the velocity dispersion profile and an extended surface density profile. Formation scenarios of GCs can lead to different explanations of these peculiarities, therefore the dynamics of stars in the outskirts of GCs are an important tool in tracing back the evolutionary history and formation of star clusters. One possible explanation for these features is that GCs are embedded in dark matter halos. Alternatively, these features are the result of a population of energetically unbound stars that can be spatially trapped within the cluster, known as potential escapers (PEs). We present a prescription for the contribution of these energetically unbound members to a family of self-consistent, distribution function-based models, which, for brevity, we call the Spherical Potential Escapers Stitched (SPES) models. We show that, when fitting to mock data of bound and unbound stars from an N-body model of a tidally-limited star cluster, the SPES models correctly reproduce the density and velocity dispersion profiles up to the Jacobi radius, and they are able to recover the value of the Jacobi radius itself to within 20%. We also provide a comparison to the number density and velocity dispersion profiles of the Galactic cluster 47 Tucanae. Such a case offers a proof of concept that an appropriate modelling of PEs is essential to accurately interpret Gaia data in the outskirts of GCs, and, in turn, to formulate meaningful present-day constraints for GC formation scenarios in the early universe. △ Less

Submitted 15 March, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 16 pages, 11 figures, submitted to MNRAS

Stellar Astrophysics and Exoplanet Science with the Maunakea Spectroscopic Explorer (MSE)

Authors: Maria Bergemann, Daniel Huber, Vardan Adibekyan, George Angelou, Daniela Barría, Timothy C. Beers, Paul G. Beck, Earl P. Bellinger, Joachim M. Bestenlehner, Bertram Bitsch, Adam Burgasser, Derek Buzasi, Santi Cassisi, Márcio Catelan, Ana Escorza, Scott W. Fleming, Boris T. Gänsicke, Davide Gandolfi, Rafael A. García, Mark Gieles, Amanda Karakas, Yveline Lebreton, Nicolas Lodieu, Carl Melis, Thibault Merle , et al. (48 additional authors not shown)

Abstract: The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fib… ▽ More The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: 31 pages, 11 figures; To appear as a chapter for the Detailed Science Case of the Maunakea Spectroscopic Explorer

The VLT-FLAMES Tarantula Survey: XXXI. Radial velocities and multiplicity constraints of red supergiant stars in 30 Doradus

Authors: L. R. Patrick, D. J. Lennon, N. Britavskiy, C. J. Evans, H. Sana, W. D. Taylor, A. Herrero, L. A. Almeida, J. S. Clark, M. Gieles, N. Langer, F. R. N. Schneider, J. Th. van Loon

Abstract: The incidence of multiplicity in cool, luminous massive stars is relatively unknown compared to their hotter counterparts. Here we present radial velocity (RV) measurements and investigate the multiplicity properties of red supergiants (RSGs) in the 30~Doradus region of the Large Magellanic Cloud. We provide absolute RV measurements for our sample and estimate line-of-sight velocities for the Hodg… ▽ More The incidence of multiplicity in cool, luminous massive stars is relatively unknown compared to their hotter counterparts. Here we present radial velocity (RV) measurements and investigate the multiplicity properties of red supergiants (RSGs) in the 30~Doradus region of the Large Magellanic Cloud. We provide absolute RV measurements for our sample and estimate line-of-sight velocities for the Hodge 301 and SL 639 clusters, which agree well with those of hot stars in the same clusters. By combining results for the RSGs with those for nearby B-type stars, we estimate systemic velocities and velocity dispersions for the two clusters, obtaining estimates for their dynamical masses of $\log (M_{\rm dyn}/M_{\odot})=$3.8$\pm$0.3 for Hodge 301, and an upper limit of $\log (M_{\rm dyn}/M_{\odot})$<3.1$\pm$0.8 for SL 639, assuming Virial equilibrium. Analysis of the multi-epoch data reveals one RV-variable, potential binary candidate (VFTS744), which is likely a semi-regular variable asymptotic giant branch star. We estimate an upper limit on the observed binary fraction for our sample of 0.3, where we are sensitive to maximum periods for individual objects in the range of 1 to 10 000 days and mass-ratios above 0.3 depending on the data quality. From simulations of the RV measurements from binary systems given the current data we conclude that systems within the parameter range q>0.3, $\log$P[days]<3.5, would be detected by our variability criteria, at the 90% confidence level. The intrinsic binary fraction, accounting for observational biases, is estimated using simulations of binary systems with an empirically defined distribution of parameters where orbital periods are uniformly distributed in the 3.3<$\log$P[days]<4.3 range. A range of intrinsic binary fractions are considered; a binary fraction of 0.3 is found to best reproduce the observed data. [Abridged] △ Less

Submitted 1 March, 2019; originally announced March 2019.

Comments: accepted in A&A

Journal ref: A&A 624, A129 (2019)

Globular cluster number density profiles using Gaia DR2

Authors: T. J. L. de Boer, M. Gieles, E. Balbinot, V. Henault-Brunet, A. Sollima, L. L. Watkins, I. Claydon

Abstract: Using data from Gaia DR2, we study the radial number density profiles of the Galactic globular cluster sample. Proper motions are used for accurate membership selection, especially crucial in the cluster outskirts. Due to the severe crowding in the centres, the Gaia data is supplemented by literature data from HST and surface brightness measurements, where available. This results in 81 clusters wi… ▽ More Using data from Gaia DR2, we study the radial number density profiles of the Galactic globular cluster sample. Proper motions are used for accurate membership selection, especially crucial in the cluster outskirts. Due to the severe crowding in the centres, the Gaia data is supplemented by literature data from HST and surface brightness measurements, where available. This results in 81 clusters with a complete density profile covering the full tidal radius (and beyond) for each cluster. We model the density profiles using a set of single-mass models ranging from King and Wilson models to generalised lowered isothermal limepy models and the recently introduced spes models, which allow for the inclusion of potential escapers. We find that both King and Wilson models are too simple to fully reproduce the density profiles, with King (Wilson) models on average underestimating(overestimating) the radial extent of the clusters. The truncation radii derived from the limepy models are similar to estimates for the Jacobi radii based on the cluster masses and their orbits. We show clear correlations between structural and environmental parameters, as a function of Galactocentric radius and integrated luminosity. Notably, the recovered fraction of potential escapers correlates with cluster pericentre radius, luminosity and cluster concentration. The ratio of half mass over Jacobi radius also correlates with both truncation parameter and PE fraction, showing the effect of Roche lobe filling. △ Less

Submitted 5 March, 2019; v1 submitted 23 January, 2019; originally announced January 2019.

Comments: 23 pages, 23 figures. Accepted in MNRAS

Linking the rotation of a cluster to the spins of its stars: The kinematics of NGC6791 and NGC6819 in 3D

Authors: Sebastian Kamann, Nate Bastian, Mark Gieles, Eduardo Balbinot, Vincent Hénault-Brunet

Abstract: The physics governing the formation of star clusters is still not entirely understood. One open question concerns the amount of angular momentum that newly formed clusters possess after emerging from their parent gas clouds. Recent results suggest an alignment of stellar spins and binary orbital spins in star clusters, which support a scenario in which clusters are born with net angular momentum c… ▽ More The physics governing the formation of star clusters is still not entirely understood. One open question concerns the amount of angular momentum that newly formed clusters possess after emerging from their parent gas clouds. Recent results suggest an alignment of stellar spins and binary orbital spins in star clusters, which support a scenario in which clusters are born with net angular momentum cascading down to stellar scales. In this paper, we combine Gaia data and published line of sight velocities to explore if NGC6791 and NGC6819, two of the clusters for which an alignment of stellar spins has been reported, rotate in the same plane as their stars. We find evidence for rotation in NGC6791 using both proper motions and line of sight velocities. Our estimate of the inclination angle is broadly consistent with the mean inclination that has been determined for its stars, but the uncertainties are still substantial. Our results identify NGC6791 as a promising follow-up candidate to investigate the link between cluster and stellar rotation. We find no evidence for rotation in NGC6819. △ Less

Submitted 20 November, 2018; originally announced November 2018.

Comments: 11 pages, 9 figures, accepted for publication in MNRAS

Black hole growth through hierarchical black hole mergers in dense star clusters: implications for gravitational wave detections

Authors: Fabio Antonini, Mark Gieles, Alessia Gualandris

Abstract: In a star cluster with a sufficiently large escape velocity, black holes (BHs) that are produced by BH mergers can be retained, dynamically form new BH binaries, and merge again. This process can repeat several times and lead to significant mass growth. In this paper, we calculate the mass of the largest BH that can be formed through repeated mergers of stellar seed BHs and determine how its value… ▽ More In a star cluster with a sufficiently large escape velocity, black holes (BHs) that are produced by BH mergers can be retained, dynamically form new BH binaries, and merge again. This process can repeat several times and lead to significant mass growth. In this paper, we calculate the mass of the largest BH that can be formed through repeated mergers of stellar seed BHs and determine how its value depends on the physical properties of the host cluster. We adopt an analytical model in which the energy generated by the black hole binaries in the cluster core is assumed to be regulated by the process of two-body relaxation in the bulk of the system. This principle is used to compute the hardening rate of the binaries and to relate this to the time-dependent global properties of the parent cluster. We demonstrate that in clusters with initial escape velocity $\gtrsim 300\rm km\ s^{-1}$ in the core and density $\gtrsim 10^5\ M_\odot\rm pc^{-3}$, repeated mergers lead to the formation of BHs in the mass range $100-10^5 \,M_\odot$, populating any upper mass gap created by pair-instability supernovae. This result is independent of cluster metallicity and the initial BH spin distribution. We show that about $10\%$ of the present-day nuclear star clusters meet these extreme conditions, and estimate that BH binary mergers with total mass $\gtrsim 100\,M _\odot$ should be produced in these systems at a maximum rate $\approx 0.05 \,\rm Gpc^{-3} yr^{-1}$, corresponding to one detectable event every few years with Advanced LIGO/VIRGO at design sensitivity. The contribution of globular clusters is likely to be negligible instead because the first BH merger remnant escapes following the relativistic kick. A possible connection of our results to the formation of massive BH seeds in galaxy nuclei and globular clusters is discussed. △ Less

Submitted 4 June, 2019; v1 submitted 8 November, 2018; originally announced November 2018.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

Mass modelling globular clusters in the Gaia era: a method comparison using mock data from an $N$-body simulation of M4

Authors: Vincent Hénault-Brunet, Mark Gieles, Antonio Sollima, Laura L. Watkins, Alice Zocchi, Ian Claydon, Elena Pancino, Holger Baumgardt

Abstract: As we enter a golden age for studies of internal kinematics and dynamics of Galactic globular clusters (GCs), it is timely to assess the performance of modelling techniques in recovering the mass, mass profile, and other dynamical properties of GCs. Here, we compare different mass-modelling techniques (distribution-function (DF)-based models, Jeans models, and a grid of N-body models) by applying… ▽ More As we enter a golden age for studies of internal kinematics and dynamics of Galactic globular clusters (GCs), it is timely to assess the performance of modelling techniques in recovering the mass, mass profile, and other dynamical properties of GCs. Here, we compare different mass-modelling techniques (distribution-function (DF)-based models, Jeans models, and a grid of N-body models) by applying them to mock observations from a star-by-star N-body simulation of the GC M 4 by Heggie. The mocks mimic existing and anticipated data for GCs: surface brightness or number density profiles, local stellar mass functions, line-of-sight velocities, and Hubble Space Telescope- and Gaia-like proper motions. We discuss the successes and limitations of the methods. We find that multimass DF-based models, Jeans, and N-body models provide more accurate mass profiles compared to single-mass DF-based models. We highlight complications in fitting the kinematics in the outskirts due to energetically unbound stars associated with the cluster ("potential escapers", not captured by truncated DF models nor by N-body models of clusters in isolation), which can be avoided with DF-based models including potential escapers, or with Jeans models. We discuss ways to account for mass segregation. For example, three-component DF-based models with freedom in their mass function are a simple alternative to avoid the biases of single-mass models (which systematically underestimate the total mass, half-mass radius, and central density), while more realistic multimass DF-based models with freedom in the remnant content represent a promising avenue to infer the total mass and the mass function of remnants. △ Less

Submitted 22 November, 2018; v1 submitted 11 October, 2018; originally announced October 2018.

Comments: 27 pages, 12 figures, MNRAS accepted

Extending the Globular Cluster System-Halo Mass Relation to the Lowest Galaxy Masses

Authors: Duncan A. Forbes, Justin I. Read, Mark Gieles, Michelle L. M. Collins

Abstract: High mass galaxies, with halo masses $M_{200} \ge 10^{10} M_{\odot}$, reveal a remarkable near-linear relation between their globular cluster (GC) system mass and their host galaxy halo mass. Extending this relation to the mass range of dwarf galaxies has been problematic due to the difficulty in measuring independent halo masses. Here we derive new halo masses based on stellar and HI gas kinemati… ▽ More High mass galaxies, with halo masses $M_{200} \ge 10^{10} M_{\odot}$, reveal a remarkable near-linear relation between their globular cluster (GC) system mass and their host galaxy halo mass. Extending this relation to the mass range of dwarf galaxies has been problematic due to the difficulty in measuring independent halo masses. Here we derive new halo masses based on stellar and HI gas kinematics for a sample of nearby dwarf galaxies with GC systems. We find that the GC system mass--halo mass relation for galaxies populated by GCs holds from halo masses of $M_{200} \sim 10^{14} M_{\odot}$ down to below $M_{200}$ $\sim 10^9 M_{\odot}$, although there is a substantial increase in scatter towards low masses. In particular, three well-studied ultra diffuse galaxies, with dwarf-like stellar masses, reveal a wide range in their GC-to-halo mass ratios. We compare our GC system--halo mass relation to the recent model of El Badry et al., finding that their fiducial model does not reproduce our data in the low mass regime. This may suggest that GC formation needs to be more efficient than assumed in their model, or it may be due to the onset of stochastic GC occupation in low mass halos. Finally, we briefly discuss the stellar mass-halo mass relation for our low mass galaxies with GCs, and we suggest some nearby dwarf galaxies for which searches for GCs may be fruitful. △ Less

Submitted 20 September, 2018; originally announced September 2018.

Comments: 16 pages, 5 figures, accepted for publication in MNRAS

From Giant H II regions and H II galaxies to globular clusters and compact dwarf ellipticals

Authors: Elena Terlevich, David Fernández-Arenas, Roberto Terlevich, Mark Gieles, Ricardo Chávez, Ana Luisa González-Morán

Abstract: Massive starforming regions like Giant HII Regions (GHIIR) and HII Galaxies (HIIG) are emission line systems ionized by compact young massive star clusters (YMC) with masses ranging from $10^4$M$_\odot$ to $10^8$M$_\odot$. We model the photometric and dynamical evolution over a Hubble time of the massive gravitationally bound systems that populate the tight relation between absolute blue magnitude… ▽ More Massive starforming regions like Giant HII Regions (GHIIR) and HII Galaxies (HIIG) are emission line systems ionized by compact young massive star clusters (YMC) with masses ranging from $10^4$M$_\odot$ to $10^8$M$_\odot$. We model the photometric and dynamical evolution over a Hubble time of the massive gravitationally bound systems that populate the tight relation between absolute blue magnitude and velocity dispersion ($M_{B}-σ$) of GHIIR and HIIG and compare the resulting relation with that one of old stellar systems: globular clusters, elliptical galaxies, bulges of spirals. After 12~Gyr of evolution their position on the $σ$ vs. M$_B$ plane coincides -- depending on the initial mass -- either with the globular clusters for systems with initial mass $M < 10^6$M$_\odot$ or with a continuation of the ellipticals, bulges of spirals and ultracompact dwarfs for YMC with $M >10^6$M$_\odot$. The slope change in the $M_{B}-σ$ and $M_B$-size relations at cluster masses around $10^6$M$_\odot$ is due to the larger impact of the dynamical evolution on the lower mass clusters. We interpret our result as an indication that the YMC that ionize GHIIR and HIIG can evolve to form globular clusters and ultra compact dwarf ellipticals in about 12 Gyr so that present day globular clusters and ultra compact dwarf ellipticals may have formed in conditions similar to those observed in today GHIIR and HIIG. △ Less

Submitted 21 August, 2018; originally announced August 2018.

Comments: 11 pages, 6 figures, accepted for publication in MNRAS