Globular Research Articles

Gravitational wave background from primordial black holes in globular clusters

Primordial black holes still represent a viable candidate for a significant fraction, if not for the totality, of dark matter. If these compact objects have masses of order tens of solar masses, their coalescence can be observed by current and future ground-based gravitational wave detectors. Therefore, finding new gravitational wave signatures associated with this dark matter candidate can either lead to their detection or help constraining their abundance. In this work we consider the phenomenology of primordial black holes in dense environments, in particular globular clusters. We model the internal structure of globular clusters in a semi-analytical fashion, and we derive the expected merger rate. We show that, if primordial black holes are present in globular clusters, their contribution to the GW background can be comparable to other well-known channels, such as early- and late-time binaries, thus enhancing the detectability prospects of primordial black holes and demonstrating that this contribution needs to be taken into account.

Testing the sources of the peculiar abundances in globular clusters

This work aims to analyze some of the polluters proposed in the self-enrichment scenarios put forward to explain the multiple populations in globular clusters (GCs), extending previous studies. Three scenarios with different polluter stars were tested: asymptotic giant branch stars (AGBs), high-mass interacting binaries (IBs), and fast rotating massive stars (FRMSs). With abundance data available from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey and ΔY estimates from precise Hubble Space Telescope (HST) photometry, twenty-six clusters were studied, increasing the number of clusters in previous studies by more than a factor of three. We also included the study of the abundances of N, C, Mg, and Al, extending previous studies that mainly focused on the abundances of He, O, and Na. In addition, we constructed an empirical model to test whether one could explain the chemical signatures of the “enriched” population of GC stars with a fixed source and dilution process based on empirical data. In agreement with work by other authors, we found that the proposed polluters can generally predict the qualitative abundance patterns in GC stars and in some cases quantitatively predict some elements, but in most cases when we compare the model yields with the observations, we find that they cannot explain the entire set of observed abundance patterns. The empirical model succeeds in reproducing the abundances of Al for a given ΔY (and vice versa), showing that there is a direct relationship between Al and He, with one increasing proportionally to the other. However, the empirical model fails to reproduce the observed abundances of Na and N, in agreement with the results of previous works. The observed decoupling between the maximum abundances of CNO-cycle elements such as N and Na with those of Al and He provides new information and constraints for future models and could take us a step closer to understanding the origin of the peculiar abundance variations of GC stars.

Extragalactic globular cluster near-infrared spectroscopy

Context. One way to constrain the evolutionary histories of galaxies is to analyse their stellar populations. In the local Universe, our understanding of the stellar population properties of galaxies has traditionally relied on the study of optical absorption and emission-line features. Aims. In order to overcome limitations intrinsic to this wavelength range, such as the age-metallicity degeneracy and the high sensitivity to dust reddening, we must use wavelength ranges beyond the optical. The near-infrared (NIR) offers a possibility to extract information on spectral signatures that are not as obvious in traditional optical bands. Moreover, with the current and forthcoming generation of instrumentation focusing on the NIR, it is mandatory to explore possibilities within this wavelength range for nearby-Universe galaxies. However, although the NIR shows great potential, we are only beginning to understand it. Widely used techniques such as a full spectral fitting and line strength indices need to be tested on systems that are as close to simple stellar populations as possible, and the result from the techniques need to be compared to the yields from a traditional optical analysis. Methods. We present a NIR spectral survey of extragalactic globular clusters (GCs). The set was composed of 21 GCs from the Centaurus A galaxy that were obtained with SOAR/TripleSpec4, which covered the ∼1.0–2.4 μm range with a spectral resolution (R = λ/Δλ) of 3500. These spectra cover Hβ equivalent widths between 0.98 Å and 4.32 Å, and [MgFe]′ between 0.24 Å and 3.76 Å. Results. This set was ideal for performing absorption band measurements and a full spectral fitting, and it can be used for kinematic studies and age and abundance measurements. With this library, we expect to be able to probe the capabilities of NIR models, as well as to further improve stellar population estimates for the GCs around the Centaurus A galaxy.

Discovery of a Millisecond Pulsar Associated with Terzan 6

Observations show that globular clusters (GCs) might be among the best places to find millisecond pulsars (MSPs). However, the GC Terzan 6 seems to be an exception without any pulsar discovered, although its high stellar encounter rate suggests that it harbors dozens of them. We report the discovery of the first radio pulsar, PSR J1750–3116A, likely associated with Terzan 6 in a search of C-band (4–8 GHz) data from the Robert C. Byrd Green Bank Telescope with a spin period of 5.33 ms and dispersion measure (DM) ≃ 383 pc cm−3. The mean flux density of this pulsar is approximately 3 μJy. The DM agrees well with predictions from the Galactic free electron density model, assuming a distance of 6.7 kpc for Terzan 6. PSR J1750–3116A is likely an isolated MSP, potentially formed through dynamical interactions, considering the core-collapsed classification and the exceptionally high stellar encounter rate of Terzan 6. This is the highest radio frequency observation that has led to the discovery of a pulsar in a GC to date. While L-band (1–2 GHz) observations of this cluster are unlikely to yield significant returns due to propagation effects, we predict that further pulsar discoveries in Terzan 6 will be made by existing radio telescopes at higher frequencies.

Impact of charge on the stability of pulsar SAX J1748.9-2021 in modified symmetric teleparallel gravity

In this paper, we explore the effect of charge on the stability of pulsar star SAX J1748.9-2021 in f(Q) gravity, where Q represents non-metricity. For this purpose, we apply the Krori-Barua metric ansatz with anisotropic fluid and use a linear f(Q) model f(Q)=ζQ, where ζ is a non-zero constant. We derive exact relativistic solutions of the corresponding field equations. Furthermore, we study its geometric and physical properties through astrophysical observations from the pulsar SAX J1748.9-2021, which is found in X-ray binary systems within globular clusters. We examine features like anisotropic pressure, the mass–radius relationship, redshift, the Zeldovich condition, energy and causality conditions, the adiabatic index, the Tolman–Oppenheimer–Volkoff equation, the equation of state parameter and compactness. Our findings align with the observational data which indicate that the pulsar SAX J1748.9-2021 is viable and stable under this modified theory of gravity.

The Discovery of Three Pulsars in the Globular Cluster M15 with FAST

We present the discovery of three pulsars in the globular cluster (GC) M15 (NGC 7078) by the Five-hundred-meter Aperture Spherical Radio Telescope. PSR J2129+1210J (M15J) is a millisecond pulsar with a spin period of 11.84 ms and a dispersion measure of 66.68 pc cm−3. Both PSR J2129+1210K and L (M15K and L) are long-period pulsars with spin periods of 1928 ms and 3961 ms, respectively. M15L is the GC pulsar with the longest spin period known. We update the timing solutions of M15A to M15H. As predicted by A. Ridolfi et al., the flux density of M15C keeps decreasing, and the latest detection in our data set was on 2022 December 20. We have also detected M15I’s signal for the first time since its discovery. Current timing suggests that it is an isolated pulsar.

The R-process Alliance: Fifth Data Release from the Search for R-process-enhanced Metal-poor Stars in the Galactic Halo with the GTC* * This paper includes data gathered with the 10.4 m Gran Telescopio Canarias located at La Palma, Canary Islands, Spain.

Understanding the abundance pattern of metal-poor stars and the production of heavy elements through various nucleosynthesis processes offers crucial insights into the chemical evolution of the Milky Way, revealing primary sites and major sources of rapid neutron-capture process (r-process) material in the Universe. In this fifth data release from the R-Process Alliance (RPA), we present the detailed chemical abundances of 41 faint (down to V = 15.8) and extremely metal-poor (down to [Fe/H] = −3.3) halo stars selected from the RPA. We obtained high-resolution spectra for these objects with the HORuS spectrograph on the Gran Telescopio Canarias. We measure the abundances of light, α, Fe-peak, and neutron-capture elements. We report the discovery of five carbon-enhanced metal-poor, one limited-r, three r-I, and four r-II stars, and six Mg-poor stars. We also identify one star of a possible globular cluster origin at an extremely low metallicity at [Fe/H] = −3.0. This adds to the growing evidence of a lower-limit metallicity floor for globular cluster abundances. We use the abundances of Fe-peak elements and the α-elements to investigate the contributions from different nucleosynthesis channels in the progenitor supernovae. We find the distribution of [Mg/Eu] as a function of [Fe/H] to have different enrichment levels, indicating different possible pathways and sites of their production. We also reveal differences in the trends of the neutron-capture element abundances of Sr, Ba, and Eu of various r-I and r-II stars from the RPA data releases, which provide constraints on their nucleosynthesis sites and subsequent evolution.

DESI dark energy time evolution is recovered by cosmologically coupled black holes

Recent baryon acoustic oscillation (BAO) measurements by the Dark Energy Spectroscopic Instrument (DESI) provide evidence that dark energy (DE) evolves with time, as parameterized by a w 0 w a equation of state. Cosmologically coupled black holes (BHs) provide a DE source that naturally evolves with time, because BH production tracks cosmic star-formation. Using DESI BAO measurements and priors informed by Big Bang Nucleosynthesis, we measure the fraction of baryonic density converted into BHs, assuming that all DE is sourced by BH production. We find that the best-fit DE density tracks each DESI best-fit w 0 w a model within 1σ, except at redshifts z ≲ 0.2, highlighting limitations of the w 0 w a parameterization. Cosmologically coupled BHs produce H 0 = (69.94 ± 0.81) km s-1 Mpc-1, with the same χ 2 as ΛCDM, and with two fewer parameters than w 0 w a . This value reduces tension with SH0ES to 2.7σ and is in excellent agreement with recent measurements from the Chicago-Carnegie Hubble Program. Because cosmologically coupled BH production depletes the baryon density established by primordial nucleosynthesis, these BHs provide a physical explanation for the “missing baryon problem” and the anomalously low sum of neutrino masses preferred by DESI.The global evolution of DE is an orthogonal probe of cosmological coupling, complementing constraints on BH mass-growth from elliptical galaxies, stellar binaries, globular clusters, the LIGO-Virgo-KAGRA merging population, and X-ray binaries.A DE density that correlates with cosmic star-formation: 1) is a natural outcome of cosmological coupling in BH populations; 2) eases tension between early and late-time cosmological probes; and 3) produces time-evolution toward a late-time ΛCDM cosmology different from Cosmic Microwave Background projections.

SETI at FAST in China

Abstract Since the commencement of the first SETI observation in 2019, China’s Search for Extraterrestrial Intelligence program has garnered momentum through domestic support and international collaborations. Several observations targeting exoplanets and nearby stars have been conducted with the Five-hundred-metre Aperture Spherical Telescope (FAST). In 2023 the introduction of the Far Neighbour Project (FNP) marked a substantial leap forward, driven by the remarkable sensitivity of the FAST telescope and some of the novel observational techniques. The FNP seeks to methodically detect technosignatures from celestial bodies, including nearby stars, exoplanetary systems, Milky Way globular clusters and more. This paper provides an overview of the progress achieved by SETI in China and offers insights into the distinct phases comprising the FNP. Additionally it underscores the significance of this project’s advancement and its potential contributions to the field.

Measuring Dwarf Galaxy Intrinsic Abundance Scatter with Mid-resolution Spectroscopic Surveys: Calibrating APOGEE Abundance Errors

The first generations of stars left their chemical fingerprints on metal-poor stars in the Milky Way and its surrounding dwarf galaxies. While instantaneous and homogeneous enrichment implies that groups of conatal stars should have the same element abundances, small amplitudes of abundance scatter are seen at fixed [Fe/H]. Measurements of intrinsic abundance scatter have been made with small high-resolution spectroscopic data sets where measurement uncertainty is small compared to this scatter. In this work, we present a method to use mid-resolution survey data, which have larger errors, to make this measurement. Using APOGEE Data Release 17, we calculate the intrinsic scatter of Al, O, Mg, Si, Ti, Ni, and Mn relative to Fe for 333 metal-poor stars across six classical dwarf galaxies around the Milky Way, and 1604 stars across 19 globular clusters (GCs). We calibrate the reported abundance errors in bins of signal-to-noise ratio and [Fe/H] using a high-fidelity halo data set. Applying these calibrated errors to the APOGEE data, we find small amplitudes of average intrinsic abundance scatter in dwarf galaxies ranging from 0.03 to 0.09 dex, with a median value of 0.047 dex. For the GCs, we find intrinsic scatters ranging from 0.01 to 0.11 dex, with particularly high scatter for Al and O. Our measurements of intrinsic abundance scatter place important upper bounds, which are limited by our calibration, on the intrinsic scatter in these systems, as well as constraints on their underlying star formation history and mixing that we can look to simulations to interpret.

JWST’s PEARLS: 119 multiply imaged galaxies behind MACS0416, lensing properties of caustic crossing galaxies, and the relation between halo mass and number of globular clusters at z = 0.4

We present a new lens model for the ɀ = 0.396 galaxy cluster MACS J0416.1 –2403 based on a previously known set of 77 spectroscop-ically confirmed, multiply imaged galaxies plus an additional set of 42 candidate multiply imaged galaxies from past HST and new JWST data. The new galaxies lack spectroscopic redshifts but have geometric and/or photometric redshift estimates that are presented here. The new model predicts magnifications and time delays for all multiple images. The full set of constraints totals 343, constituting the largest sample of multiple images lensed by a single cluster to date. Caustic-crossing galaxies lensed by this cluster are especially interesting. Some of these galaxies show transient events, most of which are interpreted as micro-lensing of stars at cosmological distances. These caustic-crossing arcs are expected to show similar events in future, deeper JWST observations. We provide time delay and magnification models for all these arcs. The time delays and the magnifications for different arcs are generally anti-correlated. In the major sub-halos of the cluster, the dark-matter mass from our lens model correlates well with the observed number of globular clusters, as expected from N-body simulations. This confirms earlier results, derived at lower redshifts, which suggest that globular clusters can be used as powerful mass proxies for the halo masses when lensing constraints are scarce or not available.

JWST imaging of the closest globular clusters

We present observations of the two closest globular clusters, NGC 6121 and NGC 6397, taken with the NIRISS detector of JWST. The combination of our new JWST data with archival Hubble Space Telescope (HST) images allows us to compute proper motions, disentangle cluster members from field objects, and probe the main sequence (MS) of the clusters down to <0.1 M⊙ as well as the brighter part of the white-dwarf sequence. We show that theoretical isochrones fall short in modeling the low-mass MS and discuss possible explanations for the observed discrepancies. Our analysis suggests that the lowest-mass members of both clusters are significantly more metal-rich and oxygen-poor than their higher-mass counterparts. It is unclear whether the difference is caused by a genuine mass-dependent chemical heterogeneity, low-temperature atmospheric processes altering the observed abundances, or systematic shortcomings in the models. We computed the present-day local luminosity and mass functions of the two clusters; our data reveal a strong flattening of the mass function indicative of a significant preferential loss of low-mass stars in agreement with previous dynamical models for these two clusters. We have made our NIRISS astro-photometric catalogs and stacked images publicly available to the community.

Massive star cluster formation

The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the TORCH framework, which uses the Astrophysical MUltipurpose Software Environment (AMUSE) to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgraded TORCH by implementing the N-body code PETAR, thereby enabling TORCH to handle massive clusters forming from 106 M⊙ clouds with ≥105 individual stars. We present results from TORCH simulations of star clusters forming from 104, 105, and 106 M⊙ turbulent spherical gas clouds (named M4, M5, M6) of radius R = 11.7 pc. We find that star formation is highly efficient and becomes more so at a higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities 2.325 × 101,2,3 M⊙ pc−2, after a free-fall time of tff = 6.7,2.1,0.67 Myr, we find that ∼30%, 40%, and 60% of the cloud mass has formed into stars, respectively. The end of simulation-integrated star formation efficiencies for M4, M5, and M6 are ϵ⋆ = M⋆/Mcloud = 36%, 65%, and 85%. Observations of nearby clusters similar in mass and size to M4 have instantaneous star formation efficiencies of ϵinst ≤ 30%, which is slightly lower than the integrated star formation efficiency of M4. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense (Σ ≥ 102 M⊙ pc−2) and massive (Mcloud ≥ 105 M⊙). In such environments, stellar feedback from winds and radiation is not strong enough to counteract the gravity from gas and stars until a majority of the gas has formed into stars.

A coherent view of Li depletion and angular momentum transport to explain the Li plateau – from Population II to Population I stars

Context. The discrepancy between the predictions of Big Bang nucleosynthesis and the lithium abundance observed in the oldest stars of our Galaxy, known as the cosmological lithium problem, has long been regarded as a challenge to the fields of both cosmology and astrophysics. Aims. In light of recent theoretical advances concerning the transport of chemicals and angular momentum in Population I low-mass stars, we re-examine the stellar depletion hypothesis to explain the lithium plateau, which spans a wide range of metallicities over a specific range of stellar effective temperature. Methods. We computed stellar evolution models with the code STAREVOL, including the same input physics that enable self-consistent reproduction of the Li depletion in the Sun and stars in open clusters, while accounting for internal rotation consistent with asteroseismic constraints. In addition to atomic diffusion and parametric turbulence, which were considered in previous studies of Li depletion along the plateau, our models include rotation-induced hydrodynamical processes and additional parametric viscosity for the transport of angular momentum as well as penetrative convection with a rotational dependence, and magnetic braking. Results. As in the case of Pop I stars, the mixing obtained with the current prescriptions for vertical and horizontal shear turbulence induced by rotation is insufficient to reproduce the Li constraints, and parametric turbulence is required. Even if the nature of the turbulence has yet to be identified, we show that the compactness of Pop II low-mass dwarf stars shall naturally lead to similar Li depletion over a large domain in the [Fe/H]–Teff plane, resulting in a plateau with little dispersion. We calibrated the efficiency of the turbulence to fit the abundance of Li in Pop II stars selected from the GALAH DR3 spectroscopic survey and from an homogeneous reanalysis of abundances from the literature. This calibration also enables the reproduction of lithium and magnesium trends in post-turnoff stars of the globular cluster NGC 6752. The same stellar structure considerations consistently explain the observed change of Li depletion and the dispersion regime for [Fe/H] above −1.5 dex, that is, at the transition in metallicity between Pop II to Pop I stars. Conclusions. Our results provide new constraints to the physical processes that transport chemicals and angular momentum in stellar interiors. They offer a comprehensive way to reproduce the observed Li patterns in low-mass dwarf stars across the entire Galactic metallicity range covered by spectroscopic surveys, including the most Fe-poor regime, as supported by the Li value in the non-CEMP star that lies on the plateau at [Fe/H] below −5.8 dex. Our careful analysis of the other very metal-poor stars with lower Li abundances supports the environmental origin of the so-called meltdown regime. Finally, the expected plateau-to-scatter transition pattern further supports the stellar solution to the cosmological problem.

Measuring the center and the dimension of our Galaxy using globular clusters

Abstract Historical measurements, such as the measurement of the size of our Galaxy and distance to its center, are based on relatively simple assumptions, that allow high-school students and undergraduate college students to explore and replicate the basic analysis of astronomical data. Here I present a practical data analysis exercise, that was developed to select the Slovenian participants to the International Olympiad on Astronomy and Astrophysics (IOAA). The exercise guides students to discover the globular cluster distribution within our Galaxy, use basic knowledge of spherical astronomy and coordinate transformation, explore the night sky and develop basic skills for the analysis of astronomical data.

Tidal mass loss in the Fornax dwarf spheroidal galaxy through N-body simulations with Gaia DR3-based orbits

Aims. The Fornax dwarf spheroidal galaxy (dSph) represents a challenge for some globular cluster (GC) formation models, because an exceptionally high fraction of its stellar mass is locked in its GC system. In order to shed light on our understanding of GC formation, we aim to constrain the amount of stellar mass that Fornax has lost via tidal interaction with the Milky Way (MW). Methods. Exploiting the flexibility of effective multi-component N-body simulations and relying on state-of-the-art estimates of Fornax’s orbital parameters, we study the evolution of the mass distribution of the Fornax dSph in observationally justified orbits in the gravitational potential of the MW over 12 Gyr. Results. We find that, though the dark-matter mass loss can be substantial, the fraction of stellar mass lost by Fornax to the MW is always negligible, even in the most eccentric orbit considered. Conclusions. We conclude that stellar-mass loss due to tidal stripping is not a plausible explanation for the unexpectedly high stellar mass of the GC system of the Fornax dSph and we discuss quantitatively the implications for GC formation scenarios.

Euclid Early Release Observations -- Globular clusters in the Fornax galaxy cluster, from dwarf galaxies to the intracluster field

Euclid Early Release Observations -- Globular clusters in the Fornax galaxy cluster, from dwarf galaxies to the intracluster field

Combined Gemini-South and HST photometric analysis of the globular cluster NGC 6558

Context. NGC 6558 is a low-galactic-latitude globular cluster projected in the direction of the Galactic bulge. Due to high reddening, this region presents challenges in deriving accurate parameters, which require meticulous photometric analysis. We present a combined analysis of near-infrared and optical photometry from multi-epoch high-resolution images collected with Gemini-South/GSAOI+GeMS (in the J and KS filters) and HST/ACS (in the F606W and F814W filters). Aims. We aim to refine the fundamental parameters of NGC 6558, utilising high-quality Gemini-South/GSAOI and HST/ACS photometries. Additionally, we intend to investigate its role in the formation of the Galactic bulge. Methods. We performed a meticulous differential reddening correction to investigate the effect of contamination from Galactic bulge field stars. To derive the fundamental parameters – age, distance, reddening, and the total-to-selective coefficient – we employed a Bayesian isochrone fitting. The results from high-resolution spectroscopy and RR Lyrae stars were implemented as priors. For the orbital parameters, we employed a barred Galactic mass model. Furthermore, we analysed the age-metallicity relation to contextualise NGC 6558 within the Galactic bulge’s history. Results. We studied the impact of two differential reddening corrections on the age derivation. When removing as much as possible of the Galactic bulge field star contamination, the isochrone fitting combined with synthetic colour-magnitude diagrams gives a distance of 8.41−0.10+0.11 kpc, an age of 13.0 ± 0.9 Gyr, and a reddening of E(B − V) = 0.34 ± 0.02. We derived a total-to-selective coefficient of RV = 3.2 ± 0.2 thanks to the simultaneous near-infrared–optical synthetic colour-magnitude diagram fitting, which, aside from errors, agrees with the commonly used value. The orbital parameters showed that NGC 6558 is confined within the inner Galaxy and it is not compatible with a bar-shape orbit, indicating that it is a bulge member. Assembling the old and moderately metal-poor ([Fe/H] ∼ −1.1) clusters in the Galactic bulge, we derived their age-metallicity relation with star formation starting at 13.6 ± 0.2 Gyr and effective yields of ρ = 0.05 ± 0.01 Z⊙. Conclusions. The old age derived for NGC 6558 is compatible with other clusters with similar metallicity and a blue horizontal branch in the Galactic bulge, which constitute the moderately metal-poor globular clusters. The age-metallicity relation shows that the starting age of star formation is compatible with the age of NGC 6558, and the chemical enrichment is ten times faster than the ex situ globular cluster branch.

Using 26Al to detect ongoing self-enrichment in young massive star clusters

Abstract Self-enrichment is one of the leading explanations for chemical anomalies in globular clusters. In this scenario, various candidate polluter stars have been proposed to eject gas with altered chemical composition during the self-enrichment process. Most of the proposed polluters will also eject radioactive 26Al into the surroundings. Hence, any detection of 26Al in young massive star clusters (YMCs) would support the self-enrichment scenario if YMCs were indeed the progenitors of globular clusters. Observations of gamma-ray data from COMPTEL and INTEGRAL, as well as detections of 26AlF molecules by the Atacama Large Millimeter-submillimeter Array (ALMA), indicate the maturing of 26Al detection methods. Detection possibilities will be enhanced in the short- to mid-term by the upcoming launch of the Compton Spectrometer and Imager (COSI). The Square Kilometer Array (SKA) could in principle also detect radio recombination lines of the positronium formed from the decay products of 26Al. Here, we show for a sample of YMCs in the nearby Universe, where self-enrichment could plausibly take place. For some nearby galaxies, this could enhance 26Al by an order of one magnitude. Detecting 26AlF with ALMA appears feasible for many candidate self-enrichment clusters, although significant challenges remain with other detection methods. The Large Magellanic Cloud, with its YMC R136, stands out as the most promising candidate. Detecting a 1.8 MeV radioactive decay line of 26Al here would require at least 15 months of targeted observation with COSI, assuming ongoing self-enrichment in R136.

Energy equipartition in multiple-population globular clusters

Abstract We present the results of Monte Carlo simulations aimed at exploring the evolution towards energy equipartition of first- (1G) and second-generation (2G) stars in multiple-population globular clusters and how this evolution is affected by the initial differences between the spatial distributions of the two populations. Our results show that these initial differences have fundamental implications for the evolution towards energy equipartition of the two populations. We find that 2G stars, which are assumed to be initially more centrally concentrated than 1G stars, are generally characterized by a more rapid evolution towards energy equipartition. The evolution towards energy equipartition depends on the velocity dispersion component and is more rapid for the tangential velocity dispersion. The extent of the present-day differences between the degree of energy equipartition of 2G and 1G stars depends on the cluster’s dynamical age and may be more significant in the tangential velocity dispersion and at intermediate distances from the cluster’s center around the half-mass radius.