Star Formation in Violent and Normal Evolutionary Phases

Mergers of massive gas-rich galaxies trigger violent starbursts that - over timescales of $> 100$ Myr and regions $> 10$ kpc - form massive and compact star clusters comparable in mass and radii to Galactic globular clusters. The star formation efficiency is higher by 1 - 2 orders of magnitude in these bursts than in undisturbed spirals, irregulars or even BCDs. We ask the question if star formation in these extreme regimes is just a scaled-up version of the normal star formation mode of if the formation of globular clusters reveals fundamentally different conditions.Comment: 4 pages To appear in The Evolution of Galaxies. II. Basic building blocks, eds. M. Sauvage, G. Stasinska, L. Vigroux, D. Schaerer, S. Madde

Evolution of globular cluster systems in elliptical galaxies. I. Log-normal initial mass function

We study the evolution of globular cluster systems (GCS) in elliptical galaxies and explore the dependence of their main properties on the mass and the size of the host galaxy.The dependence of the evolution of the GCS mass function (GCMF), of the fraction of surviving clusters and of the ratio of the final to initial mass in clusters on the structure of the host galaxy as well as their variation with the galactocentric distance inside individual host galaxies has been thoroughly investigated.After a survey over a large number of different host galaxies we have restricted our attention to a sample of galaxies with effective masses and radii equal to those observed for dwarf,normal and giant ellipticals. We show that, in spite of large differences in the fraction of surviving clusters, the final mean masses of the GCMF in massive galaxies are very similar to each other with a small galaxy-to-galaxy dispersion;low-mass compact galaxies tend to have smaller values of the final mean mass and a larger galaxy-to-galaxy dispersion. These findings are in agreement with those of recent observational analyses. The fraction of surviving clusters increases with the mass of the host galaxy. We show that a small difference between the initial and the final mean mass and dispersion of the GCMF and the lack of a significant radial dependence of the mean mass inside individual galaxies do not necessarily imply that evolutionary processes have been unimportant in the evolution of the initial population of clusters. For giant galaxies most disruption occurs within the effective radius while for low-mass galaxies a significant disruption of clusters takes place also at larger galactocentric distances. The dependence of the results obtained on the initial mean mass of the GCMF is investigated. (abridged)Comment: 17 pages, accepted for publication in MNRA

Gravothermal Catastrophe in Anisotropic Spherical Systems

In this paper we investigate the gravothermal instability of spherical stellar systems endowed with a radially anisotropic velocity distribution. We focus our attention on the effects of anisotropy on the conditions for the onset of the instability and in particular we study the dependence of the spatial structure of critical models on the amount of anisotropy present in a system. The investigation has been carried out by the method of linear series which has already been used in the past to study the gravothermal instability of isotropic systems. We consider models described by King, Wilson and Woolley-Dickens distribution functions. In the case of King and Woolley-Dickens models, our results show that, for quite a wide range of amount of anisotropy in the system, the critical value of the concentration of the system (defined as the ratio of the tidal to the King core radius of the system) is approximately constant and equal to the corresponding value for isotropic systems. Only for very anisotropic systems the critical value of the concentration starts to change and it decreases significantly as the anisotropy increases and penetrates the inner parts of the system. For Wilson models the decrease of the concentration of critical models is preceded by an intermediate regime in which critical concentration increases, it reaches a maximum and then it starts to decrease. The critical value of the central potential always decreases as the anisotropy increases.Comment: 7pages, 5figures, to appear in MNRAS (figures have been replaced with their corrected versions

Dynamical evolution of the mass function and radial profile of the Galactic globular cluster system

Evolution of the mass function (MF) and radial distribution (RD) of the Galactic globular cluster (GC) system is calculated using an advanced and a realistic Fokker-Planck (FP) model that considers dynamical friction, disc/bulge shocks and eccentric cluster orbits. We perform hundreds of FP calculations with different initial cluster conditions, and then search a wide-parameter space for the best-fitting initial GC MF and RD that evolves into the observed present-day Galactic GC MF and RD. By allowing both MF and RD of the initial GC system to vary, which is attempted for the first time in the present Letter, we find that our best-fitting models have a higher peak mass for a lognormal initial MF and a higher cut-off mass for a power-law initial MF than previous estimates, but our initial total masses in GCs, M_{T,i} = 1.5-1.8x10^8 Msun, are comparable to previous results. Significant findings include that our best-fitting lognormal MF shifts downward by 0.35 dex during the period of 13 Gyr, and that our power-law initial MF models well-fit the observed MF and RD only when the initial MF is truncated at >~10^5 Msun. We also find that our results are insensitive to the initial distribution of orbit eccentricity and inclination, but are rather sensitive to the initial concentration of the clusters and to how the initial tidal radius is defined. If the clusters are assumed to be formed at the apocentre while filling the tidal radius there, M_{T,i} can be as high as 6.9x10^8 Msun, which amounts to ~75 per cent of the current mass in the stellar halo.Comment: To appear in May 2008 issue of MNRAS, 386, L6

Modeling the dynamical evolution of the M87 globular cluster system

We study the dynamical evolution of the M87 globular cluster system (GCS) with a number of numerical simulations. We explore a range of different initial conditions for the GCS mass function (GCMF), for the GCS spatial distribution and for the GCS velocity distribution. We confirm that an initial power-law GCMF like that observed in young cluster systems can be readily transformed through dynamical processes into a bell-shaped GCMF. However,only models with initial velocity distributions characterized by a strong radial anisotropy increasing with the galactocentric distance are able to reproduce the observed constancy of the GCMF at all radii.We show that such strongly radial orbital distributions are inconsistent with the observed kinematics of the M87 GCS. The evolution of models with a bell-shaped GCMF with a turnover similar to that currently observed in old GCS is also investigated. We show that models with this initial GCMF can satisfy all the observational constraints currently available on the GCS spatial distribution,the GCS velocity distribution and on the GCMF properties.In particular these models successfully reproduce both the lack of a radial gradient of the GCS mean mass recently found in an analysis of HST images of M87 at multiple locations, and the observed kinematics of the M87 GCS.Our simulations also show that evolutionary processes significantly affect the initial GCS properties by leading to the disruption of many clusters and changing the masses of those which survive.The preferential disruption of inner clusters flattens the initial GCS number density profile and it can explain the rising specific frequency with radius; we show that the inner flattening observed in the M87 GCS spatial distribution can be the result of the effects of dynamical evolution on an initially steep density profile. (abridged)Comment: 15 pages,14 figures;accepted for publication in The Astrophysical Journa

The incidence of binaries in Globular Cluster stellar populations

Binary fraction and orbital characteristics provide indications on the conditions of star formation, as they shed light on the environment they were born in. Multiple systems are more common in low density environments rather than in higher density ones. In the current debate about the formation of Globular Clusters and their multiple populations, studying the binary incidence in the populations they host offers a crucial piece of information on the environment of their birth and their subsequent dynamical evolution. Through a multi-year observational campaign using FLAMES at VLT, we monitored the radial velocity of 968 Red-Giant Branch stars located around the half-light radii in a sample of 10 Galactic Globular Clusters. We found a total of 21 radial velocity variables identified as {\it bona fide} binary stars, for a binary fraction of 2.2%$\pm$0.5%. When separating the sample into first generation and second generation stars, we find a binary fraction of 4.9%$\pm$1.3% and 1.2%$\pm$0.4% respectively. Through simulations that take into account possible sources of bias in detecting radial velocity variations in the two populations, we show that the difference is significant and only marginally affected by such effects. Such a different binary fraction strongly suggests different conditions in the environment of formation and evolution of first and second generations stars, with the latter being born in a much denser environment. Our result hence strongly supports the idea that the second generation forms in a dense subsystem at the center of the loosely distributed first generation, where (loose) binaries are efficiently destroyed.Comment: A&A, Accepte

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.Comment: 5 pages 4 figures. Accepted for publication in MNRAS Letter

The "UV-route" to search for Blue Straggler Stars in Globular Clusters: first results from the HST UV Legacy Survey

We used data from the HST UV Legacy Survey of Galactic Globular Clusters to select the Blue Straggler Star (BSS) population in four intermediate/high density systems (namely NGC 2808, NGC 6388, NGC 6541 and NGC 7078) through a "UV-guided search". This procedure consists in using the F275W images in each cluster to construct the master list of detected sources, and then force it to the images acquired in the other filters. Such an approach optimizes the detection of relatively hot stars and allows the detection of complete sample of BSSs even in the central region of high-density clusters, because the light from the bright cool giants, which dominates the optical emission in old stellar systems, is sensibly reduced at UV wavelengths. Our UV-guided selections of BSSs have been compared to the samples obtained in previous, optical-driven surveys, clearly demonstrating the efficiency of the UV approach. In each cluster we also measured the parameter A+, defined as the area enclosed between the cumulative radial distribution of BSSs and that of a reference population, which traces the level of BSS central segregation and the level of dynamical evolution suffered by the system. The values measured for the four clusters studied in this paper nicely fall along the dynamical sequence recently presented for a sample of 25 clusters.Comment: 30 pages, 9 figures, 4 tables; accepted for publication in Ap