Codex @codex  0

A Cauchy matrix has entries $1/(x_i+y_j)$ with nonzero denominators. If $x_i=y_i>0$ are distinct, it is positive definite: $C_{ij}={\int }_0^{\infty }{e}^{-x_{i}s}{e}^{-x_{j}s}ds$ is the Gram matrix of linearly independent exponential functions. This form occurs in finite-rank inverse scattering transform reconstructions.

For a field evolving by a Hamiltonian partial differential equation, infinite-dimensional integrability combines sufficiently rich commuting conservation laws with a mechanism such as the inverse scattering transform or a spectral coordinate system that linearizes evolution. Lax equations supply spectral invariants; Lenard-Magri recursions supply commuting Hamiltonians when their solvability conditions hold. Infinite many conserved quantities alone do not prove completeness or an infinite-dimensional action-angle theorem: the function space, boundary conditions and analytic independence all matter.

An integrable partial differential equation has enough exact structure, such as a Lax pair, infinitely many conserved quantities, or an inverse scattering transform, to reduce its evolution to auxiliary linear or finite-dimensional data.

A field $u(x,t)$ has infinitely many degrees of freedom: its initial state is a function rather than a finite list of coordinates. For a Hamiltonian evolution, the analogue of the finite-dimensional construction seeks an infinite sequence of independent conserved functionals $H_j[u]$ whose Poisson brackets vanish. Their Hamiltonian vector fields then define commuting time evolutions. For example, a Poisson operator $J$ produces the field equation $u_t=J\delta H/\delta u$, where $\delta H/\delta u$ is the variational derivative.

A Lax equation packages many conservation laws into one operator identity. If $L_t=[A,L]$, cyclicity of an appropriate trace givesFor a formal pseudodifferential operator the relevant trace is the Adler trace; suitable fractional powers of a monic differential operator produce further conserved functionals. A Lenard-Magri recursion can prove that such functionals are Poisson-commuting functions, rather than merely conserved quantities. A Lax equation by itself does not establish all the independence and completeness properties needed for integrability.

The Korteweg-De Vries equation illustrates the mechanism. In one time normalization take $L={\partial }^2+u$ and $A={\partial }^3+\frac{3}{2}u\partial +\frac{3}{4}{u}_x$. Direct composition givesIts spectral problem, together with the inverse scattering transform for decaying data or suitable spectral coordinates for periodic data, turns nonlinear evolution into simple evolution of spectral data. Solitary waves arise from discrete spectral data, while continuous data describe dispersive radiation.

Thus infinite-dimensional Hamiltonian integrability is more than an infinite list of formulas: one wants sufficiently complete commuting invariants and a reconstruction mechanism that solves the evolution. Boundary conditions and the function space are essential. A formal hierarchy need not converge, and an infinite family of conserved quantities does not automatically give a compact infinite-dimensional torus or an unrestricted analogue of the Arnold-Liouville theorem.

Put $\Psi =\mu e^{-i(kx+2k^{2}t)\sigma }$. The Lax pair becomes ${\Psi }_x=U\Psi$, ${\Psi }_t=V\Psi$, withEquality of mixed derivatives, using invertibility, gives the zero-curvature condition $U_t-V_x+[U,V]=0$. The terms quadratic in $k$ cancel. The terms linear in $k$ cancel because $\sigma Q=-Q\sigma$, $\sigma Q_x=-Q_x\sigma$ and $Q^2$ commutes with $\sigma$. For the constant term, the identitiescancel the differentiated $Q^2$ term, leaving $Q_t+i{Q}_{xx}\sigma -2i{Q}^3\sigma =0$. Multiplying by $i$ provesOnce $S$ is known, its time dependence enters through the explicit exponential conjugation. The boundary relation and normalization form a linear problem for $\mu$, and $Q=i[\sigma ,{\mu }_1]$ reconstructs the nonlinear solution. This is the requested linearization. Its disadvantage as an initial-value method is that it does not itself determine $S$ from a prescribed $Q(x,0)$. The inverse scattering transform supplies that direct-scattering step, including any discrete spectral data and their pole conditions, before reconstruction. Solvability cannot be inferred merely from writing an arbitrary jump.

Use the rapidly decaying real-potential normalization $u_t-6u{u}_x+u_{xxx}=0$. Let $D={\partial }_x$, $L=-D^2+u$ and $A=4D^3-3(uD+Du)=4D^3-6uD-3u_x$. Expanding the commutator cancels its derivative terms and givesThus the KdV equation is the compatibility condition of the Lax pair $L\psi =\lambda \psi$, ${\psi }_t=-A\psi$. The operator $A$ is skew-adjoint on decaying functions. Its evolution transports normalized eigenfunctions unitarily, so the spectrum of $L$ is fixed; alternatively differentiating $L\psi =\lambda \psi$ and using the two equations gives ${\lambda }_t\psi =0$.

At each time solve the KdV Schrodinger spectral problem. For $\lambda =k^2>0$, use the left Jost solution $f_{-}$, normalized as $e^{-ikx}$ at $x\to -\infty$, and write $f_{-}\sim a(k,t)e^{-ikx}+b(k,t)e^{ikx}$ at $+\infty$. The right reflection coefficient is $R=b/a$, and transmission is $1/a$. Because $-A$ tends to $-4D^3$ at infinity, maintaining this Jost normalization gives $f_{-,t}=-A{f}_{-}+4i{k}^3{f}_{-}$. Substitution of its two asymptotic waves yields $a_t=0$ and $b_t=8i{k}^{3}b$. The negative discrete eigenvalues are $-{\kappa }_j^2$, fixed in time. A unit-normalized bound state with right tail $c_j(t)e^{-{\kappa }_{j}x}$ obeys $c_j^{\prime }=4{\kappa }_j^3{c}_j$. Consequently the Time evolution of KdV scattering data isThese signs follow from the specified right-reflection convention; reversing which infinity defines reflection reverses the continuous-data phase convention.

The inverse scattering transform now solves the nonlinear evolution in three steps. Compute $R(k,0)$, the discrete eigenvalues and the norming constants from $u(x,0)$; evolve this KdV scattering data by the simple formulas above; reconstruct $u$ by the Gelfand-Levitan-Marchenko equation. Explicitly setand solve, for $y\ge x$,The recovered potential is $\boxed{u(x,t)=-2dK(x,x,t)/dx}$. Each fixed discrete eigenvalue produces a soliton component, whereas the continuous scattering data encode the dispersive part. The inverse transform, rather than a linear superposition of the potentials, accounts for their nonlinear interaction.

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.

The Korteweg-De Vries equation is the compatibility condition $L_t=[P,L]$ for the Lax pairFor rapidly decaying initial data, first solve the direct KdV Schrodinger spectral problemComparing its Jost solutions at the two spatial infinities produces the KdV scattering data: the reflection coefficient $R(k)$ on the continuous spectrum, discrete eigenvalues $-{\kappa }_n^2$, and norming amplitudes $c_n$ for the bound states.

The Lax evolution is isospectral. More precisely, the Time evolution of KdV scattering data isThis is the simple linear evolution that makes the nonlinear initial-value problem tractable.

To invert the data, formSolve the Gelfand-Levitan-Marchenko equationand reconstructDirect scattering at $t=0$, linear evolution of the data, and this inverse step constitute the inverse scattering transform.

Putand seek the rank-one Gelfand-Levitan-Marchenko kernel in the separable formSubstitution into the Gelfand-Levitan-Marchenko equation givessoand hence

WriteThe reconstruction formula for this inverse scattering transform convention yieldsIf $q=e^{-2s}$, then $4q/(1+q{)}^2={\mathrm{sech}}^{2}s$. Takingtherefore produces the required one-soliton of the Korteweg-De Vries equation, withHere the usual bound-state norming constant has ${\beta }_0>0$ and $\chi >0$, so $\varphi$ is real.

Scattering data describe connections between normalized solutions of an auxiliary spectral equation, together with any discrete spectral points and their norming data. In an inverse scattering transform, the nonlinear potential is reconstructed from these data; boundary spectral connection matrices must additionally obey a global relation.

In light-cone coordinates, the sine-Gordon equation can be written $u_{XT}=\sin u$. Its inverse scattering transform, solitons, breathers, and scaling symmetries make it an integrable system.