The globular clusters–stellar haloes connection in early-type galaxies

This paper explores if, and to what an extent, the stellar populations of early-type galaxies can be traced through the colour distribution of their globular cluster (GC) systems. The analysis, based on a galaxy sample from the Virgo Advanced Camera for Surveys data, is an extension of a previous approach that has been successful in the cases of the giant ellipticals NGC 1399 and NGC 4486, and assumes that the two dominant GC populations form along diffuse stellar populations sharing the cluster chemical abundances and spatial distributions. The results show that (a) integrated galaxy colours can be matched to within the photometric uncertainties and are consistent with a narrow range of ages; (b) the inferred mass to luminosity ratios and stellar masses are within the range of values available in the literature; (c) most GC systems occupy a thick plane in the volume space defined by the cluster formation efficiency, total stellar mass and projected surface mass density. The formation efficiency parameter of the red clusters shows a dependency with projected stellar mass density that is absent for the blue globulars. In turn, the brightest galaxies appear clearly detached from that plane as a possible consequence of major past mergers; (d) the stellar mass–metallicity relation is relatively shallow but shows a slope change at   M _*≈ 10¹⁰ M_⊙  . Galaxies with smaller stellar masses show predominantly unimodal GC colour distributions. This result may indicate that less massive galaxies are not able to retain chemically enriched interstellar matter.     

Galaxy evolution in the high-redshift, colour-selected cluster RzCS 052 at z= 1.02

We present deep I and z ' imaging of the colour-selected cluster RzCS 052 and study the colour–magnitude relation of this cluster, its scatter, the morphological distribution on the red sequence, the luminosity and stellar mass functions of red galaxies and the cluster blue fraction. We find that the stellar populations of early-type galaxies in this cluster are uniformly old and that their luminosity function does not show any sign of evolution other than the passive evolution of their stellar populations. We rule out a significant contribution from mergers in the build-up of the red sequence of RzCS 052. The cluster has a large (∼30 per cent) blue fraction and we infer that the evolution of the blue galaxies is faster than an exponentially declining star formation model and that these objects have probably experienced starburst episodes. Mergers are unlikely to be the driver of the observed colour evolution, because of the measured constancy of the mass function, as derived from near-infrared photometry of 32 clusters, including RzCS 052, presented in a related paper. Mechanisms with clustercentric radial dependent efficiencies are disfavoured as well, because of the observed constant blue fraction with clustercentric distance.     

Principal component analysis for spectral indices of stellar populations

We apply the method of principal component analysis to a sample of simple stellar populations to select some age-sensitive spectral indices. Besides the well-known age-sensitive index, H β , we find that some other spectral indices have great potential to determine the age of stellar populations, such as G4300, Fe4383, C₂4668, and Mg b . In addition, we find that the sensitivity to age of these spectral indices depends on the metallicity of the simple stellar population (SSP): H β and G4300 are more suited to determine the age of the low-metallicity stellar population, C₂4668 and Mg b are more suited to the high-metallicity stellar population. The results suggest that the principal component analysis method provides a more objective and informative alternative to diagnostics by individual spectral lines.     

Evolutionary history of the elliptical galaxy NGC 1052

We have obtained Keck spectra for 16 globular clusters (GCs) associated with the merger remnant elliptical NGC 1052, as well as a long-slit spectrum of the galaxy. We derive ages, metallicities and abundance ratios from simple stellar population models using the recently published methods of Proctor & Sansom , applied to extragalactic GCs for the first time. A number of GCs indicate the presence of strong blue horizontal branches that are not fully accounted for in the current stellar population models. We find all of the GCs to be ∼13 Gyr old according to simple stellar populations, with a large range of metallicities. From the galaxy spectrum we find NGC 1052 to have a luminosity-weighted central age of ∼2 Gyr and metallicity of  [Fe/H]∼+0.6  . No strong gradients in either age or metallicity were found to the maximum radius measured  (0.3  r _e≃ 1 kpc)  . However, we do find a strong radial gradient in α-element abundance, which reaches a very high central value. The young central starburst age is consistent with the age inferred from the H  i tidal tails and infalling gas of ∼1 Gyr. Thus, although NGC 1052 shows substantial evidence for a recent merger and an associated starburst, it appears that the merger did not induce the formation of new GCs, perhaps suggesting that little recent star formation occurred. This interpretation is consistent with 'frosting' models for early-type galaxy formation.     

Inclusion of binaries in evolutionary population synthesis

Using evolutionary population synthesis we present integrated colours, integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as the Lick/IDS) system, for an extensive set of instantaneous-burst binary stellar populations with and without binary interactions. The ages of the populations are in the range 1–15 Gyr and the metallicities are in the range 0.0001–0.03. By comparing the results for populations with and without binary interactions we show that the inclusion of binary interactions makes the integrated U – B , B – V , V – R and R – I colours and all Lick/IDS spectral absorption indices (except for H_β) substantially smaller. In other words, binary evolution makes a population appear bluer. This effect raises the derived age and metallicity of the population.
We calculate several sets of additional solar-metallicity binary stellar populations to explore the influence of the binary evolution algorithm input parameters (the common-envelope ejection efficiency and the stellar wind mass-loss rate) on the resulting integrated colours. We also look at the dependence on the choice of distribution functions used to generate the initial binary population. The results show that variations in the choice of input model parameters and distributions can significantly affect the results. However, comparing the discrepancies that exist between the colours of various models, we find that the differences are less than those produced between the models with and those without binary interactions. Therefore it is very necessary to consider binary interactions in order to draw accurate conclusions from evolutionary population synthesis work.     

On the self-consistency of evolutionary synthesis models

Evolutionary synthesis models have been used to study the physical properties of unresolved populations in a wide range of scenarios. Unfortunately, their self-consistency is difficult to test and there are some theoretical open questions without an answer: (1) The change of the homology relations assumed in the computation of isochrones due to the effect of stellar winds (or rotation) and the discontinuities in the stellar evolution are not considered. (2) There is no consensus about how the isochrones must be integrated. (3) The discreteness of the stellar populations (that produce an intrinsic statistical dispersion) usually are not taken into account, and model results are interpreted in a deterministic way instead of a statistical one. The objective of this contribution is to present some inconsistencies in the computation and some cautions in the application of the results of such codes. This revised version was published online in August 2006 with corrections to the Cover Date.     

The influence of binary interactions on infrared passbands of populations

In our evolutionary population synthesis models, samples of binaries are reproduced by a ' patched ' Monte Carlo simulation and the stellar masses, integrated   J ,  H ,  K ,  L ,  L 2  and M magnitudes, mass-to-light ratios and broad colours involving infrared bands are presented for an extensive set of instantaneous-burst binary stellar populations. In addition, the fluctuations in the integrated colours, which have been given by Zhang et al., are reduced.
By comparing the results for binary stellar populations with (Model A) and without (Model B) binary interactions, we show that the inclusion of binary interactions makes the stellar mass of a binary stellar population smaller (  ∼3.6–4.5  per cent during the past 15 Gyr), magnitudes greater (except   U , ∼ 0.18 mag  at the most), colours bluer (∼0.15 mag for   V − K   at the most) and mass-to-light ratios greater (∼0.06 for K band) except those in the U and B passbands at higher metallicities. Binary interactions make the V magnitude less sensitive to age, and R and I magnitudes more sensitive to metallicity.
Given an age, the absolute values of the differences in the stellar mass, magnitudes and mass-to-light ratios (except those in the U and B bands) between Models A and B reach a maximum at   Z = 0.0001  , i.e. the effects of binary interactions on these parameters reach a maximum, while the differences in some colours reach a maximum at   Z ∼ 0.01–0.0004  . In contrast, the absolute value of the difference in the stellar mass is minimal at   Z = 0.03  ; those in the   U ,  B ,  V   magnitudes and the mass-to-light ratios in the U and B bands reach a minimum at   Z ∼ 0.01–0.004  .     

How special are brightest group and cluster galaxies?

We use the Sloan Digital Sky Survey (SDSS) to construct a sample of 625 brightest group and cluster galaxies (BCGs) together with control samples of non-BCGs matched in stellar mass, redshift and colour. We investigate how the systematic properties of BCGs depend on stellar mass and on their privileged location near the cluster centre. The groups and clusters that we study are drawn from the C4 catalogue of Miller et al. but we have developed improved algorithms for identifying the BCG and for measuring the cluster velocity dispersion. Since the SDSS photometric pipeline tends to underestimate the luminosities of large galaxies in dense environments, we have developed a correction for this effect which can be readily applied to the published catalogue data. We find that BCGs are larger and have higher velocity dispersions than non-BCGs of the same stellar mass, which implies that BCGs contain a larger fraction of dark matter. In contrast to non-BCGs, the dynamical mass-to-light ratio of BCGs does not vary as a function of galaxy luminosity. Hence BCGs lie on a different Fundamental Plane than ordinary elliptical galaxies. BCGs also follow a steeper Faber–Jackson relation than non-BCGs, as suggested by models in which BCGs assemble via dissipationless mergers along preferentially radial orbits. We find tentative evidence that this steepening is stronger in more massive clusters. BCGs have similar mean stellar ages and metallicities to non-BCGs of the same mass, but they have somewhat higher α/Fe ratios, indicating that star formation may have occurred over a shorter time-scale in the BCGs. Finally, we find that BCGs are more likely to host radio-loud active galactic nuclei than other galaxies of the same mass, but are less likely to host an optical active galactic nucleus (AGN). The differences we find are more pronounced for the less massive BCGs, i.e. they are stronger at the galaxy group level.     

The diverse nature of optical emission lines in brightest cluster galaxies: IFU observations of the central kiloparsec

We present integral field spectroscopy of the nebular line emission in a sample of nine brightest cluster galaxies (BCGs). The sample was chosen to probe both cooling flow and non-cooling flow clusters, as well as a range of cluster X-ray luminosities. The line emission morphology and velocity gradients suggest a great diversity in the properties of the line emitting gas. While some BCGs show evidence for filamentary or patchy emission (Abell 1060, Abell 1668 and MKW 3s), others have extended emission (Abell 1204, Abell 2199), while still others have centrally concentrated emission (Abell 2052). We examine diagnostic line ratios to determine the dominant ionization mechanisms in each galaxy. Most of the galaxies show regions with active galactic nucleus like spectra, however, for two BCGs, Abell 1060 and Abell 1204, the emission line diagnostics suggest regions which can be described by the emission from young stellar populations. The diversity of emission-line properties in our sample of BCGs suggests that the emission mechanism is not universal, with different ionization processes dominating different systems. Given this diversity, there is no evidence for a clear distinction of the emission-line properties between cooling flow and non-cooling flow BCGs. It is not always cooling flow BCGs which show emission (or young stellar populations), and non-cooling flow BCGs which do not.     

Ultraviolet dust attenuation in star-forming galaxies – II. Calibrating the A(UV) versus LTIR/LUV relation

We investigate the dependence of the total-IR (TIR) to ultraviolet (UV) luminosity ratio method for calculating the UV dust attenuation   A (UV)  from the age of the underlying stellar populations by using a library of spectral energy distributions for galaxies with different star formation histories. Our analysis confirms that the TIR/UV versus   A (UV)  relation varies significantly with the age of the underlying stellar population: that is, for the same TIR/UV ratio, systems with low specific star formation rate (SSFR) suffer a lower UV attenuation than starbursts. Using a sample of nearby field and cluster spiral galaxies, we show that the use of a standard (i.e. age-independent) TIR/UV versus   A (UV)  relation leads to a systematic overestimate up to 2 mag of the amount of UV dust attenuation suffered by objects with low SSFR and in particular H  i -deficient star-forming cluster galaxies. This result points out that the age-independent TIR/UV versus   A (UV)  relation cannot be used to study the UV properties of large samples of galaxies including low star-forming systems and passive spirals. Therefore, we give some simple empirical relations from which the UV attenuation can be estimated taking into account its dependence on the age of the stellar populations, providing a less biased view of UV properties of galaxies.     

Age and metallicity of compact stellar systems in Virgo and Fornax

From high signal-to-noise ratio GMOS-N and AAOmega spectra, we have measured line-strength indices in the Lick/IDS system for several luminous and spatially dispersed compact stellar systems (CSSs) located in the Virgo and Fornax galaxy clusters. We estimate their [α/Fe] ratios, ages and metallicities using a simple (single-burst) stellar population model. We confirm that the Virgo core region luminous CSSs (     ) contain ancient stellar populations with subsolar total metallicities, suggesting that they comprise the bright tail of M87's GC distribution. The two Virgo intracluster globular clusters have ages and metallicities consistent with the cluster core CSSs. Two Fornax luminous CSSs also have ancient stellar populations but are at the upper end of the Virgo CSS metallicity range, while the third (UCD3) appears to be relatively young, metal-rich and with a core + halo radial profile. Our results suggest that Fornax may contain an extra population of luminous CSSs formed more recently than the ancient GC-like systems found in both clusters.     

Spatially resolved kinematics and stellar populations of brightest cluster and group galaxies

We present an examination of the kinematics and stellar populations of a sample of three brightest group galaxies (BGGs) and three brightest cluster galaxies (BCGs) in X-ray groups and clusters. We have obtained high signal-to-noise ratio Gemini/Gemini South Multi-Object Spectrograph (GMOS) long-slit spectra of these galaxies and use Lick indices to determine ages, metallicities and α-element abundance ratios out to at least their effective radii. We find that the BGGs and BCGs have very uniform masses, central ages and central metallicities. Examining the radial dependence of their stellar populations, we find no significant velocity dispersion, age, or α-enhancement gradients. However, we find a wide range of metallicity gradients, suggesting a variety of formation mechanisms. The range of metallicity gradients observed is surprising, given the homogeneous environment these galaxies probe and their uniform central stellar populations. However, our results are inconsistent with any single model of galaxy formation and emphasize the need for more theoretical understanding of both the origins of metallicity gradients and galaxy formation itself. We postulate two possible physical causes for the different formation mechanisms.     

An optimized Hβ index for disentangling stellar population ages

We have defined a new Hβ absorption index definition,  Hβ_o  , which has been optimized as an age indicator for old and intermediate age stellar populations. Rather than using stellar spectra, we employed for this purpose a library of stellar population spectral energy distributions of different ages and metallicities at moderately high resolution.  Hβ_o  provides us with improved abilities for lifting the age–metallicity degeneracy affecting the standard Hβ Lick index definition. The new index, which has also been optimized against photon noise and velocity dispersion, is fully characterized with wavelength shift, spectrum shape, dust extinction and [α/Fe] abundance ratio effects.  Hβ_o  requires spectra of similar qualities as those commonly used for measuring the standard Hβ Lick index definition. Aiming at illustrating the use and capabilities of  Hβ_o  as an age indicator we apply it to Milky Way globular clusters and to a well selected sample of early-type galaxies covering a wide range in mass. The results shown here are particularly useful for applying this index and understand the involved uncertainties.     

Direct N-body modelling of stellar populations: blue stragglers in M67

We present a state-of-the-art N -body code which includes a detailed treatment of stellar and binary evolution as well as the cluster dynamics. This code is ideal for investigating all aspects relating to the evolution of star clusters and their stellar populations. It is applicable to open and globular clusters of any age. We use the N -body code to model the blue straggler population of the old open cluster M67. Preliminary calculations with our binary population synthesis code show that binary evolution alone cannot explain the observed numbers or properties of the blue stragglers. On the other hand, our N -body model of M67 generates the required number of blue stragglers and provides formation paths for all the various types found in M67. This demonstrates the effectiveness of the cluster environment in modifying the nature of the stars it contains, and highlights the importance of combining dynamics with stellar evolution. We also perform a series of N =10 000 simulations in order to quantify the rate of escape of stars from a cluster subject to the Galactic tidal field.     

Recurrent gas accretion by massive star clusters, multiple stellar populations and mass thresholds for spheroidal stellar systems

We explore the gravitational influence of pressure-supported stellar systems on the internal density distribution of a gaseous environment. We conclude that compact massive star clusters with masses  ≳10⁶ M_⊙  act as cloud condensation nuclei and are able to accrete gas recurrently from a warm interstellar medium which may cause further star formation events and account for multiple stellar populations in the most massive globular and nuclear star clusters. The same analytical arguments can be used to decide whether an arbitrary spherical stellar system is able to keep warm or hot interstellar material or not. These mass thresholds coincide with transition masses between pressure supported galaxies of different morphological types.     

The environmental dependence of the chemical properties of star-forming galaxies

We use a  0.040 < z < 0.085  sample of 37 866 star-forming galaxies from the Fourth Data Release of the Sloan Digital Sky Survey to investigate the dependence of gas-phase chemical properties on stellar mass and environment. The local density, determined from the projected distances to the fourth and fifth nearest neighbours, is used as an environment indicator. Considering environments ranging from voids, i.e.  log Σ≲−0.8  , to the periphery of galaxy clusters, i.e.  log Σ≈ 0.8  , we find no dependence of the relationship between galaxy stellar mass and gas-phase oxygen abundance, along with its associated scatter, on local galaxy density. However, the star-forming gas in galaxies shows a marginal increase in the chemical enrichment level at a fixed stellar mass in denser environments. Compared with galaxies of similar stellar mass in low-density environments, they are enhanced by a few per cent for massive galaxies to about 20 per cent for galaxies with stellar masses  ≲10^9.5 M_⊙  . These results imply that the evolution of star-forming galaxies is driven primarily by their intrinsic properties and is largely independent of their environment over a large range of local galaxy density.     

Intragroup diffuse light in compact groups of galaxies – II. HCG 15, 35 and 51

This continuing study of intragroup light in compact groups of galaxies aims to establish new constraints to models of formation and evolution of galaxy groups, specially of compact groups, which are a key part in the evolution of larger structures, such as clusters. In this paper we present three additional groups (HCG 15, 35 and 51) using deep wide-field B - and R -band images observed with the LAICA camera at the 3.5-m telescope at the Calar Alto observatory (CAHA). This instrument provides us with very stable flat-fielding, a mandatory condition for reliably measuring intragroup diffuse light. The images were analysed with the ov_wav package, a wavelet technique that allows us to uncover the intragroup component in an unprecedented way. We have detected that 19, 15 and 26 per cent of the total light of HCG 15, 35 and 51, respectively, are in the diffuse component, with colours that are compatible with old stellar populations and with mean surface brightness that can be as low as 28.4 B mag arcsec⁻². Dynamical masses, crossing times and mass-to-light ratios were recalculated using the new group parameters. Also tidal features were analysed using the wavelet technique.     

How well do starlab and nbody4 compare? I. Simple models

N -body simulations are widely used to simulate the dynamical evolution of a variety of systems, among them star clusters. Much of our understanding of their evolution rests on the results of such direct N -body simulations. They provide insight in the structural evolution of star clusters, as well as into the occurrence of stellar exotica. Although the major pure N -body codes starlab/kira and nbody4 are widely used for a range of applications, there is no thorough comparison study yet.
Here, we thoroughly compare basic quantities as derived from simulations performed either with starlab/kira or nbody4 .
We construct a large number of star cluster models for various stellar mass function settings (but without stellar/binary evolution, primordial binaries, external tidal fields, etc.), evolve them in parallel with starlab/kira and nbody4 , analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison, we develop a bootstrap algorithm for functional dependencies.
We find an overall excellent agreement between the codes, both for the clusters' structural and energy parameters as well as for the properties of the dynamically created binaries. However, we identify small differences, like in the energy conservation before core collapse and the energies of escaping stars, which deserve further studies.
Our results reassure the comparability and the possibility to combine results from these two major N -body codes, at least for the purely dynamical models (i.e. without stellar/binary evolution) we performed. Further detailed comparison studies for more complex systems, e.g. including stellar/binary evolution, are required.     

New age estimates of M31 globular clusters from multicolour photometry

The large majority of extragalactic star cluster studies performed to date essentially use multicolour photometry, combined with theoretical stellar synthesis models, to derive ages, masses, extinction estimates and metallicities. M31 offers a unique laboratory for studies of globular cluster (GC) systems. In this paper, we obtain new age estimates for 91 M31 GCs, based on improved photometric data, updated theoretical stellar synthesis models and sophisticated new fitting methods. In particular, we used photometric measurements from the Two Micron All Sky Survey (2MASS), which, in combination with optical photometry, can partially break the well-known age–metallicity degeneracy operating at ages in excess of a few Gyr. We show robustly that previous age determinations based on photometric data were affected significantly by this age–metallicity degeneracy. Except for one cluster, the ages of our other sample GCs are all older than 1 Gyr. Their age distribution shows populations of young- and intermediate-age GCs, peaking at ∼3 and 8 Gyr, respectively, as well as the 'usual' complement of well-known old GCs, i.e. GCs of similar age as the majority of the Galactic GCs. Our results also show that although there is significant scatter in metallicity at any age, there is a notable lack of young metal-poor and old metal-rich GCs, which might be indicative of an underlying age–metallicity relationship among the M31 GC population.     

The mass–metallicity relation in galaxy clusters: the relative importance of cluster membership versus local environment

Using a large (14 857), homogenously selected sample of cluster galaxies identified in the Sloan Digital Sky Survey Data Release 4, we investigate the impact of cluster membership and local density on the stellar mass–gas phase metallicity relation (MZR). We show that stellar metallicities are not suitable for this work, being relatively insensitive to subtle changes in the MZR. Accurate nebular abundances can be obtained for 1318 cluster galaxies in our sample and we show that these galaxies are drawn from clusters that are fully representative of the parent sample in terms of mass, size, velocity dispersion and richness. By comparing the MZR of the cluster galaxies with a sample of control galaxies matched in mass, redshift, fibre covering fraction and rest-frame   g − r   colour cluster galaxies are found to have, on average, higher metallicities by up to 0.04 dex. The magnitude of this offset does not depend strongly on galactic half-light radius or cluster properties such as velocity dispersion or cluster mass. The effect of local density on the MZR is investigated, using the presence of a near neighbour and both two- and three-dimensional density estimators. For all three metrics, it is found that the cluster galaxies in locally rich environments have higher median metallicities by up to ∼0.05 dex than those in locally poor environments (or without a near neighbour). Control (non-cluster) galaxies at locally high densities exhibit similar metal enhancements. Taken together, these results show that galaxies in clusters are, on average, slightly more metal rich than the field, but that this effect is driven by local overdensity and not simply cluster membership.