Modelling gaseous and stellar kinematics in the disc galaxies NGC 772, 3898 and 7782

We present V -band surface photometry and major-axis kinematics of stars and ionized gas of three early-type spiral galaxies, namely NGC 772, 3898 and 7782. For each galaxy we present a self-consistent Jeans model for the stellar kinematics, adopting the light distribution of bulge and disc derived by means of a two-dimensional parametric photometric decomposition. This allows us to investigate the presence of non-circular gas motions, and derive the mass distribution of luminous and dark matter in these objects.
NGC 772 and 7782 have apparently normal kinematics with the ionized gas tracing the gravitational equilibrium circular speed. This is not true in the innermost region (| r |≲8 arcsec) of NGC 3898, where the ionized gas is rotating more slowly than the circular velocity predicted by dynamical modelling. This phenomenon is common in the bulge-dominated galaxies for which dynamical modelling enables us to make the direct comparison between the gas velocity and the circular speed, and it poses questions about the reliability of galaxy mass distributions derived by the direct decomposition of the observed ionized-gas rotation curve into the contributions of luminous and dark matter.     

The circumnuclear material in the Galactic Centre: a clue to the accretion process

On the basis of 'sticky particle' calculations, it is argued that the gas features observed within 10 pc of the Galactic Centre — the circumnuclear disc (CND) and the ionized gas filaments — as well as the newly formed stars in the inner 1 pc can be understood in terms of tidal capture and disruption of gas clouds on low angular momentum orbits in a potential containing a point mass. The calculations demonstrate that a dissipative component forms a 'dispersion ring', an asymmetric elliptical torus precessing counter to the direction of rotation, and that this shape can be maintained for many orbital periods. For a range of plausible initial conditions, such a structure can explain the morphology and kinematics of the CND and of the most conspicuous ionized filament. While forming the dispersion ring, a small cloud with low specific angular momentum is drawn into a long filament which repeatedly collides with itself at high velocity. The compression in strong shocks is likely to lead to star formation even in the near tidal field of the point mass. This process may have general relevance to accretion on to massive black holes in normal and active galactic nuclei.     

The SAURON project – VII. Integral-field absorption and emission-line kinematics of 24 spiral galaxy bulges

We present observations of the stellar and gas kinematics for a representative sample of 24 Sa galaxies obtained with our custom-built integral-field spectrograph SAURON operating on the William Herschel Telescope. The data have been homogeneously reduced and analysed by means of a dedicated pipeline. All resulting data cubes were spatially binned to a minimum mean signal-to-noise ratio of 60 per spatial and spectral resolution element. Our maps typically cover the bulge-dominated region. We find a significant fraction of kinematically decoupled components (12/24), many of them displaying central velocity dispersion minima. They are mostly aligned and co-rotating with the main body of the galaxies, and are usually associated with dust discs and rings detected in unsharp-masked images. Almost all the galaxies in the sample (22/24) contain significant amounts of ionized gas which, in general, is accompanied by the presence of dust. The kinematics of the ionized gas are consistent with circular rotation in a disc co-rotating with respect to the stars. The distribution of mean misalignments between the stellar and gaseous angular momenta in the sample suggests that the gas has an internal origin. The [O  iii ]/Hβ ratio is usually very low, indicative of current star formation, and shows various morphologies (ring-like structures, alignments with dust lanes or amorphous shapes). The star formation rates (SFRs) in the sample are comparable with that of normal disc galaxies. Low gas velocity dispersion values appear to be linked to regions of intense star formation activity. We interpret this result as stars being formed from dynamically cold gas in those regions. In the case of NGC 5953, the data suggest that we are witnessing the formation of a kinematically decoupled component from cold gas being acquired during the ongoing interaction with NGC 5954.     

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 dynamical model for the extraplanar gas in spiral galaxies

Recent H  i observations reveal that the discs of spiral galaxies are surrounded by extended gaseous haloes. This extraplanar gas reaches large distances (several kiloparsecs) from the disc and shows peculiar kinematics (low rotation and inflow). We have modelled the extraplanar gas as a continuous flow of material from the disc of a spiral galaxy into its halo region. The output of our models is pseudo data cubes that can be directly compared to the H  i data. We have applied these models to two spiral galaxies (NGC 891 and NGC 2403) known to have a substantial amount of extraplanar gas. Our models are able to reproduce accurately the vertical distribution of extraplanar gas for an energy input corresponding to a small fraction (<4 per cent) of the energy released by supernovae. However, they fail in two important aspects: (1) they do not reproduce the right gradient in rotation velocity; (2) they predict a general outflow of the extraplanar gas, contrary to what is observed. We show that neither of these difficulties can be removed if clouds are ionized and invisible at 21 cm as they leave the disc but become visible at some point on their orbits. We speculate that these failures indicate the need for accreted material from the intergalactic medium that could provide the low angular momentum and inflow required.     

Neutral hydrogen absorption observations of the central region of NGC 5929

We present 0.15-arcsec (25-pc) resolution MERLIN observations of neutral hydrogen absorption detected towards the nuclear region of the type 2 Seyfert galaxy NGC 5929. Absorption is detected only towards the north-eastern radio component with a column density of (6.5 ± 0.6) × 10²¹ cm⁻². Based on comparison with an HST WFPC2 continuum image, we propose that the absorption is caused by a 1.5-arcsec structure of neutral gas and dust offset 0.3 arcsec south-east of the nucleus and running NE–SW. A separate cloud of dust is apparent 1.5 arcsec to the south-west of the nucleus in the HST image. A comparison of the centroid velocity (2358 ± 5 km s⁻¹) and full width at half-maximum (43 ± 6 km s⁻¹) of the absorbing gas with previous [O  III ] observations suggests that both the neutral and ionized gas are undergoing galactic rotation towards the observer in the north-east and away from the observer in the south-west. The main structure is consistent with an inclined ring of gas and dust encircling the active galactic nucleus (AGN); alternatively it may be a bar or inner spiral arm. We do not detect neutral hydrogen absorption or dust obscuration against the radio nucleus (column density < 3.1 × 10²¹ cm⁻²) expected by a torus of neutral gas and dust in unified models of AGNs for a type 2 Seyfert galaxy.     

Late-type galaxies observed with SAURON: two-dimensional stellar and emission-line kinematics of 18 spirals

We present the stellar and gas kinematics of a sample of 18 nearby late-type spiral galaxies (Hubble types ranging from Sb to Sd), observed with the integral-field spectrograph SAURON at the 4.2-m William Herschel Telescope. SAURON covers the spectral range 4800–5380 Å, allowing us to measure the Hβ, Fe, Mg b absorption features and the emission in the Hβ line and the  [O  iii ]λλ4959, 5007 Å  and  [N  i ]λλ5198, 5200 Å  doublets over a  33 × 41-arcsec²  field of view. The maps cover the nuclear region of these late-type galaxies and in all cases include the entire bulge. In many cases the stellar kinematics suggests the presence of a cold inner region, as visible from a central drop in the stellar velocity dispersion. The ionized gas is almost ubiquitous and behaves in a complicated fashion: the gas velocity fields often display more features than the stellar ones, including wiggles in the zero-velocity lines, irregular distributions, ring-like structures. The line ratio [O  iii ]/Hβ often takes on low values over most of the field, probably indicating a wide-spread star formation.     

The Planetary Nebula Spectrograph elliptical galaxy survey: the dark matter in NGC 4494

We present new Planetary Nebula Spectrograph observations of the ordinary elliptical galaxy NGC 4494, resulting in positions and velocities of 255 planetary nebulae out to seven effective radii (25 kpc). We also present new wide-field surface photometry from MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and kinematical distributions of the planetary nebulae agree with the field stars in the region of overlap. The mean rotation is relatively low, with a possible kinematic axis twist outside  1 R _e  . The velocity dispersion profile declines with radius, though not very steeply, down to  ∼70 km s⁻¹  at the last data point.
We have constructed spherical dynamical models of the system, including Jeans analyses with multi-component Λ cold dark matter (CDM) motivated galaxies as well as logarithmic potentials. These models include special attention to orbital anisotropy, which we constrain using fourth-order velocity moments. Given several different sets of modelling methods and assumptions, we find consistent results for the mass profile within the radial range constrained by the data. Some dark matter (DM) is required by the data; our best-fitting solution has a radially anisotropic stellar halo, a plausible stellar mass-to-light ratio and a DM halo with an unexpectedly low central density. We find that this result does not substantially change with a flattened axisymmetric model.
Taken together with other results for galaxy halo masses, we find suggestions for a puzzling pattern wherein most intermediate-luminosity galaxies have very low concentration haloes, while some high-mass ellipticals have very high concentrations. We discuss some possible implications of these results for DM and galaxy formation.     

Radio planetary nebulae in the Magellanic Clouds

We report the extragalactic radio-continuum detection of 15 planetary nebulae (PNe) in the Magellanic Clouds (MCs) from recent Australia Telescope Compact Array+Parkes mosaic surveys. These detections were supplemented by new and high-resolution radio, optical and infrared observations which helped to resolve the true nature of the objects. Four of the PNe are located in the Small Magellanic Cloud (SMC) and 11 are located in the Large Magellanic Cloud (LMC). Based on Galactic PNe the expected radio flux densities at the distance of the LMC/SMC are up to ∼2.5 and ∼2.0 mJy at 1.4 GHz, respectively. We find that one of our new radio PNe in the SMC has a flux density of 5.1 mJy at 1.4 GHz, several times higher than expected. We suggest that the most luminous radio PN in the SMC (N S68) may represent the upper limit to radio-peak luminosity because it is approximately three times more luminous than NGC 7027, the most luminous known Galactic PN. We note that the optical diameters of these 15 Magellanic Clouds (MCs) PNe vary from very small (∼0.08 pc or 0.32 arcsec; SMP L47) to very large (∼1 pc or 4 arcsec; SMP L83). Their flux densities peak at different frequencies, suggesting that they may be in different stages of evolution. We briefly discuss mechanisms that may explain their unusually high radio-continuum flux densities. We argue that these detections may help solve the 'missing mass problem' in PNe whose central stars were originally  1–8 M_⊙  . We explore the possible link between ionized haloes ejected by the central stars in their late evolution and extended radio emission. Because of their higher than expected flux densities, we tentatively call this PNe (sub)sample –'Super PNe'.     

The shape of the dark matter halo in the early-type galaxy NGC 2974

We present H  i observations of the elliptical galaxy NGC 2974, obtained with the Very Large Array. These observations reveal that the previously detected H  i disc in this galaxy is in fact a ring. By studying the harmonic expansion of the velocity field along the ring, we constrain the elongation of the halo and find that the underlying gravitational potential is consistent with an axisymmetric shape.
We construct mass models of NGC 2974 by combining the H  i rotation curve with the central kinematics of the ionized gas, obtained with the integral-field spectrograph SAURON. We introduce a new way of correcting the observed velocities of the ionized gas for asymmetric drift, and hereby disentangle the random motions of the gas caused by gravitational interaction from those caused by turbulence. To reproduce the observed flat rotation curve of the H  i gas, we need to include a dark halo in our mass models. A pseudo-isothermal sphere provides the best model to fit our data, but we also tested an NFW halo and modified Newtonian dynamics, which fit the data marginally worse.
The mass-to-light ratio M / L _I increases in NGC 2974 from 4.3 M_⊙/L_⊙, _I at one effective radius to 8.5 M_⊙/L_{⊙, I} at 5  R _e. This increase of M / L already suggests the presence of dark matter: we find that within 5  R _e at least 55 per cent of the total mass is dark.     

The kinematics of lopsided galaxies

Lopsidedness is a common feature in galaxies, both in the distribution of light and in the kinematics. We investigate the kinematics of a model for lopsided galaxies that consists of a disc lying off-centre in a dark halo, and circling around the halo centre. We search for families of stable, closed, non-crossing orbits, and assume that gas in our galaxies moves on these orbits. Several of our models show strong lopsided gas kinematics, especially those in which the disc spins around its axis in a retrograde sense compared with its motion around the halo centre. We are able to reproduce the H  i velocity map of the kinematically lopsided galaxy NGC 4395.
The lopsidedness in our models is most pronounced in the models where the halo provides a relatively large fraction of the total mass at small radii. This may explain why the gas shows lopsidedness more frequently in late-type galaxies, which are dominated by dark matter. Surfaces of section show large regions of irregular orbits in the models where the halo density is low. This may indicate that these models are unstable.     

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.     

A catalogue of potential adaptive optics survey fields from the UKIRT archive

We present a multicolour catalogue of faint galaxies situated close to bright stars,   V ≲ 15  , with the aim of identifying high-redshift galaxies suitable for study with adaptive optics-equipped near-infrared imagers and spectrographs. The catalogue is constructed from archival calibration observations of the United Kingdom Infrared Telescope (UKIRT) Faint Standard stars with the UKIRT Fast Track Imager (UFTI) camera on UKIRT. We have analysed the deepest 16 fields from the archive to provide a catalogue of galaxies brighter than   K ∼ 20.3  lying between 3 and 25 arcsec of the guide stars. We identify 111 objects in a total survey area of  8.7 arcmin²  . Of these, 87 are classified as galaxies based on their light profiles in our ∼0.5 arcsec median seeing K -band images. 12 of the galaxies have  ( J − K ) ≥ 2.0  consistent with them lying at high redshifts,   z ≳ 2  . These 12 very red galaxies have K -band magnitudes of   K = 18.1–20.1  and separations from the guide stars of 4–20 arcsec and hence are very well suited to adaptive optics studies to investigate their morphologies and spectral properties on sub-kpc scales. We provide coordinates and JHK photometry for all catalogued objects.     

Dark matter in dwarf spheroidals – II. Observations and modelling of Draco

We present stellar radial velocity data for the Draco dwarf spheroidal (dSph) galaxy obtained using the AF2/WYFFOS instrument combination on the William Herschel Telescope. Our data set consists of 186 member stars, 159 of which have good quality velocities, extending to a magnitude   V ≈19.5  with a mean velocity precision of ≈2 km s⁻¹. As this survey is based on a high-precision photometric target list, it contains many more Draco members at large radii. For the first time, this allows a robust determination of the radial behaviour of the velocity dispersion in a dSph.
We find statistically strong evidence of a rising velocity dispersion consistent with a dark matter halo that has a gently rising rotation curve. There is a <2 σ signature of rotation about the long axis, inconsistent with tidal disruption as the source of the rising dispersion. By comparing our data set with earlier velocities, we find that Draco probably has a binary distribution and fraction comparable to those in the solar neighbourhood.
We apply a novel maximum likelihood algorithm and fit the velocity data to a two parameter spherical model with an adjustable dark matter content and velocity anisotropy. Draco is best fit by a weakly tangentially anisotropic distribution of stellar orbits in a dark matter halo with a very slowly rising rotation law  ( v _circ∝ r ^0.17)  . We are able to rule out both a mass-follows-light distribution and an extended halo with a harmonic core at the 2.5 to 3 σ significance level, depending on the details of our assumptions about Draco's stellar binary population. Our modelling lends support to the idea that the dark matter in dwarf spheroidals is distributed in the form of massive, nearly isothermal haloes.     

Two-dimensional spectroscopy of double-barred galaxies

We describe the results of our spectroscopy for a sample of barred galaxies whose inner regions exhibit an isophotal twist commonly interpreted as a secondary bar. The line-of-sight velocity fields of the ionized gas and stars and the light-of-sight velocity dispersion fields of the stars were constructed from two-dimensional spectroscopy with the 6-m Special Astrophysical Observatory telescope. We detected various types of noncircular motions of ionized gas: radial flows within large-scale bars, counterrotation of the gas and stars at the center of NGC 3945, a polar gaseous disk in NGC 5850, etc. Our analysis of the optical and near-infrared galaxy images (both ground-based and those from the Hubble Space Telescope) revealed circumnuclear minispirals in five objects. The presence of an inner (secondary) bar in the galaxy images is shown to have no effect on the circumnuclear kinematics of the gas and stars. Thus, contrary to popular belief, the secondary bar is not a dynamically decoupled galactic structure. We conclude that the so-called double-barred galaxies are not a separate type of galaxy but are a combination of objects with distinctly different morphologies of their circumnuclear regions.     

Planetary nebula velocities in the disc and bulge of M31

We present radial velocities for a sample of 723 planetary nebulae in the disc and bulge of M31, measured using the WYFFOS fibre spectrograph on the William Herschel Telescope. Velocities are determined using the [O  iii ] λ5007 emission line. Rotation and velocity dispersion are measured to a radius of 50 arcmin (11.5 kpc), the first stellar rotation curve and velocity dispersion profile for M31 to such a radius. Our kinematics are consistent with rotational support at radii well beyond the bulge effective radius of 1.4 kpc, although our data beyond a radius of 5 kpc are limited. We present tentative evidence for kinematic substructure in the bulge of M31 to be studied fully in a later work. This paper is part of an ongoing project to constrain the total mass, mass distribution and velocity anisotropy of the disc, bulge and halo of M31.     

The tidal stripping of satellites

We present an improved analytic calculation for the tidal radius of satellites and test our results against N -body simulations.
The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star within the satellite . We demonstrate that this last point is critical and suggest using three tidal radii to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power-law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically.
Over short times (≲1–2 Gyr ∼1 satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.     

Mass distribution in nearby Abell clusters

We study the mass distribution in six nearby  ( z < 0.06)  relaxed Abell clusters of galaxies A0262, A0496, A1060, A2199, A3158 and A3558. Given the dominance of dark matter in galaxy clusters, we approximate their total density distribution by the Navarro, Frenk & White (NFW) formula characterized by virial mass and concentration. We also assume that the anisotropy of galactic orbits is reasonably well described by a constant and that galaxy distribution traces that of the total density. Using the velocity and position data for 120–420 galaxies per cluster we calculate, after removal of interlopers, the profiles of the lowest order even velocity moments, dispersion and kurtosis. We then reproduce the velocity moments by jointly fitting the moments to the solutions of the Jeans equations. Including the kurtosis in the analysis allows us to break the degeneracy between the mass distribution and anisotropy and constrain the anisotropy as well as the virial mass and concentration. The method is tested in detail on mock data extracted from the N -body simulations of dark matter haloes. We find that the best-fitting Galactic orbits are remarkably close to isotropic in most clusters. Using the fitted pairs of mass and concentration parameters for the six clusters, we conclude that the trend of decreasing concentration for higher masses found in the cosmological N -body simulations is consistent with the data. By scaling the individual cluster data by mass, we combine them to create a composite cluster with 1465 galaxies and perform a similar analysis on such sample. The estimated concentration parameter then lies in the range  1.5 < c < 14  and the anisotropy parameter in the range  −1.1 < β < 0.5  at the 95 per cent confidence level.     

Structure and dynamics of galaxies with a low surface-brightness disc – I. The stellar and ionized-gas kinematics

Photometry and long-slit spectroscopy are presented for a sample of six galaxies with a low surface-brightness stellar disc and a bulge. The characterizing parameters of the bulge and disc components were derived by means of a two-dimensional photometric decomposition of the images of the sample galaxies. Their surface-brightness distribution was assumed to be the sum of the contribution of a Sérsic bulge and an exponential disc, with each component being described by elliptical and concentric isophotes of constant ellipticity and position angle. The stellar and ionized-gas kinematics were measured along the major and minor axes in half of the sample galaxies, whereas the other half was observed only along two diagonal axes. Spectra along two diagonal axes were obtained also for one of the objects with major and minor axis spectra. The kinematic measurements extend in the disc region out to a surface-brightness level  μ _R ≈ 24  mag arcsec⁻², reaching in all cases the flat part of the rotation curve. The stellar kinematics turns out to be more regular and symmetric than the ionized-gas kinematics, which often shows the presence of non-circular, off-plane and non-ordered motions. This raises the question about the reliability of the use of the ionized gas as the tracer of the circular velocity in the modelling of the mass distribution, in particular in the central regions of low surface-brightness galaxies.