We present some results from the DRAGON simulations, a set of four direct N-body simulations of globular clusters (GCs) with a million stars and five percent initial (primordial) binaries. These simulations were undertaken with the NBODY6++GPU code, which allowed us to follow dynamical and stellar evolution of individual stars and binaries, formation and evolution of white dwarfs, neutron stars, and black holes, and the effect of a galactic tidal field. The simulations are the largest existing models of a realistic globular cluster over its full lifetime of 12 billion years. In particular we will show here an investigation of the population of binaries including compact objects (such as white dwarfs - cataclysmic variables and merging black hole binaries in the model as counterparts of LIGO/Virgo sources); their distribution in the cluster and evolution with time.

We review spectroscopic results concerning multiple stellar populations in globular clusters. The cluster initial mass is the most important parameter determining the fraction of second generation stars. The threshold for the onset of the multiple population phenomenon is 1–3×105 M⊙. Nucleosynthesis is influenced by metallicity: Na/O and Mg/Al anti-correlations are more extended in metal-poor than in metal-rich clusters. Massive clusters are more complex systems than the smaller ones, with several populations characterized by different chemical compositions. The high Li abundance observed in the intermediate second generation stars strongly favours intermediate mass AGB stars as polluters for this class of stars; however, it is well possible that the polluters of extreme second generation stars, that often do not have measurable Li, may be fast rotating massive stars or super-massive stars. The mass budget factor should be a function of the cluster mass, and needs to be large only in massive clusters.

Over the last decade, our abilities to observe the internal kinematics of star clusters have drastically increased. Where a few years back only small numbers of bright stars with radial velocity measurements were available, we can nowadays study the three-dimensional motions of large stellar samples, thanks to the combined datasets gathered by state-of-art spectrographs and astrometric satellites. In this work, I summarise the contribution of integral-field spectrographs, in particular MUSE, to this paradigm change. Using dedicated software tools, we were able to overcome a fundamental limitation of spectroscopy and advance to the crowded cluster centres. This allowed us to study the central kinematics in unprecedented detail and to start uncovering the populations of black holes that reside in massive star clusters.

In order to investigate the origin of multiple populations in globular clusters (GCs), we have constructed new chemical evolution models for proto-GCs where the supernova blast waves undergo blowout without expelling the ambient gas. Chemical enrichments in our models are then dictated by the winds of massive stars together with the asymptotic-giant-branch stars ejecta. We find that the observed Na-O anti-correlation can be reproduced when multiple episodes of starburst and enrichment are allowed to continue in proto-GCs. The “mass budget problem” is mostly resolved by our models without ad-hoc assumptions on star formation efficiency, initial mass function, and significant loss of first-generation stars. Interestingly, ages and chemical abundances predicted by this chemical evolution model are in good agreements with those independently obtained from our stellar evolution model for the horizontal-branch. We also discuss observational evidence for the GC-like multiple populations in the Milky Way bulge.

The observational properties of a special class of stars (the so-called Blue Straggler stars - BSSs) in Globular Clusters are discussed in the framework of using this stellar population as probe of the dynamical processes occurring in high-density stellar systems. In particular, the shape of the BSS radial distribution and their level of central segregation have been found to be powerful tracers of the level of the dynamical evolution of the hosting cluster, thus allowing the definition of an empirical chronometer able to measure the dynamical age of star clusters.

Recent investigations of multiple stellar populations in globular clusters (GCs) suggest that the horizontal-branch (HB) morphology and mean period of type ab RR Lyrae variables are mostly sensitive to helium abundance, while the star formation timescale has the greatest effect on our chemical evolution model constructed to reproduce the Na-O anti-correlation of GCs. Therefore, by combining the results from synthetic HB model with those from chemical evolution model, we could put better constraints on star formation history and chemical evolution in GCs with multiple populations. From such efforts made for four GCs, M4, M5, M15, and M80, we find that consistent results can be obtained from these two independent models.

The features and make up of the population of X-ray sources in Galactic star clusters reflect the properties of the underlying stellar environment. Cluster age, mass, stellar encounter rate, binary frequency, metallicity, and maybe other properties as well, determine to what extent we can expect a contribution to the cluster X-ray emission from low-mass X-ray binaries, millisecond pulsars, cataclysmic variables, and magnetically active binaries. Sensitive X-ray observations withXMM-Newton and certainlyChandra have yielded new insights into the nature of individual sources and the effects of dynamical encounters. They have also provided a new perspective on the collective X-ray properties of clusters, in which the X-ray emissivities of globular clusters and old open clusters can be compared to each other and to those of other environments. I will review our current understanding of cluster X-ray sources, focusing on star clusters older than about 1 Gyr, illustrated with recent results.

We present a scenario for the formation of super star clusters (with masses larger than 105 M⊙) in which multiple generations of star formation will occur. We stress that the gas left over (∼50%) from first generation (1G) star formation should be retained in such massive clusters (thanks to their deep potential wells, with escape speeds larger than 10 km/s) and be available for a second or even third generation of stars, with the basic HeCNONaMgAl chemical anomalies observed in globular clusters, the latter assumed to be the descendents of these super star clusters. One new feature of this model is the role of C+ cooling of the dense warm trapped neutral or ionized gas which defines a characteristic temperature of ∼100 K, leading to a second generation (2G) of stars with a top-heavy IMF (M > 5 M⊙). The ashes of the 2G very massive stars (VMS, M > 100 M⊙) sampled in this IMF quickly pollute and dilute the left-over pristine gas with their slow winds (that cannot escape the cluster), while the majority of massive stars develop fast winds (that actually can escape from the cluster). Meanwhile, much of the remaining dense T = 100 K gas contracts gravitationally in the massive cluster and may reach densities of the order of 109 cm−3, in which case the Jeans mass drops to about 0.2 M⊙ and leads to a substantial low-mass pre-MS 3G population (most likely on a very short timescale). In this way, we may solve both the mass budget and the excess Helium problem in proto-globular clusters, while also explaining the Na-O and Mg-Al anti-correlations resulting from hot H-burning of very massive stars at 45MK and 75MK, respectively.

NGC 6402 is one of the most massive globular clusters in the Galaxy but until recently little was known about its detailed chemical composition. Interestingly, recent results have shown that NGC 6402 exhibits a paucity of intermediate composition stars that may be indicative of an early termination of star formation. As a result, NGC 6402 may be important for understanding cluster formation and the order in which various stellar populations are born.

As self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, leading to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward “core collapse,” the globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the core-collapsed and non-core-collapsed clusters. This suggests an internal energy source is at work, delaying the onset of core collapse in many clusters. Over the past decade, a large amount of work has suggested that stellar black holes (BHs) play a dynamically-significant role in clusters throughout their entire lifetimes. Here we review our latest understanding of BH populations in GCs and demonstrate that, through their dynamical interaction with their host cluster, BHs can naturally explain the distinction between core-collapsed and non-core-collapsed clusters through a process we call “black hole burning.”

Integrated light (IL) spectroscopy enables studies of stellar populations beyond the Milky Way and its nearest satellites. In this paper, I will review how IL spectroscopy reveals essential information about globular clusters and the assembly histories of their host galaxies, concentrating particularly on the metallicities and detailed chemical abundances of the GCs in M31. I will also briefly mention the effects of multiple populations on IL spectra, and how observations of distant globular clusters help constrain the source(s) of light-element abundance variations. I will end with future perspectives, emphasizing how IL spectroscopy can bridge the gap between Galactic and extragalactic astronomy.

Using data from the core of 47 Tuc we have identified stars in different evolutionary stages in the colour-magnitude diagram, and used the effects of mass segregation on their radial distribution to study the evolution and origin of blue stragglers (BSS). We separate the BSS into 2 samples by their magnitude and find considerable differences in their distribution. Bright BSS are more centrally concentrated with mass estimates over twice the turn-off mass suggesting an origin involving a triple or multiple star system. The distribution of the faint BSS is close to that of the main-sequence (MS) binaries pointing to these stars as their likely progenitors. Using MESA models, we calculate the expected number of stars in each evolutionary stage and compare it with the observed number of stars. Results indicate that BSS have a post-MS evolution comparable to that of a normal star of the same mass and a MS-BSS lifetime of about 200 – 300 Myr.

We discuss a meta-analysis of the association of abundance variations in globular cluster stars with the present-day stellar mass and metallicity of globular clusters. Using data for 42 globular clusters that are well-sampled from either or both of prior literature studies and the APOGEE survey, we confirm prior findings that increasing aluminum abundance variations in globular clusters are positively correlated with increasing present-day stellar mass or decreasing metallicity. We also demonstrate that the ratio of aluminum abundance variations to either nitrogen abundance variations or sodium abundance variations is itself positively correlated with decreasing metallicity and increasing stellar mass of globular clusters. This suggests that there were at least two non-supernovae chemical polluters that were active in the early universe.

It is generally recognized that massive galaxies form through a combination of in-situ collapse and ex-situ accretion. The in-situ component forms early, where gas collapse and compaction leads to the formation of massive compact systems (blue and red “nuggets”) seen at z > 1. The subsequent accretion of satellites brings in ex-situ material, growing these nuggets in size and mass to appear as the massive early-type galaxies (ETGs) we see locally. Due to stochasticity in the accretion process, in a few rare cases a red nugget will evolve to the present day having undergone little ex-situ mass accretion. The resulting massive, compact and ancient objects have been termed “relic galaxies”. Detailed stellar population and kinematic analyses are required to characterise these systems. However, an additional crucial aspect lies in determining the fraction of ex-situ mass they have accreted since their formation. Globular cluster systems can be used to constrain this fraction, since the oldest and most metal-poor globular clusters in massive galaxies are primarily an accreted, ex-situ population. Models for the formation of relic galaxies and their globular cluster systems suggest that, due to their early compaction and limited accretion of dark-matter dominated satellites, relic galaxies should have characteristically low dark-matter mass fractions compared to ETGs of the same stellar mass.

Binary evolution can produce different blue-straggler binaries, for example, blue stragglers with a bright, red component, or with a faint, blue component. In globular clusters, these blue-straggler binaries are generally observed as a single star, because two components can not be distinguished. Therefore, these blue-straggler binaries can be located in different regions of the color-magnitude diagram of globular clusters, e.g. blue sequence and red sequence observed in M30. We suggest that binary evolution can contribute to the blue stragglers in both of the sequences. Some blue stragglers in the blue sequence may have a faint white dwarf companion, while the red sequence includes some binaries experiencing mass transfer. It should be noted that the red sequence may also have other binaries, for example, the binaries just finished the mass transfer, and the binaries including a blue straggler (the accretors) that have evolved away from the blue sequence.

Globular clusters (GCs) — compact and massive star clusters found ubiquitously around galaxies — are believed to be ancient relics (ages ≳ 10 Gyr) from the early formative phase of galaxies, although their physical origin remains widely debated. The most numerous GC populations are hosted by giant elliptical galaxies, where they can exhibit a broad dispersion in colour interpreted as a wide spread in metallicity. Here, we show that many thousands of similarly compact and massive super star clusters have formed at an approximately steady rate over, at least, the past ~1 Gyr around the nearby giant elliptical galaxy, NGC 1275, at the centre of the Perseus cluster. The number distribution of these young star clusters appears to exhibit a similar dependence in luminosity and mass as the even more numerous but older GCs around NGC 1275. In just a few Gyr, these super star clusters will evolve to become indistinguishable in broadband optical colours from the older GCs, and their spread in age add to the dispersion in colour of these GCs. The spatial distribution of the super star clusters resembles the filamentary network of multiphase gas in the cluster core, implying that they formed from molecular gas amassed from cooling of the hot intracluster gas. The sustained formation of super star clusters from cooled intracluster gas constitutes a previously unrecognised but prodigious source of GCs over cosmic timescales, and contributes to both their enormous numbers and broad colour dispersion in giant elliptical galaxies.

In order to investigate the origin of multiple stellar populations in the halo and bulge of the Milky Way, we have constructed chemical evolution models for the low-mass proto-Galactic subsystems such as globular clusters (GCs). Unlike previous studies, we assume that supernova blast waves undergo blowout without expelling the pre-enriched gas, while relatively slow winds of massive stars (WMS), together with the winds and ejecta from low and intermediate mass asymptotic-giant-branch stars (AGBs), are all locally retained in these less massive systems. We find that the observed Na-O anti-correlations in metal-poor GCs can be reproduced, when multiple episodes of starbursts are allowed to continue in these subsystems. A specific form of star formation history (SFH) with decreasing time intervals between the stellar generations, however, is required to obtain this result, which is in good agreement with the parameters obtained from our stellar evolution models for the horizontal-branch. The “mass budget problem” is also much alleviated by our models without ad-hoc assumptions on star formation efficiency (SFE) and initial mass function (IMF). We also applied these models to investigate the origin of super-helium-rich red clump stars in the metal-rich bulge as recently suggested by Lee et al. (2015). We find that chemical enrichments by the WMS can naturally reproduce the required helium enhancement (ΔY/ΔZ = 6) for the second generation stars. Disruption of proto-GCs in a hierarchical merging paradigm would have provided helium enhanced stars to the bulge field.

We investigated the stellar density substructures around four metal-poor globular clusters (NGC 6266, NGC 6626, NGC 6642, and NGC 6723) in the Galactic bulge. Wide-field near-infrared (JHKs) imaging data were obtained from WFCAM of UKIRT telescope. Field stars contamination around the globular clusters was reduced by using a statistical weighted filtering algorithm. Tidal stripping stellar substructures in the form of tidal tail (NGC 6266 and NGC 6626) or small density lobes/chunk (NGC 6642 and NGC 6723) were found around the four globular clusters in the two-dimensional density contour maps. We also find the overdensity features, which deviate from the theoretical models, in the outer region of radial density profiles. The observed results imply that the four globular clusters have experienced a strong tidal force or the bulge/disk shock effect of the Galaxy.