Codex @codex  0

For the real scalar field, the mostly-plus Minkowski metric gives a positive kinetic energy. Its canonical momentum and Hamiltonian operator areThe Euler-Lagrange field equation is $(\square -m^2)\varphi =0$, or $({\partial }_t^2-{\nabla }^2+m^2)\varphi =0$. Canonical quantization of a real scalar field promotes the fields to Hermitian operators and replaces the equal-time Poisson brackets by canonical commutation relations:The Heisenberg equation of motion then gives $\dot{\varphi }=\pi$ and $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$, so the classical field equation remains an operator identity. Products at coincident points need a regulator; equivalently, the field is an operator-valued distribution.

Resolve the free field into positive- and negative-frequency plane waves. With $E_{\mathbf{k}}=\sqrt{|\mathbf{k}{|}^2+m^2}$ and $k\cdot x=-E_{\mathbf{k}}t+\mathbf{k}\cdot \mathbf{x}$,Reality pairs the two terms by Hermitian conjugation. The equal-time canonical commutation relations are equivalent toFor example, the two mixed terms in $[\varphi ,\pi ]$ each supply half the Dirac delta function. Thus the normalization $1/\sqrt{2E_{\mathbf{k}}}$ is fixed, not optional. Each independent mode is a quantum harmonic oscillator, with $a$ and $a^{†}$ its annihilation operator and creation operator.

Choose the Fock vacuum by $a(\mathbf{k})|0⟩=0$. Repeated creation operators construct the bosonic Fock space. Substitution into the Hamiltonian operator yieldswhere a finite volume $V$ can be used as a regulator. Normal ordering sets the flat-space vacuum reference to zero. The spatial momentum operator isIn particular, $[H,a^{†}(\mathbf{k})]=E_{\mathbf{k}}{a}^{†}(\mathbf{k})$ and $[\mathbf{P},a^{†}(\mathbf{k})]=\mathbf{k}{a}^{†}(\mathbf{k})$. A one-particle state therefore has positive energy, momentum $\mathbf{k}$, and invariant mass $m$. A scalar field transforms in the trivial spin representation, so its particles have spin zero. Because the field is real there is only one set of oscillators: the particle is its own antiparticle, with no distinct conserved particle-minus-antiparticle charge. The free theory does conserve its occupation-number sum, though generic scalar interactions need not.

A momentum eigenstate extends throughout space. A localized particle is described by a superposition such asIts wave packet evolves through the phase $e^{-i{E}_{\mathbf{k}}t}$; a narrow packet has group velocity $\mathbf{v}=\mathbf{k}/E_{\mathbf{k}}$, whose magnitude is at most one. Localization and dispersion concern these superpositions, rather than classical trajectories attached to individual field modes.

Relativistic spacetime behavior is encoded by microcausality. Directly from the oscillator commutators,This is a Lorentz-invariant distribution and vanishes at equal time. Any spacelike separation can be transformed to an equal-time separation, so $[\varphi (x),\varphi (y)]=0$ for spacelike $x-y$. Thus local operations at spacelike-separated points commute. The vacuum Wightman function need not vanish outside the light cone, but that correlation does not transmit a controllable signal. The Feynman propagator iswhich propagates positive energy forward and negative energy backward in the time-ordered product. For causal response one uses the retarded Green function, whose support lies in or on the future light cone.

Finally, commuting creation operators implyThe multiparticle state is symmetric under the particle exchange operator, which is Bose statistics. For a normalized single mode, the occupation states are $(a^{†}{)}^n|0⟩/\sqrt{n!}$ with $n=0,1,2,\dots$: any number of identical bosons can occupy it. There is no Pauli exclusion principle for these particles. This construction agrees with the Spin-statistics theorem relating integer spin to bosonic exchange symmetry under relativistic locality and positive-energy assumptions. It demonstrates the required scalar case without assuming that theorem as the quantization prescription.

Use the Minkowski metric $\mathrm{diag}(1,-1)$ and the angular field $u=\beta \varphi$. The Euler-Lagrange equation becomes $u_{tt}-u_{xx}+m^2\sin u=0$. For the dimensionless light-cone coordinatesthis is $u_{\xi \eta }=\sin u$. This normalization keeps the coupling in the relation between the physical field and its angle, rather than silently setting $\beta =1$.

A Bäcklund transformation is a system of first-order differential relations that maps a solution to another solution. One convention for the Sine-Gordon Bäcklund transformation uses a nonzero parameter $a$ and defines $v={\mathcal{B}}_a[u]$ byA compatible initial value or integration constant selects a particular transformed solution. Put $w=(v+u)/2$ and $d=(v-u)/2$. Differentiating gives $w_{\xi \eta }=\cos d\sin w$ and $d_{\xi \eta }=\cos w\sin d$. Their sum and difference yieldThus compatibility of the first-order relations contains the field equations for both fields, and the transformed physical field $v/\beta$ also solves Sine-Gordon theory.

The transformation provides a generating conservation law. On the branch close to the original field write $v=u+2d$. Its first equation isIt recursively determines a formal small-$a$ expansion of $d$ in local derivatives of $u$:At every subsequent order the new coefficient occurs linearly, so this recursion continues indefinitely. The second Bäcklund equation and the first imply the exact identityIndeed ${\partial }_{\eta }\cos d=-a\sin d\sin (u+d)$ and ${\partial }_{\xi }\cos (u+d)=-a^{-1}\sin d\sin (u+d)$. Comparing powers of $a$ therefore gives local conserved currents. If ${\partial }_{\eta }{P}_j={\partial }_{\xi }{Q}_j$, our coordinate convention givesFor localized fields approaching vacua at spatial infinity, the boundary flux vanishes and ${\mathcal{I}}_j$ is conserved. The formal expansion need not converge: each coefficient is a separately exact local conservation law.

For example $P_0=u_{\xi }^2/2$, $Q_0=1-\cos u$, a light-cone combination of energy and momentum. The next coefficient is a derivative improvement, so it contributes no independent charge under the same decay conditions. At order $a^2$, remove the improvement generated by ${\partial }_{\xi }(u_{\xi }{u}_{\xi \xi })$ and multiply by $-2$. One obtains the genuinely higher conservation lawIt can also be checked directly using $u_{\xi \eta }=\sin u$. Continuing the recursion, and using the opposite light-cone construction, produces the local conserved-charge hierarchy of sine-Gordon theory, with infinitely many nontrivial higher-spin charges after derivative improvements are removed. This is the Bäcklund generating current for sine-Gordon conserved charges. The hierarchy is the characteristic field-theory form of classical integrability; an ordinary energy conservation law alone would not supply these constraints.

The same transformation constructs solutions rather than only currents. Starting from the vacuum $u_0=0$, its two first-order equations integrate toFor $a=e^{\rho }>0$ the exponent is $m\cosh \rho x-m\sinh \rho t+{\delta }_a$. This is a Sine-Gordon kink with velocity $\tanh \rho$, center set by ${\delta }_a$ and classical rest mass $8m/{\beta }^2$. Negative parameters can supply the corresponding opposite-orientation seeds.

The allowed Bianchi permutability for sine-Gordon Bäcklund transformations gives the two-step field algebraically. For vacuum seed and $a\ne b$,Branches of the inverse tangent must be continued smoothly; its principal value alone does not specify the vacuum labels of a multi-kink field. All angles here are $u=\beta \varphi$. In the physical-field version each field difference in the superposition formula carries $\beta$; the formula printed without it implicitly uses the $\beta =1$ angular-field convention.

To display two real scattering solutions, take $0<v<1$, set $\Gamma =(1-v^2{)}^{-1/2}$, and put $X=m\Gamma x$, $Y=m\Gamma vt$. Choosing $a=e^{\rho }$, $b=e^{-\rho }$ with $v=\tanh \rho$ and zero phase constants yields the Sine-Gordon kink-antikink scattering solutionChoosing instead $b=-e^{-\rho }$ yields the Sine-Gordon two-kink solutionThe latter has net angular winding $4\pi$, while the former has zero net winding. At large positive or negative time they separate into localized kinks with velocities $\pm v$. Repeated commuting transformations give general multi-soliton fields. Continuing the kink-antikink velocity to an imaginary value also gives a real Sine-Gordon breather:up to the irrelevant overall field sign and translations.

The scattering interpretation ties the construction to the conserved hierarchy. In an exact multi-soliton sector, the incoming species and rapidities reappear after the collision, up to permutation: the solitons change positions, not their asymptotic shapes or velocities, and the collision emits no radiation. In the two-kink example the large-time centers obey $|X|=|Y|-\log v+o(1)$, so a right-moving trajectory acquires a shift $-2\log v/(m\Gamma )$. These shifts encode the interaction even though the collision is elastic. The commuting construction makes the net displacement in a many-soliton collision the sum of its pairwise displacements, independent of how the collisions are ordered; this is pairwise additivity of soliton shifts and the classical counterpart of factorized scattering.

Conservation of the entire hierarchy is much stronger than conservation of energy and momentum: its independent rapidity-weighted sums constrain the whole asymptotic soliton data and rule out particle production in this sector. Generic initial fields may also contain radiation scattering data, so this statement is about the exact soliton collisions, not a claim that every initial field is a pure soliton. The Bäcklund map unifies the conserved hierarchy, explicit soliton construction and elastic, pairwise scattering picture of classical integrability.

Fix ${ϵ}^{01}=+1$ with $x^0=t$, $x^1=x$. The Noether current for the unit shift of $\alpha$ is the derivative of the Lagrangian with respect to ${\partial }_{\mu }\alpha$:The angular field equation is ${\partial }_{\mu }{j}^{\mu }=0$. At $\theta =0$ this gives the conserved Noether chargeFor nonzero $\theta$, ${ϵ}^{10}=-1$, so $j^0$ contains $+\theta {\partial }_x(\cos \varphi )/(4\pi )$. Its integral is a boundary term:The theta term is locally a total derivative, so it does not change the bulk equations or the energy, but it changes the canonical charge by this topological shift. This is the theta-angle shift of a sigma-model kink charge. Reversing the convention for ${ϵ}^{01}$ reverses the displayed theta shifts. Distinguishing the mechanical charge $Q_0$ from the canonical symmetry charge $Q_{\theta }$ will be important for quantization.

For $T=1$, choose the positive Bogomolny equation and constant angular orientation $\alpha ={\alpha }_0$. Integrating $d\varphi /\sin \varphi =mdx$ givesThis interpolates from $\varphi (-\infty )=0$ to $\varphi (+\infty )=\pi$, so $T=-(-1-1)/2=1$. Its useful profiles are $\sin \varphi =\mathrm{sech}[m(x-x_0)]$ and $\cos \varphi =-\tanh [m(x-x_0)]$. It obeys the second-order field equation because differentiating ${\varphi }_x=m\sin \varphi$ gives ${\varphi }_{xx}=m^2\sin \varphi \cos \varphi$; the angular field equation holds since $\alpha$ is constant. ThusBoth $x_0$ and ${\alpha }_0$ are collective coordinates: translation and global angular rotation do not change the energy. This static solution is the zero-charge member of the rotating charged kink in a spherical sigma model.

Take $r>0$ and $m>0$, the positive-energy sigma-model regime, and keep $\theta =0$ for the energy calculation. For a static field the energy isFor either sign $s=\pm 1$, complete the square:Here $\int {\varphi }_x\sin \varphi dx=-[\cos \varphi {]}_{-\infty }^{+\infty }=2T$. Choose $s$ to have the sign of $T$. The Bogomolny bound is thereforeSaturation requires ${\varphi }_x=sm\sin \varphi$ and $\sin \varphi {\alpha }_x=0$. Finite-energy vacua have $\sin \varphi =0$, so their endpoint cosine values are $\pm 1$ and this particular topological charge is $T\in \{-1,0,1\}$. This is the static kink bound in a mass-deformed spherical sigma model. Negative $r$ would instead make the derivative energy unbounded below, so positivity is an essential physical convention.

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 )$.

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 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.

For the pairwise moment-of-inertia identity, expand the squared pair distances:The last term vanishes in the centre of mass frame, givingLet $C={\sum }_{j<k}{m}_j{m}_k$, and label the two smallest masses $m_{(1)},m_{(2)}$. If $r=\min r_{jk}$, every term in $U$ is at most $G{m}_j{m}_k/r$, so $U\le GC/r$. One pair attains $r$, and its mass product is at least $m_{(1)}{m}_{(2)}$, so $U\ge G{m}_{(1)}{m}_{(2)}/r$. HenceThese constants show the minimum-separation binding-energy bounds: inverse binding energy is comparable to the nearest-pair distance, independently of configuration.

Similarly, every separation is at most $R=\max r_{jk}$ and at least one pair attains $R$. The pairwise identity gives the maximum-separation inertia boundsBoth lower bounds rely on actual extremal separations, rather than an arbitrary numerical lower or upper estimate for all the pair distances.

Negative energy does not imply a positive lower bound on the minimum separation. The printed request is false if it is meant to exclude close approaches or collisions; the trivial bound $r\ge 0$ says nothing of that kind. In fact $E<0$ and $T\ge 0$ give $U=T-E\ge |E|$, which, combined with $rU\le A_0$, proves the negative-energy minimum-separation upper boundThis ensures at least one close pair, not confinement of every particle or prevention of collision.

An explicit negative-energy gravitational collision is a pair initially at rest with separation $a$. Its total energy is $E=-G{m}_1{m}_2/a<0$. For reduced mass $\mu =m_1{m}_2/(m_1+m_2)$, its relative radial equation isThe separation decreases to zero in the finite timeBefore that time the motion is a regular Newtonian solution, yet its separation has no positive infimum. Even at fixed negative energy, bound Kepler orbits with eccentricity approaching one have arbitrarily small pericentre distance. Thus a missing angular momentum or collision-exclusion hypothesis cannot be supplied by the energy sign alone.

For $E>0$, the instantaneous virial theorem gives $\ddot{I}=T+E\ge E$. Integrating twice from any regular time $t_0$ givesUsing $I\le A_1{R}^2$, obtain the positive-energy linear diameter growth boundwhenever the numerator is nonnegative. Consequently, for a solution existing for arbitrarily large future times,This is the precise at-least-linear expansion statement. It does not assert that $R$ is monotone at every instant, nor that all individual particles escape. The large-time conclusion presupposes continued existence of the trajectory; the inequality itself holds on every nonsingular time interval.