From direct N-body simulations we find that the dynamical evolution of star clusters is strongly influenced by the Roche volume filling factor. We present a parameter study of the dissolution of open star clusters with different Roche volume filling factors and different particle numbers. We study both Roche volume underfilling and overfilling models and compare with the Roche volume filling case. We find that in the Roche volume overfilling limit of our simulations two-body relaxation is no longer the dominant dissolution mechanism but the changing cluster potential. We call this mechnism "mass-loss driven dissolution" in contrast to "two-body relaxation driven dissolution" which occurs in the Roche volume underfilling regime. We have measured scaling exponents of the dissolution time with the two-body relaxation time. In this experimental study we find a decreasing scaling exponent with increasing Roche volume filling factor. The evolution of the escaper number in the Roche volume overfilling limit can be described by a log-logistic differential equation. We report the finding of a resonance condition which may play a role for the evolution of star clusters and may be calibrated by the main periodic orbit in the large island of retrograde quasiperiodic orbits in the Poincar\'e surfaces of section. We also report on the existence of a stability curve which may be of relevance with respect to the structure of star clusters.

We investigate the dynamics in the close vicinity of and within the critical area in a 2D effective galactic potential with a bar of Zotos. We have calculated Poincar\'e surfaces of section and the basins of escape. In both the Poincar\'e surfaces of section and the basins of escape we find numerical evidence for the existence of a separatrix which hinders orbits from escaping out of the bar region. We present numerical evidence for the similarity between spiral arms of barred spiral galaxies and tidal tails of star clusters.

Context: Using the recent observational data of R\"oser et al. we present $N$-body simulations of the Hyades open cluster. Aims: We make an attempt to determine initial conditions of the Hyades cluster at the time of its formation in order to reproduce the present-day cumulative mass profile, stellar mass and luminosity function (LF). Methods: We performed direct $N$-body simulations of the Hyades in an analytic Milky Way potential that account for stellar evolution and include primordial binaries in a few models. Furthermore, we applied a Kroupa (2001) IMF and used extensive ensemble-averaging. Results: We find that evolved single-star King initial models with King parameters $W_0 = 6-9$ and initial particle numbers $N_0 = 3000$ provide good fits to the observational present-day cumulative mass profile within the Jacobi radius. The best-fit King model has an initial mass of $1721\ M_\odot$ and an average mass loss rate of $-2.2 \ M_\odot/\mathrm{Myr}$. The K-band LFs of models and observations show a reasonable agreement. Mass segregation is detected in both observations and models. If 33% primordial binaries are included the initial particle number is reduced by 5% as compared to the model without primordial binaries. Conclusions: The present-day properties of the Hyades can be well reproduced by a standard King or Plummer initial model when choosing appropriate initial conditions. The degeneracy of good-fitting models can be quite high due to the large dimension of the parameter space. More simulations with different Roche-lobe filling factors and primordial binary fractions are required to explore this degeneracy in more detail.

We examine the ratios $r_h/r_J$ of projected half-mass and Jacobi radius as well as $r_t/r_J$ of tidal and Jacobi radius for open and globular clusters in the Milky Way using data of both observations and simulations. We applied an improved calculation of $r_J$ for eccentric orbits of globular clusters. A sample of 236 open clusters of Piskunov et al. within the nearest kiloparsec around the Sun has been used. For the Milky Way globular clusters, data are taken from the Harris catalogue. We particularly use the subsample of 38 Milky Way globular clusters for which orbits have been integrated by Dinescu et al. We aim to quantify the differences between open and globular clusters and to understand, why they form two intrinsically distinct populations. We find under certain assumptions, or, in other words, in certain approximations, (i) that globular clusters are presently Roche volume underfilling and (ii) with at least $3\sigma$ confidence that the ratio $r_h/r_J$ of half-mass and Jacobi radius is $3 - 5$ times larger at present for an average open cluster in our sample than for an average globular cluster in our sample and (iii) that a significant fraction of globular clusters may be Roche volume overfilling at pericentre with $r_t > r_J$. Another aim of this paper is to throw light on the underlying theoretical reason for the existence of the van den Bergh correlation between half-mass and galactocentric radius.

Dynamical friction leads to an orbital decay of massive objects like young compact star clusters or Massive Black Holes in central regions of galaxies. The dynamical friction force can be well approximated by Chandrasekhar's standard formula, but recent investigations show, that corrections to the Coulomb logarithm are necessary. With a large set of N-body simulations we show that the improved formula for the Coulomb logarithm fits the orbital decay very well for circular and eccentric orbits. The local scale-length of the background density distribution serves as the maximum impact parameter for a wide range of power-law indices of -1 ... -5. For each type of code the numerical resolution must be compared to the effective minimum impact parameter in order to determine the Coulomb logarithm. We also quantify the correction factors by using self-consistent velocity distribution functions instead of the standard Maxwellian often used. These factors enter directly the decay timescale and cover a range of 0.5 ... 3 for typical orbits. The new Coulomb logarithm combined with self-consistent velocity distribution functions in the Chandrasekhar formula provides a significant improvement of orbital decay times with correction up to one order of magnitude compared to the standard case. We suggest the general use of the improved formula in parameter studies as well as in special applications.

Context: A recent new approach to apply a simple dynamical mass estimate of tidally limited star clusters is based on the identification of the tidal radius in a King profile with the dynamical Jacobi radius. The application to an unbiased open cluster catalogue yields significantly higher cluster masses compared to the classical methods. Aims: We quantify the bias in the mass determination as function of projection direction and cluster age by analysing a simulated star cluster. Methods: We use direct $N$-body simulations of a star cluster including stellar evolution in an analytic Milky Way potential and apply a best fit to the projected number density of cluster stars. Results: We obtain significantly overestimated star cluster masses which depend strongly on the viewing direction. The overestimation is typically in the range of 10-50 percent and reaches a factor of 3.5 for young clusters. Mass segregation reduces the derived limiting radii systematically.

We present N-body simulations of dissolving star clusters close to galactic centres. For this purpose, we developed a new N-body program called nbody6gc based on Aarseth's series of N-body codes. We describe the algorithm in detail. We report about the density wave phenomenon in the tidal arms which has been recently explained by Kuepper et al. (2008). Standing waves develop in the tidal arms. The wave knots or clumps develop at the position, where the emerging tidal arm hits the potential wall of the effective potential and is reflected. The escaping stars move through the wave knots further into the tidal arms. We show the consistency of the positions of the wave knots with the theory in Just et al. (2009). We also demonstrate a simple method to study the properties of tidal arms. By solving many eigenvalue problems along the tidal arms, we construct numerically a 1D coordinate system whose direction is always along a principal axis of the local tensor of inertia. Along this coordinate system, physical quantities can be evaluated. The half-mass or dissolution times of our models are almost independent of the particle number which indicates that two-body relaxation is not the dominant mechanism leading to the dissolution. This may be a typical situation for many young star clusters. We propose a classification scheme which sheds light on the dissolution mechanism.

Tidal tails of star clusters are not homogeneous but show well defined clumps in observations as well as in numerical simulations. Recently an epicyclic theory for the formation of these clumps was presented. A quantitative analysis was still missing. We present a quantitative derivation of the angular momentum and energy distribution of escaping stars from a star cluster in the tidal field of the Milky Way and derive the connection to the position and width of the clumps. For the numerical realization we use star-by-star $N$-body simulations. We find a very good agreement of theory and models. We show that the radial offset of the tidal arms scales with the tidal radius, which is a function of cluster mass and the rotation curve at the cluster orbit. The mean radial offset is 2.77 times the tidal radius in the outer disc. Near the Galactic centre the circumstances are more complicated, but to lowest order the theory still applies. We have also measured the Jacobi energy distribution of bound stars and showed that there is a large fraction of stars (about 35%) above the critical Jacobi energy at all times, which can potentially leave the cluster. This is a hint that the mass loss is dominated by a self-regulating process of increasing Jacobi energy due to the weakening of the potential well of the star cluster, which is induced by the mass loss itself.

The dissolution process of star clusters is rather intricate for theory. We investigate it in the context of chaotic dynamics. We use the simple Plummer model for the gravitational field of a star cluster and treat the tidal field of the Galaxy within the tidal approximation. That is, a linear approximation of tidal forces from the Galaxy based on epicyclic theory in a rotating reference frame. The Poincar\'e surfaces of section reveal the effect of a Coriolis asymmetry. The system is non-hyperbolic which has important consequences for the dynamics. We calculated the basins of escape with respect to the Lagrangian points $L_1$ and $L_2$. The longest escape times have been measured for initial conditions in the vicinity of the fractal basin boundaries. Furthermore, we computed the chaotic saddle for the system and its stable and unstable manifolds. The chaotic saddle is a fractal structure in phase space which has the form of a Cantor set and introduces chaos into the system.

We have studied the dynamical evolution of rotating globular clusters with direct $N$-body models. Our initial models are rotating King models; we obtained results for both equal-mass systems and systems composed out of two mass components. Previous investigations using a Fokker-Planck solver have revealed that rotation has a noticeable influence on stellar systems like globular clusters, which evolve by two-body relaxation. In particular, it accelerates their dynamical evolution through the gravogyro instability. We have validated the occurence of the gravogyro instability with direct $N$-body models. In the case of systems composed out of two mass components, mass segregation takes place, which competes with the rotation in the acceleration of the core collapse. The "accelerating" effect of rotation has not been detected in our isolated two-mass $N$-body models. Last, but not least, we have looked at rotating $N$-body models in a tidal field within the tidal approximation. It turns out that rotation increases the escape rate significantly. A difference between retrograde and prograde rotating star clusters occurs with respect to the orbit of the star cluster around the Galaxy, which is due to the presence of a ``third integral'' and chaotic scattering, respectively.

We present new models of the evolution and dissolution of star clusters evolving under the combined influence of internal relaxation and external tidal fields, using the anisotropic gaseous model based on the Fokker-Planck approximation, and a new escaper loss cone model. This model borrows ideas from loss cones of stellar distributions near massive black holes, and describes physical processes related to escaping stars by a simple model based on two timescales and a diffusion process. We compare our results with those of direct $N$-body models and of direct numerical solutions of the orbit-averaged Fokker-Planck equation. For this comparative study we limit ourselves to idealized single point mass star clusters, in order to present a detailed study of the physical processes determining the rate of mass loss, core collapse and other features of the system's evolution. With the positive results of our study the path is now open in the future to use the computationally efficient gaseous models for future studies with more realism (mass spectrum, stellar evolution).