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.     

N-body simulations for testing the stability of triaxial galaxies in MOND

We perform a stability test of triaxial models in Modified Newtonian Dynamics (MOND) using N -body simulations. The triaxial models considered here have densities that vary with   r ⁻¹  in the centre and   r ⁻⁴  at large radii. The total mass of the model varies from 10⁸ to  10¹⁰ M_⊙  , representing the mass scale of dwarfs to medium-mass elliptical galaxies, respectively, from deep MOND to quasi-Newtonian gravity. We build triaxial galaxy models using the Schwarzschild technique, and evolve the systems for 200 Keplerian dynamical times (at the typical length-scale of 1.0 kpc). We find that the systems are virial overheating, and in quasi-equilibrium with the relaxation taking approximately 5 Keplerian dynamical times (1.0 kpc). For all systems, the change of the inertial (kinetic) energy is less than 10 per cent (20 per cent) after relaxation. However, the central profile of the model is flattened during the relaxation and the (overall) axis ratios change by roughly 10 per cent within 200 Keplerian dynamical times (at 1.0 kpc) in our simulations. We further find that the systems are stable once they reach the equilibrium state.     

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.     

Post-Newtonian N-body simulations

We report on the first fully consistent conventional cluster simulation which includes terms up to the third-order post-Newtonian approximation. Numerical problems for treating extremely energetic binaries orbiting a single massive object are circumvented by employing the special 'wheel-spoke' regularization method of Zare which has not been used in large- N simulations before. Idealized models containing   N = 1 × 10⁵  particles of mass  1 M_⊙  with a central black hole (BH) of  300 M_⊙  have been studied on GRAPE-type computers. An initial half-mass radius of   r _h≃ 0.1  pc is sufficiently small to yield examples of relativistic coalescence. This is achieved by significant binary shrinkage within a density cusp environment, followed by the generation of extremely high eccentricities which are induced by Kozai cycles and/or resonant relaxation. More realistic models with white dwarfs and 10 times larger half-mass radii also show evidence of general relativity effects before disruption. An experimentation with the post-Newtonian terms suggests that reducing the time-scales for activating the different orders progressively may be justified for obtaining qualitatively correct solutions without aiming for precise predictions of the final gravitational radiation wave form. The results obtained suggest that the standard loss-cone arguments underestimate the swallowing rate in globular clusters containing a central BH.     

Effects of ram pressure on the gas distribution and star formation in the Large Magellanic Cloud

We use high-resolution N -body/smoothed particle hydrodynamics (SPH) simulations to study the hydrodynamical interaction between the Large Magellanic Cloud (LMC) and the hot halo of the Milky Way. We investigate whether ram pressure acting on the satellite's interstellar medium can explain the peculiarities observed in the H  i distribution and the location of the recent star formation activity.
Due to the present nearly edge-on orientation of the disc with respect to the orbital motion, compression at the leading edge can explain the high density region observed in H  i at the south-east border. In the case of a face-on disc (according to Mastropietro the LMC was moving almost face-on before the last perigalactic passage), ram pressure directed perpendicular to the disc produces a clumpy structure characterized by voids and high density filaments that resemble those observed by the Parkes H  i survey. As a consequence of the very recent edge-on motion, the Hα emission is mainly concentrated on the eastern side where 30 Doradus and most of the supergiant shells are located, although some Hα complexes form a patchy distribution on the entire disc. In this scenario, only the youngest stellar complexes show a progression in age along the leading border of the disc.     

Final stages of N-body star cluster encounters

We performed numerical simulations of star cluster encounters with Hernquist's treecode on a CRAY YMP-2E computer. We used different initial conditions (relative positions and velocities, cluster sizes, masses and concentration degrees) with the total number of particles per simulation ranging from 4608 to 20 480. Long-term interaction stages (up to 1 Gyr) when the pair coalesces into a single cluster are compared with isolated LMC clusters. Evidence is found that, when seen in a favourable plane, these resulting clusters show elliptical shapes as a result of the disruption of one of the companions. These elliptical shapes are essentially time-independent, but they do depend on the initial structural parameters of the pair components. We also analysed the fraction of stars that are ejected to the field by the interaction. We found that this fraction can be almost 50 per cent for the disrupted cluster. These simulations can represent a possible mechanism with which to explain the ellipticity observed in several star clusters in the Magellanic Clouds.     

The spatial evolution of stellar structures in the Large Magellanic Cloud

We present an analysis of the spatial distribution of various stellar populations within the Large Magellanic Cloud (LMC). We combine mid-infrared selected young stellar objects, optically selected samples with mean ages between ∼9 and ∼1000 Myr and existing stellar cluster catalogues to investigate how stellar structures form and evolve within the LMC. For the analysis we use Fractured Minimum Spanning Trees, the statistical Q parameter and the two-point correlation function. Restricting our analysis to young massive (OB) stars, we confirm our results obtained for M33, namely that the luminosity function of the groups is well described by a power law with index −2, and that there is no characteristic length-scale of star-forming regions. We find that stars in the LMC are born with a large amount of substructure, consistent with a two-dimensional fractal distribution with dimension     and evolve towards a uniform distribution on a time-scale of ∼175 Myr. This is comparable to the crossing time of the galaxy, and we suggest that stellar structure, regardless of spatial scale, will be eliminated in a crossing time. This may explain the smooth distribution of stars in massive/dense young clusters in the Galaxy, while other, less massive, clusters still display large amounts of structure at similar ages. By comparing the stellar and star cluster distributions and evolving time-scales, we show that infant mortality of clusters (or 'popping clusters') has a negligible influence on the galactic structure. Finally, we quantify the influence of the elongation, differential extinction and contamination of a population on the measured Q value.