Cosmological perturbation theory and the spherical collapse model — III. The velocity divergence field and the Ω dependence

Cosmological perturbation theory (PT) is a useful tool to study the cumulants of the density and velocity fields in the large-scale structure of the Universe. In Papers I and II of this series we saw that the spherical collapse (SC) model provides the exact solution to PT at tree-level and gives a good approximation to the loop corrections (next-to-leading orders), indicating negligible tidal effects. Here, we derive predictions for the (smoothed) cumulants of the velocity divergence field θ ≡ ▽ ⊙  v for an irrotational fluid in the SC model. By comparing these with the exact analytic results of Scoccimarro &38; Frieman, it is shown that, at least for the unsmoothed case, the loop corrections to the cumulants of θ are dominated by tidal effects. However, most of the tidal contribution seems to cancel out when computing the hierarchical ratios, T _J  = 〈θ ^J 〉 / 〈θ²〉  ^{J −1}. We also extend the work presented in Papers I and II to give predictions for the cumulants of the density and velocity divergence fields in non-flat spaces. In particular, we show the equivalence between the spherically symmetric solution to the equations of motion in the SC model (given in terms of the density) and that of the Lagrangian PT approach (given in terms of the displacement field). It is shown that the Ω dependence is very weak for both cosmic fields even at one loop (a 10 per cent effect at most), except for the overall factor f (Ω) that couples to the velocity divergence.     

Monte-Carlo Simulations of Cosmic-ray and Internal Radiation Effects on ISOCAM on Board ISO

One of the main limitations to the sensitivity of the infrared camera ISOCAM on-board the Infrared Space Observatory (ISO) comes from responsivity variations and glitches caused by the impacts of charged particles in photo-detectors. Glitch rate measurements, glitch properties and removal methods have already been addressed during the first ISO detector workshop(Madrid, 1998) and published in a special issue of Experimental Astronomy. It appeared that glitch rate and most of glitch properties could be reproduced by Monte-Carlo simulations. This is very interesting in order to predict before launch the effect of charged particles in photo-detectors operated in space. This paper presents results of Monte-Carlo simulations of radiation effects on ISOCAM detectors. Glitch rates, spatial and energetic properties of glitches have been computed and are compared with measured values. This revised version was published online in July 2006 with corrections to the Cover Date.     

Evolution of binary stars and the effect of tides on binary populations

We present a rapid binary-evolution algorithm that enables modelling of even the most complex binary systems. In addition to all aspects of single-star evolution, features such as mass transfer, mass accretion, common-envelope evolution, collisions, supernova kicks and angular momentum loss mechanisms are included. In particular, circularization and synchronization of the orbit by tidal interactions are calculated for convective, radiative and degenerate damping mechanisms. We use this algorithm to study the formation and evolution of various binary systems. We also investigate the effect that tidal friction has on the outcome of binary evolution. Using the rapid binary code, we generate a series of large binary populations and evaluate the formation rate of interesting individual species and events. By comparing the results for populations with and without tidal friction, we quantify the hitherto ignored systematic effect of tides and show that modelling of tidal evolution in binary systems is necessary in order to draw accurate conclusions from population synthesis work. Tidal synchronism is important but, because orbits generally circularize before Roche lobe overflow, the outcome of the interactions of systems with the same semilatus rectum is almost independent of eccentricity. It is not necessary to include a distribution of eccentricities in population synthesis of interacting binaries; however, the initial separations should be distributed according to the observed distribution of semilatera recta rather than periods or semimajor axes.     

The formation of ultracompact dwarf galaxies

Recent spectroscopic observations of galaxies in the Fornax Cluster reveal nearly unresolved 'star-like' objects with redshifts appropriate to the Fornax Cluster. These objects have intrinsic sizes of ≈100 pc and absolute B -band magnitudes in the range  −14< M _B<−11.5 mag  and lower limits for the central surface brightness   μ _B≳23 mag arcsec⁻²  , and so appear to constitute a new population of ultracompact dwarf galaxies (UCDs). Such compact dwarfs were predicted to form from the amalgamation of stellar superclusters (by Kroupa) , which are rich aggregates of young massive star clusters (YMCs) that can form in collisions between gas-rich galaxies. Here we present the evolution of superclusters in a tidal field. The YMCs merge on a few supercluster crossing times. Superclusters that are initially as concentrated and massive as knot S in the interacting Antennae galaxies evolve to merger objects that are long-lived and show properties comparable to the newly discovered UCDs. Less massive superclusters resembling knot 430 in the Antennae may evolve to ω Cen-type systems. Low-concentration superclusters are disrupted by the tidal field, dispersing their surviving star clusters while the remaining merger objects rapidly evolve into the   μ _B− M _B  region populated by low-mass Milky Way dSph satellites.     

The Monitor project: JW 380 – a 0.26-, 0.15-M⊙, pre-main-sequence eclipsing binary in the Orion nebula cluster

We present a general recipe for constructing N -body realizations of galaxies comprising near spherical and disc components. First, an exact spherical distribution function for the spheroids (halo and bulge) is determined, such that it is in equilibrium with the gravitational monopole of the disc components. Second, an N -body realization of this model is adapted to the full disc potential by growing the latter adiabatically from its monopole. Finally, the disc is sampled with particles drawn from an appropriate distribution function, avoiding local-Maxwellian approximations. We performed test simulations and find that the halo and bulge radial density profile very closely match their target model, while they become slightly oblate due to the added disc gravity. Our findings suggest that vertical thickening of the initially thin disc is caused predominantly by spiral and bar instabilities, which also result in a radial re-distribution of matter, rather than scattering off interloping massive halo particles.     

Mapping on the HEALPix grid

The natural spherical projection associated with the Hierarchical Equal-Area and Isolatitude Pixelization (HEALPix) is described and shown to be one of a hybrid class that combines the cylindrical equal-area and Collignon projections, not previously documented in the cartographic literature. Projection equations are derived for the class in general and are used to investigate its properties. It is shown that the HEALPix projection suggests a simple method of (i) storing and (ii) visualizing data sampled on the grid of the HEALPix pixelization, and also suggests an extension of the pixelization that is better suited for these purposes. Potentially useful properties of other members of the class are described, and new triangular and hexagonal pixelizations are constructed from them. Finally, the standard formalism is defined for representing the celestial coordinate system for any member of the class in the FITS data format.     

grommet: an N-body code for high-resolution simulations of individual galaxies

This paper presents a fast, economical particle-multiple-mesh N -body code optimized for large- N modelling of collisionless dynamical processes, such as black hole wandering or bar–halo interactions, occurring within isolated galaxies. The code has been specially designed to conserve linear momentum. Despite this, it also has variable softening and an efficient block-time-step scheme: the force between any pair of particles is calculated using the finest mesh that encloses them both (respecting Newton's third law) and is updated only on the longest time-step of the two (which conserves momentum). For realistic galaxy models with   N ≳ 10⁶  , it is faster than the fastest comparable tree code by factors ranging from ∼2 (using single time-steps) to ∼10 (multiple time-steps in a concentrated galaxy).