Context: Using the recent observational data of Roser 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.
(abridged) In this paper we derive observed and modelled shape parameters (apparent ellipticity and orientation of the ellipse) of 650 Galactic open clusters identified in the ASCC-2.5 catalogue. We provide the observed shape parameters of Galactic o
pen clusters, computed with the help of a multi-component analysis. For the vast majority of clusters these parameters are determined for the first time. High resolution (star by star) N-body simulations are carried out with the specially developed $phi$GRAPE code providing models of clusters of different initial masses, Galactocentric distances and rotation velocities. The comparison of models and observations of about 150 clusters reveals ellipticities of observed clusters which are too low (0.2 vs. 0.3), and offers the basis to find the main reason for this discrepancy. The models predict that after $approx 50$ Myr clusters reach an oblate shape with an axes ratio of $1.65:1.35:1$, and with the major axis tilted by an angle of $q_{XY} approx 30^circ$ with respect to the Galactocentric radius due to differential rotation of the Galaxy. Unbiased estimates of cluster shape parameters require reliable membership determination in large cluster areas up to 2-3 tidal radii where the density of cluster stars is considerably lower than the background. Although dynamically bound stars outside the tidal radius contribute insignificantly to the cluster mass, their distribution is essential for a correct determination of cluster shape parameters. In contrast, a restricted mass range of cluster stars does not play such a dramatic role, though deep surveys allow to identify more cluster members and, therefore, to increase the accuracy of the observed shape parameters.
SUPERBOX is known as a very efficient particle-mesh code with highly-resolving sub-grids. Nevertheless, the height of a typical galactic disc is small compared to the size of the whole system. Consequently, the numerical resolution in z-direction, i.
e. vertically with respect to the plane of the disc, remains poor. Here, we present a new version of SUPERBOX that allows for a considerably higher resolution along z. The improved code is applied to investigate disc heating by the infall of a galaxy satellite. We describe the improvement and communicate our results. As an important application we discuss the disruption of a dwarf galaxy within a disc-bulge-halo galaxy that consists of some 10^6 particles.
We present some preliminary results from recent numerical simulations that model the evolution of super-massive black hole (SMBH) binaries in galactic nuclei. Including the post-Newtonian terms for the binary system and adopting appropriate models fo
r the galaxies allows us, for the first time, to follow the evolution of SMBH binaries from kpc scales down to the coalescence phase. We use our results to make predictions of the detectability of such events with the gravitational wave detector LISA.
