Properties Of Globular Clusters Research Articles

JWST Imaging of the Closest Globular Clusters. II. Discovery of Brown Dwarfs in NGC 6397 and Measurement of Age from the Brown Dwarf Cooling Sequence, Using SANDee—A New Grid of Model Isochrones across the Hydrogen-burning Limit

Globular clusters contain vast repositories of metal-poor stars that represent some of the oldest stellar generations in the Universe. The archaeological footprint of early Galactic evolution may be retained in the measurable properties of globular clusters, such as their ages, mass functions, and chemical abundances. Until recently, all photometric studies of globular clusters were restricted to stellar members. Now, the sensitivity of JWST can extend this analysis to the substellar regime. If detected in sufficient numbers, brown dwarf members can provide tight constraints on the properties of their parent population. We present SANDee—a new grid of stellar models that accurately represent the color–magnitude diagrams of globular clusters across the hydrogen-burning limit at a wide range of metallicities. Using JWST NIRCam photometry and the new models, we identify three brown dwarfs in the globular cluster NGC 6397 with T eff = 1300–1800 K, confirmed by both proper motion and model fitting. We use the observed luminosities of discovered brown dwarfs to obtain the first age estimate of a globular cluster from its substellar cooling sequence: 13.4 ± 3.3 Gyr. We also derive the local mass function of the cluster across the hydrogen-burning limit and find it to be top heavy, suggesting extensive dynamical evolution. We expect that the constraints on both age and mass function of NGC 6397 derived in this work can be greatly improved by a second epoch of NIRCam imaging in the same field.

Galaxy assembly revealed by globular clusters

Many observable properties of globular clusters (GCs) provide valuable insights for unveiling the hierarchical assembly of their host galaxy. For the Milky Way (MW) in particular, GCs from different accreted satellite galaxies show distinct chemical, spatial, kinematic, and age distributions. Here we examine such clustering features for model GC populations in simulated galaxies, which are carefully selected to match various observational constraints of the MW assembly. We evaluate several widely used clustering, dimensionality reduction, and supervised classification methods on these model GCs, using 10 properties that are observable in the MW. We can categorize in-situ and ex-situ formed GCs with about 90% accuracy, based solely on their clustering features in these 10 variables. The methods are also effective in distinguishing the last major merger in MW analogs with similar accuracy. Although challenging, we still find it possible to identify one, and only one, additional smaller satellite. We develop a new technique to classify the progenitors of MW GCs by combining several methods and weighting them by the validated accuracy. According to this technique, about 60% of GCs belong to the in-situ group, 20% are associated with the Gaia-Sausage/Enceladus event, and 10% are associated with the Sagittarius dwarf galaxy. The remaining 10% of GCs cannot be reliably associated with any single accretion event.

Constrain the Dark-matter Distribution of Ultra-diffuse Galaxies with Globular-cluster Mass Segregation: A Case Study with NGC5846-UDG1

The properties of globular clusters (GCs) contain valuable information of their host galaxies and dark-matter halos. In the remarkable example of ultra-diffuse galaxy, NGC5846-UDG1, the GC population exhibits strong radial mass segregation, indicative of dynamical-friction-driven orbital decay, which opens the possibility of using imaging data alone to constrain the dark-matter content of the galaxy. To explore this possibility, we develop a semianalytical model of GC evolution, which starts from the initial mass, structural, and spatial distributions of the GC progenitors, and follows the effects of dynamical friction, tidal evolution, and two-body relaxation. Using Markov Chain Monte Carlo, we forward-model the GCs in a UDG1-like potential to match the observed GC statistics, and to constrain the profile of the host halo and the origin of the GCs. We find that, with the assumptions of zero mass segregation when the star clusters were born, UDG1 is relatively dark-matter-poor compared to what is expected from stellar-to-halo–mass relations, and its halo concentration is lower than the cosmological average, irrespective of having a cuspy or a cored profile. Its GC population has an initial spatial distribution more extended than the smooth stellar distribution. We discuss the results in the context of scaling laws of galaxy–halo connections, and warn against naively using the GC-abundance–halo–mass relation to infer the halo mass of ultra-diffuse galaxies. Our model is generally applicable to GC-rich dwarf galaxies, and is publicly available.

The accretion history of the Milky Way – II. Internal kinematics of globular clusters and of dwarf galaxies

ABSTRACT We study how structural properties of globular clusters and dwarf galaxies are linked to their orbits in the Milky Way halo. From the inner to the outer halo, orbital energy increases and stellar-systems gradually move out of internal equilibrium: in the inner halo, high-surface brightness globular clusters are at pseudo-equilibrium, while further away, low-surface brightness clusters and dwarfs appear more tidally disturbed. Dwarf galaxies are the latest to arrive into the halo as indicated by their large orbital energies and pericentres, and have no time for more than one orbit. Their (gas-rich) progenitors likely lost their gas during their recent arrival in the Galactic halo. If dwarfs are at equilibrium with their dark matter (DM) content, the DM density should anticorrelate with pericentre. However, the transformation of DM dominated dwarfs from gas-rich rotation-supported into gas-poor dispersion-supported systems is unlikely accomplished during a single orbit. We suggest instead that the above anticorrelation is brought by the combination of ram-pressure stripping and of Galactic tidal shocks. Recent gas removal leads to an expansion of their stellar content caused by the associated gravity loss, making them sufficiently fragile to be transformed near pericentre passage. Out of equilibrium dwarfs would explain the observed anticorrelation of kinematics-based DM density with pericentre without invoking DM density itself, questioning its previous estimates. Ram-pressure stripping and tidal shocks may contribute to the dwarf velocity dispersion excess. It predicts the presence of numerous stars in their outskirts and a few young stars in their cores.

Black Holes in Globular Clusters

Black holes are the laboratory for extreme physics research in the universe, and the detection of black holes is one of the most cutting-edge international scientific goals at present. Globular clusters are old, self-gravity-bound object systems in the Milky Way that contain millions of stars with a high stellar densities inside. Hundreds of pulsar objects have been observed in globular clusters, but the existence of black holes in globular clusters remains an unsolved mystery. A globular cluster may have an intermediatemass black hole at its centre, with masses ranging from 1, 000 to 10, 000 times solar masses. At the same time, there may be a large number of stellar black holes in globular clusters, which will affect the evolution of globular clusters and form. some special objects in globular clusters, which can be observed and confirmed by astronomical observation equipment. This paper will focus on the possible observational features of the presence of black holes in globular clusters and their effects on the properties of globular clusters. Black holes are coupled with their host globular clusters through dynamics, which are important sites for the formation of black hole X-ray binaries and gravitational wave sources, and black holes change the dynamical structure and evolutionary processes of their host clusters. Therefore, searching for X-ray sources, detecting gravitational wave signals, and studying the characteristics of globular clusters are important astronomical evidence for detecting black holes in globular clusters.

How the dynamical properties of globular clusters impact their γ-ray and X-ray emission

ABSTRACT The X-ray and γ-ray emission of globular clusters (GCs) is attributed to their large fraction of compact binary systems, especially those with millisecond pulsars (MSPs). We analyse a population of 124 Galactic GCs to investigate how their dynamical properties affect the formation and evolution of compact binary systems and how this can be translated into the clusters’ observed X-ray and γ-ray emission. We use mainly Chandra X-ray Observatory and Fermi Large Area Telescope observations to achieve our goals and start by detecting 39 GCs in γ rays, seven of which are not listed in previous Fermi-LAT catalogs. Additionally, we find that the total number of X-ray sources within a GC and its γ-ray luminosity are linearly correlated with the stellar encounter rate, indicating that compact binary systems are mainly formed via close stellar encounters. We also find an unexpected rise in the number of X-ray sources for GCs with low rates of stellar encounters, suggesting that there is a dynamical threshold where the formation of X-ray sources is dominated by stellar encounters. Furthermore, we use the Heggie-Hills law to find that subsequent stellar encounters in these compact binaries will, on average, make the binaries even harder, with basically no possibility of binary ionization. Finally, we find that all GCs are point-like sources in γ rays, indicating that the MSPs are concentrated in the clusters’ cores, likely due to dynamical friction.

Investigating dynamical properties of globular clusters through a family of lowered isothermal models

ABSTRACT To investigate the dynamical properties of globular clusters, surface brightness and kinematic data were collected and fitted to a family of lowered isothermal models called limepy models. For the 18 globular clusters studied, the amounts of concentration, truncation, and anisotropy were determined. In addition, the cluster mass, half-mass radius, distance, and mass-to-light ratio were also obtained. In general, limepy models can describe these clusters well. Among these 18 clusters, NGC 5139, 6388, and 7078 have been claimed to be candidates for hosting intermediate-mass black holes in the literature. The models could not appropriately fit the central proper-motion velocity dispersion of NGC 5139, nor the slope of the proper-motion velocity dispersion profile of NGC 6388. Thus, more dedicated models with intermediate-mass black holes or a group of stellar-mass black holes at cluster centres may need to be considered. Regarding NGC 7078, our model with some degree of anisotropy can fit the data. Finally, a strong concentration–truncation anticorrelation and a truncation–semimajor-axis correlation were revealed, which could be the observational imprint of the dynamical evolution of globular clusters.

Properties of Globular Clusters in Galaxy Clusters: Sensitivity from the Formation and Evolution of Globular Clusters

We investigate the properties of globular clusters (GCs) in a galaxy cluster, using the particle tagging method with a semianalytical approach in a cosmological context. We assume GCs form from dark matter halo mergers and their metallicity is assigned based on the stellar mass of the host dark matter halos and the formation redshift of GCs. Dynamical evolution and disruption of GCs are considered using semianalytical approaches, controlled by several free parameters. In this paper, we investigate how our results are changed by the choice of free parameters. We compare our fiducial results with representative observations, including the mass ratio between the GC system and its host galaxy, the GC occupancy, the number fraction of blue GCs, and the metallicity gradient with the GC mass. Because we can know the positions of GCs with time, comparison with additional observations is possible, e.g., the median radii of the GC system in individual galaxies, the mean projected density profiles of intracluster GCs, and the metallicity and age gradients of GCs with a clustercentric radius. We also find that the specific mass of the GC system in each galaxy is different with a clustercentric radius.

Properties of globular clusters formed in dark matter mini-halos

We seek to differentiate dynamical and morphological attributes between globular clusters (GCs) that were formed inside their own dark matter (DM) mini-halo from those who were not. We employed high-resolution full N-body simulations on a Graphics Processing Unit (GPU) of the GCs with and without a DM mini-halo, orbiting a Fornax-like dwarf galaxy. For GCs with DM, we observed that this dark extra mass triggers a tidal radius growth that allows the mini-halo to act as a protective shield against tidal stripping, being itself stripped beforehand. We demonstrate that this shielding effect becomes negligible when the tidal radius is smaller than the half-mass radius of the mini-halo. Contrary to previous predictions, we found that the inflation of outer stellar velocity dispersion profiles is expected for GCs with and without a mini-halo, as a result of the host’s tidal field. Moreover, we observed that GCs with a DM mini-halo should have, in general, relatively more radial outer velocity anisotropy profiles throughout all of their orbits, smaller degrees of internal rotation, and as a consequence of the latter, smaller ellipticities for their stellar distribution. Due to dynamical friction, we observed a clear bimodal evolutionary distribution of GCs with and without DM in the integrals of motion space and show that for GCs originally embedded in DM, this method is not reliable for association with previous accretion events. Finally, we provide parametric mass profiles of disrupted DM mini-halos from GCs that are to be used in Jeans modelling and orbital integration studies.

Integrated Photometry of Multiple Stellar Populations in Globular Clusters

Abstract Evidence that multiple populations (MPs) are common properties of globular clusters (GCs) has accumulated over the past decades from clusters in the Milky Way and in its satellites. This finding has revived research into GCs, and suggested that their formation at high redshift must have been a much more complex phenomenon than imagined before. However, most information on MPs is limited to nearby GCs. The main limitation is that most studies of MPs rely on resolved stars, posing a major challenge to the investigation of the MP phenomenon in distant galaxies. Here we search for integrated colors of old GCs that are sensitive to the MP phenomenon. To do this, we exploit integrated magnitudes of simulated GCs with MPs, and multiband Hubble Space Telescope photometry of 56 Galactic GCs, where MPs are widely studied, and characterized as part of the UV Legacy Survey of Galactic GCs. We find that both integrated C F275W,F336W,F438W and m F275W − m F814W colors strongly correlate with the iron abundance of the host GC. To second order, the pseudo two-color diagram built with these integrated colors is sensitive to the MP phenomenon. In particular, once the dependence on cluster metallicity is removed, the color residuals depend on the maximum internal helium variation within GCs and on the fraction of second-generation stars. This diagram, which we define here for Galactic GCs, has the potential to detect and characterize MPs from integrated photometry of old GCs, thus providing the possibility to extend their investigation outside the Local Group.

Joint constraints on the field-cluster mixing fraction, common envelope efficiency, and globular cluster radii from a population of binary hole mergers via deep learning

The recent release of the second Gravitational-Wave Transient Catalog (GWTC-2) has increased significantly the number of known GW events, enabling unprecedented constraints on formation models of compact binaries. One pressing question is to understand the fraction of binaries originating from different formation channels, such as isolated field formation versus dynamical formation in dense stellar clusters. In this paper, we combine the $\texttt{COSMIC}$ binary population synthesis suite and the $\texttt{CMC}$ code for globular cluster evolution to create a mixture model for black hole binary formation under both formation scenarios. For the first time, these code bodies are combined self-consistently, with $\texttt{CMC}$ itself employing $\texttt{COSMIC}$ to track stellar evolution. We then use a deep-learning enhanced hierarchical Bayesian analysis to constrain the mixture fraction $f$ between formation models, while simultaneously constraining the common envelope efficiency $\alpha$ assumed in $\texttt{COSMIC}$ and the initial cluster virial radius $r_v$ assumed in $\texttt{CMC}$. Under specific assumptions about other uncertain aspects of isolated binary and globular cluster evolution, we report the median and $90\%$ confidence interval of three physical parameters $(f,\alpha,r_v)=(0.20^{+0.32}_{-0.18},2.26^{+2.65}_{-1.84},2.71^{+0.83}_{-1.17})$. This simultaneous constraint agrees with observed properties of globular clusters in the Milky Way and is an important first step in the pathway toward learning astrophysics of compact binary formation from GW observations.

Chemical Composition of Globular Clusters of Milky Way Subsystems from Gaia DR2 Data

Results of a study of the kinematic and chemical properties of globular clusters of the Milky Way based on data from the Gaia DR2 catalog and meaurements with the Hubble Space Telescope are presented. A new method for dividing globular clusters into Galatic subsystems based on the elements of their Galactic orbits is proposed. Samples of globular clusters belonging to the bar/bulge, thick disk, and halo of the Milky Way are obtained. The mean metallicities of the globular clusters in various subsystems are calculated. The mean metallicities of globular clusters of the thick disk and halo display statistically significant differences. At the same time, no statistically significant differences are present between the mean metallicities of halo globular clusters moving in the direction of rotation of the Galactic disk and those moving in the retrograde direction. This argues against the suggestion that retrograde and prograde globular clusters have different origins.

Black hole mergers from globular clusters observable by LISA II. Resolved eccentric sources and the gravitational wave background

In paper I of this series we showed that a large percentage of the binary black hole (BBH) mergers that form through dynamical interactions in globular clusters will have significant eccentricity in the ~10^{-3}-10^{-1} Hz LISA band. In this work we quantify the evolution of these highly eccentric binaries through the LISA and LIGO bands, and compute the stochastic gravitational wave background from the merging, eccentric population. We find that the population of BBHs that merge in-between three-body encounters inside their cluster (~50% of all cluster-formed BBH mergers) will have measurable eccentricity for their entire lifetime in the LISA band. The population of BBHs that merge during three-body encounters (~5% of all cluster-formed BBH mergers), will be detectable by LIGO with eccentricities of e~0.1. The gravitational wave background from dynamically assembled BBHs encodes a characteristic bump due to the high initial eccentricities of these systems. The location and amplitude of this bump depends on globular cluster properties.

Binary black hole mergers from globular clusters: the impact of globular cluster properties

The dense environment of globular clusters (GCs) can facilitate the formation of binary black holes (BBHs), some of which can merge with gravitational waves (GW) within the age of the Universe. We have performed a survey of Monte-Carlo simulations following the dynamical evolution of GCs with different masses, sizes and binary fractions and explored the impact of the host GC properties on the formation of BBH mergers. We find that the number of BBH mergers from GCs is determined by the GC's initial mass, size and primordial binary fraction. We identify two groups of BBH mergers: a primordial group whose formation does not depend on cluster's dynamics and a dynamical group whose formation is driven by the cluster's dynamical evolution. We show how the BBH origin affects the BBH mergers' main properties such as the chirp mass and merging time distributions. We provide analytic expressions for the dependence of the number of BBH mergers from individual GCs on the main cluster's structural properties and the time evolution of the merger rates of these BBHs. These expressions provide an essential ingredient for a general framework allowing to estimate the merger rate density. Using the relations found in our study, we find a local merger rate density of 0.18$-$1.8 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for primordial BBH mergers and 0.6$-$18 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for dynamical BBH mergers, depending on the GC mass and size distributions, initial binary fraction and the number density of GCs in the Universe.

The properties of bright globular clusters, ultra-compact dwarfs and dwarf nuclei in the Virgo core: hints on origin of ultra-compact dwarf galaxies (UCDs)

AbstractBased on the data from the Next Generation Virgo cluster Survey (NGVS), we statistically study the photometric properties of globular clusters (GCs), ultra-compact dwarfs (UCDs) and dwarf nuclei in the Virgo core (M87) region. We found an obvious negative color (g - z) gradient in GC system associate with M87, i.e. GCs in the outer regions are bluer. However, such color gradient does not exist in UCD system, neither in dwarf nuclei system around M87. In addition, we found that many UCDs are surrounded by extended, low surface brightness envelopes. The dwarf nuclei and UCDs show different spatial distributions from GCs, with dwarf nuclei and UCDs (especially for the UCDs with visible envelopes) lying at larger distances to the Virgo center. These results support the view that UCDs (at least for a fraction of UCDs) are more tied to dwarf nuclei than to GCs.

How Black Holes Shape Globular Clusters: Modeling NGC 3201

Abstract Numerical simulations have shown that black holes (BHs) can strongly influence the evolution and present-day observational properties of globular clusters (GCs). Using a Monte Carlo code, we construct GC models that match the Milky Way cluster NGC 3201, the first cluster in which a stellar-mass BH was identified through radial velocity measurements. We predict that NGC 3201 contains ≳200 stellar-mass BHs. Furthermore, we explore the dynamical formation of main-sequence–BH binaries and demonstrate that systems similar to the observed BH binary in NGC 3201 are produced naturally. Additionally, our models predict the existence of bright blue straggler–BH binaries that are unique to core-collapsed clusters, which otherwise retain few BHs.

Globular cluster formation with multiple stellar populations: self-enrichment in fractal massive molecular clouds

Internal chemical abundance spreads are one of fundamental properties of globular clusters (GCs) in the Galaxy. In order to understand the origin of such abundance spreads, we numerically investigate GC formation from massive molecular clouds (MCs) with fractal structures using our new hydrodynamical simulations with star formation and feedback effects of supernovae (SNe) and asymptotic giant branch (AGB) stars. We particularly investigate star formation from gas chemically contaminated by SNe and AGB stars within MCs with different initial conditions and environments. The principal results are as follows. GCs with multiple generation of stars can be formed from merging of hierarchical star cluster complexes that are developed from high-density regions of fractal MCs. Feedback effects of SNe and AGB stars can control the formation efficiencies of stars formed from original gas of MCs and from gas ejected from AGB stars. The simulated GCs have radial gradients of helium abundances within the central 3 pc. The original MC masses need to be as large as 10^7 Msun for a canonical initial stellar mass function (IMF) so that the final masses of stars formed from AGB ejecta can be 10^5 Msun. Since star formation from AGB ejecta is rather prolonged (10^8 yr), their formation can be strongly suppressed by SNe of the stars themselves. This result implies that the so-called mass budget problem is much more severe than ever thought in the self-enrichment scenario of GC formation. and thus that IMF for the second generation of stars should be `top-light'.

A novel look at energy equipartition in globular clusters

Two-body interactions play a major role in shaping the structural and dynamical properties of globular clusters (GCs) over their long-term evolution. In particular, GCs evolve toward a state of partial energy equipartition that induces a mass-dependence in their kinematics. By using a set of Monte Carlo cluster simulations evolved in quasi-isolation, we show that the stellar mass dependence of the velocity dispersion $\sigma(m)$ can be described by an exponential function $\sigma^2\propto \exp(-m/m_\mathrm{eq})$, with the parameter $m_\mathrm{eq}$ quantifying the degree of partial energy equipartition of the systems. This simple parametrization successfully captures the behaviour of the velocity dispersion at lower as well as higher stellar masses, that is, the regime where the system is expected to approach full equipartition. We find a tight correlation between the degree of equipartition reached by a GC and its dynamical state, indicating that clusters that are more than about 20 core relaxation times old, have reached a maximum degree of equipartition. This equipartition$-$dynamical state relation can be used as a tool to characterize the relaxation condition of a cluster with a kinematic measure of the $m_\mathrm{eq}$ parameter. Vice versa, the mass-dependence of the kinematics can be predicted knowing the relaxation time solely on the basis of photometric measurements. Moreover, any deviations from this tight relation could be used as a probe of a peculiar dynamical history of a cluster. Finally, our novel approach is important for the interpretation of state-of-the-art Hubble Space Telescope proper motion data, for which the mass dependence of kinematics can now be measured, and for the application of modeling techniques which take into consideration multi-mass components and mass segregation.

Constraining ultracompact dwarf galaxy formation with galaxy clusters in the local universe

We compare the predictions of a semi-analytic model for ultra-compact dwarf galaxy (UCD) formation by tidal stripping to the observed properties of globular clusters (GCs) and UCDs in the Fornax and Virgo clusters. For Fornax we find the predicted number of stripped nuclei agrees very well with the excess number of GCs$+$UCDs above the GC luminosity function. GCs$+$UCDs with masses $>10^{7.3}$ M$_\odot$ are consistent with being entirely formed by tidal stripping. Stripped nuclei can also account for Virgo UCDs with masses $>10^{7.3}$ M$_\odot$ where numbers are complete by mass. For both Fornax and Virgo, the predicted velocity dispersions and radial distributions of stripped nuclei are consistent with that of UCDs within $\sim$50-100 kpc but disagree at larger distances where dispersions are too high and radial distributions too extended. Stripped nuclei are predicted to have radially biased anisotropies at all radii, agreeing with Virgo UCDs at clustercentric distances larger than 50 kpc. However, ongoing disruption is not included in our model which would cause orbits to become tangentially biased at small radii. We find the predicted metallicities and central black hole masses of stripped nuclei agree well with the metallicities and implied black hole masses of UCDs for masses $>10^{6.5}$ M$_\odot$. The predicted black hole masses also agree well with that of M60-UCD1, the first UCD with a confirmed central black hole. These results suggest that observed GC$+$UCD populations are a combination of genuine GCs and stripped nuclei, with the contribution of stripped nuclei increasing toward the high-mass end.

MOCCA code for star cluster simulations – IV. A new scenario for intermediate mass black hole formation in globular clusters

We discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the influence of an intermediate mass black hole on the host global globular cluster properties. We discuss two scenarios for intermediate mass black hole formation. The SLOW and FAST scenarios. They occur later or earlier in the cluster evolution and require smaller or extremely large central densities, respectively. In our simulations, the formation of intermediate mass black holes is highly stochastic. In general, higher formation probabilities follow from larger cluster concentrations (i.e. central densities). We further discuss possible observational signatures of the presence of intermediate mass black holes in globular clusters that follow from our simulations. These include the spatial and kinematic structure of the host cluster, possible radio, X-ray and gravitational wave emissions due to dynamical collisions or mass-transfer and the creation of hypervelocity main sequence escapers during strong dynamical interactions between binaries and an intermediate mass black hole. All simulations discussed in this paper were performed with the MOCCA Monte Carlo code. MOCCA accurately follows most of the important physical processes that occur during the dynamical evolution of star clusters but, as with other dynamical codes, it approximates the dissipative processes connected with stellar collisions and binary mergers.