We present fifth order Runge-Kutta-Nystr\"om methods, where we allow the timestep coefficients to assume complex values. Among the methods with complex timesteps, we focus on the ones with the coefficients that have positive real parts. This property makes them suitable for problems where a negative coefficient is not acceptable. In addition, the leading order terms in the error expansion of these methods are purely imaginary, effectively increasing the order of the methods by one for real variables.

In this paper, we propose a method to study the nature of resonant relaxation in near-Keplerian systems. Our technique is based on measuring the fractal dimension of the angular momentum trails and we use it to analyze the outcome of N-body simulations. With our method, we can reliably determine the timescale for resonant relaxation, as well as the rate of change of angular momentum in this regime. We find that growth of angular momentum is more rapid than random walk, but slower than linear growth. We also determine the presence of long term correlations, arising from the bounds on angular momentum growth. We develop a toy model that reproduces all essential properties of angular momentum evolution.

The orbits of stars close to a massive black hole are nearly Keplerian ellipses. Such orbits exert long term torques on each other, which lead to an enhanced angular momentum relaxation known as resonant relaxation. Under certain conditions, this process can modify the angular momentum distribution and affect the interaction rates of the stars with the massive black hole more efficiently than non-resonant relaxation. The torque on an orbit exerted by the cluster depends on the eccentricity of the orbit. In this paper, we calculate this dependence and determine the resonant relaxation timescale as a function of eccentricity. In particular, we show that the component of the torque that changes the magnitude of the angular momentum is linearly proportional to eccentricity, so resonant relaxation is much more efficient on eccentric orbits than on circular orbits.

We study the evolution of the mass function in young and dense star clusters by means of direct N-body simulations. Our main aim is to explain the recent observations of the relatively flat mass function observed near the centre of the Arches star cluster. In this region, the power law index of the mass function for stars more massive than about 5-6 solar mass, is larger than the Salpeter value by about unity; whereas further out, and for the lower mass stars, the mass function resembles the Salpeter distribution. We show that the peculiarities in the Arches mass function can be explained satisfactorily without primordial mass segregation. We draw two conclusions from our simulations: 1) The Arches initial mass function is consistent with a Salpeter slope down to ~1 solar mass, 2) The cluster is about half way towards core collapse. The cores of other star clusters with characteristics similar to those of the Arches are expected to show similar flattening in the mass functions for the high mass (>5 solar mass) stars.

Recent theoretical work has solidified the viability of the collisional runaway scenario in young dense star clusters for the formation of very massive stars (VMSs), which may be precursors to intermediate-mass black holes (IMBHs). We present first results from a numerical study of the collisional runaway process in dense star clusters containing primordial binaries. Stellar collisions during binary scattering encounters offer an alternate channel for runaway growth, somewhat independent of direct collisions between single stars. We find that clusters with binary fractions >~10% yield two VMSs via collisional runaways, presenting the exotic possibility of forming IMBH--IMBH binaries in star clusters. We discuss the implications for gravitational wave observations, and the impact on cluster structure.

We present a new study of the collisional runaway scenario to form an intermediate-mass black hole (IMBH, MBH > 100 Msun) at the centre of a young, compact stellar cluster. The first phase is the formation of a very dense central core of massive stars (Mstar =~ 30-120 Msun) through mass segregation and gravothermal collapse. Previous work established the conditions for this to happen before the massive stars evolve off the main sequence (MS). In this and a companion paper, we investigate the next stage by implementing direct collisions between stars. Using a Monte Carlo stellar dynamics code, we follow the core collapse and subsequent collisional phase in more than 100 models with varying cluster mass, size, and initial concentration. Collisions are treated either as ideal, ``sticky-sphere'' mergers or using realistic prescriptions derived from 3-D hydrodynamics computations. In all cases for which the core collapse happens in less than the MS lifetime of massive stars (~3 Myr), we obtain the growth of a single very massive star (VMS, Mstar =~ 400-4000 Msun) through a runaway sequence of mergers. Mass loss from collisions, even for velocity dispersions as high as sigma1D ~ 1000 km/s, does not prevent the runaway. The region of cluster parameter space leading to runaway is even more extended than predicted in previous work because, in clusters with sigma1D > 300 km/s, collisions accelerate (and, in extreme cases, drive) core collapse. Although the VMS grows rapidly to > 1000 Msun in models exhibiting runaway, we cannot predict accurately its final mass. This is because the termination of the runaway process must eventually be determined by a complex interplay between stellar dynamics, hydrodynamics, and the stellar evolution of the VMS. [abridged]

We present results from dynamical Monte Carlo simulations of dense star clusters near the Galactic center. These clusters sink toward the center of the Galaxy by dynamical friction. During their inspiral, they may undergo core collapse and form an intermediate-mass black hole through runaway collisions. Such a cluster can then reach within a parsec of the Galactic center before it completely disrupts, releasing many young stars in this region. This scenario provides a natural explanation for the presence of the young stars observed near the Galactic center. Here we determine the initial conditions for this scenario to work and we derive the mass distribution of cluster stars as a function of distance from the Galactic center. For reasonable initial conditions, we find that clusters massive enough for rapid inspiral would include more massive stars (m >~ 30 Msun) than currently observed in the inspiral region. We point out several possible explanations for this apparent discrepancy.

Using Monte Carlo codes, we follow the collisional evolution of clusters in a variety of scenarios. We consider the conditions under which a cluster of main sequence stars may undergo rapid core collapse due to mass segregation, thus entering a phase of runaway collisions, forming a very massive star (VMS, M >= 1000 Msun) through repeated collisions between single stars. Although collisional mass loss is accounted for realistically, we find that a VMS forms even in proto-galactic nuclei models with a high velocity dispersion (many 100 km/s). Such a VMS may be a progenitor for an intermediate-mass black hole (M >= 100 Msun). In contrast, in galactic nuclei hosting a central massive black hole, collisions are found to be disruptive. The stars which are subject to collisions are progressively ground down by high-velocity collisions and a merger sequence appears impossible.

We consider spherical stellar clusters with a broad mass function and a relaxation time short enough so that the segregation of massive stars toward the centre occurs before they have time to evolve off the main sequence. The relaxational and collisional dynamics of model clusters is followed with a high-resolution Monte Carlo code. Stellar collisions are treated in a realistic way, through the use of a outcome of a very large set of SPH simulations. We find that, even in proto-galactic nuclei models with high velocity dispersions, run-away growth of a very massive star (VMS, M>100 M_sun) occurs in all cases when the core collapse time is shorter than the MS life time of massive stars, i.e. 3 Myrs. The VMS is a likely progenitor for an intermediate-mass or massive black hole (IMBH/MBH).

We study the early dynamical evolution of young, dense star clusters using Monte Carlo simulations for systems with up to N~10^7 stars. Rapid mass segregation of massive main-sequence stars and the development of the Spitzer instability can drive these systems to core collapse in a small fraction of the initial half-mass relaxation time. If the core collapse time is less than the lifetime of the massive stars, all stars in the collapsing core may then undergo a runaway collision process leading to the formation of a massive black hole. Here we study in detail the first step in this process, up to the occurrence of core collapse. We have performed about 100 simulations for clusters with a wide variety of initial conditions, varying systematically the cluster density profile, stellar IMF, and number of stars. We also considered the effects of initial mass segregation and stellar evolution mass loss. Our results show that, for clusters with a moderate initial central concentration and any realistic IMF, the ratio of core collapse time to initial half-mass relaxation time is typically ~0.1, in agreement with the value previously found by direct N-body simulations for much smaller systems. Models with even higher central concentration initially, or with initial mass segregation (from star formation) have even shorter core-collapse times. Remarkably, we find that, for all realistic initial conditions, the mass of the collapsing core is always close to ~10^-3 of the total cluster mass, very similar to the observed correlation between central black hole mass and total cluster mass in a variety of environments. We discuss the implications of our results for the formation of intermediate-mass black holes in globular clusters and super star clusters, ultraluminous X-ray sources, and seed black holes in proto-galactic nuclei.

We review possible dynamical formation processes for central massive black holes in dense star clusters. We focus on the early dynamical evolution of young clusters containing a few thousand to a few million stars. One natural formation path for a central seed black hole in these systems involves the development of the Spitzer instability, through which the most massive stars can drive the cluster to core collapse in a very short time. The sudden increase in the core density then leads to a runaway collision process and the formation of a very massive merger remnant, which must then collapse to a black hole. Alternatively, if the most massive stars end their lives before core collapse, a central cluster of stellar-mass black holes is formed. This cluster will likely evaporate before reaching the highly relativistic state necessary to drive a runaway merger process through gravitational radiation, thereby avoiding the formation of a central massive black hole. We summarize the conditions under which these different paths will be followed, and present the results of recent numerical simulations demonstrating the process of rapid core collapse and runaway collisions between massive stars.