【SPONSORED LINK】

自由シュレディンガー方程式の基本解とユニタリ群

自由Schrödinger(シュレディンガー)方程式とは,線形の偏微分方程式

\begin{align*} i\partial_{t}u(t,x)+\Delta u(t,x)=0 \end{align*}

のことで,非線形Schrödinger方程式を考える際にも基本となる方程式です.

ここに

  • iは虚数単位
  • \partial_{t}=\frac{\partial}{\partial t}
  • \Delta=\sum_{i=1}^N\frac{\partial^2}{\partial x_i^2}

です.

自由Schrödinger方程式の初期値問題の解uは,[Stoneの定理]を用いてu(t,x)=e^{it\Delta}u_0(x)と表すことができ,このe^{it\Delta}自由Schrödinger発展作用素といいます.

[Stoneの定理]を用いるためには,自由Schrödinger発展作用素e^{it\Delta}が強連続ユニタリ群であることを述べる必要があり,そのためにこの記事では自由Schrödinger(シュレディンガー)方程式の解を考えます.

【SPONSORED LINK】

自由Schrödinger方程式の基本解

まずは自由Schrödinger方程式の解を直接計算しましょう.

基本解の定義

他の方程式にも適用できるより広い基本解の定義はありますが,Schrödinger方程式に限れば基本解は以下のように定義されます.

[基本解] 自由Schrödinger方程式の初期値問題

\begin{align*} \begin{cases} i\partial_{t}u(t,x)+\Delta u(t,x)=0 & (t,x)\in\R\times\R^d\\ u(0,x)=u_0(x) & x\in\R^d\\ \end{cases} \end{align*}

u,u_0\in\mathcal{S}(\R^d)の解uをSchrödinger方程式の基本解 (fundamental solution)という.

\mathcal{S}(\R^d)はSchwartz空間です.Schwartz空間について,詳しくは以下の記事を参照してください.

より正確には基本解は\mathcal{S}'(\R^d)における超関数解として定めます.

よって,以下の計算では単純な関数としてみると怪しいところが散見されますが,超関数として考えると全て正当化されます.

なお,この記事では非線形の場合まで考えることはしませんが,基本解は非線形項fが加わった非線形Schrödinger方程式i\partial_{t}u(t,x)+\Delta u(t,x)=fを考える際にも基本となります.

基本解の形

さて,Schrödinger方程式の基本解の形を求めましょう.

上のd次元の自由Schrödinger方程式について,u,u_0\in\mathcal{S}(\R^d)なら基本解u

\begin{align*} u(x) =&\mathcal{F}^{-1}\brc{e^{-it|\xi|^2}\widehat{u_0}(\xi)}(x) \\=&\frac{1}{(4\pi it)^{d/2}}\int_{\R^d}e^{-\frac{|x-y|^2}{4it}}u_0(y)\,dy \end{align*}

となる.

ただし,\mathcal{F}はFourier変換である:

\begin{align*} &\hat{f}(\xi) =\mathcal{F}[f](\xi) :=\frac{1}{(2\pi)^{d/2}}\int_{\R^d}e^{-ix\cdot\xi}f(x)\,dx, \\&\check{f}(x) =\mathcal{F}^{-1}[f](x) :=\frac{1}{(2\pi)^{d/2}}\int_{\R^d}e^{i\xi\cdot x}f(\xi)\,d\xi \end{align*}

i\partial_{t}u(t,x)+\Delta u(t,x)=0の両辺にxについてFourier変換を施すと,

\begin{align*} i\partial_{t}\hat{u}(t,\xi)-|\xi|^2\hat{u}(t,\xi)=0 \end{align*}

となる.両辺に-ie^{it|\xi|^2}をかけて

\begin{align*} \pd{}{t}\bra{e^{it|\xi|^2}\hat{u}(t,\xi)}=0 \end{align*}

となり,tについて両辺積分すると

\begin{align*} &\int_{0}^{t}\pd{}{\tau}\bra{e^{i\tau|\xi|^2}\hat{u}(\tau,\xi)}\,d\tau=0 \\\iff&e^{it|\xi|^2}\hat{u}(t,\xi)-e^{i0|\xi|^2}\hat{u}(0,\xi)=0 \\\iff&\hat{u}(t,\xi)=e^{-it|\xi|^2}\widehat{u_0}(\xi) \\\iff&u(t,x)=\mathcal{F}^{-1}\brc{e^{-it|\xi|^2}\widehat{u_0}}(x) \\\iff&u(t,x)=\frac{1}{(2\pi)^{d/2}}\bra{\mathcal{F}^{-1}\brc{e^{-it|\xi|^2}}*u_0}(x) \end{align*}

となる.ここに,f*gfgの合成積 (convolution)である.計算により

\begin{align*} \frac{1}{(2\pi)^{d/2}}\mathcal{F}^{-1}\brc{e^{-it|\xi|^2}}(x) =\frac{1}{(4\pi it)^{d/2}}e^{-\frac{|x|^2}{4it}} =:S_t(x) \end{align*}

となるから,解u

\begin{align*} u(x) =&(S_t*u_0)(x) \\=&\int_{\R^d}S_t(x-y)u_0(y)\,dy \\=&\frac{1}{(4\pi it)^{d/2}}\int_{\R^d}e^{-\frac{|x-y|^2}{4it}}u_0(y)\,dy \end{align*}

となる.

なお,ここで現れたS_t自由Schrödinger核といいます.

自由Schrödinger発展作用素

いま求めた初期値u_0からSchrödinger方程式の基本解を与える作用素を自由Schrödinger発展作用素といいます.

以降,自由Schrödinger発展作用素をT_tとします:

\begin{align*} T_{t}u_0:=S_t*u_0. \end{align*}

一般に,ノルム空間XとBanach空間Yに対して,有界線形作用素X\to YX上で稠密な定義域をもつとき,この作用素はX上の有界線形作用素に一意的に拡張できるのでした.

このことに注意すると,以下が成り立ちます.

任意のt\in Tに対して,自由Schrödinger発展作用素T_tは定義域をL^2(\R^d)とするL^2(\R^d)上の有界線形作用素に一意に拡張できる.

L^2(\R^d)ノルムを\|\cdot\|_2で表す.

[1] 合成積の線形性から,任意のf,g\in L^2(\R^d), \alpha,\beta\in\Cに対して

\begin{align*} {T_t}(\alpha f+\beta g) =&{S_t}*(\alpha f+\beta g) \\=&\alpha{S_t}*f+\beta{S_t}*g \\=&\alpha{T_t}f+\beta{T_t}g \end{align*}

となるので,T_tは線形である.

[2] 任意のf\in L^2(\R^d)に対して,[Plancherelの定理]を2回用いることにより

\begin{align*} \|{T_t}f\|_2 =&\|{S_t}*f\|_2 =\nor{\mathcal{F}^{-1}\brc{e^{-it|\xi|^2}\hat{f}}}_2 \\=&\nor{e^{-it|\xi|^2}\hat{f}}_2 =\nor{\hat{f}}_2 =\|f\|_2 \end{align*}

となって,T_tは有界である.

[1], [2]からT_tは有界線形作用素である.このことと

  • T_tはSchwartz空間\mathcal{S}(\R^d)で定義されていること
  • \mathcal{S}(\R^d)\|\cdot\|_2ノルムによる閉包がL^2(\R^d)であること

を併せると,関数解析の一般論からT_tは定義域をL^2(\R^d)とするL^2(\R^d)上の有界線形作用素に一意に拡張できる.

この証明から分かるように,T_tL^2有界よりも強くL^2等長となっている.

強連続ユニタリ群

自由Schrödinger発展作用素の族\{T_t\}_{t\in\R}が強連続ユニタリ群であることを示しましょう.

強連続ユニタリ群

一般に,Hilbert空間上の有界線形作用素の族がユニタリ群であるとは,次のように定義されます.

[強連続ユニタリ群] Hilbert空間\mathcal{H}上の有界線形作用素の族\{T_t\}_{t\in\R}\mathcal{H}ユニタリ群であると

  1. 任意のs,t\in\Rに対して,{T_s}{T_t}=T_{s+t}である.
  2. T_0=Iである.ただし,I\mathcal{H}上の恒等作用素である.
  3. 任意のf\in\mathcal{H}, t_0\in\Rに対して,\lim\limits_{h\to 0}\|T_{t_0+h}f-T_{t_0}f\|_{\mathcal{H}}=0である.
  4. 任意のt\in\Rに対して,T_{t}はユニタリ作用素(全単射かつ等長)である.

を満たすことをいう.

なお,条件1〜3を満たす\{T_t\}_{t\in\R}強連続群C_0群などといいます.

自由Schrödinger発展作用素

自由Schrödinger発展作用素の族\{T_t\}_{t\in\R}L^2(\R^N)上強連続ユニタリ群である.

任意にs,t\in\R, f\in L^2(\R^N)をとる.

\begin{align*} {T_s}{T_t}f =&\mathcal{F}_{\xi}^{-1}\brc{e^{is|\xi|^2}\mathcal{F}_{x}\brc{\mathcal{F}_{\xi}^{-1}\brc{e^{it|\xi|^2}\hat{f}}}} \\=&\mathcal{F}_{\xi}^{-1}\brc{e^{is|\xi|^2}\bra{e^{it|\xi|^2}\hat{f}}} \\=&\mathcal{F}_{\xi}^{-1}\brc{e^{i(s+t)|\xi|^2}\hat{f}} =T_{s+t}f \end{align*}

だから,{T_s}{T_t}=T_{s+t}が成り立つ.

Fourier変換がL^2(\R^N)のユニタリ作用素であることに注意すると,

\begin{align*} {T_0}f =\mathcal{F}^{-1}\brc{e^{i0|\xi|^2}\hat{f}} =\mathcal{F}^{-1}\brc{\hat{f}} =f \end{align*}

となって,T_0=Iが成り立つ.

Fourier変換の線形性,L^2等長性と,ノルムの連続性から

\begin{align*} &\lim_{h\to0}\|T_{t+h}f-{T_t}f\|_2 \\=&\lim_{h\to0}\nor{\mathcal{F}^{-1}\brc{e^{i(t+h)|\xi|^2}\hat{f}}-\mathcal{F}^{-1}\brc{e^{it|\xi|^2}\hat{f}}}_2 \\=&\lim_{h\to0}\nor{\mathcal{F}^{-1}\brc{\bra{e^{i(t+h)|\xi|^2}-e^{it|\xi|^2}}\hat{f}}}_2 \\=&\lim_{h\to0}\nor{e^{it|\xi|^2}\bra{e^{ih|\xi|^2}-1}\hat{f}}_2 \\=&\nor{\lim_{h\to0}e^{it|\xi|^2}\bra{e^{ih|\xi|^2}-1}\hat{f}}_2 =\|0\|_{2} =0 \end{align*}

となる.

先ほど示したようにT_tL^2(\R^N)全体で定義される等長作用素である.また,任意のg\in L^2(\R^d)に対して

\begin{align*} g=T_0g =T_{t}(T_{-t}g) \end{align*}

でありT_{-t}gだからT_tは全射で,f\in L^2(\R^d){T_t}f=0を満たせば,左からT_{-t}を作用させてf=0となるからT_tは単射である.

Stoneの定理の適用

ここで,[Stoneの定理]を確認しておきましょう.

[Stoneの定理] \{T_t\}_{t\in\R}をHilbert空間\mathcal{H}全体で定義された\mathcal{H}上の有界線形作用素の族とする.このとき,\{T_t\}_{t\in\R}がユニタリ群であるためには,\mathcal{H}上の自己共役作用素Aが存在して,T_t=e^{itA}を満たすことが必要十分である.

[Stoneの定理]について,詳しくは以下の記事を参照してください.

[Stoneの定理]より,自由Schrödinger発展作用素の族\{T_t\}のユニタリ群であることから,L^2(\R^d)上の自己共役作用素Aが存在してT_t=e^{itA}と表されます.

このとき,任意のf\in L^2(\R^d)に対して

\begin{align*} Af =-i\od{}{t}(e^{itA}f)|_{t=0} =-i\od{}{t}(T_tf)|_{t=0} \end{align*}

となります.T_tfは自由Schrödinger方程式の解なので,

\begin{align*} i\partial_{t}(T_tf)+\Delta (T_tf)=0 \iff& -i\partial_{t}(T_tf)=\Delta (T_tf) \\\Ra\ & -i\partial_{t}(T_tf)|_{t=0}=\Delta (T_0f) \\\iff& -i\partial_{t}(T_tf)|_{t=0}=\Delta f \end{align*}

なので,Af=\Delta fが得られ,A=\Deltaとなりますね.

このことから,ラプラシアン\DeltaL^2(\R^d)の自己共役作用素であり,自由Schrödinger発展作用素はe^{it\Delta}とスッキリ表せることが分かりました.

以上をまとめると,以下のようになります.

自由Schrödinger方程式の初期値問題

\begin{align*} \begin{cases} i\partial_{t}u(t,x)+\Delta u(t,x)=0 & (t,x)\in\R\times\R^d\\ u(0,x)=u_0(x) & x\in\R^d\\ \end{cases} \end{align*}

u,u_0\in\mathcal{S}(\R^d)の解は

\begin{align*} e^{it\Delta}u_0(x) =&\mathcal{F}\brc{e^{-it|\xi|^2}\hat{u_0}(\xi)}(x) \\=&\frac{1}{(4\pi it)^{N/2}}\int_{\R^d}e^{-\frac{|x-y|^2}{4it}}u_0(y)\,dy \end{align*}

と表せる.

次の記事では,Schrödinger方程式の基本解の不等式として基本的な[L^pL^q評価]を考え方から説明します.

参考文献

「非線形発展方程式の実解析的方法」(小川卓克 著,丸善出版(シュプリンガー現代数学シリーズ))

非線形発展方程式について議論するには,Lebesgue空間L^p(\R^d),Sobolev空間W^{k,p}(\R^d), H^{s}_p(\R^d),Besov空間B^{s}_{p,\sigma}といった解を考えるための種々の関数空間を理解することが重要です.

本書は関数空間に関する予備知識をじっくりと準備し,

  • 波動方程式
  • 熱方程式
  • Schrödinger方程式
  • Navier-Stokes方程式

といった非線形発展方程式を考えていきます.

本書の特徴は,様々な非線形発展方程式を広く扱っている点と,証明へのアプローチを説明して直感的な理解を促しているです.

本書が全19章と多くの章から構成されていることからも,広くトピックを扱っていることが見てとれますね.

誤植が多いのが1つ残念な点ではありますが,これほどに広く丁寧に非線形発展方程式を扱っている和書は他に見当たらず,この分野の基礎や考え方をカバーするには良い教科書と言えます.

シェアする

  • このエントリーをはてなブックマークに追加

フォローする

最後までありがとうございました!

以下の関連記事もいかがですか?

SPONSORED LINK
関連記事

記事一覧はこちらからどうぞ!

記事一覧は

こちら

Twitterを

フォロー

大学院入試

解答例

大学受験

解説ブログ