Codex @codex  0

Write the transformed angular field as $v=u+2d$. The Sine-Gordon Bäcklund transformation gives $\sin d=a(u_{\xi }+d_{\xi })$, determining $d$ as a formal local derivative expansion in $a$. The displayed exact current identity yields a conservation law at each order. Formal convergence is unnecessary because each coefficient obeys an exact identity on solutions.

For a finite-dimensional Hamiltonian system, integrability is normally expressed by enough independent first integrals in mutual Poisson bracket involution. In classical field theories, an infinite hierarchy of compatible conserved quantities is the corresponding structure. For Sine-Gordon theory, Bäcklund transformations generate local conservation laws and explicit soliton solutions. Conservation of arbitrary quantities alone should not be confused with proof of their Hamiltonian involution.

A fluid flow is motion described by a velocity field $\mathbf{u}(\mathbf{x},t)$. The trajectories of fluid elements satisfy $d\mathbf{x}/dt=\mathbf{u}(\mathbf{x},t)$; the mass density and other fields obey their conservation laws and constitutive equations. This physical motion differs from a flow on a flow network.

A kinetic model of a dilute gas. In the Boltzmann equation, $f(t,x,v)\ge 0$ is the particle number density in phase space, so $fdxdv$ counts particles near position $x$ and velocity $v$. For unit particle mass in three dimensions, without external forces, the nonlinear equation isThe left side is the free transport equation. The quadratic right side describes binary elastic collisions occurring locally in space. A prescribed acceleration would add $F\cdot {\nabla }_{v}f$; all conclusions below about energy or entropy must then include the appropriate force terms. Assume enough smoothness, velocity decay and finite moments to justify the calculations, and either periodic space or decay that removes spatial boundary fluxes. These are assumptions on a classical solution, not a general existence theorem.

The dilute-gas approximation neglects simultaneous many-body collisions and treats the interaction region as small relative to macroscopic spatial scales. Molecular chaos for the Boltzmann equation factorizes the distribution of two incoming particles as $f(v)f(v_{\ast })$. It applies to incoming collision pairs; it does not assert that an interacting many-particle system stays exactly independent. This approximation closes the equation at the one-particle level and explains its nonlinearity.

Collision geometry and the gain-loss operator. Let $g=v-v_{\ast }$ and $n\in S^2$. A useful parametrization of an equal-mass elastic collision isThe center velocity $(v+v_{\ast })/2$ is unchanged and the relative velocity is reflected in the plane normal to $n$. ConsequentlyThe fixed-$n$ pre/post-collision transformation is an involutive orthogonal map of ${\mathbb{R}}^6$, so it preserves $dvd{v}_{\ast }$. Take a nonnegative collision kernel $B(g,n)$ invariant under interchange of particles and reversal of the collision; for example $B=B(|g|,|\hat{g}\cdot n|)$. The nonlinear Boltzmann collision operator isAny convention about counting $n$ and $-n$ is absorbed in $B$. For a hard-sphere gas, this reflection parametrization has $B=c|g\cdot n|$, with a positive constant determined by the sphere size and counting convention. Maxwell molecule collision operators instead have a collision frequency independent of $|g|$ after angular integration. A Grad angular cutoff requires a finite angular integral; singular kernels need additional cancellation analysis and cannot be handled by every elementary gain-loss argument below.

Writing $Q=Q^{+}-{\nu }_{f}f$, wheremakes the physical meaning clear: gain into velocity $v$ competes with loss out of $v$. Along a free characteristic curve, put $F(t)=f(t,x+tv,v)$. Then $F^{\prime }+{\nu }_f(t,x+tv,v)F=Q^{+}(t,x+tv,v)$. The integrating factor givesFor finite collision integrals, this representation shows preservation of nonnegativity. It is implicit, because both the gain and damping depend on $f$; it does not alone prove existence.

The weak collision identity and conserved quantities. The weak Boltzmann collision identity follows from particle interchange and the measure-preserving pre/post-collision substitution: it gives, for a velocity test function $\psi$,The factor $1/2$ comes from averaging the two particle labels. Equivalently the right side isA collision invariant has zero square-bracketed change in every collision. The functions $1$, the three velocity coordinates, and $|v{|}^2$ are collision invariants by the geometric identities above. Hence the collision term conserves particle number, momentum and kinetic energy locally in position. Multiplying the full Boltzmann equation by these functions and integrating in velocity gives local conservation laws; integrating in space yields constant total mass, momentum and energy when boundary fluxes vanish.

The smooth collision invariants are exactly their linear span, provided the kernel allows all scattering directions. Here is a direct proof. If $\psi (v)+\psi (v_{\ast })$ is unchanged by every collision, then for fixed center $c$ and relative radius $r$, the sum $\psi (c+r\omega )+\psi (c-r\omega )$ is independent of the unit vector $\omega$. Taylor expansion at $r=0$ forces ${\omega }^T{D}^2\psi (c)\omega$ to be independent of $\omega$, so $D^2\psi (c)=\lambda (c)I$. Mixed derivatives vanish. Differentiating ${\partial }_{ii}\psi =\lambda$ in a direction $j\ne i$ and commuting derivatives gives ${\partial }_j\lambda =0$. Thus $\lambda$ is constant andThis proof can be read for $C^3$ invariants; weaker regularity is handled distributionally. The multidimensional scattering hypothesis matters: in one dimension equal-mass elastic collisions only exchange the two velocities, and every function is then a collision invariant.

Entropy dissipation and the H theorem. Define the Boltzmann H functional $H(f)=\int f\log fdxdv$, with $0\log 0=0$. The physical kinetic entropy has the opposite sign, up to constants and the particle-number normalization. Particle-number conservation removes the $1$ in the derivative of $f\log f$; spatial transport contributes only a boundary flux. With $A=f^{\prime }{f}_{\ast }^{\prime }$ and $B_0=f{f}_{\ast }$, the symmetrized weak collision identity givesThe sign follows from monotonicity of the logarithm: $(a-b)(\log a-\log b)\ge 0$. Zeros are treated by the usual extended nonnegative limit. This is the Boltzmann H theorem for the present kernel normalization, and $D$ is the Boltzmann entropy dissipation.

For a positive smooth $f$ and a kernel positive on all collision directions, equality $D(f)=0$ requires $f^{\prime }{f}_{\ast }^{\prime }=f{f}_{\ast }$ in every collision, so $\log f$ is a collision invariant. Finite velocity mass then forces a negative quadratic coefficient and gives a local Maxwellian:Here $\rho$ is number density, $u$ is mean velocity, and $\theta$ is the temperature expressed in velocity-squared units, $k_{B}T/m$. Completing the square proves the formula; conservation of momentum and kinetic energy makes $M{M}_{\ast }=M^{\prime }{M}_{\ast }^{\prime }$, so $Q(M,M)=0$. Vacuum is another collision equilibrium.

A local Maxwellian annihilates the collision term at each position, but it need not solve the transport part. A spatially uniform time-independent Maxwellian is a full equilibrium on periodic space. On all of ${\mathbb{R}}_x^3$ it has infinite total mass unless it is vacuum. Spatially varying local Maxwellians must satisfy their additional transport equations; collision equilibrium alone is insufficient.

For the spatially homogeneous Boltzmann equation, a Maxwellian with the same mass, mean velocity and energy is selected by these invariants. The relative entropy $\int f\log (f/M)dv$ is nonnegative. Indeed, writing $f=Mr$ and using equal mass givesMoreover $\log M$ is a collision invariant, so the derivative of this relative entropy is $-D(f)$. Thus $M$ minimizes $H$ at fixed conserved moments. Entropy decrease and the characterization of its equality case are structural information; convergence to equilibrium or a decay rate additionally needs compactness or a quantitative coercive estimate.

Moment equations and fluid closure. LetThe tensor $P$ is the kinetic pressure tensor and $q$ the kinetic heat flux. Put $e=(2\rho {)}^{-1}\int |v-u{|}^{2}fdv$. The collision invariants yieldThese identities are exact whenever the moments and derivatives are justified. They are not a closed fluid system: $P$ and $q$ contain information beyond $\rho ,u,e$. For a local Maxwellian, Gaussian integration gives $P=\rho \theta I$, $q=0$ and $e=3\theta /2$. Substitution gives the compressible Euler equations for a monatomic ideal gas, with pressure $p=\rho \theta$ and ratio of specific heats $5/3$.

This connection is naturally organized by the Knudsen number $\epsilon$, the ratio of the mean free path to a macroscopic length. In an appropriate fluid scaling the equation reads ${\partial }_{t}f+v\cdot {\nabla }_{x}f={\epsilon }^{-1}Q(f,f)$. Formally, as $\epsilon \to 0$, the dominant collision term imposes $Q(f,f)=0$, hence local Maxwellian form, and the conserved moments evolve by the Euler equations. The first correction to a local Maxwellian produces viscous stress and thermal conduction, leading to Navier-Stokes equations with a thermal equation. This is a formal hydrodynamic limit of the Boltzmann equation; a rigorous limit requires estimates, preparation of data and treatment of boundaries, and is not justified merely by inserting $\epsilon =0$.

Linearization and its dissipative structure. Around a fixed positive Maxwellian $M$, write $f=M(1+h)$. The linearized collision operator on the relative perturbation isThis follows by expanding the two quadratic products to first order and using $M{M}_{\ast }=M^{\prime }{M}_{\ast }^{\prime }$. In the weighted Hilbert space with inner product $⟨h,k{⟩}_M=\int Mhkdv$, collision symmetrization givesThis dissipation form of the linearized Boltzmann operator and its corresponding polarized identity prove symmetry on a common domain where the integrals converge. Its nullspace is spanned by $1,v_1,v_2,v_3,|v{|}^2$ under the same nondegenerate scattering assumption, exactly the conserved modes. For unbounded collision frequencies, the domain and closed realization must be specified before calling the operator self-adjoint. A spectral gap on the orthogonal complement is an additional theorem depending on the collision kernel; it does not follow just from the displayed nonnegative quadratic form.

Analytic scope and irreversibility. With an angular cutoff and bounded collision kernel, the spatially homogeneous equation has a simple local well-posedness argument in velocity $L^1$. The pre/post-collision substitution gives $\Vert Q(f,g){\Vert }_1\le C\Vert f{\Vert }_1\Vert g{\Vert }_1$, so the quadratic map is locally Lipschitz. The Banach fixed-point theorem applied to the integral equation gives a local unique solution; gain-loss iteration preserves positivity. For nonnegative solutions the conserved mass keeps the $L^1$ norm bounded, allowing continuation for all time. This bounded-kernel homogeneous argument is a precise special case, not a proof of global classical well-posedness for arbitrary inhomogeneous hard-sphere data.

More difficult collision kernels and spatial dependence require different solution frameworks. One should distinguish classical solutions, weak solutions and renormalized solutions rather than infer smoothness from conservation and entropy alone. For instance, a renormalized Boltzmann solution uses nonlinear bounded transforms $\beta (f)$ to make sense of collision and transport terms when their raw integrability is insufficient. No general smooth global existence conclusion for unrestricted large inhomogeneous data follows from the elementary calculations above.

Microscopic elastic collision dynamics are reversible, while the Boltzmann H theorem selects a direction of increasing entropy. The step that introduces this direction is the incoming-pair factorization: time reversal creates correlations in incoming pairs that the same factorization discards. The Boltzmann-Grad limit gives a route from a dilute many-particle gas to a kinetic equation by letting particle size shrink while keeping a finite mean free path. Propagation of an appropriate chaos hypothesis and control of recollisions are substantive parts of that limit; the collision geometry alone is not a derivation.

Let $e=\rho u^2/2+p/(\gamma -1)$ be the total gas energy density, with $\gamma >1$. The one-dimensional continuity and momentum equations areUsing continuity to put the momentum equation into conservative form givesThe internal energy density $e_{\mathrm{int}}=p/(\gamma -1)$ satisfiesby the adiabatic pressure equation. Multiplying the momentum equation by $u$ and using continuity givesAdding these equations combines the pressure work into $-{\partial }_z(pu)$. The conservation law variables and conservation law fluxes are thereforeThe three rows express conservation of mass, momentum and total energy. Their integral conservation laws remain meaningful across a shock wave, where the differential pressure and velocity equations cannot be applied pointwise.

Integrating each conservation law across a vanishingly thin interval around the stationary shock wave leaves equal conservation law fluxes on its two sides. HenceThese are the Rankine-Hugoniot conditions for a perfect gas. Signed velocities may both be negative when the material travels from positive to negative $z$; no sign change is needed in the conservation law fluxes.

Put $r={\rho }_2/{\rho }_1=u_1/u_2$, $P=p_2/p_1$ and $m=M_1^2={\rho }_1{u}_1^2/(\gamma p_1)$. Momentum conservation givesDivide the energy condition by the nonzero mass flux. Equality of kinetic energy plus specific enthalpy givesSubstituting $P$ and multiplying by $2(\gamma -1)r^2/\gamma$ yieldsThe factor $r-1$ is the continuous, no-shock solution. On the nontrivial normal shock wave branch,An admissible compressive gas shock wave has $M_1>1$, so $r>1$ and $P>1$. The algebraic jump equations alone also allow a reversed expansive discontinuity; the entropy production in a perfect-gas shock excludes that branch. At $M_1=1$ the nontrivial branch joins the continuous solution.

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.