期刊
ASTROPHYSICAL JOURNAL
卷 540, 期 2, 页码 969-982出版社
IOP PUBLISHING LTD
DOI: 10.1086/309350
关键词
celestial mechanics, stellar dynamics; globular clusters : general; methods : n-body simulations; methods : numerical
We present a new parallel supercomputer implementation of the Monte Carlo method for simulating the dynamical evolution of globular star clusters. Our method is based on a modified version of Henon's Monte Carlo algorithm for solving the Fokker-Planck equation. Our code allows us to follow the evolution of a cluster containing up to 5 x 10(5) stars to core collapse in less than or similar to 40 hours of computing time. In this paper we present the results of test calculations for clusters with equal-mass stars, starting from both Plummer and King model initial conditions. We consider isolated as well as tidally truncated clusters. Our results are compared to those obtained from approximate, self-similar analytic solutions, from direct numerical integrations of the Fokker-Planck equation, and from direct N-body integrations performed on a GRAPE-4 special-purpose computer with N = 16384. In all cases we find excellent agreement with other methods, establishing our new code as a robust tool for the numerical study of globular cluster dynamics using a realistic number of stars.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据