Star Clusters Research Articles

Resolved properties of a luminous hinge clump in the compact group of galaxies NGC 6845

Context. Compact groups of galaxies are unique places where galaxy-galaxy interactions play a mayor role on the evolution of its members. These strong gravitational encounters can induce star formation bursts. Aims. We study the properties of one of the most luminous so-called hinge clumps that is located in the compact group of galaxies NGC 6845. Methods. Using integral field spectroscopy from GMOS/Gemini complemented with archival MUSE data, we obtained oxygen abundances, ages, star formation rates, and velocity fields, and we also modeled a single stellar population to understand the star formation history of the hinge clump in NGC 6845. Results. We found that the hinge clump sits in a tailthat has a star formation rate of 3.4 M⊙ yr−1, which is comparable with a few other extreme cases, for instance, the star clusters in the Antennae galaxy and other reported hinge clumps in the literature. This clump represents ∼15% of total star formation rate of NGC 6845A. Large-scale modeling of the observed velocity field of NGC 6845A rules out the scenario according to which this hinge clump was a satellite galaxy. Its kinematics is compatible with that of the galactic disk of NGC 6845A. Its abundance with a mean value of 0.4 Z⊙ is also consistent with the metallicity gradient of the galaxy. Conclusions. Our analysis suggest that the hinge clump is formed by multiple stellar populations and not by a single burst. This causes the wide age range. We found that the central clump is encompassed by a ring-like structure, which might suggest a second generation of star formation. In addition, the analysis of the diagnostic diagram indicates that this central region might also be ionized by shocks from stellar and supernova winds. Finally, the derived star formation rate density Σ = 9.7 M⊙ yr−1 kpc−2 of the central clump places it in starburst regime, where gas inflows should provide gas to maintain the star formation. This work shows a resolved example of an extreme localized starburst in a compact group of galaxies.

The 33 M⊙ black hole Gaia BH3 is part of the disrupted ED-2 star cluster

Context. The Gaia Collaboration has recently reported the detection of a 33 M⊙ black hole in a wide binary system located in the solar neighbourhood. Aims. Here we explore the relationship between this black hole, known as Gaia BH3, and the nearby ED-2 halo stellar stream. Methods. We studied the orbital characteristics of the Gaia BH3 binary and present measurements of the chemical abundances of ED-2 member stars derived from high-resolution spectra obtained with the VLT. Results. We find that the Galactic orbit of the Gaia BH3 system and its metallicity are entirely consistent with being part of the ED-2 stream. The characteristics of the stream, particularly its negligible spread in metallicity and in other chemical elements, as well as its single stellar population, suggest that it originated from a disrupted star cluster of low mass. Its age is comparable to that of the globular cluster M92 that has been estimated to be as old as the Universe. Conclusions. This is the first black hole unambiguously associated with a disrupted star cluster. We infer the plausible mass range for the cluster to be relatively narrow, between 2 × 103 M⊙ and 5.2 × 104 M⊙. This implies that the black hole could have formed directly from the collapse of a massive very metal-poor star, but that the alternative scenario of binary interactions inside the cluster environment also deserves to be explored.

Metallicity relations in LMC and SMC from the slope of red giant branch stars in globular clusters

Metallicity relations in LMC and SMC from the slope of red giant branch stars in globular clusters

Bound star clusters observed in a lensed galaxy 460 Myr after the Big Bang

The Cosmic Gems arc is among the brightest and highly magnified galaxies observed at redshift z ≈ 10.2 (ref. 1). However, it is an intrinsically ultraviolet faint galaxy, in the range of those now thought to drive the reionization of the Universe2–4. Hitherto the smallest features resolved in a galaxy at a comparable redshift are between a few hundreds and a few tens of parsecs (pc)5,6. Here we report JWST observations of the Cosmic Gems. The light of the galaxy is resolved into five star clusters located in a region smaller than 70 pc. They exhibit minimal dust attenuation and low metallicity, ages younger than 50 Myr and intrinsic masses of about 106M⊙. Their lensing-corrected sizes are approximately 1 pc, resulting in stellar surface densities near 105M⊙ pc−2, three orders of magnitude higher than typical young star clusters in the local Universe7. Despite the uncertainties inherent to the lensing model, they are consistent with being gravitationally bound stellar systems, that is, proto-globular clusters. We conclude that star cluster formation and feedback likely contributed to shaping the properties of galaxies during the epoch of reionization.

Massive Interacting Binaries as an Enrichment Source for Multiple Populations in Star Clusters

We present a suite of binary evolution models with massive primaries (10 ≤ M 1 ≤ 40 M ⊙) and periods and mass ratios chosen such that the systems undergo nonconservative mass transfer while the primaries have helium cores. We track the total mass and chemical composition of the ejecta from these systems. This material shows the abundance signatures of hot hydrogen burning that are needed to explain the abundance patterns seen in multiple populations in massive star clusters. We then calculate the total yield of a population of binary stars with masses, mass ratios, and periods consistent with their distribution in a field population. We show that the overall abundance of this material is enriched in helium, nitrogen, sodium, and aluminum, and depleted in carbon, oxygen, and magnesium, by amounts that are consistent with observations. We also show that such a population of binaries will return approximately 25% of its mass in this ejecta (compared to 4% if all the stars were single), over a characteristic timescale of about 12 Myr. We argue that massive binaries must be seriously considered as a contributor to the source of enriched material needed to explain the multiple populations in massive clusters, since essentially all massive stars are formed in binaries or higher-order multiples, massive binaries are primarily formed in clusters, and massive binaries naturally produce material of the right composition.

T-ReX: The Tarantula—Revealed by X-Rays

The Tarantula Nebula (30 Doradus, 30 Dor) is the most important star-forming complex in the Local Group, offering a microscope on starburst astrophysics. At its heart lies the exceptionally rich young stellar cluster R136, containing the most massive stars known. Stellar winds and supernovae have carved 30 Dor into an amazing display of arcs, pillars, and bubbles. We present first results and advanced data-processing products from the 2 Ms Chandra X-ray Visionary Project, “The Tarantula—Revealed by X-rays” (T-ReX). The 3615 point sources in the T-ReX catalog include massive stars, compact objects, binaries, bright pre-main-sequence stars, and compact young stellar (sub)clusters in 30 Dor. After removing point sources and excluding the exceptionally bright supernova remnant N157B (30 Dor B), the global diffuse X-ray maps reveal hot plasma structures resolved at 1–10 pc scales, with an absorption-corrected total-band (0.5–7 keV) X-ray luminosity of 2.110 × 1037 erg s−1. Spatially resolved spectral modeling provides evidence for emission lines enhanced by charge-exchange processes at the interfaces. We identify a candidate for the oldest X-ray pulsar detected to date in 30 Dor, PSR J0538-6902, inside a newly resolved arcuate X-ray wind nebula, the Manta Ray. The long temporal baseline of T-ReX allowed monitoring of dozens of massive stars, several showing periodic variability tied to binary orbital periods, and captured strong flares from at least three low-mass Galactic foreground stars.

Stochastic gravitational wave background from binary black hole mergers dynamically assembled in dense star clusters

Stochastic gravitational wave background from binary black hole mergers dynamically assembled in dense star clusters

The VISCACHA survey – X. A fast method to build completeness maps for Magellanic Clouds star clusters

ABSTRACT Photometric completeness affects the photometry of stars in crowded regions such as the cores of star clusters. Some analysis such as deriving the structural parameters of star clusters using radial density profile is heavily affected by photometric completeness and the classical techniques to map this completeness in a given field are very expensive computationally. In most surveys, for example, the large quantity of data makes it impracticable to estimate the completeness using the traditional method for the whole sample due to time and computational requirements. In this work, we present a new method that is significantly faster and results in similar completeness curves and maps as the traditional approach, providing a great first-step completeness estimator for a large sample of data. Using the completeness corrected data for each cluster we built radial density profiles improving significantly the inner portion of the profile; we also fitted the King model to them, determining the clusters’ structural parameters based on a more realistic cluster profile. In this preliminary analysis, we derived structural parameters for nine selected clusters covering a range of core radii (5–40 arcesc) and tidal radii (40–180 arcesc) and discuss how the photometric completeness affects the determination of these parameters when we count stars to trace the radial profile of a star cluster.

Dynamics of star clusters with tangentially anisotropic velocity distribution

Recent high-precision observations with HST and Gaia enabled new investigations of the internal kinematics of star clusters (SCs) and the dependence of kinematic properties on the stellar mass. These studies raised new questions about the dynamical evolution of self-gravitating stellar systems. We aim to develop a more complete theoretical understanding of how the various kinematical properties of stars affect the global dynamical development of their host SCs. We perform $N$-body simulations of globular clusters with isotropic, radially anisotropic, and tangentially anisotropic initial velocity distributions. We also study the effect of an external Galactic tidal field. We obtain three main results. First, compared to the conventional, isotropic case, the relaxation processes are accelerated in the tangentially anisotropic models and, in agreement with our previous investigations, are slower in the radially anisotropic ones. This leads to, for example, more rapid mass segregation in the central regions of the tangential models or their earlier core collapse. Second, although all SCs become isotropic in the inner regions after several relaxation times, we observe differences in the anisotropy profile evolution in the outer cluster regions --- all tidally filling models gain tangential anisotropy there, while the underfilling models become radially anisotropic. Third, we observe different rates of evolution towards energy equipartition (EEP). While all SCs evolve towards EEP in their inner regions (regardless of the filling factor), the outer regions of the tangentially anisotropic and isotropic models are evolving to an `inverted' EEP (i.e.\ with the high-mass stars having higher velocity dispersion than the low-mass ones). The extent (both spatial and temporal) of this inversion can be attributed to the initial velocity anisotropy --- it grows with increasing tangential anisotropy and decreases as the radial anisotropy rises.

Evolution of Stellar Rotation in Open Clusters

Understanding stellar rotation and the role of magnetic braking remains a major outstanding problem in astrophysics. In this paper, stellar rotation of the young open cluster Blanco 1 is investigated and compared with other clusters with ages ranging from ∼35 Myr to ∼950 Myr and star masses in the range 0.2 < M⊙ < 1.4. It is proposed that rotation rates of stars in young open clusters are determined by the early angular momentum acquired during formation and not by magnetic braking effects, which operate on longer time scales. On the convective C Sequence, for lighter stars, the early angular velocity is related directly to the moment of inertia. However, for the interface or I Sequence, differential rotation exists and observed rotation rates are essentially the rotation rates of the outer convection zone. It is estimated that the convection zone can rotate at over 10× the rate of the inner radiative zone for heavier stars. This differential rotation process can drive saturated interface dynamos that lead to large rotational mass dependences which can be described by a simple energy transfer model. Comparing with clusters of different ages magnetic braking outlines the evolution of older clusters from young clusters. For old clusters such as Praesepe the majority of stars have experienced spin down, that can be described by a simple model based on an additional momentum loss proportional to the angular momentum of the early cluster, which for the I Sequence is consistent with the well known Skumanich relation.

Formation of the inner Oort cloud in the presence of an early solar binary

ABSTRACT Consistent with the notion that most Sun-like stars form in multistellar systems, this study explores the impact of a temporarily bound stellar binary companion on the early dynamical evolution of the Solar system. Using N-body simulations, we illustrate how such a companion markedly enhances the trapping of scattered bodies on inner Oort cloud-like orbits, with perihelion distances exceeding $q \gt 40$ au. We further find that the orbital geometry of the Sun-binary system plays a central role in regulating the efficiency of small-body implantation on to high-perihelion orbits, and demonstrate that this process is driven by the von Zeipel–Kozai–Lidov mechanism. Incorporating the transiency of stellar clusters and the eventual Sun-binary pair dissociation due to passing stars, we show how the binary can be stripped away by an approximately solar-mass ejector star, with only a modest impact on the generated inner Oort cloud population. Collectively, our results highlight a previously underappreciated process that could have contributed to the formation of the inner Oort cloud.

Detection of open cluster members inside and beyond tidal radius by machine learning methods based on Gaia DR3

ABSTRACT In our previous work, we introduced a method that combines two unsupervised algorithms: Density-based spatial clustering of applications with noise (DBSCAN) and Gaussian mixture model (GMM). We applied this method to 12 open clusters based on Gaia Early Data Release 3 (EDR3) data, demonstrating its effectiveness in identifying reliable cluster members within the tidal radius. However, for studying cluster morphology, we need a method capable of detecting members both inside and outside the tidal radius. By incorporating a supervised algorithm into our approach, we successfully identified members beyond the tidal radius. In our current work, we initially applied DBSCAN and GMM to identify reliable members of cluster stars. Subsequently, we trained the random forest algorithm using DBSCAN- and GMM-selected data. Leveraging the random forest, we can identify cluster members outside the tidal radius and observe cluster morphology across a wide field of view. Our method was then applied to 15 open clusters based on Gaia DR3, which exhibit a wide range of metallicity, distances, members, and ages. Additionally, we calculated the tidal radius for each of the 15 clusters using the King profile and detected stars both inside and outside this radius. Finally, we investigated mass segregation and luminosity distribution within the clusters. Overall, our approach significantly improved the estimation of the tidal radius and detection of mass segregation compared to the previous work. We found that in Collinder 463, low-mass stars do not segregate in comparison to high-mass and intermediate-mass stars. Additionally, we detected a peak of luminosity in the clusters, some of which were located far from the centre, beyond the tidal radius.

Orbital evolution of LIGO/Virgo binaries in stellar clusters driven by cluster tides, stellar encounters, and general relativity

ABSTRACT Origin of LIGO/Virgo gravitational wave events may involve production of binaries with relativistic components in dense stellar systems – globular or nuclear star clusters – and their subsequent evolution towards merger. Orbital parameters of these binaries (the inner orbit) and their motion inside the cluster (the outer orbit) evolve due to both external agents – random encounters with cluster stars and cluster tides due to the smooth cluster potential – and the internal ones – various sources of dissipation and precession within the binary. We present a numerical framework – binary evolution in stellar clusters (besc) – that follows the evolution of the binary inner and outer orbits accounting for all these effects simultaneously, enabling efficient Monte Carlo studies. The secular effect of cluster tides is computed in the singly averaged approximation, without averaging over the outer binary orbit. As to stellar encounters, we include the effects of both close and distant flybys on the inner and outer orbits of the binary, respectively. Also, given our focus on the LIGO/Virgo sources, we include the general relativistic precession (which suppresses cluster tides at high eccentricities) and the gravitational wave emission (shrinking the binary orbit). We use besc to illustrate a number of characteristic binary evolutionary outcomes and discuss relative contributions of different physical processes. In particular, we find that stellar encounters often dominate over the cluster tides in the evolution of compact binaries. besc can also be used to study other objects in clusters, e.g. blue stragglers, hot Jupiters, X-ray binaries, etc.

Most nearby young star clusters formed in three massive complexes.

Efforts to unveil the structure of the local interstellar medium and its recent star-formation history have spanned the past 70 years (refs. 1-6). Recent studies using precise data from space astrometry missions have revealed nearby, newly formed star clusters with connected origins7-12. Nonetheless, mapping young clusters across the entire sky back to their natal regions has been hindered by a lack of clusters with precise radial-velocity data. Here we show that 155 out of 272 (57%) high-quality young clusters13,14 within 1 kiloparsec of the Sun arise from three distinct spatial volumes. This conclusion is based on the analysis of data from the third Gaia release15 and other large-scale spectroscopic surveys. At present, dispersed throughout the solar neighbourhood, their past positions more than 30 million yearsago reveal that these families of clusters each formed in one of three compact, massive star-forming complexes. One of these families includes all of the young clusters near the Sun-the Taurus and Scorpius-Centaurus star-forming complexes16,17. We estimate that more than 200 supernovae were produced from these families and argue that these clustered supernovae produced both the Local Bubble18 and the largest nearby supershell GSH 238+00+09 (ref. 19), both of which are clearly visible in modern three-dimensional dust maps20-22.

A baseline on the relation between chemical patterns and the birth stellar cluster

A baseline on the relation between chemical patterns and the birth stellar cluster

CO in the Draco nebula: the atomic–molecular transition

ABSTRACT This paper presents maps of the J=2-1 transition of CO towards the Draco nebula intermediate-velocity cloud (IVC). The maps cover 8500 square arcmin with a velocity resolution of 0.33 km s$^{-1}$ and angular resolution of 38 ${\rm arcsec}$, or 0.11 pc at the cloud distance of 600 pc. The mapped area includes all the emission detected by the Herschel satellite with 250 $\mu$m intensity >5 MJy sr−1. Previously published observations of the far-infrared (far-IR) emission and the 21 cm line of HI are used to derive the column density distribution of H$_2$ and the abundance ratio CO/H$_2$, as well as the distribution of the molecular fraction of hydrogen, which approaches 90 per cent over much of the brighter parts of the nebula. The CO emission is highly clumpy and closely resembles the structures seen in far-IR images. The kinematics of the CO show supersonic motions between clumps but near-thermal to trans-sonic motions within clumps, consistent with model predictions that the scale length for dissipation of supersonic turbulence should be $\sim 0.1$ pc, mediated by kinematic viscosity and/or ambipolar diffusion. Different parts of the nebula show evidence for a spread of molecular formation time-scales of a few 10$^5$ yr, comparable with the dynamical time-scale of the infalling gas. The IVC will likely merge with the Galactic interstellar medium in $\sim 10^7$ yr, and the densest clumps may form an unbound cluster of low-mass stars.

Rescuer: cosmological K-corrections for star clusters

ABSTRACT The advent of JWST now allows entire star cluster populations to be imaged in galaxies at cosmologically significant redshifts, bringing with it the need to apply K-corrections to their magnitudes and colour indices. Since the stellar populations within star clusters can be well approximated by a single age and metallicity, their spectral energy distributions are very different from those of galaxies or supernovae, and their K-corrections behave differently. We derive the photometric K-corrections versus redshift for model star clusters that cover a wide range of ages and metallicities, illustrating the results particularly for the broad-band filters on the Hubble Space Telescope/Advanced Camera for Surveys and the JWST/Near-Infrared Camera that are most commonly being used for imaging of populations of star clusters in distant galaxies. In an appendix, we introduce a webtool called rescuer (REdshifted Star ClUstERs) that can generate K-values and their uncertainties for any user-defined combination of cluster properties.

Deriving physical parameters of unresolved star clusters. VIII. Limits of aperture photometry for star cluster studies

Recently, it has been noticed that the discrepancies in the integrated colour indices (CIs) between star clusters and models are mostly due to the projection of bright stars in the apertures. In order to reduce this problem, the method of adaptive aperture photometry has been proposed. This method has been applied to star clusters from the M\,31 Panchromatic $Hubble$ Andromeda Treasury (PHAT) survey, and studies show that the adaptive aperture photometry performs better than the conventional approach. The aim of this study is to determine the best achievable limits on the accuracy and applicability of the aperture photometry method for studying star clusters in the local Universe. We computed a large network of artificial 3D star clusters spanning the parameter space of the M\,31 clusters. We then simulated images of these clusters by projecting each onto a 2D plane from 100 directions. Star cluster images were generated in six passbands to match the PHAT survey. To investigate the limiting accuracy of aperture photometry and the limits of its applicability to star cluster studies, we measured the simulated images and performed parameter determination tests. We demonstrate that star clusters with and without post-main-sequence stars have significant photometric differences. We show that in order to obtain reliable physical parameters of star clusters, the CIs must be measured using an aperture with a radius larger than the cluster’s half-light radius. Furthermore, we demonstrate that the parameter determination of young clusters (sim 10\,Myr) is problematic regardless of the aperture size used. Therefore, it is advisable to determine the parameters of these clusters using colour-magnitude diagram fitting methods, when possible. We also show that the randomness of the viewing angle can lead to a CI uncertainty of up to 0.1\,mag, depending on cluster parameters and aperture size.

Great Balls of FIRE

Context. Despite the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Aims. Here we focus on understanding the dynamical formation of merging BBHs within massive clusters in galaxies of different masses. Methods. To this end, we applied a new framework to consistently model the formation and evolution of massive star clusters in zoom-in cosmological simulations of galaxies. Each simulation, taken from the FIRE project, provides a realistic star formation environment, with a unique star formation history, that hosts realistic giant molecular clouds that constitute the birthplace of star clusters. Combined with the code for star cluster evolution CMC, we are able to produce populations of dynamically formed merging BBHs across cosmic time in different environments. Results. As the most massive star clusters preferentially form in dense massive clouds of gas, we find that, despite their low metallicities favouring the creation of black holes, low-mass galaxies contain few massive clusters and therefore make a limited contribution to the global production of dynamically formed merging BBHs. Furthermore, we find that massive clusters can host hierarchical BBH mergers with clear, identifiable physical properties. Looking at the evolution of the BBH merger rate in different galaxies, we find strong correlations between BBH mergers and the most extreme episodes of star formation. Finally, we discuss the implications for future LIGO-Virgo-KAGRA gravitational wave observations.

A Comprehensive Kinematic Model of the Large Magellanic Cloud Disk from Star Clusters and Field Stars using Gaia DR3: Tracing the Disk Characteristics, Rotation, Bar, and Outliers

The internal kinematics of the Large Magellanic Cloud (LMC) disk have been modeled by several studies using different tracers with varying coverage, resulting in a range of parameters. Here, we model the LMC disk using 1705 star clusters and field stars, based on a robust Markov Chain Monte Carlo method, using Gaia DR3 data. The dependency of the model parameters on the age, coverage, and strength of the clusters are also presented. This is the first comprehensive 2D kinematic study using star clusters. Red clump (RC) stars and young main-sequence stars are also modeled for comparison. The clusters and field stars are found to have distinctly different kinematic centers, disk inclination, position angle of the line of nodes, and scale radius. We also note a significant radial variation of the disk parameters. Clusters and young stars are found to have a large residual proper motion and a relatively large velocity dispersion when compared to the RC field population, which could be due to perturbation from the bar and spiral arms. We trace the presence of the large residual proper motion and noncircular motion among clusters likely to be due to the bar and detect a decrease in the scale radius as a result of the possible evolution of the bar. The kinematically deviant clusters point to a spatiotemporal disturbance in the LMC disk, matching with the expected impact factor and time of the recent collision between the LMC and the Small Magellanic Cloud.