Theoretical astrophysics and plasma physics at RPC

Abstract:

It is logically possible that early two-body relaxation in simulations of cosmological clustering influences the final structure of massive clusters. Convergence studies in which mass and spatial resolution are simultaneously increased, cannot eliminate this possibility. We test the importance of two-body relaxation in cosmological simulations with simulations in which there are two species of particles. The cases of two mass ratios, sqrt(2):1 and 4:1, are investigated. Simulations are run with both a spatially fixed softening length and adaptive softening using the publicly available codes GADGET and MLAPM, respectively. The effects of two-body relaxation are detected in both the density profiles of halos and the mass function of halos. The effects are more pronounced with a fixed softening length, but even in this case they are not so large as to suggest that results obtained with one mass species are significantly affected by two-body relaxation. The simulations that use adaptive softening are less affected by two-body relaxation and produce slightly higher central densities in the largest halos. They run about three times faster than the simulations that use a fixed softening length.

Abstract:

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.

Abstract:

We present a computer code written in C that is designed to simulate structure formation from collisionless matter. The code is purely grid-based and uses a recursively refined Cartesian grid to solve Poisson's equation for the potential, rather than obtaining the potential from a Green's function. Refinements can have arbitrary shapes and in practice closely follow the complex morphology of the density field that evolves. The timestep shortens by a factor two with each successive refinement. It is argued that an appropriate choice of softening length is of great importance and that the softening should be at all points an appropriate multiple of the local inter-particle separation. Unlike tree and P3M codes, multigrid codes automatically satisfy this requirement. We show that at early times and low densities in cosmological simulations, the softening needs to be significantly smaller relative to the inter-particle separation than in virialized regions. Tests of the ability of the code's Poisson solver to recover the gravitational fields of both virialized halos and Zel'dovich waves are presented, as are tests of the code's ability to reproduce analytic solutions for plane-wave evolution. The times required to conduct a LCDM cosmological simulation for various configurations are compared with the times required to complete the same simulation with the ART, AP3M and GADGET codes. The power spectra, halo mass functions and halo-halo correlation functions of simulations conducted with different codes are compared.

Abstract:

For two decades the steady-state cooling-flow model has dominated the literature of cluster and elliptical-galaxy X-ray sources. For ten years this model has been in severe difficulty from a theoretical point of view, and it is now coming under increasing pressure observationally. For two decades the steady-state cooling-flow model has dominated the literature of cluster and elliptical-galaxy X-ray sources. For ten years this model has been in severe difficulty from a theoretical point of view, and it is now coming under increasing pressure observationally. A small number of enthusiasts have argued for a radically different interpretation of the data, but had little impact on prevailing opinion because the unsteady heating picture that they advocate is extremely hard to work out in detail. Here I explain why it is difficult to extract robust observational predictions from the heating picture. Major problems include the variability of the sources, the different ways in which a bi-polar flow can impact on X-ray emission, the weakness of synchrotron emission from sub-relativistic flows, and the sensitivity of synchrotron emission to a magnetic field that is probably highly localized.