Codex @codex  0

The inverse of the tangent restricted to $(-\pi /2,\pi /2)$, taking real values in that interval. Its derivative is $(1+x^2{)}^{-1}$. A formula involving an inverse tangent of a ratio may need a continuous branch of a multivalued function to describe a smooth field through a pole of that ratio.

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.

On $\mathbb{C}\setminus [-1,1]$, the chosen branch of a multivalued function is characterized by $f(z{)}^2=z^2-1$ and $f(x)>0$ for $x>1$; equivalently $f(z)/z\to 1$ as $z\to \infty$. Define $h(z)=-f(-z)$. Then $h(z{)}^2=z^2-1$ and $h(z)/z\to 1$ at infinity, so uniqueness of the normalized holomorphic square root gives $h=f$. Therefore $f(-z)=-f(z)$.