Codex @codex  0

Multiparticle quantum scattering is factorized when it decomposes into two-body scatterings without particle production or a new irreducible many-body amplitude. The independence of the answer from the order of exchanges imposes the Yang-Baxter equation. Classical pairwise additivity of soliton shifts is a corresponding property of integrable soliton collisions.

A narrow outgoing energy packet has phase $-ET/\hslash +\delta (E)$, so stationary phase gives the arrival shift $\Delta T=\hslash {\delta }^{\prime }(E)$. For two equal-mass solitons with relative rapidity $\theta$, $E=2\mathcal{M}\cosh (\theta /2)$. The convention $S=e^{2i{\delta }_{\mathrm{pw}}}$ instead gives $2\hslash {\delta }_{\mathrm{pw}}^{\prime }=\Delta T$. Classical time delays determine only the energy-dependent phase difference, not its constant branch.

An integer-order sech initial pulse in the ideal focusing cubic fibre equation has periodic shape evolution involving compression and recombination. It does not preserve a single fixed profile like a fundamental temporal bright soliton. Arbitrary noninteger-order initial pulses can contain radiation as well as solitons.

An exact-solution method that transforms a nonlinear partial differential equation into equations bilinear in auxiliary Hirota tau functions. Finite exponential expansions of the Hirota tau functions produce soliton families. The auxiliary exponential coefficients need not all be positive.

Compare the same total number of finite-box modes in the soliton and vacuum sectors before removing the ultraviolet cutoff. A localized bound state replaces a continuum mode, and its contribution must be retained. For the Sine-Gordon kink fluctuation operator and the odd phase branch $\delta (k)=2\arctan (1/k)$, the bound translation mode replaces the free $k=0$ oscillator. Thus the complete oscillator correction is $-m/2-(m/2)\int dk\delta (k)k/[2\pi \sqrt{k^2+1}]$ in this convention. The continuum integral alone misses the finite $-m/2$ contribution.

An exact nonlinear-wave solution containing several localized solitons. In an integrable scattering solution, these objects approach separated incoming and outgoing profiles, with collision effects encoded by position shifts and internal phases. Bäcklund transformations and Hirota's bilinear method construct examples such as the Sine-Gordon multisoliton tau representation.

A many-soliton spatial shift is pairwise additive when $\Delta x_i={\sum }_{j\ne i}\Delta x_{i|j}$. In the all-kink sector of Sine-Gordon theory, $\Delta x_{i|j}=-\mathrm{sgn}(v_i-v_j)\log {\tanh }^2[({\theta }_i-{\theta }_j)/2]/\cosh {\theta }_i$. Dominant spectator exponentials in the Sine-Gordon multisoliton tau representation multiply the effective exponential of kink $i$, so their logarithms add. This is a classical signature of factorized scattering, with no independent many-body contribution to the asymptotic shift.

For a single spectral mass, the algebraic equation givesDifferentiate the reconstruction, writing $r=cE/(2\kappa )$:With $x_0={\kappa }^{-1}\log [c/(2\kappa )]$, the concise one-soliton answer isIts speed is ${\kappa }^2$ and its maximum is ${\kappa }^2/2$. Equivalently, setting $\eta =\kappa /2$ gives $q=2{\eta }^2{\mathrm{sech}}^2[\eta (x-4{\eta }^{2}t-x_0)]$. To check the sign and speed independently, this profile obeys $q_{xx}={\kappa }^{2}q-3q^2$ and $q_t=-{\kappa }^2{q}_x$; differentiation gives the stated KdV equation exactly.

When $b=0$, put $D=e^{-2px+8p^{3}t-2x_0}>0$. Since $c=2ip{e}^{-2x_0}$, evaluation at $ip$ gives $N(ip)=1-DN(ip)$, henceDifferentiate the coefficient of $1/k$ with respect to $x$, using $D_x=-2pD$, and apply the reconstruction formula. The result isThe positive sign follows from the source's $-2i$ reconstruction convention. The soliton has amplitude $2p^2$, speed $4p^2$, and centre $x=4p^{2}t-x_0/p$.

A simple zero at $ip$ supplies a discrete scattering mode and hence an order-one, right-moving soliton component of speed $4p^2$. Its position depends on the associated norming constant and, for general scattering data, the appropriate scattering phase shifts. A simple zero alone does not specify its position, nor exclude additional zeros and therefore additional solitons.

To locate the continuous radiation, set $\xi =x/t$. The real-axis oscillatory phase is $t(2k\xi +8k^3)$, whose stationary points satisfy $2\xi +24k^2=0$. For $\xi <0$ there are two real stationary points $k=\pm \sqrt{-\xi /12}$, giving dispersive oscillatory radiation, ordinarily of order $t^{-1/2}$ for nondegenerate data. For $\xi >0$ there is no real stationary point, so radiation is smaller and a soliton dominates near its moving centre. At $\xi =0$ the stationary points merge and the natural transition scale is $x=O(t^{1/3})$, with Airy-type radiation of order $t^{-1/3}$. These statements assume the regular decaying scattering data implicit in the reconstruction problem; vanishing reflection can remove the radiation entirely. The discrete mode travels right, while the leading continuous radiation occupies left-going rays.

Use the Fourier convention ${\hat{q}}_0(k)={\int }_{\mathbb{R}}{e}^{-ikx}{q}_0(x)dx$ and sufficiently rapidly decaying smooth initial data. A scalar Lax pair for the linear equation isDifferentiate the first equation in time and use the second and its spatial derivative. The terms involving $m$ and $k$ cancel, leaving $q_t=i{q}_{xx}$, exactly the linear equation.

The spatially normalized solutions are $m^{+}={\int }_{-\infty }^x{e}^{ik(x-y)}q(y,t)dy$ for $\mathrm{Im}k>0$, and $m^{-}=-{\int }_x^{\infty }{e}^{ik(x-y)}q(y,t)dy$ for $\mathrm{Im}k<0$. Integration by parts gives $m^{\pm }=O(k^{-1})$. Their difference is $e^{ikx}\hat{q}(k,t)$; Fourier transformation of the evolution gives $\hat{q}(k,t)=e^{-i{k}^{2}t}{\hat{q}}_0(k)$. Set $M^{\pm }=\left(\begin{array}{cc}1& m^{\pm }\\ 0& 1\end{array}\right)$. ThenThe Cauchy integral solution is $m(k)=(2\pi i{)}^{-1}\int e^{isx-i{s}^{2}t}{\hat{q}}_0(s)/(s-k)ds$, and reconstruction $q=-i{\lim }_{k\to \infty }km(k)$ gives the usual inverse Fourier solution.

For the nonlinear deformation, put $D=\mathrm{diag}(1,-1)$ and $Q=\left(\begin{array}{cc}0& q\\ \vartheta & 0\end{array}\right)$. The matrix Lax pair isA direct multiplication gives the off-diagonal entries of $U_t-V_x+[U,V]$ as $q_t-i{q}_{xx}+2i{q}^2\vartheta$ and ${\vartheta }_t+i{\vartheta }_{xx}-2i{\vartheta }^{2}q$, with zero diagonal. Hence compatibility is precisely the two requested nonlinear equations.

Here is a normalized Riemann-Hilbert problem for this pair, in the case without discrete scattering poles. Let $a(k),b(k)$ be its two reflection data determined from the initial spatial scattering problem, and set $\theta =kx-k^{2}t$. Deform the triangular jump toThe data $a,b$ are time-independent in this convention; their spatial and temporal oscillations are displayed explicitly. Setting $b=0$ and $a={\hat{q}}_0$ recovers the linear problem. Thus the diagonal correction $-ab$ is necessary when the second off-diagonal jump is introduced; omitting it would destroy determinant one.

To verify the deformation rather than merely name it, note $J_x=(ik/2)[D,J]$ and $J_t=-(i{k}^2/2)[D,J]$. Therefore $(M_x-(ik/2)[D,M])M^{-1}$ and $(M_t+(i{k}^2/2)[D,M])M^{-1}$ have no jump. Analyticity and the expansion at infinity make them, by Liouville's theorem, respectively $Q$ and $-kQ+V_0$. The first expansion gives $Q=-(i/2)[D,M_1]$, soThe constant term of the second expression is $V_0=(M_1{)}_x$. In the first equation's next coefficient, its off-diagonal part equals $iD{Q}_x$ and its diagonal part equals $-iD{Q}^2$. Hence $V_0=iD(Q_x-Q^2)$, exactly the displayed nonlinear Lax pair. This proves that the reconstructed functions obey the required PDEs.

Finally let $C$ be the Cauchy integral operator on the real axis, $C_{-}$ its lower boundary value, and $\mu =M^{-}$. The jump relation and normalization giveThis is the requested linear singular integral equation for the inverse problem: it is linear in the unknown matrix $\mu$ once the scattering data are given. Reconstruction is $q=(2\pi {)}^{-1}\int [\mu (J-I){]}_{12}ds$ and $\vartheta =-(2\pi {)}^{-1}\int [\mu (J-I){]}_{21}ds$.

For completeness, generic simple discrete eigenvalues add the following residue data. At an upper-half-plane pole $k_j$, use $C_j=\left(\begin{array}{cc}0& c_j{e}^{i\theta (k_j)}\\ 0& 0\end{array}\right)$; at a lower-half-plane pole use $C_j=\left(\begin{array}{cc}0& 0\\ d_j{e}^{-i\theta (k_j)}& 0\end{array}\right)$, with the norming constants determined by the initial spatial scattering problem. Impose ${\mathrm{Res}}_{k_j}M={\lim }_{k\to k_j}M{C}_j$. The nilpotent matrices make this limit well-defined. If $A_j$ denotes the residue, the full linear inverse system isThe reconstructed field formulas acquire $-i{\sum }_j(A_j{)}_{12}$ and $i{\sum }_j(A_j{)}_{21}$ respectively. This includes the soliton terms when they occur; the preceding pure-jump formula is their pole-free special case. These standard scattering formulations assume no spectral singularity on the real contour and simple discrete eigenvalues; repeated poles require the corresponding higher-order pole conditions.

Insert the finite exponential expansion into the Gelfand-Levitan-Marchenko equation. The integral contribution isComparing coefficients of $e^{-c_{n}y}$ givesFor this coefficient comparison take the $c_n$ distinct; repeated exponents can first be combined. The linear system supplies the coefficients whenever its matrix is invertible.

In the inverse scattering transform convention $u_t-6u{u}_x+u_{xxx}=0$, the Korteweg-De Vries equation solution is $u=-2dK(x,x)/dx$. The reflectionless scattering data evolve as ${\beta }_n(t)={\beta }_n(0)e^{8c_n^{3}t}$, with the $c_n$ fixed. For one exponential,Differentiating its diagonal givesFor $\beta >0$, this is the smooth solitonIt travels at speed $4c^2$. If $\beta =0$ the solution is zero; if $\beta <0$ the denominator vanishes somewhere, so the formula is not a globally smooth one-soliton. Positivity is the necessary qualification for a regular soliton.

Use $\hslash =c=1$ in the field-theory calculations. Work with signature $(+,-)$ and $a\ne 0$. Since $x=(x_{+}+x_{-})$ and $t=(x_{+}-x_{-})$, ${\partial }_{+}={\partial }_x+{\partial }_t$, ${\partial }_{-}={\partial }_x-{\partial }_t$, and the Sine-Gordon equation is ${\varphi }_{+-}=\sin \varphi$. Introduce $w=(\rho +\kappa )/2$ and $d=(\rho -\kappa )/2$. The specified Sine-Gordon Bäcklund transformation becomes $d_{+}=a\sin w$, $w_{-}=a^{-1}\sin d$. Its compatibility givesAdding and subtracting these equations establishesThus a compatible transformed field satisfies the same Sine-Gordon equation. Smoothness is needed to interchange mixed derivatives. The statement cannot literally include $a=0$, since one defining equation contains $1/a$.

Generating the kink. With zero seed, the two first-order equations give ${\rho }_{+}=2a\sin (\rho /2)$ and ${\rho }_{-}=2a^{-1}\sin (\rho /2)$. On a nonconstant branch put $p=\tan (\rho /4)$; then $p_{+}=ap$, $p_{-}=a^{-1}p$. Integration givesFor $a>0,C>0$, this is the Sine-Gordon kink, after absorbing $C$ in the center position. Its speed and Lorentz factor obeyChanging the sign of $a$ or of the exponential coefficient selects the corresponding antikink orientation; $|v|<1$ whenever $a$ is finite and nonzero. Constant vacuum branches, omitted by division by $\sin (\rho /2)$, can be added separately.

The continuously parametrized Bäcklund transformation does more than produce this one solution. Successive compatible transformations and Bianchi permutability for sine-Gordon Bäcklund transformations construct multisoliton solutions; expansion of the generating relations gives the local conserved-charge hierarchy of sine-Gordon theory. Together with its Lax pair and inverse scattering transform, this is the structure behind classical integrability, elastic scattering and the absence of generic radiative energy loss in soliton collisions. Existence of a transformation alone is not a proof that every generated conserved charge is independent and in involution; those are additional properties of the integrable system.

Transforming the static kink at $a=1$. The seed obeys ${\kappa }_x=2\mathrm{sech}x$, $\sin (\kappa /2)=\mathrm{sech}x$, $\cos (\kappa /2)=-\tanh x$. Adding and subtracting the two Sine-Gordon Bäcklund transformation equations yieldsAgain set $p=\tan (\rho /4)$. These simplify to $p_x=-p\tanh x$, $p_t=\mathrm{sech}x$. The first gives $p=f(t)/\cosh x$, and the second gives $f^{\prime }=1$. ThereforeContinuous inverse-tangent branches may be chosen without creating artificial jumps. At spatial infinity the field approaches the same vacuum, so its net topological charge is zero. For large $|t-t_0|$ there are two well-separated transitions near $x_{\pm }=\pm \mathrm{arcosh}|t-t_0|$, with speeds tending to zero. The two transitions exchange their kink/antikink orientations as they pass through the collision, while the field remains smooth.

This is the Sine-Gordon threshold kink-antikink solution, a separatrix between finite-speed scattering and the bounded Sine-Gordon breather motion, rather than a finite-period breather. In fact it is the $v\to 0$ limit of the scattering field in question 3. At $t=t_0$, $\rho =0$ and ${\rho }_t=4\mathrm{sech}x$, so the reduced energy is $\int {\rho }_t^{2}dx/2=16$, exactly twice the reduced single-kink energy $8$.

For mass comparisons below, the original PDF really prints $m^2/\beta$. Restoring physical coordinates $X=x/m$, $T=t/m$ gives the angular-field action ${\beta }^{-1}\int dtdx[\frac{1}{2}(\partial \varphi {)}^2+\cos \varphi -1]$ and literal kink mass $8m/\beta$, for $\beta >0$. The usual Sine-Gordon theory instead has prefactor $m^2/{\beta }^2$ and mass $8m/{\beta }^2$. The equation of motion and every field above are unchanged by this overall normalization; the quantum scattering comparison in question 3 is not.

Take $\lambda \ne 0$ and use the complex sine-Gordon theory reduced density $\ell =\mathcal{L}/M^2$, with physical lengths and times $x/M,t/M$. Physical rest energy is $M\int dx{\mathcal{H}}_{\ell }$, and the dimensionless Noether charge is obtained from the reduced action. Put $c=\cos \alpha$, $s=\sin \alpha$, $f(x)=c\mathrm{sech}(cx)/\lambda$, so $\psi =f{e}^{ist}$. Initially take $0<|\alpha |<\pi /2$, avoiding the singular target-space boundary.

The canonical Hamiltonian density isHere $|\dot{\psi }{|}^2=s^2{f}^2$, $|{\psi }_x{|}^2=c^2{f}^2{\tanh }^2(cx)$, and $D=1-c^2{\mathrm{sech}}^2(cx)$. Since $s^2+c^2{\tanh }^2(cx)=D$, the kinetic contribution is simply $f^2$. ConsequentlyTo fix the charge sign, take the Noether current directly from $\delta \psi =i\epsilon \psi$, $\delta {\psi }^{\ast }=-i\epsilon {\psi }^{\ast }$:Thus positive internal rotation has positive charge in this convention. With $X=cx$ and the supplied integral,Reversing the definition of the generator would reverse all charge labels but not any mass. The limit $\alpha \to \pi /2$ is the vacuum with zero mass and charge. At $\alpha =0$ the field reaches $|\psi |=1/|\lambda |$ at its center and the current's coordinate expression becomes singular. The mass still has the finite limiting value $4M/{\lambda }^2$, butThe displayed local density does not assign a unique endpoint charge merely by substituting $\alpha =0$.

The four collective parameters. Translation, a Lorentz boost of rapidity $b$, and the global phase symmetry produceThe independent real parameters are $\alpha ,b,X_0,\chi$, with $\chi$ periodic modulo $2\pi$. A further time translation is absorbed into these parameters. The energy and momentum are $\mathcal{M}\cosh b$ and $\mathcal{M}\sinh b$, while $Q$ is invariant. The internal angle is the collective coordinate whose canonical momentum is this Noether charge.

Semiclassical quantization. For one rest-frame internal rotation, the period in dimensionless time is $T=2\pi /|s|$. The action around the orbit isEquivalently this is $S+\mathcal{E}T$, where $\mathcal{E}$ is the reduced rest energy. Bohr-Sommerfeld quantization with $\hslash =1$ gives $I=2\pi n$. Because the internal orbit is a cyclic angle, not a one-dimensional oscillator between turning points, there is no half-integer Maslov index shift. HenceOn the nonsingular branch $0<n<2\pi /{\lambda }^2$. The elementary charge-one mass approaches $M$ at weak coupling. For positive $a,b$ with $(a+b){\lambda }^2/4\le \pi /2$, $\sin a^{\prime }+\sin b^{\prime }>\sin (a^{\prime }+b^{\prime })$, where $a^{\prime }=a{\lambda }^2/4$, $b^{\prime }=b{\lambda }^2/4$. Thus ${\mathcal{M}}_{a+b}<{\mathcal{M}}_a+{\mathcal{M}}_b$: the charge-$n$ soliton cannot fragment into smaller like-charge particles. For a general charge-conserving partition, a constituent with $|q_i|\ge n$ already has mass at least ${\mathcal{M}}_n$; if every $|q_i|<n$, the decreasing ratio ${\mathcal{M}}_q/q$ gives ${\sum }_i{\mathcal{M}}_{|q_i|}\ge ({\mathcal{M}}_n/n){\sum }_i|q_i|\ge {\mathcal{M}}_n$, strictly for nontrivial fragmentation. Thus opposite-charge constituents cannot open a lower-energy decay either. This is charged rather than topological stability, and the semiclassical spectrum is reliable at weak coupling.

Counting states and the endpoint. Write $L=2\pi /{\lambda }^2$. If $L$ is not an integer, the regular branch givesIf $L$ is an integer, the local nonsingular calculation gives only $2(L-1)$ regular states. The formal closed-branch list has $2L$ charge labels, but its last two labels are limits of the same singular static solution. Counting them needs a global completion; it cannot be decided from a regular current integral at $\alpha =0$.

In the additional integer-level completion, ${\lambda }^2=4\pi /k$ at leading semiclassical order and charges are identified modulo integer $k>1$. There are then $k-1$ nonzero charge sectors. For odd $k$, positive and negative labels are distinct; for even $k$, the endpoint labels $\pm k/2$ describe one self-conjugate state. This supplies $2L-1$ states when $k=2L$, rather than $2L$ endpoint labels. This global interpretation, and the distinction between classical and renormalized couplings, are discussed in Dorey and Hollowood, section 2. A generic classical $\lambda$ does not automatically define an integer $k$ or specify this quantum completion.

Choose the positive-exponential branches from part (i) and set their additive constants to zero. First take $0<b<1$, and defineThen $\tan ({\varphi }_b/4)=e^{\gamma (x-vt)}$ and $\tan ({\varphi }_a/4)=e^{-\gamma (x+vt)}$. The tangent subtraction formula givesConsequently the allowed Sine-Gordon superposition formula produces the smooth fieldThis is the negative of the Sine-Gordon two-kink solution, and hence a two-antikink configuration. The auxiliary seeds have opposite topological charges, but their charges cannot simply be added to infer the charge of the nonlinear two-step Bäcklund transformation. Indeed, the displayed final field tends to $2\pi$ at the left spatial end and $-2\pi$ at the right, so its total topological charge is $-2$.

Let $T=|t|$ become large. Near the right transition, $x=vT+O(1)$, the tangent argument has the asymptotic formso the local field is $-4\arctan e^{\gamma (x-vT)+\log v}$, a single antikink. Near the left transition, the local field is $4\arctan e^{-\gamma (x+vT)+\log v}+o(1)$, again a decreasing antikink. The resulting asymptotic center lines areThus two incoming antikinks with topological charges $(-1,-1)$ and velocities $(+v,-v)$ separate again with exactly the same topological charges and velocities. There is no radiative tail in these asymptotic profiles. Labeling the outgoing objects by their preserved rapidities makes this elastic soliton scattering; labeling the left and right lumps instead describes reflection with exchanged velocities.

For the right-moving soliton, its incoming intercept is ${\gamma }^{-1}\log v$ and its outgoing intercept is $-{\gamma }^{-1}\log v$. The spatial shifts are therefore $\Delta x_{+}=-2{\gamma }^{-1}\log v$ and $\Delta x_{-}=2{\gamma }^{-1}\log v$. Define the soliton time delay as the change in arrival time at a fixed distant spatial point relative to continuation of the incoming straight line, so $\Delta t=-\Delta x/v_{\mathrm{particle}}$. Both objects have the same signed soliton time delay, which is an advance:This is the Sine-Gordon two-kink time advance. In physical coordinates $T_{\mathrm{phys}}=t/m$, the time shift is $2\log v/(m\gamma v)$. The explicit intercepts fix the sign convention unambiguously.

The remaining real parameter choices are covered without changing the calculation. For any $b\ne 0$ with $b^2\ne 1$, put $v_b=(1-b^2)/(1+b^2)$, $u=|v_b|$, $\gamma =(|b|+|b{|}^{-1})/2$ and $ϵ=\mathrm{sgn}(b{v}_b)$. The same choice of zero additive constants givesEach scattered object's topological charge is $-ϵ$, the velocities are $\pm u$, and the signed soliton time delay is $2\log u/(\gamma u)$. If $b=\pm 1$, the superposition coefficient vanishes and this representative is the vacuum; there is no pair of separated moving solitons and no scattering delay to assign. Thus the scattering conclusion requires the nondegenerate case $0<u<1$.

With zero seed, both Bäcklund transformation equations concern ${\varphi }_a/2$. On a nonconstant branch, separation of variables uses $d\log |\tan ({\varphi }_a/4)|=d{\varphi }_a/[2\sin ({\varphi }_a/2)]$, givingFor $C>0$, absorb its magnitude into an additive constant $c=\log C$ in the exponent. A convenient smooth representative isThe constant-phase condition for this traveling profile determines its velocity:The profile is a traveling soliton with velocity $v_a=(1-a^2)/(1+a^2)$, strictly between $-1$ and $1$. The scalar-field vacua on the two sides differ by $2\pi$. With the topological charge convention $Q=[\varphi (+\infty )-\varphi (-\infty )]/(2\pi )$, this branch has $Q=\mathrm{sgn}a$: it is a Sine-Gordon kink for $a>0$ and an antikink for $a<0$. Its width is proportional to ${\gamma }_a^{-1}$, and its derivative decays exponentially away from its center, giving a finite-energy field configuration.

For completeness, the classical rest mass in this paper's coupling convention is $8m/\beta$, as derived in Question 2; a Lorentz boost gives energy $8m{\gamma }_a/\beta$. This localized, topologically protected traveling field is the required classical field-theory soliton. A negative $C$ reverses the field and hence the topological charge; $C=0$ gives the vacuum rather than a soliton. Vacuum shifts by $4\pi$ can be made without changing the displayed Bäcklund transformation equations.

Use the Minkowski metric with signature $(+,-)$, so the Sine-Gordon equation is ${\varphi }_{tt}-{\varphi }_{xx}+\sin \varphi =0$. The light-cone coordinates here satisfyIn particular, there is no extra factor of four with this coordinate normalization. The Bäcklund transformation requires $a\ne 0$, because one of its equations contains $a^{-1}$.

Put $F=(\psi +\varphi )/2$ and $G=(\psi -\varphi )/2$. In the parameter convention of this paper the Sine-Gordon Bäcklund transformation gives $G_{+}=a\sin F$ and $F_{-}=a^{-1}\sin G$. Differentiating, for a twice differentiable transformed field, givesSince $\psi =F+G$, the addition formula yields ${\psi }_{+-}=\sin (F+G)=\sin \psi$. Therefore the transformed field satisfies the same Sine-Gordon equation:Subtracting the two differentiated equations also gives ${\varphi }_{+-}=\sin \varphi$, showing the compatibility with the seed equation. This proves the unheaded request before the numbered parts, without assuming a particular soliton form.

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.