§E13.33Hill 方程式

最終更新

周期的な係数をもつ微分方程式では、Floquet 理論によって一周期ごとの解の変化を記述する。しかし、直線全体で微分作用素を考えると、一周期の解の性質が作用素のスペクトルにどのように反映されるかが問題となる。周期的な実ポテンシャルをもつ Hill 方程式は、この関係を調べる基本的な対象であり、スペクトルは閉区間の帯の和集合として現れる。たとえば、Mathieu 方程式−u′′+2hcos⁡(2x)u=λu-u''+2h\cos(2x)u=\lambda uでは、十分小さい非零のhhに対し、帯の間にスペクトルを含まない有限の開区間が生じる。本記事では、Hill 方程式の帯スペクトルとその基本的な性質を扱う。

1 周期作用素と左右の Weyl 解

定義 1.1.T>0T>0とし、q∈Lloc1(R;R)q\in L^1_{\mathrm{loc}}(\R;\R)はほとんど至る所でq(x+T)=q(x)q(x+T)=q(x)を満たすとする。Hilbert 空間L2(R;C)L^2(\R;\C)上の作用素AqA_qを

D(Aq)={u∈L2(R):u,u′∈ACloc(R), −u′′+qu∈L2(R)},Aqu=−u′′+qu\begin{split} D(A_q)&=\{u\in L^2(\R):u,u'\in AC_{\mathrm{loc}}(\R),\ -u''+qu\in L^2(\R)\},\\ A_qu&=-u''+qu \end{split}

によって定める。AqA_qをポテンシャルqqの 周期 Schrödinger 作用素 (periodic Schrödinger operator) という。局所絶対連続な関数はその連続な代表元で扱う。

定理 1.2.T>0T>0、q∈Lloc1(R;R)q\in L^1_{\mathrm{loc}}(\R;\R)とし、qqはTT周期であるとする。このとき、周期 Schrödinger 作用素AqA_qについて、次が成り立つ。

  1. AqA_qは下に有界な自己共役作用素である。対応する閉形式は D(aq)=H1(R),aq(u,v)=∫R(u′v′‾+quv‾)dxD(\mathfrak a_q)=H^1(\R),\qquad \mathfrak a_q(u,v)=\int_\R\bigl(u'\overline{v'}+qu\overline v\bigr)dx であり、Hill 方程式−u′′+qu=zu-u''+qu=zuの両端±∞\pm\inftyは極限点型である。
  2. z∈C∖Rz\in\C\setminus\Rとする。このとき、Hill 方程式のモノドロミー行列M(z)M(z)の二つの乗数は相異なり、一方ρ+(z)\rho_+(z)は∣ρ+(z)∣<1|\rho_+(z)|<1、他方ρ−(z)=ρ+(z)−1\rho_-(z)=\rho_+(z)^{-1}は∣ρ−(z)∣>1|\rho_-(z)|>1を満たす。乗数ρ+\rho_+の Floquet 解は右端の Weyl 解であり、乗数ρ−\rho_-の Floquet 解は左端の Weyl 解である。特にΔ(z)∉[−2,2]\Delta(z)\notin[-2,2]である。

証明.Q±=∫0Tq±(x) dxQ_\pm=\int_0^Tq_\pm(x)\,dxと置く。長さ一の任意の区間は、高々NT=⌈1/T⌉+2N_T=\lceil1/T\rceil+2個の区間[nT,(n+1)T][nT,(n+1)T]に含まれる。周期性から各区間上のq±q_\pmの積分はQ±Q_\pmなので、

sup⁡x∈R∫xx+1q±(t) dt≤NTQ±<∞\sup_{x\in\R}\int_x^{x+1}q_\pm(t)\,dt\le N_TQ_\pm<\infty

である。§E13.31 補題 3.2をr=q+r=q_+に適用すると、すべてのu∈H1(R)u\in H^1(\R)について∫q+∣u∣2<∞\int q_+|u|^2<\inftyとなる。従って§E13.31 定理 3.5の自然形式域はH1(R)H^1(\R)に等しい。同定理によりaq\mathfrak a_qは稠密に定義された下に有界な閉形式である。§E13.31 定理 3.3により両端は極限点型であり、§E13.31 定理 3.6により形式に対応する自己共役作用素はAqA_qに一致する。

主張 (2)を示す。∣ρ∣=1|\rho|=1の乗数が存在すると仮定し、対応する非零 Floquet 解をuuとする。この解は

(u(T),u′(T))=ρ(u(0),u′(0))(u(T),u'(T))=\rho(u(0),u'(0))

を満たす。ρ=eiη\rho=e^{i\eta}、0≤η<π0\le\eta<\piならばR=IR=Iとし、ρ=−eiη\rho=-e^{i\eta}、0≤η<π0\le\eta<\piならばR=−IR=-Iとする。いずれも§E13.30 定理 3.4の結合型条件である。uuの[0,T][0,T]への制限はこの自己共役実現の域に属し、非実数zzの固有関数となる。これは、自己共役作用素のスペクトルが実軸に含まれるという§E12.15 命題 3.1 (1)に反する。

§E13.13 命題 5.2 (2)により乗数の積は一であるから、一方だけが単位円の内側にある。重根ならば乗数は11または−1-1となるので、二根は相異なる。§E13.13 命題 5.2 (3)により、内側の乗数に対応する解とその導関数は右端で指数的に減衰する。外側の乗数に対応する解は左端で指数的に減衰する。各端点は極限点型であるから、§E13.32 補題 1.2により、これらはそれぞれの端点の Weyl 解である。Δ(z)∈[−2,2]\Delta(z)\in[-2,2]ならば乗数の二次方程式の根は単位円上にあるので、Δ(z)∉[−2,2]\Delta(z)\notin[-2,2]である。▨

補題 1.3.T>0T>0、q∈Lloc1(R;R)q\in L^1_{\mathrm{loc}}(\R;\R)とし、qqはTT周期であるとする。このとき、任意のλ∈R\lambda\in\Rに対して、Hill 方程式−u′′+qu=λu-u''+qu=\lambda uの非零解はL2(R)L^2(\R)に属さない。従って周期 Schrödinger 作用素AqA_qは固有値をもたない。

証明. 非零解uuの初期値をv=(u(0),u′(0))T≠0v=(u(0),u'(0))^{\mathsf T}\ne0と置く。初期値(1,0)T,(0,1)T(1,0)^{\mathsf T},(0,1)^{\mathsf T}で標準化した基本解をθ,φ\theta,\varphiとし、ξ=(ξ1,ξ2)T∈C2\xi=(\xi_1,\xi_2)^{\mathsf T}\in\C^2に対して

Q(ξ)=∫0T∣θ(x,λ)ξ1+φ(x,λ)ξ2∣2dxQ(\xi)=\int_0^T|\theta(x,\lambda)\xi_1+\varphi(x,\lambda)\xi_2|^2dx

と定める。Q(ξ)=0Q(\xi)=0ならば連続性によって被積分関数は区間全体で零となる。従って、その解の初期値ξ\xiは零である。従ってQQは正定値であり、定数c>0c>0が存在して、すべてのξ∈C2\xi\in\C^2に対しQ(ξ)≥c∥ξ∥2Q(\xi)\ge c\|\xi\|^2が成り立つ。§E13.13 補題 1.2 (3)により、すべてのn∈Zn\in\Zについて

∫nT(n+1)T∣u(x)∣2dx=Q(M(λ)nv)\int_{nT}^{(n+1)T}|u(x)|^2dx=Q(M(\lambda)^nv)

となる。

∣Δ(λ)∣<2|\Delta(\lambda)|<2ならば、M(λ)M(\lambda)は相異なる絶対値一の固有値をもつ。固有ベクトルを列とする可逆行列をSSとすると、S−1MnvS^{-1}M^nvのノルムは∥S−1v∥>0\|S^{-1}v\|>0に等しい。従って∥Mnv∥\|M^nv\|はn∈Zn\in\Zによらない正の下界をもつ。従って、一周期ごとの質量をn∈Zn\in\Zについて足した和は発散する。

∣Δ(λ)∣>2|\Delta(\lambda)|>2ならば、乗数は実数ρ,ρ−1\rho,\rho^{-1}で∣ρ∣>1|\rho|>1としてよい。固有ベクトルの座標でvvの第一成分が零でなければ、∥Mnv∥\|M^nv\|はn→+∞n\to+\inftyで指数的に増大する。第一成分が零ならば第二成分は零でなく、n→−∞n\to-\inftyで指数的に増大する。従ってこの場合も質量の和は発散する。

Δ(λ)=2ε\Delta(\lambda)=2\varepsilon、ε∈{1,−1}\varepsilon\in\{1,-1\}ならば、二行二列の特性方程式から

M=ε(I+N),N2=0,Mnv=εn(v+nNv)(n∈Z)M=\varepsilon(I+N),\qquad N^2=0,\qquad M^nv=\varepsilon^n(v+nNv)\quad(n\in\Z)

と書くことができる。Nv=0Nv=0ならば∥Mnv∥=∥v∥>0\|M^nv\|=\|v\|>0である。Nv≠0Nv\ne0ならば∥v+nNv∥≥∣n∣∥Nv∥−∥v∥\|v+nNv\|\ge |n|\|Nv\|-\|v\|である。いずれの場合も一周期の質量の和は発散する。

判別式の三つの場合のいずれでも一周期の質量の和は発散するので、u∉L2(R)u\notin L^2(\R)である。AqA_qの固有関数は定義域の条件により全直線上の二乗可積分解であり、自己共役性により固有値は実数なので、固有値は存在しない。▨

2 Green 核と Weyl 行列

補題 2.1.T>0T>0、q∈Lloc1(R;R)q\in L^1_{\mathrm{loc}}(\R;\R)とし、qqはTT周期であるとする。Hill 方程式の判別式についてλ∈R\lambda\in\R、∣Δ(λ)∣>2|\Delta(\lambda)|>2とする。u+u_+を絶対値が一より小さい乗数の Floquet 解、u−u_-をその逆数の乗数の Floquet 解とし、W=u−u+′−u−′u+W=u_-u_+'-u_-'u_+と置く。このときW≠0W\ne0であり、

Gλ(x,t)=−u−(min⁡{x,t})u+(max⁡{x,t})WG_\lambda(x,t)=-\frac{u_-(\min\{x,t\})u_+(\max\{x,t\})}{W}

はL2(R)L^2(\R)上の有界積分作用素を定める。この作用素は(Aq−λI)−1(A_q-\lambda I)^{-1}に等しく、従ってλ∈ρ(Aq)\lambda\in\rho(A_q)である。

証明.∣ρ∣<1|\rho|<1の乗数をρ\rhoとし、a=−T−1log⁡∣ρ∣>0a=-T^{-1}\log|\rho|>0と置く。§E13.13 命題 5.2 (3)により二解は一次独立であるからW≠0W\ne0である。同命題の指数評価と左向きの整数周期移動から、定数C0>0C_0>0が存在して、すべてのx∈Rx\in\Rに対し

∣u+(x)∣≤C0e−ax,∣u−(x)∣≤C0eax(x∈R)|u_+(x)|\le C_0e^{-ax},\qquad |u_-(x)|\le C_0e^{ax}\qquad(x\in\R)

である。従ってC=C02/∣W∣C=C_0^2/|W|とすると

∣Gλ(x,t)∣≤Ce−a∣x−t∣,sup⁡x∫R∣Gλ(x,t)∣dt≤2Ca,sup⁡t∫R∣Gλ(x,t)∣dx≤2Ca.|G_\lambda(x,t)|\le Ce^{-a|x-t|},\qquad \sup_x\int_\R|G_\lambda(x,t)|dt\le\frac{2C}{a},\qquad \sup_t\int_\R|G_\lambda(x,t)|dx\le\frac{2C}{a}.

各xxについてGλ(x,⋅)∈L2(R)G_\lambda(x,\cdot)\in L^2(\R)なので、f∈L2(R)f\in L^2(\R)に対する核積分Rf(x)=∫Gλ(x,t)f(t)dtRf(x)=\int G_\lambda(x,t)f(t)dtは絶対収束する。Cauchy–Schwarz の不等式により

∣Rf(x)∣2≤(∫R∣Gλ(x,t)∣dt)(∫R∣Gλ(x,t)∣∣f(t)∣2dt).|Rf(x)|^2\le\left(\int_\R|G_\lambda(x,t)|dt\right) \left(\int_\R|G_\lambda(x,t)||f(t)|^2dt\right).

xxについて積分し、非負関数の積分順序を交換すると、∥Rf∥2≤(2C/a)∥f∥2\|Rf\|_2\le(2C/a)\|f\|_2を得る。

f∈L2(R)f\in L^2(\R)がコンパクトな台をもつとする。核積分を

Rf(x)=−1W{u+(x)∫−∞xu−(t)f(t)dt+u−(x)∫x∞u+(t)f(t)dt}Rf(x)=-\frac1W\left\{ u_+(x)\int_{-\infty}^xu_-(t)f(t)dt+ u_-(x)\int_x^\infty u_+(t)f(t)dt\right\}

と書く。積分の被積分関数は局所可積分である。微分するとf(x)f(x)を含む二項が相殺され、二度目の微分ではu−u+′−u−′u+=Wu_-u_+'-u_-'u_+=Wによって

(Rf)′′=(q−λ)Rf−f(Rf)''=(q-\lambda)Rf-f

を得る。従ってRf,(Rf)′∈ACloc(R)Rf,(Rf)'\in AC_{\mathrm{loc}}(\R)である。既にRf∈L2Rf\in L^2を示したので、Rf∈D(Aq)Rf\in D(A_q)、(Aq−λI)Rf=f(A_q-\lambda I)Rf=fとなる。

一般のf∈L2f\in L^2に対し、コンパクト台のfnf_nをffへL2L^2収束させると、Rfn→RfRf_n\to RfおよびAqRfn=fn+λRfn→f+λRfA_qRf_n=f_n+\lambda Rf_n\to f+\lambda RfがL2L^2で成り立つ。AqA_qは自己共役で閉作用素なので、Rf∈D(Aq)Rf\in D(A_q)、(Aq−λI)Rf=f(A_q-\lambda I)Rf=fである。補題 1.3によりAq−λIA_q-\lambda Iは単射であるから、R=(Aq−λI)−1R=(A_q-\lambda I)^{-1}となる。▨

命題 2.2.T>0T>0、q∈Lloc1(R;R)q\in L^1_{\mathrm{loc}}(\R;\R)とし、qqはTT周期であるとする。Hill 方程式−u′′+qu=zu-u''+qu=zuの基本解を

θ(0,z)=φ′(0,z)=1,θ′(0,z)=φ(0,z)=0\theta(0,z)=\varphi'(0,z)=1,\qquad \theta'(0,z)=\varphi(0,z)=0

で標準化し、モノドロミー行列を

M(z)=(a(z)b(z)c(z)d(z)),Δ(z)=a(z)+d(z)M(z)=\begin{pmatrix}a(z)&b(z)\\c(z)&d(z)\end{pmatrix},\qquad \Delta(z)=a(z)+d(z)

と書く。z∈C∖Rz\in\C\setminus\Rに対し、ρ+\rho_+を絶対値が一より小さい乗数、ρ−\rho_-を他方の乗数とし、s=ρ+−ρ−s=\rho_+-\rho_-と置く。このとき、左右の Weyl 解をψ±=θ+m±φ\psi_\pm=\theta+m_\pm\varphiと正規化すると、

m±(z)=ρ±(z)−a(z)b(z),s(z)2=Δ(z)2−4m_\pm(z)=\frac{\rho_\pm(z)-a(z)}{b(z)},\qquad s(z)^2=\Delta(z)^2-4

である。基本解(θ,φ)(\theta,\varphi)に対応する Weyl 行列をM\mathcal Mと書くと、

M(z)=−1s(z)(b(z)(d(z)−a(z))/2(d(z)−a(z))/2−c(z)).\mathcal M(z)=-\frac1{s(z)} \begin{pmatrix} b(z)&(d(z)-a(z))/2\\ (d(z)-a(z))/2&-c(z) \end{pmatrix}.

実点λ\lambdaが∣Δ(λ)∣≠2|\Delta(\lambda)|\ne2を満たすならば、上半平面からのM\mathcal Mはλ\lambdaの複素近傍へ正則に延長される。∣Δ(λ)∣<2|\Delta(\lambda)|<2のときには、上半平面からの境界値について

m+(λ+i0)=α(λ)+iβ(λ),m−(λ+i0)=α(λ)−iβ(λ),β(λ)>0m_+(\lambda+i0)=\alpha(\lambda)+i\beta(\lambda),\qquad m_-(\lambda+i0)=\alpha(\lambda)-i\beta(\lambda),\qquad \beta(\lambda)>0

であり、α,β\alpha,\betaは各開いた帯の内部で連続である。このとき

ℑM(λ+i0)=12β(λ)(1α(λ)α(λ)α(λ)2+β(λ)2)\Im\mathcal M(\lambda+i0) =\frac1{2\beta(\lambda)} \begin{pmatrix} 1&\alpha(\lambda)\\ \alpha(\lambda)&\alpha(\lambda)^2+\beta(\lambda)^2 \end{pmatrix}

は正定値である。

証明.z∈C∖Rz\in\C\setminus\Rとし、b(z)=0b(z)=0と仮定する。解φ(⋅,z)\varphi(\cdot,z)はφ(0,z)=φ(T,z)=0\varphi(0,z)=\varphi(T,z)=0を満たすので、§E13.30 定理 3.4の Dirichlet 自己共役実現の非実固有関数となる。これは矛盾である。従ってb(z)≠0b(z)\ne0である。

定理 1.2 (2)により二つの乗数は相異なり、ψ±\psi_\pmは乗数ρ±\rho_\pmの Floquet 解である。初期値(1,m±)T(1,m_\pm)^{\mathsf T}に対する固有ベクトルの第一行はa+bm±=ρ±a+bm_\pm=\rho_\pmなので、表示されたm±m_\pmの式を得る。ρ++ρ−=a+d\rho_++\rho_-=a+d、ρ+ρ−=ad−bc=1\rho_+\rho_-=ad-bc=1より

m+−m−=sb,m++m−=d−ab,m−m+=−cb.m_+-m_-=\frac{s}{b},\qquad m_++m_-=\frac{d-a}{b},\qquad m_-m_+=-\frac cb.

§E13.32 定義 5.1の式へ代入すると、M\mathcal Mの表示を得る。

g=Δ2−4g=\Delta^2-4と置き、g(z0)≠0g(z_0)\ne0とする。§E13.13 系 5.4によりa,b,c,d,ga,b,c,d,gは整関数である。r02=g(z0)r_0^2=g(z_0)を満たすr0≠0r_0\ne0を一つ取り、z0z_0の円板を∣g(z)/g(z0)−1∣<1/2|g(z)/g(z_0)-1|<1/2となるように取る。この円板上で

r(z)=r0exp⁡(12Log⁡g(z)g(z0))r(z)=r_0\exp\left(\frac12\operatorname{Log}\frac{g(z)}{g(z_0)}\right)

と置く。§E5.3 命題 2.10と§E5.3 命題 1.8によりrrは正則であり、r2=gr^2=gである。従って二次方程式の二根(Δ+r)/2(\Delta+r)/2、(Δ−r)/2(\Delta-r)/2は相異なる正則関数となる。非実点では絶対値が一より小さい根の選択が局所で一定なので、ρ±\rho_\pmは正則である。実点λ\lambdaでもΔ(λ)2−4≠0\Delta(\lambda)^2-4\ne0ならば、上半円板で選んだ二根を円板全体へ正則に延長することができる。ssはこの円板を小さくすれば零にならず、M\mathcal Mの式にはbbによる除算が残っていない。従ってM\mathcal Mは円板全体へ正則に延長される。

∣Δ(λ)∣<2|\Delta(\lambda)|<2とする。もしb(λ)=0b(\lambda)=0ならば、実数a,da,dはad=1ad=1を満たすので、∣a+d∣≥2|a+d|\ge2となる。従ってb(λ)≠0b(\lambda)\ne0である。二つの乗数は互いに共役な非実数なので、m−(λ+i0)=m+(λ+i0)‾m_-(\lambda+i0)=\overline{m_+(\lambda+i0)}であり、ℑm+(λ+i0)≠0\Im m_+(\lambda+i0)\ne0である。上半平面でのℑm+>0\Im m_+>0は§E13.32 命題 1.5によるから、連続な境界値の虚部は正である。m+=α+iβm_+=\alpha+i\beta、m−=α−iβm_-=\alpha-i\betaを§E13.32 定義 5.1へ代入すると、主張の虚部の式を得る。右辺の第一対角成分は1/(2β)>01/(2\beta)>0、行列式は1/4>01/4>0なので、正定値である。▨

注意 2.3.命題 2.2の基本解とモノドロミー行列を用い、実点λ\lambdaがΔ(λ)2≠4\Delta(\lambda)^2\ne4を満たすとする。同命題の証明により、λ\lambdaの複素近傍で二つの乗数ρ±\rho_\pmは相異なる正則関数として取ることができる。この近傍では

P±(z)=M(z)−ρ∓(z)Iρ±(z)−ρ∓(z)P_\pm(z)=\frac{M(z)-\rho_\mp(z)I}{\rho_\pm(z)-\rho_\mp(z)}

は各固有直線への正則な射影である。実際、相異なる固有値に対応する固有ベクトルの基底において、P±P_\pmは対応する成分を保ち、他方の成分を零にする。

P±(λ)v≠0P_\pm(\lambda)v\ne0となる定数ベクトルvvを選ぶと、P±(z)vP_\pm(z)vはλ\lambdaの近傍で零にならない正則な固有ベクトルである。第一成分が非零の点では、第一成分で割ることによって対応する固有解をu=θ+mφu=\theta+m\varphiと正規化することができる。

実点λ\lambdaでb(λ)=0b(\lambda)=0となり、乗数に対応する非零の固有解uuがu(0)=0u(0)=0を満たす場合には、初期値の一意性によりu′(0)≠0u'(0)\ne0である。その固有直線に対応するP±P_\pmと、P±(λ)v≠0P_\pm(\lambda)v\ne0を満たす定数ベクトルvvを取ると、P±(z)vP_\pm(z)vの第二成分はλ\lambdaの近傍で非零である。したがって、第二成分で割ることによって、対応する固有解を

u=φ+nθ,u′(0)=1,n=u(0)u′(0)u=\varphi+n\theta,\qquad u'(0)=1,\qquad n=\frac{u(0)}{u'(0)}

と正規化することができ、nnは正則である。二つの正規化をともに用いることができる点ではn=1/mn=1/mである。従って、mmの極は固有直線の座標をnnへ替えることによって扱うことができる。命題 2.2のbbを分母に含まない表示ではs=ρ+−ρ−≠0s=\rho_+-\rho_-\ne0であり、この座標の変更は Weyl 行列M\mathcal Mに特異点を生じない。

命題 2.4.T>0T>0、q∈Lloc1(R;R)q\in L^1_{\mathrm{loc}}(\R;\R)とし、qqはTT周期であるとする。周期 Schrödinger 作用素AqA_qの基本解を00で標準化し、対応するスペクトル行列をΩ\Omega、Weyl 行列をM\mathcal Mとする。このとき、開集合

O={λ∈R:∣Δ(λ)∣<2}O=\{\lambda\in\R:|\Delta(\lambda)|<2\}

上では

dΩ(λ)=1πℑM(λ+i0) dλd\Omega(\lambda)=\frac1\pi\Im\mathcal M(\lambda+i0)\,d\lambda

であり、この密度は連続かつ正定値である。特に、OOに含まれる空でない開区間IIではΩ(I)≠0\Omega(I)\ne0である。

証明.命題 2.2により、OOの各点の近傍でM\mathcal Mは上半平面から正則に延長され、その虚部は実軸を含めて連続である。§E13.32 命題 7.1を適用すると、主張の密度を得る。正定値性も命題 2.2による。IIに含まれる非退化なコンパクト区間JJと、非零の定数ベクトルv∈C2v\in\C^2を取ると、

v∗Ω(J)v=1π∫Jv∗ℑM(λ+i0)v dλ>0v^*\Omega(J)v=\frac1\pi\int_Jv^*\Im\mathcal M(\lambda+i0)v\,d\lambda>0

であり、Ω(I)≠0\Omega(I)\ne0である。▨

3 帯スペクトルと絶対連続性

定義 3.1. 複素 Hilbert 空間HH上の自己共役作用素AAと、その射影値測度EEを取る。すべてのf∈Hf\in Hについて、有限正測度

μf(B)=⟨E(B)f,f⟩(B⊆R は Borel 集合)\mu_f(B)=\langle E(B)f,f\rangle\qquad(B\subseteq\R\text{ は Borel 集合})

が Lebesgue 測度に関して絶対連続であるとき、AAのスペクトルは 純絶対連続 (purely absolutely continuous) であるという。§E10.14 定理 2.3により、すべてのf∈Hf\in Hに対して非負可積分関数wfw_fが存在し、すべての Borel 集合BBについてμf(B)=∫Bwf(λ)dλ\mu_f(B)=\int_Bw_f(\lambda)d\lambdaが成り立つことは、AAのスペクトルが純絶対連続であるための必要十分条件である。

定理 3.2.T>0T>0、q∈Lloc1(R;R)q\in L^1_{\mathrm{loc}}(\R;\R)とし、qqはTT周期であるとする。このとき、周期 Schrödinger 作用素AqA_qのスペクトルは

σ(Aq)={λ∈R:∣Δ(λ)∣≤2}\sigma(A_q)=\{\lambda\in\R:|\Delta(\lambda)|\le2\}

であり、固有値をもたず、純絶対連続である。00で標準化した基本解に対するスペクトル行列Ω\Omegaも、実軸全体で Lebesgue 測度に関して絶対連続である。

一周期[0,T][0,T]の周期型固有値をλ0,λ1,…\lambda_0,\lambda_1,\ldots、反周期型固有値をλ0′,λ1′,…\lambda'_0,\lambda'_1,\ldotsとし、それぞれ重複度を込めて並べると、

λ0<λ0′≤λ1′<λ1≤λ2<λ2′≤λ3′<⋯ ,\lambda_0<\lambda'_0\le\lambda'_1<\lambda_1\le\lambda_2<\lambda'_2\le\lambda'_3<\cdots,σ(Aq)=⋃m∈N≥0([λ2m,λ2m′]∪[λ2m+1′,λ2m+1]).\sigma(A_q)=\bigcup_{m\in\N} \left([\lambda_{2m},\lambda'_{2m}]\cup[\lambda'_{2m+1},\lambda_{2m+1}]\right).

隣り合う帯端λ2m′,λ2m+1′\lambda'_{2m},\lambda'_{2m+1}またはλ2m+1,λ2m+2\lambda_{2m+1},\lambda_{2m+2}が一致すると、その間のギャップは空になる。帯端の集合{λ∈R:Δ(λ)2=4}\{\lambda\in\R:\Delta(\lambda)^2=4\}は実軸上に集積点をもたない。

証明.λ∈R\lambda\in\R、∣Δ(λ)∣>2|\Delta(\lambda)|>2ならば、補題 2.1によりλ∈ρ(Aq)\lambda\in\rho(A_q)である。AqA_qは自己共役なので、

σ(Aq)⊆{λ∈R:∣Δ(λ)∣≤2}\sigma(A_q)\subseteq\{\lambda\in\R:|\Delta(\lambda)|\le2\}

を得る。

EEをAqA_qの射影値測度とする。§E13.32 定理 5.5により、L2(R)L^2(\R)からL2(Ω)L^2(\Omega)へのユニタリ写像UUが存在してUE(I)U−1UE(I)U^{-1}は1I\mathbf1_Iによる乗算である。有界 Borel 集合IIについて、Ω(I)=0\Omega(I)=0ならばこの乗算は零である。逆にΩ(I)≠0\Omega(I)\ne0ならば、非負行列Ω(I)\Omega(I)のある対角成分は正であるから、対応する座標ベクトルをeje_jとして1Iej\mathbf1_Ie_jはL2(Ω)L^2(\Omega)の非零元となる。従ってE(I)≠0E(I)\ne0である。

∣Δ(λ)∣<2|\Delta(\lambda)|<2とする。λ\lambdaの任意の近傍は、{∣Δ∣<2}\{|\Delta|<2\}に含まれる有界な空でない開区間IIを含む。命題 2.4によりΩ(I)≠0\Omega(I)\ne0であり、従ってE(I)≠0E(I)\ne0である。§E12.15 命題 3.1 (1)によりλ∈σ(Aq)\lambda\in\sigma(A_q)となる。

§E13.28 定理 5.4と§E13.28 系 5.6 (1)をp=w=1p=w=1に適用すると、主張の固有値の配列と閉区間の和を得る。§E13.28 系 5.6 (2)により各帯の内部では∣Δ∣<2|\Delta|<2であり、配列の狭義不等号から各帯の長さは正である。従って各帯端は帯内部の点の極限である。スペクトルは閉集合なので、すべての帯端もスペクトルに属する。

§E13.13 系 5.4によりΔ2−4\Delta^2-4は整関数である。定理 1.2 (2)により非実点では零にならないので、恒等的に零ではない。§E5.8 命題 2.2により零点は孤立している。実軸上の有界閉区間に無限個の零点があれば、その区間内の集積点も連続性によって零点となり、孤立性に反する。従って帯端は各有界区間に有限個しかなく、全体でも可算集合である。

補題 1.3によりAqA_qは固有値をもたない。§E12.15 命題 3.1 (3)により、すべてのλ∈R\lambda\in\RでE({λ})=0E(\{\lambda\})=0である。UE({λ})U−1UE(\{\lambda\})U^{-1}が1{λ}\mathbf1_{\{\lambda\}}による乗算であることから、Ω({λ})=0\Omega(\{\lambda\})=0も成り立つ。ギャップでは§E12.15 命題 3.1 (1)によりEEが零であり、ユニタリ表示からΩ\Omegaも零である。帯内部では命題 2.4によりΩ\Omegaは絶対連続である。残りの帯端は可算集合で、各点の測度が零なので、帯端全体のΩ\Omega測度も零である。従ってΩ\Omegaは実軸全体で絶対連続である。

Lebesgue 零 Borel 集合N⊆RN\subseteq\Rを取る。Ω(N∩[−r,r])=0\Omega(N\cap[-r,r])=0なので、ユニタリ表示によりE(N∩[−r,r])=0E(N\cap[-r,r])=0である。r→∞r\to\inftyとすると、任意のf∈L2(R)f\in L^2(\R)に対し⟨E(N)f,f⟩=0\langle E(N)f,f\rangle=0を得る。従ってAqA_qのスペクトルは純絶対連続である。▨

例 3.3.q(x)=x2q(x)=x^2は実数値の滑らかな関数で下に有界であるが、周期関数ではない。§E13.32 例 8.4により、L2(R)L^2(\R)上の調和振動子A=−d2/dx2+x2A=-d^2/dx^2+x^2は、固有値2n+12n+1、n∈N≥0n\in\Nに対応する完全正規直交固有系をもつ。従ってAAのスペクトルは純点スペクトルである。実数値で滑らかかつ下に有界なポテンシャルでも、周期性を仮定しなければ、固有値が存在してスペクトルが純絶対連続でない場合がある。

4 自由方程式と Mathieu 方程式

例 4.1.T>0T>0とし、q=0q=0とする。このとき、自由方程式−u′′=λu-u''=\lambda uの判別式は

Δ(λ)=2∑n∈N≥0(−1)nT2nλn(2n)!=2cos⁡(Tλ)\Delta(\lambda)=2\sum_{n\in\N}\frac{(-1)^nT^{2n}\lambda^n}{(2n)!} =2\cos(T\sqrt\lambda)

である。右辺の余弦は平方根の符号に依存せず、λ=0\lambda=0でも級数によって定まる。周期 Schrödinger 作用素A0A_0はσ(A0)=[0,∞)\sigma(A_0)=[0,\infty)をもち、正のエネルギーにあるすべてのギャップは閉じている。一方、0∈σ(A0)0\in\sigma(A_0)であっても、解u(x)=xu(x)=xは実軸上で有界ではない。

証明.§E13.13 例 5.6のh=0h=0の基本解をx=Tx=Tで評価すると、λ=k2>0\lambda=k^2>0ではΔ(λ)=2cos⁡(kT)\Delta(\lambda)=2\cos(kT)、λ=0\lambda=0ではΔ(0)=2\Delta(0)=2、λ=−κ2<0\lambda=-\kappa^2<0ではΔ(λ)=2cosh⁡(κT)>2\Delta(\lambda)=2\cosh(\kappa T)>2を得る。余弦の冪級数により、表示されたΔ\Deltaの級数も得る。従って{∣Δ∣≤2}=[0,∞)\{|\Delta|\le2\}=[0,\infty)であり、定理 3.2によりスペクトルは[0,∞)[0,\infty)である。

同じ基本解とその導関数をx=Tx=Tで評価すると、n∈N≥1n\in\NNに対してM((nπ/T)2)=(−1)nIM((n\pi/T)^2)=(-1)^nIである。§E13.28 定理 5.4 (2)により、各正の帯端は二重であり、対応する二つの端点は一致する。従って正のギャップはすべて閉じる。λ=0\lambda=0の解u=xu=xは非有界なので、閉じた帯の全点で全解が有界であるとは限らない。▨

例 4.2.h∈Rh\in\Rに対して、周期π\piの Mathieu 方程式

−u′′+2hcos⁡(2x)u=λu-u''+2h\cos(2x)u=\lambda u

を考える。このとき、δ>0\delta>0とC2C^2級関数ae,ao:(−δ,δ)→Ra_e,a_o:(-\delta,\delta)\to\Rが存在し、偶関数の反周期解に対応する固有値ae(h)a_e(h)と、奇関数の反周期解に対応する固有値ao(h)a_o(h)は

ae(h)=1+h+O(h2),ao(h)=1−h+O(h2)(h→0)a_e(h)=1+h+O(h^2),\qquad a_o(h)=1-h+O(h^2)\qquad(h\to0)

を満たす。0<∣h∣<δ0<|h|<\deltaでは、最初の有限ギャップは

(min⁡{ae(h),ao(h)}, max⁡{ae(h),ao(h)})\bigl(\min\{a_e(h),a_o(h)\},\ \max\{a_e(h),a_o(h)\}\bigr)

であり、その長さは2∣h∣+O(h2)2|h|+O(h^2)である。

証明.L=π/2L=\pi/2とし、初期値

θ(0;h,λ)=1,θ′(0;h,λ)=0,φ(0;h,λ)=0,φ′(0;h,λ)=1\theta(0;h,\lambda)=1,\quad\theta'(0;h,\lambda)=0,\qquad \varphi(0;h,\lambda)=0,\quad\varphi'(0;h,\lambda)=1

をもつ解を取る。一次系の係数行列は(h,λ)(h,\lambda)について滑らかであり、有界なパラメータ集合上の各偏導関数は[0,π][0,\pi]上で一様に有界である。§E13.8 定理 4.1をr=2r=2で適用すると、

Fe(h,λ)=θ(L;h,λ),Fo(h,λ)=φ′(L;h,λ)F_e(h,\lambda)=\theta(L;h,\lambda),\qquad F_o(h,\lambda)=\varphi'(L;h,\lambda)

はC2C^2級である。h=0h=0では両者ともcos⁡(Lλ)\cos(L\sqrt\lambda)に等しいので、

Fe(0,1)=Fo(0,1)=0,∂λFe(0,1)=∂λFo(0,1)=−π4≠0.F_e(0,1)=F_o(0,1)=0,\qquad \partial_\lambda F_e(0,1)=\partial_\lambda F_o(0,1)=-\frac\pi4\ne0.

§E4.8 系 3.1により、零点を与えるC1C^1級の枝ae,aoa_e,a_oが存在する。各枝は

aj′(h)=−∂hFj(h,aj(h))∂λFj(h,aj(h))(j=e,o)a_j'(h)=-\frac{\partial_hF_j(h,a_j(h))}{\partial_\lambda F_j(h,a_j(h))}\qquad(j=e,o)

を満たす。近傍を小さくすれば分母は零でなく、FjF_jはC2C^2級、aja_jはC1C^1級なので、右辺はC1C^1級である。従ってae,aoa_e,a_oはC2C^2級である。

ue(x,h)=θ(x;h,ae(h))u_e(x,h)=\theta(x;h,a_e(h))、uo(x,h)=φ(x;h,ao(h))u_o(x,h)=\varphi(x;h,a_o(h))と置く。ポテンシャルはx=0x=0について偶であるから、初期値問題の一意性によりueu_eは偶、uou_oは奇である。またcos⁡(2(π−x))=cos⁡2x\cos(2(\pi-x))=\cos2xであり、ue(L,h)=0u_e(L,h)=0なので、−ue(π−x,h)-u_e(\pi-x,h)とue(x,h)u_e(x,h)はx=Lx=Lで同じ値と導関数をもつ。従って両者は一致する。同様にuo′(L,h)=0u_o'(L,h)=0からuo(π−x,h)=uo(x,h)u_o(\pi-x,h)=u_o(x,h)を得る。よって

(uj(π,h),uj′(π,h))=−(uj(0,h),uj′(0,h))(j=e,o)(u_j(\pi,h),u_j'(\pi,h))=-(u_j(0,h),u_j'(0,h))\qquad(j=e,o)

であり、各h∈(−δ,δ)h\in(-\delta,\delta)に対して、ae(h),ao(h)a_e(h),a_o(h)は反周期型固有値である。

j∈{e,o}j\in\{e,o\}を固定し、u=uj(⋅,h)u=u_j(\cdot,h)、u˙=∂huj(⋅,h)\dot u=\partial_hu_j(\cdot,h)と置く。§E13.8 定理 3.8のパラメータ微分の方程式から

(−d2dx2+2hcos⁡2x−aj(h))u˙=(aj′(h)−2cos⁡2x)u\left(-\frac{d^2}{dx^2}+2h\cos2x-a_j(h)\right)\dot u =\bigl(a_j'(h)-2\cos2x\bigr)u

である。uuとu˙\dot uは[0,L][0,L]の両端で同じ混合分離型条件を満たす。§E13.26 定理 1.4 (1)と§E13.26 定理 1.4 (2)により境界項は零となり、

aj′(h)∫0Luj(x,h)2dx=∫0L2cos⁡2x uj(x,h)2dxa_j'(h)\int_0^Lu_j(x,h)^2dx =\int_0^L2\cos2x\,u_j(x,h)^2dx

を得る。ue(x,0)=cos⁡xu_e(x,0)=\cos x、uo(x,0)=sin⁡xu_o(x,0)=\sin xであり、

∫0π/2cos⁡2x dx=∫0π/2sin⁡2x dx=π4,\int_0^{\pi/2}\cos^2x\,dx=\int_0^{\pi/2}\sin^2x\,dx=\frac\pi4,∫0π/22cos⁡2xcos⁡2x dx=π4,∫0π/22cos⁡2xsin⁡2x dx=−π4.\int_0^{\pi/2}2\cos2x\cos^2x\,dx=\frac\pi4,\qquad \int_0^{\pi/2}2\cos2x\sin^2x\,dx=-\frac\pi4.

従ってae′(0)=1a_e'(0)=1、ao′(0)=−1a_o'(0)=-1である。二階導関数は00の近傍で有界なので、積分剰余の表示からae(h)=1+h+O(h2)a_e(h)=1+h+O(h^2)、ao(h)=1−h+O(h2)a_o(h)=1-h+O(h^2)を得る。特に、十分小さいh≠0h\ne0では二つの値は異なる。

[0,π][0,\pi]上の反周期自己共役実現の固有値を、重複度を込めてγ1(h)≤γ2(h)≤⋯\gamma_1(h)\le\gamma_2(h)\le\cdotsと並べる。§E13.30 定理 4.3にX={(z,−z):z∈C}X=\{(z,-z):z\in\C\}、B=0B=0を代入すると、対応する形式は

D(ah)={u∈H1(0,π):u(π)=−u(0)},ah(u,u)=∫0π(∣u′∣2+2hcos⁡2x ∣u∣2)dxD(\mathfrak a_h)=\{u\in H^1(0,\pi):u(\pi)=-u(0)\},\qquad \mathfrak a_h(u,u)=\int_0^\pi\bigl(|u'|^2+2h\cos2x\,|u|^2\bigr)dx

である。形式域はhhに依存せず、

∣ah(u,u)−a0(u,u)∣≤2∣h∣∥u∥22|\mathfrak a_h(u,u)-\mathfrak a_0(u,u)|\le2|h|\|u\|_2^2

なので、§E13.30 定理 6.1により∣γn(h)−γn(0)∣≤2∣h∣|\gamma_n(h)-\gamma_n(0)|\le2|h|である。自由方程式の反周期固有値は1,1,9,9,…1,1,9,9,\ldotsである。従ってhhが十分小さいとき、11の近傍にある反周期型固有値は重複度を込めて最初の二つだけである。ae(h),ao(h)a_e(h),a_o(h)は、h→0h\to0のとき11に収束する。また、十分小さいh≠0h\ne0に対して相異なる反周期型固有値であるから、最初の二つの反周期型固有値に一致する。

§E13.28 定理 5.4 (4)と定理 3.2により、最初の二つの反周期型固有値の間は最初の有限ギャップである。その長さは

∣ae(h)−ao(h)∣=∣2h+O(h2)∣=2∣h∣+O(h2)|a_e(h)-a_o(h)|=|2h+O(h^2)|=2|h|+O(h^2)

となる。▨

命題 4.3.T>0T>0、q∈Lloc1(R;R)q\in L^1_{\mathrm{loc}}(\R;\R)とし、qqはTT周期であるとする。このとき、任意の定数c∈Rc\in\Rに対して

D(Aq+c)=D(Aq),Aq+c=Aq+cI,Δq+c(λ)=Δq(λ−c)(λ∈C)D(A_{q+c})=D(A_q),\qquad A_{q+c}=A_q+cI,\qquad \Delta_{q+c}(\lambda)=\Delta_q(\lambda-c)\quad(\lambda\in\C)

である。従ってσ(Aq+c)=c+σ(Aq)\sigma(A_{q+c})=c+\sigma(A_q)であり、周期型・反周期型の帯端とギャップもccだけ平行移動する。

証明は演習とする(問題 5.1)。

5 演習

問題 5.1.命題 4.3の証明を完成させよ。また、一定ポテンシャルq=cq=cの周期 Schrödinger 作用素のスペクトルと、周期TTに対する帯端を求めよ。

解答.

u∈L2(R)u\in L^2(\R)、u,u′∈ACloc(R)u,u'\in AC_{\mathrm{loc}}(\R)とする。cu∈L2(R)cu\in L^2(\R)なので、−u′′+qu∈L2-u''+qu\in L^2であることは、−u′′+(q+c)u∈L2-u''+(q+c)u\in L^2であるための必要十分条件である。従って定義域は一致し、Aq+c=Aq+cIA_{q+c}=A_q+cIである。方程式

−u′′+(q+c)u=λu-u''+(q+c)u=\lambda u

は−u′′+qu=(λ−c)u-u''+qu=(\lambda-c)uに等しい。標準化初期値は共通なので、初期値問題の一意性によりMq+c(λ)=Mq(λ−c)M_{q+c}(\lambda)=M_q(\lambda-c)、従って判別式の等式を得る。定理 3.2へ代入すると、スペクトル・帯端・ギャップがccだけ平行移動する。

q=0q=0に適用すると、例 4.1によりσ(Ac)=[c,∞)\sigma(A_c)=[c,\infty)である。周期型固有値はccとc+(2mπ/T)2c+(2m\pi/T)^2(m∈N≥1m\in\NN)、反周期型固有値はc+((2m+1)π/T)2c+((2m+1)\pi/T)^2(m∈N≥0m\in\N)である。最下端ccは一重であり、残りの帯端は二重である。特にc<0c<0ならば、負のエネルギーの区間[c,0)[c,0)も純絶対連続スペクトルに含まれる。▨

参考文献

  1. Gerald Teschl, Mathematical Methods in Quantum Mechanics: With Applications to Schrödinger Operators, 2nd ed., Graduate Studies in Mathematics 157, American Mathematical Society, 2014.周期 Schrödinger 作用素の Floquet 乗数、Weyl 関数、帯スペクトルと絶対連続性を扱う参考文献である。

前提記事

14 本の記事・単元を表示