Codex @codex  0

A Dirac field theory invariant under $\psi \mapsto e^{i\alpha }\psi$, $\bar{\psi }\mapsto e^{-i\alpha }\bar{\psi }$ has a conserved global charge, subject to the quantum theory preserving the symmetry. In particle language this charge counts particles minus antiparticles. When the S-matrix commutes with it, a scattering amplitude between states of different charge vanishes. A neutral scalar and a Yukawa interaction preserve the charge, so a process with one incoming Dirac particle and only one outgoing Dirac antiparticle plus neutral particles is forbidden. This net charge is distinct from a sum of positive occupation numbers for all particle and antiparticle modes.

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.

For particles in the vector representation of the orthogonal group, the two-body S-matrix is a linear combination of the identity $I$, permutation $P$ and trace contraction $K$: $S=S_{2}I+S_{3}P+S_{1}K$. Here $K(e_i\otimes e_j)={\delta }_{ij}{\sum }_a{e}_a\otimes e_a$, $P^2=I$, $PK=KP=K$ and $K^2=NK$. The orthogonal-invariant scattering channels diagonalize these three operators simultaneously.

There are exactly two tree Feynman diagrams: the incident scalar can be absorbed before the final scalar is emitted, or after it. The internal fermion four-momenta are respectively $r_s=p+k=q+l$ and $r_u=p-l=q-k$. A scalar exchange diagram would require an absent scalar self-interaction, so there is no additional tree channel.

Define the S-matrix convention $S_{fi}={\delta }_{fi}+i(2\pi {)}^4{\delta }^{(4)}(p+k-q-l)T_{s^{\prime }s}$. The two tree scattering amplitudes areUsing the Dirac propagator from the preceding Feynman rules,The two terms add with the same relative sign: neither diagram exchanges identical external fermions. Their interference must be retained when squaring the complete scattering amplitude. In the phase convention above, ${\gamma }^5(M-i/r){\gamma }^5=-(M+i/r)$, which gives an equivalent simplified numerator. An overall phase depends on the S-matrix convention and does not change a relativistic scattering cross-section.

The Dirac field has the global symmetry $\psi \mapsto e^{i\alpha }\psi$, $\bar{\psi }\mapsto e^{-i\alpha }\bar{\psi }$, while the real scalar is unchanged. Both the free Dirac action and the pseudoscalar Yukawa interaction preserve this symmetry. Therefore Dirac fermion number conservation holds: its charge counts particles minus antiparticles.

The initial state has charge $+1$, whereas a final antiparticle and a neutral scalar have charge $-1$. Since the S-matrix commutes with that charge,This holds at every order, not just tree level. In the Feynman rules, a continuous fermion arrow cannot connect these specified external states. Replacing only the outgoing $\bar{u}$ by a $v$ in the previous expression would not give a physical amplitude. Changing both external fermions to antiparticles would instead produce an allowed process with reversed fermion flow and the appropriate $v$ spinors; moving a leg between initial and final states is a different operation governed by crossing symmetry.

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.

Use the momentum-space Feynman rules with the scalar Feynman propagator $i/(p^2-m^2+i0)$. The four-leg vertex of a factorial-normalized scalar interaction has weight $i\lambda$ for the positive interaction sign printed here. There are $4!$ Wick contractions assigning four external legs to its four fields, cancelling the factorial in its coefficient. A negative interaction sign would give $-i\lambda$; its squared tree amplitude is the same.

At each vertex include $(2\pi {)}^4{\delta }^4(\sum p)$ with all incident momenta taken incoming. Assign an internal momentum to each line and integrate each independent loop with $\int d^4\ell /(2\pi {)}^4$. Divide a graph by its Feynman-diagram symmetry factor, sum the graphs at the chosen order, and omit disconnected vacuum graphs from normalized amplitudes. For an S-matrix element, amputate external propagators and put the external momenta on shell as in the LSZ reduction formula; the external one-particle residues are one at tree level.

Define the invariant amplitude by the relativistically normalized matrix elementwith $⟨p|p^{\prime }⟩=2E_{\mathbf{p}}(2\pi {)}^3{\delta }^3(\mathbf{p}-{\mathbf{p}}^{\prime })$. The lowest-order connected four-point graph is one contact vertex:There is no exchange graph at this order because there is no three-field interaction.

For the elastic scattering from a quartic scalar contact interaction, write $\sqrt{s}$ for the total centre-of-mass energy and $E_p=\sqrt{s}/2$ for the energy of each incoming particle. The incoming and outgoing spatial momentum magnitudes both equal $k=\sqrt{E_p^2-m^2}$, with $E_p>m$. The invariant incident flux isThe Lorentz-invariant phase-space measure for two outgoing particles isIn the centre-of-mass frame, the spatial delta function sets ${\mathbf{p}}_4=-{\mathbf{p}}_3$, while the energy delta function has radial derivative $2k/E_p$. Therefore the relativistic two-body phase space satisfiesThe two outgoing real-scalar particles are identical. Integrating over the full solid angle counts each unordered pair twice, so include the identical final-state symmetry factor $1/2!$. This givesThis is isotropic. When $E$ denotes each particle's energy, $s=4E^2$ and the full-sphere event density is ${\lambda }^2/(512{\pi }^2{E}^2)$. If $E$ denotes the total energy of the pair, $s=E^2$ and it is ${\lambda }^2/(128{\pi }^2{E}^2)$. Stating the answer in $s$ removes that energy-label ambiguity.

An equally valid angular convention selects one outgoing particle in a hemisphere, so each event is represented once. In that convention omit $1/2!$ and use $d\sigma /d\Omega ={\lambda }^2/(64{\pi }^{2}s)$ on the hemisphere, or ${\lambda }^2/(256{\pi }^2{E}_p^2)$. Both conventions give ${\sigma }_{\mathrm{event}}={\lambda }^2/(32\pi s)$ at this order. The identical-state factor concerns counting final states and is separate from the vertex factorial.

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.

For one spatial dimension and $c=1$, a massive particle has $v=\tanh \theta$, $E=M\cosh \theta$ and $p=M\sinh \theta$. A Lorentz boost adds a constant to all rapidities, so a rapidity difference is invariant. Rapidity is useful for factorized scattering, where two-body S-matrices depend only on that difference.

The S-matrix maps incoming asymptotic states to outgoing asymptotic states. Its connected transition matrix elements contain an overall energy-momentum Dirac delta function and an invariant scattering amplitude. In quantum field theory, the LSZ reduction formula obtains these elements from amputated field correlators.