Most recent progress in understanding the dynamical evolution of star clusters relies on direct N-body simulations. Owing to the computational demands, and the desire to model more complex and more massive star clusters, hardware calculational accele
rators, such as GRAPE special-purpose hardware or, more recently, GPUs (i.e. graphics cards), are generally utilised. In addition, simulations can be accelerated by adjusting parameters determining the calculation accuracy (i.e. changing the internal simulation time step used for each star). We extend our previous thorough comparison (Anders et al. 2009) of basic quantities as derived from simulations performed either with STARLAB/KIRA or NBODY6. Here we focus on differences arising from using different hardware accelerations (including the increasingly popular graphic card accelerations/GPUs) and different calculation accuracy settings. We use the large number of star cluster models (for a fixed stellar mass function, without stellar/binary evolution, primordial binaries, external tidal fields etc) already used in the previous paper, evolve them with STARLAB/KIRA (and NBODY6, where required), analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison, we apply the bootstrap algorithm for functional dependencies developed in our previous study. In general we find very high comparability of the simulation results, independent of the used computer hardware (including the hardware accelerators) and the used N-body code. For the tested accuracy settings we find that for reduced accuracy (i.e. time step at least a factor 2.5 larger than the standard setting) most simulation results deviate significantly from the results using standard settings. The remaining deviations are comprehensible and explicable.
Evolutionary synthesis models are the prime method to construct models of stellar populations, and to derive physical parameters from observations. One of the assumptions for such models so far has been the time-independence of the stellar mass funct
ion. However, dynamical simulations of star clusters in tidal fields have shown the mass function to change due to the preferential removal of low-mass stars from clusters. Here we combine the results from dynamical simulations of star clusters in tidal fields with our evolutionary synthesis code GALEV to extend the models by a new dimension: the total cluster disruption time. We reanalyse the mass function evolution found in N-body simulations of star clusters in tidal fields, parametrise it as a function of age and total cluster disruption time and use this parametrisation to compute GALEV models as a function of age, metallicity and the total cluster disruption time. We study the impact of cluster dissolution on the colour (generally, they become redder) and magnitude (they become fainter) evolution of star clusters, their mass-to-light ratios (off by a factor of ~2 -- 4 from standard predictions), and quantify the effect on the cluster age determination from integrated photometry (in most cases, clusters appear to be older than they are, between 20 and 200%). By comparing our model results with observed M/L ratios for old compact objects in the mass range 10^4.5 -- 10^8 Msun, we find a strong discrepancy for objects more massive than 10^7 Msun (higher M/L). This could be either caused by differences in the underlying stellar mass function or be an indication for the presence of dark matter in these objects. Less massive objects are well represented by the models. The models for a range of total cluster disruption times are available online. (shortened)
N-body simulations are widely used to simulate the dynamical evolution of a variety of systems, among them star clusters. Much of our understanding of their evolution rests on the results of such direct N-body simulations. They provide insight in the
structural evolution of star clusters, as well as into the occurrence of stellar exotica. Although the major pure N-body codes STARLAB/KIRA and NBODY4 are widely used for a range of applications, there is no thorough comparison study yet. Here we thoroughly compare basic quantities as derived from simulations performed either with STARLAB/KIRA or NBODY4. We construct a large number of star cluster models for various stellar mass function settings (but without stellar/binary evolution, primordial binaries, external tidal fields etc), evolve them in parallel with STARLAB/KIRA and NBODY4, analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison we develop a bootstrap algorithm for functional dependencies. We find an overall excellent agreement between the codes, both for the clusters structural and energy parameters as well as for the properties of the dynamically created binaries. However, we identify small differences, like in the energy conservation before core collapse and the energies of escaping stars, which deserve further studies. Our results reassure the comparability and the possibility to combine results from these two major N-body codes, at least for the purely dynamical models (i.e. without stellar/binary evolution) we performed. (abridged)
The study of young star cluster (YSC) systems, preferentially in starburst and merging galaxies, has seen great interest in the recent past, as it provides important input to models of star formation. However, even some basic properties (like the lum
inosity function [LF]) of YSC systems are still under debate. Here we study the photometric properties of the YSC system in the nearest major merger system, the Antennae galaxies. We find evidence for the existence of a statistically significant turnover in the LF.
Two aspects of our recent N-body studies of star clusters are presented: (1) What impact does mass segregation and selective mass loss have on integrated photometry? (2) How well compare results from N-body simulations using NBODY4 and STARLAB/KIRA?
