We present Sapporo, a library for performing high-precision gravitational N-body simulations on NVIDIA Graphical Processing Units (GPUs). Our library mimics the GRAPE-6 library, and N-body codes currently running on GRAPE-6 can switch to Sapporo by a
simple relinking of the library. The precision of our library is comparable to that of GRAPE-6, even though internally the GPU hardware is limited to single precision arithmetics. This limitation is effectively overcome by emulating double precision for calculating the distance between particles. The performance loss of this operation is small (< 20%) compared to the advantage of being able to run at high precision. We tested the library using several GRAPE-6-enabled N-body codes, in particular with Starlab and phiGRAPE. We measured peak performance of 800 Gflop/s for running with 10^6 particles on a PC with four commercial G92 architecture GPUs (two GeForce 9800GX2). As a production test, we simulated a 32k Plummer model with equal mass stars well beyond core collapse. The simulation took 41 days, during which the mean performance was 113 Gflop/s. The GPU did not show any problems from running in a production environment for such an extended period of time.
We consider the effect of mass segregation on the observable integrated properties of star clusters. The measurable properties depend on a combination of the dynamical age of the cluster and the physical age of the stars in the cluster. To investigat
e all possible combinations of these two quantities we propose an analytical model for the mass function of segregated star clusters that agrees with the results of N-body simulations, in which any combination can be specified. For a realistic degree of mass segregation and a fixed density profile we find with increasing age an increase in the measured core radii and a central surface brightness that decreases in all filters more rapidly than what is expected from stellar evolution alone. Within a Gyr the measured core radius increases by a factor of two and the central surface density in all filters of a segregated cluster will be overestimated by a similar factor when not taking into account mass segregation in the conversion from light to mass. We find that the $V-I$ colour of mass segregated clusters decreases with radius by about 0.1-0.2 mag, which could be observable. From recent observations of partially resolved extra-galactic clusters a decreasing half-light radius with increasing wavelength was observed, which was attributed to mass segregation. These observations can not be reproduced by our models. We find that the differences between measured radii in different filters are always smaller than 5%.
The early evolution of dense star clusters is possibly dominated by close interactions between stars, and physical collisions between stars may occur quite frequently. Simulating a stellar collision event can be an intensive numerical task, as detail
ed calculations of this process require hydrodynamic simulations in three dimensions. We present a computationally inexpensive method in which we approximate the merger process, including shock heating, hydrodynamic mixing and mass loss, with a simple algorithm based on conservation laws and a basic qualitative understanding of the hydrodynamics of stellar mergers. The algorithm relies on Archimedes principle to dictate the distribution of the fluid in the stable equilibrium situation. We calibrate and apply the method to mergers of massive stars, as these are expected to occur in young and dense star clusters. We find that without the effects of microscopic mixing, the temperature and chemical composition profiles in a collision product can become double-valued functions of enclosed mass. Such an unphysical situation is mended by simulating microscopic mixing as a post-collision effect. In this way we find that head-on collisions between stars of the same spectral type result in substantial mixing, while mergers between stars of different spectral type, such as type B and O stars ($sim$10 and $sim$40msun respectively), are subject to relatively little hydrodynamic mixing.
