Normalized to: Ernst, A.

**11** article(s) in total. **13** co-authors, from **1** to **9** common article(s). Median position in authors list is **1,0**.

[1] oai:arXiv.org:1506.02862 [pdf] - 1270658

Roche volume filling and the dissolution of open star clusters

Comments: 14 pages, 10+1 figures, accepted by Astronomische Nachrichten

Submitted: **2015-06-09**

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.

[2] oai:arXiv.org:1407.0509 [pdf] - 1270594

Fractal basins of escape and the formation of spiral arms in a galactic
potential with a bar

Comments: 12 pages, 10 figures, accepted by MNRAS

Submitted: **2014-07-02**

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.

[3] oai:arXiv.org:1110.1274 [pdf] - 1084682

Simulations of the Hyades

Comments: 14 pages, 16+1 figures, hopefully final version, contains a note
added in proof

Submitted: **2011-10-06**, last modified: **2013-09-29**

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.

[4] oai:arXiv.org:1212.2021 [pdf] - 1158390

Roche volume filling of star clusters in the Milky Way

Comments: 8 pages, 5+1 figures, accepted by MNRAS

Submitted: **2012-12-10**

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.

[5] oai:arXiv.org:1009.2455 [pdf] - 298418

Dynamical friction of massive objects in galactic centres

Comments: 22 pages, 28 figures, accepted by MNRAS

Submitted: **2010-09-13**

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.

[6] oai:arXiv.org:1009.0710 [pdf] - 1034687

Calibration of radii and masses of open clusters with a simulation

Comments: 9 pages, 10+1 figures, accepted by Astronomy and Astrophysics

Submitted: **2010-09-03**

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.

[7] oai:arXiv.org:0906.4459 [pdf] - 1002674

On the dissolution of star clusters in the Galactic centre. I. Circular
orbits

Comments: 18 pages, 20 figures; accepted by MNRAS

Submitted: **2009-06-24**

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.

[8] oai:arXiv.org:0808.3293 [pdf] - 315035

Quantitative analysis of clumps in the tidal tails of star clusters

Comments: 14 pages, 17 figures; accepted by MNRAS

Submitted: **2008-08-25**, last modified: **2008-10-16**

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.

[9] oai:arXiv.org:0710.4485 [pdf] - 314897

Escape from the vicinity of fractal basin boundaries of a star cluster

Comments: Accepted by MNRAS, Figures have lower quality

Submitted: **2007-10-24**

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.

[10] oai:arXiv.org:astro-ph/0702206 [pdf] - 89260

N-body Models of Rotating Globular Clusters

Comments: 16 pages, 17 figures, accepted by MNRAS

Submitted: **2007-02-07**, last modified: **2007-03-26**

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.

[11] oai:arXiv.org:astro-ph/0412698 [pdf] - 70102

Anisotropic gaseous models of tidally limited star clusters --
comparison with other methods

Comments: 13 pages, 14 figures, added two figures, updated conclusion and
discussion, accepted for publication in MNRAS

Submitted: **2004-12-30**, last modified: **2005-10-28**

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).