Rotational Brownian motion (astronomy)

In astronomy, rotational Brownian motion is the random walk in orientation of a binary star's orbital plane, induced by gravitational perturbations from passing stars.
Theory.
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 formula_1 and formula_2 be the masses of the two components of the binary whose total mass is formula_3. A field star that approaches the binary with impact parameter formula_4 and velocity formula_5 passes a distance formula_6 from the binary, where
formula_7
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 formula_8, is approximately formula_9 where formula_10 and formula_11 are the number density and velocity dispersion of the field stars and formula_12 is the semi-major axis of the binary.
As it passes near the binary, the field star experiences a change in velocity of order 
formula_13, 
where formula_14 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, formula_15, is then Δ"l" ≈ "a" "V"bin. Conservation of angular momentum implies that the binary's angular momentum changes by Δ"l"bin ≈ -(m/μ12)Δ"l" where "m" is the mass of a field star and μ12 is the binary reduced mass. Changes in the magnitude of "l"bin correspond to changes in the binary's orbital eccentricity via the relation "e" = 1 - "l"b2/"GM"12μ12"a". Changes in the direction of "l"bin correspond to changes in the orientation of the binary, leading to rotational diffusion. The rotational diffusion coefficient is
formula_16
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 
formula_17
If <Δξ2>, "a", ρ and σ are constant in time, this becomes
formula_18
where μ = cos θ and τ is the time in units of the relaxation time "t"rel, where
formula_19
The solution to this equation states that the expectation value of μ decays with time as
formula_20
Hence, "t"rel is the time constant for the binary's orientation to be randomized by torques from field stars.
Applications.
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
formula_21
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.