BRIDGE: A Direct-Tree Hybrid N -Body Algorithm for Fully Self-Consistent Simulations of Star Clusters and Their Parent Galaxies

We developed a new direct-tree hybrid N -body algorithm for fully self-consistent N -body simulations of star clusters in their parent galaxies. In such simulations, star clusters need high accuracy, while galaxies need a fast scheme because of the large number of particles required to model it. In our new algorithm, the internal motion of the star cluster is calculated accurately using the direct Hermite scheme with individual timesteps, and all other motions are calculated using the tree code with a second-order leapfrog integrator. The direct and tree schemes are combined using an extension of the mixed variable symplectic (MVS) scheme. Thus, the Hamiltonian corresponding to everything other than the internal motion of the star cluster is integrated with the leapfrog, which is symplectic. Using this algorithm, we performed fully self-consistent N -body simulations of star clusters in their parent galaxy. The internal and orbital evolutions of the star cluster agreed well with those obtained using the direct scheme. We also performed fully self-consistent N -body simulation for large-N models (N = 2 � 10 6 ). In this case, the calculation speed was seven-times faster than what would be if the direct scheme was used.

[1]  Jessica R. Lu,et al.  Orbits and Origins of the Young Stars in the Central Parsec , 2006 .

[2]  Simon F. Portegies Zwart,et al.  The Runaway Growth of Intermediate-Mass Black Holes in Dense Star Clusters , 2002, astro-ph/0201055.

[3]  Stephen L. W. McMillan,et al.  An O(N log N) integration scheme for collisional stellar systems , 1993 .

[4]  Hiroshi Nakai,et al.  Symplectic integrators and their application to dynamical astronomy , 1990 .

[5]  Stuart L. Shapiro,et al.  Random Gravitational Encounters and the Evolution of Spherical Systems. III. Halo , 1971 .

[6]  Kiel,et al.  The efficiency of the spiral- in of a black hole to the Galactic centre , 2002, astro-ph/0212494.

[7]  Simon Portegies Zwart,et al.  The Origin of IRS 16: Dynamically Driven In-Spiral of a Dense Star Cluster to the Galactic Center? , 2003, astro-ph/0303599.

[8]  H. Cohn,et al.  Late core collapse in star clusters and the gravothermal instability , 1980 .

[9]  Piet Hut,et al.  A hierarchical O(N log N) force-calculation algorithm , 1986, Nature.

[10]  J. Wisdom,et al.  Symplectic maps for the N-body problem. , 1991 .

[11]  Junichiro Makino,et al.  A Fast Parallel Treecode with GRAPE , 2004 .

[12]  J. Makino,et al.  Dynamical Friction on Satellite Galaxies , 2005, astro-ph/0511651.

[13]  T. Paumard,et al.  The Two Young Star Disks in the Central Parsec of the Galaxy: Properties, Dynamics, and Formation , 2006, astro-ph/0601268.

[14]  Piet Hut,et al.  Use of Supercomputers in Stellar Dynamics , 1986 .

[15]  J. Makino,et al.  GRAPE-6A: A Single-Card GRAPE-6 for Parallel PC-GRAPE Cluster Systems , 2005, astro-ph/0504407.

[16]  Alfred Krabbe,et al.  The Nuclear Cluster of the Milky Way: Star Formation and Velocity Dispersion in the Central 0.5 Parsec , 1995 .

[17]  T. Paumard,et al.  New results on the helium stars in the galactic center using BEAR spectro-imagery , 2000, astro-ph/0011215.

[18]  L. Hernquist,et al.  TREESPH: A Unification of SPH with the Hierarchical Tree Method , 1989 .

[19]  Mark Morris,et al.  Dynamical Friction on Star Clusters near the Galactic Center , 2003, astro-ph/0307271.

[20]  Subrahmanyan Chandrasekhar,et al.  Dynamical friction. I. General considerations: the coefficient of dynamical friction , 1943 .

[21]  Naoto Kobayashi,et al.  Object 17: Another Cluster of Emission-Line Stars Near the Galactic Center , 1995 .

[22]  G. Neugebauer,et al.  The First Measurement of Spectral Lines in a Short-Period Star Bound to the Galaxy’s Central Black Hole: A Paradox of Youth , 2003 .

[23]  J. Makino,et al.  To Circularize or Not To Circularize?—Orbital Evolution of Satellite Galaxies , 2002, astro-ph/0208452.

[24]  Piet Hut,et al.  Core radius and density measurements in N-body experiments Connections with theoretical and observational definitions , 1985 .

[25]  Toshiyuki Fukushige,et al.  GRAPE-6: Massively-Parallel Special-Purpose Computer for Astrophysical Particle Simulations , 2003, astro-ph/0310702.

[26]  Yuri Levin,et al.  Stellar Disk in the Galactic Center: A Remnant of a Dense Accretion Disk? , 2003, astro-ph/0303436.

[27]  Frederic A. Rasio,et al.  The Disruption of stellar clusters containing massive black holes near the Galactic Center , 2005 .