Source author record

Enrico Vesperini

Enrico Vesperini appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2022arXiv

Early dynamical evolution of rotating star clusters in a tidal field

In order to explore how the early internal rotational properties of star clusters are affected by the external potential of their host galaxies, we have run a suite of $N$-body simulations following the early dynamical evolution and violent relaxation of rotating star clusters embedded in a tidal field. Our study focuses on models for which the cluster's rotation axis has a generic orientation relative to the torque of the tidal field. The interaction between the violent relaxation process, angular momentum of the cluster, and the external torque creates a complex kinematic structure within the cluster, most prominently a radial variation in the position of the rotation axis along both the polar and azimuthal directions. We also examine the cluster's velocity dispersion anisotropy and show that the projected anisotropy may be affected by the variation of the rotation axis directions within the cluster; the combination of projection effects and the complex kinematical features may result in the measurement of tangential anisotropy in the cluster's inner regions. We also characterize the structural properties of our clusters as a function of their initial rotation and virial ratio and find that clusters may develop a triaxial morphology and a radial variation of the minor axis not necessarily aligned with the rotation axis. Finally, we examine the long-term evolution of these complex kinematic features.

preprint2022arXiv

HST observations of the globular cluster NGC 6402 (M14) and its peculiar multiple populations

We present Hubble Space Telescope (HST) photometric results for NGC 6402, a highly reddened very luminous Galactic globular cluster (GC). Recent spectroscopic observations of its red giant stars have shown a quite peculiar behavior in the chemistry of its multiple populations. These results have prompted UV and optical HST observations aimed at obtaining the cluster's "Chromosome map" (ChM), an efficient tool to classify GCs and characterize their multiple populations. We find that the discontinuity in the abundance distributions of O, Mg, Al and Na inferred from spectroscopy is more nuanced in the ChM, which is mostly sensitive to nitrogen. Nevertheless, photometry in optical bands reveals a double main sequence, indicating a discontinuity in the helium content of the populations. The population with the largest chemical anomalies (extreme) peaks at a helium mass fraction Y~0.31. This helium content is consistent with results from the analysis of the distribution of horizontal-branch stars and the spectrophotometry of the red giants. The ChM and the color magnitude diagrams are compared with those in NGC 2808, a prototype GC with helium abundances up to Y > 0.35, and both confirm that NGC 6402 does not host stellar populations with such extreme helium content. Further, the ChM reveals the presence of a group of stars with larger metallicity, thus indicating that NGC 6402 is a Type II cluster. The modalities of formation of the multiple populations in NGC 6402 are briefly surveyed, with main attention on the Asymptotic Giant Branch and Supermassive star models, and on possible clusters' merging.

preprint2022arXiv

Long-term evolution of multimass rotating star clusters

We investigate the long-term dynamical evolution of the internal kinematics of multimass rotating star clusters. We have performed a set of N-body simulations to follow the internal evolution of clusters with different degrees of initial rotation and have explored the evolution of the rotational velocity, the degree of energy equipartition, and anisotropy in the velocity distribution. Our simulations show that: 1) as the cluster evolves, the rotational velocity develops a dependence on the stellar mass with more massive stars characterised by a more rapid rotation and a peak in the rotation curve closer to the cluster centre than low-mass stars; 2) the degree of energy equipartition in the cluster's intermediate and outer regions depends on the component of the velocity dispersion measured; for more rapidly rotating clusters, the evolution towards energy equipartition is more rapid in the direction of the rotational velocity; 3) the anisotropy in the velocity distribution is stronger for massive stars; 4) both the degree of mass segregation and energy equipartition are characterised by spatial anisotropy; they have a dependence on both $R$ and $z$, correlated with the flattening in the spatial variation of the cluster's density and velocity dispersion, as shown by 2D maps of the mass segregation and energy equipartition on the ($R$-$z$) meridional plane.

preprint2022arXiv

Mass segregation and dynamics of primordial binaries in star clusters with a radially anisotropic velocity distribution

This paper is the third in a series investigating, by means of $N$-body simulations, the implications of an initial radially anisotropic velocity distribution on the dynamics of star clusters. Such a velocity distribution may be imprinted during a cluster's early evolutionary stages and several observational studies have found examples of old globular clusters in which radial anisotropy is still present in the current velocity distribution. Here we focus on its influence on mass segregation and the dynamics of primordial binary stars (disruptions, ejections, and component exchanges). The larger fraction of stars on radial/highly eccentric orbits in the outer regions of anisotropic clusters lead to an enhancement in the dynamical interactions between inner and outer stars that affects both the process of mass segregation and the evolution of primordial binaries. The results of our simulations show that the time scale of mass segregation of the initially anisotropic cluster is longer in the core and shorter in the outer regions, when compared to the initially isotropic system. The evolution of primordial binaries is also significantly affected by the initial velocity distribution and we find that the rate of disruptions, ejections and exchange events affecting the primordial binaries in the anisotropic clusters is higher than in the isotropic ones.

preprint2022arXiv

The role of rotation on the formation of second generation stars in globular clusters

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.

preprint2020arXiv

Identifying Multiple Populations in M71 Using CN

We have observed the CN features at $\sim$3800 ~Å~ and 4120~ Å~ as well as the CH band at $\sim$4300~Å~ for 145 evolved stars in the Galactic globular cluster M71 using the multi-object spectrograph, Hydra, on the WIYN-3.5m telescope. We use these measurements to create two $δ$CN indices finding that both distributions are best fit by two populations, a CN-enhanced and CN-normal. We find that 42 $\pm$ 4\% of the RGB stars in our sample are CN-enhanced. The percentage of CN-enhanced is 40 $\pm$ 13\% for the AGB and 33 $\pm$ 9\% for the HB stars, which suggests there are no missing second generation stars at these stages of stellar evolution. The two generations also separate in magnitude and color on the HB, which allows us to find the difference in He abundance between the two populations by fitting appropriate ZAHBs. The broad range of distances from the cluster's center covered by our sample allows us to study the dependence of the ratio of the number of first to second population stars on the distance from the cluster's center, and we find that this ratio does not vary radially and that the two populations are spatially mixed. Finally, we compare our identification of multiple populations with the classification based on the Na-O anti-correlation and the HST UV photometry, and we find good agreement with both methods.

preprint2020arXiv

Substructure in the Globular Cluster Populations of the Virgo Cluster Elliptical Galaxies M84 and M86

We have carried out a search for substructure within the globular cluster systems of M84 (NGC 4374) and M86 (NGC 4406), two giant elliptical galaxies in the Virgo Cluster. We use wide-field (36 arcmin x 36 arcmin), multi-color broadband imaging to identify globular cluster candidates in these two galaxies as well as several other nearby lower-mass galaxies. Our analysis of the spatial locations of the globular cluster candidates reveals several substructures, including: a peak in the projected number density of globular clusters in M86 that is offset from the system center and may be at least partly due to the presence of the dwarf elliptical galaxy NGC 4406B; a bridge that connects the M84 and M86 globular cluster systems; and a boxy iso-density contour along the southeast side of the M86 globular cluster system. We divide our sample into red (metal-rich) and blue (metal-poor) globular cluster candidates to look for differences in the spatial distributions of the two populations and find that the blue cluster candidates are the dominant population in each of the substructures we identify. We also incorporate the measurements from two radial velocity surveys of the globular clusters in the region and find that the bridge substructure is populated by globular clusters with a mix of velocities that are consistent with either M86 and M84, possibly providing further evidence for interaction signatures between the two galaxies.

preprint2016arXiv

Dynamical Formation of Cataclysmic Variables in Globular Clusters

The formation and evolution of X-ray sources in globular clusters is likely to be affected by the cluster internal dynamics and the stellar interactions in the cluster dense environment.Several observational studies have revealed a correlation between the number of X-ray sources and the stellar encounter rate and provided evidence of the role of dynamics in the formation of X-ray binaries. We have performed a survey of Monte-Carlo simulations aimed at exploring the connection between the dynamics and formation of cataclysmic variables (CVs) and the origin of the observed correlation between the number of these objects, $N_{\rm cv}$, and the stellar encounter rate, $Γ$.The results of our simulations show a correlation between $N_{\rm cv}$ and $Γ$ as found in observational data, illustrate the essential role played by dynamics, and shed light on the dynamical history behind this correlation. CVs in our simulations are more centrally concentrated than single stars with masses close to those of turn-off stars, although this trend is stronger for CVs formed from primordial binaries undergoing exchange encounters, which include a population of more massive CVs absent in the group of CVs formed from binaries not suffering any component exchange.

preprint2016arXiv

Evolution of Binary Stars in Multiple-Population Globular Clusters - II. Compact Binaries

We present the results of a survey of N-body simulations aimed at exploring the evolution of compact binaries in multiple-population globular clusters.We show that as a consequence of the initial differences in the structural properties of the first-generation (FG) and the second-generation (SG) populations and the effects of dynamical processes on binary stars, the SG binary fraction decreases more rapidly than that of the FG population. The difference between the FG and SG binary fraction is qualitatively similar to but quantitatively smaller than that found for wider binaries in our previous investigations.The evolution of the radial variation of the binary fraction is driven by the interplay between binary segregation, ionization and ejection. Ionization and ejection counteract in part the effects of mass segregation but for compact binaries the effects of segregation dominate and the inner binary fraction increases during the cluster evolution. We explore the variation of the difference between the FG and the SG binary fraction with the distance from the cluster centre and its dependence on the binary binding energy and cluster structural parameters. The difference between the binary fraction in the FG and the SG populations found in our simulations is consistent with the results of observational studies finding a smaller binary fraction in the SG population.

preprint2016arXiv

Investigating mass segregation process in globular clusters with Blue Straggler Stars: the impact of dark remnants

We present the results of a set of N-body simulations aimed at exploring how the process of mass segregation (as traced by the spatial distribution of blue straggler stars, BSSs) is affected by the presence of a population of heavy dark remnants (as neutron stars and black holes). To this end, clusters characterized by different initial concentrations and different fractions of dark remnants have been modeled. We find that an increasing fraction of stellar-mass black holes significantly delays the mass segregation of BSSs and the visible stellar component. In order to trace the evolution of BSS segregation, we introduce a new parameter ($A^+$) that can be easily measured when the cumulative radial distribution of these stars and a reference population are available. Our simulations show that $A^+$ might also be used as an approximate indicator of the time remaining to the core collapse of the visible component.

preprint2016arXiv

Kinematical evolution of tidally limited star clusters: the role of retrograde stellar orbits

The presence of an external tidal field often induces significant dynamical evolutionary effects on the internal kinematics of star clusters. Previous studies investigating the restricted three-body problem with applications to star cluster dynamics have shown that unbound stars on retrograde orbits (with respect to the direction of the cluster's orbit) are more stable against escape than prograde orbits, and predicted that a star cluster might acquire retrograde rotation through preferential escape of stars on prograde orbits. In this study we present evidence of this prediction, but we also illustrate that there are additional effects that cannot be accounted for by the preferential escape of prograde orbits alone. Specifically, in the early evolution, initially underfilling models increase their fraction of retrograde stars without losing significant mass, and acquire a retrograde angular velocity. We attribute this effect to the development of preferentially eccentric/radial orbits in the outer regions of star clusters as they are expanding into their tidal limitation. We explore the implications of the evolution of the fraction of prograde and retrograde stars for the evolution of the cluster internal rotation, and its dependence on the initial structural properties. Although all the systems studied here evolve towards an approximately solid-body internal rotation with angular velocity equal to about half of the angular velocity of the cluster orbital motion around the host galaxy, the evolutionary history of the radial profile of the cluster internal angular velocity depends on the cluster initial structure.

preprint2016arXiv

On the Link Between Energy Equipartition and Radial Variation in the Stellar Mass Function of Star Clusters

We make use of $N$-body simulations to determine the relationship between two observable parameters that are used to quantify mass segregation and energy equipartition in star clusters. Mass segregation can be quantified by measuring how the slope of a cluster's stellar mass function $α$ changes with clustercentric distance r, and then calculating $δ_α= \frac{d α(r)}{d ln(r/r_m)}$ where $r_m$ is the cluster's half-mass radius. The degree of energy equipartition in a cluster is quantified by $η$, which is a measure of how stellar velocity dispersion $σ$ depends on stellar mass m via $σ(m) \propto m^{-η}$. Through a suite of $N$-body star cluster simulations with a range of initial sizes, binary fractions, orbits, black hole retention fractions, and initial mass functions, we present the co-evolution of $δ_α$ and $η$. We find that measurements of the global $η$ are strongly affected by the radial dependence of $σ$ and mean stellar mass and the relationship between $η$ and $δ_α$ depends mainly on the cluster's initial conditions and the tidal field. Within $r_m$, where these effects are minimized, we find that $η$ and $δ_α$ initially share a linear relationship. However, once the degree of mass segregation increases such that the radial dependence of $σ$ and mean stellar mass become a factor within $r_m$, or the cluster undergoes core collapse, the relationship breaks down. We propose a method for determining $η$ within $r_m$ from an observational measurement of $δ_α$. In cases where $η$ and $δ_α$ can be measured independently, this new method offers a way of measuring the cluster's dynamical state.

preprint2016arXiv

Radial Variation in the Stellar Mass Functions of Star Clusters

A number of recent observational studies of Galactic globular clusters have measured the variation in the slope of a cluster's stellar mass function $α$ with clustercentric distance $r$. In order to gather a deeper understanding of the information contained in such observations, we have explored the evolution of $α(r)$ for star clusters with a variety of initial conditions using a large suite of $N$-body simulations. We have specifically studied how the time evolution of $α(r)$ is affected by initial size, mass, binary fraction, primordial mass segregation, black hole retention, an external tidal field, and the initial mass function itself. Previous studies have shown that the evolution of $α_G$ is closely related to the amount of mass loss suffered by a cluster. Hence for each simulation we have also followed the evolution of the slope of the cluster's global stellar mass function, $α_G$, and have shown that clusters follow a well-defined track in the $α_G$-$dα(r)/d(ln(r/r_m))$ plane. The location of a cluster on the $α_G-dα(r)/d(ln(r/r_m))$ plane can therefore constrain its dynamical history and, in particular, constrain possible variations in the stellar initial mass function. The $α_G$-$dα(r)/d(ln(r/r_m))$ plane thus serves as a key tool for fully exploiting the information contained in wide field studies of cluster stellar mass functions.

preprint2016arXiv

The Hubble Space Telescope UV Legacy Survey of Galactic Globular Clusters. X. The radial distribution of stellar populations in NGC 2808

Due to their extreme helium abundance, the multiple stellar populations of the globular cluster NGC 2808 have been widely investigated from a photometric, spectroscopic, and kinematic perspective. The most striking feature of the color-magnitude diagram of NGC 2808 is the triple main sequence (MS), with the red MS corresponding to a stellar population with primordial helium, and the middle and the blue MS being enhanced in helium up to Y$\sim$0.32 and $\sim$0.38, respectively. A recent study has revealed that this massive cluster hosts at least five distinct stellar populations (A, B, C, D, and E). Among them populations A, B, and C correspond to the red MS, while populations C and D are connected to the middle and the blue MS. In this paper we exploit Hubble-Space-Telescope photometry to investigate the radial distribution of the red, the middle and the blue MS from the cluster center out to about 8.5 arcmin. Our analysis shows that the radial distribution of each of the three MSs is different. In particular, as predicted from multiple-population formation models, both the blue MS and the middle MS appears to be more concentrated than the red MS with a significance level for this result wich is above 3σ.

preprint2015arXiv

Chemical Abundances in NGC 5053: A Very Metal-Poor and Dynamically Complex Globular Cluster

NGC 5053 provides a rich environment to test our understanding of the complex evolution of globular clusters (GCs). Recent studies have found that this cluster has interesting morphological features beyond the typical spherical distribution of GCs, suggesting that external tidal effects have played an important role in its evolution and current properties. Additionally, simulations have shown that NGC 5053 could be a likely candidate to belong to the Sagittarius dwarf galaxy (Sgr dSph) stream. Using the Wisconsin-Indiana-Yale-NOAO-Hydra multi-object spectrograph, we have collected high quality (signal-to-noise ratio $\sim$ 75-90), medium-resolution spectra for red giant branch stars in NGC 5053. Using these spectra we have measured the Fe, Ca, Ti, Ni, Ba, Na, and O abundances in the cluster. We measure an average cluster [Fe/H] abundance of -2.45 with a standard deviation of 0.04 dex, making NGC 5053 one of the most metal-poor GCs in the Milky Way (MW). The [Ca/Fe], [Ti/Fe], and [Ba/Fe] we measure are consistent with the abundances of MW halo stars at a similar metallicity, with alpha-enhanced ratios and slightly depleted [Ba/Fe]. The Na and O abundances show the Na-O anti-correlation found in most GCs. From our abundance analysis it appears that NGC 5053 is at least chemically similar to other GCs found in the MW. This does not, however, rule out NGC 5053 being associated with the Sgr dSph stream.

preprint2015arXiv

Evolution of Binary Stars in Multiple-Population Globular Clusters

The discovery of multiple stellar populations in globular clusters has implications for all the aspects of the study of these stellar systems. In this paper, by means of N-body simulations, we study the evolution of binary stars in multiple-population clusters and explore the implications of the initial differences in the spatial distribution of different stellar populations for the evolution and survival of their binary stars. Our simulations show that initial differences between the spatial distribution of first-generation (FG) and second-generation (SG) stars can leave a fingerprint in the current properties of the binary population. SG binaries are disrupted more efficiently than those of the FG population resulting in a global SG binary fraction smaller than that of the FG. As for surviving binaries, dynamical evolution produces a difference between the SG and the FG binary binding energy distribution with the SG population characterized by a larger fraction of high binding energy (more bound) binaries. We have also studied the dependence of the binary properties on the distance from the cluster centre. Although the global binary fraction decreases more rapidly for the SG population, the local binary fraction measured in the cluster inner regions may still be dominated by SG binaries. The extent of the differences between the surviving FG and SG binary binding energy distribution also varies radially within the cluster and is larger in the cluster inner regions.

preprint2015arXiv

Light-Element Abundances of Giant Stars in the Globular Cluster M71 (NGC 6838)

Aluminum is the heaviest light element displaying large star--to--star variations in Galactic globular clusters (GCs). This element may provide additional insight into the origin of the multiple populations now known to be a common place in GCs, and also the nature of the first-generation stars responsible for a cluster's chemical inhomogeneities. In a previous analysis we found that, unlike more metal-poor GCs, 47 Tuc did not exhibit a strong Na-Al correlation, which motivates a careful study of the similar metallicity but less massive GC M71. We present chemical abundances of O, Na, Al, and Fe for 33 giants in M71 using spectra obtained with the WIYN-Hydra spectrograph. Our spectroscopic analysis finds that, similar to 47 Tuc and in contrast with more metal-poor GCs, M71 stars do not exhibit a strong Na-Al correlation and span a relatively narrow range in [Al/Fe], which are characteristics that GC formation models must reproduce.

preprint2015arXiv

Rapidly rotating second-generation progenitors for the blue hook stars of ω Cen

Horizontal Branch stars belong to an advanced stage in the evolution of the oldest stellar galactic population, occurring either as field halo stars or grouped in globular clusters. The discovery of multiple populations in these clusters, that were previously believed to have single populations gave rise to the currently accepted theory that the hottest horizontal branch members (the blue hook stars, which had late helium-core flash ignition, followed by deep mixing) are the progeny of a helium-rich "second generation" of stars. It is not known why such a supposedly rare event (a late flash followed by mixing) is so common that the blue hook of ω Cen contains \sim 30% of horizontal branch stars 10 , or why the blue hook luminosity range in this massive cluster cannot be reproduced by models. Here we report that the presence of helium core masses up to \sim 0.04 solar masses larger than the core mass resulting from evolution is required to solve the luminosity range problem. We model this by taking into account the dispersion in rotation rates achieved by the progenitors, whose premain sequence accretion disc suffered an early disruption in the dense environment of the cluster's central regions where second-generation stars form. Rotation may also account for frequent late-flash-mixing events in massive globular clusters.

preprint2015arXiv

Velocity anisotropy in tidally limited star clusters

We explore the long-term evolution of the anisotropy in the velocity space of star clusters starting with different structural and kinematical properties. We show that the evolution of the radial anisotropy strength and its radial variation within a cluster contain distinct imprints of the cluster initial structural properties, dynamical history, and of the external tidal field of its host galaxy. Initially isotropic and compact clusters with small initial values of the ratio of the half-mass to Jacobi radius, $r_h/r_J$, develop a strong radial anisotropy during their long-term dynamical evolution. Many clusters, if formed with small values of $r_h/r_J$, should now be characterized by a significant radial anisotropy increasing with the distance from the cluster centre, reaching its maximum at a distance between 0.2 $r_J$ and 0.4 $r_J$, and then becoming more isotropic or mildly tangentially anisotropic in the outermost regions. A similar radial variation of the anisotropy can also result from an early violent relaxation phase. In both cases, as a cluster continues its evolution and loses mass, the anisotropy eventually starts to decrease and the system evolves toward an isotropic velocity distribution. However, in order to completely erase the strong anisotropy developed by these compact systems during their evolution, they must be in the advanced stages of their evolution and lose a large fraction of their initial mass. Clusters that are initially isotropic and characterized by larger initial values of $r_h/r_J$, on the other hand, never develop a significant radial anisotropy.

preprint2014arXiv

Kinematical fingerprints of star cluster early dynamical evolution

We study the effects of the external tidal field on the violent relaxation phase of star clusters dynamical evolution, with particular attention to the kinematical properties of the equilibrium configurations emerging at the end of this phase.We show that star clusters undergoing the process of violent relaxation in the tidal field of their host galaxy can acquire significant internal differential rotation and are characterized by a distinctive radial variation of the velocity anisotropy. These kinematical properties are the result of the symmetry breaking introduced by the external tidal field in the collapse phase and of the action of the Coriolis force on the orbit of the stars. The resulting equilibrium configurations are characterized by differential rotation, with a peak located between one and two half-mass radii. As for the anisotropy, similar to clusters evolving in isolation, the systems explored in this Letter are characterized by an inner isotropic core, followed by a region of increasing radial anisotropy. However for systems evolving in an external tidal field the degree of radial anisotropy reaches a maximum in the cluster intermediate regions and then progressively decreases, with the cluster outermost regions being characterized by isotropy or a mild tangential anisotropy. Young or old but less-relaxed dynamically young star clusters may keep memory of these kinematical fingerprints of their early dynamical evolution.

preprint2014arXiv

Pre-main sequence accretion and the formation of multiple populations in Globular Clusters

We investigate the viability of a model in which the chemical anomalies among Globular Cluster stars are due to accretion of gas onto the protostellar discs of low mass stars. This model has been suggested as a way to reduce the large initial cluster masses required by other models for the formation of multiple stellar generations. We numerically follow the evolution of the accreting stars, and we show that the structure of the seed star does not remain fully convective for the whole duration of the accretion phase. Stellar populations showing discrete abundances of helium in the core, that seem to be present in some clusters, might be formed with this mechanism only if accretion occurs before the core of the stars become radiative (within 2-3Myr) or if a thermohaline instability is triggered, to achieve full mixing after the accretion phase ends. We also show that the lithium abundances in accreted structures may vary by orders of magnitude in equal masses obtained by accreting different masses. In addition, the same thermohaline mixing which could provide a homogeneous helium distribution down to the stellar center, would destroy any lithium surviving in the envelope, so that both helium homogeneity and lithium survival require that the accretion phase be limited to the first couple of million years of the cluster evolution. Such a short accretion phase strongly reduces the amount of processed matter available, and reintroduces the requirement of an extremely large initial mass for the protocluster.

preprint2013arXiv

Simulating Star Clusters with the AMUSE Software Framework: I. Dependence of Cluster Lifetimes on Model Assumptions and Cluster Dissolution Modes

We perform a series of simulations of evolving star clusters using AMUSE (the Astrophysical Multipurpose Software Environment), a new community-based multi-physics simulation package, and compare our results to existing work. These simulations model a star cluster beginning with a King model distribution and a selection of power-law initial mass functions, and contain a tidal cut-off. They are evolved using collisional stellar dynamics and include mass loss due to stellar evolution. After determining that the differences between AMUSE results and prior publications are understood, we explored the variation in cluster lifetimes due to the random realization noise introduced by transforming a King model to specific initial conditions. This random realization noise can affect the lifetime of a simulated star cluster by up to 30%. Two modes of star cluster dissolution were identified: a mass evolution curve that contains a run-away cluster dissolution with a sudden loss of mass, and a dissolution mode that does not contain this feature. We refer to these dissolution modes as "dynamical" and "relaxation" dominated respectively. For Salpeter-like initial mass functions, we determined the boundary between these two modes in terms of the dynamical and relaxation time scales.

preprint2013arXiv

The MICA Experiment: Astrophysics in Virtual Worlds

We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. We conducted scientific collaboration meetings, professional seminars, a workshop, classroom instruction, public lectures, informal discussions and gatherings, and experiments in immersive, interactive visualization of high-dimensional scientific data. Perhaps the most successful of these was our program of popular science lectures, illustrating yet again the great potential of immersive VR as an educational and outreach platform. While the members of our research groups and some collaborators found the use of immersive VR as a professional telepresence tool to be very effective, we did not convince a broader astrophysics community to adopt it at this time, despite some efforts; we discuss some possible reasons for this non-uptake. On the whole, we conclude that immersive VR has a great potential as a scientific and educational platform, as the technology matures and becomes more broadly available and accepted.

preprint2012arXiv

Rapid Mass Segregation in Massive Star Clusters

Several dynamical scenarios have been proposed that can lead to prompt mass segregation on the crossing time scale of a young cluster. They generally rely on cool and/or clumpy initial conditions, and are most relevant to small systems. As a counterpoint, we present a novel dynamical mechanism that can operate in relatively large, homogeneous, cool or cold systems. This mechanism may be important in understanding the assembly of large mass-segregated clusters from smaller clumps.

preprint2012arXiv

Space-Based UV/Optical Wide-Field Imaging and Spectroscopy: Near-Field Cosmology and Galaxy Evolution Using Globular Clusters in Nearby Galaxies

Star formation plays a central role in the evolution of galaxies and of the Universe as a whole. Studies of star-forming regions in the local universe have shown that star formation typically occurs in a clustered fashion. Building a coherent picture of how star clusters form and evolve is therefore critical to our overall understanding of the star formation process. Most clusters disrupt after they form, thus contributing to the field star population. However, the most massive and dense clusters remain bound and survive for a Hubble time. These globular clusters provide unique observational probes of the formation history of their host galaxies. In particular, the age and metallicity can be determined for each globular cluster individually, allowing the distribution of ages and metallicities within host galaxies to be constrained. We show how space-based UV-to-near-IR imaging covering a wide field of view (>= 20 arcmin per axis) and deep UV/Optical multi-object spectroscopy of globular cluster systems in nearby galaxies would allow one to place important new constraints on the formation history of early-type galaxies and their structural subcomponents (e.g., bulge, halo).

preprint2012arXiv

Super AGB and AGB evolution and the chemical inventory in NGC 2419

We follow the scenario of formation of second generation stars in globular clusters by matter processed by hot bottom burning (HBB) in massive asymptotic giant branch (AGB) stars and Super AGB stars (SAGB). In the cluster NGC 2419 we assume the presence of an extreme population directly formed from the AGB and SAGB ejecta, so we can directly compare the yields for a metallicity Z=0.0003 with the chemical inventory of the cluster NGC 2419. At such a low metallicity, the HBB temperatures (well above 108K) allow a very advanced nucleosynthesis. Masses of about 6Mo deplete Mg and synthesize Si, going beyond Al, so this latter element results only moderately enhanced; sodium can not be enhanced. The models are consistent with the observations, although the predicted Mg depletion is not as strong as in the observed stars. We predict that the oxygen abundance must be depleted by a huge factor in the Mg poor stars. The HBB temperatures are close to the region where other p capture reactions on heavier nuclei become possible. We show that high potassium abundance found in Mg poor stars can be achieved during HBB, by p captures on the argon nuclei, if the relevant cross section(s) are larger than listed in the literature or if the HBB temperature is higher. Finally, we speculate that some calcium production is occurring owing to proton capture on potassium. We emphasize the importance of a strong effort to measure a larger sample of abundances in this cluster.

preprint2012arXiv

The role of super-asymptotic giant branch ejecta in the abundance patterns of multiple populations in globular clusters

In order to account for the chemical composition of a stellar second generation (SG), Globular Clusters (GCs) evolution models based on the asymptotic giant branch (AGB) scenario so far included only the yields available for the massive AGB stars, while the possible role of super-AGB ejecta was either extrapolated or not considered. In this work, we explore the role of super-AGB ejecta using yields recently calculated by Ventura and D'Antona. Models of clusters showing extended Na-O anticorrelations, like NGC 2808, indicate that a SG formation history similar to that outlined in our previous work is required: formation of an Extreme population with very large helium content from the pure ejecta of super-AGB stars, followed by formation of an Intermediate population by dilution of stellar ejecta with pristine gas. The very O-poor Na-rich Extreme stars can be accounted for once deep-mixing is assumed in SG giants forming in a gas with helium abundance Y> 0.34, which significantly reduces the atmospheric oxygen content, while preserving the sodium abundance. On the other hand, for clusters showing a mild O-Na anticorrelation, like M 4, the use of the new yields broadens the range of SG formation routes leading to abundance patterns consistent with observations. It is shown that models in which SG stars form only from super-AGB ejecta promptly diluted with pristine gas can reproduce the observations. We discuss the variety of small helium variations occurring in this model and its relevance for the horizontal branch morphology. In some of these models the duration of the SG formation episode can be as short as \sim10 Myr; the formation time of the SG is thus compatible with the survival of a cooling flow in the GC, previous to the explosion of the SG core collapse supernovae. We also explore models with formation of multiple populations in individual bursts, each lasting no longer than \sim10 Myr.

preprint2011arXiv

Binary Star Disruption in Globular Clusters with Multiple Stellar Populations

The discovery of multiple stellar populations in globular clusters raises fundamental questions concerning the formation and dynamical history of these systems. In a previous study aimed at exploring the formation of second-generation (SG) stars from the ejecta of first-generation (FG) AGB stars, and the subsequent dynamical evolution of the cluster, we showed that SG stars are expected to form in a dense subsystem concentrated in the inner regions of the FG cluster. In this paper we explore the implications of the structural properties of multiple-population clusters, and in particular the presence of the inner SG subsystem, for the disruption of binary stars. We quantify the enhancement of the binary disruption rate due to the presence of the central SG subsystem for a number of different initial conditions. Our calculations show that SG binaries, which are assumed to be more concentrated in the cluster inner regions, are disrupted at a substantially larger rate than FG binaries. Assuming a similar initial fraction of FG and SG binaries, our dynamical study indicates that the SG population is now expected to contain a significantly smaller binary fraction than the FG population.

preprint2011arXiv

Formation of Multiple Populations in Globular Clusters: constraints on the dilution by pristine gas

The star-to-star differences in the abundance of light elements observed in the globular clusters (GCs) can be explained assuming that a second generation (SG) of stars form in the gas ejected by the asymptotic giant branch (AGB) stars belonging to a first stellar generation. However, while Na and O appear to be anticorrelated in the cluster stars, from the stellar models they turn out to be correlated into the AGB ejecta. In order to reconcile the stellar theory with the observational findings, all the GC models invoke an early dilution of AGB ejecta with pristine gas occurring during the SG formation. Despite a vast consensus about the occurrence of such dilution, the physical process behind it is still unknown. In the present paper we set some general constraints on the pristine gas dynamics and on the possible amount of pristine gas involved in the SG formation, making use of a one zone chemical model. We find that such a dilution is a necessary ingredient in the SG star formation to explain the observed abundance patterns. We confirm the conclusion of our previous works showing that clusters must have been initially much more massive. We also show that models assuming that clusters had an initial mass similar to their current one, and adopting a large fraction of pristine gas to form SG stars, fail to reproduce the observed Na-O anticorrelation and are not viable. We finally show that the dilution event should be restricted in time, rather than extended for the all duration of the SG formation.

preprint2010arXiv

Abundance patterns of multiple populations in Globular Clusters: a chemical evolution model based on yields from AGB ejecta

A large number of spectroscopic studies have provided evidence of the presence of multiple populations in globular clusters by revealing patterns in the stellar chemical abundances. This paper is aimed at studying the origin of these abundance patterns. We explore a model in which second generation (SG) stars form out of a mix of pristine gas and ejecta of the first generation of asymptotic giant branch stars. We first study the constraints imposed by the spectroscopic data of SG stars in globular clusters on the chemical properties of the asymptotic and super asymptotic giant branch ejecta. With a simple one-zone chemical model, we then explore the formation of the SG population abundance patterns focussing our attention on the Na-O, Al-Mg anticorrelations and on the helium distribution function. We carry out a survey of models and explore the dependence of the final SG chemical properties on the key parameters affecting the gas dynamics and the SG formation process. Finally, we use our chemical evolution framework to build specific models for NGC 2808 and M4, two Galactic globular clusters which show different patterns in the Na-O and Mg-Al anticorrelation and have different helium distributions. We find that the amount of pristine gas involved in the formation of SG stars is a key parameter to fit the observed O-Na and Mg-Al patterns. The helium distribution function for these models is in general good agreement with the observed one. Our models, by shedding light on the role of different parameters and their interplay in determining the final SG chemical properties, illustrate the basic ingredients, constraints and problems encountered in this self-enrichment scenario which must be addressed by more sophisticated chemical and hydrodynamic simulations.

preprint2010arXiv

Intermediate-Mass Black Holes in Early Globular Clusters

Spectroscopic and photometric observations show that many globular clusters host multiple stellar populations, challenging the common paradigm that globular clusters are "simple stellar populations" composed of stars of uniform age and chemical composition. The chemical abundances of second-generation (SG) stars constrain the sources of gas out of which these stars must have formed, indicating that the gas must contain matter processed through the high-temperature CNO cycle. First-generation massive Asymptotic Giant Branch (AGB) stars have been proposed as the source of this gas. In a previous study, by means of hydrodynamical and N-body simulations, we have shown that the AGB ejecta collect in a cooling flow in the cluster core, where the gas reaches high densities, ultimately forming a centrally concentrated subsystem of SG stars. In this Letter we show that the high gas density can also lead to significant accretion onto a pre-existing seed black hole. We show that gas accretion can increase the black hole mass by up to a factor of 100. The details of the gas dynamics are important in determining the actual black hole growth. Assuming a near-universal seed black hole mass and small cluster-to-cluster variations in the duration of the SG formation phase, the outcome of our scenario is one in which the present intermediate-mass black hole (IMBH) mass may have only a weak dependence on the current cluster properties. The scenario presented provides a natural mechanism for the formation of an IMBH at the cluster center during the SG star-formation phase.

preprint2010arXiv

Widespread presence of shallow cusps in the surface-brightness profile of globular clusters

Surface brightness profiles of globular clusters with shallow central cusps (Sigma ~ R^v with -0.3<~ v <~ -0.05) have been associated by several recent studies with the presence of a central intermediate mass black hole (IMBH). Such shallow slopes are observed in several globular clusters thanks to the high angular resolution of Hubble Space Telescope imaging. In this Letter we evaluate whether shallow cusps are a unique signature of a central IMBH by analyzing a sample of direct N-body simulations of star clusters with and without a central IMBH. We ``observe'' the simulations as if they were HST images. Shallow cusps are common in our simulation sample: star clusters without an IMBH have v >~ -0.3 in the pre-core-collapse and core-collapse phases. Post-core-collapse clusters without an IMBH transition to steeper cusps, -0.7<~ v <~ -0.4, only if the primordial binary fraction is very small, f_{bin}< 3 per cent, and if there are few stellar-mass black holes remaining. Otherwise v values overlap the range usually ascribed to the presence of an IMBH throughout the entire duration of the simulations. In addition, measuring v is intrinsically prone to significant uncertainty, therefore typical measurement errors may lead to v > -0.3 even when <v> <~ -0.4. Overall our analysis shows that a shallow cusp is not an unequivocal signature of a central IMBH and casts serious doubts on the usefulness of measuring v in the context of the hunt for IMBHs in globular clusters.

preprint2009arXiv

Effects of Primordial Mass Segregation on the Dynamical Evolution of Star Clusters

In this paper we use N-body simulations to study the effects of primordial mass segregation on the early and long-term evolution of star clusters. Our simulations show that in segregated clusters early mass loss due to stellar evolution triggers a stronger expansion than for unsegregated clusters. Tidally limited, strongly segregated clusters may dissolve rapidly as a consequence of this early expansion, while segregated clusters initially underfilling their Roche lobe can survive the early expansion and have a lifetime similar to that of unsegregated clusters. Long-lived initially segregated clusters tend to have looser structure and reach core collapse later in their evolution than initially unsegregated clusters. We have also compared the effects of dynamical evolution on the global stellar mass function (MF) of low-mass main sequence stars. In all cases the MF flattens as the cluster loses stars. The amount of MF flattening induced by a given amount of mass loss in a rapidly dissolving initially segregated cluster is less than for an unsegregated cluster. The evolution of the MF of a long-lived segregated cluster, on the other hand, is very similar to that of an initially unsegregated cluster.

preprint2009arXiv

Evidence for two populations of Galactic globular clusters from the ratio of their half-mass to Jacobi radii

We investigate the ratio between the half-mass radii r_h of Galactic globular clusters and their Jacobi radii r_J given by the potential of the Milky Way and show that clusters with galactocentric distances R_{GC}>8 kpc fall into two distinct groups: one group of compact, tidally-underfilling clusters with r_h/r_J<0.05 and another group of tidally filling clusters which have 0.1 < r_h/r_J<0.3. We find no correlation between the membership of a particular cluster to one of these groups and its membership in the old or younger halo population. Based on the relaxation times and orbits of the clusters, we argue that compact clusters and most clusters in the inner Milky Way were born compact with half-mass radii r_h < 1 pc. Some of the tidally-filling clusters might have formed compact as well, but the majority likely formed with large half-mass radii. Galactic globular clusters therefore show a similar dichotomy as was recently found for globular clusters in dwarf galaxies and for young star clusters in the Milky Way. It seems likely that some of the tidally-filling clusters are evolving along the main sequence line of clusters recently discovered by Kuepper et al. (2008) and are in the process of dissolution.

preprint2009arXiv

Star cluster dynamics

Dynamical evolution plays a key role in shaping the current properties of star clusters and star cluster systems. A detailed understanding of the effects of evolutionary processes is essential to be able to disentangle the properties which result from dynamical evolution from those imprinted at the time of cluster formation. In this review, we focus our attention on globular clusters and review the main physical ingredients driving their early and long-term evolution, describe the possible evolutionary routes and show how cluster structure and stellar content are affected by dynamical evolution.

preprint2009arXiv

Tidal disruption, global mass function and structural parameters evolution in star clusters

[abridged] We present a unified picture for the evolution of star clusters on the two-body relaxation timescale. We use direct N-body simulations of star clusters in a galactic tidal field starting from different multi-mass King models, up to 10% of primordial binaries and up to Ntot=65536 particles. An additional run also includes a central Intermediate Mass Black Hole. We find that for the broad range of initial conditions we have studied the stellar mass function of these systems presents a universal evolution which depends only on the fractional mass loss. The structure of the system, as measured by the core to half mass radius ratio, also evolves toward a universal state, which is set by the efficiency of heating on the visible population of stars induced by dynamical interactions in the core of the system. Interactions with dark remnants are dominant over the heating induced by a moderate population of primordial binaries (3-5%), especially under the assumption that most of the neutron stars and black holes are retained in the system. All our models without primordial binaries undergo a deep gravothermal collapse in the radial mass profile. However their projected light distribution can be well fitted by medium concentration King models (with parameter W0 ~ 8), even though there tends to be an excess over the best fit for the innermost points of the surface brightness. This excess is consistent with a shallow cusp in the surface brightness (mu(R) ~ R^{-v} with v ~ 0.4-0.7), like it has been observed for many globular clusters from high-resolution HST imaging. Classification of core-collapsed globular clusters based on their surface brightness profile is likely to fail in systems that have already bounced back to lower concentrations.

preprint2006arXiv

A dynamical origin for early mass segregation in young star clusters

Some young star clusters show a degree of mass segregation that is inconsistent with the effects of standard two-body relaxation from an initially unsegregated system without substructure, in virial equilibrium, and it is unclear whether current cluster formation models can account for this degree of initial segregation in clusters of significant mass. In this Letter we demonstrate that mergers of small clumps that are either initially mass segregated, or in which mass segregation can be produced by two-body relaxation before they merge, generically lead to larger systems which inherit the progenitor clumps' segregation. We conclude that clusters formed in this way are naturally mass segregated, accounting for the anomalous observations and suggesting that this process of prompt mass segregation due to initial clumping should be taken fully into account in constructing cluster dynamical models.