Theoretical astrophysics and plasma physics at RPC

Abstract:

The GAIA mission will produce an extraordinary database from which we should be able to deduce not only the Galaxy's current structure, but also much of its history, and thus cast a powerful light on the way in which galaxies in general are made up of components, and of how these formed. The database can be fully exploited only by fitting to it a sophisticated model of the entire Galaxy. Steady-state models are of fundamental importance even though the Galaxy cannot be in a steady state. A very elaborate model of the Galaxy will be required to reproduce the great wealth of detail that GAIA will reveal. A systematic approach to model-building will be required if such a model is to be successfully constructed, however. The natural strategy is to proceed through a series of models of ever increasing elaborateness, and to be guided in the specification of the next model by mismatches between the data and the current model. An approach to the dynamics of systems with steady gravitational potentials that we call the `torus programme' promises to provide an appropriate framework within which to carry out the proposed modelling programme. The basic principles of this approach have been worked out in some detail and are summarized here. Some extensions will be required before the GAIA database can be successfully confronted. Other modelling techniques that might be employed are briefly examined.

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 time-step shortens by a factor of 2 with each successive refinement. Competing approaches to N-body simulation are discussed from the point of view of the basic theory of N-body simulation. It is argued that an appropriate choice of softening length ∈ is of great importance and that ∈ should be at all points an appropriate multiple of the local interparticle separation. Unlike tree and P3M codes, multigrid codes automatically satisfy this requirement. We show that at early times and low densities in cosmological simulations, ∈ needs to be significantly smaller relative to the interparticle separation than in virialized regions. Tests of the ability of the code's Poisson solver to recover the gravitational fields of both virialized haloes 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 ACDM 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. The code is available from http://www-thphys.physics.ox.ac.uk/users/MLAPM.

Abstract:

Following the suggestion of Black that some massive extrasolar planets may be associated with the tail of the distribution of stellar companions, we investigate a scenario in which 5 ≤ N ≤ 100 planetary mass objects are assumed to form rapidly through a fragmentation process occuring in a disc or protostellar envelope on a scale of 100 au. These are assumed to have formed rapidly enough through gravitational instability or fragmentation that their orbits can undergo dynamical relaxation on a time-scale of ∼100 orbits. Under a wide range of initial conditions and assumptions, the relaxation process ends with either (i) one potential 'hot Jupiter' plus up to two 'external' companions, i.e. planets orbiting near the outer edge of the initial distribution; (ii) one or two 'external' planets or even none at all; (iii) one planet on an orbit with a semi-major axis of 10 to 100 times smaller than the outer boundary radius of the inital distribution together with an 'external' companion. Most of the other objects are ejected and could contribute to a population of free-floating planets. Apart from the potential 'hot Jupiters', all the bound objects are on orbits with high eccentricity, and also with a range of inclination with respect to the stellar equatorial plane. We found that, apart from the close orbiters, the probability of ending up with a planet orbiting at a given distance from the central star increases with the distance. This is because of the tendency of the relaxation process to lead to collisions with the central star. The scenario we envision here does not impose any upper limit on the mass of the planets. We discuss the application of these results to some of the more massive extrasolar planets.