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

2020 ◽  
Vol 500 (4) ◽  
pp. 4578-4596
Author(s):  
Madeleine McKenzie ◽  
Kenji Bekki

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.

Author(s):  
Yingtian Chen ◽  
Hui Li ◽  
Mark Vogelsberger

Abstract We perform a suite of hydrodynamic simulations to investigate how initial density profiles of giant molecular clouds (GMCs) affect their subsequent evolution. We find that the star formation duration and integrated star formation efficiency of the whole clouds are not sensitive to the choice of different profiles but are mainly controlled by the interplay between gravitational collapse and stellar feedback. Despite this similarity, GMCs with different profiles show dramatically different modes of star formation. For shallower profiles, GMCs first fragment into many self-gravitation cores and form sub-clusters that distributed throughout the entire clouds. These sub-clusters are later assembled ‘hierarchically’ to central clusters. In contrast, for steeper profiles, a massive cluster is quickly formed at the center of the cloud and then gradually grows its mass via gas accretion. Consequently, central clusters that emerged from clouds with shallower profiles are less massive and show less rotation than those with the steeper profiles. This is because 1) a significant fraction of mass and angular momentum in shallower profiles is stored in the orbital motion of the sub-clusters that are not able to merge into the central clusters 2) frequent hierarchical mergers in the shallower profiles lead to further losses of mass and angular momentum via violent relaxation and tidal disruption. Encouragingly, the degree of cluster rotations in steeper profiles is consistent with recent observations of young and intermediate-age clusters. We speculate that rotating globular clusters are likely formed via an ‘accretion’ mode from centrally-concentrated clouds in the early Universe.


2019 ◽  
Vol 14 (S351) ◽  
pp. 524-527
Author(s):  
Maria A. Tiongco ◽  
Enrico Vesperini ◽  
Anna Lisa Varri

AbstractWe present several results of the study of the evolution of globular clusters’ internal kinematics, as driven by two-body relaxation and the interplay between internal angular momentum and the external Galactic tidal field. Via a large suite of N-body simulations, we explored the three-dimensional velocity space of tidally perturbed clusters, by characterizing their degree of velocity dispersion anisotropy and their rotational properties. These studies have shown that a cluster’s kinematical properties contain distinct imprints of the cluster’s initial structural properties, dynamical history, and tidal environment. Building on this fundamental understanding, we then studied the dynamics of multiple stellar populations in globular clusters, with attention to the largely unexplored role of angular momentum.


2019 ◽  
Vol 622 ◽  
pp. A53 ◽  
Author(s):  
K. Bekki

Context. One of the currently favored scenarios for the formation of globular clusters (GCs) with multiple stellar populations is that an initial massive stellar system forms (“first generation”, FG), subsequently giving rise to gaseous ejecta which is converted into a second-generation (SG) of stars to form a GC. How such GCs with such FG and SG populations form and evolve, however, remains unclear. Aims. We therefore investigate, for the first time, the sequential formation processes of both FG and SG stars from star-forming massive gas clumps in gas-rich dwarf disk galaxies. Methods. We adopt a novel approach to resolve the two-stage formation of GCs in hydrodynamical simulations of dwarf galaxies. In the new simulations, new gas particles that are much less massive than their parent star particle are generated around each new star particle when the new star enters into the asymptotic giant branch (AGB) phase. Furthermore, much finer maximum time step width (~105 yr) and smaller softening length (~2 pc) are adopted for such AGB gas particles to properly resolve the ejection of gas from AGB stars and AGB feedback effects. Therefore, secondary star formation from AGB ejecta can be properly investigated in galaxy-scale simulations. Results. An FG stellar system can first form from a massive gas clump developing due to gravitational instability within its host gas-rich dwarf galaxy. Initially the FG stellar system is not a single massive cluster, but instead is composed of several irregular stellar clumps (or filaments) with a total mass larger than 106 M⊙. While the FG system is dynamically relaxing, gaseous ejecta from AGB stars can be gravitationally trapped by the FG system and subsequently converted into new stars to form a compact SG stellar system within the FG system. Interestingly, about 40% of AGB ejecta is from stars that do not belong to the FG system (“external gas accretion”). FG and SG stellar systems have different amplitudes of internal rotation and V∕σ. The mass-density (MSG−ρSG) relation for SG stellar systems can be approximated as ρSG ∝ MSG1.5. There can be a threshold total mass of GC host galaxies (Mth = [5 − 23] × 109 M⊙) beyond which the formation of GCs with compact SG stellar systems is possible. Both the initial baryonic mass fraction and the gas mass fraction in dwarfs are crucial parameters that determine whether or not GCs can contain multiple stellar populations. GCs with compact SG stellar systems are more likely to form in dwarf disks with larger gas mass fractions and higher surface mass densities. Formation of binary GCs with SGs and the subsequent GC merging are clearly seen in some models. The derived external gas-accretion process in FG systems initially consisting of stellar clumps will need to be investigated further in more sophisticated simulations.


2006 ◽  
Vol 2 (14) ◽  
pp. 436-437
Author(s):  
Alison I. Sills ◽  
Jonathan M. Downing

AbstractWe investigate some implications of having two star formation episodes in globular clusters, rather than the traditional single-burst approximation. Evidence for more than one stellar generation is accumulating in observations of abundances of elements lighter than iron in globular cluster stars, and is thought to imply some self-enrichment of the globular cluster gas. In particular, we explore models based on the assumption that the self-enrichment comes from an early generation of asymptotic giant branch (AGB) stars.


2019 ◽  
Vol 491 (1) ◽  
pp. 440-454 ◽  
Author(s):  
Long Wang ◽  
Pavel Kroupa ◽  
Koh Takahashi ◽  
Tereza Jerabkova

ABSTRACT Many possible scenarios for the formation of multiple stellar populations (MSPs) in globular clusters (GCs) have been discussed so far, including the involvement of asymptotic giant branch stars, fast-rotating main-sequence stars, very massive main-sequence stars and mass-transferring massive binaries based on stellar evolution modelling. But self-consistent, dynamical simulations of very young GCs are usually not considered. In this work, we perform direct N-body modelling of such systems with total masses up to 3.2 × 105 M⊙, taking into account the observationally constrained primordial binary properties, and discuss the stellar mergers driven both by binary stellar evolution and dynamical evolution of GCs. The occurrence of stellar mergers is enhanced significantly in binary-rich clusters such that stars forming from the gas polluted by merger-driven ejection/winds would appear as MSPs. We thus emphasize that stellar mergers can be an important process that connects MSP formation with star cluster dynamics, and that multiple MSP formation channels can naturally work together. The scenario studied here, also in view of a possible top-heavy initial mass function, may be particularly relevant for explaining the high mass fraction of MSPs (the mass budget problem) and the absence of MSPs in young and low-mass star clusters.


2020 ◽  
Vol 495 (2) ◽  
pp. 2247-2264
Author(s):  
Evelyn J Johnston ◽  
Thomas H Puzia ◽  
Giuseppe D’Ago ◽  
Paul Eigenthaler ◽  
Gaspar Galaz ◽  
...  

ABSTRACT Clues to the formation and evolution of nuclear star clusters (NSCs) lie in their stellar populations. However, these structures are often very faint compared to their host galaxy, and spectroscopic analysis of NSCs is hampered by contamination of light from the rest of the system. With the introduction of wide-field integral field unit (IFU) spectrographs, new techniques have been developed to model the light from different components within galaxies, making it possible to cleanly extract the spectra of the NSCs and study their properties with minimal contamination from the light of the rest of the galaxy. This work presents the analysis of the NSCs in a sample of 12 dwarf galaxies in the Fornax Cluster observed with the Multi-Unit Spectroscopic Explorer (MUSE). Analysis of the stellar populations and star formation histories reveal that all the NSCs show evidence of multiple episodes of star formation, indicating that they have built up their mass further since their initial formation. The NSCs were found to have systematically lower metallicities than their host galaxies, which is consistent with a scenario for mass assembly through mergers with infalling globular clusters, whilst the presence of younger stellar populations and gas emission in the core of two galaxies is indicative of in-situ star formation. We conclude that the NSCs in these dwarf galaxies likely originated as globular clusters that migrated to the core of the galaxy that have built up their mass mainly through mergers with other infalling clusters, with gas-inflow leading to in-situ star formation playing a secondary role.


2021 ◽  
Vol 503 (4) ◽  
pp. 5997-6004
Author(s):  
Hagai B Perets ◽  
Paz Beniamini

ABSTRACT Environments of supernovae (SNe) and gamma-ray bursts (GRBs) link their progenitors to the underlying stellar population, providing critical clues for their origins. However, various transients including Ca-rich SNe and short-GRBs, appear to be located at remote locations, far from the stellar population of their host galaxy, challenging our understanding of their origin and/or physical evolution. These findings instigated models suggesting that either large velocity-kicks were imparted to the transient progenitors, allowing them to propagate to large distances and attain their remote locations; or that they formed in dense globular-clusters residing in the haloes. Here we show that instead, large spatial-offsets of such transients are naturally explained by observations of highly extended stellar populations in (mostly early-type) galaxy haloes, typically missed since they can only be identified through ultra-deep/stacked images. Consequently, no large velocity kicks, nor halo globular–cluster environments are required in order to explain the origin of these transients. These findings support thermonuclear explosions on white-dwarfs, for the origins of Ca-rich SNe progenitors, and the existence of small (or zero) kick-velocities given to short-GRB progenitors. Furthermore, since stacked/ultra-deep imaging show that early-type galaxies are more extended than late-type galaxies, studies of transients’ offset-distribution (e.g. type Ia SNe or FRBs) should account for host galaxy-type. Since early-type galaxies contain older stellar populations, transient arising from older stellar populations would have larger fractions of early-type hosts, and consequently larger fractions of large-offset transients. In agreement with our results for short-GRBs and Ca-rich SNe showing different offset distributions in early versus late-type galaxies.


2009 ◽  
Vol 50 (52) ◽  
pp. 67-80 ◽  
Author(s):  
Patrick Heimbach ◽  
Véronique Bugnion

AbstractWe extend the application of control methods to a comprehensive three-dimensional thermomechanical ice-sheet model, SICOPOLIS (SImulation COde for POLythermal Ice Sheets). Lagrange multipliers, i.e. sensitivities, are computed with an exact, efficient adjoint model that has been generated from SICOPOLIS by rigorous application of automatic differentiation. The case study uses the adjoint model to determine the sensitivity of the total Greenland ice volume to various control variables over a 100 year period. The control space has of the order 1.2 × 106 elements, consisting of spatial fields of basal flow parameters, surface and basal forcings and initial conditions. Reliability of the adjoint model was tested through finite-difference perturbation calculations for various control variables and perturbation regions, ascertaining quantitative inferences of the adjoint model. As well as confirming qualitative aspects of ice-sheet sensitivities (e.g. expected regional variations), we detect regions where model sensitivities are seemingly unexpected or counter-intuitive, albeit ‘real’ in the sense of actual model behavior. An example is inferred regions where sensitivities of ice-sheet volume to basal sliding coefficient are positive, i.e. where a local increase in basal sliding parameter increases the ice-sheet volume. Similarly, positive (generally negative) ice temperature sensitivities in certain parts of the ice sheet are found, the detection of which seems highly unlikely if only conventional perturbation experiments had been used. The object of this paper is largely a proof of concept. Available adjoint-code generation tools now open up a variety of novel model applications, notably with regard to sensitivity and uncertainty analyses and ice-sheet state estimation or data assimilation.


2016 ◽  
Vol 12 (S323) ◽  
pp. 227-230
Author(s):  
Muhammad Akashi ◽  
Noam Soker

AbstractWe conduct three-dimensional hydrodynamical simulations of two opposite jets launched obliquely to the orbital plane around an asymptotic giant branch (AGB) star and within its dense wind, and demonstrate the formation of a ‘messy’ planetary nebula (PN), namely, a PN lacking any type of symmetry (highly irregular). In building the initial conditions we assume that a tight binary system orbits the AGB star, and that the orbital plane of the tight binary system is inclined to the orbital plane of binary system and the AGB star. We further assume that the accreted mass onto the tight binary system forms an accretion disk around one of the stars, and that the plane of the disk is in between the two orbital planes. The highly asymmetrical lobes that we obtain support the notion that messy PNe might be shaped by triple stellar systems.


2020 ◽  
Vol 492 (4) ◽  
pp. 4874-4883 ◽  
Author(s):  
Duncan A Forbes ◽  
Adebusola Alabi ◽  
Aaron J Romanowsky ◽  
Jean P Brodie ◽  
Nobuo Arimoto

ABSTRACT Ultra-diffuse galaxies (UDGs) reveal extreme properties. Here, we compile the largest study to date of 85 globular cluster (GC) systems around UDGs in the Coma cluster, using new deep ground-based imaging of the known UDGs and existing imaging from the Hubble Space Telescope of their GC systems. We find that the richness of GC systems in UDGs generally exceeds that found in normal dwarf galaxies of the same stellar mass. These GC-rich UDGs imply haloes more massive than expected from the standard stellar mass–halo mass relation. The presence of such overly massive haloes presents a significant challenge to the latest simulations of UDGs in cluster environments. In some exceptional cases, the mass in the GC system is a significant fraction of the stellar content of the host galaxy. We find that rich GC systems tend to be hosted in UDGs of lower luminosity, smaller size, and fainter surface brightness. Similar trends are seen for normal dwarf galaxies in the Coma cluster. A toy model is presented in which the GC-rich UDGs are assumed to be ‘failed’ galaxies within massive haloes that have largely old, metal-poor, alpha-element-enhanced stellar populations. On the other hand, GC-poor UDGs are more akin to normal, low surface brightness dwarfs that occupy less massive dark matter haloes. Additional data on the stellar populations of UDGs with GC systems will help to further refine and test this simplistic model.


Sign in / Sign up

Export Citation Format

Share Document