The colour–magnitude relation for galaxies in the Coma cluster

We present a new photometric catalogue of the Coma galaxy cluster in the Johnson U and V bands. We cover an area of 3360 arcmin² of sky, to a depth of     in a 13-arcsec diameter aperture, and produce magnitudes for ∼1400 extended objects in metric apertures from 8.8- to 26-arcsec diameters. The mean internal rms scatter in the photometry is 0.014 mag in V , and 0.026 mag in U , for     .
We place new limits on the levels of scatter in the colour–magnitude relation (CMR) in the Coma cluster, and investigate how the slope and scatter of the CMR depend on galaxy morphology, luminosity and position within the cluster. As expected, the lowest levels of scatter are found in the elliptical galaxies, while the late-type galaxies have the highest numbers of galaxies bluewards of the CMR. We investigate whether the slope of the CMR is an artefact of colour gradients within galaxies, and show that it persists when the colours are measured within a diameter that scales with galaxy size. Looking at the environmental dependence of the CMR, we find a trend of systematically bluer galaxy colours with increasing projected radius from the centre of the cluster. Surprisingly, this is accompanied by a decreased scatter of the CMR. We investigate whether this gradient could be caused by dust in the cluster potential, however the reddening required would produce too large a scatter in the colours of the central galaxies. The gradient appears to be better reproduced by a gradient in the mean galactic ages with projected radius.     

Dissipationless mergers of elliptical galaxies and the evolution of the fundamental plane

We carry out numerical simulations of dissipationless major mergers of elliptical galaxies using initial galaxy models that consist of a dark matter haloes and a stellar bulge with properties consistent with the observed fundamental plane. By varying the density profile of the dark matter haloes [standard Navarro, Frenk & White (NFW) profile versus adiabatically contracted NFW profile], the global stellar to dark matter mass ratio and the orbit of the merging galaxies, we are able to assess the impact of each of these factors on the structure of the merger remnant. Our results indicate that the properties of the remnant bulge depend primarily on the angular momentum and energy of the orbit; for a cosmologically motivated orbit, the effective radius and velocity dispersion of the remnant bulge remain approximately on the fundamental plane. This indicates that the observed properties of elliptical galaxies are consistent with significant growth via late dissipationless mergers. We also find that the dark matter fraction within the effective radius of our remnants increases after the merger, consistent with the hypothesis that the tilt of the fundamental plane from the virial theorem is due to a varying dark matter fraction as a function of galaxy mass.     

Minor axis rotation and the intrinsic shape of the shell elliptical NGC 3923

We present kinematic observations of NGC 3923, one of the best-known examples of an elliptical galaxy with shells. Although NGC 3923 has no rotation on its major axis, it has minor axis rotation with an amplitude of some 20 km s⁻¹ out to 25 arcsec radius. This is possibily the result of a kinematically decoupled core, but we present arguments against this hypothesis, and if the core were formed by a merger, this must have been a different merger from that which gave rise to the shells. The minor axis rotation suggests that NGC 3923 is prolate or triaxial, at least in the inner region. A prolate geometry for the whole galaxy would favour the merger model over the interaction model for the shell formation mechanism. There is other evidence, however, that the shape of the underlying galaxy changes with radius, in which case it is possible that the minor axis rotation could be confined to the galaxy core. The kinematics could reflect a triaxial mass distribution, in which the long and short axes are aligned with the photometric major and minor axes, and the axial ratios change such that the galaxy is near-prolate in the inner regions and near-oblate at larger radii. Measuring the rotation curve to about double this radius along the minor axis, as well as offset from the nucleus along the major axis at the radius of the innermost shells, would enable the true extent of the minor axis rotation to be determined.     

-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.     

Investigating the Andromeda stream – III. A young shell system in M31

Published maps of red giant stars in the halo region of M31 exhibit a giant stellar stream to the south of this galaxy, as well as a giant 'shelf' to the northeast of M31's centre. Using these maps, we find that there is a fainter shelf of comparable size on the western side as well. By choosing appropriate structural and orbital parameters for an accreting dwarf satellite within the accurate M31 potential model of Geehan et al., we produce a very similar structure in an N -body simulation. In this scenario, the tidal stream produced at pericentre of the satellite's orbit matches the observed southern stream, while the forward continuation of this tidal stream makes up two orbital loops, broadened into fan-like structures by successive pericentric passages; these loops correspond to the north-eastern and western shelves. The tidal debris from the satellite also reproduces a previously observed 'stream' of counterrotating planetary nebulae and a related stream seen in red giant stars. The debris pattern in our simulation resembles the shell systems detected around many elliptical galaxies, though this is the first identification of a shell system in a spiral galaxy and the first in any galaxy close enough to allow measurements of stellar velocities and relative distances. We discuss the physics of these partial shells, highlighting the role played by spatial and velocity caustics in the observations. We show that kinematic surveys of the tidal debris will provide a sensitive measurement of M31's halo potential, while quantifying the surface density of debris in the shelves will let us reconstruct the original mass and time of disruption of the progenitor satellite.     

BV photometry of five shell galaxies

Current views consider shell structures as bona fide signatures of a recent minor/major merging event though also weak interaction models (WIM) could produce long lasting shells on host galaxies possessing a stellar thick disc.We present a B V band photometric study of a sample of 5 shell galaxies belonging to the Malin & Carter (1983) compilation. The structural properties and colors of the galaxies, as well as the colors of their shells are examined in detail. We did not find signatures of the presence of double nuclei. NGC 7585 is the only E galaxy in the sample and has a moderately boxy structure. The other galaxies have either a discy structure or are mixed E/S0 type galaxies. NGC 474 is a true lenticular. NGC 6776 shows a diffuse asymmetric outer structure and a system of tails of the the same color of the galaxy body; but not clear shells. In general, the color of the shells in our sample is similar or slightly redder than that of the host galaxy, whose color, in turn, is typical of the early‐type morphological class. One of the outer shells of NGC 474 is significantly bluer than the body of the galaxy. Since NGC 474 appears to be interacting with NGC 470, the color of this one shell could be explained as result of a recent acquisition of material through tidal interaction. The WIM hypothesis could explain both the red and the blue shells of NGC 474, this latter acquired from the fly‐by of the nearby companion NGC 470, but the lack of the constancy of shell surface brightness as a ratio of the underlying galaxy brightness argues against WIM. We speculate about evidence, which also comes from different observations, that suggests a merging/accretion origin of the shells. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Galaxy density profiles and shapes – I. Simulation pipeline for lensing by realistic galaxy models

Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, which we aim to investigate in a consistent way. In this first paper in a series, we develop a flexible but efficient pipeline to simulate lensing by realistic galaxy models. These galaxy models have separate stellar and dark matter components, each with a range of density profiles and shapes representative of early-type, central galaxies without significant contributions from other nearby galaxies. We use Fourier methods to calculate the lensing properties of galaxies with arbitrary surface density distributions, and Monte Carlo methods to compute lensing statistics such as point-source lensing cross-sections. Incorporating a variety of magnification bias modes lets us examine different survey limitations in image resolution and flux. We rigorously test the numerical methods for systematic errors and sensitivity to basic assumptions. We also determine the minimum number of viewing angles that must be sampled in order to recover accurate orientation-averaged lensing quantities. We find that for a range of non-isothermal stellar and dark matter density profiles typical of elliptical galaxies, the combined density profile and corresponding lensing properties are surprisingly close to isothermal around the Einstein radius. The converse implication is that constraints from strong lensing and/or stellar kinematics, which are indeed consistent with isothermal models near the Einstein radius, cannot trivially be extrapolated to smaller and larger radii.     

QSO lensing magnification associated with galaxy groups

We simulated both the matter and light (galaxy) distributions in a wedge of the Universe and calculated the gravitational lensing magnification caused by the mass along the line-of-sight of galaxies and galaxy groups identified in sky surveys. A large volume redshift cone containing cold dark matter particles mimics the expected cosmological matter distribution in a flat universe with low matter density and a cosmological constant. We generate a mock galaxy catalogue from the matter distribution and identify thousands of galaxy groups in the luminous sky projection. We calculate the expected magnification around galaxies and galaxy groups and then the induced quasi-stellar object (QSO)–lens angular correlation due to magnification bias. This correlation is observable and can be used both to estimate the average mass of the lens population and to make cosmological inferences. We also use analytical calculations and various analyses to compare the observational results with theoretical expectations for the cross-correlation between faint QSOs from the 2dF Survey and nearby galaxies and groups from the Automated Plate Measurement and Sloan Digital Sky Survey Early Data Release. The observed QSO–lens anticorrelations are stronger than the predictions for the cosmological model used. This suggests that there could be unknown systematic errors in the observations and data reduction, or that the model used is not adequate. If the observed signal is assumed to be solely due to gravitational lensing, then the lensing is stronger than expected, due to more massive galactic structures or more efficient lensing than simulated.     

On the dynamics of the satellite galaxies in NGC 5044

The NGC 5044 galaxy group is dominated by a luminous elliptical galaxy that is surrounded by ∼160 dwarf satellites. The projected number density profile of this dwarf population deviates within ∼1/3 of the virial radius from a projected Navarro, Frenk and White (NFW) profile, which is assumed to approximate the underlying total matter distribution. By means of a semi-analytic model, we demonstrate that the interplay between gravitation, dynamical friction and tidal mass loss and destruction can explain the observed number density profile. We use only two parameters in our models: the total to stellar mass fraction of the satellite haloes and the disruption efficiency. The disruption efficiency is expressed by a minimum radius. If the tidal radius of a galaxy (halo) falls below this radius, it is assumed to become unobservable. The preferred parameters are an initial total to stellar mass fraction of ∼20 and a disruption radius of  4 kpc  . In that model, about 20 per cent of all the satellites are totally disrupted on their orbits within the group environment. Dynamical friction is less important in shaping the inner slope of the number density profile because the reduction in mass by tidal forces lowers the impact of the friction term. The main destruction mechanism is tide. In the preferred model, the total B -band luminosity of all disrupted galaxies is about twice the observed luminosity of the central elliptical galaxy, indicating that a significant fraction of stars are scattered into the intragroup medium. Dwarf galaxy satellites closer to the centre of the NGC 5044 group may exhibit optical evidence of partial tidal disruption. If dynamical friction forces the satellite to merge with the central elliptical, the angular momentum of the satellite tends to be removed at the apocentre passage. Afterwards, the satellite drops radially towards the centre.     

Density–potential pairs for spherical stellar systems with Sérsic light profiles and (optional) power-law cores

Popular models for describing the luminosity-density profiles of dynamically hot stellar systems (e.g. Jaffe, Hernquist, Dehnen) were constructed with the desire to match the deprojected form of an   R ^1/4  light profile. Real galaxies, however, are now known to have a range of different light-profile shapes that scale with mass. Consequently, although highly useful, the above models have implicit limitations, and this is illustrated here through their application to a number of real galaxy density profiles. On the other hand, the analytical density profile given by Prugniel & Simien closely matches the deprojected form of Sérsic   R ^{1/ n}   light profiles – including deprojected exponential light profiles. It is thus applicable for describing bulges in spiral galaxies, dwarf elliptical galaxies, and both ordinary and giant elliptical galaxies. Moreover, the observed Sérsic quantities define the parameters of the density model. Here we provide simple equations, in terms of elementary and special functions, for the gravitational potential and force associated with this density profile. Furthermore, to match galaxies with partially depleted cores, and better explore the supermassive black hole/galaxy connection, we have added a power-law core to this density profile and derived similar expressions for the potential and force of this hybrid profile. Expressions for the mass and velocity dispersion, assuming isotropy, are also given. These spherical models may also prove appropriate for describing the dark matter distribution in haloes formed from ΛCDM cosmological simulations.     

From kinematics to dynamics in thin galactic discs

We present a method for recovering the distribution functions of edge-on thin axisymmetric discs directly from their observable kinematic properties. The most generally observable properties of such a stellar system are the line-of-sight velocity distributions of the stars at different projected radii along the galaxy. If the gravitational potential is known, then the general two-integral distribution function can be reconstructed using the shapes of the high-velocity tails of these line-of-sight distributions. If the wrong gravitational potential is adopted, then a distribution function can still be constructed using this technique, but the low-velocity parts of the observed velocity distributions will not be reproduced by the derived dynamical model. Thus, the gravitational potential is also tightly constrained by the observed kinematics.     

Gas dynamic stripping and X-ray emission of cluster elliptical galaxies

Detailed three-dimensional numerical simulations of an elliptical galaxy orbiting in a gas-rich cluster of galaxies indicate that gas dynamic stripping is less efficient than the results from previous, simpler calculations by Takeda et al. and Gaetz et al. implied. This result is consistent with X-ray data for cluster elliptical galaxies. Hydrodynamic torques and direct accretion of orbital angular momentum can result in the formation of a cold gaseous disc, even in a non-rotating galaxy. The gas lost by cluster galaxies via the process of gas dynamic stripping tends to produce a colder, chemically enriched cluster gas core. A comparison of the models with the available X-ray data of cluster galaxies shows that the X-ray luminosity distribution of cluster galaxies may reflect hydrodynamic stripping, but also that a purely hydrodynamic treatment is inadequate for the cooler interstellar medium near the centre of the galaxy.     

The nature of the dwarf population in Abell 868

We present the results of a study of the morphology of the dwarf galaxy population in Abell 868, a rich, intermediate-redshift     cluster which has a galaxy luminosity function (LF) with a steep faint-end slope     . A statistical background subtraction method is employed to study the     colour distribution of the cluster galaxies. This distribution suggests that the galaxies contributing to the faint-end of the measured cluster LF can be split into three populations: dwarf irregular galaxies (dIrrs) with     dwarf elliptical galaxies (dEs) with     and contaminating background giant ellipticals (gEs) with     . The removal of the contribution of the background gEs from the counts only marginally lessens the faint-end slope     . However, the removal of the contribution of the dIrrs from the counts produces a flat LF     . The dEs and the dIrrs have similar spatial distributions within the cluster, except that the dIrrs appear to be totally absent within a central projected radius of about 0.2 Mpc     . The number densities of both dEs and dIrrs appear to fall off beyond a projected radius of ≃ 0.35 Mpc. We suggest that the dE and dIrr populations of A868 have been associated with the cluster for similar time-scales, but evolutionary processes such as 'galaxy harassment' tend to fade the dIrr galaxies while having a much smaller effect on the dE galaxies. The harassment would be expected to have the greatest effect on dwarfs residing in the central parts of the cluster.     

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.     

Imaging of the protoelliptical NGC 1700 and its globular cluster system

An excellent candidate for a young elliptical, or 'protoelliptical' galaxy is NGC 1700. Here we present new B -, V - and I -band imaging using the Keck telescope, and reanalyse existing V - and I -band images from the Hubble Space Telescope . After subtracting a model of the galaxy from the Keck images, NGC 1700 reveals two symmetric tidal tail-like structures. If this interpretation is correct, it suggests a past merger event involving two spiral galaxies. These tails are largely responsible for the 'boxiness' of the galaxy isophotes observed at a radius of ∼13 kpc.
We also show that the B − I colour distribution of the globular cluster system is bimodal. The mean colour of the blue population is consistent with that of old Galactic globular clusters. Relative to this old, metal-poor population, we find that the red population is younger and more metal-rich. This young population has an age and metallicity similar to that inferred for the central stars, suggesting that both populations are associated with an episode of star formation triggered by the merger that may have formed the galaxy. We find that, although they have large errors, the majority of the age estimates of NGC 1700 are reasonably consistent and we adopt a 'best estimate' for the age of 3.0±1.0 Gyr. This relatively low age places NGC 1700 within the age range where there is a notable lack of obvious candidates for protoellipticals. The total globular cluster specific frequency is rather low for a typical elliptical, even after taking into account the fading of the galaxy over the next 10 Gyr. We speculate that NGC 1700 will eventually form a relatively 'globular cluster poor' elliptical galaxy.     

The cosmological evolution of colour gradients in spheroids

The analysis of the four-colour maps of galaxies in the Hubble Deep Field (HDF) has revealed, in the sample of 0.4< z <1 early-type field galaxies, the existence of ellipticals with a predominantly old coeval stellar population. However, there is another, unexpected, category of HDF early-type galaxies, in which the galaxy core is significantly bluer than the outer regions. We demonstrate that these colour gradients are predicted by the multizone chemodynamical model for the evolution of elliptical galaxies.
We suggest that the colour gradient could be used as a chronometer for the evolution of elliptical galaxies: galaxies younger than a few Gyr exhibit cores bluer than the surrounding galaxy as a result of ongoing star formation, while more evolved galaxies have redder cores because of metallicity gradients increasing toward the centre.     

Galaxy aggregates in the Coma cluster

We present evidence for a new morphologically defined form of small-scale substructure in the Coma cluster, which we call galaxy aggregates. Aggregates are dominated by a central galaxy, which is on average 5 mag brighter than the smaller aggregate members, nearly all of which lie to one side of the central galaxy. We have found three such galaxy aggregates: two dominated by the S0 galaxies RB 55 and RB 60, and one by the starbursting SBb NGC 4858.   RB 55 and 60 are both equidistant between the two dominant D galaxies NGC 4874 and 4889, while NGC 4858 is located near the large E0 galaxy NGC 4860. All three central galaxies have redshifts consistent with Coma cluster membership. We describe the spatial structures of these unique objects, and suggest several possible mechanisms to explain their origin. These include: chance superpositions from background galaxies, interactions between other galaxies and with the cluster gravitational potential, and ram pressure. We conclude that the most probable scenario of creation is an interaction with the cluster through its gravitational potential.     

Morphology of flows and buoyant bubbles in the Virgo cluster

There is growing evidence that the active galactic nuclei (AGN) associated with the central elliptical galaxy in clusters of galaxies are playing an important role in the evolution of the intracluster medium (ICM) and clusters themselves. We use high-resolution three-dimensional simulations to study the interaction of the cavities created by AGN outflows (bubbles) with the ambient ICM. The gravitational potential of the cluster is modelled using the observed temperature and density profiles of the Virgo cluster. We demonstrate the importance of the hydrodynamical Kutta–Zhukovsky forces associated with the vortex ring structure of the bubbles, and discuss possible effects of diffusive processes on their evolution.     

The effect of dwarf galaxy disruption in semi-analytic models

We present results for a galaxy formation model that includes a simple treatment for the disruption of dwarf galaxies by gravitational forces and galaxy encounters within galaxy clusters. This is implemented a posteriori in a semi-analytic model by considering the stability of cluster dark matter subhaloes at   z = 0  . We assume that a galaxy whose dark matter substructure has been disrupted will itself disperse, while its stars become part of the population of intracluster stars responsible for the observed intracluster light. Despite the simplicity of this assumption, our results show a substantial improvement over previous models and indicate that the inclusion of galaxy disruption is indeed a necessary ingredient of galaxy formation models. We find that galaxy disruption suppresses the number density of dwarf galaxies by about a factor of 2. This makes the slope of the faint end of the galaxy luminosity function shallower, in agreement with observations. In particular, the abundance of faint, red galaxies is strongly suppressed. As a result, the luminosity function of red galaxies and the distinction between the red and the blue galaxy populations in colour–magnitude relationships are correctly predicted. Finally, we estimate a fraction of intracluster light comparable to that found in clusters of galaxies.     

Stellar velocity profiles and line strengths out to four effective radii in the early-type galaxies NGC 3379 and 821

We use the integral-field spectrograph SAURON to measure the stellar line-of-sight velocity distribution and absorption line strengths out to four effective radii ( R _e) in the early-type galaxies NGC 3379 and 821. With our newly developed observing technique, we can now probe these faint regions in galaxies that were previously not accessible with traditional long-slit spectroscopy. We make optimal use of the large field-of-view and high throughput of the spectrograph: by adding the signal of all ∼1400 lenslets into one spectrum, we obtain sufficient signal-to-noise in a few hours of observing time to reliably measure the absorption line kinematics and line strengths out to large radius.
We find that the line strength gradients previously observed within 1 R _e remain constant out to at least 4 R _e, which puts constraints on the merger histories of these galaxies. The stellar halo populations are old and metal poor. By constructing orbit-based Schwarzschild dynamical models, we find that dark matter is necessary to explain the observed kinematics in NGC 3379 and 821, with 30–50 per cent of the total matter being dark within 4 R _e. The radial anisotropy in our best-fitting halo models is less than in our models without halo, due to differences in orbital structure. The halo also has an effect on the  Mg  b – V _esc  relation: its slope is steeper when a dark matter halo is added to the model.