Massive dark matter haloes around bright isolated galaxies in the 2dFGRS

We identify a large sample of isolated bright galaxies and their fainter satellites in the 2dF Galaxy Redshift Survey (2dFGRS). We analyse the dynamics of ensembles of these galaxies selected according to luminosity and morphological type by stacking the positions of their satellites and estimating the velocity dispersion of the combined set. We test our methodology using realistic mock catalogues constructed from cosmological simulations. The method returns an unbiased estimate of the velocity dispersion provided that the isolation criterion is strict enough to avoid contamination and that the scatter in halo mass at fixed primary luminosity is small. Using a maximum likelihood estimator that accounts for interlopers, we determine the satellite velocity dispersion within a projected radius of 175  h ⁻¹ kpc. The dispersion increases with the luminosity of the primary and is larger for elliptical galaxies than for spiral galaxies of similar b _J luminosity. Calibrating the mass–velocity dispersion relation using our mock catalogues, we find a dynamical mass within 175  h ⁻¹ kpc of     for elliptical galaxies and     for spiral galaxies. Finally, we compare our results with recent studies and investigate their limitations using our mock catalogues.     

Black hole demographics from the relation

We analyse a sample of 32 galaxies for which a dynamical estimate of the mass of the hot stellar component, M _bulge, is available. For each of these galaxies, we calculate the mass of the central black hole, M _•, using the tight empirical correlation between M _• and bulge stellar velocity dispersion. The frequency function     is reasonably well described as a Gaussian with     and standard deviation ∼0.45; the implied mean ratio of black hole mass to bulge mass is a factor of ∼5 smaller than generally quoted in the literature. We present marginal evidence for a lower, average black hole mass fraction in more massive galaxies. The total mass density in black holes in the local Universe is estimated to be ∼     consistent with that inferred from high-redshift     active galactic nuclei.     

Late superhumps and the stream–disc impact in IY UMa

We present two simple dynamical models for Sagittarius based on N -body simulations of the progressive disruption of a satellite galaxy orbiting for 12.5 Gyr within a realistic Galactic potential. In both models the satellite initially has observable properties similar to those of current outlying dwarfs; in one case it is purely stellar while in the other it is embedded in an extended massive halo. The purely stellar progenitor is a King model with a total velocity dispersion of 18.9 km s⁻¹, a core radius of 0.44 kpc and a tidal radius of 3 kpc. The initial stellar distribution in the other case follows a King profile with the same core radius, a slightly larger total velocity dispersion and similar extent. Both these models are consistent with all published data on the current Sagittarius system, they match not only the observed properties of the main body of Sagittarius, but also those reported for unbound debris at larger distances.     

Self-consistent response of a galactic disc to vertical perturbations

We study the self-consistent, linear response of a galactic disc to vertical perturbations, as induced, say, by a tidal interaction. We calculate the self-gravitational potential corresponding to a non-axisymmetric, self-consistent density response of the disc using the Green's function method. The response potential is shown to oppose the perturbation potential because the self-gravity of the disc resists the imposed potential, and this resistance is stronger in the inner parts of a galactic disc. For the   m = 1  azimuthal wavenumber, the disc response opposes the imposed perturbation up to a radius that spans a range of 4–6 disc scalelengths, so that the disc shows a net warp only beyond this region. This physically explains the well known but so far unexplained observation that warps typically set in beyond this range of radii. We show that the inclusion of a dark matter halo in the calculation only marginally changes (by ∼10 per cent) the radius for the onset of warps. For perturbations with higher azimuthal wavenumbers, the net signature of the vertical perturbations can only be seen at larger radii – for example, beyond 7 exponential disc scalelengths for   m = 10  . Also, for the high- m cases, the magnitude of the negative disc response due to the disc self-gravity is much smaller. This is shown to result in corrugations of the mid-plane density, which explains the puzzling scalloping with   m = 10  detected in H  i in the outermost regions ∼30 kpc in the Galaxy.     

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.     

Investigating the Andromeda stream – II. Orbital fits and properties of the progenitor

We construct test-particle orbits and simple N -body models that match the properties of the giant stellar stream observed to the south of M31, using the model of M31's potential derived in the companion paper by Geehan et al. We introduce a simple approximation to account for the difference in position between the stream and the orbit of the progenitor; this significantly affects the best-fitting orbits. The progenitor orbits we derive have orbital apocentre  ∼60 kpc  and pericentre  ∼3 kpc  , though these quantities vary somewhat with the current orbital phase of the progenitor which is as yet unknown. Our best combined fit to the stream and galaxy properties implies a mass within 125 kpc of M31 of  (7.4 ± 1.2) × 10¹¹ M_⊙  . Based on its length, width, luminosity, and velocity dispersion, we conclude that the stream originates from a progenitor satellite with mass   M _s∼ 10⁹ M_⊙  , and at most modest amounts of dark matter; the estimate of M _s is again correlated with the phase of the progenitor. M31 displays a large number of faint features in its inner halo which may be progenitors or continuations of the stream. While the orbital fits are not constrained enough for us to conclusively identify the progenitor, we can identify several plausible candidates, of which a feature in the planetary nebula distribution found by Merrett et al. is the most plausible, and rule out several others. We make predictions for the kinematic properties of the successful candidates. These may aid in observational identification of the progenitor object, which would greatly constrain the allowed models of the stream.     

The extended rotation curve and the dark matter halo of M33

We present the 21-cm rotation curve of the nearby galaxy M33 out to a galactocentric distance of 16 kpc (13 disc scalelengths). The rotation curve keeps rising out to the last measured point and implies a dark halo mass ≳5×10¹⁰ M_⊙. The stellar and gaseous discs provide virtually equal contributions to the galaxy gravitational potential at large galactocentric radii, but no obvious correlation is found between the radial distribution of dark matter and the distribution of stars or gas.
Results of the best fit to the mass distribution in M33 picture a dark halo which controls the gravitational potential from 3 kpc outward, with a matter density which decreases radially as R ^−1.3. The density profile is consistent with the theoretical predictions for structure formation in hierarchical clustering cold dark matter (CDM) models, and favours lower mass concentrations than those expected in the standard cosmogony.     

The mass of the Andromeda galaxy

This paper argues that the Milky Way galaxy is probably the largest member of the Local Group. The evidence comes from estimates of the total mass of the Andromeda galaxy (M31) derived from the three-dimensional positions and radial velocities of its satellite galaxies, as well as the projected positions and radial velocities of its distant globular clusters and planetary nebulae. The available data set comprises 10 satellite galaxies, 17 distant globular clusters and nine halo planetary nebulae with radial velocities. We find that the halo of Andromeda has a mass of together with a scalelength of 90 kpc and a predominantly isotropic velocity distribution. For comparison, our earlier estimate for the Milky Way halo is Although the error bars are admittedly large, this suggests that the total mass of M31 is probably less than that of the Milky Way . We verify the robustness of our results to changes in the modelling assumptions and to errors caused by the small size and incompleteness of the data set.
Our surprising claim can be checked in several ways in the near future. The numbers of satellite galaxies, planetary nebulae and globular clusters with radial velocities can be increased by ground-based spectroscopy, while the proper motions of the companion galaxies and the unresolved cores of the globular clusters can be measured using the astrometric satellites Space Interferometry Mission ( SIM ) and Global Astrometric Interferometer for Astrophysics ( GAIA ). Using 100 globular clusters at projected radii 20 R 50 kpc with both radial velocities and proper motions, it will be possible to estimate the mass within 50 kpc to an accuracy of 20 per cent. Measuring the proper motions of the companion galaxies with SIM and GAIA will reduce the uncertainty in the total mass caused by the small size of the data set to 22 per cent.     

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.     

Tracing spiral density waves in M81

We use Spitzer IRAC 3.6 and 4.5 μm near-infrared data from the Spitzer Infrared Nearby Galaxies Survey (SINGS), optical B, V and I and Two-Micron All-Sky Survey K _s-band data to produce mass surface density maps of M81. The IRAC 3.6- and 4.5-μm data, whilst dominated by emission from old stellar populations, are corrected for small-scale contamination by young stars and polycyclic aromatic hydrocarbon emission. The I -band data are used to produce a mass surface density map by a   B − V   colour correction, following the method of Bell and de Jong. We fit a bulge and exponential disc to each mass map, and subtract these components to reveal the non-axisymmetric mass surface density. From the residual mass maps, we are able to extract the amplitude and phase of the density wave, using azimuthal profiles. The response of the gas is observed via dust emission in the 8-μm IRAC band, allowing a comparison between the phase of the stellar density wave and gas shock. The relationship between this angular offset and radius suggests that the spiral structure is reasonably long-lived and allows the position of corotation to be determined.     

Chandra measurements of the distribution of mass in the luminous lensing cluster Abell 2390

We present spatially resolved X-ray spectroscopy of the luminous lensing cluster Abell 2390, using observations made with the Chandra observatory. The temperature of the X-ray gas rises with increasing radius within the central ∼ 200 kpc of the cluster, and then remains approximately isothermal, with kT =11.5_−1.6^+1.5 keV , out to the limits of the observations at r ∼1.0 Mpc . The total mass profile determined from the Chandra data has a form in good agreement with the predictions from numerical simulations. Using the parametrization of Navarro, Frenk and White, we measure a scale radius r _s∼0.8 Mpc and a concentration parameter c ∼3 . The best-fitting X-ray mass model is in good agreement with independent gravitational lensing results and optical measurements of the galaxy velocity dispersion in the cluster. The X-ray gas to total mass ratio rises with increasing radius with f _gas∼21 per cent at r =0.9 Mpc . The azimuthally averaged 0.3–7.0 keV surface brightness profile exhibits a small core radius and a clear 'break' at r ∼500 kpc , where the slope changes from S _X ∼^∝  r ^−1.5 to S _X ∼^∝  r ^−3.6 . The data for the central region of the cluster indicate the presence of a cooling flow with a mass deposition rate of 200–300 M_⊙ yr⁻¹ and an effective age of 2–3 Gyr .     

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.     

Hubble Space Telescope photometry of globular clusters in the Sombrero galaxy

We explore the rich globular cluster (GC) system of the nearby Sa galaxy M104, the 'Sombrero' (NGC 4594), using archive Wide Field Planetary Camera 2 data. The GC colour distribution is found to be bimodal at the >99 per cent confidence level, with peaks at     and     . The inferred metallicities are very similar to those of GCs in our Galaxy and M31. However, the Sombrero reveals a much enhanced number of red (metal-rich) GCs compared to other well-studied spirals. Because the Sombrero is dominated by a huge bulge and only has a modest disc, we associate the two subpopulations with the halo and bulge components, respectively. Thus our analysis supports the view that the metal-rich GCs in spirals are associated with the bulge rather than with the disc. The Sombrero GCs have typical effective (half-light) radii of ∼2 pc with the red ones being ∼30 per cent smaller than the blue ones. We identify many similarities between the GC system of the Sombrero and those of both late-type spirals and early-type galaxies. Thus both the GC system and the Hubble type of the Sombrero galaxy appear to be intermediate in their nature.     

GMOS spectroscopy of the S0 galaxy NGC 3115

We present Gemini Multi-Object Spectrograph longslit spectroscopy of the isolated S0 galaxy NGC 3115. We have determined kinematical data and Lick/IDS absorption line-strength indices for the major axis out to around 9 kpc and for the minor axis out to around 5 kpc (around 2R _e ). Using stellar population models which include the effects of variable [α/Fe] ratios, we derive metallicities, abundance ratios and ages for the stellar population of NGC 3115. We find that [α/Fe] remains fairly constant with increasing radius at around  [α/Fe]= 0.17  for the major axis but increases rapidly for the minor axis to around  [α/Fe]= 0.3  . We also find that to first order, this behaviour can be explained by a simple spheroid + disc model, where the spheroid has  [α/Fe]= 0.3  and the disc shows close to solar abundance ratios. The disc also appears considerably younger than the spheroid, having an age of around 6 Gyr compared to 12 Gyr for the spheroid. We compare these results to those previously presented for the globular cluster system of NGC 3115.     

Modelling the bar in the centre of the starburst galaxy M82

We present VLA A-array 21-cm atomic hydrogen (H  i ) absorption observed against the central region of the starburst galaxy M82 with an angular resolution of ∼1.3 arcsec (≃20 pc). These observations, together with MERLIN H  i absorption measurements, are compared with the molecular (CO) and ionized ([Ne  ii ]) gas distributions and are used to constrain the dynamics and structure of the ionized, neutral and molecular gas in this starburst.
A position–velocity diagram of the H  i distribution reveals an unusual 'hole' feature which, when previously observed in CO, has been interpreted as an expanding superbubble contained within a ring of gas in solid body rotation. However, we interpret this feature as a signature of a nearly edge-on barred galaxy. In addition, we note that the CO, H  i and [Ne  ii ] position–velocity diagrams reveal two main velocity gradients, and we interpret these as gas moving on x₁- and x₂-orbits within a bar potential. We find the best fit to the data to be produced using a bar potential with a flat rotation curve velocity v _b=140 km s⁻¹ and a total length of 1 kpc, a non-axisymmetry parameter q =0.9, an angular velocity of the bar Ω_b=217 km s⁻¹ arcsec⁻¹, a core radius R _c=25 pc, an inclination angle i =80° and a projected angle between the bar and the major axis of the galaxy φ '=4°. We also discuss the orientation of the disc and bar in M82.     

A new superwind galaxy: XMM–Newton observations of NGC 6810

We present the first imaging X-ray observation of the highly inclined  ( i = 78°)  Sab Seyfert 2 galaxy NGC 6810 using XMM–Newton , which reveals soft X-ray emission that extends out to a projected height of ∼7 kpc away from the plane of the galaxy. The soft X-ray emission beyond the optical disc of the galaxy is most plausibly extraplanar, although it could instead come from large galactic radius. This extended X-ray emission is spatially associated with diffuse Hα emission, in particular with a prominent 5-kpc-long Hα filament on the north-west of the disc. A fraction ≲35 per cent of the total soft X-ray emission of the galaxy arises from projected heights  | z | ≥ 2 kpc  . Within the optical disc of the galaxy the soft X-ray emission is associated with the star-forming regions visible in ground-based Hα and XMM–Newton optical monitor near-UV imaging. The temperature, supersolar α-element-to-iron abundance ratio, soft X-ray/Hα correlation, and X-ray to far-infrared (FIR) flux ratio of NGC 6810 are all consistent with local starbursts with winds, although the large base radius of the outflow would make NGC 6810 one of the few 'disc-wide' superwinds currently known. Hard X-ray emission from NGC 6810 is weak, and the total   E = 2–10 keV  luminosity and spectral shape are consistent with the expected level of X-ray binary emission from the old and young stellar populations. The X-ray observations provide no evidence of any active galactic nucleus activity. We find that the optical, IR and radio properties of NGC 6810 are all consistent with a starburst galaxy, and that the old classification of this galaxy as a Seyfert 2 galaxy is probably incorrect.     

A radio study of neutral gas and dust in the radio galaxy 3C 293

We have observed broad H  i absorption in the radio galaxy 3C 293 using Multi-Element Radio Linked Interferometric Network (MERLIN) at 0.2-arcsec angular resolution and the Giant Meterwavelength Radio Telescope (GMRT) at arcsec resolution. Extensive H  i is found in absorption across the centre of this peculiar radio galaxy, allowing a detailed study of the dynamics of the neutral gas on linear scales down to ∼160 pc. In optical depth position–velocity diagrams across the central few kpc we detect a distinct velocity gradient of 179 km s⁻¹ arcsec⁻¹ associated with the broad absorption. This is interpreted as a ring of neutral gas rotating around the suspected position of the active galactic nucleus (AGN) . The radius of this high velocity gradient ring is found to be >0.74 arcsec (600 pc), implying an upper limit upon the enclosed mass of     , assuming a near edge-on disc with an inclination of i . The optical depth of H  i is mapped across the entire central region of 3C 293 showing enhancements of a factor of 4 in the areas that are co-spatial with dust lanes seen in Hubble Space Telescope ( HST ) imaging of this galaxy.     

Unravelling the mystery of the M31 bar

The inclination of M31 is too close to edge-on for a bar component to be easily recognized and is not sufficiently edge-on for a boxy/peanut bulge to protrude clearly out of the equatorial plane. Nevertheless, a sufficient number of clues allow us to argue that this galaxy is barred. We use fully self-consistent N -body simulations of barred galaxies and compare them with both photometric and kinematic observational data for M31. In particular, we rely on the near-infrared photometry presented in a companion paper. We compare isodensity contours to isophotal contours and the light profile along cuts parallel to the galaxy major axis and offset towards the north, or the south, to mass profiles along similar cuts on the model. All these comparisons, as well as position–velocity diagrams for the gaseous component, give us strong arguments that M31 is barred. We compare four fiducial N -body models to the data and thus set constraints on the parameters of the M31 bar, as its strength, length and orientation. Our 'best' models, although not meant to be exact models of M31, reproduce in a very satisfactory way the main relevant observations. We present arguments that M31 has both a classical and a boxy/peanut bulge. Its pseudo-ring-like structure at roughly 50 arcmin is near the outer Lindblad resonance of the bar and could thus be an outer ring, as often observed in barred galaxies. The shape of the isophotes also argues that the vertically thin part of the M31 bar extends considerably further out than its boxy bulge, that is, that the boxy bulge is only part of the bar, thus confirming predictions from orbital structure studies and from previous N -body simulations. It seems very likely that the backbone of M31's boxy bulge is families of periodic orbits, members of the x₁-tree and bifurcating from the x₁ family at its higher order vertical resonances, such as the x1v3 or x1v4 families.     

Near-infrared bulge–disc correlations of lenticular galaxies

We consider the luminosity and environmental dependence of structural parameters of lenticular galaxies in the near-infrared K band. Using a 2D galaxy image decomposition technique, we extract bulge and disc structural parameters for a sample of 36 lenticular galaxies observed by us in the K band. By combining data from the literature for field and cluster lenticulars with our data, we study correlations between parameters that characterize the bulge and the disc as a function of luminosity and environment. We find that scaling relations such as the Kormendy relation, photometric plane and other correlations involving bulge and disc parameters show a luminosity dependence. This dependence can be explained in terms of galaxy formation models in which faint lenticulars  ( M _T > −24.5)  formed via secular formation processes that likely formed the pseudo-bulges of late-type disc galaxies, while brighter lenticulars  ( M _T < −24.5)  formed through a different formation mechanism most likely involving major mergers. On probing variations in lenticular properties as a function of environment, we find that faint cluster lenticulars show systematic differences with respect to faint field lenticulars. These differences support the idea that the bulge and disc components fade after the galaxy falls into a cluster, while simultaneously undergoing a transformation from spiral to lenticular morphologies.