Use a common invariant operator domain on which the stated differentiations and commutator are valid. For a differentiable normalized eigenfunction, the Hellmann-Feynman calculation and self-adjointness give
The Lax equation preserves the spectrum. Differentiating now gives , so is in the same eigenspace, or is zero. For a nondegenerate eigenvalue it equals . Multiplying by produces with while preserving its eigenvalue.
Write and . The spatial eigenvalue equation gives , so . For the two first-order systems are therefore
The second row of follows by differentiating the first time equation in and using the spatial equation once more. Conversely these matrix equations imply both scalar equations, so the formulations are equivalent.
Solve , . Since , differentiation gives . To see why the Lax equation preserves the spectrum, direct differentiation gives
Thus is an orthogonal similarity, so all eigenvalues and all coefficients of its characteristic polynomial are first integrals.
A -dimensional Hamiltonian system is Liouville integrable when it has functionally independent first integrals, including the Hamiltonian, whose pairwise Poisson brackets vanish, on an open dense regular region.
For the three-particle periodic Toda lattice, take . The proposed Hamiltonian has
which are exactly the equations of motion. Differentiating the new variables gives
The transformation is not canonical. For example
rather than canonical coordinate brackets. Moreover , so the map forgets the common translation of all the and is not an invertible coordinate change on the six-dimensional phase space.
Multiplication of the displayed matrices gives diagonal entries and off-diagonal entries in . Therefore they form the required Lax pair. The three characteristic coefficients are
They are independent first integrals, as permitted without proof. Equivalently one may take
Indeed , and . Their mutual Poisson commutation also follows directly: generates the common translation, and because is conserved.