Energy ejection in the collapse of a cold spherical self-gravitating cloud

When an open system of classical point particles interacting by Newtonian gravity collapses and relaxes violently, an arbitrary amount of energy may, in principle, be carried away by particles which escape to infinity. We investigate here, using numerical simulations, how this released energy and other related quantities (notably the binding energy and size of the virialized structure) depend on the initial conditions, for the one-parameter family of starting configurations given by randomly distributing N cold particles in a spherical volume. Previous studies have established that the minimal size reached by the system scales approximately as   N ^1/3  , a behaviour which follows trivially when the growth of perturbations (which regularize the singularity of the cold collapse in the   N  →∞  limit) is assumed to be unaffected by the boundaries. Our study shows that the energy ejected grows approximately in proportion to   N ^1/3  , while the fraction of the initial mass ejected grows only very slowly with N , approximately logarithmically, in the range of N simulated. We examine in detail the mechanism of this mass and energy ejection, showing explicitly that it arises from the interplay of the growth of perturbations with the finite size of the system. A net lag of particles compared to their uniform spherical collapse trajectories develops first at the boundaries and then propagates into the volume during the collapse. Particles in the outer shells are then ejected as they scatter through the time-dependent potential of an already re-expanding central core. Using modified initial configurations, we explore the importance of fluctuations at different scales and discreteness (i.e. non-Vlasov) effects in the dynamics.     

Multimass spherical structure models for N-body simulations

We present a simple and efficient method to set up spherical structure models for N -body simulations with a multimass technique. This technique reduces by a substantial factor the computer run time needed in order to resolve a given scale as compared to single-mass models. It therefore allows to resolve smaller scales in N -body simulations for a given computer run time. Here, we present several models with an effective resolution of up to  1.68 × 10⁹  particles within their virial radius which are stable over cosmologically relevant time-scales. As an application, we confirm the theoretical prediction by Dehnen that in mergers of collisionless structures like dark matter haloes always the cusp of the steepest progenitor is preserved. We model each merger progenitor with an effective number of particles of approximately 10⁸ particles. We also find that in a core–core merger the central density approximately doubles whereas in the cusp–cusp case the central density only increases by approximately 50 per cent. This may suggest that the central regions of flat structures are better protected and get less energy input through the merger process.     

An iterative method for constructing equilibrium phase models of stellar systems

We present a new method for constructing equilibrium phase models for stellar systems, which we call the iterative method. It relies on constrained, or guided evolution, so that the equilibrium solution has a number of desired parameters and/or constraints. This method is very powerful, to a large extent due to its simplicity. It can be used for mass distributions with an arbitrary geometry and a large variety of kinematical constraints. We present several examples illustrating it. Applications of this method include the creation of initial conditions for N -body simulations and the modelling of galaxies from their photometric and kinematic observations.     

Constrained simulations of the local universe – II. The nature of the local Hubble flow

Using a suite of N -body simulations in different cold dark matter (CDM) scenarios, with cosmological constant (ΛCDM) and without (OCDM, SCDM), we study the Hubble flow (σ_H) in Local Volumes (LV) around Local Group (LG) like objects found in these simulations, and compare the numerical results with the most recent observations. We show that ΛCDM and OCDM models exhibit the same behaviour of σ_H. Hence, we demonstrate that the observed coldness of the Hubble flow is not likely to be a manifestation of the dark energy, contrary to previous claims. The coldness does not constitute a problem by itself but it poses a problem to the standard ΛCDM model only if the mean density within the LV is greater than twice the mean matter cosmic density. The lack of blueshifted galaxies in the LV, outside of the LG can be considered as another manifestation of the coldness of the flow. Finally, we show that the main dynamical parameter that affects the coldness of the flow is the relative isolation of the LG, and the absence of nearby Milky Way like objects within a distance of about  3 Mpc  .     

An energy-conserving formalism for adaptive gravitational force softening in smoothed particle hydrodynamics and N-body codes

In this paper, we describe an adaptive softening length formalism for collisionless N -body and self-gravitating smoothed particle hydrodynamics (SPH) calculations which conserves momentum and energy exactly. This means that spatially variable softening lengths can be used in N -body calculations without secular increases in energy. The formalism requires the calculation of a small additional term to the gravitational force related to the gradient of the softening length. The extra term is similar in form to the usual SPH pressure force (although opposite in direction) and is therefore straightforward to implement in any SPH code at almost no extra cost. For N -body codes, some additional cost is involved as the formalism requires the computation of the density through a summation over neighbouring particles using the smoothing kernel. The results of numerical tests demonstrate that, for homogeneous mass distributions, the use of adaptive softening lengths gives a softening which is always close to the 'optimal' choice of fixed softening parameter, removing the need for fine-tuning. For a heterogeneous mass distribution (as may be found in any large-scale N -body simulation), we find that the errors on the least-dense component are lowered by an order of magnitude compared to the use of a fixed softening length tuned to the densest component. For SPH codes, our method presents a natural and an elegant choice of softening formalism which makes a small improvement to both the force resolution and the total energy conservation at almost zero additional cost.     

Observable consequences of planet formation models in systems with close-in terrestrial planets

To date, two planetary systems have been discovered with close-in, terrestrial-mass planets     . Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form 'hot Earths' and their observable predictions. Models include: (1) in situ accretion; (2) formation at larger orbital distance followed by inward 'type 1' migration; (3) formation from material being 'shepherded' inward by a migrating gas giant planet; (4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disc; (5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances and (6) photoevaporative mass-loss of a close-in giant planet. Models 1–4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets     with very close-in perihelion distances (≲0.025 au), and even then the net inward movement in orbital distance is at most only 0.1–0.15 au. For planets of less than     , photoevaporation can remove the planet's envelope and leave behind the solid core on a Gyr time-scale, but only for planets inside 0.025–0.05 au. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet's radius).     

nmagic: a fast parallel implementation of a χ2-made-to-measure algorithm for modelling observational data

We describe a made-to-measure (M2M) algorithm for constructing N -particle models of stellar systems from observational data (χ²M2M), extending earlier ideas by Syer & Tremaine. The algorithm properly accounts for observational errors, is flexible, and can be applied to various systems and geometries. We implement this algorithm in a parallel code nmagic and carry out a sequence of tests to illustrate its power and performance. (i) We reconstruct an isotropic Hernquist model from density moments and projected kinematics and recover the correct differential energy distribution and intrinsic kinematics. (ii) We build a self-consistent oblate three-integral maximum rotator model and compare how the distribution function is recovered from integral field and slit kinematic data. (iii) We create a non-rotating and a figure rotating triaxial stellar particle model, reproduce the projected kinematics of the figure rotating system by a non-rotating system of the same intrinsic shape, and illustrate the signature of pattern rotation in this model. From these tests, we comment on the dependence of the results from χ²M2M on the initial model, the geometry, and the amount of available data.     

High-inclination planets and asteroids in multistellar systems

The Kozai mechanism often destabilizes high-inclination orbits. It couples changes in the eccentricity and inclination, and drives high inclination, circular orbits to low inclination, eccentric orbits. In a recent study of the dynamics of planetesimals in the quadruple star system HD 98800, there were significant numbers of stable particles in circumbinary polar orbits about the inner binary pair which are apparently able to evade the Kozai instability.
Here, we isolate this feature and investigate the dynamics through numerical and analytical models. The results show that the Kozai mechanism of the outer star is disrupted by a nodal libration induced by the inner binary pair on a shorter time-scale. By empirically modelling the period of the libration, a criteria for determining the high-inclination stability limits in general triple systems is derived. The nodal libration feature is interesting and, although affecting inclination and node only, shows many parallels to the Kozai mechanism. This raises the possibility that high-inclination planets and asteroids may be able to survive in multistellar systems.     

Multimass schemes for collisionless N-body simulations

We present a general scheme for constructing Monte Carlo realizations of equilibrium, collisionless galaxy models with known distribution function (DF) f ₀. Our method uses importance sampling to find the sampling DF f _s that minimizes the mean-square formal errors in a given set of projections of the DF f ₀. The result is a multimass N -body realization of the galaxy model in which 'interesting' regions of phase space are densely populated by lots of low-mass particles, increasing the effective N there, and less interesting regions by fewer, higher mass particles.
As a simple application, we consider the case of minimizing the shot noise in estimates of the acceleration field for an N -body model of a spherical Hernquist model. Models constructed using our scheme easily yield a factor of ∼100 reduction in the variance at the central acceleration field when compared to a traditional equal-mass model with the same number of particles. When evolving both models with a real N -body code, the diffusion coefficients in our model are reduced by a similar factor. Therefore, for certain types of problems, our scheme is a practical method for reducing the two-body relaxation effects, thereby bringing the N -body simulations closer to the collisionless ideal.     

Are mergers responsible for universal halo properties?

N -body simulations of cold dark matter (CDM) have shown that, in this hierarchical structure formation model, dark matter halo properties, such as the density profile, the phase-space density profile, the distribution of axial ratio, the distribution of spin parameter and the distribution of internal specific angular momentum, follow 'universal' laws or distributions. Here, we study the properties of the first generation of haloes in a hot dark matter (HDM) dominated universe, as an example of halo formation through monolithic collapse. We find all these universalities to be present in this case also. Halo density profiles are very well fit by the Navarro, Frenk & White profile over two orders of magnitude in mass. The concentration parameter depends on mass as   c ∝ M ^0.2  , reversing the dependence found in a hierarchical CDM universe. However, the concentration–formation time relation is similar in the two cases: earlier forming haloes tend to be more concentrated than their later forming counterparts. Halo formation histories are also characterized by two phases in the HDM case: an early phase of rapid accretion followed by slower growth. Furthermore, there is no significant difference between the HDM and CDM cases concerning the statistics of other halo properties: the phase-space density profile; the velocity anisotropy profile; the distribution of shape parameters; the distribution of spin parameter and the distribution of internal specific angular momentum are all similar in the two cases. Only substructure content differs dramatically. These results indicate that mergers do not play a pivotal role in establishing the universalities, thus contradicting models which explain them as consequences of mergers.