-
Internal Dynamics of Multiple Populations in 28 Galactic Globular Clusters: A Wide-Field study with Gaia and the Hubble Space Telescope
Authors:
Giacomo Cordoni,
Luca Casagrande,
Antonino Milone,
Emanuele Dondoglio,
Alessandra Mastrobuono-Battisti,
Sohee Jang,
Anna Marino,
Edoardo Lagioia,
Maria Vittoria Legnardi,
Tuila Ziliotto,
Fabrizio Muratore,
Vernica Mehta,
Elena Lacchin,
Marco Tailo
Abstract:
We present the first comprehensive analysis of the internal dynamics of multiple stellar populations (MPs) in 28 Galactic Globular Clusters (GCs) across a wide field of view, extending from the innermost regions to the clusters' outskirts. Using astro-photometric catalogs from ground-based observations, Gaia, and the Hubble Space Telescope (HST), we identify first- (1P) and second-population (2P)…
▽ More
We present the first comprehensive analysis of the internal dynamics of multiple stellar populations (MPs) in 28 Galactic Globular Clusters (GCs) across a wide field of view, extending from the innermost regions to the clusters' outskirts. Using astro-photometric catalogs from ground-based observations, Gaia, and the Hubble Space Telescope (HST), we identify first- (1P) and second-population (2P) stars, and study the internal dynamics of MPs using high-precision Gaia DR3 and HST proper motions. Our results reveal that while the 1P transitions from isotropy to slight tangential anisotropy toward the outer regions, 2P stars become increasingly radially anisotropic beyond the half-light radius. We also explore the connection between the dynamics of MPs and the clusters' structural and dynamical properties, finding statistically significant differences in the anisotropy profiles of dynamically young and non-relaxed clusters, particularly beyond the 1-2 half-light radii. In these regions, 1P stars transition from isotropic to slightly tangentially anisotropic motion, while 2P stars become more radially anisotropic. In contrast, dynamically older clusters, with mixed MPs, exhibit weaker relative differences. Furthermore, clusters with orbits closer to the Galactic center exhibit larger dynamical differences between 1P and 2P stars than those with larger peri-Galactic radii. These findings are consistent with a scenario where 2P stars form in a more centrally concentrated environment, where the interaction with the Milky Way tidal field plays a crucial role in the dynamical evolution of MPs, especially of 1P.
△ Less
Submitted 3 September, 2024;
originally announced September 2024.
-
The boring history of Gaia BH3 from isolated binary evolution
Authors:
Giuliano Iorio,
Stefano Torniamenti,
Michela Mapelli,
Marco Dall'Amico,
Alessandro A. Trani,
Sara Rastello,
Cecilia Sgalletta,
Stefano Rinaldi,
Guglielmo Costa,
Bera A. Dhal-Lahtinen,
Gaston J. Escobar,
Erika Korb,
M. Paola Vaccaro,
Elena Lacchin,
Benedetta Mestichelli,
Ugo Niccolò di Carlo,
Mario Spera,
Manuel Arca Sedda
Abstract:
Gaia BH3 is the first observed dormant black hole (BH) with a mass of $\approx{30}$ M$_\odot$ and represents the first confirmation that such massive BHs are associated with metal-poor stars. Here, we explore the isolated binary formation channel for Gaia BH3 focusing on the old and metal-poor stellar population of the Milky Way halo. We use the MIST stellar models and our open-source population s…
▽ More
Gaia BH3 is the first observed dormant black hole (BH) with a mass of $\approx{30}$ M$_\odot$ and represents the first confirmation that such massive BHs are associated with metal-poor stars. Here, we explore the isolated binary formation channel for Gaia BH3 focusing on the old and metal-poor stellar population of the Milky Way halo. We use the MIST stellar models and our open-source population synthesis code SEVN to evolve $5.6 \times 10^8$ binaries exploring 20 sets of parameters. We find that systems like Gaia BH3 form preferentially from binaries initially composed of a massive star ($40-60$ M$_\odot$) and a low mass companion ($<1$ M$_\odot$) in a wide ($P>10^3$ days) and eccentric orbit ($e>0.6$). Such progenitor binary stars do not undergo any Roche-lobe overflow episode during their entire evolution, so that the final orbital properties of the BH-star system are determined at the core collapse of the primary star. Low natal kicks ($\lesssim$ 10~km/s) significantly favour the formation of Gaia BH3-like systems, but high velocity kicks up to $\approx 220$ km/s are also allowed. We estimate the formation efficiency for Gaia BH3-like systems in old ($t>10$ Gyr) and metal-poor ($Z<0.01$) populations to be $\sim 4 \times 10^{-8}$ M$_\odot^{-1}$ (for our fiducial model), representing $\sim 3\%$ of the whole simulated BH-star population. We expect up to $\approx 4000$ BH-star systems in the Galactic halo formed through isolated evolution, of which $\approx 100$ are compatible with Gaia BH3-like. Given the density profile of the Galactic halo we do not expect more than one at the observed distance of Gaia BH3. Our models show that, even if it was born inside a stellar cluster, Gaia BH3 is compatible with a primordial binary star that escaped from its parent cluster without experiencing significant dynamical interactions.
△ Less
Submitted 5 August, 2024; v1 submitted 26 April, 2024;
originally announced April 2024.
-
Multiple stellar population mass loss in massive Galactic globular clusters
Authors:
Elena Lacchin,
Alessandra Mastrobuono-Battisti,
Francesco Calura,
Carlo Nipoti,
Antonino P. Milone,
Massimo Meneghetti,
Eros Vanzella
Abstract:
The degree of mass loss, i.e. the fraction of stars lost by globular clusters, and specifically by their different populations, is still poorly understood. Many scenarios of the formation of multiple stellar populations, especially the ones involving self-enrichment, assume that the first generation (FG) was more massive at birth than now to reproduce the current mass of the second generation (SG)…
▽ More
The degree of mass loss, i.e. the fraction of stars lost by globular clusters, and specifically by their different populations, is still poorly understood. Many scenarios of the formation of multiple stellar populations, especially the ones involving self-enrichment, assume that the first generation (FG) was more massive at birth than now to reproduce the current mass of the second generation (SG). This assumption implies that, during their long-term evolution, clusters lose around $90\%$ of the FG. We have tested whether such strong mass loss could take place in a massive globular cluster orbiting the Milky Way at $4\ {\rm kpc}$ from the centre and composed of two generations. We perform a series of $N$-body simulations for ${12\ \rm Gyr}$ to probe the parameter space of internal cluster properties. We have derived that, for an extended FG and a low-mass second one, the cluster loses almost $98\%$ of its initial FG mass and the cluster mass can be as much as 20 times lower after a Hubble time. Furthermore, under these conditions, the derived fraction of SG stars, $f_{\rm enriched}$, falls in the range occupied by observed clusters of similar mass ($\sim 0.6-0.8$). In general, the parameters that affect the most the degree of mass loss are the presence or not of primordial segregation, the depth of the central potential, $W_{0,FG}$, the initial mass of the SG, $M^{ini}_{SG}$, and the initial half-mass radius of the SG, $r_{h,SG}$. Higher $M^{ini}_{SG}$ have not been found to imply higher final $f_{\rm enriched}$ due to the deeper cluster potential well which slows down mass loss.
△ Less
Submitted 26 September, 2023;
originally announced September 2023.
-
Shaping the unseen: the influence of baryons and environment on low-mass, high-redshift dark matter haloes in the SIEGE simulations
Authors:
R. Pascale,
F. Calura,
A. Lupi,
J. Rosdahl,
E. Lacchin,
M. Meneghetti,
C. Nipoti,
E. Vanzella,
E. Vesperini,
A. Zanella
Abstract:
We use zoom-in, hydrodynamical, cosmological $N$-body simulations tracing the formation of the first stellar clumps from the SImulating the Environments where Globular clusters Emerged (SIEGE) project, to study key structural properties of dark matter haloes when the Universe was only $0.92$ Gyr old. The very high-resolution (maximum physical resolution 0.3 h$^{-1}$ pc at $z=6.14$, smallest dark-m…
▽ More
We use zoom-in, hydrodynamical, cosmological $N$-body simulations tracing the formation of the first stellar clumps from the SImulating the Environments where Globular clusters Emerged (SIEGE) project, to study key structural properties of dark matter haloes when the Universe was only $0.92$ Gyr old. The very high-resolution (maximum physical resolution 0.3 h$^{-1}$ pc at $z=6.14$, smallest dark-matter particle mass $164\,M_{\odot}$) allows us to reach the very low mass end of the stellar-to-halo mass relation ($M_{\rm vir}=10^{7.5-9.5}\,M_{\odot}$) to study the processes that mould dark matter haloes during the first stages of structure formation. We investigate the role of baryonic cooling and stellar feedback, modeled from individual stars, in shaping haloes, and of environmental effects as accretion of dark matter along cosmic filaments and mergers. We find that the onset of star formation (typically for $\log M_{\rm vir}/M_{\odot}\simeq7.6$) causes the inner cusp in the haloes density profile to flatten into a core with constant density and size proportionally to the halo virial mass. Even at these mass scales, we confirm that baryons make haloes that have formed stars rounder in the central regions than haloes that have not formed stars yet, with median minor-to-major $\langle q \rangle$ and intermediate-to-major $\langle s \rangle$ axes 0.66 and 0.84, respectively. Our morphological analysis shows that, at $z=6.14$, haloes are largely prolate in the outer parts, with the major axis aligned along filaments of the cosmic web or towards smaller sub-haloes, with the degree of elongation having no significant dependence on the halo mass.
△ Less
Submitted 4 September, 2023;
originally announced September 2023.
-
The role of rotation on the formation of second generation stars in globular clusters
Authors:
Elena Lacchin,
Francesco Calura,
Enrico Vesperini,
Alessandra Mastrobuono-Battisti
Abstract:
By means of 3D hydrodynamic simulations, we explore the effects of rotation in the formation of second-generation (SG) stars in globular clusters (GC). Our simulations follow the SG formation in a first-generation (FG) internally rotating GC; SG stars form out of FG asymptotic giant branch (AGB) ejecta and external pristine gas accreted by the system. We have explored two different initial rotatio…
▽ More
By means of 3D hydrodynamic simulations, we explore the effects of rotation in the formation of second-generation (SG) stars in globular clusters (GC). Our simulations follow the SG formation in a first-generation (FG) internally rotating GC; SG stars form out of FG asymptotic giant branch (AGB) ejecta and external pristine gas accreted by the system. We have explored two different initial rotational velocity profiles for the FG cluster and two different inclinations of the rotational axis with respect to the direction of motion of the external infalling gas, whose density has also been varied. For a low (10^-24 g cm^-3) external gas density, a disk of SG helium-enhanced stars is formed. The SG is characterized by distinct chemo-dynamical phase space patterns: it shows a more rapid rotation than the FG with the helium-enhanced SG subsystem rotating more rapidly than the moderate helium-enhanced one. In models with high external gas density (10^-23 g cm^-3), the inner SG disc is disrupted by the early arrival of external gas and only a small fraction of highly enhanced helium stars preserves the rotation acquired at birth. Variations in the inclination angle between the rotation axis and the direction of the infalling gas and the velocity profile can slightly alter the extent of the stellar disc and the rotational amplitude. No significant variation has been found in the timespan of our simulations when changing the inclination angle between the rotation axis and the direction of the infalling gas, while different velocity profiles can slightly alter the extent of the stellar disc and the rotational amplitude. The results of our simulations illustrate the complex link between dynamical and chemical properties of multiple populations and provide new elements for the interpretation of observational studies and future investigations of the dynamics of multiple-population GCs.
△ Less
Submitted 12 September, 2022;
originally announced September 2022.
-
Sub-parsec resolution cosmological simulations of star-forming clumps at high redshift with feedback of individual stars
Authors:
F. Calura,
A. Lupi,
J. Rosdahl,
E. Vanzella,
M. Meneghetti,
P. Rosati,
E. Vesperini,
E. Lacchin,
R. Pascale,
R. Gilli
Abstract:
We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass f…
▽ More
We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified 'delayed cooling' feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n~10^3-10^5 cm^{-3}). The sites where star formation ignites are characterised by maximum densities of the order of 10^5 cm^{-3} and gravitational pressures P/k>10^7 K/cm^3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z=6.14 is 3.4x10^7 M_sun, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of ~10^6 M_sun and half-mass sizes of ~100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one order of magnitude lower than the observed one. In the size-mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies.
△ Less
Submitted 16 August, 2022; v1 submitted 27 June, 2022;
originally announced June 2022.
-
On the role of Type Ia supernovae in the second generation star formation in globular clusters
Authors:
Elena Lacchin,
Francesco Calura,
Enrico Vesperini
Abstract:
By means of 3D hydrodynamic simulations, we study how Type Ia supernovae (SNe) explosions affect the star formation history and the chemical properties of second generation (SG) stars in globular clusters (GC). SG stars are assumed to form once first generation asymptotic giant branch (AGB) stars start releasing their ejecta; during this phase, external gas is accreted by the system and SNe Ia beg…
▽ More
By means of 3D hydrodynamic simulations, we study how Type Ia supernovae (SNe) explosions affect the star formation history and the chemical properties of second generation (SG) stars in globular clusters (GC). SG stars are assumed to form once first generation asymptotic giant branch (AGB) stars start releasing their ejecta; during this phase, external gas is accreted by the system and SNe Ia begin exploding, carving hot and tenuous bubbles. Given the large uncertainty on SNe Ia explosion times, we test two different values for the 'delay time'. We run two different models for the external gas density: in the low-density scenario with short delay time, the explosions start at the beginning of the SG star formation, halting it in its earliest phases. The external gas hardly penetrates the system, therefore most SG stars present extreme helium abundances (Y > 0.33). The low-density model with delayed SN explosions has a more extended SG star formation epoch and includes SG stars with modest helium enrichment. On the contrary, the high-density model is weakly affected by SN explosions, with a final SG mass similar to the one obtained without SNe Ia. Most of the stars form from a mix of AGB ejecta and pristine gas and have a modest helium enrichment. We show that gas from SNe Ia may produce an iron spread of $\sim 0.14$ dex, consistent with the spread found in about 20% of Galactic GCs, suggesting that SNe Ia might have played a key role in the formation of this sub-sample of GCs.
△ Less
Submitted 16 July, 2021;
originally announced July 2021.
-
Chemical evolution of ultra-faint dwarf galaxies: testing the IGIMF
Authors:
E. Lacchin,
F. Matteucci,
F. Vincenzo,
M. Palla
Abstract:
We test the integrated galactic initial mass function (IGIMF) on the chemical evolution of 16 ultra-faint dwarf (UFD) galaxies discussing in detail the results obtained for three of them: Boötes I, Boötes II and Canes Venatici I, taken as prototypes of the smallest and the largest UFDs. These objects have very small stellar masses ($\sim 10^3-10^4 \mathrm{M_{\odot}}$) and quite low metallicities (…
▽ More
We test the integrated galactic initial mass function (IGIMF) on the chemical evolution of 16 ultra-faint dwarf (UFD) galaxies discussing in detail the results obtained for three of them: Boötes I, Boötes II and Canes Venatici I, taken as prototypes of the smallest and the largest UFDs. These objects have very small stellar masses ($\sim 10^3-10^4 \mathrm{M_{\odot}}$) and quite low metallicities ([Fe/H]$<-1.0$ dex). We consider four observational constraints: the present-day stellar mass, the [$α$/Fe] vs. [Fe/H] relation, the stellar metallicity distribution function and the cumulative star formation history. Our model follows in detail the evolution of several chemical species (H, He, $α$-elements and Fe). We take into account detailed nucleosynthesis and gas flows (in and out). Our results show that the IGIMF, coupled with the very low star formation rate predicted by the model for these galaxies ($\sim 10^{-4}-10^{-6}\ \mathrm{M_{\odot}yr^{-1}}$), cannot reproduce the main chemical properties, because it implies a negligible number of core-collapse SNe and even Type Ia SNe, the most important polluters of galaxies. On the other hand, a constant classical Salpeter IMF gives the best agreement with data, but we cannot exclude that other formulations of the IGIMF could reproduce the properties of these galaxies. Comparing with Galaxy data we suggest that UFDs could not be the building blocks of the entire Galactic halo, although more data are necessary to draw firmer conclusions.
△ Less
Submitted 26 February, 2020; v1 submitted 19 November, 2019;
originally announced November 2019.
-
The influence of a top-heavy integrated galactic IMF and dust on the chemical evolution of high-redshift starbursts
Authors:
Marco Palla,
Francesco Calura,
Francesca Matteucci,
XiLong Fan,
Fiorenzo Vincenzo,
Elena Lacchin
Abstract:
We study the effects of the integrated galactic initial mass function (IGIMF) and dust evolution on the abundance patterns of high redshift starburst galaxies. In our chemical models, the rapid collapse of gas clouds triggers an intense and rapid star formation episode, which lasts until the onset of a galactic wind, powered by the thermal energy injected by stellar winds and supernova explosions.…
▽ More
We study the effects of the integrated galactic initial mass function (IGIMF) and dust evolution on the abundance patterns of high redshift starburst galaxies. In our chemical models, the rapid collapse of gas clouds triggers an intense and rapid star formation episode, which lasts until the onset of a galactic wind, powered by the thermal energy injected by stellar winds and supernova explosions. Our models follow the evolution of several chemical elements (C, N, $α$-elements and Fe) both in the gas and dust phases. %The most recent stellar yield and dust prescriptions are adopted. We test different values of $β$, the slope of the embedded cluster mass function for the IGIMF, where lower $β$ values imply a more top-heavy initial mass function (IMF). The computed abundances are compared to high-quality abundance measurements obtained in lensed galaxies and from composite spectra in large samples of star-forming galaxies in the redshift range $2 \lesssim z \lesssim 3$. The adoption of the IGIMF causes a sensible increase of the rate of star formation with respect to a standard Salpeter IMF, with a strong impact on chemical evolution. We find that in order to reproduce the observed abundance patterns in these galaxies, either we need a very top-heavy IGIMF ($β< 2$) or large amounts of dust. In particular, if dust is important, the IGIMF should have $β\ge 2$, which means an IMF slightly more top-heavy than the Salpeter one. The evolution of the dust mass with time for galaxies of different mass and IMF is also computed, highlighting that the dust amount increases with a top-heavier IGIMF.
△ Less
Submitted 24 March, 2020; v1 submitted 19 August, 2019;
originally announced August 2019.
-
On the Origin of Dust in Galaxy Clusters at Low to Intermediate Redshift
Authors:
Eda Gjergo,
Francesca Matteucci,
Marco Palla,
Andrea Biviano,
Elena Lacchin,
XiLong Fan
Abstract:
Stacked analyses of galaxy clusters at low-to-intermediate redshift show signatures attributable to dust, but the origin of this dust is uncertain. We test the hypothesis that the bulk of cluster dust derives from galaxy ejecta. To do so, we employ dust abundances obtained from detailed chemical evolution models of galaxies. We integrate the dust abundances over cluster luminosity functions (one-s…
▽ More
Stacked analyses of galaxy clusters at low-to-intermediate redshift show signatures attributable to dust, but the origin of this dust is uncertain. We test the hypothesis that the bulk of cluster dust derives from galaxy ejecta. To do so, we employ dust abundances obtained from detailed chemical evolution models of galaxies. We integrate the dust abundances over cluster luminosity functions (one-slope and two-slope Schechter functions). We consider both a hierarchical scenario of galaxy formation and an independent evolution of the three main galactic morphologies: elliptical/S0, spiral and irregular. We separate the dust residing within galaxies from the dust ejected in the intracluster medium. To the latter, we apply thermal sputtering. The model results are compared to low-to-intermediate redshift observations of dust masses. We find that in any of the considered scenarios, elliptical/S0 galaxies contribute negligibly to the present-time intracluster dust, despite producing the majority of gas-phase metals in galaxy clusters. Spiral galaxies, instead, provide both the bulk of the spatially-unresolved dust and of the dust ejected into the intracluster medium. The total dust-to-gas mass ratio in galaxy clusters amounts to $10^{-4}$, while the intracluster medium dust-to-gas mass ratio amounts to $10^{-6}$ at most. These dust abundances are consistent with the estimates of cluster observations at $0.2 < z <1$. We propose that galactic sources, spiral galaxies in particular, are the major contributors to the cluster dust budget.
△ Less
Submitted 13 February, 2020; v1 submitted 12 July, 2019;
originally announced July 2019.