Theoretical astrophysics and plasma physics at RPC

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.

Abstract:

We analyze the nuclear stellar dynamics of the SBO galaxy NGC 1023, utilizing observational data both from the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope and from the ground. The stellar kinematics measured from these long-slit spectra show rapid rotation (V ≈ 70 km s^-1 at a distance of O″. 1 = 4.9 pc from the nucleus) and increasing velocity dispersion toward the nucleus (where σ = 295 ± 30 km s^-1). We model the observed stellar kinematics assuming an axisymmetric mass distribution with both two and three integrals of motion. Both modeling techniques point to the presence of a central dark compact mass (which presumably is a supermassive black hole) with confidence greater than 99%. The isotropic two-integral models yield a best-fitting black hole mass of (6.0 ± 1.4) x 10⁷ M⊙ and mass-to-light ratio (M/LV) of 5.38 ± 0.08, and the goodness of fit (χ²) is insensitive to reasonable values for the galaxy's inclination. The three-integral models, which nonparametrically fit the observed line-of-sight velocity distribution as a function of position in the galaxy, suggest a black hole mass of (3.9 ± 0.4) x 10⁷ M⊙ and M/LV of 5.56 ± 0.02 (internal errors), and the edge-on models are vastly superior fits over models at other inclinations. The internal dynamics in NGC 1023 as suggested by our best-fit three-integral model shows that the velocity distribution function at the nucleus is tangentially anisotropic, suggesting the presence of a nuclear stellar disk. The nuclear line-of-sight velocity distribution has enhanced wings at velocities ≥ 600 km s^-1 from systemic, suggesting that perhaps we have detected a group of stars very close to the central dark mass.