The Star Cluster Population of M51: II. Age distribution and relations among the derived parameters

We use archival Hubble Space Telescope observations of broad-band images from the ultraviolet (F255W-filter) through the near infrared (NICMOS F160W-filter) to study the star cluster population of the interacting spiral galaxy M 51. We obtain age, mass, extinction, and effective radius estimates for 1152 star clusters in a region of ~7.3 × 8.1 kpc centered on the nucleus and extending into the outer spiral arms. In this paper we present the data set and exploit it to determine the age distribution and relationships among the fundamental parameters (i.e. age, mass, effective radius). We show the critical dependence of the age distribution on the sample selection, and confirm that using a constant mass cut-off, above which the sample is complete for the entire age range of interest, is essential. In particular, in this sample we are complete only for masses above 5× 104~M? for the last 1 Gyr. Using this dataset we find: i) that the cluster formation rate seems to have had a large increase ~50-70 Myr ago, which is coincident with the suggested second passage of its companion, NGC 5195; ii) a large number of extremely young (<10 Myr) star clusters, which we interpret as a population of unbound clusters of which a large majority will disrupt within the next ~10 Myr; and iii) that the distribution of cluster sizes can be well approximated by a power-law with exponent, -? = -2.2 ± 0.2, which is very similar to that of Galactic globular clusters, indicating that cluster disruption is largely independent of cluster radius. In addition, we have used this dataset to search for correlations among the derived parameters. In particular, we do not find any strong trends between the age and mass, mass and effective radius, nor between the galactocentric distance and effective radius. There is, however, a strong correlation between the age of a cluster and its extinction, with younger clusters being more heavily reddened than older clusters

Hierarchical Star-Formation in M33: Fundamental properties of the star-forming regions

Star-formation within galaxies appears on multiple scales, from spiral structure, to OB associations, to individual star clusters, and often sub-structure within these clusters. This multitude of scales calls for objective methods to find and classify star-forming regions, regardless of spatial size. To this end, we present an analysis of star-forming groups in the local group spiral galaxy M33, based on a new implementation of the Minimum Spanning Tree (MST) method. Unlike previous studies which limited themselves to a single spatial scale, we study star-forming structures from the effective resolution limit (~20pc) to kpc scales. We find evidence for a continuum of star-forming group sizes, from pc to kpc scales. We do not find a characteristic scale for OB associations, unlike that found in previous studies, and we suggest that the appearance of such a scale was caused by spatial resolution and selection effects. The luminosity function of the groups is found to be well represented by a power-law with an index, -2, similar to that found for clusters and GMCs. Additionally, the groups follow a similar mass-radius relation as GMCs. The size distribution of the groups is best described by a log-normal distribution and we show that within a hierarchical distribution, if a scale is selected to find structure, the resulting size distribution will have a log-normal distribution. We find an abrupt drop of the number of groups outside a galactic radius of ~4kpc, suggesting a change in the structure of the star-forming ISM, possibly reflected in the lack of GMCs beyond this radius. (abridged)Comment: 12 pages, 16 figures, accepted MNRA

Monte Carlo Simulations of Globular Cluster Evolution. V. Binary Stellar Evolution

We study the dynamical evolution of globular clusters containing primordial binaries, including full single and binary stellar evolution using our Monte Carlo cluster evolution code updated with an adaptation of the single and binary stellar evolution codes SSE/BSE from Hurley et. al (2000, 2002). We describe the modifications we have made to the code. We present several test calculations and comparisons with existing studies to illustrate the validity of the code. We show that our code finds very good agreement with direct N-body simulations including primordial binaries and stellar evolution. We find significant differences in the evolution of the global properties of the simulated clusters using stellar evolution compared to simulations without any stellar evolution. In particular, we find that the mass loss from stellar evolution acts as a significant energy production channel simply by reducing the total gravitational binding energy and can significantly prolong the initial core contraction phase before reaching the binary-burning quasi steady state of the cluster evolution as noticed in Paper IV. We simulate a large grid of clusters varying the initial cluster mass, binary fraction, and concentration and compare properties of the simulated clusters with those of the observed Galactic globular clusters (GGCs). We find that our simulated cluster properties agree well with the observed GGC properties. We explore in some detail qualitatively different clusters in different phases of their evolution, and construct synthetic Hertzprung-Russell diagrams for these clusters.Comment: 46 preprint pages, 18 figures, 3 tables, submitted to Ap

Understanding the Dynamical State of Globular Clusters: Core-Collapsed vs Non Core-Collapsed

We study the dynamical evolution of globular clusters using our H\'enon-type Monte Carlo code for stellar dynamics including all relevant physics such as two-body relaxation, single and binary stellar evolution, Galactic tidal stripping, and strong interactions such as physical collisions and binary mediated scattering. We compute a large database of several hundred models starting from broad ranges of initial conditions guided by observations of young and massive star clusters. We show that these initial conditions very naturally lead to present day clusters with properties including the central density, core radius, half-light radius, half-mass relaxation time, and cluster mass, that match well with those of the old Galactic globular clusters. In particular, we can naturally reproduce the bimodal distribution in observed core radii separating the "core-collapsed" vs the "non core-collapsed" clusters. We see that the core-collapsed clusters are those that have reached or are about to reach the equilibrium "binary burning" phase. The non core-collapsed clusters are still undergoing gravo-thermal contraction.Comment: 42 pages, 12 figures, 1 table, submitted to MNRA

Lifetimes of tidally limited star clusters with different radii

We study the escape rate, dN/dt, from clusters with different radii in a tidal field using analytical predictions and direct N-body simulations. We find that dN/dt depends on the ratio R=r_h/r_j, where r_h is the half-mass radius and r_j the radius of the zero-velocity surface. For R>0.05, the "tidal regime", there is almost no dependence of dN/dt on R. To first order this is because the fraction of escapers per relaxation time, t_rh, scales approximately as R^1.5, which cancels out the r_h^1.5 term in t_rh. For R<0.05, the "isolated regime", dN/dt scales as R^-1.5. Clusters that start with their initial R, Ri, in the tidal regime dissolve completely in this regime and their t_dis is insensitive to the initial r_h. We predicts that clusters that start with Ri<0.05 always expand to the tidal regime before final dissolution. Their t_dis has a shallower dependence on Ri than what would be expected when t_dis is a constant times t_rh. For realistic values of Ri, the lifetime varies by less than a factor of 1.5 due to changes in Ri. This implies that the "survival" diagram for globular clusters should allow for more small clusters to survive. We note that with our result it is impossible to explain the universal peaked mass function of globular cluster systems by dynamical evolution from a power-law initial mass function, since the peak will be at lower masses in the outer parts of galaxies. Our results finally show that in the tidal regime t_dis scales as N^0.65/w, with w the angular frequency of the cluster in the host galaxy. [ABRIDGED

The spatial distribution of star and cluster formation in M51

Aims. We study the connection between spatially resolved star formation and young star clusters across the disc of M51. Methods. We combine star cluster data based on B, V, and I-band Hubble Space Telescope ACS imaging, together with new WFPC2 U-band photometry to derive ages, masses, and extinctions of 1580 resolved star clusters using SSP models. This data is combined with data on the spatially resolved star formation rates and gas surface densities, as well as Halpha and 20cm radio-continuum (RC) emission, which allows us to study the spatial correlations between star formation and star clusters. Two-point autocorrelation functions are used to study the clustering of star clusters as a function of spatial scale and age. Results. We find that the clustering of star clusters among themselves decreases both with spatial scale and age, consistent with hierarchical star formation. The slope of the autocorrelation functions are consistent with projected fractal dimensions in the range of 1.2-1.6, which is similar to other galaxies, therefore suggesting that the fractal dimension of hierarchical star formation is universal. Both star and cluster formation peak at a galactocentric radius of 2.5 and 5 kpc, which we tentatively attribute to the presence of the 4:1 resonance and the co-rotation radius. The positions of the youngest (<10 Myr) star clusters show the strongest correlation with the spiral arms, Halpha, and the RC emission, and these correlations decrease with age. The azimuthal distribution of clusters in terms of kinematic age away from the spiral arms indicates that the majority of the clusters formed 5-20 Myr before their parental gas cloud reached the centre of the spiral arm.Comment: 14 pages, 21 figures, accepted for publication in A&

Wide-field HST/ACS images of M81: The Population of Compact Star Clusters

We study the population of compact stellar clusters (CSCs) in M81, using the HST/ACS images in the filters F435W, F606W and F814W covering, for the first time, the entire optical extent of the galaxy. Our sample contains 435 clusters of FWHM less than 10 ACS pixels (9 pc). The sample shows the presence of two cluster populations, a blue group of 263 objects brighter than B=22 mag, and a red group of 172 objects, brighter than B=24 mag. Based on the analysis of colour magnitude diagrams and making use of simple stellar population models, we find the blue clusters are younger than 300 Myr with some clusters as young as few Myr, and the red clusters are as old as globular clusters. The luminosity function of the blue group follows a power-law distribution with an index of 2.0, typical value for young CSCs in other galaxies. The power-law shows unmistakable signs of truncation at I=18.0 mag (M_I=-9.8 mag), which would correspond to a mass-limit of 4x10^4 M_solar if the brightest clusters are younger than 10 Myr. The red clusters have photometric masses between 10^5 to 2x10^7 M_solar for the adopted age of 5 Gyr and their luminosity function resembles very much the globular cluster luminosity function in the Milky Way. The brightest GC in M81 has M_B^0=-10.3 mag, which is ~0.9 mag brighter than w-Cen, the most massive GC in the Milky Way.Comment: Accepted by MNRAS. The paper contains 10 figures and 3 tables. Table 3 will be published in full online onl