Lax pairs and isospectral flows

[10/80%]
Show
LaTeX

Let LL and KK be two n×nn \times n matrices. Assume that the entries depend on a parameter tt and are differentiable with respect to tt. Assume that K−1(t)K^{-1}(t) exists for all tt. Assume that the time-evolution of LL is given by

L(t)=K(t)L(0)K−1(t). L(t) = K(t)L(0)K^{-1}(t) .
?
(i)

Show that L(t)L(t) satisfies the matrix differential equation

dLdt=[L,B](t) \frac{dL}{dt} = [L, B](t)

where [ , ][\, , \, ] denotes the commutator and

B=−dKdtK−1(t). B = -\frac{dK}{dt}K^{-1}(t) .
(ii)

Show that if L(t)L(t) is hermitian and K(t)K(t) is unitary, then the matrix B(t)B(t) is skew-hermitian.

Consider the linear operators LL and MM defined by

Lψ(x,t,λ)=(i∂∂x+U(x,t,λ))ψ(x,t,λ) L\psi (x,t,\lambda ) = \left( i\frac{\partial }{\partial x} + U(x,t,\lambda ) \right) \psi (x,t,\lambda ) Mψ(x,t,λ)=(i∂∂t+V(x,t,λ))ψ(x,t,λ). M\psi (x,t,\lambda ) = \left( i\frac{\partial }{\partial t} + V(x,t,\lambda ) \right) \psi (x,t,\lambda ) .

Find the condition on LL and MM such that [L,M]=0[L,M] = 0, where [ , ][\, , \, ] denotes the commutator. The potentials U(x,t,λ)U(x,t,\lambda ) and V(x,t,λ)V(x,t,\lambda ) (square matrices) are typically chosen as elements of some semisimple Lie algebra.

?

Consider the differential operators

L:=d2dx2+f(x),B:=d3dx3+(32f+12x)ddx+34dfdx. L := \frac{d^2}{dx^2} + f(x), \quad B := \frac{d^3}{dx^3} + \left( \frac{3}{2}f + \frac{1}{2}x \right)\frac{d}{dx} + \frac{3}{4}\frac{df}{dx} .

From the condition [L,B]ϕ(x)=Lϕ(x)[L,B]\phi (x) = L\phi (x) find the ordinary differential equation that the function ff obeys.

?

Consider the linear equations

(∂∂u−U)ϕ(u,v,λ)=0,(∂∂v−V)ϕ(u,v,λ)=0 \left( \frac{\partial }{\partial u} - U \right)\phi (u,v,\lambda ) = 0, \quad \left( \frac{\partial }{\partial v} - V \right)\phi (u,v,\lambda ) = 0

where U,VU, V are functions of u,v,λu, v, \lambda. The linear operators

L=∂∂u−U,M=∂∂v−V L = \frac{\partial }{\partial u} - U, \quad M = \frac{\partial }{\partial v} - V

are the Lax pair. The integrability condition for ϕ\phi is the zero curvature condition on the Lax pair

[L,M]ϕ=0. [L,M]\phi = 0 .

Let g(u,v)g(u,v) be an invertible 2×22 \times 2 matrix and let

U=g−1∂g∂u,V=g−1∂g∂v. U = g^{-1}\frac{\partial g}{\partial u}, \quad V = g^{-1}\frac{\partial g}{\partial v} .

Find the nonlinear differential equation for [L,M]=0[L,M] = 0.

?

Consider the n×nn \times n tridiagonal matrices

A=(a1b10…0b1a2b2…0⋮⋮⋮bn−10…0bn−1an),B=(0b10…0−b10b2…0⋮⋮⋮bn−10…0−bn−10) A = \begin{pmatrix} a_1 & b_1 & 0 & \ldots & 0 \\ b_1 & a_2 & b_2 & \ldots & 0 \\ \vdots & \vdots & \vdots & & b_{n-1} \\ 0 & \ldots & 0 & b_{n-1} & a_n \end{pmatrix}, \quad B = \begin{pmatrix} 0 & b_1 & 0 & \ldots & 0 \\ -b_1 & 0 & b_2 & \ldots & 0 \\ \vdots & \vdots & \vdots & & b_{n-1} \\ 0 & \ldots & 0 & -b_{n-1} & 0 \end{pmatrix}

where the entries aja_j and bjb_j are real analytic functions of tt. Thus AA is symmetric over R\mathbb {R} and BB is skew-symmetric over R\mathbb {R}.

?
(i)

Let n=2n=2. Find the equations of motion (Lax representation)

dAdt=[A,B](t). \frac{dA}{dt} = [A,B](t) .

Find a first integral.

(ii)

Let n=3n=3. Find the equations of motion

dAdt=[A,B](t). \frac{dA}{dt} = [A,B](t) .

Find a first integral.

?
(i)

Show that the nonlinear autonomous system of first order differential equations

du1dt=u2,du2dt=u1u3,du3dt=−u1u2 \frac{du_1}{dt} = u_2, \quad \frac{du_2}{dt} = u_1u_3, \quad \frac{du_3}{dt} = -u_1u_2

admits the Lax representation dL/dt=[B,L](t)dL/dt = [B,L](t) with

L=(u3(u2−12u12)/2(u2+12u12)/212u12),B=(u1−u1/2u1/20). L = \begin{pmatrix} u_3 & (u_2-\frac{1}{2}u_1^2)/\sqrt{2} \\ (u_2+\frac{1}{2}u_1^2)/\sqrt{2} & \frac{1}{2}u_1^2 \end{pmatrix}, \quad B = \begin{pmatrix} u_1 & -u_1/\sqrt{2} \\ u_1/\sqrt{2} & 0 \end{pmatrix} .
(ii)

Find the first integrals from \tr[L]\tr [L] and 12\tr[L2]\frac{1}{2}\tr [L^2].

Many soliton equations in 1+11+1 dimensions can be written as an integrability condition for a set of linear equations called the linear scattering problem

dv→=Ωv→ d\overrightarrow {v} = \Omega \overrightarrow {v}

where Ω\Omega is an N×NN \times N matrix consisting of a one-parameter ζ\zeta (spectral parameter) family of differential one-forms. The integrability condition is Θ=0\Theta = 0 (zero curvature equation), where

Θ:=dΩ−Ω∧Ω \Theta := d\Omega - \Omega \wedge \Omega

dd is the exterior derivative and ∧\wedge is the exterior product. Assume that

Ω=(ζD+P(x,t))dx+Q(ζ,x,t)dt,\tr[Ω]=0 \Omega = (\zeta D + P(x,t))dx + Q(\zeta ,x,t)dt, \quad \tr [\Omega ] = 0

where ζ\zeta is the complex spectral parameter, DD is a diagonal N×NN \times N matrix with djjd_{jj} complex and pairwise distinct. The N×NN \times N matrix P(x,t)P(x,t) is off-diagonal and it is assumed to satisfy the boundary conditions

Pjk(x,t),∂Pjk(x,t)∂x,…,(∂∂x)mPjk(x,t),⋯→0as∣x∣→∞. P_{jk}(x,t), \quad \frac{\partial P_{jk}(x,t)}{\partial x}, \quad \ldots , \quad \left( \frac{\partial }{\partial x} \right)^m P_{jk}(x,t), \quad \cdots \to 0 \quad \text{as} \quad \left|x\right| \to \infty .
?
(i)

Write down the linear problem dv→=Ωv→d\overrightarrow {v} = \Omega \overrightarrow {v} for this case.

(ii)

Give the zero-curvature equation dΩ−Ω∧Ω=0d\Omega - \Omega \wedge \Omega = 0 for this case.

A number of Hamilton systems can be written in the form

dLdt=[A,L](t) \frac{dL}{dt} = [A,L](t)

where LL and AA are time-dependent n×nn \times n matrices. This is called the Lax representation of the Hamilton system. We find that

dLkdt=[A,Lk](t) \frac{dL^k}{dt} = [A,L^k](t)

and the that \tr[Lk]\tr [L^k] (k=1,2,…k=1,2,\ldots) are first integrals. Consider the Hamilton function (Toda lattice)

H(p→,q→)=12(p12+p22+p32)+exp⁡(q1−q2)+exp⁡(q2−q3)+exp⁡(q3−q1). H(\overrightarrow {p},\overrightarrow {q}) = \frac{1}{2}(p_1^2+p_2^2+p_3^2) + \exp (q_1-q_2) + \exp (q_2-q_3) + \exp (q_3-q_1) .

Introducing the quantities

aj:=12exp⁡(12(qj−qj+1)),bj:=12pj a_j := \frac{1}{2}\exp \left( \frac{1}{2}(q_j-q_{j+1}) \right), \quad b_j := \frac{1}{2}p_j

and cyclic boundary conditions (i.e., q4≡q1q_4 \equiv q_1) we find that the Hamilton equations of motion take the form (with b3=0b_3 = 0)

dajdt=aj(bj−bj+1),db1dt=−2a12,db2dt=2(a12−a22),db3dt=2a22 \frac{da_j}{dt} = a_j(b_j-b_{j+1}), \quad \frac{db_1}{dt} = -2a_1^2, \quad \frac{db_2}{dt} = 2(a_1^2-a_2^2), \quad \frac{db_3}{dt} = 2a_2^2

where j=1,2j=1,2. Introducing the matrices (Lax pair)

L:=(b1a10a1b2a20a2b3),A:=(0−a10a10−a20a20) L := \begin{pmatrix} b_1 & a_1 & 0 \\ a_1 & b_2 & a_2 \\ 0 & a_2 & b_3 \end{pmatrix}, \quad A := \begin{pmatrix} 0 & -a_1 & 0 \\ a_1 & 0 & -a_2 \\ 0 & a_2 & 0 \end{pmatrix}

the equations of motion can be written as Lax representation. From LL we find the first integral as \tr[Ln]\tr [L^n], where n=1,2,…n = 1,2,\ldots, where \tr\tr denotes the trace. We obtain

\trL=b1+b2+b3,\tr[L2]=b12+b22+b32+2a12+2a22. \tr L = b_1+b_2+b_3, \quad \tr [L^2] = b_1^2+b_2^2+b_3^2+2a_1^2+2a_2^2 .

Write a SymbolicC++ program lax.cpp that finds [L,A][L,A] and shows that \tr(L)\tr (L), \tr(L2)\tr (L^2) and the determinant of LL are first integrals.

?

Let ∂:=∂/∂x\partial := \partial /\partial x. The two pseudodifferential operators

L:=∂+∑j=1∞uj+1∂−j,W:=1+∑j=1∞wj∂−j L := \partial + \sum _{j=1}^{\infty } u_{j+1}\partial ^{-j}, \quad W := 1+\sum _{j=1}^{\infty } w_j\partial ^{-j}

are called the Lax operator and gauge operator, respectively. With ∂−1∂=∂∂−1=1\partial^{-1}\partial = \partial \partial^{-1} = 1 the generalized Leibniz rule is satisfied

∂kf⋅=∑j=0∞(kj)(∂jf)∂k−j. \partial ^k f\cdot = \sum _{j=0}^{\infty } \binom {k}{j} (\partial ^j f)\partial ^{k-j} .

for any k∈Zk \in \mathbb {Z} and L=W∂W−1L = W\partial W^{-1}. The KP hierarchy is given by the Lax equations

∂nL=[Bn,L] \partial _n L = [B_n,L]

where ∂n:=∂/∂tn\partial_n := \partial /\partial t_n and Bn:=L+nB_n := L_+^n is the differential part of Ln=L+n+L−nL^n = L_+^n+L_-^n with

L+n=∑j=0∞cjn∂j,L−n=∑j=−∞−1cjn∂j. L_+^n = \sum _{j=0}^{\infty } c_j^n \partial ^j, \quad L_-^n = \sum _{j=-\infty }^{-1} c_j^n \partial ^j .

Derive the first three equations of the KP hierarchy.

?

Let v→(x,t)=(v1(x,t),v2(x,t))T\overrightarrow {v}(x,t) = (v_1(x,t),v_2(x,t))^T. Consider the linear equations

∂v→∂x=(i∣u∣2−2iζ2i∂u/∂t+2ζui∂u∗/∂t−2ζ∗u−i∣u∣2+2iζ2)v→ \frac{\partial \overrightarrow {v}}{\partial x} = \begin{pmatrix} i\left|u\right|^2-2i\zeta ^2 & i\partial u/\partial t+2\zeta u \\ i\partial u^*/\partial t-2\zeta ^*u & -i\left|u\right|^2+2i\zeta ^2 \end{pmatrix} \overrightarrow {v} ∂v→∂t=(−iζu−u∗iζ)v→. \frac{\partial \overrightarrow {v}}{\partial t} = \begin{pmatrix} -i\zeta & u \\ -u^* & i\zeta \end{pmatrix} \overrightarrow {v} .

Show that the compatibility condition of these two equations with ∂ζ/∂t=0\partial \zeta /\partial t = 0 provides the unstable nonlinear Schrodinger equation

i∂u∂x+∂2u∂t2+2∣u∣2u=0. i\frac{\partial u}{\partial x} + \frac{\partial ^2 u}{\partial t^2} + 2\left|u\right|^2u = 0 .
?

References

[StHa17] Steeb W., Hardy Y. Problems and Solutions in Introductory and Advanced Matrix Calculus. 2nd ed., World Scientific Publishing, 2017.
[STH12] Steeb W., Tanski I., Hardy Y. Problems and Solutions for Groups, Lie Groups, Lie Algebras with Applications. World Scientific Publishing, 2012.