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.     

The peculiar rotation curve of NGC 157

We present the results of a new H i , optical, and Hα interferometric study of the nearby spiral galaxy NGC 157. Our combined C- and D-array observations with the VLA show a large-scale, ring-like structure in the neutral hydrogen underlying the optical disc, together with an extended, low surface density component going out to nearly twice the Holmberg radius. Beginning just inside the edge of the star-forming disc, the line of nodes in the gas disc commences a 60° warp, while at the same time, the rotation velocity drops by almost half its peak value of 200 km s⁻¹, before levelling off again in the outer parts. While a flat rotation curve in NGC 157 cannot be ruled out, supportive evidence for an abrupt decline comes from the ionized gas kinematics, the optical surface photometry, and the global H i profile. A standard 'maximum-disc' mass model predicts comparable amounts of dark and luminous matter within NGC 157. Alternatively, a model employing a disc truncated at 2 disc scalelengths could equally well account for the unusual form of the rotation curve in NGC 157.     

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.     

Supersonic motions of galaxies in clusters

We study motions of galaxies in galaxy clusters formed in the concordance Λ cold dark matter cosmology. We use high-resolution cosmological simulations that follow the dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing the motions of galaxies and the properties of intracluster gas in a sample of eight simulated clusters at z = 0, we study the velocity dispersion profiles of the dark matter, gas and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ≈1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move slightly faster than the dark matter particles. The magnitude of the velocity bias,   b _v ≈ 1.1  , is, however, smaller than the bias estimated for subhaloes in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ≈20–30 per cent of the gas sound speed. These random motions provide about 10 per cent of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.     

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.     

Orientation dependency of broad-line widths in quasars and consequences for black hole mass estimation

We have identified two new galaxies with gas counter-rotation (NGC 1596 and 3203) and have confirmed similar behaviour in another one (NGC 128), this using results from separate studies of the ionized-gas and stellar kinematics of a well-defined sample of 30 edge-on disc galaxies. Gas counter-rotators thus represent 10 ± 5 per cent of our sample, but the fraction climbs to 21 ± 11 per cent when only lenticular (S0) galaxies are considered and to 27 ± 13 per cent for S0 galaxies with detected ionized gas only. Those fractions are consistent with but slightly higher than previous studies. A compilation from well-defined studies of S0 galaxies in the literature yields fractions of 15 ± 4 and 23 ± 5 per cent, respectively. Although mainly based on circumstantial evidence, we argue that the counter-rotating gas originates primarily from minor mergers and tidally induced transfer of material from nearby objects. Assuming isotropic accretion, twice those fractions of objects must have undergone similar processes, underlining the importance of (minor) accretion for galaxy evolution. Applications of gas counter-rotators to barred galaxy dynamics are also discussed.     

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.     

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.     

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.     

Dark matter in early-type galaxies: dynamical modelling of IC 1459, IC 3370, NGC 3379 and NGC 4105

We analyse long-slit spectra of four early-type galaxies which extend from ∼1 to 3 effective radii: IC 1459; IC 3370; NGC 3379 and NGC 4105. We have extracted the full line-of-sight velocity distribution (in the case of NGC 3379 we also used data from the literature), which we model using the two-integral approach. Using two-integral modelling, we find no strong evidence for dark haloes, but the fits suggest that three-integral modelling is necessary. We also find that the inferred constant mass-to-light ratio in all the four cases is typical for early-type galaxies. Finally, we also discuss the constraints on the mass-to-light ratio, which can be obtained using X-ray haloes in the case of IC 1459, NGC 3379 and NGC 4105, and compare the estimated values with the predictions from the dynamical modelling.     

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.     

Using slitless spectroscopy to study the kinematics of the planetary nebula population in M94

The planetary nebula populations of relatively nearby galaxies can be easily observed and provide both a distance estimate and a tool with which dynamical information can be obtained. Usually the requisite radial velocities are obtained by multi-object spectroscopy once the planetary nebulae have been located by direct imaging. Here we report on a technique for measuring planetary nebula kinematics using the double-beam ISIS spectrograph at the William Herschel Telescope in a novel slitless mode, which enables the detection and radial velocity measurements to be combined into a single step. The results on our first target, the Sab galaxy NGC 4736, allow the velocity dispersion of the stellar population in a disc galaxy to be traced out to four scalelengths for the first time and are consistent with a simple isothermal sheet model.     

Tidal evolution of discy dwarf galaxies in the Milky Way potential: the formation of dwarf spheroidals

We conduct high-resolution collisionless N -body simulations to investigate the tidal evolution of dwarf galaxies on an eccentric orbit in the Milky Way (MW) potential. The dwarfs originally consist of a low surface brightness stellar disc embedded in a cosmologically motivated dark matter halo. During 10 Gyr of dynamical evolution and after five pericentre passages, the dwarfs suffer substantial mass loss and their stellar component undergoes a major morphological transformation from a disc to a bar and finally to a spheroid. The bar is preserved for most of the time as the angular momentum is transferred outside the galaxy. A dwarf spheroidal (dSph) galaxy is formed via gradual shortening of the bar. This work thus provides a comprehensive quantitative explanation of a potentially crucial morphological transformation mechanism for dwarf galaxies that operates in groups as well as in clusters. We compare three cases with different initial inclinations of the disc and find that the evolution is fastest when the disc is coplanar with the orbit. Despite the strong tidal perturbations and mass loss, the dwarfs remain dark matter dominated. For most of the time, the one-dimensional stellar velocity dispersion, σ, follows the maximum circular velocity, V _max, and they are both good tracers of the bound mass. Specifically, we find that   M _bound∝ V ^3.5_max  and     in agreement with earlier studies based on pure dark matter simulations. The latter relation is based on directly measuring the stellar kinematics of the simulated dwarf, and may thus be reliably used to map the observed stellar velocity dispersions of dSphs to halo circular velocities when addressing the missing satellites problem.     

The association between gas and galaxies – I. CFHT spectroscopy and pair analysis

We investigate the relative distribution of the gaseous contents of the Universe (as traced by a sample of Lyα absorbers), and the luminous baryonic matter (as traced by a redshift survey of galaxies in the same volume searched for Lyα absorbers), along 16 lines of sight (LOS) between redshifts 0 and 1. Our galaxy redshift survey was made with the multi-object spectrograph on the Canada–France–Hawaii Telescope and, when combined with galaxies from the literature in the same LOS, gives us a galaxy sample of 636 objects. By combining this with an absorption-line sample of 406 absorbing systems drawn from published works, we are able to study the relationship between gas and galaxies over the latter half of the age of the Universe. A correlation between absorbers and galaxies is detected out to separation of 1.5 Mpc. This correlation is weaker than the galaxy–galaxy correlation. There is also some evidence that the absorbing systems seen in C  iv are more closely related to galaxies, although this correlation could be with column density rather than metallicity. The above results are all consistent with the absorbing gas and the galaxies coexisting in dark matter filaments and knots as predicted by current models where the column density of the absorbing gas is correlated with the underlying matter density.     

Collisional baryonic dark matter haloes

If dark haloes are composed of dense gas clouds, as has recently been inferred, then collisions between clouds lead to galaxy evolution. Collisions introduce a core in an initially singular dark matter distribution, and can thus help to reconcile scale-free initial conditions – such as are found in simulations – with observed haloes, which have cores. A pseudo-Tully–Fisher relation, between halo circular speed and visible mass (not luminosity), emerges naturally from the model: M _vis∝ V ^7/2.
Published data conform astonishingly well to this theoretical prediction. For our sample of galaxies, the mass–velocity relationship has much less scatter than the Tully–Fisher relation, and holds as well for dwarf galaxies (where diffuse gas makes a sizeable contribution to the total visible mass) as it does for giants. It seems very likely that this visible-mass/velocity relationship is the underlying physical basis for the Tully–Fisher relation, and this discovery in turn suggests that the dark matter is both baryonic and collisional.     

Long‐slit échelle spectroscopy of galactic outflows: The case ofNGC 4449

At high redshift the ubiquity of outflows and winds in strongly star‐forming galaxies has been demonstrated using rest frame UV absorption lines. In the cases with optical emission lines, the studies mostly had to rely on low and intermediate dispersion spectra. This implies that for detailed studies of galactic wind physics we have to use local objects. In particular, dwarf galaxies are well suited to extrapolation to high redshift protogalaxies. Several kinematic studies of strongly starforming dwarf galaxies using Fabry‐Pérot and IFU spectrographs exist. Unfortunately, similar as for high redshift galaxies the employed spectral resolution is often significantly higher that the thermal line width. As a result faint high velocity features and details of the turbulent motion are hidden or unresolved. Here we will present an analysis of the ionized gas kinematics of the prototypical star‐forming irregular galaxy NGC 4449 using long‐slit, high‐dispersion échelle spectra. (© 2014 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.     

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.     

A minor-merger model for NGC 7479

The overall morphology of the barred spiral galaxy NGC 7479 is modelled in numerical simulations of a minor merger. Special attention is paid to the morphology and velocity field of the asymmetric spiral structure and the strong stellar bar. The mass of the satellite galaxy is 1/10 of the mass of the primary disc, or 1/30 of the total mass of the primary. The satellite is placed initially in a circular prograde orbit at six disc scalelengths from the centre of the primary. We follow the evolution of the merger until the secondary galaxy reaches the nuclear region of the primary. A comparison between the modelled and observed morphologies of the stellar and the ionized and neutral gas distributions and velocity fields supports the hypothesis that the transient look of NGC 7479 is a result of a minor merger. We vary several of the initial parameters of the merger and discuss their effects on the resulting morphology. The merging satellite galaxy is likely to lie within the bar of NGC 7479. We identify a possible candidate in the observational data. We discuss briefly the most probable future evolution of NGC 7479 in the light of our minor-merger simulations, and conclude that NGC 7479 is likely to evolve toward an earlier Hubble type.