ラックス・ペア

ラックス・ペア (英: Lax pair)は、数学の可積分系の理論における用語であり、ある微分方程式 (時刻発展型偏微分方程式) がラックス方程式を用いて書き換え可能な場合に、その中で使われる時刻に依存する作用素の対を指す。このような場合、元の微分方程式は、そのラックス・ペアを持つと表現される。これらはアメリカ合衆国の数学者であるピーター・ラックスによって連続媒体中の解を論ずるために導入された。時刻発展型偏微分方程式をラックス方程式で書き換えることにより、逆散乱法（English版）を用いて微分方程式の解を求めることができるようになる。この方法により、従来の方法では解くことができなかった多数の非線形時刻発展型偏微分方程式の厳密解が得られており、それらは、ほとんどの場合ソリトン解を持つことが知られている。

定義

1組のラックス・ペア L(t) 、A(t) は、あるヒルベルト空間 [math]\mathbb{H}[/math]に作用する、時刻 t に依存する作用素で、次のラックス方程式を満たすものである。

[math]\frac{\partial L(t)}{\partial t} = [A(t),L(t)]\ [/math]

ここで [math][A,L] = AL - LA[/math] は交換子である。

ラックス方程式を微分方程式の解法に用いる場合、[math]\mathbb{H}[/math] は通常は関数空間である。以後その空間方向の独立変数を x で表すことにする。 x は1次元の場合も、もっと多次元の場合も有り得る。[math]\mathbb{I}[/math] を時刻t の定義域とする場合、[math]\mathbb{I}[/math] から [math]\mathbb{H}[/math] への写像を考えて、それを [math]\psi(t)[/math] などと表すことにする。個々のt ごとに[math]\psi(t)[/math] は[math]\mathbb{H}[/math] の元であり、変数x の関数でもある。このため [math]\psi[/math] をt、x の関数と考えて [math]\psi(t,x )[/math] と表すことも有り得ることにする。

L(t) 、A(t) も、実際には x の関数でもあるので、同様に L(t,x) 、A(t,x) と表すことも、あるいは t も x も省略して L 、A と表すことも有り得ることにする。ラックス方程式の左辺で偏微分記号を用いているのは、時刻による微分であることをはっきりさせるためである。

また、 L と A は決定すべき未知関数 u(t,x) を内部に含み、u がある微分方程式を満たすことが、ラックス方程式が成立する条件となる。この場合、u についてのその微分方程式がラックス方程式で書き換えられたと呼ぶのである。

左辺の時刻微分について注意しておく。これは作用素の時刻微分であり、次のように定義される。

[math]\frac{\partial L(t)}{\partial t} = \lim_{h\to 0} \frac{L(t+h)-L(t)}{h}[/math]

u が x のみの関数であれば、時刻と共に変動しないので、当然その時刻微分は 0 である。これと同様に、[math]\frac{\partial}{\partial x}[/math] あるいは [math] \frac{\partial^n}{\partial x^n}[/math] は時刻と共に変動しない作用素であり、その時刻微分は 0 である。

L(t) 、A(t) は、[math]\mathbb{H}[/math] のスカラー体の上の任意の時刻関数 (これを[math]\lambda(t)[/math] とする) と交換する。つまり、

[math][A(t),\lambda(t)] = [L(t),\lambda(t)] = 0\ [/math]

である。これは L(t) 、A(t) が時刻微分作用素 [math]\frac{\partial}{\partial t}[/math] をその内部に含んでいないことを意味する (これが、L(t) 、A(t) は [math]\mathbb{H}[/math] に作用するという表現の暗黙的な意味である。 [math]\frac{\partial}{\partial t}[/math] は [math]\mathbb{H}[/math] への作用素ではなく、 [math]\mathbb{I}\times\mathbb{H}[/math] への作用素である)。

ラックス方程式の形式は不変のままで L 、A の形式を変えることにより、未知関数 u が満たすべき微分方程式を様々な形式に変化させることができる。下記の例のように、ラックス方程式から導かれる微分方程式は大抵の場合、非線形偏微分方程式となる。

例

次の微分方程式をKdV方程式(Korteweg–de Vries equation)と呼ぶ。この方程式は、その数値解析において初めてソリトン解が発見されたことで有名である。KdV方程式を例に、時刻発展型偏微分方程式がラックス方程式で書き換え可能であるとはどういう意味かをもう少し具体的に説明する。

[math]\frac {\partial u}{\partial t} = 6u\frac {\partial u}{\partial x} - \frac {\partial^3 u}{\partial x^3}[/math]

ここで、解 u は時刻 t および空間方向について1次元の独立変数x の関数である (つまり[math]u = u (t,x )[/math] )。

実は、ラックス・ペアを次のように取ると、ラックス方程式からKdV方程式が導かれる。

[math]L = -\partial_x^2 + u[/math]

[math]A = -4\partial_x^3 + 3u\partial_x + 3\partial_x u = -4\partial_x^3 + 6u\partial_x + 3u_x [/math]

ここで、作用素 [math]\partial_t[/math] および [math]\partial_x[/math] はそれぞれ、[math]\frac{\partial}{\partial t}[/math]、[math]\frac{\partial}{\partial x}[/math] を表すものとする。

この場合のL はスツルム＝リウヴィル型作用素と呼ばれる形式になっているが、これは量子力学におけるシュレーディンガー方程式のハミルトニアンと同じ形をしていることを注意しておく (実はこの類似性から、量子力学で用いられていた逆散乱法の応用による解法が可能となったのである)。

また、記号上の注意であるが、上の式では、L、A を構成する各作用素は、「作用対象に右側のものから順に作用する」という規則に従う。例えば、上の式の中にある u は、作用対象である[math]\mathbb{H}[/math] の任意の元を [math]\psi[/math] とすると、作用の結果 [math]u\psi[/math] を生成する作用素である。ここまでは当然のように思えるが、[math]\partial_x u[/math] という作用素を考えると、直観とはやや異なった結果になる。この規則の下に[math]\partial_x u[/math] を [math]\psi[/math] に実際に作用させてみると、結果は次のようになる。

[math]\partial_x u \psi = \partial_x (u \psi) = (\partial_x u)\psi + u(\partial_x \psi)[/math]

[math]\psi[/math] を用いずに作用素間のみの関係として見れば、

[math]\partial_x u = u\partial_x + u_x[/math]

である。以降、上式で用いたように、u の x または t についての偏微分(偏導関数)を表すには、u に [math]u_t, u_x, u_{xx}\ [/math] などのように偏微分した変数を添えることにする。

上式の表現によれば [math]\partial_x^2 u[/math]、[math]\partial_x^3 u[/math] については次のようになる。

[math]\partial_x^2 u = u\partial_x^2 + 2u_x\partial_x + u_{xx}[/math]

[math]\partial_x^3 u = u\partial_x^3 + 3u_x\partial_x^2 + 3u_{xx}\partial_x + u_{xxx}[/math]

以上の関係を交換子を用いて表すと次のようになる。

[math][\partial_x, u] = u_x[/math]

[math][\partial_x^2, u] = 2u_x\partial_x + u_{xx}[/math]

[math][\partial_x^3, u] = 3u_x\partial_x^2 + 3u_{xx}\partial_x + u_{xxx}[/math]

以上で準備ができたので、実際にラックス方程式 [math]\frac{\partial L}{\partial t} = [A,L][/math] からKdV方程式を導いてみる。ラックス方程式の左辺が [math]u_t\ [/math] となることはすぐ分かる。 [A,L] の計算は以下の通りである。

[math][A,L] = [-4\partial_x^3 + 6u\partial_x + 3u_x, -\partial_x^2 + u] \ [/math]

[math]= [-4\partial_x^3, u] + [6u\partial_x, -\partial_x^2] + [6u\partial_x, u] + [3u_x, -\partial_x^2]\ [/math]

[math]= -4(3u_x\partial_x^2 + 3u_{xx}\partial_x + u_{xxx}) + 6(2u_x\partial_x + u_{xx})\partial_x + 6uu_x + 3(2u_{xx}\partial_x + u_{xxx}) \ [/math]

[math]= -12u_x\partial_x^2 - 12u_{xx}\partial_x -4u_{xxx} + 12u_x\partial_x^2 + 6u_{xx}\partial_x + 6uu_x + 6u_{xx}\partial_x + 3u_{xxx}\ [/math]

[math]= 6uu_x - u_{xxx}\ [/math]

従って、

[math]u_t = 6uu_x - u_{xxx}\ [/math]

となって、確かに KdV 方程式と一致する。上のラックス方程式においては [math]\partial_x[/math]、 [math]\partial_t[/math] を最も右側に持つ項はお互いに打ち消しあって最終的には式に現れない。このような場合 L と A は準可換 (semicommutative) であると呼ぶ。

等スペクトル性

[math]\lambda(t)[/math]、[math]\psi(t)[/math] をそれぞれ時刻 t における L(t) の1つの固有値、およびその固有値を持つ固有ベクトルの1つとすれば [math]L(t) \psi(t) = \lambda(t) \psi(t)[/math] である。少々天下り的ではあるが、この [math]\lambda(t)[/math] が時刻と共に変動しない条件を考えてみる。まずこの式の両辺を時刻 t で偏微分すれば、

[math] \frac{\partial L}{\partial t} \psi + L \frac{\partial \psi}{\partial t} = \frac{\partial \lambda}{\partial t} \psi + \lambda \frac{\partial \psi}{\partial t}[/math]

である。[math] \frac{\partial L}{\partial t}[/math] をラックス方程式を用いて書き換え、さらに [math] \frac{\partial \lambda}{\partial t} = 0[/math] とすれば、

[math] (AL-LA) \psi + L \frac{\partial \psi}{\partial t} - \lambda \frac{\partial \psi}{\partial t} = 0[/math]

λ(t) と A(t) が交換可能なことを利用して、さらに整理すれば、

[math] (L- \lambda) (A - \frac{\partial }{\partial t}) \psi = 0[/math]

従って、

[math]\frac{\partial \psi}{\partial t}=A \psi[/math]

であれば、[math]\lambda(t)[/math] は時刻と共に変動しないことが分かる (ただしこれは十分条件であって必要条件ではないことを注意しておく)。このような場合、行列 (または作用素) である L(t) は t の変動に関して等スペクトル的（English版）であると表現される。

一方、ラックス方程式は、実は量子力学におけるハイゼンベルクの運動方程式の特別な場合(観測可能量を表す作用素が時刻t を陽に含まない場合) と全く同じ形式をしている。ハイゼンベルク方程式においては、系のハミルトニアンを H(t) とすると、A(t) に相当するのは、[math]- \frac{i}{\hbar}H(t)[/math] である。また位置や運動量についてのハイゼンベルク形式の観測可能量(オブザーバブル)が L(t) に相当する。ハイゼンベルク方程式との類推から、ラックス方程式の解は、

[math]L(t)=U(t,t_0) L(t_0) U(t,t_0)^{-1} \ [/math]

と表される。ここで t、[math]t_0[/math] は任意の時刻を表し、[math]U(t,t_0)[/math] は次の方程式の解であり、時間推進作用素(time evolution operator)と呼ばれる。

[math] \frac{\partial}{\partial t} U(t,t_0) = A(t) U(t,t_0), \qquad U(t_0,t_0) = I[/math]

I は単位行列を表す。このような[math]U(t,t_0)[/math] は常に存在することが証明でき、特に任意の時刻 [math]t_1[/math]、[math]t_2[/math] で、[math]A(t_1)[/math] と [math]A(t_2)[/math] が可換であれば、

[math] U(t,t_0) = e^{\int_{t_0}^{t} A(\tau) d\tau}[/math]

である。また [math]t_1[/math] を任意の時刻として、

[math] U(t,t_1) U(t_1,t_0) = U(t,t_0)\ [/math]

という関係が常に成り立つ。従って、

[math]U(t,t_0)^{-1} = U(t_0,t) \ [/math]

である。なお、もし A(t) が歪エルミートであれば、[math]U(t,t_0)[/math] はユニタリとなることを注意しておく。

さて、 [math]U(t,t_0)[/math] を用いれば、[math]\psi(t)[/math] は、

[math]\psi(t) = U(t,t_0) \psi(t_0) \ [/math]

と表現される。実際これに L(t) を作用させれば、

[math]L(t) \psi(t) = U(t,t_0) L(t_0) U(t,t_0)^{-1} U(t,t_0) \psi(t_0) \ [/math]

[math]= U(t,t_0) L(t_0) \psi(t_0) = U(t,t_0) \lambda(t_0) \psi(t_0) \ [/math]

[math]= \lambda(t_0) \psi(t) \ [/math]

となって、確かに L(t) の固有値 [math]\lambda[/math] は時刻と共に変動しないことが分かる。

以上をまとめると、あるラックス・ペアにおいて、時刻 [math]t_0[/math] で [math]L(t_0)[/math] 、[math]A(t_0)[/math] についての初期条件が与えられ、その時刻における固有値問題 [math]L(t_0) \psi(t_0) = \lambda(t_0) \psi(t_0)[/math] の解が得られれば、任意の時刻 t における固有値問題 [math]L(t) \psi(t) = \lambda(t) \psi(t)[/math] の解は、次の式で与えられるということである。

[math]\lambda(t)=\lambda(t_0) \ [/math] (固有値またはスペクトルは不変)

[math]\psi(t) = U(t,t_0) \psi(t_0) \ [/math]

逆散乱法とのリンク

上記の性質は逆散乱法のための基礎となる。この方法においては、時刻 [math]t_0[/math] で [math]u (t_0,x)[/math] は初期条件として与えられており、 [math]u(t_0,x)[/math] は [math]|x| \to \infty[/math] で [math]|u(t_0,x)| \to 0 [/math] を満たすものと仮定する (以降、[math]|u(t_0,x)|[/math] が十分小さい x についての領域を散乱領域と呼ぶことにする)。

この方法は、次のような概略にて進められる：

[math]L(t_0)[/math] のスペクトルを計算し、 [math]\lambda[/math] と [math]\psi(t_0,x)[/math] を得る。
散乱領域においては[math]A[/math] は既知と見なせるので、初期条件 [math]\psi(t_0,x)[/math] の下で、 [math]\psi[/math] を [math]\psi(t) = U(t,t_0) \psi(t_0)[/math] を用いて時間発展させる。
散乱領域における [math]\psi(t)[/math] が分かったので、これから [math]u(t,x)[/math] を逆散乱法で計算する。