Vector Resonant Relaxation Research Articles

Constraining intermediate-mass black holes from the stellar disc of SgrA*

ABSTRACT Stars evolving around a supermassive black hole see their orbital orientations diffuse efficiently, a process called ‘vector resonant relaxation’. In particular, stars within the same disc, i.e. neighbours in orientations, will slowly diffuse away from one another through this stochastic process. We use jointly (i) detailed kinetic predictions for the efficiency of this dilution and (ii) the recent observation of a stellar disc around SgrA*, the supermassive black hole at the centre of the Milky Way, to constrain SgrA*’s unobserved stellar cluster. Notably, we investigate quantitatively the impact of a population of intermediate-mass black holes on the survivability of the stellar disc.

Black hole discs and spheres in galactic nuclei – exploring the landscape of vector resonant relaxation equilibria

ABSTRACT Vector resonant relaxation (VRR) is known to be the fastest gravitational process that shapes the geometry of stellar orbits in nuclear star clusters. This leads to the realignment of the orbital planes on the corresponding VRR time-scale tVRR of a few million years, while the eccentricity e and semimajor axis a of the individual orbits are approximately conserved. The distribution of orbital inclinations reaches an internal equilibrium characterized by two conserved quantities, the total potential energy among stellar orbits, Etot, and the total angular momentum, Ltot. On time-scales longer than tVRR, the eccentricities and semimajor axes change slowly, and the distribution of orbital inclinations are expected to evolve through a series of VRR equilibria. Using a Monte Carlo Markov Chain method, we determine the equilibrium distribution of orbital inclinations in the microcanonical ensemble with fixed Etot and Ltot for isolated nuclear star clusters with a power-law distribution of a, e, and m, where m is the stellar mass. We explore the possible equilibria for nine representative Etot–Ltot pairs that cover the possible parameter space. For all cases, the equilibria show anisotropic mass segregation, where the distribution of more massive objects is more flattened than that for lighter objects. Given that stellar black holes are more massive than the average main-sequence stars, these findings suggest that black holes reside in disc-like structures within nuclear star clusters for a wide range of initial conditions.

A numerical study of stellar discs in galactic nuclei

ABSTRACT We explore the dynamics of stellar discs in the close vicinity of a supermassive black hole (SMBH) by means of direct N-body simulations. We show that an isolated nuclear stellar disc exhibits anisotropic mass segregation meaning that massive stars settle to lower orbital inclinations and more circular orbits than the light stars. However, in systems in which the stellar disc is embedded in a much more massive isotropic stellar cluster, anisotropic mass segregation tends to be suppressed. In both cases, an initially thin stellar disc becomes thicker, especially in the inner parts due to the fluctuating anisotropy in the spherical component. We find that vector resonant relaxation is quenched in the disc by nodal precession, but it is still the most efficient relaxation process around SMBHs of mass $10^6\, {\rm M}_{\odot }$ and above. Two-body relaxation may dominate for less massive SMBHs found in dwarf galaxies. Stellar discs embedded in massive isotropic stellar clusters ultimately tend to become isotropic on the local two-body relaxation time-scale. Our simulations show that the dynamics of young stars at the centre of the Milky Way is mostly driven by vector resonant relaxation leading to an anticorrelation between the scatter of orbital inclinations and distance from the SMBH. If the S-stars formed in a disc less than 10 Myr ago, they may coexist with a cusp of stellar mass black holes or an intermediate mass black hole with mass up to $1000\, {\rm M}_{\odot }$ to reproduce the observed scatter of angular momenta.

Orbital alignment and mass segregation in galactic nuclei via vector resonant relaxation

ABSTRACT Supermassive black holes dominate the gravitational potential in galactic nuclei. In these dense environments, stars follow nearly Keplerian orbits and see their orbital planes relax through the potential fluctuations generated by the stellar cluster itself. For typical astrophysical galactic nuclei, the most likely outcome of this vector resonant relaxation is that the orbital planes of the most massive stars spontaneously self-align within a narrow disc. We present a maximum entropy method to systematically determine this long-term distribution of orientations and use it for a wide range of stellar orbital parameters and initial conditions. The heaviest stellar objects are found to live within a thin equatorial disc. The thickness of this disk depends on the stars’ initial mass function, and on the geometry of the initial cluster. This work highlights a possible (indirect) novel method to constrain the distribution of intermediate mass black holes in galactic nuclei.

Young stellar cluster dilution near supermassive black holes: the impact of vector resonant relaxation on neighbour separation

ABSTRACT We investigate the rate of orbital orientation dilution of young stellar clusters in the vicinity of supermassive black holes. Within the framework of vector resonant relaxation, we predict the time evolution of the two-point correlation function of the stellar orbital plane orientations as a function of their initial angular separation and diversity in orbital parameters (semimajor axis, eccentricity). As expected, the larger the spread in initial orientations and orbital parameters, the more efficient the dilution of a given set of co-eval stars, with a characteristic time-scale set up by the coherence time of the background potential fluctuations. A Markovian prescription that matches numerical simulations allows us to efficiently probe the underlying kinematic properties of the unresolved nucleus when requesting consistency with a given dilution efficiency, imposed by the observed stellar disc within the 1 arcsec of Sgr A*. As a proof of concept, we compute maps of constant dilution times as a function of the semimajor axis cusp index and fraction of intermediate-mass black holes in the old background stellar cluster. This computation suggests that vector resonant relaxation should prove useful in this context since it impacts orientations on time-scales comparable to the stars’ age.

Statistical mechanics of gravitational systems with regular orbits: rigid body model of vector resonant relaxation

I consider a self-gravitating, N-body system assuming that the N constituents follow regular orbits about the center of mass of the cluster, where a central massive object may be present. I calculate the average over a characteristic timescale of the full, N-body Hamiltonian including all kinetic and potential energy terms. The resulting effective system allows for the identification of the orbital planes with N rigid, disk-shaped tops, that can rotate about their fixed common centre and are subject to mutual gravitational torques. The time-averaging imposes boundaries on the canonical generalized momenta of the resulting canonical phase space. I investigate the statistical mechanics induced by the effective Hamiltonian on this bounded phase space and calculate the thermal equilibrium states. These are a result of the relaxation of spins’ directions, identified with orbital planes’ orientations, which is called vector resonant relaxation. I calculate the dependence of spins’ angular velocity dispersion on temperature and calculate the velocity distribution functions. I argue that the range of validity of the gravitational phase transitions, identified in the special case of zero kinetic term by Roupas, Kocsis and Tremaine, is expanded to non-zero values of the ratio of masses between the cluster of N-bodies and the central massive object. The relevance with astrophysics is discussed focusing on stellar clusters. The same analysis performed on an unbounded phase space accounts for continuous rigid tops.

Vector Resonant Relaxation of Stars around a Massive Black Hole

Abstract In the vicinity of a massive black hole, stars move on precessing Keplerian orbits. The mutual stochastic gravitational torques between the stellar orbits drive a rapid reorientation of their orbital planes, through a process called vector resonant relaxation. We derive, from first principles, the correlation of the potential fluctuations in such a system, and the statistical properties of random walks undergone by the stellar orbital orientations. We compare this new analytical approach with numerical simulations. We also provide a simple scheme to generate the random walk of a test star’s orbital orientation using a stochastic equation of motion. We finally present quantitative estimations of this process for a nuclear stellar cluster such as that of the Milky Way.

Resonant Relaxation in Globular Clusters

Abstract Resonant relaxation has been discussed as an efficient process that changes the angular momenta of stars orbiting around a central supermassive black hole due to the fluctuating gravitational field of the stellar cluster. Other spherical stellar systems, such as globular clusters, exhibit a restricted form of this effect where enhanced relaxation rate only occurs in the directions of the angular momentum vectors, but not in their magnitudes; this is called vector resonant relaxation (VRR). To explore this effect, we performed a large set of direct N-body simulations, with up to 512k particles (where k =1024) and ∼500 dynamical times. Contrasting these simulations, which naturally include the collective effects, with Spitzer-style Monte Carlo simulations, which by design only exhibit two-body relaxation, we show that the temporal behavior of the angular momentum vectors in N-body simulations cannot be explained by two-body relaxation alone. VRR operates efficiently in globular clusters with N > 104. The fact that VRR operates in globular clusters may open a way to use powerful tools in statistical physics for their description. In particular, since the distribution of orbital planes relaxes much more rapidly than the distribution of the magnitude of angular momentum and the radial action, the relaxation process reaches an internal statistical equilibrium in the corresponding part of phase space while the whole cluster is generally out of equilibrium, in a state of quenched disorder. We point out the need to include effects of VRR in Monte Carlo simulations of globular clusters.

The Impact of Vector Resonant Relaxation on the Evolution of Binaries near a Massive Black Hole: Implications for Gravitational-wave Sources

Abstract Binaries within the sphere of influence of a massive black hole (MBH) in galactic nuclei are susceptible to the Lidov–Kozai (LK) mechanism, which can drive orbits to high eccentricities and trigger strong interactions within the binary such as the emission of gravitational waves (GWs) and mergers of compact objects. These events are potential sources for GW detectors such as Advanced LIGO and VIRGO. The LK mechanism is only effective if the binary is highly inclined with respect to its orbit around the MBH (within a few degrees of 90°), implying low rates. However, close to an MBH, torques from the stellar cluster give rise to the process of vector resonant relaxation (VRR). VRR can bring a low-inclination binary into an “active” LK regime in which high eccentricities and strong interactions are triggered in the binary. Here, we study the coupled LK–VRR dynamics, with implications for LIGO and VIRGO GW sources. We carry out Monte Carlo simulations and find that the merger fraction enhancement due to LK–VRR dynamics is up to a factor of ∼10 for the lower end of assumed MBH masses ( ) and decreases sharply with increasing M •. We find that, even in our most optimistic scenario, the baseline BH–BH merger rate is small, and the enhancement by LK–VRR coupling is not large enough to increase the rate to well above the LIGO/VIRGO lower limit, . For the Galactic Center, the LK–VRR-enhanced rate is ∼100 times lower than the LIGO/VIRGO limit, and for , the rate barely reaches .

Black Hole Disks in Galactic Nuclei.

Gravitational torques among objects orbiting a supermassive black hole drive the rapid reorientation of orbital planes in nuclear star clusters (NSCs), a process known as vector resonant relaxation. In this Letter, we determine the statistical equilibrium of systems with a distribution of masses, semimajor axes, and eccentricities. We average the interaction over the apsidal precession time and construct a MonteCarlo Markov chain method to sample the microcanonical ensemble of the NSC. We examine the case of NSCs formed by 16 episodes of star formation or globular cluster infall. We find that the massive stars and stellar mass black holes form a warped disk, while low mass stars resemble a spherical distribution with a possible net rotation. This explains the origin of the clockwise disk in the Galactic center and predicts a population of black holes (BHs) embedded within this structure. The rate of mergers among massive stars, tidal disruption events of massive stars by BHs, and BH-BH mergers are highly increased in such disks. The first two may explain the origin of the observed G1 and G2 clouds, the latter may be important for gravitational wave detections with LIGO and VIRGO. More generally, black holes are expected to settle in disks in all dense spherical stellar systems assembled by mergers of smaller systems including globular clusters.

Isotropic–Nematic Phase Transitions in Gravitational Systems. II. Higher Order Multipoles

Abstract The gravitational interaction among bodies orbiting in a spherical potential leads to the rapid relaxation of the orbital planes’ distribution, a process called vector resonant relaxation. We examine the statistical equilibrium of this process for a system of bodies with similar semimajor axes and eccentricities. We extend the previous model of Roupas et al. by accounting for the multipole moments beyond the quadrupole, which dominate the interaction for radially overlapping orbits. Nevertheless, we find no qualitative differences between the behavior of the system with respect to the model restricted to the quadrupole interaction. The equilibrium distribution resembles a counterrotating disk at low temperature and a spherical structure at high temperature. The system exhibits a first-order phase transition between the disk and the spherical phase in the canonical ensemble if the total angular momentum is below a critical value. We find that the phase transition erases the high-order multipoles, i.e., small-scale structure in angular momentum space, most efficiently. The system admits a maximum entropy and a maximum energy, which lead to the existence of negative temperature equilibria.

Isotropic–Nematic Phase Transitions in Gravitational Systems

Abstract We examine dense self-gravitating stellar systems dominated by a central potential, such as nuclear star clusters hosting a central supermassive black hole. Different dynamical properties of these systems evolve on vastly different timescales. In particular, the orbital-plane orientations are typically driven into internal thermodynamic equilibrium by vector resonant relaxation before the orbital eccentricities or semimajor axes relax. We show that the statistical mechanics of such systems exhibit a striking resemblance to liquid crystals, with analogous ordered-nematic and disordered-isotropic phases. The ordered phase consists of bodies orbiting in a disk in both directions, with the disk thickness depending on temperature, while the disordered phase corresponds to a nearly isotropic distribution of the orbit normals. We show that below a critical value of the total angular momentum, the system undergoes a first-order phase transition between the ordered and disordered phases. At a critical point, the phase transition becomes second order, while for higher angular momenta there is a smooth crossover. We also find metastable equilibria containing two identical disks with mutual inclinations between 90° and 180°.

A numerical study of vector resonant relaxation

Stars bound to a supermassive black hole interact gravitationally. Persistent torques acting between stellar orbits lead to the rapid resonant relaxation of the orbital orientation vectors ("vector" resonant relaxation) and slower relaxation of the eccentricities ("scalar" resonant relaxation), both at rates much faster than two-body or non-resonant relaxation. We describe a new parallel symplectic integrator, N-ring, which follows the dynamical evolution of a cluster of N stars through vector resonant relaxation, by averaging the pairwise interactions over the orbital period and periapsis-precession timescale. We use N-ring to follow the evolution of clusters containing over 10^4 stars for tens of relaxation times. Among other results, we find that the evolution is dominated by torques among stars with radially overlapping orbits, and that resonant relaxation can be modelled as a random walk of the orbit normals on the sphere, with angular step size ranging from 0.5-1 radian. The relaxation rate in a cluster with a fixed number of stars is proportional to the RMS mass of the stars. The RMS torque generated by the cluster stars is reduced below the torque between Kepler orbits due to apsidal precession and declines weakly with the eccentricity of the perturbed orbit. However since the angular momentum of an orbit also decreases with eccentricity, the relaxation rate is approximately eccentricity-independent for e<0.7 and grows rapidly with eccentricity for e>0.8. We quantify the relaxation using the autocorrelation function of the spherical multipole moments; this decays exponentially and the e-folding time may be identified with the vector resonant relaxation timescale.

Resonant relaxation and the warp of the stellar disc in the Galactic Centre

Observations of the spatial distribution and kinematics of young stars in the Galactic centre can be interpreted as showing that the stars occupy one, or possibly two, discs of radii ~0.05-0.5 pc. The most prominent (`clockwise') disc exhibits a strong warp: the normals to the mean orbital planes in the inner and outer third of the disc differ by ~60 deg. Using an analytical model based on Laplace-Lagrange theory, we show that such warps arise naturally and inevitably through vector resonant relaxation between the disc and the surrounding old stellar cluster.