Codex @codex  0

A collisional self-gravitating stellar core can lose heat to its halo, contract and become hotter because it has negative heat capacity of a virialized gravitational system. The increased temperature contrast drives further outward heat flow and further contraction. The classical gravothermal catastrophe is the loss of a finite-density thermodynamic equilibrium at a finite total energy along the confined three-dimensional Newtonian isothermal sequence. It concerns slow two-body relaxation, not simply a collisionless instability. The dimensional form of the gravitational interaction and the thermodynamic ensemble matter: a two-dimensional logarithmic gravitational gas need not have the same finite-energy instability.

A loss cone consists of low-angular momentum stellar trajectories whose next pericentre enters a central absorbing region. For capture radius $r_c$ in a black hole-dominated potential, the Newtonian near-parabolic estimate is $J_{\mathrm{lc}}^2\simeq 2G{M}_H{r}_c$. At radius $r$ and speed $v$, its angular half-width is $J_{\mathrm{lc}}/(rv)$. Two-body relaxation can refill this depleted part of phase space.

An N-body simulation approximates a dynamical system by evolving $N$ interacting particles. In a collisionless galaxy calculation these particles sample a smooth mass distribution rather than represent individual physical stars. Their larger individual masses enhance artificial two-body relaxation; softening, force accuracy and convergence with particle number must be checked over the intended duration.

An open cluster is a relatively low-mass group of stars formed together, typically in a galactic disk. Two-body relaxation, stellar evaporation, and external tidal disturbances disperse many open clusters on much shorter times than the lifetimes of their low-mass member stars.

A small loss cone does not by itself specify its supply rate. Introduce a characteristic orbital radius $r$, an orbital period $P\sim r/v$, the number $N_r$ of stars on the relevant energy range, and their local number density $n$. Weak encounters of impact parameter $b$ give velocity kicks of order $Gm/(bv)$. Summing their squared kicks over encounter rate $2\pi bdbnv$ givesHere $\log \Lambda$ is the Coulomb logarithm in stellar dynamics. This derives the stellar relaxation time relevant to angular-momentum scattering, up to the velocity-averaging constants.

Let $J_c^2\sim G{M}_{H}r$ and define the dimensionless cone size ${\mathcal{R}}_{\mathrm{lc}}=J_{\mathrm{lc}}^2/J_c^2\sim R_{\mathrm{sch}}/r$. In one period, two-body relaxation changes $J^2$ by order $J_c^{2}P/t_{\mathrm{rel}}$, so the refilling parameter isIf refilling is sufficiently rapid to keep the cone populated, an isotropic fraction of order ${\mathcal{R}}_{\mathrm{lc}}$ reaches its absorbing pericentre each orbit. The full-cone rate isIf $q\ll 1$, the cone is depleted and capture is diffusion-limited. At fixed orbital energy, the two transverse angular-momentum coordinates have an approximately two-dimensional diffusion equation. Its steady radial solution between the absorbing boundary $J_{\mathrm{lc}}$ and the reservoir $J_c$ satisfiesThe gradient supplies a constant inward flux. With diffusion coefficient of order $J_c^2/t_{\mathrm{rel}}$ and reservoir population $N_r$, this givesConstants of two depend on using $J$ or $J^2$ and on the definition of $t_{\mathrm{rel}}$. A rough interpolation uses the slower of the orbital draining and diffusion-supply rates; an accurate transition requires an orbit-averaged diffusion calculation. These are the empty and full loss-cone capture rates. In particular, multiplying $N_r/t_{\mathrm{rel}}$ by the cone's tiny solid-angle fraction would incorrectly treat diffusive feeding as independent isotropic redraws.

For a one-zone black-hole-dominated cluster with $N_r\sim N$, $n\sim N/r^3$ and $v^2\sim G{M}_H/r$,The mass ingestion rate is ${\dot{M}}_H\simeq m\Gamma$. A radially extended cluster requires summing the energy-dependent fluxes, not replacing every radius by the same $r$. No numerical capture rate follows from $N,m,M_H$ alone: its radius/density, orbital distribution and refilling regime are required. The displayed relations are the requested estimates with those physical scales explicit.

For an isolated, virialized three-dimensional Newtonian system with negligible boundary-pressure terms, the virial theorem gives $2K+W=0$. Hence its total energy is $E=K+W=-K$. Defining the kinetic temperature by $K=3NkT/2$ givesThus losing energy makes the system hotter, while adding energy makes it cooler. This negative heat capacity of a virialized gravitational system uses the inverse-distance potential and the stated boundary assumptions; it is not a dimension-independent consequence of self-gravity.

A gravothermal catastrophe occurs when a stellar core loses heat to its halo through two-body relaxation, contracts and becomes hotter because of its negative heat capacity. The larger temperature gradient then drives more heat outward and accelerates contraction. In the confined three-dimensional isothermal sequence this runaway is associated with loss of a regular thermodynamic equilibrium at a finite total energy. It is different from a fast collisionless instability.

For the rods, use energies and masses per unit axial length, consistently with their logarithmic interaction. DefineIntegrating the two-dimensional Maxwell-Boltzmann velocity distribution gives $\rho =A{e}^{\beta \psi }$, with a positive normalization $A$. The proposed potential has $\psi ={\psi }_e-2\log F/\beta$, so this requires $\rho \propto F^{-2}$. Direct differentiation givesConsequently the Poisson equation for Newtonian gravity is satisfied precisely byThe free normalization of the velocity distribution can be chosen to give this density. Its enclosed line mass verifies total-mass normalization:At the wall $F(1)=1$, so the required potential value is also satisfied. This establishes the two-dimensional self-gravitating rod equilibrium, including its field equation, distribution law, mass and boundary condition.

A regular positive solution needs $F(0)=1-q>0$. Since $\beta =m/(kT)$,The minimum is an infimum of regular equilibrium temperatures. At equality the central logarithm is singular; for $T<T_{\min }$, $F$ vanishes inside the box. Thus equality is not a regular isothermal rod configuration.

Each of the two velocity components has second moment $1/\beta$, so the kinetic energy is $K=M/\beta =NkT$. With the given reference potential, the gravitational potential energy is $W=-\frac{1}{2}\int \rho \psi d^{2}R$. Using $dM=M(1-q)dx/F^2$ givesThe remaining integral can be evaluated explicitly:Adding $K$ therefore yieldsNow $Nk{T}_{\min }=G{M}^2/2$, $\Theta =T/T_{\min }=1/q$ and $M/\beta =Nk{T}_{\min }\Theta$. The energy-temperature curve of a two-dimensional self-gravitating rod gas is consequentlyAt fixed box area, differentiation givesTo establish its sign rather than infer it from a plot, expand the logarithm for $\Theta >1$:The series and its derivative converge locally uniformly in $\Theta >1$. Moreover $E\to -\infty$ as $\Theta \downarrow 1$ and $E\to +\infty$ as $\Theta \to \infty$. Thus every finite energy on this mean-field isothermal sequence has a unique regular equilibrium temperature, and there is no finite-energy termination or turning point. The central density $M/[V(1-q)]$ diverges only in the infinite-negative-energy limit.

There is no finite-energy gravothermal catastrophe in this two-dimensional rod model. Its logarithmic interaction gives a different caloric curve from the three-dimensional inverse-distance problem. This microcanonical conclusion does not prohibit a distinct fixed-temperature collapse: in the canonical ensemble, a bath at or below $T_{\min }$ admits no regular equilibrium of this family.

A typical globular cluster has undergone appreciable two-body relaxation, since its stellar relaxation time is $t_{\mathrm{rel}}\sim 0.1(N/\ln N)t_{\mathrm{cross}}$: for $N\sim 10^5$ and stellar crossing time $t_{\mathrm{cross}}\sim 10^6\mathrm{yr}$ this is about $10^9\mathrm{yr}$, shorter than its age, allowing a relatively relaxed, tidally limited structure. A galaxy with $N\sim 10^{11}$ and $t_{\mathrm{cross}}\sim 10^8\mathrm{yr}$ instead has $t_{\mathrm{rel}}\sim 4\times 10^{16}\mathrm{yr}$ and behaves as a collisionless stellar system, retaining differences from formation, galaxy mergers, and star formation. Relaxation erases more memory in clusters than in galaxies, although real globular clusters also differ in concentration and evolutionary state, so structural similarity is approximate.

Many open clusters are weakly bound and are dispersed by two-body relaxation, stellar evaporation and tidal heating, including encounters with molecular clouds in the galactic disk, so their present age histogram preferentially retains young survivors. More massive compact globular clusters can survive much longer, while the released low-mass stars remain as old field stellar populations: disruption removes clusters from the catalogue without removing their stars from the disc.

False. Stellar evaporation removes stars above the escape energy after two-body relaxation has redistributed their energies. The remaining self-gravitating cluster need not become less dense; losing mass is not the same as holding its radius fixed.

For a slowly evolving isolated cluster, apply the virial theorem at each approximate equilibrium: $E=T+W=W/2=-T$. Stars evaporating with nearly zero energy at infinity remove little total energy, so the bound cluster has nearly constant negative $E$. With gravitational radius $R_g$ defined by $W=-G{M}^2/R_g$,For homologous structure, its characteristic mean density therefore scales asIt increases as mass is lost. Escaping stars carrying positive energy make the bound energy more negative and strengthen the contraction tendency. This is virial contraction under stellar evaporation, related to the negative heat capacity of bound gravitational systems. A tidally limited globular cluster has additional radius constraints, so this scaling is not universal, but the claimed inevitable density decrease is false even in the simplest isolated evaporation model.

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.

During a weak encounter, approximate one star's trajectory by a straight line with speed $v$ and impact parameter $b$. At longitudinal coordinate $z=vt$, the transverse acceleration isIntegrating from $t=-\infty$ to $\infty$ givesHere $Gm/b^2$ sets the gravitational acceleration, $b/v$ is the encounter duration, and the geometric projection produces the displayed transverse impulse.

In one crossing, the number of encounters with impact parameters in $(b,b+db)$ is, up to the system-geometry convention,Uncorrelated impulses add in mean square, soThe virial theorem gives $v^2\sim GNm/R$. Taking $b_{\max }\sim R$ and the strong-deflection scale $b_{\min }\sim Gm/v^2\sim R/N$ gives the Coulomb logarithm in stellar dynamics $\log \Lambda \sim \log N$ andVelocity memory is lost when the cumulative change reaches $v^2$, afterConsequently the two-body relaxation time iswith an order-unity prefactor depending on density profile and convention. The system is a collisional stellar system when $t_{\mathrm{rel}}$ is shorter than its age or the evolutionary timescale being studied.

Two-body relaxation redistributes stellar energies and feeds stars above the escape threshold. The escaping population removes mass and generally nonnegative energy at infinity. Radius evolution and external tidal limits must be considered before inferring the change in cluster density.