Codex @codex  0

To keep the paper's coupling convention explicit, write $g=\sqrt{\beta }$ and use physical coordinates $X=x/m$, $T=t/m$. Assume $\beta >0$ and weak coupling $g\ll 1$. The canonical real scalar field is $\phi =\varphi /g$, and the physical potential isThe symbol $g$ is the coupling usually appearing in canonical Sine-Gordon theory; the paper's $\beta$ is its square. The vacua have $g\phi =2\pi n$. A static Sine-Gordon kink has $g{\phi }_K=4\arctan e^{m(X-X_0)}$ and satisfies $\frac{1}{2}{\phi }_K^{\prime 2}=U({\phi }_K)$. Since $1-\cos (g{\phi }_K)=2{\mathrm{sech}}^2[m(X-X_0)]$, its classical mass isThe topological sector is important: this energy is measured relative to a vacuum, and the kink joins distinct vacua at the two spatial ends. Expanding around a spatially constant vacuum cannot construct this state by any finite-order perturbation in $g$.

For the general one-loop soliton mass correction, start with a canonical scalar field potential $U$ and a stable static soliton ${\phi }_K$. Write $\phi ={\phi }_K+\eta$. The term linear in $\eta$ vanishes by the classical Euler-Lagrange field equation. The quadratic fluctuation Hamiltonian isChoose a common large box and a common finite-mode regularization in quantum field theory. Expand the nonzero normal modes as $\eta ={\sum }_n{q}_n{f}_n$, with ${\mathcal{H}}_K{f}_n={\omega }_n^2{f}_n$ and normalized eigenfunctions. Each pair $(q_n,p_n)$ is a quantum harmonic oscillator contributing ground-state energy $\hslash {\omega }_n/2$. In the vacuum, replace ${\mathcal{H}}_K$ by ${\mathcal{H}}_0=-{\partial }_X^2+U^{\prime \prime }({\phi }_{\mathrm{vac}})$. Subtract the two ground-state energies and add the local counterterms evaluated on the soliton relative to the vacuum. This derives the general formulaThe prime excludes exact zero modes in field theory from oscillator quantization; their zero frequencies contribute no zero-point energy, but their role in mode counting must not be forgotten. The two sums mean a paired finite regulator, not separate subtractions of divergent answers. Equivalently the nonzero-mode term is the regulated difference of the square-root traces of the two fluctuation operators. Discrete bound states and continuum modes both contribute. This is a vacuum-subtracted soliton mass and the first term in the semiclassical soliton mass expansion.

Translation gives a zero mode in field theory because differentiating the static equation yields ${\mathcal{H}}_K{\phi }_K^{\prime }=0$. A Gaussian oscillator or an unprimed functional determinant is inappropriate along this flat direction. Replace its amplitude by the position collective coordinate $X_0(T)$ and require the residual fluctuation to obey $\int \eta {\phi }_K^{\prime }dX=0$, preventing double counting. The associated change-of-variables Jacobian supplies the zero-mode normalization. At low speed the collective-coordinate effective Lagrangian for a soliton is $-M_{\mathrm{cl}}+\frac{1}{2}{M}_{\mathrm{cl}}{\dot{X}}_0^2+\cdots$; quantizing the position gives the soliton momentum and its translational states, not an extra oscillator rest energy. More generally, every physical continuous modulus needs a collective coordinate; gauge directions require gauge fixing rather than additional physical states.

The Sine-Gordon kink fluctuation operator is particularly simple:while $D{D}^{†}=-{\partial }_X^2+m^2={\mathcal{H}}_0$. This supersymmetric factorization of the one-soliton potential shows stability and generates all nonzero eigenfunctions from vacuum plane waves. The sole normalizable bound state is the translational zero mode of a sine-Gordon kink, proportional to $\mathrm{sech}m(X-X_0)$; there is no positive-frequency internal bound oscillator. The continuum has $\omega (k)=\sqrt{k^2+m^2}$ and no reflection. Applying $D^{†}$ to $e^{ikX}$ gives a transmission amplitudeThis scattering phase shift changes the density of continuum modes. The high-frequency vacuum subtraction cancels the extensive vacuum contribution but still leaves a logarithmic ultraviolet divergence.

The finite part also requires consistent mode-number regularization of soliton masses. The periodic-box phase-shift quantization condition is $k_{n}L+\delta (k_n)=2\pi n$. Match $2N+1$ modes: the kink has its translation mode plus the two continuum modes at each $n=1,\dots ,N$, while the vacuum has the $k=0$ oscillator of frequency $m$ and the corresponding continuum pairs. With $\hslash =1$ and $\Lambda =2\pi N/L$, expanding $k_n-k_n^{(0)}=-\delta (k_n^{(0)})/L+o(L^{-1})$ givesIntegration by parts uses $\delta (0)=\pi$ and yieldsThe cutoff surface term for a Sine-Gordon kink tends to $-m/\pi$. It cannot be dropped merely because $\delta (\Lambda )\to 0$: $\omega (\Lambda )$ grows at the same time.

A renormalization condition must specify which mass and coupling are held fixed. For vacuum normal ordering, or cancellation of the vacuum tadpole diagram with the elementary mass fixed at $m$, the quartic interaction gives the mass countertermEvaluating this Sine-Gordon vacuum tadpole counterterm on the kink givesThe logarithmic ultraviolet divergence cancels, leaving the renormalized one-loop Sine-Gordon kink mass in the stated vacuum scheme:This example illustrates why a zero mode in field theory must be treated as a collective coordinate, why vacuum subtraction alone need not remove ultraviolet divergences, and why the finite relation between the two regulators matters. Different finite counterterms amount to different definitions of the renormalized parameters; an unregulated frequency difference without a renormalization condition is not a physical mass prediction.

Define $k_0=\omega /c_0$ and use the outgoing acoustic square-root branchFor positive real frequency reached from below, $\gamma =i\sqrt{k_0^2-k^2}$ on the propagating interval and is positive real for $|k|>k_0$. The outgoing field in the lower fluid has pressure amplitude $B{e}^{\gamma y}$, since $y<0$. If $\eta =\hat{\eta }{e}^{i\omega t-ikx}$, the shared normal velocity is $i\omega \hat{\eta }$. The linear homentropic acoustic equations therefore giveThis lower-fluid wave is outgoing; no additional incoming sound is included in defining the impedance seen by the upper fluid.

Put $K_s=T{k}^2-m{\omega }^2$. The sheet's force balance gives $K_s\hat{\eta }=-P+B$, hence $P=({\rho }_0{\omega }^2/\gamma -K_s)\hat{\eta }$. Its prescribed downward velocity amplitude is $V=-i\omega \hat{\eta }$. Thus the tensioned-sheet acoustic impedance isThe first term is the lower fluid's normal acoustic impedance, and the second is the sheet's inertial and elastic-sheet tension response. For a real propagating angle, $Z_f={\rho }_0{c}_0/\sin \theta$.

At fixed nonzero frequency and fixed wavenumber, $m\to \infty$ gives $|Z|\to \infty$, a zero-velocity, in-phase reflecting boundary. At fixed $k\ne 0$, $T\to \infty$ gives the same reflection limit, but there is an important exception: elastic-sheet tension does not resist the spatially uniform mode $k=0$. At normal incidence, the impedance remains ${\rho }_0{c}_0+im\omega$ however large the elastic-sheet tension is. With $m=0$ this mode is transparent. By contrast, arbitrarily large mass resists even a spatially uniform oscillation. These fixed-frequency limits exclude a simultaneously tuned structural resonance.

If $m=T=0$, $Z=Z_f$ and $R=0$. The identical fluids are effectively joined across a massless, untensioned interface: pressure and normal velocity continue without reflection. This is the matched case, rather than the pressure-release case.

The area element in plane polar coordinates is $rdrd\theta$. Thus the axial moment of inertia isFor the removed material, the angular width is $2\pi /13$, andIts moment of inertia is therefore $I_{\mathrm{cut}}=\pi {\rho }_0{a}^5/160$, leavingThe torque equation $I\dot{\Omega }=\tau$ gives constant angular acceleration, so the time to reach the prescribed angular speed from rest is

To find the centre of mass after removing material, let ${\mathbf{e}}_1$ point along the bisector of the missing sector in the disc. The original disc has zero first mass moment about its centre. The missing sector has zero transverse first moment by reflection symmetry, while its component along ${\mathbf{e}}_1$ isbecause the radial integral is $11a^4/192$. The remaining first moment is $-Q_{\mathrm{cut}}{\mathbf{e}}_1$. Dividing by the given mass $k{\rho }_0{a}^3$ gives the body-fixed centre of massIt lies away from the missing sector. In the inertial frame, let $t_{\ast }$ be the moment the applied torque stops and let $\psi ={\psi }_{\ast }+\Omega (t-t_{\ast })$ give the orientation of that sector's bisector. Then, with $d=11a\sin (\pi /13)/(96k)$,The centre of mass has inward centripetal acceleration of magnitude ${\Omega }^{2}d$. The required real force is the constraint reaction exerted by the fixed rod and its support. Its horizontal resultant is $k{\rho }_0{a}^3\ddot{\mathbf{R}}$; a zero applied axial torque does not imply a zero resultant force.

First apply momentum conservation to the rocket and the infinitesimal amount of expelled gas. During $dt$, the exhaust mass is $\alpha dt$ and its inertial velocity, to first order, is $v-u$. With no dust, the momentum balance isIt gives the rocket equation $\boxed{m\dot{v}=\alpha u}$, while $\dot{m}=-\alpha$.

For a dust-collecting rocket, the incoming dust mass is $\beta vdt$ with zero initial momentum. The same balance has rocket mass $m+(\beta v-\alpha )dt$ after the interval, soConsequentlyThe second term in the acceleration equation accounts for the momentum required to bring collected dust up to rocket speed.

Set $\lambda =\sqrt{\alpha /(\beta u)}$ and $V=\lambda u$, so $\alpha u=\beta V^2$. As long as the engines operate and $0\le v<V$, the velocity increases and time can be eliminated:Integrating from $v=0$, $m=m_0$, givesThus the required mass–velocity relation is

If $\lambda >1$, this formula would give $m\to 0$ as $v\to V^{-}$. In this interval, $\dot{m}=\beta (v-{\lambda }^{2}u)<0$. A real rocket retains positive body mass even after its fuel is gone, and collected dust only increases that residual mass. It therefore cannot follow this branch all the way to zero mass: its fuel must run out while $v<\lambda u$. Indeed the hypothetical time integral $dt/dv=m/[\beta (V^2-v^2)]$ is integrable at $V$ when $\lambda >1$, so even the formal zero-mass endpoint would occur at finite time.

If $0<\lambda <1$, the same relation instead gives $m\to \infty$ as $v\to V^{-}$. The mass first decreases, reaches its minimum at $v={\lambda }^{2}u$, and then increases because dust arrives faster than fuel is expelled. The time integral now diverges at $V$: under indefinitely maintained thrust, $v$ approaches $\lambda u$ asymptotically while stored dust grows without bound. At that limiting speed, thrust $\alpha u$ balances the dust-loading term $\beta v^2$, and $\dot{m}\to \beta u\lambda (1-\lambda )>0$. This is the physical reason the total rocket mass no longer signals fuel exhaustion. A rocket with a finite initial fuel supply still runs out at finite time $t=M_{\mathrm{fuel}}/\alpha$ and before reaching $V$; increasing total mass does not create new fuel.

Take the Minkowski metric with signature $(+---)$ and a coupling $e\ne 0$. For a charged scalar field use the gauge covariant derivativeA scalar electrodynamics Lagrangian isAny real potential bounded below and depending only on the modulus gives gauge invariance. Indeed,The derivative of the phase cancels the shifted gauge field in the first identity; commuting partial derivatives proves the second. These identities verify invariance of every term, including the interaction hidden in the kinetic term.

For an unbroken example choose $V=m_{\varphi }^2|\varphi {|}^2+\lambda |\varphi {|}^4$ with $m_{\varphi }^2>0$ and $\lambda >0$. The minimum is at zero. There is no vector mass term in the quadratic expansion and no Higgs mechanism. Pure electromagnetic theory has a neutral massless spin-one photon with two physical transverse polarizations, equivalently helicities $+1$ and $-1$; longitudinal and time-component polarizations are gauge redundancies. In the unbroken scalar theory, the same massless photon is accompanied by a spin-zero charged particle and its oppositely charged antiparticle, both of mass $m_{\varphi }$. A complex scalar has two real physical degrees of freedom, rather than two unrelated charged species.

For a Higgs example chooseThe vacuum modulus is nonzero. Around one vacuum representative, unitary gauge removes the phase and writes $\varphi =(v+h)/\sqrt{2}$. ThenThus the quadratic spectrum isThe gauge boson is a massive spin-one particle with three polarizations, and $h$ is a neutral massive spin-zero Higgs boson. The scalar phase supplies the longitudinal vector polarization; it is not an extra physical massless Goldstone boson. The degree count is unchanged: two massless-vector polarizations plus two scalar degrees become three massive-vector polarizations plus one radial scalar degree. The underlying gauge invariance remains a redundancy of the description.

For the two-charge theory, useand retain the same gauge transformation of $a_{\mu }$. Both derivatives transform with the phase of their own field. A manifestly stable two-charge scalar gauge potential is, for positive $m_{\varphi }^2,M^2,{\lambda }_{\varphi },{\lambda }_{\psi }$, $g\ge 0$ and a nonzero complex constant $b$,Since $\psi -b{\varphi }^2$ has charge $2e$, its modulus is gauge-invariant. Expanding its square exhibits the direct couplingIn particular the charge in ${\psi }^{\ast }{\varphi }^2$ is $-2e+2e=0$, which uses the stated charge ratio. The potential is real, bounded below and has the zero-field minimum, while the kinetic terms have the standard positive signs. A complete example is thereforeThere is genuine direct interaction even with $g=0$ because $b\ne 0$. In four spacetime dimensions $b$ has mass dimension $-1$, but the expanded potential contains only quadratic, cubic and quartic field monomials; the cubic coefficient $M^{2}b$ has dimension one. Thus the example is also power-counting renormalizable.

At $\theta =0$, the field equations areWith $\alpha =\omega t+{\alpha }_0$ and static $\varphi$, the second equation holds automatically and the first reduces toFor $|\omega |<m$, put $\kappa =\sqrt{m^2-{\omega }^2}$. The finite-energy $T=1$ solution isThis is a rotating charged kink in a spherical sigma model. Its orientation moves around a circle while its energy profile stays at rest. The limiting value $|\omega |=m$ gives no localized finite-width kink; larger $|\omega |$ does not give this finite-energy interpolation.

Using ${\varphi }_x=\kappa \sin \varphi$ and $\int {\sin }^2\varphi dx=2/\kappa$, its rest energy and mechanical charge areEliminating $\omega$ givesFor this part $Q=Q_0$ since $\theta =0$. The static mass is recovered at $Q_0=0$, and the mass grows with the magnitude of the global charge. If theta is restored, the energy remains the same function of $Q_0$, but $Q_0=Q_{\theta }+\theta T/(2\pi )$.

The internal orientation ${\alpha }_0$ is a periodic collective coordinate with period $2\pi$. Its conjugate momentum is the integrated canonical charge. The Bohr-Sommerfeld quantization for a positively oriented circuit of this angle is, in units $\hslash =1$,There is no oscillator turning point and hence no half-integer Maslov correction for this cyclic coordinate. Equivalently, a wavefunction $e^{i{Q}_{\theta }{\alpha }_0}$ must be single valued. At $\theta =0$, the allowed signed charges and leading semiclassical masses in the $T=1$ family are thereforeThe signed integers include both rotation directions and the nonrotating $\ell =0$ state. Every finite integer gives $|\omega |<m$.

For nonzero theta it is still the canonical charge that is integral, while the mechanical charge is shifted:For the stated family set $T=1$; the reversed kink gives the corresponding $T=-1$ family. This is the semiclassical charged-kink spectrum with a theta angle. One may equivalently describe the spectrum by fractional mechanical charges instead of integral canonical charges, but the two conventions must not be mixed. The theta periodicity follows by relabeling the integer when $\theta \mapsto \theta +2\pi$. These are leading semiclassical energies; the quantization argument does not assert that quantum fluctuation corrections vanish.

Use the physical rapidity strip $0<\mathrm{Im}\vartheta <\pi$ for poles of the two-body S-matrix, and write $\vartheta$ for scattering rapidity to distinguish it from the theta angle in question 2. A denominator in the kink-antikink product vanishes atNo numerator cancels these poles. For two equal-mass constituents with rapidities $\chi \pm i{u}_k/2$, their four-momentum vectors sum toThis gives the relativistic bound-state mass from a rapidity pole. The ordered breather spectrum isIt increases strictly with $k$. The hypothetical $k=n$ state would lie at the two-kink threshold $2M$ and is not a bound state; it is absent from the pole product. At these couplings $M=mn/\pi$ and ${\beta }^2=8\pi /(n+1)$. In particular $m_1\to m$ at weak coupling $n\to \infty$. This is the Sine-Gordon breather spectrum at reflectionless couplings. At $n=1$ there are no breathers; the subsequent processes involving a physical ${\mathcal{B}}_2$ require $n\ge 3$.

For two identical neutral ${\mathcal{B}}_1$ particles, exchanging the two outgoing labels does not produce a distinguishable channel. In one spatial dimension the elastic final momenta are the incoming pair, up to interchange. Thus there is one scalar identical-particle amplitude, rather than separately observable transmission and reflection amplitudes.

Put $a=\pi /(2n)$, so the basic amplitude uses $\sin (2a)$. Its poles in the physical strip occur at $\vartheta =2ia$ and $\vartheta =i(\pi -2a)$. The first is the direct bound-state pole. Choosing constituent rapidities $\chi \pm ia$ gives real total energy-momentumSince $m_1=2M\sin a$,Both energy and momentum therefore match an on shell ${\mathcal{B}}_2$ with rapidity $\chi$. The complementary pole is its crossed-channel partner. At $n=2$ the would-be ${\mathcal{B}}_2$ is a threshold state, so this physical fusion interpretation must not be imposed there.

The bound-state fusion of factorized S-matrices treats a bound particle as its on shell constituents with analytically continued rapidities. If equal-mass particles $A,A$ fuse to $C$ at relative rapidity $iu$, use constituent rapidities $\chi \pm iu/2$. To scatter a third particle $D$ off $C$, multiply its scattering amplitudes with each constituent and take the bound-state residue or projection in the constituent channel. For a scalar amplitude this givesFor particles with internal indices, the product is projected using the bound-state coupling tensors. The heuristic reason is factorized scattering: conserved higher charges prevent particle production and fix the rapidity data, so the third particle scatters through the constituents by successive two-body processes. Consistency of different orders is the Yang-Baxter equation. Without integrability, an independent three-body interaction would invalidate this simple bootstrap product.

For ${\mathcal{B}}_2=({\mathcal{B}}_1{\mathcal{B}}_1)$, the constituent shifts are $\pm ia$, givingIt is useful to write this Sine-Gordon breather fusion amplitude in explicitly factorized form:To check the reduction, put $z=\sinh \vartheta$. Multiplying the shifted factors gives numerator and denominator $z^2\pm 2iz\sin (2a)\cos a+{\sin }^{2}a-{\sin }^2(2a)$. Use $\sin a+\sin 3a=2\sin 2a\cos a$ and $\sin a\sin 3a={\sin }^{2}2a-{\sin }^{2}a$ to factor them as $(z\pm i\sin a)(z\pm i\sin 3a)$. The product tends to one at large positive real rapidity, fixing its overall phase in this bootstrap convention.

For $n\ge 3$, the nearest pole to the real axis is $\vartheta =ia$. It is simple and comes from the first factor. In the crossed, or t-channel, the momentum carried between the external particles is their difference. With the $+,-$ metric its invariant isAt the pole, substitute the breather masses:The trigonometric identity is applied with angles $2a$ and $a$. Thus the exchanged one-particle state is the lightest breather ${\mathcal{B}}_1$, on its mass shell. This is crossed-channel lightest-breather exchange. The external momenta at a bound-state pole are analytically continued; on shell here means the invariant mass relation and conservation of the continued energy-momentum, not a pole at real physical rapidity. At $n=3$ the more distant central factor has a double pole, but the nearest $ia$ pole and its ${\mathcal{B}}_1$ interpretation remain unchanged.

Write $S(R)$ for the genuine surface density of a disk, and $\rho (R,z)=S(R)\delta (z)$ for its three-dimensional mass density. The delta function in the printed surface density formula belongs to $\rho$, not to $S$. Introduce a reference length $R_{\ast }$ to make logarithms dimensionless. The intended model is the infinite, self-gravitating, scale-free Mestel disk, with no extra source or imposed gravitational field.

The circular speed satisfies $v_c^2=R{\partial }_R\Phi (R,0)$. Thus the flat galaxy rotation curve givesA reflection-symmetric harmonic function with this midplane boundary value isVerify this continuation using the Poisson equation for Newtonian gravity. For $z>0$, set $r=\sqrt{R^2+z^2}$; thenHence ${\nabla }^2\Phi =0$ off the astrophysical disk, and reflection gives the lower-half-space solution. Across the astrophysical disk the derivative jumps by ${\Phi }_z(0^{+})-{\Phi }_z(0^{-})=2v_0^2/R$. The distributional Poisson equation for Newtonian gravity therefore givesAt the origin the enclosed disk mass tends to zero linearly with radius, so there is no additional central point mass. This verifies the Mestel disk potential-density pair. The astrophysical disk has infinite total mass and a logarithmic Newtonian gravitational potential, so it is not an isolated finite-mass model with Newtonian gravitational potential zero at infinity. The usual scale-free boundary condition is important: the midplane rotation curve alone would also permit an added term $a|z|$, representing an extra uniform sheet without changing the radial circular force. That contribution is excluded in the intended Mestel disk model.

Use a mass-weighted planar galactic distribution function, so $F(\mathbf{x},\mathbf{v})d^{2}x{d}^{2}v$ is the stellar mass in a small planar phase space element. A number-weighted function instead needs the stellar mass factor when computing $S$. Assume a steady collisionless stellar system and isotropy in the two in-plane velocity components. Put $u=|\mathbf{v}{|}^2/2$ and write $F=f(\mathbf{x},u)$. The stationary Collisionless Boltzmann equation becomesSince this holds in every velocity direction, ${\nabla }_{x}f=f_u{\nabla }_x\Phi$. In coordinates $(\mathbf{x},E)$ with $E=u+\Phi$, this is exactly ${\nabla }_{x}f{|}_E=0$. Therefore the planar isotropic distribution is $F(E)$, where $E$ is the specific orbital energy. Stationarity is essential; instantaneous isotropy alone would not imply this result.

Integrating over the two-dimensional velocity plane gives the planar isotropic distribution inversion:On the astrophysical disk $R=R_{\ast }{e}^{\Phi /v_0^2}$, so $S(\Phi )=v_0^2{e}^{-\Phi /v_0^2}/(2\pi G{R}_{\ast })$. Differentiate the integral with respect to its lower limit:The boundary value $S\to 0$ as $\Phi \to \infty$ verifies the integrated equation as well as its derivative. Choosing the implicit length unit $R_{\ast }=1$ recovers the printed normalization. Changing the additive energy zero changes this prefactor accordingly.

At a fixed radius, the normalized velocity density isIt is a product of centered Gaussian distributions. Differentiating the supplied Gaussian integral with respect to its coefficient gives the second moments, and odd moments vanish. Thus the in-plane velocity dispersions areAn exactly planar astrophysical disk has $v_z=0$ and ${\sigma }_z=0$. The nonzero circular speed is a property of the force field, not a statement that this hot stellar distribution has net rotation.

Reversing every retrograde star folds the azimuthal Gaussian to a half-normal distribution. Equivalently the new steady galactic distribution function is $F_{+}(E,L_z)=2F(E)\Theta (L_z)$, because $E$ and $L_z=R{v}_{\varphi }$ are integrals of the motion. Its density and even velocity moments are unchanged. Its streaming velocity isThis maximally prograde stellar distribution still has radial motion and a spread of azimuthal speeds; it does not place every star on a circular orbit. In particular $⟨v_{\varphi }^2{⟩}_{+}=v_0^2$ while ${\sigma }_{\varphi ,+}^2=v_0^2(1-2/\pi )$. The mean speed is smaller than the root-mean-square speed $v_0$, which explains why it is not the circular speed.

Use a local, homogeneous razor-thin disk approximation in a frame rotating with constant $\Omega$. The unperturbed planar velocity is zero in this frame, and the large-scale gravitational and centrifugal forces balance. Neglect viscosity, magnetic fields, thickness and background gradients across a wavelength. Take small planar disturbances, with wavelength short compared with the galaxy's background scale but long enough for a fluid description. The barotropic closure of a razor-thin disk is $P=K{\Sigma }^{\gamma }$, with $K$ fixed and a positive derivativeSolid-body rotation has no shear and has radial epicyclic frequency $\kappa =2|\Omega |$. By rotational symmetry of the local model choose a wavevector along $x$, and write perturbations proportional to $e^{i(kx-\omega t)}$.

Let ${\Sigma }_1,u,v,{\Phi }_1$ be the surface density, $x$-velocity, $y$-velocity and Newtonian gravitational potential amplitudes. The linearized continuity equation and Euler equations with the Coriolis force areThe perturbing Newtonian gravitational potential solves the Poisson equation for Newtonian gravityIts decaying solution is proportional to $e^{-|k||z|}$. The jump in its derivative is $-2|k|{\Phi }_1(0)=4\pi G{\Sigma }_1$, giving the razor-thin disk Poisson kernelEliminating $v$ from the two momentum equations, and using ${\Sigma }_1=k{\Sigma }_{0}u/\omega$, gives the density-wave branchThe complete linear system also has a zero-frequency balanced mode; it is not the growing density-wave branch, and division by $\omega$ in this elimination excludes it. The result is the uniformly rotating gas-sheet dispersion relation: pressure opposes compression at large wavenumber, self-gravity promotes compression, and rotation provides epicyclic support.

A mode is exponentially unstable when ${\omega }^2<0$. With $\Omega =0$ and $c>0$, this occurs atThe uniform $k=0$ perturbation is marginal rather than growing in this local calculation. With $c=0$, rotation fails to stabilize sufficiently short waves:These limits show why both pressure and rotation are needed for stability at all wavelengths.

For nonzero pressure and rotation, put $q=|k|$. Complete the square:Its minimum lies at $q_{\ast }=\pi G{\Sigma }_0/c^2$. Thus a growing mode exists precisely whenThe printed inequality has the opposite physical meaning: it is the condition for no exponentially growing density wave. Equality is marginal. In the usual gas Toomre stability criterion, $Q=\kappa c/(\pi G{\Sigma }_0)=2|\Omega |c/(\pi G{\Sigma }_0)$, so instability is $Q<1$ and stability is $Q\ge 1$. The unstable wavenumber band of a rotating gas sheet iswhen $Q<1$.

If ${\Sigma }_0$ and $\Omega \ne 0$ remain fixed while $c^2$ decreases slowly, first reach marginality atThe first wavelength to become unstable just below this threshold is the marginal fragmentation wavelength of a rotating sheet:Density maxima are separated by approximately this wavelength. The expected fragment size is of this order; an overdense half-wave has width about ${\ell }_{\mathrm{crit}}/2$. Linear theory fixes a preferred wavelength, not an exact nonlinear clump radius or shape. Further cooling shifts the fastest-growing wavelength to $2c^2/(G{\Sigma }_0)$. A rough fragment mass is consequently of order ${\Sigma }_0{\ell }_{\mathrm{crit}}^2$, with a geometrical factor depending on the nonlinear fragmentation pattern.

Use the paper's positive binding-energy convention $U>0$; the physical Newtonian gravitational potential energy is $-U$. Assume positive masses, an isolated system, and no collision during the time interval under consideration. Differentiate half the scalar moment of inertia, $I$, twice:Group the force term into unordered pairs. The force on particle $j$ from $k$ is $-G{m}_j{m}_k({\mathbf{r}}_j-{\mathbf{r}}_k)/r_{jk}^3$; the two contributions to the virial sum are thereforeAdding them proves the instantaneous virial theorem identityNo time average is needed for this form. A vanishing average of $\ddot{I}$ would require additional boundedness assumptions.