1 ユニタリ群と弱解
定義 1.1 (Schrödinger 方程式の弱解).Hを複素 Hilbert 空間、AをH上の自己共役作用素、I⊆Rを0を含む区間とする。内積(u,v)は第一変数について線形とする。u0∈H、f∈Lloc1(I;H)に対して、u∈C(I;H)が初期値u0をもつ方程式iut=Au+fの 弱解 (weak solution) であるとは、u(0)=u0であり、任意のv∈D(A)について(u(⋅),v)が局所絶対連続であって
idtd(u(t),v)=(u(t),Av)+(f(t),v)がほとんどすべてのt∈Iで成り立つことをいう。外力の空間は Bochner 可積分性によって定める。
定理 1.2.Hを複素 Hilbert 空間、AをH上の自己共役作用素とし、u0∈Hとする。このとき
u(t)=U(t)u0,U(t)=e−itA(t∈R)はiut=Auの唯一のC(R;H)弱解であり、∥u(t)∥=∥u0∥を満たす。さらにu0∈D(A)ならば、グラフノルムに関して
u∈C(R;D(A))∩C1(R;H),iu′(t)=Au(t)が成り立つ。
証明.§E12.25 命題 3.2を−Aに適用すると、Uは生成作用素−iAをもつ強連続ユニタリ群である。v∈D(A)に対して
(U(t)u0,v)=(u0,U(−t)v)を微分すると、第二変数の共役線形性から
dtd(U(t)u0,v)=−i(U(t)u0,Av)を得る。したがってU(t)u0は弱解であり、ユニタリ性からノルムを保存する。
EをAのスペクトル測度、PN=E([−N,N])、N∈N≥1とする。初期値零の弱解wに対して、任意のh∈HについてPNh∈D(A)であるから
idtd(PNw(t),h)=(APNw(t),h).右辺はtの連続関数であり、APNは有界である。スカラーの積分等式をすべてのhについて用いると
PNw(t)=−i∫0tAPNw(s)dsとなる。ゆえにPNwはHのノルムに関して連続微分可能である。有界自己共役作用素A∣PNHの指数関数を掛けて微分すると
dtd(eitA∣PNHPNw(t))=0.初期値が零なのでPNw(t)=0である。PN→1の強収束からw(t)=0を得る。
u0∈D(A)の場合、生成作用素の定義と群の法則からu′(t)=−iU(t)Au0となる。§E12.15 定理 1.2の積の公式によりAU(t)u0=U(t)Au0であり、uとAuの強連続性からグラフノルムについての連続性が従う。▨
系 1.3.Hを複素 Hilbert 空間、Aを自己共役作用素とし、(ej)j∈N≥1をAej=λjejを満たすHの完全正規直交系とする。このとき、任意のu0∈Hについて
e−itAu0=j≥1∑e−itλj(u0,ej)ejがHで成り立ち、級数はt∈Rについて一様に収束する。
証明. スペクトル汎関数計算からU(t)ej=e−itλjejである。初期値の正規直交展開に有界作用素U(t)を作用させると表示を得る。第N部分和との差のノルムの二乗は∑j>N∣(u0,ej)∣2であり、tによらず零へ収束する。▨
例 1.4. 区間(0,π)の Dirichlet Laplace 作用素AD=−d2/dx2と初期値u0(x)=sinx+sin2xに対して、解は
u(t,x)=e−itsinx+e−4itsin2xである。実際、各項を微分するとiut=ADu、u(t,0)=u(t,π)=0、u(0,x)=u0(x)を得る。定理 1.2の一意性から、この関数はユニタリ群による解に等しい。二つの正弦関数は直交し、それぞれのノルムの二乗がπ/2なので、∥u(t)∥22=πである。
2 定義域とエネルギーの保存
命題 2.1.Aを複素 Hilbert 空間H上の自己共役作用素、m≥1を整数とする。このとき、すべてのt∈Rについて
e−itAD(Am)=D(Am),Ame−itAv=e−itAAmv(v∈D(Am))が成り立つ。v∈D(Am)に対してt↦e−itAvはグラフノルムに関して連続であり、∥Ame−itAv∥=∥Amv∥である。したがってグラフノルム(∥v∥2+∥Amv∥2)1/2も保存される。
証明.EをAのスペクトル測度とする。汎関数計算の積の公式§E12.15 定理 1.2と∣e−itλ∣=1により
∫R∣λ∣2md(E(λ)e−itAv,e−itAv)=∫R∣λ∣2md(E(λ)v,v).よって定義域への所属とAmのノルムは保存される。同じ積の公式から可換性を得る。時刻−tの作用で逆包含を得る。e−itAvとe−itAAmvはともにHで連続なので、グラフノルムに関する連続性も従う。▨
証明. スペクトル積分における定義域の条件は∫[0,∞)λd(E(λ)v,v)<∞である。絶対値一の乗数e−itλはこの積分を変えない。時刻−tによって全射性も従う。形式との一致は§E15.24 命題 1.1、全空間の勾配ノルムとの一致は§E15.24 命題 1.4による。▨
系 2.3.Ω⊆Rnを有界C∞開集合、ADを Dirichlet Laplace 作用素、m≥1を整数とする。このとき、
u0∈{v∈H2m(Ω):TΔjv=0 (0≤j<m)}ならば、u(t)=e−itADu0はすべてのt∈Rで同じ空間に属し、supt∈R∥u(t)∥H2m≤Cm∥u0∥H2mを満たす。定数CmはΩ,mのみに依存する。
証明.§E15.23 系 6.1によって表示された空間はD(ADm)に等しく、そのH2mノルムはグラフノルムと同値である。命題 2.1を適用すると所属と評価を得る。▨
命題 2.4.n≥1、s∈Rとする。Fourier 変換をv(ξ)=∫e−2πix⋅ξv(x)dxとし、超関数は共役を入れない複素線形双対とする。このとき、Hs(Rn)上で
eitΔv=F−1(e−4π2it∣ξ∣2v(ξ))と定めると、eitΔは強連続等長群である。s=0の場合はL2上のe−itA、A=−Δに一致する。
証明.§E15.13 定理 1.2の等長同型Jsにより、表示された作用素はL2上の絶対値一の乗数と共役である。群の法則と等長性は乗数の積と絶対値から従う。t→t0のとき、乗数の差の絶対値は2以下で各点零へ収束するので、優収束定理から強連続性を得る。§E15.24 命題 1.4によるAの Fourier 表示と汎関数計算からs=0の一致が従う。▨
例 2.5.(0,π)でej(x)=(2/π)1/2sin(jx)とし、u0=∑j≥1j−1ejとする。このとき、直交性と∑jj−2<∞により、この級数はL2で収束する。EをADのスペクトル測度とすると、ADej=j2ejなので
∫[0,∞)λd(E(λ)u0,u0)=j≥1∑j2j−2=∞.したがってu0∈/D(AD1/2)=H01(0,π)である。時刻tの解の係数はe−itj2/jであるから、同じ発散がすべての時刻で成立する。一般にも、ある時刻tでe−itADu0∈D(AD1/2)ならば、命題 2.2を時刻−tに適用してu0∈D(AD1/2)を得る。よってこの時間発展はL2の初期値を正の時刻で一律にH01へ写さない。
3 非斉次問題
系 3.1.Aを複素 Hilbert 空間H上の自己共役作用素とし、T>0、u0∈H、f∈L1(0,T;H)とする。このとき
u(t)=e−itAu0−i∫0te−i(t−s)Af(s)ds(0≤t≤T)はiut=Au+fの唯一のC([0,T];H)弱解であり、
∥u(t)∥≤∥u0∥+∫0t∥f(s)∥dsを満たす。積分は Bochner 積分であり、fの代表元によらない。
証明.§E15.24 定理 2.1を Schrödinger の場合(β=i)に外力g=−ifとして適用する。u′+iAu=−ifはiu′=Au+fと同じ方程式であり、表示されたuは同定理の積分表示に等しい。したがってu∈C([0,T];H)であり、uはすべてのv∈D(A)とt∈[0,T]について
(u(t),v)=(u0,v)−i∫0t((u(s),Av)+(f(s),v))dsを満たす唯一の連続なH値関数である。
連続関数w:[0,T]→Hとv∈D(A)について、s↦(w(s),Av)+(f(s),v)は[0,T]上で可積分である。したがって、uをwに替えた上の等式がすべてのt∈[0,T]で成り立つことは、(w(0),v)=(u0,v)であり、(w(⋅),v)が絶対連続であって、ほとんどすべてのtでidtd(w(t),v)=(w(t),Av)+(f(t),v)が成り立つことと同値である。D(A)はHで稠密であるから、すべてのv∈D(A)について(w(0),v)=(u0,v)であることはw(0)=u0と同値である。ゆえにuは唯一のC([0,T];H)弱解である。ほとんど至る所で等しい二つの可積分関数の Bochner 積分は等しいので、表示はfの代表元によらない。
0<t≤Tとし、同定理を区間[0,t]上でデータ(u0,−if)と(0,0)に適用する。後者の軟解は零であるから、データの差に対する評価により∥u(t)∥≤∥u0∥+∫0t∥f(s)∥dsを得る。t=0では評価は自明である。▨
系 3.2.Aを複素 Hilbert 空間H上の自己共役作用素、T>0、m≥1を整数、u0∈D(Am)、f∈L1(0,T;D(Am))とし、外力の可積分性をグラフノルムについて課す。このとき、系 3.1の解はC([0,T];D(Am))に属し、
Amu(t)=e−itAAmu0−i∫0te−i(t−s)AAmf(s)ds,∥Amu(t)∥≤∥Amu0∥+∫0t∥Amf(s)∥dsを満たす。さらにu0∈D(A)であり、f∈C1([0,T];H)またはf∈C([0,T];D(A))ならば、解はC1([0,T];H)∩C([0,T];D(A))に属する古典解である。
証明.命題 2.1によりe−itAはグラフノルムをもつ Banach 空間D(Am)上の強連続等長群でもある。§E12.26 定理 2.1をこの空間に適用すると、積分表示はD(Am)で連続となる。Am:D(Am)→Hは有界なので Bochner 積分と交換し、表示と評価を得る。
古典解の主張は、§E15.24 系 2.2をβ=iと外力g=−ifに適用したものである。gはfと同じ正則性をもち、u′+iAu=−ifはiu′=Au+fと同じ方程式である。▨
例 3.3.Aを複素 Hilbert 空間上の自己共役作用素とし、Ae=λe、∥e∥=1とする。このとき、初期値零、一定外力f(t)=eに対する解は
u(t)=⎩⎨⎧λe−itλ−1e,−ite,λ=0,λ=0.実際、積分表示の被積分関数はe−i(t−s)λeなので、sについて積分すると表示を得る。λ=0の場合は∥u(t)∥=2∣sin(tλ/2)∣/∣λ∣、λ=0の場合は∥u(t)∥=tである。したがって外力があるときにはHのノルムは一般に保存されない。
4 自由方程式の核と分散評価
定理 4.1.n≥1、t∈R∖{0}とする。このとき、Rn上の自由 Schrödinger 群eitΔの Fourier 乗数はe−4π2it∣ξ∣2である。任意のu0∈S(Rn)に対して、解は各点で
eitΔu0(x)=∫RnKt(x−y)u0(y)dy,Kt(x)=(4π∣t∣)−n/2e−isgn(t)nπ/4ei∣x∣2/(4t)と表される。
証明. 乗数の表示は命題 2.4による。ℜz>0とし、z−n/2はこの右半平面で正実軸上の正の値から定める。§E15.3 命題 4.1と Fourier 反転§E15.3 定理 2.1により、正実数zについて
∫Rne−4π2z∣ξ∣2e2πix⋅ξdξ=(4πz)−n/2e−∣x∣2/(4z).左辺は右半平面の任意のコンパクト集合上で、被積分関数とz微分が可積分な Gauss 関数に支配される。したがって左辺はzについて正則である。右辺も正則なので、正則関数の一致の定理から等式はℜz>0で成立する。
ε>0とし、z=ε+itを代入した核をKε,tとおく。u0∈Sならば、Fourier 反転と Fubini の定理により
F−1(e−4π2(ε+it)∣ξ∣2u0)(x)=∫RnKε,t(x−y)u0(y)dy.実際、Fubini の適用に要る絶対積分は∥u0∥1∫e−4π2ε∣ξ∣2dξ<∞である。ε↓0のとき、左辺はu0∈L1による優収束で各点eitΔu0(x)へ収束し、Plancherel の定理によってL2でも同じ極限をもつ。
右辺ではKε,t(x)→Kt(x)であり、
∣Kε,t(x)∣=(4π)−n/2(ε2+t2)−n/4exp(−4(ε2+t2)ε∣x∣2)≤(4π∣t∣)−n/2.したがって∣u0∣による優収束から核表示を得る。複素べきの極限はarg(ε+it)→sgn(t)π/2なので、表示された位相係数となる。▨
定理 4.2.n≥1、t=0とする。このとき、核積分はL1(Rn)から有界連続関数への複素線形作用素を定め、
∥eitΔu0∥∞≤(4π∣t∣)−n/2∥u0∥1を満たす。この作用素はL1∩L2上でL2のユニタリ群に一致する。さらに1≤p≤2、1/p+1/p′=1に対して、eitΔはLp(Rn)からLp′(Rn)への有界作用素として定まり、
∥eitΔu0∥p′≤(4π∣t∣)−n(1/p−1/2)∥u0∥pを満たす。
証明.∣Kt(x−y)∣=(4π∣t∣)−n/2なので、u0∈L1に対する積分は絶対収束し、端点評価を満たす。x→x0のとき被積分関数は各点収束し、定数倍の∣u0∣に支配されるので、積分は連続である。
u0∈L1∩L2に対して、uj∈SをL1とL2の両方でu0に収束するように取る。この同時近似は、u0を有界な球と有界な値で切断してから滑らかな近似恒等核で畳み込むことで得られる。定理 4.1により二つの作用素はuj上で一致する。端点評価から核積分の像は一様収束し、ユニタリ性から群の像はL2で収束する。L2収束列のほとんど至る所収束する部分列を取ると、二つの極限は一致する。
二つの作用素が共通部分で一致するので、L1+L2上にTt(g+h)=Kt∗g+eitΔhと定めることができる。二つの分解の差はL1∩L2に属するので、値は分解によらない。Ttの端点ノルムはL1→L∞で(4π∣t∣)−n/2以下、L2→L2で一である。Lebesgue 測度は半有限なので§E12.36 定理 2.1を適用することができる。1/p=1−θ/2とおくと1/p′=θ/2となり、補間定数は
((4π∣t∣)−n/2)1−θ=(4π∣t∣)−n(1/p−1/2)である。p=1,2は端点評価そのものである。▨
例 4.3.Rn上の初期値u0(x)=e−∣x∣2に対する自由解は
u(t,x)=(1+4it)−n/2exp(−1+4it∣x∣2)である。複素べきはt=0の値一から連続に定める。実際、§E15.3 命題 4.1によりu0(ξ)=πn/2e−π2∣ξ∣2である。この式に自由乗数を掛け、定理 4.1の証明で得た複素 Gauss 積分を用いると表示を得る。絶対値は
∣u(t,x)∣=(1+16t2)−n/4exp(−1+16t2∣x∣2)なので、∥u(t)∥∞=(1+16t2)−n/4であり、指数関数の空間的な尺度は(1+16t2)1/2である。
例 4.4.Ω⊂Rnを空でない有界C∞開集合とし、e=0を Dirichlet Laplace 作用素ADの固有関数、ADe=λeとする。このとき、自由空間の分散評価と同じ形の
∥e−itADv∥∞≤C∣t∣−n/2∥v∥1(t=0)という全時刻の評価は成立しない。実際、ADe=λeを繰り返し用いると、任意の整数r≥1についてe∈D(ADr)である。§E15.23 系 6.1によりe∈H2r(Ω)である。2r>n/2となるrを取り、§E15.10 系 4.2を適用するとeは有界な連続代表をもつ。またΩの測度は有限なのでe∈L1である。e−itADe=e−itλeだから、仮定した評価から0<∥e∥∞≤C∣t∣−n/2∥e∥1を得る。右辺は∣t∣→∞で零へ収束して矛盾する。たとえばΩ=(0,π)、e(x)=sinxでは∥e∥∞=1、∥e∥1=2であり、任意の有限なCに対して十分大きい∣t∣で1≤2C∣t∣−1/2は破れる。