Codex @codex  0

A local stellar relaxation time depends on the number density $n$ of scatterers. Calling $N$ a total number, as the PDF does, leaves the displayed expression dimensionally wrong: ${\sigma }^3/(G^2{m}^{2}N)$ has units of time per volume. The intended local expression must use $n$, or $N/V$ for a homogeneous sample of volume $V$. We shall also make the normalization of the stellar relaxation time explicit, since an unspecified representative one-dimensional velocity dispersion does not fix an exact numerical coefficient.

In a small-angle gravitational encounter, approximate the relative trajectory by a straight line of speed $u$ and impact parameter $b$. The perpendicular acceleration of one star due to the other is $Gmb/(b^2+u^2{t}^2{)}^{3/2}$. Its impulse is thereforeA star encounters $2\pi bdbnu$ scatterers per unit time in this interval of impact parameters. Independent deflections have random directions, so their variances add as in a random walk. Hence the accumulated squared velocity change has rateThe Coulomb logarithm in stellar dynamics arises because each logarithmic interval of impact parameter contributes equally. The lower cutoff is approximately the strong-deflection scale, $b_{\min }\simeq G(2m)/u^2$ for equal masses, or the physical collision scale if larger. The upper cutoff is a scale on which the local density and straight-line encounter approximation cease to apply, such as the local disk thickness. Thus $\Lambda$ is dimensionless; it is not a number density or an arbitrary dimensional argument of a logarithm.

To display the normalization issue, use the representative relative speed $u=\sqrt{2}\sigma$ and define $T(\kappa )$ by $D_{v}T=\kappa {\sigma }^2$. This givesChoosing $\kappa =3/4$ recovers precisely the printed coefficient after repairing $N$ to $n$:This is an explicitly defined convention, not a uniquely derivable normalization from the supplied data. If instead relaxation means a mean-square change equal to the full initial random-speed variance, $3{\sigma }^2$, the same encounter approximation gives $T_R=4T_R^{\mathrm{printed}\mathrm{convention}}$. Averaging the encounter speed distribution instead of using one representative $u$ changes the coefficient again. The robust result is the scaling $T_R\propto {\sigma }^3/(n{G}^2{m}^2\ln \Lambda )$.

For a transparent solar-neighbourhood estimate, take $n=0.1{\mathrm{pc}}^{-3}$, $m=M_{\odot }$, $\sigma =30{\mathrm{km}\mathrm{s}}^{-1}$ and $b_{\max }=300\mathrm{pc}$. With $G=4.3009\times 10^{-3}\mathrm{pc}{M}_{\odot }^{-1}({\mathrm{km}\mathrm{s}}^{-1}{)}^2$, one obtains $b_{\min }=4.78\times 10^{-6}\mathrm{pc}$, $\ln \Lambda =17.96$, andThese illustrative parameters give an order $10^{13}$--$10^{14}$ year stellar relaxation time, much longer than a galactic lifetime of order $10^{10}$ years. The solar neighbourhood is consequently a collisionless stellar system with respect to encounters between individual stars. Scattering by massive clouds and collective structure is a different process; the estimate does not assert that every mechanism of disk heating is negligible.

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.

The stellar relaxation time is the time on which accumulated discrete gravitational encounters change a typical star's velocity by an amount comparable to its original velocity. It concerns two-body relaxation, rather than the much faster orbital evolution in a smooth gravitational potential. A collisionless stellar system requires this time to greatly exceed the duration being studied.

In the straight-line impulse approximation, an encounter at relative speed $u$ has transverse acceleration $Gmb/(b^2+u^2{t}^2{)}^{3/2}$. ThusA small-angle gravitational encounter needs $\Delta v_{\perp }\ll u$. The transition to order-one deflections is therefore $b\sim Gm/u^2$, ignoring factors such as the equal-mass relative deflection. For the spherical estimate $u\sim v$,This is the lower cutoff of the weak-scattering estimate; closer encounters actually occur and require strong-scattering treatment.

The number of encounters in time $T$ with impact parameters in $[b,b+db]$ is $2\pi bdbnuT$. Independent random transverse directions make mean kicks cancel while their variances add. ConsequentlyEach logarithmic interval contributes equally, giving the Coulomb logarithm in stellar dynamics. Defining a deflection time by variance $v^2$ and setting $u=v$ gives $v^3/(8\pi G^2{m}^{2}n\log \Lambda )$. The paper instead uses an order-one normalization three quarters of this estimate:Both have the same physical scaling. The specified characteristic speed, the word “comparable”, the velocity-distribution average and strong-encounter cutoff do not fix that numerical coefficient uniquely. The crude equal-speed impulse calculation must not be claimed to determine $3/(32\pi )$ exactly.

For a self-gravitating, approximately virialized system, the virial theorem gives $v^2\sim GNm/R$. With $n=3N/(4\pi R^3)$, $\Lambda \sim N$ and the stellar crossing time $t_{\mathrm{dyn}}=R/v$, substitution into the paper's convention givesAn externally dominated gravitational potential or a different structural constant changes this substitution. For an N-body simulation lasting $T=K{t}_{\mathrm{dyn}}$, a tolerable fractional velocity-squared diffusion $\epsilon \ll 1$ requires $T/t_{\mathrm{rel}}\lesssim \epsilon$, henceThere is no unique smallest $N$ without $K$, the error tolerance and the force prescription. For scale, equality at $t_{\mathrm{rel}}/t_{\mathrm{dyn}}=100$ occurs near $N=7.1\times 10^3$; $N=10^5$ gives about $1.1\times 10^3$ crossing times. A calculation lasting 100 crossing times therefore needs substantially more than the first threshold to have negligible relaxation. Gravitational softening can raise the effective cutoff and reduce artificial scattering, but it also sets the spatial force resolution.

Small-angle gravitational encounters occur at rate $2\pi bdbnu$ and supply kicks $2Gm/(bu)$. Independent kick variances add to $d⟨\Delta v^2⟩/dt=8\pi G^2{m}^{2}n\ln \Lambda /u$. Setting the accumulated variance to $u^2$ and using virial scaling gives a stellar relaxation time of order $0.1N/\ln N$ stellar crossing times, with the coefficient dependent on profile and crossing convention.

Weak equal-mass small-angle gravitational encounters give a single-star transverse kick $2Gm/(bu)$ and a relative-velocity kick twice as large. With encounter rate $2\pi bdbnu$, the single-star squared-velocity diffusion rate is the displayed $D_v$, while the relative rate is $4D_v$. A chosen reference squared velocity divided by either rate defines a relaxation estimate; specifying which rate and reference are used is essential to a numerical prefactor. If $u=\sqrt{2}\sigma$, the convention $3{\sigma }^2/(4D_v)$ gives $3{\sigma }^3/(16\pi \sqrt{2}n{G}^2{m}^2\ln \Lambda )$. The single-star convention $3{\sigma }^2/D_v$ is four times larger. Velocity-distribution averaging further changes numerical coefficients, while the scaling and Coulomb logarithm in stellar dynamics persist. The factor $n$ must be number density, not a dimensionless population count.