Simulations of globular clusters within their parent galaxies: multiple stellar populations and internal kinematics

TL;DR

This study uses 3D simulations to explore how multiple stellar populations form in globular clusters within their host galaxies, highlighting the influence of environment and initial conditions on their properties.

Contribution

It introduces a novel simulation approach modeling GC formation from molecular clouds and gas accretion, linking environment to stellar population characteristics.

Findings

The fraction of second-generation stars depends on the cluster's environment.

Simulations reproduce observed abundance patterns and kinematics.

A positive correlation between 2G star fraction and initial cluster mass is found.

Abstract

Using three-dimensional smoothed particle hydrodynamics simulations, we investigate the formation of multiple stellar populations (MSPs) in globular clusters (GCs) within the context of their parent galaxies. In our scenario, the second generation (2G) of stars originate from both asymptotic giant branch (AGB) polluters and pristine gas accreted from the host galaxy. Previous theoretical and numerical studies have demonstrated that this 'AGB with dilution' model has the potential to alleviate several problems faced by the classical AGB scenario. However, the accretion of pristine gas on to the GC has yet to be investigated within the context of the parent galaxy. This paper presents the preliminary results from our original simulation code which models GC formation from giant molecular clouds in a host galaxy, and subsequent gas accretion on to the GC. By simulating the genesis of the 2G over a 370 Myr time frame, we demonstrate that the fraction of 2G stars are inextricably linked to the GC's environment. Our simulations rationalize the wide variety of abundance patterns, kinematics, and 2G concentrations by altering the initial conditions of both the GC progenitor and the host galaxy itself. Most notably, we reproduce a positive correlation between the fraction of 2G stars and the initial mass of the cluster. We discuss the physical implications of our scenario and compare our simulations with observations of the Galactic GC 47 Tucanae (47 Tuc). Finally, we present scaling relations that encompass the wider GC population and serve as a reference for future observations.