Accelerating NBODY6 with a GPU-Enabled Particle-Particle Particle-Tree Scheme

We describe a modified version of the NBODY6 code for simulating star clusters which greatly improves computational efficiency while sacrificing little in the way of accuracy. The distant force calculator is replaced by a GPU-enabled Barnes-Hut code, and integration is done with a standard leap frog scheme. Short-range forces continue to use the CPU-based fourth-order Hermite predictor-corrector scheme of NBODY6. Our code outperforms NBODY6 for systems with more than $3\times {10}^5$ particles and runs more than a factor 2 faster for systems of ${10}^6$ particles with similar energy conservation. Our code should be useful for simulating realistic dense stellar clusters, such as globular clusters or galactic nuclei.