Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published October 4, 2002 | Published + Submitted
Journal Article Open

A new Monte Carlo code for star cluster simulations II. Central black hole and stellar collisions

Abstract

We have recently written a new code to simulate the long term evolution of spherical clusters of stars. It is based on the pioneering Monte Carlo scheme proposed by Hénon in the 70's. Unlike other implementations of this numerical method which were successfully used to investigate the dynamics of globular clusters, our code has been devised in the specific goal to treat dense galactic nuclei. In a previous paper, we described the basic version of our code which includes 2-body relaxation as the only physical process. In the present work, we go on and include further physical ingredients that are mostly relevant to galactic nuclei, namely the presence of a central (growing) black hole (BH) and collisions between (main sequence) stars. Stars that venture too close to the BH are destroyed by the tidal field. We took particular care of this process because of its importance, both as a channel to feed the BH and a way to produce accretion flares from otherwise quiescent galactic nuclei. Collisions between stars have often been proposed as another mechanism to drive stellar matter into the central BH. Furthermore, non disruptive collisions may create peculiar stellar populations which are of great observational interest in the case of the central cluster of our Galaxy. To get the best handle on the role of this process in galactic nuclei, we include it with unpreceded realism through the use of a set of more than 10 000 collision simulations carried out with a SPH (Smoothed Particle Hydrodynamics) code. Stellar evolution has also been introduced in a simple way, similar to what has been done in previous dynamical simulations of galactic nuclei. To ensure that this physics is correctly simulated, we realized a variety of tests whose results are reported here. This unique code, featuring most important physical processes, allows million particle simulations, spanning a Hubble time, in a few CPU days on standard personal computers and provides a wealth of data only rivalized by N-body simulations.

Additional Information

© 2002 ESO. Received: 17 April 2002. Accepted: 31 July 2002. Stellar models to initiate SPH collision simulations where kindly provided by the Geneva stellar evolution group, with precious help from Georges Meynet, and by Isabelle Baraffe and Corinne Charbonel. M.F. wants to thank Rainer Spurzem for interesting discussions and Gerald Quinlan and Roeland van der Marel for providing pycode, the BH adiabatic growth code. Comments by the anonymous referee helped to clarify the paper. Most simulations have been realized on the "Beowulfs'' clusters GRAVITOR at Geneva Observatory[*], and ISIS at Bern University. This work has been supported by the Swiss National Science Foundation. The writing of this paper has been finished at Caltech, with partial support from NASA under grant NAG5-10707.

Attached Files

Published - aa2593.pdf

Submitted - 0204292.pdf

Files

0204292.pdf
Files (3.8 MB)
Name Size Download all
md5:08dd1e7756579534f2e68b83837a4251
1.8 MB Preview Download
md5:e08e67e20e0c63c74cda6ddc79afa505
2.0 MB Preview Download

Additional details

Created:
August 19, 2023
Modified:
October 18, 2023