Rotational Brownian motion (astronomy)
Consider a binary that consists of two massive objects (stars, black holes etc.) and that is embedded in a stellar system containing a large number of stars. Let and be the masses of the two components of the binary whose total mass is . A field star that approaches the binary with impact parameter and velocity passes a distance from the binary, where
the latter expression is valid in the limit that gravitational focusing dominates the encounter rate. The rate of encounters with stars that interact strongly with the binary, i.e. that satisfy , is approximately where and are the number density and velocity dispersion of the field stars and is the semi-major axis of the binary.
As it passes near the binary, the field star experiences a change in velocity of order
where is the relative velocity of the two stars in the binary. The change in the field star's specific angular momentum with respect to the binary, , is then Δl ≈ a Vbin. Conservation of angular momentum implies that the binary's angular momentum changes by Δlbin ≈ -(m/μ12)Δl where m is the mass of a field star and μ12 is the binary reduced mass. Changes in the magnitude of lbin correspond to changes in the binary's orbital eccentricity via the relation e = 1 - lb2/GM12μ12a. Changes in the direction of lbin correspond to changes in the orientation of the binary, leading to rotational diffusion. The rotational diffusion coefficient is
where ρ = mn is the mass density of field stars.
Let F(θ,t) be the probability that the rotation axis of the binary is oriented at angle θ at time t. The evolution equation for F is 
If <Δξ2>, a, ρ and σ are constant in time, this becomes
where μ = cos θ and τ is the time in units of the relaxation time trel, where
The solution to this equation states that the expectation value of μ decays with time as
Hence, trel is the time constant for the binary's orientation to be randomized by torques from field stars.
Rotational Brownian motion was first discussed in the context of binary supermassive black holes at the centers of galaxies. Perturbations from passing stars can alter the orbital plane of such a binary, which in turn alters the direction of the spin axis of the single black hole that forms when the two coalesce.
Rotational Brownian motion is often observed in N-body simulations of galaxies containing binary black holes. The massive binary sinks to the center of the galaxy via dynamical friction where it interacts with passing stars. The same gravitational perturbations that induce a random walk in the orientation of the binary, also cause the binary to shrink, via the gravitational slingshot. It can be shown that the rms change in the binary's orientation, from the time the binary forms until the two black holes collide, is roughly
In a real galaxy, the two black holes would eventually coalesce due to emission of gravitational waves. The spin axis of the coalesced hole will be aligned with the angular momentum axis of the orbit of the pre-existing binary. Hence, a mechanism like rotational Brownian motion that affects the orbits of binary black holes can also affect the distribution of black hole spins. This may explain in part why the spin axes of supermassive black holes appear to be randomly aligned with respect to their host galaxies.
- Debye, P. (1929). Polar Molecules. Dover.
- Merritt, D. (2002), Rotational Brownian Motion of a Massive Binary, The Astrophysical Journal, 568, 998-1003.
- Löckmann, U. and Baumgardt, H. (2008), Tracing intermediate-mass black holes in the Galactic Centre, Monthly Notices of the Royal Astronomical Society, 384, 323-330.
- Matsubayashi, T., Makino, J. and Ebisuzaki, T. (2007), Evolution of an IMBH in the Galactic Nucleus with a Massive Central Black Hole, The Astrophysical Journal, 656, 879-896
- Kinney, A. et al. (2000), Jet Directions in Seyfert Galaxies, The Astrophysical Journal, 537, 152-177
- Gravitational Scattering Review article on the dynamics of encounters between binaries and single stars.