Theoretical astrophysics and plasma physics at RPC

Abstract:

We consider long-slit emission-line spectra of galactic nuclei when the slit is wider than the instrumental point spread function, and the target has large velocity gradients. The finite width of the slit generates complex distributions of brightness at a given spatial point in the measured spectrum, which can be misinterpreted as coming from additional physically distinct nuclear components. We illustrate this phenomenon for the case of a thin disc in circular motion around a nuclear black hole (BH). We develop a new method for estimating the mass of the BH that exploits a feature in the spectrum at the outer edge of the BH's sphere of influence, and therefore gives higher sensitivity to BH detection than traditional methods. Moreover, with this method we can determine the BH mass and the inclination of the surrounding disc separately, whereas the traditional approach to BH estimation requires two long-slit spectra to be taken. We show that, with a given spectrograph, the detectability of a BH depends on the sense of rotation of the nuclear disc. We apply our method to estimate the BH mass in M84 from a publicly available spectrum, and recover a value four times lower than that published previously from the same data.

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.