Anisotropic and inhomogeneous S-type Riemann ellipsoids inside spheroidal halos. II

Classical S-type Riemann ellipsoids are generalized taking a Ferrers-type inhomogeneity in the mass distribution, a three-dimensional anisotropy in the dispersion in the velocities, and the gravitation of a spherical halo into account. A Ferrers inhomogeneity in the mass distribution does not affect the equilibrium conditions for the ellipsoids, but only changes the numerical coefficients in the equilibrium parameters. An anisotropy in the dispersion in the velocities does change the equilibrium conditions, expanding or contracting the region in which embedded ellipsoids exist. A halo extends this region in all cases. __________ Translated from Astrofizika, Vol. 49, No. 3, pp.359–373 (August 2006).     

Triaxial haloes and particle dark matter detection

This paper presents the properties of a family of scale-free triaxial haloes. We adduce arguments to suggest that the velocity ellipsoids of such models are aligned in conical coordinates. We provide an algorithm to find the set of conically aligned velocity second moments that support a given density against the gravity field of the halo. The case of the logarithmic ellipsoidal model – the simplest triaxial generalization of the familiar isothermal sphere – is examined in detail. The velocity dispersions required to hold up the self-consistent model are analytic. The velocity distribution of the dark matter can be approximated as a triaxial Gaussian with semiaxes equal to the velocity dispersions.
There are roughly 20 experiments worldwide that are searching for evidence of scarce interactions between weakly interacting massive-particle dark matter (WIMP) and detector nuclei. The annual modulation signal, caused by the Earth's rotation around the Sun, is a crucial discriminant between WIMP events and the background. The greatest rate is in June, the least in December. We compute the differential detection rate for energy deposited by the rare WIMP–nucleus interactions in our logarithmic ellipsoidal halo models. Triaxiality and velocity anisotropy change the total rate by up to ∼40 per cent, and have a substantial effect on the amplitude of the annual modulation signal. The overall rate is greatest, but the amplitude of the modulation is weakest, in our radially anisotropic halo models. Even the sign of the signal can be changed. Restricting attention to low energy events, the models predict that the maximum rate occurs in December, and not in June.     

The haloes of merger remnants

We perform collisionless N -body simulations of 1:1 galaxy mergers, using models which include a galaxy halo, disc and bulge, focusing on the behaviour of the halo component. The galaxy models are constructed without recourse to a Maxwellian approximation. We investigate the effect of varying the galaxies' orientation, their mutual orbit and the initial velocity anisotropy or cusp strength of the haloes upon the remnant halo density profiles and shape, as well as on the kinematics. We observe that the halo density profile (determined as a spherical average, an approximation we find appropriate) is exceptionally robust in mergers, and that the velocity anisotropy of our remnant haloes is nearly independent of the orbits or initial anisotropy of the haloes. The remnants follow the halo anisotropy – local density slope (β–γ) relation suggested by Hansen & Moore in the inner parts of the halo, but β is systematically lower than this relation predicts in the outer parts. Remnant halo axis ratios are strongly dependent on the initial parameters of the haloes and on their orbits. We also find that the remnant haloes are significantly less spherical than those described in studies of simulations which include gas cooling.     

Anisotropy and rotation in homeoidally striated Jacobi ellipsoids

In this paper a unified theory of systematically rotating and peculiar motions is developed for homeoidally striated Jacobi ellipsoids, where both real and imaginary rotations are considered. The effect of positive or negative residual motion excess along the equatorial plane is considered to be equivalent either to an additional real or an imaginary rotation, respectively. The principle results consist of (i) the discovery that homeoidally striated Jacobi ellipsoids always admit an adjoint configuration i.e. a classical Jacobi ellipsoid of equal mass and axes; (ii) the establishment of further constraints on the amount of residual velocity anisotropy along the principal axes for triaxial configurations; (iii) the finding that bifurcation points from axisymmetric to triaxial configurations occur as in classical Jacobi ellipsoids, contrary to earlier findings. An interpretation of recent results from numerical simulations on stability is provided in the light of the model. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Satellite kinematics – I. A new method to constrain the halo mass–luminosity relation of central galaxies

Satellite kinematics can be used to probe the masses of dark matter haloes of central galaxies. In order to measure the kinematics with sufficient signal-to-noise ratio, one uses the satellite galaxies of a large number of central galaxies stacked according to similar properties (e.g. luminosity). However, in general, the relation between the luminosity of a central galaxy and the mass of its host halo will have non-zero scatter. Consequently, this stacking results in combining the kinematics of satellite galaxies in haloes of different masses, which complicates the interpretation of the data. In this paper, we present an analytical framework to model satellite kinematics, properly accounting for this scatter and for various selection effects. We show that in the presence of scatter in the halo mass–luminosity relation, the commonly used velocity dispersion of satellite galaxies can not be used to infer a unique halo mass–luminosity relation. In particular, we demonstrate that there is a degeneracy between the mean and the scatter of the halo mass–luminosity relation. We present a new technique that can break this degeneracy, and which involves measuring the velocity dispersions using two different weighting schemes: host weighting (each central galaxy gets the same weight) and satellite weighting (each central galaxy gets a weight proportional to its number of satellites). The ratio between the velocity dispersions obtained using these two weighting schemes is sensitive to the scatter in the halo mass–luminosity relation, and can thus be used to infer a unique relation between light and mass from the kinematics of satellite galaxies.     

Riemann Ellipsoids with Spherical Halos

Possible ellipsoidal figures of equilibrium are obtained for a rotating, gravitating fluid mass with internal mass flows of constant vorticity, embedded inside a homogeneous gravitating sphere. The classical ellipsoidal figures of equilibrium are generalized and new S-ellipsoids and ellipsoids with oblique rotation are obtained. The stability of embedded S-ellipsoids is investigated and the criterion for their stability is obtained. The existence of an ellipsoid with oblique rotation of type II inside a relatively dense halo becomes impossible.     

Homeoidally striated density profiles: sequences of virial equilibrium configurations with constant anisotropy parameters

The formulation of the tensor virial equations is generalized to unrelaxed configurations, where virial equilibrium does not coincide with dynamical (or hydrostatic) equilibrium. Homeoidally striated, Jacobi ellipsoids, which generalize classical Jacobi ellipsoids, are studied in detail. Further investigation is devoted to the generation of sequences of virial equilibrium configurations where the anisotropy parameters are left unchanged, including both flattened and elongated, triaxial configurations, and the determination of the related bifurcation points. An application is made to dark matter haloes hosting giant galaxies (M ≈ 10¹² m_⊙), with regard to assigned initial and final configuration, following and generalizing to many respects a procedure conceived by Thuan & Gott (1975). The dependence of the limiting axis ratios, below which no configuration is allowed for the sequence under consideration, on the change in mass, total energy, and angular momentum, during the evolution, is illustrated in some representative situations. The dependence of the axis ratios, ε₃₁ and ε₂₁, on a parameter, related to the initial conditions of the density perturbation, is analysed in connection with a few special cases. The same is done for the rotation parameters. Within the range of the rotation parameter, λ, deduced from high‐resolution numerical simulations, the shape of dark matter haloes is mainly decided by the amount of anisotropy in residual velocity distribution. On the other hand, the contribution of rotation has only a minor effect on the meridional plane, and no effect on the equatorial plane, as bifurcation points occur for larger values of λ. To this respect, dark matter haloes are found to resemble giant elliptical galaxies. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The origin and formation of cuspy density profiles through violent relaxation of stellar systems

It is shown that the cuspy density distributions observed in the cores of elliptical galaxies can be realized by dissipationless gravitational collapse. The initial models consist of power-law density spheres such as ρ ∝ r ⁻¹ with anisotropic velocity dispersions. Collapse simulations are carried out by integrating the collisionless Boltzmann equation directly, on the assumption of spherical symmetry. From the results obtained, the extent of constant density cores, formed through violent relaxation, decreases as the velocity anisotropy increases radially, and practically disappears for extremely radially anisotropic models. As a result, the relaxed density distributions become more cuspy with increasing radial velocity anisotropy. It is thus concluded that the velocity anisotropy could be a key ingredient for the formation of density cusps in a dissipationless collapse picture. The velocity dispersions increase with radius in the cores according to the nearly power-law density distributions. The power-law index, n , of the density profiles, defined as ρ ∝ r ^{− n} , changes from n ≈2.1 at intermediate radii to a shallower power than n ≈2.1 toward the centre. This density bend can be explained from our postulated local phase-space constraint that the phase-space density accessible to the relaxed state is determined at each radius by the maximum phase-space density of the initial state.     

Three-integral models for axisymmetric galactic discs

We present new equilibrium component distribution functions that depend on three analytic integrals in a Stäckel potential, and that can be used to model stellar discs of galaxies. These components are generalizations of two-integral ones and can thus provide thin discs in the two-integral approximation. Their most important properties are the partly analytical expression for their moments, the disc-like features of their configuration space densities (exponential decline in the galactic plane and finite extent in the vertical direction) and the anisotropy of their velocity dispersions. We further show that a linear combination of such components can fit a van der Kruit disc.     

Kinematics and metallicity analysis for nearby F,G and K stars

A sample containing 1 026 stars of spectral types F, G, and K, mainly dwarfs, from the solar neighbourhood with available space velocities and metallicities is treated. The treatment comprises a statistical analysis of the metallicity and velocity data and calculation of galactocentric orbits. Sample stars identified as members of the galactic halo are detached from the rest of the sample based on the values of their metallicities, velocity components and galactocentric orbits. In identifying halo stars a new, kinematical, criterion is proposed. Except one, these halo stars are the metal‐poorest ones in the sample. Besides, they have very high velocities with respect to LSR. On the other hand, the separation between the thin disc and thick one is done statistically based on LSR space velocities, membership probability (Schwarzschild distribution with assumed parameters) and galactocentric orbits. In the metallicity these two groups are not much different. For each of the three subsamples the mean motion and velocity ellipsoid are calculated. The elements of the velocity ellipsoids agree well with the values found in the literature, especially for the thin disc. The fractions of the subsystems found for the present sample are: thin disc 93%, thick disc 6%, halo 1%. The sample stars established to be members of the thin disc are examined for existence of star streams. Traces of both, known and unknown, star streams are not found (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Exploring disc galaxy dynamics using integral field unit data

In order to test the basic equations believed to dictate the dynamics of disc galaxies, we present and analyse deep two-dimensional spectral data obtained using the PPAK integral field unit for the early-type spiral systems NGC 2273, NGC 2985, NGC 3898 and NGC 5533. We describe the care needed to obtain and process such data to a point where reliable kinematic measurements can be obtained from these observations, and a new more optimal method for deriving the rotational motion and velocity dispersions in such disc systems. The data from NGC 2273 and NGC 2985 show systematic variations in velocity dispersion with azimuth, as one would expect if the shapes of their velocity ellipsoids are significantly anisotropic, while the hotter discs in NGC 3898 and NGC 5533 appear to have fairly isotropic velocity dispersions. Correcting the rotational motion for asymmetric drift using the derived velocity dispersions reproduces the rotation curves inferred from emission lines reasonably well, implying that this correction is quite robust, and that the use of the asymmetric drift equation is valid. NGC 2985 is sufficiently close to face on for the data, combined with the asymmetric drift equation, to determine all three components of the velocity ellipsoid. The principal axes of this velocity ellipsoid are found to be in the ratio  σ _z : σ_φ : σ _R ≈ 0.7 : 0.7 : 1  , which shows unequivocally that this disc distribution function respects a third integral of motion. The ratio is also consistent with the predictions of epicyclic theory, giving some confidence in the application of this approximation to even fairly early-type disc galaxies.     

Effects of merging histories on angular momentum distribution of dark matter haloes

The effects of merging histories of proto-objects on the angular momentum distributions of the present-time dark matter haloes are analysed. An analytical approach to the analysis of the angular momentum distributions assumes that the haloes are initially homogeneous ellipsoids and that the growth of the angular momentum of the haloes halts at their maximum expansion time. However, the maximum expansion time cannot be determined uniquely, because in the hierarchical clustering scenario each progenitor, or subunit, of the halo has its own maximum expansion time. Therefore the merging history of the halo may be important in estimating its angular momentum. Using the merger tree model by Rodrigues &38; Thomas, which takes into account the spatial correlations of the density fluctuations, we have investigated the effects of the merging histories on the angular momentum distributions of dark matter haloes. It was found that the merger effects, that is, the effects of the inhomogeneity of the maximum expansion times of the progenitors which finally merge together into a halo, do not strongly affect the final angular momentum distributions, so that the homogeneous ellipsoid approximation happens to be good for the estimation of the angular momentum distribution of dark matter haloes. This is because the effect of the different directions of the angular momenta of the progenitors cancels out the effect of the inhomogeneity of the maximum expansion times of the progenitors.   The contribution of the orbital angular momentum to the total angular momentum when two or more pre-existing haloes merge together was also investigated. It is shown that this contribution is more important than that of the angular momentum of diffuse accreting matter to the total angular momentum, especially when the mergers occur many times.     

Stability of a Riemann S Ellipsoid with a Spheroidal Halo

The stability of Riemann S ellipsoids inside an oblate halo with respect to the second form of oscillations is investigated. It is shown that some ellipsoids with reverse internal circulation of matter, which are stable inside a spherical halo or in its absence, become unstable with respect to second odd forms of oscillation inside an oblate halo. Here there is asymmetry between conjugate ellipsoids from the standpoint of their stability. Only those conjugate ellipsoids that correspond to higher frequencies of reverse circulation of matter than their corresponding angular rotation rates are unstable. The domains of instability of light and heavy conjugate embedded ellipsoids are obtained as a function of the oblateness measure and relative density of the halo.     

Velocity dispersions of dwarf spheroidal galaxies: dark matter versus MOND

We present predictions for the line-of-sight velocity dispersion profiles of dwarf spheroidal galaxies and compare them to observations in the case of the Fornax dwarf. The predictions are made in the framework of standard dynamical theory of spherical systems with different velocity distributions. The stars are assumed to be distributed according to Sérsic laws with parameters fitted to observations. We compare predictions obtained assuming the presence of dark matter haloes (with density profiles adopted from N -body simulations) with those resulting from Modified Newtonian Dynamics (MOND). If the anisotropy of velocity distribution is treated as a free parameter, observational data for Fornax are reproduced equally well by models with dark matter and with MOND. If stellar mass-to-light ratio of 1 M_⊙/L_⊙ is assumed, the required mass of the dark halo is     , two orders of magnitude larger than the mass in stars. The derived MOND acceleration scale is     . In both cases a certain amount of tangential anisotropy in the velocity distribution is needed to reproduce the shape of the velocity dispersion profile in Fornax.     

Testing fundamental physics with distant star clusters: theoretical models for pressure-supported stellar systems

We investigate the mean velocity dispersion and the velocity dispersion profile of stellar systems in modified Newtonian dynamics (MOND), using the N -body code n-mody , which is a particle-mesh-based code with a numerical MOND potential solver developed by Ciotti, Londrillo & Nipoti. We have calculated mean velocity dispersions for stellar systems following Plummer density distributions with masses in the range of 10⁴ to  10⁹ M_⊙  and which are either isolated or immersed in an external field. Our integrations reproduce previous analytic estimates for stellar velocities in systems in the deep MOND regime  ( a _i, a _e≪ a ₀)  , where the motion of stars is either dominated by internal accelerations  ( a _i≫ a _e)  or constant external accelerations  ( a _e≫ a _i)  . In addition, we derive for the first time analytic formulae for the line-of-sight velocity dispersion in the intermediate regime  ( a _i∼ a _e∼ a ₀)  . This allows for a much-improved comparison of MOND with observed velocity dispersions of stellar systems. We finally derive the velocity dispersion of the globular cluster Pal 14 as one of the outer Milky Way halo globular clusters that have recently been proposed as a differentiator between Newtonian and MONDian dynamics.     

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.     

Dark matter subhaloes in numerical simulations

We use cosmological Λ cold dark matter (CDM) numerical simulations to model the evolution of the substructure population in 16 dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10¹¹ to  10¹⁴ M_⊙  at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation  d n /d v = (1/8)( v _cmax/ v _cmax,host)⁻⁴  . Halo to halo variance in the VDF is a factor of roughly 2 to 4. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for   r ≳ 0.3 r _vir  , where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease towards the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias towards more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power-law cosmologies of   P ∝ kⁿ   , and demonstrate that a steeper spectral index, n , results in significantly less substructure.     

The graininess of dark matter haloes

We use the recently completed one billion particle Via Lactea II Λ cold dark matter simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as the structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions, the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo cannot be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.