Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMAMPOSSt: Modelling Anisotropy and Mass Profiles of Observed Spherical Systems. I. Gaussian 3D velocities
97
0
0.0
(
0
)
Authors
Gary A. Mamon n IAP
Ask ChatGPT about the research
No Arabic abstract
Mass modelling of spherical systems through internal motions is hampered by the mass/velocity anisotropy (VA) degeneracy inherent in the Jeans equation, as well as the lack of techniques that are both fast and adaptable to realistic systems. A new fast method, called MAMPOSSt, which performs a maximum likelihood fit of the distribution of observed tracers in projected phase space, is developed and thoroughly tested. MAMPOSSt assumes a shape for the gravitational potential, but instead of postulating a shape for the distribution function in terms of energy and angular momentum, or supposing Gaussian line-of-sight velocity distributions, MAMPOSSt assumes a VA profile and a shape for the 3D velocity distribution, here Gaussian. MAMPOSSt requires no binning, differentiation, nor extrapolation of the observables. Tests on cluster-mass haloes from LambdaCDM cosmological simulations show that, with 500 tracers, MAMPOSSt is able to jointly recover the virial radius, tracer scale radius, dark matter scale radius and outer or constant VA with small bias (<10% on scale radii and <2% on the two other quantities) and inefficiencies of 10%, 27%, 48% and 20%, respectively. MAMPOSSt does not perform better when some parameters are frozen, and even worse when the virial radius is set to its true value, which appears to be the consequence of halo triaxiality. The accuracy of MAMPOSSt depends weakly on the adopted interloper removal scheme, including an efficient iterative Bayesian scheme that we introduce here, which can directly obtain the virial radius with as good precision as MAMPOSSt. Our tests show that MAMPOSSt with Gaussian 3D velocities is very competitive with, and up to 1000x faster than other methods. Hence, MAMPOSSt is a very powerful and rapid tool for the mass and anisotropy modeling of systems such as clusters and groups of galaxies, elliptical and dwarf spheroidal galaxies.
rate research
Read More
Traditionally, the mass / velocity anisotropy degeneracy (MAD) inherent in the spherical, stationary, non-streaming Jeans equation has been handled by assuming a mass profile and fitting models to the observed kinematical data. Here, the opposite approach is considered: the equation of anisotropic kinematic projection is inverted for known arbitrary anisotropy to yield the space radial velocity dispersion profile in terms of an integral involving the radial profiles of anisotropy and isotropic dynamical pressure. Then, through the Jeans equation, the mass profile is derived in terms of double integrals of observable quantities. Single integral formulas for both deprojection and mass inversion are provided for several simple anisotropy models (isotropic, radial, circular, general constant, Osipkov-Merritt, Mamon-Lokas and Diemand-Moore-Stadel). Tests of the mass inversion on NFW models with these anisotropy models yield accurate results in the case of perfect observational data, and typically better than 70% (in 4 cases out of 5) accurate mass profiles for the sampling errors expected from current observational data on clusters of galaxies. For the NFW model with mildly increasing radial anisotropy, the mass is found to be insensitive to the adopted anisotropy profile at 7 scale radii and to the adopted anisotropy radius at 3 scale radii. This anisotropic mass inversion method is a useful complementary tool to analyze the mass and anisotropy profiles of spherical systems. It provides the practical means to lift the MAD in quasi-spherical systems such as globular clusters, round dwarf spheroidal and elliptical galaxies, as well as groups and clusters of galaxies, when the anisotropy of the tracer is expected to be linearly related to the slope of its density.
We investigate the connection between collisionless equilibria and the phase-space relation between density $rho$ and velocity dispersion $sigma$ found in simulations of dark matter halo formation, $F=psd propto r^{-alpha}$. Understanding this relation will shed light on the physics relevant to collisionless collapse and on the subsequent structures formed. We show that empirical density profiles that provide good fits to N-body halos also happen to have nearly scale-free psd distributions when in equilibrium. We have also done a preliminary investigation of variables other than $r$ that may match or supercede the correlation with $F$. In the same vein, we show that $rho/sigma^m$, where $m=3$ is the most appropriate combination to use in discussions of the power-law relationship. Since the mechanical equilibrium condition that characterizes the final systems does not by itself lead to power-law $F$ distributions, our findings prompt us to posit that dynamical collapse processes (such as violent relaxation) are responsible for the radial power-law nature of the psd distributions of virialized systems.
We present a parametric analysis of the intracluster medium and gravitating mass distribution of a statistical sample of 20 galaxy clusters using the phenomenological cluster model of Ascasibar and Diego. We describe an effective scheme for the estimation of errors on model parameters and derived quantities using bootstrap resampling. We find that the model provides a good description of the data in all cases and we quantify the mean fractional intrinsic scatter about the best-fit density and temperature profiles, finding this to have median values across the sample of 2 and 5 per cent, respectively. In addition, we demonstrate good agreement between r500 determined directly from the model and that estimated from a core-excluded global spectrum. We compare cool core and non-cool core clusters in terms of the logarithmic slopes of their gas density and temperature profiles and the distribution of model parameters and conclude that the two categories are clearly separable. In particular, we confirm the effectiveness of the logarithmic gradient of the gas density profile measured at 0.04 r500 in differentiating between the two types of cluster.
We apply four different mass modelling methods to a suite of publicly available mock data for spherical stellar systems. We focus on the recovery of the density and velocity anisotropy as a function of radius, using either line-of-sight velocity data only, or adding proper motion data. All methods perform well on isotropic and tangentially anisotropic mock data, recovering the density and velocity anisotropy within their 95% confidence intervals over the radial range 0.25 < R/Rhalf < 4, where Rhalf is the half light radius. However, radially-anisotropic mocks are more challenging. For line-of-sight data alone, only methods that use information about the shape of the velocity distribution function are able to break the degeneracy between the density profile and the velocity anisotropy to obtain an unbiased estimate of both. This shape information can be obtained through directly fitting a global phase space distribution function, by using higher order Virial Shape Parameters, or by assuming a Gaussian velocity distribution function locally, but projecting it self-consistently along the line of sight. Including proper motion data yields further improvements, and in this case, all methods give a good recovery of both the radial density and velocity anisotropy profiles.
We study the comoving relative velocities, v12, of model isolated galaxy pairs at z=0.5. For this purpose, we use the predictions from the GALFORM semi-analytical model of galaxy formation and evolution based on a Lambda cold dark matter cosmology consistent with the results from WMAP7. In real space, we find that isolated pairs of galaxies are predicted to form an angle t with the line-of-sight that is uniformily distributed as expected if the Universe is homogeneous and isotropic. We also find that isolated pairs of galaxies separated by a comoving distance between 1 and 3 Mpc/h are predicted to have <v12>=0. For galaxies in this regime, the distribution of the angle t is predicted to change minimally from real to redshift space, with a change smaller than 5% in <sin^2 t>. However, the distances defining the comoving regime strongly depends on the applied isolation criteria.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
