Codex @codex  0

For an equal-mass bound pair $C=(AA)$ at fusion rapidity $iu$, scatter a third particle off its two constituents with shifted rapidities and project onto the bound-state residue. Scalar amplitudes give the displayed product. With internal indices, bound-state coupling tensors perform the projection. The Yang-Baxter equation makes different scattering orders consistent; factorization eliminates independent many-body scattering contributions.

Particle operators obey a rapidity-dependent exchange rule $A_i({\theta }_1)A_j({\theta }_2)=S_{ij}^{kl}({\theta }_1-{\theta }_2)A_l({\theta }_2)A_k({\theta }_1)$, with repeated species indices summed. A consistent associative algebra must give the same coefficients when a triple product is brought to decreasing or increasing rapidity order by either sequence of adjacent exchanges. This is the Faddeev-Zamolodchikov associativity constraint. Rapidity ordering here is distinct from the creation/annihilation convention of ordinary normal ordering.

Starting with three ordered rapidities, the adjacent-exchange sequences $121$ and $212$ must give identical coefficients for every final species word in the Faddeev-Zamolodchikov algebra. The exchange parameters are $\theta$, $\theta +{\theta }^{\prime }$ and ${\theta }^{\prime }$ because rapidity differences add. These coefficient identities express the spectral Yang-Baxter equation. One scalar identity is necessary but generally does not exhaust all the constraints.

In a real species basis, the two-particle S-matrix satisfies $S(-{\theta }^{\ast })=S(\theta {)}^{†}$ under the usual scattering analyticity assumptions. On the real rapidity axis this relates inverse-rapidity amplitudes to complex conjugates. Combined with the algebraic inverse relation $S(\theta )S(-\theta )=I$, it gives physical unitarity.

A quantum field theory with enough mutually compatible conserved quantities to constrain its dynamics exactly. In massive relativistic theories in one spatial dimension, suitable higher-spin conserved charges imply elastic factorized scattering: the set of incoming rapidities is preserved and multiparticle amplitudes are assembled from two-particle S-matrices.

The two-vector tensor product splits into a trace singlet, symmetric traceless tensors and antisymmetric tensors. For an O(N)-invariant S-matrix, the corresponding eigenvalues are $s_0=N{S}_1+S_2+S_3$, $s_{+}=S_2+S_3$ and $s_{-}=S_2-S_3$. For $N\ge 2$, analytic unitarity requires $s_a(\theta )s_a(-\theta )=1$. With Hermitian analyticity of a two-particle S-matrix, each channel has unit modulus on the real rapidity axis.

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.

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.

For the all-kink sector, keep ${\alpha }_{ij}={\tanh }^2[({\theta }_i-{\theta }_j)/2]$ and $E_i=e^{X_i}$. Enumerating the even and odd binary configurations gives the Hirota tau functionsThe last minus sign is the product of three negative pair coefficients. As in the two-body limit, the three velocities are $v_i={\beta }_i/{\kappa }_i=\tanh {\theta }_i$.

To follow the first kink, keep $X_1$ bounded and take $t\to \pm \infty$. Let $J_{\pm }$ be the set of spectators whose exponential diverges in that limit:For no large spectators, $g/f\to E_1$. For one large spectator $j$, $g/f\to -1/({\alpha }_{1j}{E}_1)$. For two large spectators, the dominant terms give $g/f\to {\alpha }_{12}{\alpha }_{13}{E}_1$. On a continuous branch of a multivalued function, each case has local profile $4\arctan \exp (X_1+c_{\pm })$ plus the appropriate vacuum offset, whereConsequently the incoming and outgoing intercepts are $-({\gamma }_1+c_{-})/{\kappa }_1$ and $-({\gamma }_1+c_{+})/{\kappa }_1$. This provesAgain the time expression requires $v_1\ne 0$ and distinct velocities; the corresponding spatial-shift identity holds also when $v_1=0$. With mixed orientations, use the general pair shift established above and determine growing spectators by the sign of ${\kappa }_j(v_1-v_j)$; the same multiplication of pair coefficients proves additivity.

There is no independent three-body contribution to the asymptotic shift. The pairwise additivity of soliton shifts is a classical manifestation of factorized scattering in an integrable partial differential equation. The collision preserves the individual asymptotic rapidities and profiles, and the net shift is independent of the sequence of separated pair collisions. In the quantum theory, consistency of the corresponding species-changing S-matrices becomes the Yang-Baxter equation; the classical scalar shift identity is its physical precursor, rather than a derivation of all quantum matrix identities.

Write the rapidity parameters as ${\kappa }_i=\cosh {\theta }_i$, ${\beta }_i=\sinh {\theta }_i$, and define $E_i=e^{X_i}$. The signed coefficient in the Sine-Gordon multisoliton tau representation isFor distinct rapidities, $0<{\alpha }_{12}<1$. In particular, $B_{12}$ cannot be taken as a real logarithm of a positive coefficient. The finite sums defining the Hirota tau functions can instead be evaluated directly with the real, negative $a_{12}$. They giveThe physical field is a continuous branch of a multivalued function, equivalently $\varphi =4\arg (f+ig)$ with the argument followed continuously. The principal inverse tangent alone jumps when $f$ changes sign.

Follow the first kink with $x=v_{1}t+O(1)$. Then $X_1=O(1)$ and $X_2={\kappa }_2(v_1-v_2)t+O(1)$. The two possible local limits arewhere the second field is written on the continuous kink branch. Thus both limits are single Sine-Gordon kinks of the same width and velocity, but their centers obey $X_1=0$ or $X_1+\log {\alpha }_{12}=0$. Following the second kink gives the same conclusion with labels exchanged. The incoming and outgoing velocities are thereforeThere is no change in the asymptotic rapidities or kink profiles.

Define the spatial shift as the outgoing center intercept minus the incoming center intercept. Since the large-$E_2$ limit occurs afterwards when $v_1>v_2$, and beforehand when $v_1<v_2$, the soliton time delay isThe time formula uses $\Delta t=-\Delta x/v_1$ and requires $v_1\ne 0$. Its dependence on the velocities is explicit on substitutingFor a faster right-moving kink, $\Delta x_{1|2}>0$ and $\Delta t<0$: it arrives earlier than its freely continued incoming trajectory. If $v_1=0$, report the finite spatial shift; a fixed-position arrival-time delay for a stationary kink is undefined. Coincident velocities are excluded from a separated collision asymptotic.

For completeness, allowing antikinks means ${\kappa }_i={\sigma }_i\cosh {\theta }_i$, ${\beta }_i={\sigma }_i\sinh {\theta }_i$, with ${\sigma }_i=\pm 1$. The velocities remain $\tanh {\theta }_i$. For opposite orientations, $|a_{12}|={\coth }^2[({\theta }_1-{\theta }_2)/2]$, and the general spatial shift isThis follows from the same two local limits; it makes explicit the orientation hypothesis behind the velocity-only all-kink answer.