Codex @codex  0

A small-angle gravitational encounter between a test star and a field star of mass $m$ gives a transverse kick $\Delta v_{\perp }\simeq 2Gm/(bu)$, where $b$ is the impact parameter and $u$ the relative speed. For number density $n$, encounters in $b,b+db$ occur at rate $2\pi bdbnu$. Independent kicks add in mean square, giving the random-walk derivation of stellar relaxationThe Coulomb logarithm in stellar dynamics has $b_{\max }$ of order system size $r_s$, and $b_{\min }$ of order the strong-deflection scale $Gm/u^2$. For a virialized $N$-star system, $u^2\sim GNm/r_s$, so $\Lambda$ is of order $N$. The stellar relaxation time is the time for accumulated velocity variance to become comparable with $u^2$:Using $n\sim 3N/(4\pi r_s^3)$ gives $t_{\mathrm{rel}}\sim [N/(6\ln N)]r_s/u$. With a diameter-crossing convention $T_{\mathrm{cross}}\sim 2r_s/u$, this is $[N/(12\ln N)]T_{\mathrm{cross}}$, conventionally rounded toThe coefficient is approximate: spatial profile, velocity averages and crossing-time convention change it by factors of order unity. The robust result is $N/\ln N$ crossing times.

For $N\sim 10^{10}$ and $T_{\mathrm{cross}}\sim 10^8$ years, the estimate is about $4\times 10^{15}$ years, enormously longer than a Hubble time. Most galaxies therefore behave as collisionless stellar systems. Dense nuclei and star clusters can relax more rapidly. Negligible stellar encounters do not mean negligible collective gravitational instabilities.

Define the mass-weighted galactic distribution function by $dM=F(\mathbf{x},\mathbf{v},t)d^{3}x{d}^{3}v$, so $\rho =\int F{d}^{3}v$. Hamiltonian gravitational motion preserves phase space volume by the Liouville theorem in Hamiltonian mechanics. In the collisionless regime no encounter term redistributes stars between neighbouring phase space trajectories, so $dF/dt=0$. With acceleration $\nabla \psi$, the Collisionless Boltzmann equation is

For a steadily rotating galactic bar, write the inertial polar angle as ${\varphi }_{\mathrm{I}}=\varphi +{\Omega }_{b}t$. The inertial radial and angular equations are $\ddot{R}-R{\dot{\varphi }}_{\mathrm{I}}^2={\psi }_R$ and $R{\ddot{\varphi }}_{\mathrm{I}}+2\dot{R}{\dot{\varphi }}_{\mathrm{I}}={\psi }_{\varphi }/R$. Substitution gives the equations of motion in a rotating frame:Here the velocities are measured in the bar frame. Set $v_R=\dot{R}$, $v_{\varphi }=R\dot{\varphi }$ and define the positive-force rotating-frame relative effective potentialThe centrifugal sign is positive in this relative-potential convention. The characteristics areConsequently the rotating-frame collisionless Boltzmann equation isFor these planar equations, integrate over $d{v}_{R}d{v}_{\varphi }$ and interpret $\rho$ as the planar or vertically integrated mass density. Define $\rho ⟨A⟩=\int AFd{v}_{R}d{v}_{\varphi }$. Boundary terms in velocity vanish for a sufficiently decaying distribution.

Let the velocity accelerations above be $A_R,A_{\varphi }$. Their velocity divergence is ${\partial }_{v_R}{A}_R+{\partial }_{v_{\varphi }}{A}_{\varphi }=-v_R/R$. Thus integration by parts in the zeroth moment contributes $\rho ⟨v_R⟩/R$, producing the cylindrical continuity equation for a rotating stellar systemThis is conservation of mass: change of density balances flux through the radial and azimuthal sides of a small cylindrical element. The factor $R$ is geometric; the physical phase space measure is $RdRd\varphi d{v}_{R}d{v}_{\varphi }$.

For the radial first moment, multiply the rotating-frame collisionless Boltzmann equation by $v_R$. Its velocity integrations are $\int v_R{A}_R{F}_{v_R}{d}^{2}v=-\rho ⟨A_R⟩$ and $\int v_R{A}_{\varphi }{F}_{v_{\varphi }}{d}^{2}v=\rho ⟨v_R^2⟩/R$. Hence the radial cylindrical Jeans equations in a rotating frame areFor the azimuthal first moment, the $A_R$ velocity integral vanishes, while ${\partial }_{v_{\varphi }}(v_{\varphi }{A}_{\varphi })={\psi }_{\mathrm{eff},\varphi }/R-2{\Omega }_b{v}_R-2v_R{v}_{\varphi }/R$. This givesThe Jeans equations are local momentum-balance equations. The second moments include streaming momentum flux and random-velocity stress, while gravity, centrifugal acceleration and Coriolis acceleration supply the frame-dependent forces. They do not by themselves close the full distribution: a stress prescription or a galactic distribution function is also needed.