The importance of interloper removal in galaxy clusters: saving more objects for the Jeans analysis

We study the effect of contamination by interlopers in kinematic samples of galaxy clusters. We demonstrate that without the proper removal of interlopers the inferred parameters of the mass distribution in the cluster are strongly biased towards higher mass and lower concentration. The interlopers are removed using two procedures previously shown to work most efficiently on simulated data. One is based on using the virial mass estimator and calculating the maximum velocity available to cluster members and the other relies on the ratio of the virial and projected mass estimators. We illustrate the performance of the methods in detail using the example of A576, a cluster with a strong uniform background contamination, and compare the case of A576 to 15 other clusters with different degree of contamination. We model the velocity dispersion and kurtosis profiles obtained for the cleaned data samples of these clusters solving the Jeans equations to estimate the mass, concentration and anisotropy parameter. We present the mass–concentration relation for the total sample of 22 clusters.     

Populating a cluster of galaxies – I. Results at

We simulate the assembly of a massive rich cluster and the formation of its constituent galaxies in a flat, low-density universe. Our most accurate model follows the collapse, the star formation history and the orbital motion of all galaxies more luminous than the Fornax dwarf spheroidal, while dark halo structure is tracked consistently throughout the cluster for all galaxies more luminous than the SMC. Within its virial radius this model contains about     dark matter particles and almost 5000 distinct dynamically resolved galaxies. Simulations of this same cluster at a variety of resolutions allow us to check explicitly for numerical convergence both of the dark matter structures produced by our new parallel N -body and substructure identification codes, and of the galaxy populations produced by the phenomenological models we use to follow cooling, star formation, feedback and stellar aging. This baryonic modelling is tuned so that our simulations reproduce the observed properties of isolated spirals outside clusters. Without further parameter adjustment our simulations then produce a luminosity function, a mass-to-light ratio, luminosity, number and velocity dispersion profiles, and a morphology–radius relation which are similar to those observed in real clusters. In particular, since our simulations follow galaxy merging explicitly, we can demonstrate that it accounts quantitatively for the observed cluster population of bulges and elliptical galaxies.     

The evolution of substructure – III. The outskirts of clusters

We present an investigation of satellite galaxies in the outskirts of galaxy clusters taken from a series of high-resolution N -body simulations. We focus on the so-called backsplash population, i.e. satellite galaxies that once were inside the virial radius of the host but now reside beyond it. We find that this population is significant in number and needs to be appreciated when interpreting the various galaxy morphology environmental relationships and decoupling the degeneracy between nature and nurture. Specifically, we find that approximately half of the galaxies with current cluster-centric distance in the interval 1–2 virial radii of the host are backsplash galaxies that once penetrated deep into the cluster potential, with 90 per cent of these entering to within 50 per cent of the virial radius. These galaxies have undergone significant tidal disruption, losing on average 40 per cent of their mass. This results in a mass function for the backsplash population different from those galaxies infalling for the first time. We further show that these two populations are kinematically distinct and should be observable within existent spectroscopic surveys.     

The concentration–velocity dispersion relation in galaxy groups

Based on results from cold dark matter N -body simulations, we develop a dynamical model for the evolution of subhaloes within group-sized host haloes. Only subhaloes more massive than 5 × 10⁸ M_⊙ are considered, because they are massive enough to possibly host luminous galaxies. On their orbits within a growing host potential the subhaloes are subject to tidal stripping and dynamical friction. At the present time  ( z = 0)  , all model hosts have equal mass  ( M _vir= 3.9 × 10¹³ M_⊙)  but different concentrations associated with different formation times. We investigate the variation of subhalo (or satellite galaxy) velocity dispersion with host concentration and/or formation time. In agreement with the Jeans equation, the velocity dispersion of subhaloes increases with the host concentration. Between concentrations of ∼5 and ∼20, the subhalo velocity dispersions increase by a factor of ∼1.25. By applying a simplified tidal disruption criterion, that is, rejection of all subhaloes with a tidal truncation radius below 3  kpc at   z = 0  , the central velocity dispersion of the 'surviving' subhalo sample increases substantially for all concentrations. The enhanced central velocity dispersions in the surviving subhalo samples are caused by a lack of slow tangential motions. Additionally, we present a fitting formula for the anisotropy parameter which does not depend on concentration if the group-centric distances are scaled by r _s, the characteristic radius of the Navarro, Frenk & White profile. Since the expected loss of subhaloes and galaxies due to tidal disruption increases the velocity dispersion of surviving galaxies, the observed galaxy velocity dispersion can substantially overestimate the virial mass.     

Enhanced mergers of galaxies in low-redshift clusters

An ensemble cluster has been formed from a data set comprising a complete magnitude-limited sample of 680 giant galaxies  ( M ⁰ _B ≲−19)  in eight low-redshift clusters, normalized by the velocity dispersions and virial radii for the early-type cluster populations. Distinct galaxy populations have been identified, including an infall population. A majority (50–70 per cent or greater) of the infall population are found to be in interacting or merging systems characterized by slow gravitational encounters. The observed enhancement of galaxy–galaxy encounters in the infall population compared to the field can be explained by gravitational shocking. It is shown that disc galaxy mergers in the infall population integrated over the estimated lifetime of the cluster (∼10 Gyr) can readily account for the present cluster S0 population.     

The universal rotation curve of spiral galaxies – II. The dark matter distribution out to the virial radius

In the current ΛCDM cosmological scenario, N -body simulations provide us with a universal mass profile, and consequently a universal equilibrium circular velocity of the virialized objects, as galaxies. In this paper we obtain, by combining kinematical data of their inner regions with global observational properties, the universal rotation curve of disc galaxies and the corresponding mass distribution out to their virial radius. This curve extends the results of Paper I, concerning the inner luminous regions of Sb–Im spirals, out to the edge of the galaxy haloes.     

Mass estimation in the outer regions of galaxy clusters

We present a technique for estimating the mass in the outskirts of galaxy clusters where the usual assumption of dynamical equilibrium is not valid. The method assumes that clusters form through hierarchical clustering and requires only galaxy redshifts and positions on the sky. We apply the method to dissipationless cosmological N -body simulations where galaxies form and evolve according to semi-analytic modelling. The method recovers the actual cluster mass profile within a factor of 2 to several megaparsecs from the cluster centre. This error originates from projection effects, sparse sampling, and contamination by foreground and background galaxies. In the absence of velocity biases, this method can provide an estimate of the mass-to-light ratio on scales ∼1–10  h ⁻¹ Mpc where this quantity is still poorly known.     

Effects of cluster galaxies on arc statistics

We present the results of a set of numerical simulations evaluating the effect of cluster galaxies on arc statistics.
We perform a first set of gravitational lensing simulations using three independent projections for each of nine different galaxy clusters obtained from N -body simulations. The simulated clusters consist of dark matter only. We add a population of galaxies to each cluster, mimicking the observed luminosity function and the spatial galaxy distribution, and repeat the lensing simulations including the effects of cluster galaxies, which themselves act as individual lenses. Each galaxy is represented by a spherical Navarro, Frenk & White density profile.
We consider the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of lengths, widths, curvature radii and length-to-width ratios of long arcs. We find some changes to the distribution of short-arc properties in the presence of cluster galaxies. The differences appear in the distribution of curvature radii for arc lengths smaller than 12 arcsec, while the distributions of lengths, widths and length-to-width ratios are significantly changed only for arcs shorter than 4 arcsec.     

Properties of spherical galaxies and clusters with an NFW density profile

Using the standard dynamical theory of spherical systems, we calculate the properties of spherical galaxies and clusters whose density profiles obey the universal form first obtained in high-resolution cosmological N -body simulations by Navarro, Frenk & White (NFW). We adopt three models for the internal kinematics: isotropic velocities, constant anisotropy and increasingly radial OsipkovMerritt anisotropy. Analytical solutions are found for the radial dependence of the mass, gravitational potential, velocity dispersion, energy and virial ratio and we test their variability with the concentration parameter describing the density profile and amount of velocity anisotropy. We also compute structural parameters, such as half-mass radius, effective radius and various measures of concentration. Finally, we derive projected quantities, the surface mass density and line-of-sight as well as aperture-velocity dispersion, all of which can be directly applied in observational tests of current scenarios of structure formation. On the mass scales of galaxies, if constant mass-to-light is assumed, the NFW surface density profile is found to fit HubbleReynolds laws well. It is also well fitted by Sérsic R ^{1/ m} laws, for     but in a much narrower range of m and with much larger effective radii than are observed. Assuming in turn reasonable values of the effective radius, the mass density profiles imply a mass-to-light ratio that increases outwards at all radii.     

Probing galactic dark matter in dense environments: on the strong lensing efficiency of galaxies in rich clusters

The recent detection by Limousin et al. of five new strong lensing events dominated by galaxy cluster members in Abell 1689, and outside the critical regime of the cluster itself, offers a way to obtain constraints on the cluster mass distribution in a region inaccessible to standard lensing analysis. In addition, modelling such systems will provide another window on the dark matter haloes of galaxies in very dense environments. Here, it is shown that the boost in image separation due to the external shear and convergence from a smooth cluster component means that more numerous, less massive galaxies have the potential to create multiple images with detectable separations, relative to isolated field galaxies. This comes in addition to a potential increase in their lensing (source plane) cross-section. To gain insight into the factors involved and as a precursor to a numerical study using N -body simulations, a simple analytic model of a cluster at   z = 0.3  lensing background galaxies at   z = 2  is considered here. The fiducial model has cluster members with isothermal density profiles and luminosities L , distributed in a Schechter function (faint-end slope  ν=−1.25  ), related to their velocity dispersions σ via the Faber–Jackson scaling L ∝σ⁴. Just outside the critical regime of the cluster, the scale of galaxy-dominated image separations is significantly increased. Folding in the fact that less massive galaxies present a lower lensing cross-section, and that the cross-section can itself be enhanced in an external field leads to a factor of a few times more detected events relative to field galaxies. These values will be higher closer to the critical curve. Given that the events in Abell 1689 were detected over a very small region of the cluster where ACS data were available, this motivates the search for such events in other clusters.     

-ray wakes as probes of galaxy cluster dynamics

If a galaxy resides in a cluster, then its passage through the pervasive intracluster medium will produce a detectable signature in the X-ray emission from the cluster. Such features have now been detected in a number of systems. The simplest kinematic information that can be extracted from this signature is the galaxy's direction of motion on the plane of the sky. This paper explores the constraints on cluster dynamics that could be derived from such information. In particular, we show that it is possible to define a projected anisotropy parameter, B ( r ), which is directly analogous to the usual orbital anisotropy parameter. We describe an estimator for this quantity, ( R ), which can be derived in a robust and straightforward manner. We present a simple dynamical model for a cluster consisting of a Michie distribution function of galaxies orbiting in a truncated singular isothermal sphere potential. Using this model, we demonstrate the ambiguity between the distribution of mass and the distribution of galaxy orbits when interpreting the traditional measures of cluster kinematics (the projected density of galaxies and their line-of-sight velocity dispersion). As an example, we show how two very different dynamical models can fit the kinematic properties of the Coma cluster. We demonstrate that the measurement of using a relatively small sample of wake directions ( N _wake≈50) would provide an effective mechanism for lifting this degeneracy. Thus, by combining X-ray measurements of wake directions with number counts and line-of-sight velocities derived from optical data, it will prove possible to measure both the orbit distribution and the form of the gravitational potential in clusters of galaxies. The requisite X-ray observations lie within reach of the soon-to-be-launched AXAF satellite.     

Spin vector orientation of galaxies in the region 15h 48m≤α(2000) ≤ 19h 28m, −68°≤δ(2000) ≤−62°

We present an analysis of the orientations of 1433 galaxies found in the region  15^h 48^m≤α(2000) ≤ 19^h 28^m, −68°≤δ(2000) ≤−62°  . In this region we investigated three Abell clusters (S0794, S0797, S0805) of richness Class 0 and the Triangulum Australis cluster. Our aim is to examine non-random effects in galaxy orientations in clusters. In addition, we classified the investigated galaxies into subsamples on the basis of their axial ratio, major diameter and morphology. The spin vector orientations of total galaxies in the investigated region is found to be random. No preferred orientation is found in the clusters. We could not note any morphological dependence of the galaxy orientations in our samples. No preferred orientations can be seen for the spiral galaxies. The morphologically unidentified galaxies, galaxies having major diameters of <47 arcsec, and the nearly edge-on galaxies  ( b / a < 0.5, 0.4 < b / a ≤ 0.5)  show anisotropy: spin vectors of galaxies tend to be oriented perpendicular to the Local Supercluster plane and spin vector projections tend to point radially with respect to the Virgo cluster centre.     

Mass profiles and anisotropies of early-type galaxies

We discuss the problem of using stellar kinematics of early-type galaxies to constrain the orbital anisotropies and radial mass profiles of galaxies. We demonstrate that compressing the light distribution of a galaxy along the line of sight produces approximately the same signature in the line-of-sight velocity profiles as radial anisotropy. In particular, fitting spherically symmetric dynamical models to apparently round, isotropic face-on flattened galaxies leads to a spurious bias towards radial orbits in the models, especially if the galaxy has a weak face-on stellar disc. Such face-on stellar discs could plausibly be the cause of the radial anisotropy found in spherical models of intermediate luminosity ellipticals such as NGC 2434, 3379 and 6703.
In the light of this result, we use simple dynamical models to constrain the outer mass profiles of a sample of 18 round, early-type galaxies. The galaxies follow a Tully–Fisher relation parallel to that for spiral galaxies, but fainter by at least 0.8 mag ( I -band) for a given mass. The most luminous galaxies show clear evidence for the presence of a massive dark halo, but the case for dark haloes in fainter galaxies is more ambiguous. We discuss the observations that would be required to resolve this ambiguity.     

Luminosity distributions within rich clusters — III. A comparative study of seven Abell/ACO clusters

We recover the luminosity distributions over a wide range of absolute magnitude (−24.5 <  M _R  < −16.5) for a sample of seven rich southern galaxy clusters. We find a large variation in the ratio of dwarf to giant galaxies, DGR: 0.8 ≤ DGR ≤ 3.1. This variation is shown to be inconsistent with a ubiquitous cluster luminosity function. The DGR shows a smaller variation from cluster to cluster in the inner regions ( r  ≲ 0.56 Mpc). Outside these regions we find the DGR to be strongly anticorrelated with the mean local projected galaxy density, with the DGR increasing towards lower densities. In addition, the DGR in the outer regions shows some correlation with Bautz–Morgan type. Radial analysis of the clusters indicates that the dwarf galaxies are less centrally clustered than the giants, and they form a significant halo around clusters. We conclude that measurements of the total cluster luminosity distribution based on the inner core alone are likely to be severe underestimates of the dwarf component, the integrated cluster luminosity and the contribution of galaxy masses to the cluster's total mass. Further work is required to quantify this. The observational evidence that the unrelaxed, lower density outer regions of clusters are dwarf-rich adds credence to the recent evidence and conjecture that the field is a predominantly dwarf-rich environment, and that the dwarf galaxies are under-represented in measures of the local field luminosity function.     

Star formation in galaxies falling into clusters along supercluster-scale filaments

With the help of a statistical parameter derived from optical spectra, we show that the current star formation rate of a galaxy, falling into a cluster along a supercluster filament, is likely to undergo a sudden enhancement before the galaxy reaches the virial radius of the cluster. From a sample of 52 supercluster-scale filaments of galaxies joining a pair of rich clusters of galaxies within the two-degree Field Redshift Survey region, we find a significant enhancement of star formation, within a narrow range between ∼2 and  3  h ⁻¹₇₀ Mpc  of the centre of the cluster into which the galaxy is falling. This burst of star formation is almost exclusively seen in the fainter dwarf galaxies  ( M _B ≥−20)  . The relative position of the peak does not depend on whether the galaxy is a member of a group or not, but non-group galaxies have on average a higher rate of star formation immediately before falling into a cluster. From the various trends, we conclude that the predominant process responsible for this rapid burst is the close interaction with other galaxies falling into the cluster along the same filament, if the interaction occurs before the gas reservoir of the galaxy gets stripped off due to the interaction with the intracluster medium.     

Dark matter haloes within clusters

We examine the properties of dark matter haloes within a rich galaxy cluster using a high-resolution simulation that captures the cosmological context of a cold dark matter universe. The mass and force resolution permit the resolution of 150 haloes with circular velocities larger than 80 km s⁻¹ within the cluster virial radius of 2 Mpc (with Hubble constant H ₀ = 50 km s⁻¹ Mpc⁻¹). This enables an unprecedented study of the statistical properties of a large sample of dark matter haloes evolving in a dense environment. The cumulative fraction of mass attached to these haloes varies from close to zero per cent at 200 kpc to 13 per cent at the virial radius. Even at this resolution the overmerging problem persists; haloes that pass within 100–200 kpc of the cluster centre are tidally disrupted. Additional substructure is lost at earlier epochs within the massive progenitor haloes. The median ratio of apocentric to pericentric radii is 6:1, so that the orbital distribution is close to isotropic, circular orbits are rare and radial orbits are common. The orbits of haloes are unbiased with respect to both position within the cluster and the orbits of the smooth dark matter background, and no velocity bias is detected. The tidal radii of surviving haloes are generally well-fitted using the simple analytic prediction applied to their orbital pericentres. Haloes within clusters have higher concentrations than those in the field. Within the cluster, halo density profiles can be modified by tidal forces and individual encounters with other haloes that cause significant mass loss —'galaxy harassment'. Mergers between haloes do not occur inside the cluster virial radius.     

Ram-pressure histories of cluster galaxies

Ram-pressure stripping can remove significant amounts of gas from galaxies that orbit in clusters and massive groups, and thus has a large impact on the evolution of cluster galaxies. In this paper, we reconstruct the present-day distribution of ram pressure and the ram-pressure histories of cluster galaxies. To this aim, we combine the Millennium Simulation and an associated semi-analytic model of galaxy evolution with analytic models for the gas distribution in clusters. We find that about one quarter of galaxies in massive clusters are subject to strong ram pressures that are likely to cause an expedient loss of all gas. Strong ram pressures occur predominantly in the inner core of the cluster, where both the gas density and the galaxy velocity are higher. Since their accretion on to a massive system, more than 64 per cent of galaxies that reside in a cluster today have experienced strong ram pressures of  >10⁻¹¹ dyn cm⁻²  which most likely led to a substantial loss of the gas.     

Satellite systems around galaxies in hydrodynamic simulations

We investigate the properties of satellite galaxies formed in N -body/SPH simulations of galaxy formation in the ΛCDM cosmology. The simulations include the main physical effects thought to be important in galaxy formation and, in several cases, produce realistic spiral discs. In total, a sample of nine galaxies of luminosity comparable to the Milky Way was obtained. At magnitudes brighter than the resolution limit,   M_V =−12  , the luminosity function of the satellite galaxies in the simulations is in excellent agreement with data for the Local Group. The radial number density profile of the model satellites, as well as their gas fractions also match observations very well. In agreement with previous N -body studies, we find that the satellites tend to be distributed in highly flattened configurations whose major axis is aligned with the major axis of the (generally triaxial) dark halo. In two out of three systems with sufficiently large satellite populations, the satellite system is nearly perpendicular to the plane of the galactic disc, a configuration analogous to that observed in the Milk Way. The discs themselves are perpendicular to the minor axis of their host haloes in the inner parts, and the correlation between the orientation of the galaxy and the shape of the halo persists even out to the virial radius. However, in one case the disc's minor axis ends up, at the virial radius, perpendicular to the minor axis of the halo. The angular momenta of the galaxies and their host halo tend to be well aligned.     

Ram pressure stripping of spiral galaxies in clusters

We use three-dimensional SPH/ N -body simulations to study ram pressure stripping of gas from spiral galaxies orbiting in clusters. We find that the analytic expectation of Gunn & Gott, relating the gravitational restoring force provided by the disc to the ram pressure force, provides a good approximation to the radius at which gas will be stripped from a galaxy. However, at small radii it is also important to consider the potential provided by the bulge component. A spiral galaxy passing through the core of a rich cluster, such as Coma, will have its gaseous disc truncated to ∼4 kpc, thus losing ∼80 per cent of its diffuse gas mass. The time-scale for this to occur is a fraction of a crossing time ∼10⁷ yr. Galaxies orbiting within poorer clusters, or inclined to the direction of motion through the intracluster medium, will lose significantly less gas. We conclude that ram pressure alone is insufficient to account for the rapid and widespread truncation of star formation observed in cluster galaxies, or the morphological transformation of Sabs to S0s that is necessary to explain the Butcher–Oemler effect.     

NoSOCS in SDSS – I. Sample definition and comparison of mass estimates

We use Sloan Digital Sky Survey (SDSS) data to investigate galaxy cluster properties of the systems first detected within Digitized Second Palomar Observatory Sky Survey. With the high-quality photometry of SDSS, we derived new photometric redshifts and estimated richness and optical luminosity. For a subset of low-redshift  ( z ≤ 0.1)  clusters, we have used SDSS spectroscopic data to identify groups in redshift space in the region of each cluster, complemented with massive systems from the literature to assure the continuous mass sampling. A method to remove interlopers is applied, and a virial analysis is performed resulting in the estimates of velocity dispersion, mass and a physical radius for each low- z system. We discuss the choice of maximum radius and luminosity range in the dynamical analysis, showing that a spectroscopic survey must be complete to at least   M *+ 1  if one wishes to obtain accurate and unbiased estimates of velocity dispersion and mass. We have measured X-ray luminosity for all clusters using archival data from ROSAT All Sky Survey. For a smaller subset (21 clusters), we selected temperature measures from the literature and estimated mass from the   M − T _X  relation, finding that they show good agreement with the virial estimate. However, these two mass estimates tend to disagree with the caustic results. We measured the presence of substructure in all clusters of the sample and found that clusters with substructure have virial masses higher than those derived from T _X. This trend is not seen when comparing the caustic and X-ray masses. That happens because the caustic mass is estimated directly from the mass profile, so it is less affected by substructure.