High-performance direct gravitational N-body simulations on graphics processing units

We present the results of gravitational direct N-body simulations using the commercial graphics processing units (GPU) NVIDIA Quadro FX1400 and GeForce 8800GTX, and compare the results with GRAPE-6Af special purpose hardware. The force evaluation of the N-body problem was implemented in Cg using the GPU directly to speed-up the calculations. The integration of the equations of motions were, running on the host computer, implemented in C using the 4th order predictor–corrector Hermite integrator with block time steps. We find that for a large number of particles (N  10⁴) modern graphics processing units offer an attractive low cost alternative to GRAPE special purpose hardware. A modern GPU continues to give a relatively flat scaling with the number of particles, comparable to that of the GRAPE. The GRAPE is designed to reach double precision, whereas the GPU is intrinsically single-precision. For relatively large time steps, the total energy of the N-body system was conserved better than to one in 10⁶ on the GPU, which is impressive given the single-precision nature of the GPU. For the same time steps, the GRAPE gave somewhat more accurate results, by about an order of magnitude. However, smaller time steps allowed more energy accuracy on the grape, around 10⁻¹¹, whereas for the GPU machine precision saturates around 10⁻⁶ For N  10⁶ the GeForce 8800GTX was about 20 times faster than the host computer. Though still about a factor of a few slower than GRAPE, modern GPUs outperform GRAPE in their low cost, long mean time between failure and the much larger onboard memory; the GRAPE-6Af holds at most 256k particles whereas the GeForce 8800GTX can hold 9 million particles in memory.     

Physical parameters of the continuum emission layer of Mira type stars. II. Radii and expansion velocities

The Pulkovo Spectrophotometric Data Base is used to determine the expansion velocities and radii of the continuum emission layers of the stars RR SCO and OMICR CET. Dependences of the radii and velocities on the phase of the variability cycle are found which agree with the known values obtained from the layers which emit line radiation. __________ Translated from Astrofizika, Vol. 50, No. 4, pp. 557–564 (November 2007).     

Virialization of cosmological structures in models with time-varying equation of state

We study the virialization of the cosmic structures in the framework of flat cosmological models where the dark energy component plays an important role in the global dynamics of the Universe. In particular, our analysis focuses on the study of the spherical matter perturbations, as the latter decouple from the background expansion, start to 'turn around' and finally collapse. We generalize this procedure, taking into account models with an equation of state which vary with time, and provide a complete formulation of the cluster virialization attempting to address the non-linear regime of structure formation. In particular, assuming that clusters have collapsed prior to the epoch of z _f≃ 1.4, in which the most distant cluster has been found, we show that the behaviour of the spherical collapse model depends on the functional form of the equation of state.