Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body simulations of such massive clusters using three different direct-summation N-body codes, exploring different Galactic orbits and particle numbers. By inspecting the rate at which a star's energy changes as it becomes energetically unbound from the cluster, we can neatly identify two populations we call kicks and sweeps, that escape through two-body encounters internal to the cluster and the external tidal field, respectively. We find that for a typical halo globular cluster on a moderately eccentric orbit, sweeps are far more common than kicks but the total mass loss rate is so low that these clusters can survive for hundreds of Hubble times. The different N-body codes give largely consistent results, but we find that numerical artifacts may arise in relation to the time-step parameter of the Hermite integration scheme, namely that the value required for convergent results is sensitive to the number of particles.

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster.

Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt's evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals' gravitational responses to analogues of the four solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of $\gtrsim 100$ km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close ($\ll 10^3$ au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs.

The hierarchical galaxy formation model predicts supermassive black hole binaries (SMBHBs) in galactic nuclei. Due to the gas poor environment and the limited spatial resolution in observations they may hide in the center of many a galaxy. However, a close encounter of a star with one of the supermassive black holes (SMBHs) may tidally disrupt it to produce a tidal disruption event (TDE) and temporarily light up the SMBH. In a previous work, we investigated with direct N-BODY simulations the evolution of TDE rates of SMBHB systems in galaxy mergers of equal mass. In this work we extend to unequal mass mergers. Our results show that, when two SMBHs are far away from each other, the TDE rate of each host galaxy is similar as in an isolated galaxy. As the two galaxies and their SMBHs separation shrinks, the TDE rate is increasing gradually until it reaches a maximum shortly after the two SMBHs become bound. In this stage, the averaged TDE rate can be enhanced by several times to an order of magnitude relative to isolated single galaxies. Our simulations show that the dependence of the TDE accretion rate on the mass ratio in this stage can be well fitted by power law relations for both SMBHs. After the bound SMBHB forms, the TDE rate decreases with its further evolution. We also find that in minor mergers TDEs of the secondary SMBH during and after the bound binary formation are mainly contributed by stars from the other galaxy.

As of August 2019, among the more than 4000 confirmed exoplanets, only one has been detected in a globular cluster (GC) M4. The scarce of exoplanet detections motivates us to employ direct $N$-body simulations to investigate the dynamical stability of planets in young massive clusters (YMCs), which are potentially the progenitors of GCs. In an $N=128{\rm k}$ cluster of virial radius 1.7 pc (comparable to Westerlund-1), our simulations show that most wide-orbit planets ($a\geq 20$~au) will be ejected within a timescale of 10 Myr. Interestingly, more than $70\%$ of planets with $a<5$~au survive in the 100 Myr simulations. Ignoring planet-planet scattering and tidal damping, the survivability at $t$ Myr as a function of initial semi-major axis $a_0$ in au in such a YMC can be described as $f_{\rm surv}(a_0, t)=-0.33 \log_{10}(a_0) \left(1 - e^{-0.0482t} \right) + 1$. Upon ejection, about $28.8\%$ of free-floating planets (FFPs) have sufficient speeds to escape from the host cluster at a crossing timescale. The other FFPs will remain bound to the cluster potential, but the subsequent dynamical evolution of the stellar system can result in the delayed ejection of FFPs from the host cluster. Although a full investigation of planet population in GCs requires extending the simulations to multi-Gyr, our results suggest that wide-orbit planets and free-floating planets are unlikely to be found in GCs.

Young stars are mostly found in dense stellar environments, and even our own Solar system may have formed in a star cluster. Here, we numerically explore the evolution of planetary systems similar to our own Solar system in star clusters. We investigate the evolution of planetary systems in star clusters. Most stellar encounters are tidal, hyperbolic, and adiabatic. A small fraction of the planetary systems escape from the star cluster within 50 Myr; those with low escape speeds often remain intact during and after the escape process. While most planetary systems inside the star cluster remain intact, a subset is strongly perturbed during the first 50 Myr. Over the course of time, 0.3 % - 5.3 % of the planets escape, sometimes up to tens of millions of years after a stellar encounter occurred. Survival rates are highest for Jupiter, while Uranus and Neptune have the highest escape rates. Unless directly affected by a stellar encounter itself, Jupiter frequently serves as a barrier that protects the terrestrial planets from perturbations in the outer planetary system. In low-density environments, Jupiter provides protection from perturbations in the outer planetary system, while in high-density environments, direct perturbations of Jupiter by neighbouring stars is disruptive to habitable-zone planets. The diversity amongst planetary systems that is present in the star clusters at 50 Myr, and amongst the escaping planetary systems, is high, which contributes to explaining the high diversity of observed exoplanet systems in star clusters and in the Galactic field

We present some preliminary results of our ongoing project about planetary systems around S-stars in the vicinity of Sgr A* black hole. Since S-stars might have migrated in the Galactic Centre (GC) from elsewhere, they probably still keep their planetary systems throughout their voyage. In this work, we study the destiny of their putative planetary systems after close interaction with the central black hole of our galaxy.

We investigate the dissolution process for dynamically evolving star clusters embedded in an external tidal field by exploring the MOCCA Survey Database I, with focus on the presence and evolution of a stellar-mass black hole subsystem. We argue that the presence of a black hole subsystem can lead to the dissolution of tidally filling star clusters and this can be regarded as a third type of cluster dissolution mechanism (in addition to well-known mechanisms connected with strong mass loss due to stellar evolution and mass loss connected with the relaxation process). This third process is characterized by abrupt cluster dissolution connected with the loss of dynamical equilibrium. The abrupt dissolution is powered by strong energy generation from a stellar-mass black hole subsystem accompanied by tidal stripping. Additionally, we argue that such a mechanism should also work for even tidally under-filling clusters with top-heavy initial mass function. Observationally, star clusters which undergo dissolution powered by the third mechanism would look as a 'dark cluster' i.e. composed of stellar mass black holes surrounded by an expanding halo of luminous stars (Banerjee & Kroupa 2011), and they should be different from 'dark clusters' harbouring intermediate mass black holes as discussed by Askar et al. (2017a). An additional observational consequence of an operation of the third dissolution mechanism should be a larger than expected abundance of free floating black holes in the Galactic halo.

We investigate the dissolution process of star clusters embedded in an external tidal field and harboring a subsystem of stellar-mass black hole. For this purpose we analyzed the MOCCA models of real star clusters contained in the Mocca Survey Database I. We showed that the presence of a stellar-mass black hole subsystem in tidally filling star cluster can lead to abrupt cluster dissolution connected with the loss of cluster dynamical equilibrium. Such cluster dissolution can be regarded as a third type of cluster dissolution mechanism. We additionally argue that such a mechanism should also work for tidally under-filling clusters with a top-heavy initial mass function.

We study the dynamics and evolution of the Milky Way nuclear star cluster performing a high resolution direct one-million-body simulation. Focusing on the interactions between such stellar system and the central supermassive black hole, we find that different stellar components adapt their overall distribution differently. After 5 Gyr, stellar mass black holes are characterized by a spatial distribution with power-slope $-1.75$, fully consistent with the prediction of Bahcall-Wolf pioneering work. Using the vast amount of data available, we infer the rate for tidal disruption events, being $4 \times 10^{-6}$ per yr, and estimate the number of objects that emit gravitational waves during the phases preceding the accretion onto the super-massive black hole, $\sim 270$ per Gyr. We show that some of these sources could form extreme mass-ratio inspirals. We follow the evolution of binary stars population, showing that the initial binary fraction of $5\%$ drops down to $2.5\%$ inside the inner parsec. Also, we explored the possible formation of binary systems containing a compact object, discussing the implications for millisecond pulsars formation and the development of Ia Supernovae.

A large fraction of galactic nuclei is expected to host supermassive black hole binaries (BHB), likely formed during the early phase of galaxies assembly and merging. In this paper, we use a large set of state-of-art numerical models to investigate the interplay between a BHB and a massive star cluster (GCs) driven toward the galactic centre by dynamical friction. Varying the BHB mass and mass ratio and the GC orbit, we show that the reciprocal feedback exerted between GCs and the BHB shapes their global properties. We show that, at GC-to-BHB mass ratios above 0.1, the GC affects notably the BHB orbital evolution, possibly boosting its coalescence. This effect is maximized if the GC moves on a retrograde orbit, and for a non-equal mass BHB. We show that the GC debris dispersed around the galactic nucleus can lead to the formation of a nuclear cluster, depending on the BHB tidal field, and that the distribution of compact remnants resulting from the GC disruption can carry information about the BHB orbital properties. We find that red giant stars delivered by the spiralling GC can be disrupted at a rate of $\simeq (0.7-7)\times 10^{-7}$ yr$^{-1}$ for BHB masses $\sim 10^7{\rm M}_\odot$, while tens to hundreds of stars can be possibly observed in the galactic halo as high-velocity stars, with velocities up to $\sim 2000$ km s$^{-1}$, depending on the BHB orbital properties.

Most stars, perhaps even all stars, form in crowded stellar environments. Such star forming regions typically dissolve within ten million years, while others remain bound as stellar groupings for hundreds of millions to billions of years, and then become the open clusters or globular clusters that are present in our Milky Way galaxy today. A large fraction of stars in the Galaxy hosts planetary companions. To understand the origin and dynamical evolution of such exoplanet systems, it is necessary to carefully study the effect of their environments. Here, we combine theoretical estimates with state-of-the-art numerical simulations of evolving planetary systems similar to our own solar system in different star cluster environments. We combine the REBOUND planetary system evolution code, and the NBODY6++GPU star cluster evolution code, integrated in the AMUSE multi-physics environment. With our study we can constrain the effect of external perturbations of different environments on the planets and debris structures of a wide variety of planetary systems, which may play a key role for the habitability of exoplanets in the Universe.

The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10^-2 yr^-1 from binary stars, 5 x 10^-5 yr^-1 from globular clusters and 10^-5 yr^-1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc^3). These models are thus in tension with the detection of GW170817 with an observed rate 1.5 yr^-1 (per 100 Mpc^3; LIGO/Virgo estimate). Our results imply that either (i) the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; or that (ii) physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that (iii) another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play.

We study the rates of tidal disruption of stars by intermediate-mass to supermassive black holes on bound to unbound orbits by using high-accuracy direct N-body experiments. The approaching stars from the star cluster to the black hole can take three types of orbit: eccentric, parabolic, and hyperbolic orbits. Since the mass fallback rate shows a different variability depending on these orbital types, we can classify tidal disruption events (TDEs) into three main categories: eccentric, parabolic, and hyperbolic TDEs. Respective TDEs are characterized by two critical values of the orbital eccentricity: the lower critical eccentricity is the one below which the stars on eccentric orbits cause the finite, intense accretion, and the higher critical eccentricity above which the stars on hyperbolic orbits cause no accretion. Moreover, we find that the parabolic TDEs are divided into three subclasses: precisely parabolic, marginally eccentric, and marginally hyperbolic TDEs. We analytically derive that the mass fallback rate of the marginally eccentric TDEs can be flatter and slightly higher than the standard fallback rate proportional to $t^{-5/3}$, whereas it can be flatter and lower for the marginally hyperbolic TDEs. We confirm by N-body experiments that only few eccentric, precisely parabolic, and hyperbolic TDEs can occur in a spherical stellar system with a single intermediate-mass to supermassive black hole. A substantial fraction of the stars approaching to the black hole would cause the marginally eccentric or marginally hyperbolic TDEs.

We perform high resolution direct $N$-body simulations to study the effect of an accretion disc on stellar dynamics in an active galactic nucleus (AGN). We show that the interaction of the nuclear stellar cluster (NSC) with the gaseous disc (AD) leads to formation of a stellar disc in the central part of the NSC. The accretion of stars from the stellar disc onto the super-massive black hole is balanced by the capture of stars from the NSC into the stellar disc, yielding a stationary density profile. We derive the migration time through the AD to be 3\% of the half-mass relaxation time of the NSC. The mass and size of the stellar disc are 0.7\% of the mass and 5\% of the influence radius of the super-massive black hole. An AD lifetime shorter than the migration time would result in a less massive nuclear stellar disc. The detection of such a stellar disc could point to past activity of the hosting galactic nucleus.

Numerical simulations were carried out to study the origin of multiple stellar populations in the intermediate-age clusters NGC 411 and NGC 1806 in the Magellanic Clouds. We performed NBODY6++ simulations based on two different formation scenarios, an ad hoc formation model where second-generation (SG) stars are formed inside a cluster of first-generation (FG) stars using the gas accumulated from the external intergalactic medium and a minor merger model of unequal mass ($M_{SG}$/$M_{FG}$ ~5-10%) clusters with an age difference of a few hundred million years. We compared our results such as the radial profile of the SG-to-FG number ratio with observations on the assumption that the SG stars in the observations are composed of cluster members, and confirmed that both the ad hoc formation and merger scenarios reproduce the observed radial trend of the SG-to-FG number ratio which shows less centrally concentrated SG than FG stars. It is difficult to constrain the formation scenario for the multiple populations by only using the spatial distribution of the SG stars. SG stars originating from the merger scenario show a significant velocity anisotropy and rotational features compared to those from the ad hoc formation scenario. Thus, observations aimed at kinematic properties like velocity anisotropy or rotational velocities for SG stars should be obtained to better understand the formation of the multiple populations in these clusters. This is, however, beyond current instrumentation capabilities.

Most stars form in star clusters and stellar associated. To understand the roles of star cluster environments in shaping the dynamical evolution of planetary systems, we carry out direct $N$-body simulations of four planetary systems models in three different star cluster environments with respectively N=2k, 8k and 32k stars. In each cluster, an ensemble of initially identical planetary systems are assigned to solar-type stars with $\sim 1 M_{\odot}$ and evolved for 50~Myr. We found that following the depletion of protoplanetary disks, external perturbations and planet-planet interactions are two driving mechanisms responsible for the destabilization of planetary systems. The planet survival rate varies from $\sim 95\%$ in the N=2k cluster to $\sim 60\%$ in the N=32k cluster, which suggests that most planetary systems can indeed survive in low-mass clusters, except in the central regions. We also find that planet ejections through stellar encounters are cumulative processes, as only $\sim 3\%$ of encounters are strong enough to excite the eccentricity by $\Delta e \geq 0.5$. Short-period planets can be perturbed through orbit crossings with long-period planets. When taking into account planet-planet interactions, the planet ejection rate nearly doubles, and therefore multiplicity contributes to the vulnerability of planetary systems. In each ensemble, $\sim 0.2\%$ of planetary orbits become retrograde due to random directions of stellar encounters. Our results predict that young low-mass star clusters are promising sites for next-generation planet surveys, yet low planet detection rates are expected in dense globular clusters such as 47 Tuc. Nevertheless, planets in denser stellar environments are likely to have shorter orbital periods, which enhances their detectability.

Supermassive black hole binaries (SMBHBs) are productions of the hierarchical galaxy formation model. There are many close connections between central SMBH and its host galaxy because the former plays very important roles on the formation and evolution of a galaxy. For this reason, the evolution of SMBHBs in merging galaxies is an essential problem. Since there are many discussions about SMBHB evolution in gas rich environment, we focus on the quiescent galaxy, using tidal disruption as a diagnostic tool. Our study is based on a series of numerical large particle number direct N-body simulations for dry major mergers. According to the simulation results, the evolution can be divided into three phases. In phase I, the tidal disruption rate for two well separated SMBHs in merging system has similar level to single SMBH in isolate galaxy. After two SMBHs getting close enough to form a bound binary in phase II, the disruption rate can be enhanced for ~ 2 order of magnitudes within a short time. This "boosted" disruption stage finishes after the SMBHB evolving to compact binary system in phase III, corresponding to a drop back of disruption rate to a level of a few times higher than that in Phase I. How to correctly extrapolate our N-body simulation results to reality, and implications of our results to observations, are discussed too.

Thousands of confirmed and candidate exoplanets have been identified in recent years. Consequently, theoretical research on the formation and dynamical evolution of planetary systems has seen a boost, and the processes of planet-planet scattering, secular evolution, and interaction between planets and gas/debris disks have been well-studied. Almost all of this work has focused on the formation and evolution of isolated planetary systems, and neglect the effect of external influences, such as the gravitational interaction with neighbouring stars. Most stars, however, form in clustered environments that either quickly disperse, or evolve into open clusters. Under these conditions, young planetary systems experience frequent close encounters with other stars, at least during the first 1-10 Myr, which affects planets orbiting at any period range, as well as their debris structures.

We examine the effect of an accretion disc on the orbits of stars in the central star cluster surrounding a central massive black hole by performing a suite of 39 high-accuracy direct N-body simulations using state-of-the art software and accelerator hardware, with particle numbers up to 128k. The primary focus is on the accretion rate of stars by the black hole (equivalent to their tidal disruption rate for black holes in the small to medium mass range) and the eccentricity distribution of these stars. Our simulations vary not only the particle number, but disc model (two models examined), spatial resolution at the centre (characterised by the numerical accretion radius) and softening length. The large parameter range and physically realistic modelling allow us for the first time to confidently extrapolate these results to real galactic centres. While in a real galactic centre both particle number and accretion radius differ by a few orders of magnitude from our models, which are constrained by numerical capability, we find that the stellar accretion rate converges for models with N > 32k. The eccentricity distribution of accreted stars, however, does not converge. We find that there are two competing effects at work when improving the resolution: larger particle number leads to a smaller fraction of stars accreted on nearly-circular orbits, while higher spatial resolution increases this fraction. We scale our simulations to some nearby galaxies and find that the expected boost in stellar accretion (or tidal disruption, which could be observed as X-ray flares) in the presence of a gas disc is about a factor of 10. Even with this boost, the accretion of mass from stars is still a factor of ~ 100 slower than the accretion of gas from the disc. Thus, it seems accretion of stars is not a major contributor to black hole mass growth.

Introducing the DRAGON simulation project, we present direct $N$-body simulations of four massive globular clusters (GCs) with $10^6$ stars and 5$\%$ primordial binaries at a high level of accuracy and realism. The GC evolution is computed with NBODY6++GPU and follows the dynamical and stellar evolution of individual stars and binaries, kicks of neutron stars and black holes, and the effect of a tidal field. We investigate the evolution of the luminous (stellar) and dark (faint stars and stellar remnants) GC components and create mock observations of the simulations (i.e. photometry, color-magnitude diagrams, surface brightness and velocity dispersion profiles). By connecting internal processes to observable features we highlight the formation of a long-lived 'dark' nuclear subsystem made of black holes (BHs), which results in a two-component structure. The inner core is dominated by the BH subsystem and experiences a core collapse phase within the first Gyr. It can be detected in the stellar (luminous) line-of-sight velocity dispersion profiles. The outer extended core - commonly observed in the (luminous) surface brightness profiles - shows no collapse features and is continuously expanding. We demonstrate how a King (1966) model fit to observed clusters might help identify the presence of post core-collapse BH subsystems. For global observables like core and half-mass radii the direct simulations agree well with Monte-Carlo models. Variations in the initial mass function can result in significantly different GC properties (e.g. density distributions) driven by varying amounts of early mass loss and the number of forming BHs.

We present the first detailed comparison between million-body globular cluster simulations computed with a H\'enon-type Monte Carlo code, CMC, and a direct $N$-body code, NBODY6++GPU. Both simulations start from an identical cluster model with $10^6$ particles, and include all of the relevant physics needed to treat the system in a highly realistic way. With the two codes "frozen" (no fine-tuning of any free parameters or internal algorithms of the codes) we find excellent agreement in the overall evolution of the two models. Furthermore, we find that in both models, large numbers of stellar-mass black holes (> 1000) are retained for 12 Gyr. Thus, the very accurate direct $N$-body approach confirms recent predictions that black holes can be retained in present-day, old globular clusters. We find only minor disagreements between the two models and attribute these to the small-$N$ dynamics driving the evolution of the cluster core for which the Monte Carlo assumptions are less ideal. Based on the overwhelming general agreement between the two models computed using these vastly different techniques, we conclude that our Monte Carlo approach, which is more approximate, but dramatically faster compared to the direct $N$-body, is capable of producing a very accurate description of the long-term evolution of massive globular clusters even when the clusters contain large populations of stellar-mass black holes.

We present $GalevNB$ (Galev for $N$-body simulations), an utility that converts fundamental stellar properties of $N$-body simulations into observational properties using the $GALEV$ (GAlaxy EVolutionary synthesis models) package, and thus allowing direct comparisons between observations and $N$-body simulations. It works by converting fundamental stellar properties, such as stellar mass, temperature, luminosity and metallicity into observational magnitudes for a variety of filters of mainstream instruments/telescopes, such as HST, ESO, SDSS, 2MASS, etc.), and into spectra that spans from far-UV (90 $\rm \AA$) to near-IR (160 $\rm \mu$m). As an application, we use $GalevNB$ to investigate the secular evolution of spectral energy distribution (SED) and color-magnitude diagram (CMD) of a simulated star cluster over a few hundred million years. With the results given by $GalevNB$ we discover an UV-excess in the SED of the cluster over the whole simulation time. We also identify four candidates that contribute to the FUV peak, core helium burning stars, thermal pulsing asymptotic giant branch (TPAGB) stars, white dwarfs and naked helium stars.

Tidal Disruption of stars by supermassive central black holes from dense rotating star clusters is modelled by high-accuracy direct N-body simulation. As in a previous paper on spherical star clusters we study the time evolution of the stellar tidal disruption rate and the origin of tidally disrupted stars, now according to several classes of orbits which only occur in axisymmetric systems (short axis tube and saucer). Compared with that in spherical systems, we found a higher TD rate in axisymmetric systems. The enhancement can be explained by an enlarged loss-cone in phase space which is raised from the fact that total angular momentum $\bf J$ is not conserved. As in the case of spherical systems, the distribution of the last apocenter distance of tidally accreted stars peaks at the classical critical radius. However, the angular distribution of the origin of the accreted stars reveals interesting features. Inside the influence radius of the supermassive black hole the angular distribution of disrupted stars has a conspicuous bimodal structure with a local minimum near the equatorial plane. Outside the influence radius this dependence is weak. We show that the bimodal structure of orbital parameters can be explained by the presence of two families of regular orbits, namely short axis tube and saucer orbits. Also the consequences of our results for the loss cone in axisymmetric galactic nuclei are presented.

We compare the performance of two very different parallel gravitational $N$-body codes for astrophysical simulations on large GPU clusters, both pioneer in their own fields as well as in certain mutual scales - NBODY6++ and Bonsai. We carry out the benchmark of the two codes by analyzing their performance, accuracy and efficiency through the modeling of structure decomposition and timing measurements. We find that both codes are heavily optimized to leverage the computational potential of GPUs as their performance has approached half of the maximum single precision performance of the underlying GPU cards. With such performance we predict that a speed-up of $200-300$ can be achieved when up to 1k processors and GPUs are employed simultaneously. We discuss the quantitative information about comparisons of two codes, finding that in the same cases Bonsai adopts larger time steps as well as relative energy errors than NBODY6++, typically ranging from $10-50$ times larger, depending on the chosen parameters of the codes. While the two codes are built for different astrophysical applications, in specified conditions they may overlap in performance at certain physical scale, and thus allowing the user to choose from either one with finetuned parameters accordingly.

We present a hybrid code combining the OpenMP-parallel tree code VINE with an algorithmic chain regularization scheme. The new code, called "rVINE", aims to significantly improve the accuracy of close encounters of massive bodies with supermassive black holes in galaxy-scale numerical simulations. We demonstrate the capabilities of the code by studying two test problems, the sinking of a single massive black hole to the centre of a gas-free galaxy due to dynamical friction and the hardening of a supermassive black hole binary due to close stellar encounters. We show that results obtained with rVINE compare well with NBODY7 for problems with particle numbers that can be simulated with NBODY7. In particular, in both NBODY7 and rVINE we find a clear N-dependence of the binary hardening rate, a low binary eccentricity and moderate eccentricity evolution, as well as the conversion of the galaxy's inner density profile from a cusp to a a core via the ejection of stars at high velocity. The much larger number of particles that can be handled by rVINE will open up exciting opportunities to model stellar dynamics close to SMBHs much more accurately in a realistic galactic context. This will help to remedy the inherent limitations of commonly used tree solvers to follow the correct dynamical evolution of black holes in galaxy scale simulations.

Astrophysical research in recent decades has made significant progress thanks to the availability of various $N$-body simulation techniques. With the rapid development of high-performance computing technologies, modern simulations have been able to take the computing power of massively parallel clusters with more than $10^5$ GPU cores. While unprecedented accuracy and dynamical scales have been achieved, the enormous amount of data being generated continuously poses great challenges for the subsequent procedures of data analysis and archiving. As an urgent response to these challenges, in this paper we propose an adaptive storage scheme for simulation data, inspired by the block time step integration scheme found in a number of direct $N$-body integrators available nowadays. The proposed scheme, namely the block time step storage scheme, works by minimizing the data redundancy with assignments of data with individual output frequencies as required by the researcher. As demonstrated by benchmarks, the proposed scheme is applicable to a wide variety of simulations. Despite the main focus of developing a solution for direct $N$-body simulation data, the methodology is transferable for grid-based or tree-based simulations where hierarchical time stepping is used.

Accurate direct $N$-body simulations help to obtain detailed information about the dynamical evolution of star clusters. They also enable comparisons with analytical models and Fokker-Planck or Monte-Carlo methods. NBODY6 is a well-known direct $N$-body code for star clusters, and NBODY6++ is the extended version designed for large particle number simulations by supercomputers. We present NBODY6++GPU, an optimized version of NBODY6++ with hybrid parallelization methods (MPI, GPU, OpenMP, and AVX/SSE) to accelerate large direct $N$-body simulations, and in particular to solve the million-body problem. We discuss the new features of the NBODY6++GPU code, benchmarks, as well as the first results from a simulation of a realistic globular cluster initially containing a million particles. For million-body simulations, NBODY6++GPU is $400-2000$ times faster than NBODY6 with 320 CPU cores and 32 NVIDIA K20X GPUs. With this computing cluster specification, the simulations of million-body globular clusters including $5\%$ primordial binaries require about an hour per half-mass crossing time.

We report on GalevNB (Galev for N-Body simulations), an integrated software solution that provides N-body users direct access to the software package GALEV (GALaxy EVolutionary synthesis models). GalevNB is developed for the purpose of a direct comparison between N-body simulations and observations. It converts the fundamental stellar properties of N-body simulations, i.e., stellar mass, temperature, stellar luminosity and metallicity, into observational magnitudes for a variety of filters of widely used instruments/telescopes (HST, ESO, SDSS, 2MASS), and into spectra that span from far-UV (90 $\rm \AA$) to near-IR (160 $\rm \mu$m).

Tidal interactions of galaxies in galaxy clusters have been proposed as one potential explanation of the morphology-density relation at low masses. The efficiency of tidal transformation is expected to depend strongly on the orbit of a galaxy within the cluster halo. The orbit determines both the strength of the cluster's global tidal field and the probability of encounters with other cluster members. Additionally, we investigate the effect of an inclination between the disk of the infalling galaxy and its orbital plane. We also compare our results to observational data. We find that galaxies that entered a cluster from the outskirts are unlikely to be significantly transformed. Closer to the cluster centre, tidal interactions are a more efficient mechanism for producing harassed galaxies. But the inclination of the disk can reduce the mass loss significantly, yet it amplifies the thickening.

In the solar neighborhood, where the typical relaxation timescale is larger than the cosmic age, at least 10\% to 15\% of Sun-like stars have planetary systems with Jupiter-mass planets. In contrast, dense star clusters, charactered by frequent close encounters, have been found to host very few planets. We carry out numerical simulations with different initial conditions to investigate the dynamical stability of planetary systems in star cluster environments.

The present work reports on a feasibility study commissioned by the Chinese Academy of Sciences of China to explore various possible mission options to detect gravitational waves in space alternative to that of the eLISA/LISA mission concept. Based on the relative merits assigned to science and technological viability, a few representative mission options descoped from the ALIA mission are considered. A semi-analytic Monte Carlo simulation is carried out to understand the cosmic black hole merger histories starting from intermediate mass black holes at high redshift as well as the possible scientific merits of the mission options considered in probing the light seed black holes and their coevolution with galaxies in early Universe. The study indicates that, by choosing the armlength of the interferometer to be three million kilometers and shifting the sensitivity floor to around one-hundredth Hz, together with a very moderate improvement on the position noise budget, there are certain mission options capable of exploring light seed, intermediate mass black hole binaries at high redshift that are not readily accessible to eLISA/LISA, and yet the technological requirements seem to within reach in the next few decades for China.

Tidal Disruption of stars by super massive central black holes from dense star clusters is modeled by high-accuracy direct $N$-body simulation. The time evolution of the stellar tidal disruption rate, the effect of tidal disruption on the stellar density profile and for the first time the detailed origin of tidally disrupted stars are carefully examined and compared with classic papers in the field. Up to 128k particles are used in simulation to model the star cluster around the super massive black hole, we use the particle number and the tidal radius of black hole as free parameters for a scaling analysis. The transition from full to empty loss-cone is analyzed in our data, the tidal disruption rate scales with the particle number $N$ in the expected way for both cases. For the first time in numerical simulations (under certain conditions) we can support the concept of a critical radius of Frank & Rees (1976), which claims that most stars are tidally accreted on highly eccentric orbits originating from regions far outside the tidal radius. Due to the consumption of stars moving on radial orbits, a velocity anisotropy is founded inside the cluster. Finally we make an estimation for the real galactic center based on our simulation results and the scaling analysis.

We present GPU implementations of two fast force calculation methods, based on series expansions of the Poisson equation. One is the Self-Consistent Field (SCF) method, which is a Fourier-like expansion of the density field in some basis set; the other is the Multipole Expansion (MEX) method, which is a Taylor-like expansion of the Green's function. MEX, which has been advocated in the past, has not gained as much popularity as SCF. Both are particle-field method and optimized for collisionless galactic dynamics, but while SCF is a "pure" expansion, MEX is an expansion in just the angular part; it is thus capable of capturing radial structure easily, where SCF needs a large number of radial terms. We show that despite the expansion bias, these methods are more accurate than direct techniques for the same number of particles. The performance of our GPU code, which we call ETICS, is profiled and compared to a CPU implementation. On the tested GPU hardware, a full force calculation for one million particles took ~ 0.1 seconds (depending on expansion cutoff), making simulations with as many as $10^8$ particles fast on a comparatively small number of nodes.

Dense stellar systems such as globular clusters and dense nuclear clusters are the breeding ground of sources of gravitational waves for the advanced detectors LIGO and Virgo. These systems deserve a close study to estimate rates and parameter distribution. This is not an easy task, since the evolution of a dense stellar cluster involves the integration of $N$ bodies with high resolution in time and space and including hard binaries and their encounters and, in the case of gravitational waves (GWs), one needs to take into account important relativistic corrections. In this work we present the first implementation of the effect of spin in mergers in a direct-summation code, NBODY6. We employ non-spinning post-Newtonian corrections to the Newtonian accelerations up to 3.5 post-Newtonian (PN) order as well as the spin-orbit coupling up to next-to-lowest order and the lowest order spin-spin coupling. We integrate spin precession and add a consistent treatment of mergers. We analyse the implementation by running a set of two-body experiments and then we run a set of 500 simulations of a relativistic stellar cluster. In spite of the large number of mergers in our tests, the application of the algorithm is robust. We find in particular the formation of a runaway star whose spin decays with the mass it wins, independently of the initial value of the spins of the stars. More remarkably, the subset of compact objects that do not undergo many mergers, and hence represent a more realistic system, has a correlation between the final absolute spin and the initial choice for the initial distribution, which could provide us with information about the evolution of spins in dense clusters once the first detections have started.

We present direct astrophysical N-body simulations with up to a few million bodies using our parallel MPI/CUDA code on large GPU clusters in China, Ukraine and Germany, with different kinds of GPU hardware. These clusters are directly linked under the Chinese Academy of Sciences special GPU cluster program in the cooperation of ICCS (International Center for Computational Science). We reach about the half the peak Kepler K20 GPU performance for our phi-GPU code [2], in a real application scenario with individual hierarchically block time-steps with the high (4th, 6th and 8th) order Hermite integration schemes and a real core-halo density structure of the modeled stellar systems. The code and hardware are mainly used to simulate star clusters [23, 24] and galactic nuclei with supermassive black holes [20], in which correlations between distant particles cannot be neglected.

The hierarchical galaxy formation picture suggests that super massive black holes (MBHs) observed in galactic nuclei today have grown from coalescence of massive black hole binaries (MBHB) after galaxy merging. Once the components of a MBHB become gravitationally bound, strong three-body encounters between the MBHB and stars dominate its evolution in a "dry" gas free environment, and change the MBHB's energy and angular momentum (semi-major axis, eccentricity and orientation). Here we present high accuracy direct N-body simulations of spherical and axisymmetric (rotating) galactic nuclei with order a million stars and two massive black holes that are initially unbound. We analyze the properties of the ejected stars due to slingshot effects from three-body encounters with the MBHB in detail. Previous studies have investigated the eccentricity and energy changes of MBHs using approximate models or Monte-Carlo three body scatterings. We find general agreement with the average results of previous semi-analytic models for spherical galactic nuclei, but our results show a large statistical variation. Our new results show many more phase space details of how the process works, and also show the influence of stellar system rotation on the process. We detect that the angle between the orbital plane of the MBHBs and that of the stellar system (when it rotates) influences the phase-space properties of the ejected stars. We also find that massive MBHB tend to switch stars with counter-rotating orbits into co-rotating orbits during their interactions.

The majority of stars form in star clusters and many are thought to have planetary companions. We demonstrate that multi-planet systems are prone to instabilities as a result of frequent stellar encounters in these star clusters much more than single-planet systems. The cumulative effect of close and distant encounters on these planetary systems are investigated using Monte Carlo scattering experiments. We consider two types of planetary configurations orbiting Sun-like stars: (i) five Jupiter-mass planets in the semi-major axis range 1-42 AU orbiting a Solar mass star, with orbits that are initially co-planar, circular, and separated by 10 mutual Hill radii, and (ii) the four gas giants of our Solar system. Planets with short orbital periods are not directly affected by encountering stars. However, secular evolution of perturbed systems may result in the ejection of the innermost planets or in physical collisions of the innermost planets with the host star, up to many thousands of years after a stellar encounter. The simulations of the Solar system-like systems indicate that Saturn, Uranus and Neptune are affected by both direct interactions with encountering stars, as well as planet-planet scattering. Jupiter, on the other hand, is almost only affected by direct encounters with neighbouring stars, as its mass is too large to be substantially perturbed by the other three planets. Our results indicate that stellar encounters can account for the apparent scarcity of exoplanets in star clusters, not only for those on wide-orbit that are directly affected by stellar encounters, but also planets close to the star which can disappear long after a stellar encounter has perturbed the planetary system.

To understand the effects of the initial rotation on the evolution of the tidally limited clusters with mass spectrum, we have performed N-body simulations of the clusters with different initial rotations and compared the results with those of the Fokker-Planck (FP) simulations. We confirmed that the cluster evolution is accelerated by not only the initial rotation but also the mass spectrum. For the slowly rotating models, the time evolutions of mass, energy and angular momentum show good agreements between N-body and FP simulations. On the other hand, for the rapidly rotating models, there are significant differences between these two approaches at the early stage of the evolutions because of the development of bar instability in N-body simulations. The shape of the cluster for N-body simulations becomes tri-axial or even prolate, which cannot be produced by the 2-dimensional FP simulations. The total angular momentum and the total mass of the cluster decrease rapidly while bar-like structure persists. After the rotational energy becomes smaller than the critical value for the bar instability, the shape of the cluster becomes nearly axisymmetric again, and follows the evolutionary track predicted by the FP equation. We have confirmed again that the energy equipartiton is not completely achieved when M2/M1(m2>/m1)^(3/2) > 0.16. By examining the angular momentum at each mass component, we found that the exchange of angular momentum between different mass components occurs, similar to the energy exchange leading to the equipartition.

The formation and evolution of protoplanetary discs remains a challenge from both a theoretical and numerical standpoint. In this work we first perform a series of tests of our new hybrid algorithm presented in Glaschke, Amaro-Seoane and Spurzem 2011 (henceforth Paper I) that combines the advantages of high accuracy of direct-summation N-body methods with a statistical description for the planetesimal disc based on Fokker-Planck techniques. We then address the formation of planets, with a focus on the formation of protoplanets out of planetesimals. We find that the evolution of the system is driven by encounters as well as direct collisions and requires a careful modelling of the evolution of the velocity dispersion and the size distribution over a large range of sizes. The simulations show no termination of the protoplanetary accretion due to gap formation, since the distribution of the planetesimals is only subjected to small fluctuations. We also show that these features are weakly correlated with the positions of the protoplanets. The exploration of different impact strengths indicates that fragmentation mainly controls the overall mass loss, which is less pronounced during the early runaway growth. We prove that the fragmentation in combination with the effective removal of collisional fragments by gas drag sets an universal upper limit of the protoplanetary mass as a function of the distance to the host star, which we refer to as the mill condition.

We investigate the dynamical interaction of a central star cluster surrounding a super-massive black hole and a central accretion disk. The dissipative force acting on stars in the disk leads to an enhanced mass flow towards the super-massive black hole and to an asymmetry in the phase space distribution due to the rotating accretion disk. The accretion disk is considered as a stationary Keplerian rotating disk, which is vertically extended in order to employ a fully self-consistent treatment of stellar dynamics including the dissipative force originating from star-gas ram pressure effects. The stellar system is treated with a direct high-accuracy N-body integration code. A star-by-star representation, desirable in N-body simulations, cannot be extended to real particle numbers yet. Hence, we carefully discuss the scaling behavior of our model with regard to particle number and tidal accretion radius. The main idea is to find a family of models for which the ratio of two-body relaxation time and dissipation time (for kinetic energy of stellar orbits) is constant, which then allows us to extrapolate our results to real parameters of galactic nuclei. Our model is derived from basic physical principles and as such it provides insight into the role of physical processes in galactic nuclei, but it should be regarded as a first step towards more realistic and more comprehensive simulations. Nevertheless, the following conclusions appear to be robust: the star accretion rate onto the accretion disk and subsequently onto the super-massive black hole is enhanced by a significant factor compared to purely stellar dynamical systems neglecting the disk. This process leads to enhanced fueling of central disks in active galactic nuclei and to an enhanced rate of tidal stellar disruptions. [Abridged]

Galaxy centers are residing places for Super Massive Black Holes (SMBHs). Galaxy mergers bring SMBHs close together to form gravitationally bound binary systems which, if able to coalesce in less than a Hubble time, would be one of the most promising sources of gravitational waves for the Laser Interferometer Space Antenna (LISA). In spherical galaxy models, SMBH binaries stall at a separation of approximately one parsec, leading to the "final parsec problem" (FPP). On the other hand, it has been shown that merger-induced triaxiality of the remnant in equal-mass mergers is capable of supporting a constant supply of stars on so-called centrophilic orbits that interact with the binary and thus avoid the FPP. In this paper, using a set of direct N-body simulations of mergers of initially spherically symmetric galaxies with different mass ratios, we show that the merger-induced triaxiality is able to drive unequal-mass SMBH binaries to coalescence. The binary hardening rates are high and depend only weakly on the mass ratios of SMBHs for a wide range of mass ratios q. The hardening rates are significantly higher for galaxies having steep cusps in comparison with those having shallow cups at centers. The evolution of the binary SMBH leads to relatively shallower inner slopes at the centers of the merger remnants. The stellar mass displaced by the SMBH binary on its way to coalescence is ~ 1-5 times the combined mass of binary SMBHs. The coalescence times for SMBH binary with mass ~ million solar masses are less than 1 Gyr and for those at the upper end of SMBH masses (~ billion solar masses) are 1-2 Gyr for less eccentric binaries whereas less than 1 Gyr for highly eccentric binaries. SMBH binaries are thus expected to be promising sources of gravitational waves at low and high redshifts.

Supermassive black hole binaries (SMBHBs) are the products of frequent galaxy mergers. The coalescence of the SMBHBs is a distinct source of gravitational wave (GW) radiation. The detections of the strong GW radiation and their possible electromagnetic counterparts are essential. Numerical relativity suggests that the post-merger supermassive black hole (SMBH) gets a kick velocity up to 4000 km/s due to the anisotropic GW radiations. Here we investigate the dynamical co-evolution and interaction of the recoiling SMBHs and their galactic stellar environments with one million direct N-body simulations including the stellar tidal disruption by the recoiling SMBHs. Our results show that the accretion of disrupted stars does not significantly affect the SMBH dynamical evolution. We investigate the stellar tidal disruption rates as a function of the dynamical evolution of oscillating SMBHs in the galactic nuclei. Our simulations show that most of stellar tidal disruptions are contributed by the unbound stars and occur when the oscillating SMBHs pass through the galactic center. The averaged disruption rate is ~10^{-6} M_\odot yr^{-1}, which is about an order of magnitude lower than that by a stationary SMBH at similar galactic nuclei. Our results also show that a bound star cluster is around the oscillating SMBH of about ~ 0.7% the black hole mass. In addition, we discover a massive cloud of unbound stars following the oscillating SMBH. We also investigate the dependence of the results on the SMBH masses and density slopes of the galactic nuclei.

NBody realizations of axisymmetric collisional galaxy cores (e.g. M32, M33, NGC205, Milky Way) with embedded growing black holes are presented. Stars which approach the disruption sphere are disrupted and accreted to the black hole. We measure the zone of influence of the black hole and disruption rates in relaxation time scales. We show that secular gravitational instabilities dominate the initial core dynamics, while the black hole is small and growing due to consumption of stars. Later, the black hole potential dominates the core, and loss cone theory can be applied. Our simulations show that central rotation in galaxies can not be neglected for relaxed systems, and compare and discuss our results with the standard theory of spherically symmetric systems.

The young star clusters we observe today are the building blocks of a new generation of stars and planets in our Galaxy and beyond. Despite their fundamental role we still lack knowledge about the conditions under which star clusters form and the impact of these often harsh environments on the evolution of their stellar and substellar members. We demonstrate the vital role numerical simulations play to uncover both key issues. Using dynamical models of different star cluster environments we show the variety of effects stellar interactions potentially have. Moreover, our significantly improved measure of mass segregation reveals that it can occur rapidly even for star clusters without substructure. This finding is a critical step to resolve the controversial debate on mass segregation in young star clusters and provides strong constraints on their initial conditions.

Investigations of mass segregation are of vital interest for the understanding of the formation and dynamical evolution of stellar systems on a wide range of spatial scales. Our method is based on the minimum spanning tree (MST) that serves as a geometry-independent measure of concentration. Compared to previous such approaches we obtain a significant refinement by using the geometrical mean as an intermediate-pass. It allows the detection of mass segregation with much higher confidence and for much lower degrees of mass segregation than other approaches. The method shows in particular very clear signatures even when applied to small subsets of the entire population. We confirm with high significance strong mass segregation of the five most massive stars in the Orion Nebula Cluster (ONC). Our method is the most sensitive general measure of mass segregation so far and provides robust results for both data from simulations and observations. As such it is ideally suited for tracking mass segregation in young star clusters and to investigate the long standing paradigm of primordial mass segregation by comparison of simulations and observations.

The formation and evolution of protoplanetary systems, the breeding grounds of planet formation, is a complex dynamical problem that involves many orders of magnitudes. To serve this purpose, we present a new hybrid algorithm that combines a Fokker-Planck approach with the advantages of a pure direct-summation N-body scheme, with a very accurate integration of close encounters for the orbital evolution of the larger bodies with a statistical model, envisaged to simulate the very large number of smaller planetesimals in the disc. Direct-summation techniques have been historically developped for the study of dense stellar systems such as open and globular clusters and, within some limits imposed by the number of stars, of galactic nuclei. The number of modifications to adapt direct-summation N-body techniques to planetary dynamics is not undemanding and requires modifications. These include the way close encounters are treated, as well as the selection process for the "neighbour radius" of the particles and the extended Hermite scheme, used for the very first time in this work, as well as the implementation of a central potential, drag forces and the adjustment of the regularisation treatment. For the statistical description of the planetesimal disc we employ a Fokker-Planck approach. We include dynamical friction, high- and low-speed encounters, the role of distant encounters as well as gas and collisional damping and then generalise the model to inhomogenous discs. We then describe the combination of the two techniques to address the whole problem of planetesimal dynamics in a realistic way via a transition mass to integrate the evolution of the particles according to their masses.

We investigate a purely stellar dynamical solution to the Final Parsec Problem. Galactic nuclei resulting from major mergers are not spherical, but show some degree of triaxiality. With $N$-body simulations, we show that massive black hole binaries (MBHB) hosted by them will continuously interact with stars on centrophilic orbits and will thus inspiral---in much less than a Hubble time---down to separations at which gravitational wave (GW) emission is strong enough to drive them to coalescence. Such coalescences will be important sources of GWs for future space-borne detectors such as the {\it Laser Interferometer Space Antenna} (LISA). Based on our results, we expect that LISA will see between $\sim 10$ to $\sim {\rm few} \times 10^2$ such events every year, depending on the particular MBH seed model as obtained in recent studies of merger trees of galaxy and MBH co-evolution. Orbital eccentricities in the LISA band will be clearly distinguishable from zero with $e \gtrsim 0.001-0.01$.

We study the compact binary population in star clusters, focusing on binaries containing black holes, using a self-consistent Monte Carlo treatment of dynamics and full stellar evolution. We find that the black holes experience strong mass segregation and become centrally concentrated. In the core the black holes interact strongly with each other and black hole-black hole binaries are formed very efficiently. The strong interactions, however, also destroy or eject the black hole-black hole binaries. We find no black hole-black hole mergers within our simulations but produce many hard escapers that will merge in the galactic field within a Hubble time. We also find several highly eccentric black hole-black hole binaries that are potential LISA sources, suggesting that star clusters are interesting targets for space-based detectors. We conclude that star clusters must be taken into account when predicting compact binary population statistics.

We present a new study of the evolution of the Carina dwarf galaxy that includes a simultaneous derivation of its orbit and star formation history. The structure of the galaxy is constrained through orbital parameters derived from the observed distance, proper motions, radial velocity and star formation history. The different orbits admitted by the large proper motion errors are investigated in relation to the tidal force exerted by an external potential representing the Milky Way (MW). Our analysis is performed with the aid of fully consistent N-body simulations that are able to follow the dynamics and the stellar evolution of the dwarf system in order to determine self-consistently the star formation history of Carina. We find a star formation history characterized by several bursts, partially matching the observational expectation. We find also compatible results between dynamical projected quantities and the observational constraints. The possibility of a past interaction between Carina and the Magellanic Clouds is also separately considered and deemed unlikely.

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.

We use a self-consistent Monte Carlo treatment of stellar dynamics to investigate black hole binaries that are dynamically ejected from globular clusters to determine if they will be gravitational wave sources. We find that many of the ejected binaries have initially short periods and will merge within a Hubble time due to gravitational wave radiation. Thus they are potential sources for ground-based gravitational wave detectors. We estimate the yearly detection rate for current and advanced ground-based detectors and find a modest enhancement over the rate predicted for binaries produced by pure stellar evolution in galactic fields. We also find that many of the ejected binaries will pass through the longer wavelength Laser Interferometer Space Antenna (LISA) band and may be individually resolvable. We find a low probability that the Galaxy will contain a binary in the LISA band during its three-year mission. Some such binaries may, however, be detectable at Mpc distances implying that there may be resolvable stellar-mass LISA sources beyond our Galaxy. We conclude that globular clusters have a significant effect on the detection rate of ground-based detectors and may produce interesting LISA sources in local group galaxies.

Dense stellar systems such as globular clusters, galactic nuclei and nuclear star clusters are ideal loci to study stellar dynamics due to the very high densities reached, usually a million times higher than in the solar neighborhood; they are unique laboratories to study processes related to relaxation. There are a number of different techniques to model the global evolution of such a system. In statistical models we assume that relaxation is the result of a large number of two-body gravitational encounters with a net local effect. We present two moment models that are based on the collisional Boltzmann equation. By taking moments of the Boltzmann equation one obtains an infinite set of differential moment equations where the equation for the moment of order $n$ contains moments of order $n+1$. In our models we assume spherical symmetry but we do not require dynamical equilibrium. We truncate the infinite set of moment equations at order $n=4$ for the first model and at order $n=5$ for the second model. The collisional terms on the right-hand side of the moment equations account for two-body relaxation and are computed by means of the Rosenbluth potentials. We complete the set of moment equations with closure relations which constrain the degree of anisotropy of our model by expressing moments of order $n+1$ by moments of order $n$. The accuracy of this approach relies on the number of moments included from the infinite series. Since both models include fourth order moments we can study mechanisms in more detail that increase or decrease the number of high velocity stars. The resulting model allows us to derive a velocity distribution function, with unprecedented accuracy, compared to previous moment models.

We present status and results of AstroGrid-D, a joint effort of astrophysicists and computer scientists to employ grid technology for scientific applications. AstroGrid-D provides access to a network of distributed machines with a set of commands as well as software interfaces. It allows simple use of computer and storage facilities and to schedule or monitor compute tasks and data management. It is based on the Globus Toolkit middleware (GT4). Chapter 1 describes the context which led to the demand for advanced software solutions in Astrophysics, and we state the goals of the project. We then present characteristic astrophysical applications that have been implemented on AstroGrid-D in chapter 2. We describe simulations of different complexity, compute-intensive calculations running on multiple sites, and advanced applications for specific scientific purposes, such as a connection to robotic telescopes. We can show from these examples how grid execution improves e.g. the scientific workflow. Chapter 3 explains the software tools and services that we adapted or newly developed. Section 3.1 is focused on the administrative aspects of the infrastructure, to manage users and monitor activity. Section 3.2 characterises the central components of our architecture: The AstroGrid-D information service to collect and store metadata, a file management system, the data management system, and a job manager for automatic submission of compute tasks. We summarise the successfully established infrastructure in chapter 4, concluding with our future plans to establish AstroGrid-D as a platform of modern e-Astronomy.

Evolution of self-gravitating rotating dense stellar systems (e.g. globular clusters, galactic nuclei) with embedded black holes is investigated. The interaction between the black hole and stellar component in differential rotating flattened systems is followed. The interplay between velocity diffusion due to relaxation and black hole star accretion is investigated together with cluster rotation using 2D+1 Fokker-Planck numerical methods. The models can reproduce the Bahcall-Wolf solution $f \propto E^{1/4}$ ($n \propto r^{-7/4}$) inside the zone of influence of the black hole. Gravo-gyro and gravothermal instabilities conduce the system to a faster evolution leading to shorter collapse times with respect to the non-rotating systems. Angular momentum transport and star accretion support the development of central rotation in relaxation time scales. We explore system dissolution due to mass-loss in the presence of an external tidal field (e.g. globular clusters in galaxies).

The chemodynamical evolution of spherical multi-component self-gravitating models for isolated dwarf galaxies is studied. We compare their evolution with and without feedback effects from star formation processes. We find that initially cuspy dark matter profiles flatten with time as a result of star formation, without any special tuning conditions. Thus the seemingly flattened profiles found in many dwarfs do not contradict the cuspy profiles predicted by cosmological models. We also calculate the chemical evolution of stars and gas, to permit comparisons with observational data.

Supermassive black holes (SMBHs) found in the centers of many galaxies have been recognized to play a fundamental active role in the cosmological structure formation process. In hierarchical formation scenarios, SMBHs are expected to form binaries following the merger of their host galaxies. If these binaries do not coalesce before the merger with a third galaxy, the formation of a black hole triple system is possible. Numerical simulations of the dynamics of triples within galaxy cores exhibit phases of very high eccentricity (as high as $e \sim 0.99$). During these phases, intense bursts of gravitational radiation can be emitted at orbital periapsis. This produces a gravitational wave signal at frequencies substantially higher than the orbital frequency. The likelihood of detection of these bursts with pulsar timing and the Laser Interferometer Space Antenna ({\it LISA}) is estimated using several population models of SMBHs with masses $\gtrsim 10^7 {\rm M_\odot}$. Assuming a fraction of binaries $\ge 0.1$ in triple system, we find that few to few dozens of these bursts will produce residuals $>1$ ns, within the sensitivity range of forthcoming pulsar timing arrays (PTAs). However, most of such bursts will be washed out in the underlying confusion noise produced by all the other 'standard' SMBH binaries emitting in the same frequency window. A detailed data analysis study would be required to assess resolvability of such sources. Implementing a basic resolvability criterion, we find that the chance of catching a resolvable burst at a one nanosecond precision level is 2-50%, depending on the adopted SMBH evolution model. On the other hand, the probability of detecting bursts produced by massive binaries (masses $\gtrsim 10^7\msun$) with {\it LISA} is negligible.

The dynamical evolution of binaries of intermediate-massive black holes (IMBHs, massive black holes with a mass ranging between $10^2$ and $10^4 M_{\odot}$) in stellar clusters has recently received an increasing amount of attention. This is at least partially due to the fact that if the binary is hard enough to evolve to the phase at which it will start emitting gravitational waves (GWs) efficiently, there is a good probability that it will be detectable by future space-borne detectors like LISA. We study this evolution in the presence of rotation in the cluster. The eccentricity is strongly connected to the initial IMBHs velocities, and values of $\sim 0.7$ up to 0.9 are reached for low initial velocities, while almost circular orbits result if the initial velocities are increased. A Monte Carlo study indicates that these sources will be detectable by a detector such as LISA with median signal to noise ratios of between 10 and 20 over a three year period, although some events had signal to noise ratios of 300 or greater. Furthermore, one should also be able to estimate the chirp-mass with median fractional errors of $10^{-4}$, reduced mass on the order of $10^{-3}$ and luminosity distance on the order of $10^{-1}$. Finally, these sources will have a median angular resolution in the LISA detector of about 3 square degrees, putting events firmly in the field of view of future electromagnetic detectors such as LSST.

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.

At least 10-15% of nearby sun-like stars have known Jupiter-mass planets. In contrast, very few planets are found in mature open and globular clusters such as the Hyades and 47 Tuc. We explore here the possibility that this dichotomy is due to the post-formation disruption of planetary systems associated with the stellar encounters in long-lived clusters. One supporting piece of evidence for this scenario is the discovery of freely floating low-mass objects in star forming regions. We use two independent numerical approaches, a hybrid Monte Carlo and a direct $N$-body method, to simulate the impact of the encounters. We show that the results of numerical simulations are in reasonable agreement with analytical determinations in the adiabatic and impulsive limits. They indicate that distant stellar encounters generally do not significantly modify the compact and nearly circular orbits. However, moderately close stellar encounters, which are likely to occur in dense clusters, can excite planets' orbital eccentricity and induce dynamical instability in systems which are closely packed with multiple planets. We compute effective cross sections for the dissolution of planetary systems and show that, for all initial eccentricities, dissolution occurs on time scales which are longer than the dispersion of small stellar associations, but shorter than the age of typical open and globular clusters. Although it is much more difficult to disrupt short-period planets, close encounters can excite modest eccentricity among them, such that subsequent tidal dissipation leads to orbital decay, tidal inflation, and even disruption of the close-in planets.

This paper studies the formation and evolution of binary supermassive black holes (SMBHs) in rotating galactic nuclei, focusing on the role of stellar dynamics. We present the first N-body simulations that follow the evolution of the SMBHs from kiloparsec separations all the way to their final relativistic coalescence, and that can robustly be scaled to real galaxies. The N-body code includes post-Newtonian (PN) corrections to the binary equations of motion up to order 2.5; we show that the evolution of the massive binary is only correctly reproduced if the conservative 1PN and 2PN terms are included. The orbital eccentricities of the massive binaries in our simulations are often found to remain large until shortly before coalescence. This directly affects not only their orbital evolution rates, but has important consequences as well for the gravitational waveforms emitted during the relativistic inspiral. We estimate gravitational wave amplitudes when the frequencies fall inside the band of the (planned) Laser Interferometer Space Antennae (LISA). We find significant contributions -- well above the LISA sensitivity curve -- from the higher-order harmonics.

We present preliminary results from self-consistent, high resolution direct {\it N}-body simulations of massive black hole binaries in mergers of galactic nuclei. The dynamics of the black hole binary includes the full Post-Newtonian corrections (up to 2.5PN) to its equations of motion. We show that massive black holes starting at separations of 100 pc can evolve down to gravitational-wave-induced coalescence in less than a Hubble time. The binaries, in our models, often form with very high eccentricity and, as a result, reach separations of 50 Schwarzschild radius with eccentricities which are clearly distinct from zero -- even though gravitational wave emission damps the eccentricity during the inspiral. These deviations from exact circular orbits, at such small separations, may have important consequences for LISA data analysis.

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

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 carried out N-body simulations for rotating star clusters with equal mass and compared the results with Fokker-Planck models. These two different approaches are found to produce fairly similar results, although there are some differences with regard to the detailed aspects. We confirmed the acceleration of the core collapse of a cluster due to an initial non-zero angular momentum and found a similar evolutionary trend in the central density and velocity dispersion in both simulations. The degree of acceleration depends on the initial angular momentum. Angular momentum is being lost from the cluster due to the evaporation of stars with a large angular momentum on a relaxation time scale.

We study the survival of substructures (clumps) within larger self-gravitating dark matter halos. Building on scaling relations obtained from N-body calculations of violent relaxation, we argue that the tidal field of galaxies and halos can only destroy substructures if spherical symmetry is imposed at formation. We explore other mechanisms that may tailor the number of halo substructures during the course of virialization. Unless the larger halo is built up from a few large clumps, we find that clump-clump encounters are unlikely to homogenize the halo on a dynamical timescale. Phase mixing would proceed faster in the inner parts and allow for the secular evolution of a stellar disk.

N-body simulations give us a rough idea of how the shape of a simulated object appears in three-dimensional space. From an observational point of view this may give us a misleading picture. The faint stars may be the most common stars in the system but the morphological information obtained by observations of an object may be dominated by the color properties of the bright stars. Due to dynamical effects, such as energy equipartition, different masses of stars may populate different regions in the object. Since stars are evolving in mass the stellar evolution may also influence the dynamics of the system. Hence, if one is interested in simulating what the morphology will look like through a telescope, one needs to simulate in addition evolving stars and weight them by their luminosity. Therefore we need to combine simulations of the dynamical evolution and a stellar population synthesis at the same time. For the dynamical evolution part we use a parallel version of a direct N-body code, NBODY6++. This code also includes stellar evolution. We link the parameters from this stellar evolution routine to the BaSeL 2.0 stellar library. For isolated clusters we found results assimilable to standard stellar population synthesis codes such as the PEGASE code. For clusters in a tidal field we found that the integrated colors become relatively bluer due to energy equipartition effects. In the time shortly before dissolution of the cluster the stellar M/L ratio becomes lower compared to isolated clusters. We compared the results of our simulations to integrated spectra of galactic globular clusters. For the cluster NGC 1851 we found good agreement between simulation and observation. For extragalactic clusters in M81 and M31 we looked at medium band spectral energy distributions and found for some clusters also a good agreement.

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 address the dynamical evolution of an isolated self--gravitating system with two stellar mass groups. We vary the individual ratio of the heavy to light bodies, $\mu$ from 1.25 to 50 and alter also the fraction of the total heavy mass $\MH$ from 5% to 40% of the whole cluster mass. Clean-cut properties of the cluster dynamics are examined, like core collapse, the evolution of the central potential, as well as escapers. We present in this work collisional $N$-body simulations, using the high--order integrator NBODY6++ with up to ${\cal N}_{\star}=2\cdot 10^4$ particles improving the statistical significancy of the lower--${\cal N}_{\star}$ simulations by ensemble averages. Equipartition slows down the gravothermal contraction of the core slightly. Beyond a critical value of $\mu \approx 2$, no equipartition can be achieved between the different masses; the heavy component decouples and collapses. For the first time the critical boundary between Spitzer--stable and --unstable systems is demonstrated in direct $N$-body models.

Globular clusters rotate significantly, and with the increasing amount of detailed morphologicaland kinematical data obtained in recent years on galactic globular clusters many interesting features show up. We show how our theoretical evolutionary models of rotating clusters can be used to obtain fits, which at least properly model the overall rotation and its implied kinematics in full 2D detail (dispersions, rotation velocities). Our simplified equal mass axisymmetric rotatingmodel provides detailed two-dimensional kinematical and morphological data for star clusters. The degree of rotation is not dominant in energy, but also non-negligible for the phase space distribution function, shape and kinematics of clusters. Therefore the models are well applicable for galactic globular clusters. Since previously published papers on that matter by us made it difficult to do detailed comparisons with observations we provide a much more comprehensive and easy-to-use set of data here, which uses as entries dynamical age and flattening of observed cluster andthen offers a limited range of applicable models in full detail. The method, data structure and some exemplary comparison with observations are presented. Future work will improve modelling anddata base to take a central black hole, a mass spectrum and stellar evolution into account.

Direct-summation N-body algorithms compute the gravitational interaction between stars in an exact way and have a computational complexity of O(N^2). Performance can be greatly enhanced via the use of special-purpose accelerator boards like the GRAPE-6A. However the memory of the GRAPE boards is limited. Here, we present a performance analysis of direct N-body codes on two parallel supercomputers that incorporate special-purpose boards, allowing as many as four million particles to be integrated. Both computers employ high-speed, Infiniband interconnects to minimize communication overhead, which can otherwise become significant due to the small number of "active" particles at each time step. We find that the computation time scales well with processor number; for 2*10^6 particles, efficiencies greater than 50% and speeds in excess of 2 TFlops are reached.

Compact object clusters are likely to exist in the centre of some galaxies because of mass segregation. The high densities and velocities reached in them deserves a better understanding. The formation of binaries and their subsequent merging by gravitational radiation emission is important to the evolution of such clusters. We address the evolution of such a system in a relativistic regime. The recurrent mergers at high velocities create an object with a mass much larger than the average. For this aim we modified the direct {\sc Nbody6}++ code to include post-Newtonian effects to the force during two-body encounters. We adjusted the equations of motion to include for the first time the effects of both periastron shift and energy loss by emission of gravitational waves and so to study the eventual decay and merger of radiating binaries. The method employed allows us to give here an accurate post-Newtonian description of the formation of a run-away compact object by successive mergers with surrounding particles, as well as the distribution of characteristic eccentricities in the events. This study should be envisaged as a first step towards a detailed, accurate study of possible gravitational waves sources thanks to the combination of the direct {\sc Nbody} numerical tool with the implementation of post-Newtonian terms on it.

Binary supermassive black holes form naturally in galaxy mergers, but their long-term evolution is uncertain. In spherical galaxies, N-body simulations show that binary evolution stalls at separations much too large for significant emission of gravitational waves (the "final parsec problem"). Here, we follow the long-term evolution of a massive binary in more realistic, triaxial and rotating galaxy models. We find that the binary does not stall. The binary hardening rates that we observe are sufficient to allow complete coalescence of binary SBHs in 10 Gyr or less, even in the absence of collisional loss-cone refilling or gas-dynamical torques, thus providing a potential solution to the final parsec problem.

The relevance of encounters on the destruction of protoplanetary discs in the Orion Nebula Cluster (ONC) is investigated by combining two different types of numerical simulation. First, star-cluster simulations are performed to model the stellar dynamics of the ONC, the results of which are used to investigate the frequency of encounters, the mass ratio and separation of the stars involved, and the eccentricity of the encounter orbits. The results show that interactions that could influence the star-surrounding disc are more frequent than previously assumed in the core of the ONC, the so-called Trapezium cluster. Second, a parameter study of star-disc encounters is performed to determine the upper limits of the mass loss of the discs in encounters. For simulation times of $\sim$ 1-2 Myr (the likely age of the ONC) the results show that gravitational interaction might account for a significant disc mass loss in dense clusters. Disc destruction is dominated by encounters with high-mass stars, especially in the Trapezium cluster, where the fraction of discs destroyed due to stellar encounters can reach 10-15%. These estimates are in accord with observations of (Lada et al. 2000) who determined a stellar disc fraction of 80-85%. Thus, it is shown that in the ONC - a typical star-forming region - stellar encounters do have a significant effect on the mass of protoplanetary discs and thus affect the formation of planetary systems.

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

We use direct-summation N-body integrations to follow the evolution of binary black holes at the centers of galaxy models with large, constant-density cores. Particle numbers as large as 400K are considered. The results are compared with the predictions of loss-cone theory, under the assumption that the supply of stars to the binary is limited by the rate at which they can be scattered into the binary's influence sphere by gravitational encounters. The agreement between theory and simulation is quite good; in particular, we are able to quantitatively explain the observed dependence of binary hardening rate on N. We do not verify the recent claim of Chatterjee, Hernquist & Loeb (2003) that the hardening rate of the binary stabilizes when N exceeds a particular value, or that Brownian wandering of the binary has a significant effect on its evolution. When scaled to real galaxies, our results suggest that massive black hole binaries in gas-poor nuclei would be unlikely to reach gravitational-wave coalescence in a Hubble time.

We investigate the dynamical decay of non-hierarchical accreting triple systems and its implications on the ejection model as Brown Dwarf formation scenario. A modified chain-regularization scheme is used to integrate the equations of motion, that also allows for mass changes over time as well as for momentum transfer from the accreted gas mass onto the bodies. We integrate an ensemble of triple systems within a certain volume with different accretion rates, assuming several prescriptions of how momentum is transferred onto the bodies. We follow their evolution until the systems have decayed. We analyze the end states and decay times of these systems and determine the fraction of Brown Dwarfs formed, their escape speeds as well as the semi-major axis distribution of the formed Brown Dwarf binaries. We find that the formation probability of Brown Dwarfs depends strongly on the assumed momentum transfer which is related to the motion of the gas. Due to ongoing accretion and consequent shrinkage of the systems, the median escape velocity is increased by a factor of 2 and the binary separations are decreased by a factor of 5 compared with non-accreting systems. Furthermore, the obtained semi-major axis distribution drops off sharply to either side of the median, which is also supported by observations. We conclude that accretion and momentum transfer of accreted gas during the dynamical decay of triple systems is able to produce the observed distribution of close binary Brown Dwarfs, making the ejection model a viable option as Brown Dwarf formation scenario.

This paper has been withdrawn temporarily by the authors. As stars close to the galactic centre have short orbital periods it has been possible to trace large fractions of their orbits in the recent years. Previously the data of the orbit of the star S2 have been fitted with Keplerian orbits corresponding to a massive black hole (MBH) with a mass of M$_{BH}$=3-4$\times10^6$M_sun implying an insignificant cusp mass.However, it has also been shown that the central black hole resides in a ~1'' diameter stellar cluster of a priori unknown mass. In a spherical potential which is neither Keplerian nor harmonic, orbits will precess resulting in inclined rosetta shaped trajectories on the sky. In this case, the assumption of non-Keplerian orbits is a more physical approach. It is also the only approach through which cusp mass information can be obtained via stellar dynamics of the cusp members. This paper presents the first exemplary modelling efforts in this direction. Using positional and radial data of star S2, we find that there could exist an unobserved extended mass component of several 10$^5$M$_{\odot}$ forming a so-called 'cusp' centred on the black hole position. Considering only the fraction of the cusp mass M$_{S2_{apo}}$ within the apocenter of the S2 orbit we find as an upper limit that M$_{S2_{apo}}$/(M$_{BH}$ + M$_{S2_{apo}}$) $\le$ 0.05. A large extended cusp mass, if present, is unlikely to be composed of sub-solar mass constituents, but could be explained rather well by a cluster of high M/L stellar remnants, which we find to form a stable configuration.

We present a clear N-body realization of the growth of a Bahcall-Wolf f ~ E^0.25 (rho ~ 1/r^1.75) density cusp around a massive object ("black hole") at the center of a stellar system. Our N-body algorithm incorporates a novel implementation of Mikkola-Aarseth chain regularization to handle close interactions between star and black hole particles. Forces outside the chain were integrated on a GRAPE-6A/8 special-purpose computer with particle numbers up to N=0.25 x 10^6. We compare our N-body results with predictions of the isotropic Fokker-Planck equation and verify that the time dependence of the density (both configuration and phase-space) predicted by the Fokker-Planck equation is well reproduced by the N-body algorithm, for various choices of N and of the black hole mass. Our results demonstrate the feasibility of direct-force integration techniques for simulating the evolution of galactic nuclei on relaxation time scales.

To be presented is a study of the secular evolution of a spherical stellar system with a central star-accreting black hole (BH) using the anisotropic gaseous model. This method solves numerically moment equations of the full Fokker-Planck equation, with Boltzmann-Vlasov terms on the left- and collisional terms on the right-hand sides. We study the growth of the central BH due to star accretion at its tidal radius and the feedback of this process on to the core collapse as well as the post-collapse evolution of the surrounding stellar cluster in a self-consistent manner. Diffusion in velocity space into the loss-cone is approximated by a simple model. The results show that the self-regulated growth of the BH reaches a certain fraction of the total mass cluster and agree with other methods. Our approach is much faster than competing ones (Monte Carlo, $N$--body) and provides detailed informations about the time and space dependent evolution of all relevant properties of the system. In this work we present the method and study simple models (equal stellar masses, no stellar evolution or collisions). Nonetheless, a generalisation to include such effects is conceptually simple and under way.

First we present a recently developed 3D chemodynamical code for galaxy evolution from the K**2 collaboration. It follows the evolution of all components of a galaxy such as dark matter, stars, molecular clouds and diffuse interstellar matter (ISM). Dark matter and stars are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse (hot) gas and a sticky particle scheme for the (cool) molecular clouds. Physical processs such as star formation, stellar death or condensation and evaporation processes of clouds interacting with the ISM are described locally. An example application of the model to a star forming dwarf galaxy will be shown for comparison with other codes. Secondly we will discuss new kinds of exotic chemodynamical processes, as they occur in dense gas-star systems in galactic nuclei, such as non-standard ``drag''-force interactions, destructive and gas producing stellar collisions. Their implementation in 1D dynamical models of galactic nuclei is presented. Future prospects to generalize these to 3D are work in progress and will be discussed.

We have studied the dynamical evolution of rotating star clusters with mass spectrum using a Fokker-Planck code. As a simplest multi-mass model, we first investigated the two-component clusters. Rotation is found to accelerate the dynamical evolution through the transfer of angular momentum outward, as well as from the high masses to the low masses. However, the degree of acceleration depends sensitively on the assumed initial mass function since dynamical friction, which generates mass segregation, also tends to accelerate the evolution, and the combined effect of both is not linear or multiplicative. As long as dynamical friction dominates in the competition with angular momentum exchange the heavy masses lose random energy and angular momentum, sink towards the centre, but their remaining angular momentum is sufficient to speed them up rotationally. This is gravo-gyro instability. As a consequence, we find that the high mass stars in the central parts rotate faster than low mass stars. This leads to the suppression of mass segregation compared to the non-rotating clusters. From the study of multi-component models, we observe similar trends to the two-component models in almost all aspects. The mass function changes less drastically for clusters with rotation. Unlike non-rotating clusters, the mass function depends on $R$ and $z$. Our models are the only ones that can predict mass function and other quantities to be compared with new observations.

In the 60s and 70s super-massive central objects (from now onwards SMOs) were thought to be the main source of active galactic nuclei (AGNs) characteristics (luminosities of $L \approx 10^{12} L_{\odot}$). The release of gravitational binding energy by the accretion of material on to an SMO in the range of $10^7 - 10^9 M_{\odot}$ has been suggested to be the primary powerhouse (Lynden-Bell 1969). That rather exotic idea in early time has become common sense nowadays. Not only our own galaxy harbours a few million-solar mass black hole (Genzel 2001) but also many of other non-active galaxies show kinematic and gas-dynamic evidence of these objects (Magorrian et al. 1998) The concept of central super-massive stars (SMSs henceforth) (${\cal M} \ge 5 \times 10^4 M_{\odot}$, where ${\cal M}$ is the mass of the SMS) embedded in dense stellar systems was suggested as a possible explanation for high- energy emissions phenomena occurring in AGNs and quasars (Vilkoviski 1976, Hara 1978), such as X-ray emissions (Bahcall and Ostriker, 1975). SMSs and super-massive black holes (SMBHs) are two possibilities to explain the nature of SMOs, and SMSs may be an intermediate step towards the formation of SMBHs (Rees 1984). In this paper we give the equations that describe the dynamics of such a dense star-gas system which are the basis for the code that will be used in a prochain future to simulate this scenario. We also briefly draw the mathematical fundamentals of the code

The concept of central super-massive stars (${\cal M} \ge 5 \times 10^4 M_{\odot}$, where ${\cal M}$ is the mass of the super-massive star) embedded in dense stellar systems was suggested as a possible explanation for high- energy emissions phenomena occurring in active galactic nuclei and quasars (Vilkoviski 1976, Hara 1978), such as X-ray emissions (Bahcall and Ostriker, 1975). SMSs and super-massive black holes are two possibilities to explain the nature of super-massive central objects, and super-massive stars may be an intermediate step towards the formation of super-massive black holes (Rees 1984). Therefore it is important to study such a dense gas-star system in detail. We address here the implementation of radiative transfer in a model which was presented in former work (Amaro-Seoane and Spurzem 2001, Amaro-Seoane et al. 2002). In this sense, we extend here and improve the work done by Langbein et al. (1990) by describing the radiative transfer in super-massive stars using previous work on this subject (Castor 1972).

Spherically symmetric equal mass star clusters containing a large amount of primordial binaries are studied using a hybrid method, consisting of a gas dynamical model for single stars and a Monte Carlo treatment for relaxation of binaries and the setup of close resonant and fly-by encounters of single stars with binaries and binaries with each other (three- and four-body encounters). What differs from our previous work is that each encounter is being integrated using a highly accurate direct few-body integrator which uses regularized variables. Hence we can study the systematic evolution of individual binary orbital parameters (eccentricity, semi-major axis) and differential and total cross sections for hardening, dissolution or merging of binaries (minimum distance) from a sampling of several ten thousands of scattering events as they occur in real cluster evolution including mass segregation of binaries, gravothermal collapse and reexpansion, binary burning phase and ultimately gravothermal oscillations. For the first time we are able to present empirical cross sections for eccentricity variation of binaries in close three- and four-body encounters. It is found that a large fraction of three-body and four-body encounters results in merging. Previous cross sections obtained by Spitzer and Gao for strong encounters can be reproduced, while for weak encounters non-standard processes like formation of hierarchical triples occur.

In this paper we present some test results of our newly developed Multi-Phase Chemo-Dynamical Smoothed Particle Hydrodynamics (MP- CD-SPH) code for galaxy evolution. At first, we present a test of the ``pure'' hydro SPH part of the code. Then we describe and test the multi-phase description of the gaseous components of the interstellar matter. In this second part we also compare our condensation and evaporation description with the results of a previous 2d multi-phase hydrodynamic mesh code.

We present our recently developed 3-dimensional chemodynamical code for galaxy evolution. It follows the evolution of all components of a galaxy such as dark matter, stars, molecular clouds and diffuse interstellar matter (ISM). Dark matter and stars are treated as collisionless $N$-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse (hot) gas and a sticky particle scheme for the (cool) molecular clouds. Additionally, the galactic components are coupled by several phase transitions like star formation, stellar death or condensation and evaporation processes within the ISM. As an example here we present the dynamical, chemical and photometric evolution of a star forming dwarf galaxy with a total baryonic mass of $2 \times 10^9 {\rm M}_\odot$.

We follow the sinking of two massive black holes in a spherical stellar system where the black holes become bound under the influence of dynamical friction. Once bound, the binary hardens by three-body encounters with surrounding stars. We find that the binary wanders inside the core, providing an enhanced supply of reaction partners for the hardening. The binary evolves into a highly eccentric orbit leading to coalescence well beyond a Hubble time. These are the first results from a hybrid ``self consistent field'' (SCF) and direct Aarseth N-body integrator (NBODY6), which combines the advantages of the direct force calculation with the efficiency of the field method. The code is designed for use on parallel architectures and is therefore applicable to collisional N-body integrations with extraordinarily large particle numbers (> 10^5). This creates the possibility of simulating the dynamics of both globular clusters with realistic collisional relaxation and stellar systems surrounding supermassive black holes in galactic nuclei.

Challenging stellar dynamical problems, such as the study of gravothermal oscillations in star clusters, have in the past initiated the very successful building of GRAPE special purpose computers. It is discussed, that present day tasks such as the formation and evolution of galactic nuclei with one or more massive black holes and the coupled stellar and gas dynamical processes in the formation of nuclei and star clusters, demand a new kind of hybrid architecture, using both GRAPE and a reconfigurable logics board called RACE. For such a system we have developed first implementations and floating point performance studies in the case of the SPH algorithm (smoothed particle hydrodynamics), which will be of great advantage for SPH modelling and also for direct $N$-body simulations using the more efficient Ahmad-Cohen neighbour schemes.

We present our recently developed 3-dimensional chemodynamical code for galaxy evolution. This code follows the evolution of different galactic components like stars, dark matter and different components of the interstellar medium (ISM), i.e. a diffuse gaseous phase and the molecular clouds. Stars and dark matter are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse gas and a sticky particle scheme for the molecular clouds. Additionally, the galactic components are coupled by several phase transitions like star formation, stellar death or condensation and evaporation processes within the ISM. As an example we show the dynamical and chemical evolution of a star forming dwarf galaxy with a total baryonic mass of $2 \cdot 10^9 M_\odot$. After a moderate collapse phase the stars and the molecular clouds follow an exponential radial distribution, whereas the diffuse gas shows a central depression as a result of stellar feedback. The metallicities of the galactic components behave quite differently with respect to their temporal evolution as well as their radial distribution. Especially, the ISM is at no stage well mixed.

We use a gaseous model and a semi-analytical approach to study the evolution of a super-massive central gaseous object (a super-massive star, (SMS from now on) in an AGN and its evolution by interactions with the surrounding stellar system. Our future work in this field is outlined, which aims at a more detailed study of energy flows in the interstellar medium, stellar evolution and the relation between QSOs and galaxy formation.

We present the first post core collapse models of initially rotating star clusters, using the numerical solution of an orbit-averaged 2D Fokker-Planck equation. Based on the code developed by Einsel & Spurzem (1999), we have improved the speed and the stability and included the steady three-body binary heating source. We have confirmed that rotating clusters, whether they are in a tidal field or not, evolve significantly faster than non-rotating ones. Consequences for observed shapes, density distribution, and kinematic properties of young and old star clusters are discussed. The results are compared with gaseous and 1D Fokker-Planck models in the non-rotating case.

We address the classical problem of star accretion onto a supermassive central gaseous object in a galactic nucleus. The resulting supermassive central gas-star object is assumed to be located at the centre of a dense stellar system for which we use a simplified model consisting of a Plummer model with an embedded density cusp using stellar point masses. From the number of stars belonging to the loss-cone, which plunge onto the central object on elongated orbits from outside, we estimate the accretion rate taking into account a possible anisotropy of the surrounding stellar distribution. The total heating rate in the supermassive star due to the loss-cone stars plunging onto it is estimated. This semi-analytical study, revisiting and expanding classical paper's work, is a starting point of future work on a more detailed study of early evolutionary phases of galactic nuclei. It merits closer examination, because it is one of the key features for the link between cosmology and galaxy formation.

This is a starting point of future work on a more detailed study of early evolutionary phases of galactic nuclei. We put into port the study of the evolution of star accretion onto a supermassive gaseous object in the central region of an active galactic nucleus, which was previously addressed semi-analytically. For this purpose, we use a gaseous model of relaxing dense stellar systems, whose equations are solved numerically. The model is shortly described and some first model calculations for the specific case of black hole mass growth in time are given. Similar concepts will be applied in future work for a system in which there is a dense interstellar medium as a possible progenitor of a black hole in the galactic centre.

Interacting galaxies like the famous Antennae (NGC 4038/4039) or Stephan's Quintet (HCG 92) show considerable star forming activity in their tidal arms. High resolution images (e.g. from HST-observations) indicate that these regions consist of up to hundreds of massive stellar clusters or tidal dwarf galaxies (TDG). In this paper we want to investigate the future fate of these clusters of massive star clusters (in this work called super-clusters). We simulate compact super-clusters in the tidal field of a host-galaxy and investigate the influence of orbital and internal parameters on the rate and timescale of the merging process. We show that it is possible that such configurations merge and build a dwarf galaxy, which could be an important mechanism of how long-lived dwarf satellite galaxies form. A detailed study of the merger object will appear in a follow-up paper.

We present Superbox, a particle-mesh code with high resolution sub-grids and an NGP (nearest grid point) force-calculation scheme based on the second derivatives of the potential. Superbox implements a fast low-storage FFT-algorithm, giving the possibility to work with millions of particles on desk-top computers. Test calculations show energy and angular momentum conservation to one part in 10^5 per crossing-time. The effects of grid and numerical relaxation remain negligible, even when these calculations cover a Hubble-time of evolution. As the sub-grids follow the trajectories of individual galaxies, the code allows a highly resolved treatment of interactions in clusters of galaxies, such as high-velocity encounters between elliptical galaxies and the tidal disruption of dwarf galaxies. Excellent agreement is obtained in a comparison with a direct-summation N-body code running on special-purpose Grape3 hardware. The orbital decay of satellite galaxies due to dynamical friction obtained with Superbox agrees with Chandrasekhar's treatment when the Coulomb logarithm is approximately 1.5.

The new approach outlined in Paper I (Spurzem \& Giersz 1996) to follow the individual formation and evolution of binaries in an evolving, equal point-mass star cluster is extended for the self-consistent treatment of relaxation and close three- and four-body encounters for many binaries (typically a few percent of the initial number of stars in the cluster). The distribution of single stars is treated as a conducting gas sphere with a standard anisotropic gaseous model. A Monte Carlo technique is used to model the motion of binaries, their formation and subsequent hardening by close encounters, and their relaxation (dynamical friction) with single stars and other binaries. The results are a further approach towards a realistic model of globular clusters with primordial binaries without using special hardware. We present, as our main result, the self-consistent evolution of a cluster consisting of 300.000 equal point-mass stars, plus 30.000 equal mass binaries over several hundred half-mass relaxation times, well into the phase where most of the binaries have been dissolved and evacuated from the core. In a self-consistent model it is the first time that such a realistically large number of binaries is evolving in a cluster with an even ten times larger number of single stars. Due to the Monte Carlo treatment of the binaries we can at every moment analyze their external and internal parameters in the cluster as in an N-body simulation.

In astrophysics numerical star cluster simulations and hydrodynamical methods like SPH require computational performance in the petaflop range. The GRAPE family of ASIC-based accelerators improves the cost-performance ratio compared to general purpose parallel computers, however with limited flexibility. The AHA-GRAPE architecture adds a reconfigurable FPGA-processor to accelerate the SPH computation. The basic equations of the algorithm consist of three parts each scaling with the order of O(N), O(N*Nn), and O(N**2) respectively, where N is in the range of 10**4 to 10**7 and Nn approximately 50. These equations can profitably be distributed across a host workstation, an FPGA processor, and a GRAPE subsystem. With the new ATLANTIS FPGA processor we expect a scalable SPH performance of 1.5 Gflops per board. The first prototype AHA-GRAPE system will be available in mid-2000. This 3-layered system will deliver an increase in performance by a factor of 10 as compared to a pure GRAPE solution.

Special high-accuracy direct force summation N-body algorithms and their relevance for the simulation of the dynamical evolution of star clusters and other gravitating N-body systems in astrophysics are presented, explained and compared with other methods. Other methods means here approximate physical models based on the Fokker-Planck equation as well as other, approximate algorithms to compute the gravitational potential in N-body systems. Questions regarding the parallel implementation of direct ``brute force'' N-body codes are discussed. The astrophysical application of the models to the theory of relaxing rotating and non-rotating collisional star clusters is presented, briefly mentioning the questions of the validity of the Fokker-Planck approximation, the existence of gravothermal oscillations and of rotation and primordial binaries.

Astrophysical Challenges which demand the solution of the one million (or more) gravitating body problem are briefly discussed for the fields of cosmology, galactic nuclei and globular star clusters. Results from the classical three-body problem are to be combined with methods of statistical mechanics or thermodynamics in order to provide a reasonable model for astrophysically relevant particle numbers of the order 10**4-10**8. For high-precision modelling of star clusters and galactic nuclei over the required hundreds and thousands of crossing times specially tailored ``brute-force'' N-body integration methods are used together with special purpose hardware. A recent implementation of such code on a general purpose parallel computer is presented. After a comparison of the relevance of different N-body algorithms a new concept for a more flexible customized special purpose computer based on a combination of GRAPE and FPGA is proposed. It is an ideal machine for all kinds of N-body simulations using neighbour schemes, as the Ahmad-Cohen direct N-body codes and smoothed particle hydrodynamics (SPH) for systems including interstellar gas dynamics.

1- Simulations of star clusters and galactic nuclei: rotation and other challenges (R. Spurzem & C. Einsel) 2- Relaxation of elliptical galaxies (C. Boily) 3- Simulations with an Hybrid Nbody-Scf code (M. Hemsendorf) These are three two-page contributions to the Rutgers galaxy dynamics held in august, 1998. Each will be included in the proceedings to be published by the Astronomical Society of the Pacific.

A database on shape evolution of direct N-body models formed out of cold, dissipationless collapse is generated using GRAPE and HARP special purpose computers. Such models are important to understand the formation of elliptical galaxies. Three dynamically distinct phases of shape evolution were found, first a fast dynamical collapse which gives rise to the radial orbit instability (ROI) and generates at its end the maximal triaxiality of the system. Subsequently, two phases of violent and two-body shape relaxation occur, which drive the system first towards axisymmetry, finally to spherical symmetry (the final state, however, is still much more concentrated than the initial model). In a sequence of models the influence of numerical and physical parameters, like particle number, softening, initial virial ratio, timestep choice, different N-body codes, are examined. We find that an improper combination of softening and particle number can produce erroneous results. Selected models were evolved on the secular timescale until they became spherically symmetric again. The secular shape relaxation time scale is shown to agree very well with the two-body relaxation time, if softening is properly taken into account for the latter. Finally, we argue, that the intermediate phase of violent shape relaxation after collapse is induced by strong core oscillations in the centre, which cause potential fluctuations, dampening out the triaxiality.

The subjects and key questions faced by computational astrophysics using N-body simulations are discussed in the fields of globular star cluster dynamics, galactic nuclei and cosmological structure formation. After a comparison of the relevance of different N-body algorithms a new concept for a more flexible customized special purpose computer based on a combination of GRAPE and FPGA is proposed. It is an ideal machine for all kinds of N-body simulations using neighbour schemes, as the Ahmad-Cohen direct N-body codes and smoothed particle hydrodynamics for systems including gas.

We review recent progress in computational methods for studying the dynamical evolution of large star clusters. After describing some highlights of a new comparison of different models, including N-body models, we mention in detail recent developments in three codes which allow models with anisotropic distributions of velocity to be studied: (i) a Fokker-Planck code using finite differences, (ii) a Monte Carlo code and (iii) a gaseous model code.

The influence of rotation on the dynamical evolution of collisional stellar systems is investigated by solving the orbit-averaged Fokker-Planck equation in (E,J_z)-space. We find that large amounts of initial rotation drive the system into a phase of strong mass loss while it is moderately contracting. The core is rotating even faster than before although angular momentum is transported outwards. At the same time the core is heating. Given these features this phase can be associated with the gravo-gyro 'catastrophe' found by Hachisu (1979). The increase in central angular momentum levels off after about 2-3 initial half-mass relaxation times indicating that the source of this `catastrophe' is depleted. Finally, the central angular velocity increases again, but with a rather small power of the central density -- the same power as for the central velocity dispersion during self-similar contraction towards the collapse. The rotation curves flatten and the ellipticity variations decrease with time, but their shapes are very similar. These results suggest the existence of a self-similar solution for a rotating cluster as well. The maximum values of rotational velocity and ellipticity occur at about the half mass radius.

A collisional $N$-body simulation using NBODY5 on a single CRAY YMP processor is followed well into the post-collapse regime. This is presently one of the largest particle numbers of all such models published, but some data for an even larger $N$ produced by using special-purpose computers have recently been presented. In contrast to previous ensemble-averaged $N$-body simulations the noise here is low enough to just compare this one single run with the expectations from statistical models based on the Fokker-Planck approximation. Agreement is as good as could be expected for the case of the evolution of the Lagrangian radii, radial and tangential velocity dispersions and various core quantities. We briefly discuss approximate models to understand number and energy of escapers and the question of gravothermal core oscillations; although the system exhibits post-collapse oscillations they turn out to be directly binary driven and we cannot prove the existence of a gravothermal expansion at this particle number. Finally in a detailed examination of the wandering of the density centre we find in contrast to some previous studies a clear long-time period of the order of about 14 half-mass crossing times.

We compare the results for the dynamical evolution of star clusters derived from anisotropic gaseous models with the data from N-body simulations of isolated and one-component systems, each having modest number of stars. The statistical quality of N-body data was improved by averaging results from many N-body runs, each with the same initial parameters but with different sequences of random numbers used to initialize positions and velocities of the particles. We study the development of anisotropy, the spatial evolution and energy generation by three-body binaries and its N-dependence. We estimate the following free parameters of anisotropic gaseous models: the time scale for collisional anisotropy decay and the coefficient in the formulae for energy generation by three-body binaries. To achieve a fair agreement between N-body and gaseous models for the core in pre- as well as in post-collapse only the energy generation by binaries had to be varied by N. We find that anisotropy has considerable influence on the spatial structure of the cluster particularly for the intermediate and outer regions.