The mass and anisotropy profiles of galaxy clusters from the projected phase-space density: testing the method on simulated data

We present a new method of constraining the mass and velocity anisotropy profiles of galaxy clusters from kinematic data. The method is based on a model of the phase-space density, which allows the anisotropy to vary with radius between two asymptotic values. The characteristic scale of transition between these asymptotes is fixed and tuned to a typical anisotropy profile resulting from cosmological simulations. The model is parametrized by two values of anisotropy, at the centre of the cluster and at infinity, and two parameters of the NFW density profile, the scale radius and the scale mass. In order to test the performance of the method in reconstructing the true cluster parameters, we analyse mock kinematic data for 20 relaxed galaxy clusters generated from a cosmological simulation of the standard Λ cold dark matter model. We use Bayesian methods of inference and the analysis is carried out following the Markov Chain Monte Carlo approach. The parameters of the mass profile are reproduced quite well, but we note that the mass is typically underestimated by 15 per cent, probably due to the presence of small velocity substructures. The constraints on the anisotropy profile for a single cluster are in general barely conclusive. Although the central asymptotic value is determined accurately, the outer one is subject to significant systematic errors caused by substructures at large clustercentric distance. The anisotropy profile is much better constrained if one performs joint analysis of at least a few clusters. In this case, it is possible to reproduce the radial variation of the anisotropy over two decades in radius inside the virial sphere.     

Resonances and bifurcations in axisymmetric scale-free potentials

We investigate an analytical treatment of bifurcations of families of resonant 'thin' tubes in axisymmetric galactic potentials. We verify that the most relevant bifurcations are due to the (1:1) resonance producing the 'inclined' orbits through two different mechanisms: from the disc orbit and from the 'thin' tube associated with the vertical oscillation. The closest resonances occurring after these are the (4:3) resonance in the oblate case and the (2:1) resonance in the prolate case. The (1:1) resonances are treated in a straightforward way using a second-order truncated normal form. The higher order resonances are instead cumbersome to investigate, because the normal form has to be truncated to a high degree and the number of terms grows very rapidly. We therefore adopt a further simplification giving analytical formulae for the values of the parameters at which bifurcations ensue and compare them with selected numerical results. Thanks to the asymptotic nature of the series involved, the predictions are reliable well beyond the convergence radius of the original series.     

Star Formation Properties of the Galaxy Cluster Abell 2199

Through the morphological classifications for 290 member galaxies in the nearby galaxy Abell 2199, the star formation rates and their relations with their morphology and related physical properties are investigated in this paper. It is found that the typical star formation rate in galaxies of this galaxy cluster is strongly correlated with the H_α equivalent width, and the degree of discontinuity of the galaxy spectrum at 4000 Å is also strongly correlated with the stellar mass included in the galaxy. It is also found that star formation activities in these galaxies do not exhibit the obvious circumstance effect. This result indicates that this galaxy cluster is still situated at the stage of the violent dynamical evolution and far from the dynamical equilibrium.     

Complexity in small‐scale dwarf spheroidal galaxies – Ludwig Biermann Award Lecture 2008

Our knowledge about the dynamics, the chemical abundances and the evolutionary histories of the more luminous dwarf spheroidal (dSph) galaxies is constantly growing. However, very little is known about the enrichment of the ultra‐faint systems recently discovered in large numbers in large sky surveys. Current low‐resolution spectroscopy and photometric data indicate that these galaxies are highly dark matter dominated and predominantly metal poor. On the other hand, recent high‐resolution abundance analyses indicate that some dwarf galaxies experienced highly inhomogeneous chemical enrichment, where star formation proceeds locally on small scales. In this article, I will review the kinematic and chemical abundance information of the Milky Way satellite dSphs that is presently available from low‐ and high resolution spectroscopy. Moreover, some of the most peculiar element and inhomogeneous enrichment patterns will be discussed and related to the question of to what extent the faintest dSph candidates could have contributed to the Galactic halo, compared to more luminous systems (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Galactic satellite systems: radial distribution and environment dependence of galaxy morphology

We have studied the radial distribution of the early (E/S0) and late (S/Irr) types of satellites around bright host galaxies. We made a volume-limited sample of 4986 satellites brighter than   M _r=−18.0  associated with 2254 hosts brighter than   M _r=−19.0  from the Sloan Digital Sky Survey Data Release 5 sample. The morphology of satellites is determined by an automated morphology classifier, but the host galaxies are visually classified. We found segregation of satellite morphology as a function of the projected distance from the host galaxy. The amplitude and shape of the early-type satellite fraction profile are found to depend on the host luminosity. This is the morphology–radius/density relation at the galactic scale. There is a strong tendency for morphology conformity between the host galaxy and its satellites. The early-type fraction of satellites hosted by early-type galaxies is systematically larger than that of late-type hosts, and is a strong function of the distance from the host galaxies. Fainter satellites are more vulnerable to the morphology transformation effects of hosts. Dependence of satellite morphology on the large-scale background density was detected. The fraction of early-type satellites increases in high-density regions for both early- and late-type hosts. It is argued that the conformity in morphology of galactic satellite system is mainly originated by the hydrodynamical and radiative effects of hosts on satellites.     

Masses for the Local Group and the Milky Way

We use the very large Millennium Simulation of the concordance Λ cold dark matter cosmogony to calibrate the bias and error distribution of Timing Argument estimators of the masses of the Local Group and of the Milky Way. From a large number of isolated spiral–spiral pairs similar to the Milky Way/Andromeda system, we find the interquartile range of the ratio of timing mass to true mass to be a factor of 1.8, while the 5 and 95 per cent points of the distribution of this ratio are separated by a factor of 5.7. Here, we define true mass as the sum of the 'virial' masses, M ₂₀₀, of the two dominant galaxies. For present best values of the distance and approach velocity of Andromeda, this leads to a median likelihood estimate of the true mass of the Local Group of  5.27 × 10¹² M_⊙  or  log  M _LG/M_⊙= 12.72  , with an interquartile range of [12.58, 12.83] and a 5–95 per cent range of [12.26, 13.01]. Thus, a 95 per cent lower confidence limit on the true mass of the Local Group is  1.81 × 10¹² M_⊙  . A timing estimate of the Milky Way's mass based on the large recession velocity observed for the distant satellite Leo I works equally well, although with larger systematic uncertainties. It gives an estimated virial mass for the Milky Way of  2.43 × 10¹² M_⊙  with a 95 per cent lower confidence limit of  0.80 × 10¹² M_⊙  .     

Modelling the dynamical evolution of the Bootes dwarf spheroidal galaxy

We investigate a wide range of possible evolutionary histories for the recently discovered Bootes dwarf spheroidal galaxy, a Milky Way satellite. By means of N -body simulations, we follow the evolution of possible progenitor galaxies of Bootes for a variety of orbits in the gravitational potential of the Milky Way. The progenitors considered cover the range from dark matter-free star clusters to massive, dark matter-dominated outcomes of cosmological simulations. For each type of progenitor and orbit, we compare the observable properties of the remnant after 10 Gyr with those of Bootes observed today. Our study suggests that the progenitor of Bootes must have been, and remains now, dark matter-dominated. In general, our models are unable to reproduce the observed high velocity dispersion in Bootes without dark matter. Our models do not support time-dependent tidal effects as a mechanism able to inflate significantly the internal velocity dispersion. As none of our initially spherical models is able to reproduce the elongation of Bootes, our results suggest that the progenitor of Bootes may have had some intrinsic flattening. Although the focus of this paper is the Bootes dwarf spheroidal, these models may be of general relevance to understanding the structure, stability and dark matter content of all dwarf spheroidal galaxies.     

Galactic H i on the 50-au scale in the direction of three extragalactic sources observed with MERLIN

We present MERLIN observations of Galactic 21-cm H  i absorption at an angular resolution of  ∼0.1–0.2  arcsec and a velocity resolution of 0.5 km s⁻¹, in the direction of three moderately low latitude  (−8° < b < −12°)  extragalactic radio sources, 3C 111, 3C 123 and 3C 161, all of which are heavily reddened. H  i absorption is observed against resolved background emission sources up to ∼2 arcsec in extent and we distinguish details of the opacity distribution within 1–1.5 arcsec regions towards 3C 123 and 3C 161. This study is the second MERLIN investigation of small-scale structure in interstellar H  i (earlier work probed Galactic H  i in the directions of the compact sources 3C 138 and 3C 147). The 0.1-arcsec scale is intermediate between H  i absorption studies made with other fixed element interferometers with resolution of 1–10 arcsec and very long baseline interferometry studies with resolutions of 10–20 mas. At a scale of 1 arcsec (about 500 au), prominent changes in Galactic H  i opacity in excess of 1–1.5 are determined in the direction of 3C 161 with a signal-to-noise ratio of at least 10σ. Possible fluctuations in the H  i opacity at the level of about 1 are detected at the  2.5–3σ  level in the direction of 3C 123.     

The distribution function of dark matter in massive haloes

We study the distribution function (DF) of dark matter particles in haloes of mass range  10¹⁴–10¹⁵ M_⊙  . In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard Λ cold dark matter (ΛCDM) model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L , a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form     . This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (β₀ and  β_∞  ) and the scale of transition between them ( L ₀). The energy part   f _E ( E )  is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase-space density   Q =ρ( r )/σ³( r )  .