Monte Carlo Stellar Dynamics near Massive Black Holes: Two-dimensional Fokker-Planck Solutions of Multiple Mass Components

In this study we present a novel Monte Carlo code, referred to as GNC, which enables the investigation of dynamical relaxation in clusters comprising multiple mass components in the vicinity of supermassive black holes at the centers of galaxies. Our method is based on two-dimensional Fokker-Planck equations in the energy and angular momentum space, and allows the evolution of multiple mass components, including stars and compact objects. The code demonstrates remarkable flexibility in incorporating additional complex dynamics. By employing a weighting method, we effectively enhance the statistical accuracy of rare particle results. In this initial publication, we present the fundamental version of our method, focusing on two-body relaxations and loss cone effects. Through comparisons with previous studies, we establish consistent outcomes in terms of relaxation processes, energy and angular momentum distributions, density profiles, and loss cone consumption rates. We consistently observe the development of tangential anisotropy within the cluster, while the outer regions tend to retain near-isotropic characteristics. GNC holds great promise for exploring a wide range of intriguing phenomena within galactic nuclei, including relativistic stellar dynamics, providing detailed and insightful outcomes.