No Arabic abstract
We derive the kinetic equation that describes the secular evolution of a large set of particles orbiting a dominant massive object, such as stars bound to a supermassive black hole or a proto-planetary debris disc encircling a star. Because the particles move in a quasi-Keplerian potential, their orbits can be approximated by ellipses whose orientations remain fixed over many dynamical times. The kinetic equation is obtained by simply averaging the BBGKY equations over the fast angle that describes motion along these ellipses. This so-called Balescu-Lenard equation describes self-consistently the long-term evolution of the distribution of quasi-Keplerian orbits around the central object: it models the diffusion and drift of their actions, induced through their mutual resonant interaction. Hence, it is the master equation that describes the secular effects of resonant relaxation. We show how it captures the phenonema of mass segregation and of the relativistic Schwarzschild barrier recently discovered in $N$-body simulations.
Our current understanding of the stellar initial mass function and massive star evolution suggests that young globular clusters may have formed hundreds to thousands of stellar-mass black holes, the remnants of stars with initial masses from $sim 20 - 100, M_odot$. Birth kicks from supernova explosions may eject some black holes from their birth clusters, but most should be retained. Using a Monte Carlo method we investigate the long-term dynamical evolution of globular clusters containing large numbers of stellar black holes. We describe numerical results for 42 models, covering a range of realistic initial conditions, including up to $1.6times10^6$ stars. In almost all models we find that significant numbers of black holes (up to $sim10^3$) are retained all the way to the present. This is in contrast to previous theoretical expectations that most black holes should be ejected dynamically within a few Gyr. The main reason for this difference is that core collapse driven by black holes (through the Spitzer mass segregation instability) is easily reverted through three-body processes, and involves only a small number of the most massive black holes, while lower-mass black holes remain well-mixed with ordinary stars far from the central cusp. Thus the rapid segregation of stellar black holes does not lead to a long-term physical separation of most black holes into a dynamically decoupled inner core, as often assumed previously. Combined with the recent detections of several black hole X-ray binary candidates in Galactic globular clusters, our results suggest that stellar black holes could still be present in large numbers in many globular clusters today, and that they may play a significant role in shaping the long-term dynamical evolution and the present-day dynamical structure of many clusters.
Chandrasekhars most important contribution to stellar dynamics was the concept of dynamical friction. I briefly review that work, then discuss some implications of Chandrasekhars theory of gravitational encounters for motion in galactic nuclei.
The angular momentum evolution of stars close to massive black holes (MBHs) is driven by secular torques. In contrast to two-body relaxation, where interactions between stars are incoherent, the resulting resonant relaxation (RR) process is characterized by coherence times of hundreds of orbital periods. In this paper, we show that all the statistical properties of RR can be reproduced in an autoregressive moving average (ARMA) model. We use the ARMA model, calibrated with extensive N-body simulations, to analyze the long-term evolution of stellar systems around MBHs with Monte Carlo simulations. We show that for a single-mass system in steady-state, a depression is carved out near an MBH as a result of tidal disruptions. Using Galactic center parameters, the extent of the depression is about 0.1 pc, of similar order to but less than the size of the observed hole in the distribution of bright late-type stars. We also find that the velocity vectors of stars around an MBH are locally not isotropic. In a second application, we evolve the highly eccentric orbits that result from the tidal disruption of binary stars, which are considered to be plausible precursors of the S-stars in the Galactic center. We find that RR predicts more highly eccentric (e > 0.9) S-star orbits than have been observed to date.
The spin angular momentum S of a supermassive black hole (SBH) precesses due to torques from orbiting stars, and the stellar orbits precess due to dragging of inertial frames by the spinning hole. We solve the coupled post-Newtonian equations describing the joint evolution of S and the stellar angular momenta Lj, j = 1...N in spherical, rotating nuclear star clusters. In the absence of gravitational interactions between the stars, two evolutionary modes are found: (1) nearly uniform precession of S about the total angular momentum vector of the system; (2) damped precession, leading, in less than one precessional period, to alignment of S with the angular momentum of the rotating cluster. Beyond a certain distance from the SBH, the time scale for angular momentum changes due to gravitational encounters between the stars is shorter than spin-orbit precession times. We present a model, based on the Ornstein-Uhlenbeck equation, for the stochastic evolution of star clusters due to gravitational encounters and use it to evaluate the evolution of S in nuclei where changes in the Lj are due to frame dragging close to the SBH and to encounters farther out. Long-term evolution in this case is well described as uniform precession of the SBH about the clusters rotational axis, with an increasingly important stochastic contribution when SBH masses are small. Spin precessional periods are predicted to be strongly dependent on nuclear properties, but typical values are 10-100 Myr for low-mass SBHs in dense nuclei, 100 Myr - 10 Gyr for intermediate mass SBHs, and > 10 Gyr for the most massive SBHs. We compare the evolution of SBH spins in stellar nuclei to the case of torquing by an inclined, gaseous accretion disk.
Interaction of a binary supermassive black hole with stars in a galactic nucleus can result in changes to all the elements of the binarys orbit, including the angles that define its orientation. If the nucleus is rotating, the orientation changes can be large, causing large changes in the binarys orbital eccentricity as well. We present a general treatment of this problem based on the Fokker-Planck equation for f, defined as the probability distribution for the binarys orbital elements. First- and second-order diffusion coefficients are derived for the orbital elements of the binary using numerical scattering experiments, and analytic approximations are presented for some of these coefficients. Solutions of the Fokker-Planck equation are then derived under various assumptions about the initial rotational state of the nucleus and the binary hardening rate. We find that the evolution of the orbital elements can become qualitatively different when we introduce nuclear rotation: 1) the orientation of the binarys orbit evolves toward alignment with the plane of rotation of the nucleus; 2) binary orbital eccentricity decreases for aligned binaries and increases for counter-aligned ones. We find that the diffusive (random-walk) component of a binarys evolution is small in nuclei with non-negligible rotation, and we derive the time-evolution equations for the semimajor axis, eccentricity and inclination in that approximation. The aforementioned effects could influence gravitational wave production as well as the relative orientation of host galaxies and radio jets.