ある量子系にわずかな外力を加えると、エネルギーの値(固有値)はどう変わるでしょうか。固有値が一つだけなら、外力の大きさの冪級数で追えます。ところが同じ値の固有値が複数重なっている場合は、外力を入れた途端に分かれたり、κ \sqrt{\kappa} κ のように動いたりしえます。この記事では、そうした分かれ方を有限の大きさの行列の問題に直して調べ、自己共役な場合には分岐が起こらないことを示します。最後に、近似値から本当の固有値を上下から挟む方法も扱います。
固有射影 P ( κ ) P(\kappa) P ( κ ) が κ \kappa κ について解析的なとき、微分方程式 U ′ = [ P ′ , P ] U U'=[P',P]U U ′ = [ P ′ , P ] U の解 U ( κ ) U(\kappa) U ( κ ) (Kato の変換関数 )は P ( κ ) U ( κ ) = U ( κ ) P ( 0 ) P(\kappa)U(\kappa)=U(\kappa)P(0) P ( κ ) U ( κ ) = U ( κ ) P ( 0 ) を満たし、実の κ \kappa κ ではユニタリ です。これで κ \kappa κ ごとに動く固有空間を固定した空間に引き戻し、縮退した固有値の問題を m × m m\times m m × m 行列の解析族に縮約できます。
解析的な係数をもつ多項式の根は、0 < ∣ κ ∣ < ρ 0<|\kappa|<\rho 0 < ∣ κ ∣ < ρ では長さ p p p の巡回をなし、各巡回は κ 1 / p \kappa^{1/p} κ 1/ p のべき級数で表されます(Puiseux 展開 )。固有値が分岐する(κ \sqrt{\kappa} κ のように振る舞う)のは、この巡回の長さが 2 2 2 以上のときです。
実の κ \kappa κ で自己共役な族では巡回の長さは必ず 1 1 1 で、固有値と固有射影は交差を含めて解析的に選べます(Rellich の定理 )。
単純固有値の Rayleigh–Schrödinger 級数 は、簡約レゾルベント S S S を使った漸化式 E n = ⟨ B φ 0 , χ n − 1 ⟩ E_n=\langle B\varphi_0,\chi_{n-1}\rangle E n = ⟨ B φ 0 , χ n − 1 ⟩ 、χ n = S ( − B χ n − 1 + ∑ k = 1 n − 1 E k χ n − k ) \chi_n=S(-B\chi_{n-1}+\sum_{k=1}^{n-1}E_k\chi_{n-k}) χ n = S ( − B χ n − 1 + ∑ k = 1 n − 1 E k χ n − k ) で全次数が決まり、収束半径の下界も具体的に書けます。Hellmann–Feynman の定理 E ′ ( κ ) = ⟨ ψ ( κ ) , B ψ ( κ ) ⟩ E'(\kappa)=\langle\psi(\kappa),B\psi(\kappa)\rangle E ′ ( κ ) = ⟨ ψ ( κ ) , B ψ ( κ )⟩ は、固有ベクトルが連続でありさえすれば成り立ちます。
試行関数 ψ \psi ψ の 2 次のずれ ε 2 = ∥ ( A − η ) ψ ∥ 2 \varepsilon^2=\|(A-\eta)\psi\|^2 ε 2 = ∥ ( A − η ) ψ ∥ 2 から、固有値を上下から挟む Kato–Temple の不等式 が得られます。部分空間の上の作用素で全空間の作用素を近似する一般化強レゾルベント収束 では、スペクトルは極限で消えません。
この記事で使う事前知識
リース射影 : 閉曲線で囲んだ孤立固有値に対応する固有空間への射影を、レゾルベントの周回積分で与えたもの。→ リース射影
直交射影 : 閉部分空間への、像と核が直交する射影(P 2 = P = P ∗ P^2=P=P^* P 2 = P = P ∗ )。→ 直交射影
自己共役作用素 : A = A ∗ A=A^* A = A ∗ を満たす作用素で、スペクトルが実数になる。→ 自己共役作用素
強レゾルベント収束 : 作用素の列のレゾルベントが各ベクトルについて収束すること。→ 強レゾルベント収束
特性多項式 : 行列 T T T に対する det ( λ − T ) \det(\lambda-T) det ( λ − T ) のこと。その根が固有値。→ 特性多項式
作用素の収束と摂動 の §5 では、自己共役作用素 A A A の孤立した有限多重度の固有値 E 0 E_0 E 0 が相対有界な摂動 A κ = A + κ B A_\kappa=A+\kappa B A κ = A + κ B のもとでどう動くかを調べ、リース射影 P ( κ ) P(\kappa) P ( κ ) の解析性、単純固有値の解析性、Rayleigh–Schrödinger 級数の 2 次までを示しました。そこで残った問いが三つあります。
一つ目は縮退した固有値 です。E 0 E_0 E 0 の多重度が m ≥ 2 m\ge 2 m ≥ 2 だと、P ( κ ) P(\kappa) P ( κ ) は解析的でも、その値域の中で m m m 個の固有値が個別にどう動くかは分かりません。非対称な 2 × 2 2\times 2 2 × 2 行列 ( 0 1 κ 0 ) \begin{pmatrix}0&1\\\kappa&0\end{pmatrix} ( 0 κ 1 0 ) の固有値 ± κ \pm\sqrt\kappa ± κ のように、固有値は一般には κ \kappa κ の解析関数になりません。どのような分岐が起こりうるのか(Puiseux 展開)、自己共役性がそれをどう禁じるのか(Rellich の定理)が本記事の前半の主題です。
二つ目は級数の高次 です。量子力学の教科書の摂動論は 2 次や 3 次で止まりますが、定理として扱うには全次数の係数と収束半径が必要です。
三つ目は数値的な検証 です。変分法(Rayleigh–Ritz)は固有値の上界しか与えません。計算機で得た近似固有値を定理の一部として使うには下界も必要で、それを与えるのが Temple と Kato の不等式です。また、有限の箱で計算した結果を全空間の作用素と比べるには、異なる Hilbert 空間上の作用素の収束を定式化する必要があります。
以下 H \mathcal{H} H は複素 Hilbert 空間です。Banach 空間(とくに B ( H ) B(\mathcal{H}) B ( H ) )に値をとる関数 F F F が円板 ∣ κ ∣ < ρ |\kappa|<\rho ∣ κ ∣ < ρ 上で解析的 であるとは、各点の近くでノルム収束するべき級数に展開できることをいいます。積・和・ノルム収束する一様な極限・径路積分は解析性を保ち、⟨ φ , F ( κ ) ψ ⟩ \langle\varphi,F(\kappa)\psi\rangle ⟨ φ , F ( κ ) ψ ⟩ を考えれば、一致の定理や Cauchy の積分公式は複素関数の場合に帰着します。
本記事を通じて、作用素の収束と摂動 の設定と記号を使います。
A A A は自己共役、B B B は対称で D ( A ) ⊂ D ( B ) D(A)\subset D(B) D ( A ) ⊂ D ( B ) 、∥ B φ ∥ ≤ a ∥ A φ ∥ + b ∥ φ ∥ \|B\varphi\|\le a\|A\varphi\|+b\|\varphi\| ∥ B φ ∥ ≤ a ∥ A φ ∥ + b ∥ φ ∥ 。A κ = A + κ B A_\kappa=A+\kappa B A κ = A + κ B (D ( A κ ) = D ( A ) D(A_\kappa)=D(A) D ( A κ ) = D ( A ) )。
E 0 E_0 E 0 は A A A の孤立固有値で多重度 m < ∞ m<\infty m < ∞ 、d = d i s t ( E 0 , σ ( A ) ∖ { E 0 } ) d=\mathrm{dist}(E_0,\sigma(A)\setminus\{E_0\}) d = dist ( E 0 , σ ( A ) ∖ { E 0 }) 、Γ = { ∣ z − E 0 ∣ = r } \Gamma=\{|z-E_0|=r\} Γ = { ∣ z − E 0 ∣ = r } (0 < r < d 0<r<d 0 < r < d )、δ = min ( r , d − r ) \delta=\min(r,d-r) δ = min ( r , d − r ) 、M = sup z ∈ Γ ∥ B ( A − z ) − 1 ∥ M=\sup_{z\in\Gamma}\|B(A-z)^{-1}\| M = sup z ∈ Γ ∥ B ( A − z ) − 1 ∥ 、κ 0 = min ( 1 / M , 1 / a ) \kappa_0=\min(1/M,1/a) κ 0 = min ( 1/ M , 1/ a ) 。
P ( κ ) = − 1 2 π i ∮ Γ ( A κ − z ) − 1 d z P(\kappa)=-\frac{1}{2\pi i}\oint_\Gamma(A_\kappa-z)^{-1}dz P ( κ ) = − 2 π i 1 ∮ Γ ( A κ − z ) − 1 d z はリース射影(定義 5.1[作用素の収束と摂動] )。
定理 5.4[作用素の収束と摂動] により、実の ∣ κ ∣ < κ 0 |\kappa|<\kappa_0 ∣ κ ∣ < κ 0 では P ( κ ) P(\kappa) P ( κ ) は階数 m m m の直交射影で A κ A_\kappa A κ のスペクトル射影です。その定理は ∣ κ ∣ < κ 0 |\kappa|<\kappa_0 ∣ κ ∣ < κ 0 で述べられていますが、証明の級数 ( A κ − z ) − 1 = ( A − z ) − 1 ∑ k ≥ 0 ( − κ ) k ( B ( A − z ) − 1 ) k (A_\kappa-z)^{-1}=(A-z)^{-1}\sum_{k\ge0}(-\kappa)^k(B(A-z)^{-1})^k ( A κ − z ) − 1 = ( A − z ) − 1 ∑ k ≥ 0 ( − κ ) k ( B ( A − z ) − 1 ) k は ∣ κ ∣ < 1 / M |\kappa|<1/M ∣ κ ∣ < 1/ M で z ∈ Γ z\in\Gamma z ∈ Γ について一様にノルム収束し、D ( A ) D(A) D ( A ) 上で A κ − z = ( 1 + κ B ( A − z ) − 1 ) ( A − z ) A_\kappa-z=(1+\kappa B(A-z)^{-1})(A-z) A κ − z = ( 1 + κ B ( A − z ) − 1 ) ( A − z ) なので Γ ⊂ ρ ( A κ ) \Gamma\subset\rho(A_\kappa) Γ ⊂ ρ ( A κ ) です(ここでは a a a による制限は不要で、1 / a 1/a 1/ a は実の κ \kappa κ での自己共役性にだけ使われます)。したがって ∣ κ ∣ < 1 / M |\kappa|<1/M ∣ κ ∣ < 1/ M で P ( κ ) P(\kappa) P ( κ ) は解析的で、P ( κ ) = ∑ k κ k P k P(\kappa)=\sum_k\kappa^kP_k P ( κ ) = ∑ k κ k P k 、∥ P k ∥ ≤ r M k / δ \|P_k\|\le rM^k/\delta ∥ P k ∥ ≤ r M k / δ が成り立ちます。これに加えて、作用素 A A A を掛けても解析性が保たれることを確認します。次の補題は、A A A が非有界でも、A P ( κ ) AP(\kappa) A P ( κ ) を有界作用素として扱えるようにするためのものです(グラフノルムは、∥ φ ∥ + ∥ A φ ∥ \|\varphi\|+\|A\varphi\| ∥ φ ∥ + ∥ A φ ∥ で測るノルムのことです)。
補題 2.1 (グラフノルムでの解析性 ) ✓
∣ κ ∣ < 1 / M |\kappa|<1/M ∣ κ ∣ < 1/ M のとき、A P ( κ ) AP(\kappa) A P ( κ ) と A κ P ( κ ) = − 1 2 π i ∮ Γ z ( A κ − z ) − 1 d z A_\kappa P(\kappa)=-\frac{1}{2\pi i}\oint_\Gamma z(A_\kappa-z)^{-1}dz A κ P ( κ ) = − 2 π i 1 ∮ Γ z ( A κ − z ) − 1 d z は H \mathcal{H} H 上の有界作用素で、κ \kappa κ について解析的である。特に P ( κ ) P(\kappa) P ( κ ) は D ( A ) D(A) D ( A ) に値をとり、κ ↦ P ( κ ) ψ \kappa\mapsto P(\kappa)\psi κ ↦ P ( κ ) ψ は D ( A ) D(A) D ( A ) のグラフノルムについて解析的である。
証明(補題 2.1) 既出の定理の証明により、z ∈ Γ z\in\Gamma z ∈ Γ で ( A κ − z ) − 1 = ( A − z ) − 1 ∑ k ≥ 0 ( − κ ) k ( B ( A − z ) − 1 ) k (A_\kappa-z)^{-1}=(A-z)^{-1}\sum_{k\ge 0}(-\kappa)^k(B(A-z)^{-1})^k ( A κ − z ) − 1 = ( A − z ) − 1 ∑ k ≥ 0 ( − κ ) k ( B ( A − z ) − 1 ) k が z z z について一様にノルム収束する。左から A A A を掛けると、∥ A ( A − z ) − 1 ∥ ≤ 1 + ( ∣ E 0 ∣ + r ) / δ \|A(A-z)^{-1}\|\le 1+(|E_0|+r)/\delta ∥ A ( A − z ) − 1 ∥ ≤ 1 + ( ∣ E 0 ∣ + r ) / δ は z ∈ Γ z\in\Gamma z ∈ Γ で一様に有界なので、A ( A κ − z ) − 1 A(A_\kappa-z)^{-1} A ( A κ − z ) − 1 も同じ級数で一様に収束する解析的な関数である。これを Γ \Gamma Γ 上で積分すれば A P ( κ ) AP(\kappa) A P ( κ ) が解析的で、A A A が閉なので積分と A A A は交換する。A κ ( A κ − z ) − 1 = 1 + z ( A κ − z ) − 1 A_\kappa(A_\kappa-z)^{-1}=1+z(A_\kappa-z)^{-1} A κ ( A κ − z ) − 1 = 1 + z ( A κ − z ) − 1 と ∮ Γ d z = 0 \oint_\Gamma dz=0 ∮ Γ d z = 0 から A κ P ( κ ) A_\kappa P(\kappa) A κ P ( κ ) の表示が出る。
∎
縮退した固有値を扱う最初の障害は、固有空間 R a n P ( κ ) \mathrm{Ran}\,P(\kappa) Ran P ( κ ) が κ \kappa κ とともに動くことです。基底を κ \kappa κ について解析的に選べれば、A κ A_\kappa A κ を R a n P ( κ ) \mathrm{Ran}\,P(\kappa) Ran P ( κ ) に制限したものを m × m m\times m m × m 行列として書けます。そのための道具が変換関数です。この節では、P ( κ ) P(\kappa) P ( κ ) は円板 ∣ κ ∣ < ρ |\kappa|<\rho ∣ κ ∣ < ρ 上の解析的な有界作用素の族で、各 κ \kappa κ で P ( κ ) 2 = P ( κ ) P(\kappa)^2=P(\kappa) P ( κ ) 2 = P ( κ ) (冪等)を満たすとします。′ ' ′ は κ \kappa κ についての微分です。
補題 3.1 (冪等族の微分 ) ✓
P P ′ P = 0 PP'P=0 P P ′ P = 0 が成り立つ。さらに Q = P ′ P − P P ′ Q=P'P-PP' Q = P ′ P − P P ′ とおくと P ′ = Q P − P Q P'=QP-PQ P ′ = QP − P Q である。
証明(補題 3.1) P 2 = P P^2=P P 2 = P を微分して P ′ P + P P ′ = P ′ P'P+PP'=P' P ′ P + P P ′ = P ′ 。左右から P P P を掛けると P P ′ P + P P ′ P = P P ′ P PP'P+PP'P=PP'P P P ′ P + P P ′ P = P P ′ P なので P P ′ P = 0 PP'P=0 P P ′ P = 0 。これを使うと Q P = P ′ P − P P ′ P = P ′ P QP=P'P-PP'P=P'P QP = P ′ P − P P ′ P = P ′ P 、P Q = P P ′ P − P P ′ = − P P ′ PQ=PP'P-PP'=-PP' P Q = P P ′ P − P P ′ = − P P ′ なので、Q P − P Q = P ′ P + P P ′ = P ′ QP-PQ=P'P+PP'=P' QP − P Q = P ′ P + P P ′ = P ′ 。
∎
証明(定理 3.2) 存在。 0 < ρ ′ < ρ 0<\rho'<\rho 0 < ρ ′ < ρ を固定し、q = sup ∣ κ ∣ ≤ ρ ′ ∥ Q ( κ ) ∥ q=\sup_{|\kappa|\le\rho'}\|Q(\kappa)\| q = sup ∣ κ ∣ ≤ ρ ′ ∥ Q ( κ ) ∥ とする。U 0 = 1 U_0=1 U 0 = 1 、U n + 1 ( κ ) = 1 + ∫ 0 κ Q ( ζ ) U n ( ζ ) d ζ U_{n+1}(\kappa)=1+\int_0^\kappa Q(\zeta)U_n(\zeta)\,d\zeta U n + 1 ( κ ) = 1 + ∫ 0 κ Q ( ζ ) U n ( ζ ) d ζ (線分に沿う積分)と定めると、帰納法で ∥ U n + 1 ( κ ) − U n ( κ ) ∥ ≤ ( q ∣ κ ∣ ) n + 1 / ( n + 1 ) ! \|U_{n+1}(\kappa)-U_n(\kappa)\|\le(q|\kappa|)^{n+1}/(n+1)! ∥ U n + 1 ( κ ) − U n ( κ ) ∥ ≤ ( q ∣ κ ∣ ) n + 1 / ( n + 1 )! が成り立つ(∣ ζ ∣ ≤ ∣ κ ∣ |\zeta|\le|\kappa| ∣ ζ ∣ ≤ ∣ κ ∣ 上で前段の評価を積分する)。円板上の解析関数の原始関数は解析的なので各 U n U_n U n は解析的で、U n U_n U n は ∣ κ ∣ ≤ ρ ′ |\kappa|\le\rho' ∣ κ ∣ ≤ ρ ′ で一様収束する。極限 U U U は解析的で U ( κ ) = 1 + ∫ 0 κ Q U d ζ U(\kappa)=1+\int_0^\kappa QU\,d\zeta U ( κ ) = 1 + ∫ 0 κ Q U d ζ を満たし、微分して U ′ = Q U U'=QU U ′ = Q U 。ρ ′ \rho' ρ ′ は任意である。
一意性。 解析的な W W W が W ′ = Q W W'=QW W ′ = Q W (あるいは後で使う W ′ = [ Q , W ] W'=[Q,W] W ′ = [ Q , W ] や W ′ = − W Q W'=-WQ W ′ = − W Q )を満たすなら、両辺を κ = 0 \kappa=0 κ = 0 で n n n 回微分した Leibniz の公式から W ( n + 1 ) ( 0 ) W^{(n+1)}(0) W ( n + 1 ) ( 0 ) は W ( 0 ) , … , W ( n ) ( 0 ) W(0),\dots,W^{(n)}(0) W ( 0 ) , … , W ( n ) ( 0 ) と Q Q Q の微分係数で決まる。したがって W ( 0 ) W(0) W ( 0 ) が同じなら Taylor 係数がすべて一致し、解は一致する。
可逆性。 同様に V ′ = − V Q V'=-VQ V ′ = − V Q 、V ( 0 ) = 1 V(0)=1 V ( 0 ) = 1 の解析的な解 V V V を取る。( V U ) ′ = − V Q U + V Q U = 0 (VU)'=-VQU+VQU=0 ( V U ) ′ = − V Q U + V Q U = 0 なので V U = 1 VU=1 V U = 1 。W = U V W=UV W = U V は W ′ = Q U V − U V Q = [ Q , W ] W'=QUV-UVQ=[Q,W] W ′ = Q U V − U V Q = [ Q , W ] 、W ( 0 ) = 1 W(0)=1 W ( 0 ) = 1 を満たし、定数 1 1 1 も同じ方程式の解なので、一意性から U V = 1 UV=1 U V = 1 。よって U − 1 = V U^{-1}=V U − 1 = V 。
P U = U P ( 0 ) PU=UP(0) P U = U P ( 0 ) 。 X = V P U X=VPU X = V P U とおくと、補題 3.1 の P ′ = Q P − P Q P'=QP-PQ P ′ = QP − P Q を使って
X ′ = − V Q P U + V ( Q P − P Q ) U + V P Q U = 0 X'=-VQPU+V(QP-PQ)U+VPQU=0 X ′ = − V QP U + V ( QP − P Q ) U + V P Q U = 0 なので X ( κ ) = P ( 0 ) X(\kappa)=P(0) X ( κ ) = P ( 0 ) 、すなわち P ( κ ) U ( κ ) = U ( κ ) P ( 0 ) P(\kappa)U(\kappa)=U(\kappa)P(0) P ( κ ) U ( κ ) = U ( κ ) P ( 0 ) 。
ユニタリ性。 実の κ \kappa κ で P ( κ ) ∗ = P ( κ ) P(\kappa)^*=P(\kappa) P ( κ ) ∗ = P ( κ ) なら、実軸に沿った微分も自己共役で P ′ ∗ = P ′ P'^*=P' P ′ ∗ = P ′ 、したがって Q ∗ = P P ′ − P ′ P = − Q Q^*=PP'-P'P=-Q Q ∗ = P P ′ − P ′ P = − Q 。実軸上で ( U ∗ U ) ′ = U ∗ Q ∗ U + U ∗ Q U = 0 (U^*U)'=U^*Q^*U+U^*QU=0 ( U ∗ U ) ′ = U ∗ Q ∗ U + U ∗ Q U = 0 なので U ∗ U = 1 U^*U=1 U ∗ U = 1 。U U U は可逆なので U ∗ = U − 1 U^*=U^{-1} U ∗ = U − 1 。
∎
U ( κ ) U(\kappa) U ( κ ) は「R a n P ( 0 ) \mathrm{Ran}\,P(0) Ran P ( 0 ) を R a n P ( κ ) \mathrm{Ran}\,P(\kappa) Ran P ( κ ) へ最短距離で回す」作用素です。生成子 Q Q Q は R a n P \mathrm{Ran}\,P Ran P と R a n ( 1 − P ) \mathrm{Ran}(1-P) Ran ( 1 − P ) を入れ替える成分だけをもつので(P Q P = 0 PQP=0 P QP = 0 、( 1 − P ) Q ( 1 − P ) = 0 (1-P)Q(1-P)=0 ( 1 − P ) Q ( 1 − P ) = 0 が 補題 3.1 から従う)、R a n P ( κ ) \mathrm{Ran}\,P(\kappa) Ran P ( κ ) の内部で余計な回転を起こしません。物理では断熱定理の「断熱的な時間発展」がこの U U U です。
系 3.3 (縮退固有値の有限行列への縮約 ) ✓
§2 の設定で、∣ κ ∣ < 1 / M |\kappa|<1/M ∣ κ ∣ < 1/ M において U ( κ ) U(\kappa) U ( κ ) を P ( κ ) P(\kappa) P ( κ ) の変換関数とし、M = R a n P ( 0 ) = K e r ( A − E 0 ) \mathcal{M}=\mathrm{Ran}\,P(0)=\mathrm{Ker}(A-E_0) M = Ran P ( 0 ) = Ker ( A − E 0 ) (m m m 次元)上の作用素
T ( κ ) = U ( κ ) − 1 A κ P ( κ ) U ( κ ) ∣ M T(\kappa)=U(\kappa)^{-1}A_\kappa P(\kappa)U(\kappa)\big|_{\mathcal{M}} T ( κ ) = U ( κ ) − 1 A κ P ( κ ) U ( κ ) M を考える。T ( κ ) T(\kappa) T ( κ ) は M \mathcal{M} M を M \mathcal{M} M に写し、M \mathcal{M} M の正規直交基底で表した m × m m\times m m × m 行列は κ \kappa κ について解析的である。実の ∣ κ ∣ < κ 0 |\kappa|<\kappa_0 ∣ κ ∣ < κ 0 では T ( κ ) T(\kappa) T ( κ ) は自己共役で、その固有値(重複度込み)は σ ( A κ ) ∩ ( E 0 − r , E 0 + r ) \sigma(A_\kappa)\cap(E_0-r,E_0+r) σ ( A κ ) ∩ ( E 0 − r , E 0 + r ) の固有値(重複度込み)に一致する。さらに φ ∈ M \varphi\in\mathcal{M} φ ∈ M について A κ U ( κ ) φ = U ( κ ) T ( κ ) φ A_\kappa U(\kappa)\varphi=U(\kappa)T(\kappa)\varphi A κ U ( κ ) φ = U ( κ ) T ( κ ) φ である。
証明(系 3.3) 補題 2.1 から A κ P ( κ ) A_\kappa P(\kappa) A κ P ( κ ) は有界解析的で、U , U − 1 U,U^{-1} U , U − 1 も解析的なので T ( κ ) T(\kappa) T ( κ ) は解析的である。U U U は M \mathcal{M} M を R a n P ( κ ) \mathrm{Ran}\,P(\kappa) Ran P ( κ ) の上に写し(定理 3.2 )、A κ P ( κ ) A_\kappa P(\kappa) A κ P ( κ ) は P ( κ ) P(\kappa) P ( κ ) と可換なので R a n P ( κ ) \mathrm{Ran}\,P(\kappa) Ran P ( κ ) を保ち、U − 1 U^{-1} U − 1 は R a n P ( κ ) \mathrm{Ran}\,P(\kappa) Ran P ( κ ) を M \mathcal{M} M に戻す。実の κ \kappa κ では U U U はユニタリ、A κ P ( κ ) A_\kappa P(\kappa) A κ P ( κ ) は自己共役(P ( κ ) P(\kappa) P ( κ ) は A κ A_\kappa A κ のスペクトル射影)なので T ( κ ) T(\kappa) T ( κ ) は自己共役で、A κ ∣ R a n P ( κ ) A_\kappa|_{\mathrm{Ran}P(\kappa)} A κ ∣ Ran P ( κ ) とユニタリ同値である。後者の固有値が σ ( A κ ) ∩ ( E 0 − r , E 0 + r ) \sigma(A_\kappa)\cap(E_0-r,E_0+r) σ ( A κ ) ∩ ( E 0 − r , E 0 + r ) であることは既出の定理の 3 による。最後の等式は U U − 1 = 1 UU^{-1}=1 U U − 1 = 1 と A κ P ( κ ) = A κ A_\kappa P(\kappa)=A_\kappa A κ P ( κ ) = A κ (R a n P ( κ ) \mathrm{Ran}\,P(\kappa) Ran P ( κ ) 上)から従う。
∎
こうして無限次元の問題は、解析的な行列族の固有値 の問題に帰着しました。次の 2 節ではこの有限次元の問題を解きます。
行列 T ( κ ) T(\kappa) T ( κ ) の固有値は特性多項式 det ( λ − T ( κ ) ) \det(\lambda-T(\kappa)) det ( λ − T ( κ )) の根なので、問題は「係数が κ \kappa κ の解析関数であるモニック多項式の根が κ \kappa κ にどう依存するか」です。κ \kappa κ の関数としての根の振る舞いは、根が重なる点(分岐点)を除けば解析的で、分岐点のまわりを一周すると根が入れ替わります。この入れ替わりの巡回の長さが分岐の次数になります。
この節の証明は長く、代数(多項式環)と位相(被覆空間)の道具を使います。定理の主張と、その後の例だけを読んでも、次の節以降の流れは追えます。補題 重根の除去 は、重根があると根の追跡が難しくなるので、重根のない多項式に取り替えておくためのものです。定理の証明は、根を κ \kappa κ の周りで一周ずつ追いかけ、一周したあとの根の入れ替わり方(置換)を調べる、という方針で進みます。
補題 4.1 (根の大きさの評価 ) ✓
λ m + c 1 λ m − 1 + ⋯ + c m = 0 \lambda^m+c_1\lambda^{m-1}+\dots+c_m=0 λ m + c 1 λ m − 1 + ⋯ + c m = 0 の根は ∣ λ ∣ ≤ max { 1 , ∑ j ∣ c j ∣ } |\lambda|\le\max\{1,\sum_j|c_j|\} ∣ λ ∣ ≤ max { 1 , ∑ j ∣ c j ∣ } を満たす。
証明(補題 4.1) ∣ λ ∣ > 1 |\lambda|>1 ∣ λ ∣ > 1 なら ∣ λ ∣ m = ∣ ∑ j c j λ m − j ∣ ≤ ∑ j ∣ c j ∣ ∣ λ ∣ m − 1 |\lambda|^m=|\sum_jc_j\lambda^{m-j}|\le\sum_j|c_j||\lambda|^{m-1} ∣ λ ∣ m = ∣ ∑ j c j λ m − j ∣ ≤ ∑ j ∣ c j ∣∣ λ ∣ m − 1 なので ∣ λ ∣ ≤ ∑ j ∣ c j ∣ |\lambda|\le\sum_j|c_j| ∣ λ ∣ ≤ ∑ j ∣ c j ∣ 。
∎
補題 4.2 (重根の除去 ) ✓
p ( κ , λ ) = λ m + c 1 ( κ ) λ m − 1 + ⋯ + c m ( κ ) p(\kappa,\lambda)=\lambda^m+c_1(\kappa)\lambda^{m-1}+\dots+c_m(\kappa) p ( κ , λ ) = λ m + c 1 ( κ ) λ m − 1 + ⋯ + c m ( κ ) の係数が κ = 0 \kappa=0 κ = 0 の近くで解析的とする。このとき ρ > 0 \rho>0 ρ > 0 と、係数が ∣ κ ∣ < ρ |\kappa|<\rho ∣ κ ∣ < ρ で解析的なモニック多項式 r ( κ , λ ) r(\kappa,\lambda) r ( κ , λ ) があり、0 < ∣ κ ∣ < ρ 0<|\kappa|<\rho 0 < ∣ κ ∣ < ρ の各 κ \kappa κ で、r ( κ , ⋅ ) r(\kappa,\cdot) r ( κ , ⋅ ) の根はすべて単根で、その集合は p ( κ , ⋅ ) p(\kappa,\cdot) p ( κ , ⋅ ) の根の集合に一致する。さらに p = ∏ i q i e i p=\prod_iq_i^{e_i} p = ∏ i q i e i 、r = ∏ i q i r=\prod_iq_i r = ∏ i q i と分解でき、q i q_i q i たちは 0 < ∣ κ ∣ < ρ 0<|\kappa|<\rho 0 < ∣ κ ∣ < ρ で共通根をもたない。
証明(補題 4.2) κ = 0 \kappa=0 κ = 0 で収束するべき級数の環 R = C { κ } R=\mathbb{C}\{\kappa\} R = C { κ } は、零でない元が κ k u \kappa^ku κ k u (u u u は定数項が 0 0 0 でない単元)と一意に書けるので一意分解整域で、その分数体を K K K とする。p p p を K [ λ ] K[\lambda] K [ λ ] でモニックな既約因子の積 ∏ i q i e i \prod_iq_i^{e_i} ∏ i q i e i (q i q_i q i は互いに異なる)に分解すると、Gauss の補題(定理 3.4[多項式環と体の拡大] )により各因子は R [ λ ] R[\lambda] R [ λ ] のモニック多項式に取れる(最高次係数の積が 1 1 1 なので単元で割ってモニックにできる)。r = ∏ i q i r=\prod_iq_i r = ∏ i q i とおく。K K K の標数は 0 0 0 なので既約多項式 q i q_i q i は自分の導関数と互いに素で、異なる q i , q j q_i,q_j q i , q j も互いに素だから、r r r と ∂ λ r \partial_\lambda r ∂ λ r は K [ λ ] K[\lambda] K [ λ ] で互いに素である。Euclid の互除法により α r + β ∂ λ r = 1 \alpha r+\beta\,\partial_\lambda r=1 α r + β ∂ λ r = 1 となる α , β ∈ K [ λ ] \alpha,\beta\in K[\lambda] α , β ∈ K [ λ ] があり、分母を払うと
α ~ ( κ , λ ) r ( κ , λ ) + β ~ ( κ , λ ) ∂ λ r ( κ , λ ) = c ( κ ) , c ∈ R , c ≠ 0 \tilde\alpha(\kappa,\lambda)\,r(\kappa,\lambda)+\tilde\beta(\kappa,\lambda)\,\partial_\lambda r(\kappa,\lambda)=c(\kappa),\qquad c\in R,\ c\ne 0 α ~ ( κ , λ ) r ( κ , λ ) + β ~ ( κ , λ ) ∂ λ r ( κ , λ ) = c ( κ ) , c ∈ R , c = 0 となる α ~ , β ~ ∈ R [ λ ] \tilde\alpha,\tilde\beta\in R[\lambda] α ~ , β ~ ∈ R [ λ ] がある。現れる有限個の級数がすべて収束し、かつ c c c が κ = 0 \kappa=0 κ = 0 以外で零点をもたない(零点の孤立性)ような ρ \rho ρ を取れば、0 < ∣ κ ∣ < ρ 0<|\kappa|<\rho 0 < ∣ κ ∣ < ρ で r ( κ , ⋅ ) r(\kappa,\cdot) r ( κ , ⋅ ) と ∂ λ r ( κ , ⋅ ) \partial_\lambda r(\kappa,\cdot) ∂ λ r ( κ , ⋅ ) は共通根をもたず、r ( κ , ⋅ ) r(\kappa,\cdot) r ( κ , ⋅ ) の根はすべて単根である。p p p と r r r の根の集合が一致するのは因子分解から明らかで、q i q_i q i と q j q_j q j の共通根は r r r の重根になるので存在しない。
∎
定理 4.3 (Puiseux 展開 ) ✓
p ( κ , λ ) p(\kappa,\lambda) p ( κ , λ ) を 補題 4.2 と同じモニック多項式とする。このとき ρ > 0 \rho>0 ρ > 0 、正の整数 p 1 , … , p s p_1,\dots,p_s p 1 , … , p s 、および ∣ t ∣ < ρ 1 / p j |t|<\rho^{1/p_j} ∣ t ∣ < ρ 1/ p j で解析的な関数 f 1 , … , f s f_1,\dots,f_s f 1 , … , f s があって、0 < ∣ κ ∣ < ρ 0<|\kappa|<\rho 0 < ∣ κ ∣ < ρ と κ \kappa κ の p j p_j p j 乗根の一つ κ 1 / p j \kappa^{1/p_j} κ 1/ p j について、
{ λ : p ( κ , λ ) = 0 } = { f j ( ω κ 1 / p j ) : 1 ≤ j ≤ s , ω p j = 1 } \{\lambda:\ p(\kappa,\lambda)=0\}=\bigl\{f_j(\omega\kappa^{1/p_j}):\ 1\le j\le s,\ \omega^{p_j}=1\bigr\} { λ : p ( κ , λ ) = 0 } = { f j ( ω κ 1/ p j ) : 1 ≤ j ≤ s , ω p j = 1 } が成り立ち、すべての f j ( ω κ 1 / p j ) f_j(\omega\kappa^{1/p_j}) f j ( ω κ 1/ p j ) は互いに異なる。すなわち p ( κ , ⋅ ) p(\kappa,\cdot) p ( κ , ⋅ ) の根は、κ 1 / p j \kappa^{1/p_j} κ 1/ p j のべき級数 ∑ k ≥ 0 a j , k κ k / p j \sum_{k\ge 0}a_{j,k}\kappa^{k/p_j} ∑ k ≥ 0 a j , k κ k / p j (Puiseux 級数)の p j p_j p j 個の値の組(巡回)に分かれる。
証明(定理 4.3) 補題 4.2 の ρ \rho ρ と r r r を取り(必要なら ρ \rho ρ を小さくして係数が閉円板 ∣ κ ∣ ≤ ρ |\kappa|\le\rho ∣ κ ∣ ≤ ρ で有界になるようにする)、n = deg r n=\deg r n = deg r とする。半平面 Ω = { s ∈ C : R e s < log ρ } \Omega=\{s\in\mathbb{C}:\mathrm{Re}\,s<\log\rho\} Ω = { s ∈ C : Re s < log ρ } は s ↦ e s s\mapsto e^s s ↦ e s で穴あき円板 0 < ∣ κ ∣ < ρ 0<|\kappa|<\rho 0 < ∣ κ ∣ < ρ に写る。
段階 1(根の局所的な解析性)。 s 1 ∈ Ω s_1\in\Omega s 1 ∈ Ω で r ( e s 1 , ⋅ ) r(e^{s_1},\cdot) r ( e s 1 , ⋅ ) の n n n 個の単根を λ 1 , … , λ n \lambda_1,\dots,\lambda_n λ 1 , … , λ n とし、互いに交わらない円板 D j = D ( λ j , ε ) D_j=D(\lambda_j,\varepsilon) D j = D ( λ j , ε ) を取る。根の連続性(定理 7.3[偏角の原理と Rouché の定理] )から、s 1 s_1 s 1 の円板近傍 W W W があって、s ∈ W s\in W s ∈ W では各 D j D_j D j に r ( e s , ⋅ ) r(e^s,\cdot) r ( e s , ⋅ ) の根がちょうど一つ μ j ( s ) \mu_j(s) μ j ( s ) あり、他に根はない。単根 μ j \mu_j μ j での z ∂ λ r / r z\,\partial_\lambda r/r z ∂ λ r / r の留数は μ j \mu_j μ j なので
μ j ( s ) = 1 2 π i ∮ ∂ D j z ∂ λ r ( e s , z ) r ( e s , z ) d z \mu_j(s)=\frac{1}{2\pi i}\oint_{\partial D_j}\frac{z\,\partial_\lambda r(e^s,z)}{r(e^s,z)}\,dz μ j ( s ) = 2 π i 1 ∮ ∂ D j r ( e s , z ) z ∂ λ r ( e s , z ) d z で、被積分関数は s s s について解析的(∂ D j \partial D_j ∂ D j 上で分母は 0 0 0 にならない)だから μ j \mu_j μ j は W W W 上で解析的である。
段階 2(大域的な根の関数)。 X = { ( s , λ ) ∈ Ω × C : r ( e s , λ ) = 0 } X=\{(s,\lambda)\in\Omega\times\mathbb{C}: r(e^s,\lambda)=0\} X = {( s , λ ) ∈ Ω × C : r ( e s , λ ) = 0 } と射影 π ( s , λ ) = s \pi(s,\lambda)=s π ( s , λ ) = s を考える。段階 1 から π − 1 ( W ) \pi^{-1}(W) π − 1 ( W ) は互いに交わらない開集合 { ( s , μ j ( s ) ) : s ∈ W } \{(s,\mu_j(s)):s\in W\} {( s , μ j ( s )) : s ∈ W } (= X ∩ ( W × D j ) =X\cap(W\times D_j) = X ∩ ( W × D j ) )の合併で、それぞれ π \pi π で W W W と同相なので、π \pi π は被覆写像である(定義 3.1[基本群と被覆空間] )。Ω \Omega Ω は凸なので単連結かつ局所道連結であり、持ち上げ定理(定理 5.1[基本群と被覆空間] )を恒等写像 Ω → Ω \Omega\to\Omega Ω → Ω に適用すると、基点 s 0 s_0 s 0 の上の各点 ( s 0 , λ j ) (s_0,\lambda_j) ( s 0 , λ j ) を通る連続な切断 s ↦ ( s , ν j ( s ) ) s\mapsto(s,\nu_j(s)) s ↦ ( s , ν j ( s )) がある。W W W 上では切断は一つの層に入る(W W W は連結)ので、ν j \nu_j ν j は局所的に段階 1 の μ k \mu_k μ k に一致し、解析的である。二つの切断が一点で一致すれば全体で一致する(補題 3.5[基本群と被覆空間] )ので、ν 1 ( s ) , … , ν n ( s ) \nu_1(s),\dots,\nu_n(s) ν 1 ( s ) , … , ν n ( s ) は各 s s s で互いに異なり、r ( e s , ⋅ ) r(e^s,\cdot) r ( e s , ⋅ ) の根のすべてである。
段階 3(一周による置換)。 s ↦ ν j ( s + 2 π i ) s\mapsto\nu_j(s+2\pi i) s ↦ ν j ( s + 2 π i ) も根の連続な切断なので、一意性からある σ ( j ) \sigma(j) σ ( j ) について ν σ ( j ) \nu_{\sigma(j)} ν σ ( j ) に一致する。σ \sigma σ は { 1 , … , n } \{1,\dots,n\} { 1 , … , n } の置換である。σ \sigma σ の長さ p p p の巡回 j , σ ( j ) , … , σ p − 1 ( j ) j,\sigma(j),\dots,\sigma^{p-1}(j) j , σ ( j ) , … , σ p − 1 ( j ) を一つ取り、0 < ∣ t ∣ < ρ 1 / p 0<|t|<\rho^{1/p} 0 < ∣ t ∣ < ρ 1/ p で
f ( t ) = ν j ( p log t ) f(t)=\nu_j(p\log t) f ( t ) = ν j ( p log t ) とおく。log t \log t log t の分枝を 2 π i ℓ 2\pi i\ell 2 π i ℓ だけ取り替えると p log t p\log t p log t は 2 π i p ℓ 2\pi ip\ell 2 π i p ℓ ずれるが、ν j ( s + 2 π i p ℓ ) = ν σ p ℓ ( j ) ( s ) = ν j ( s ) \nu_j(s+2\pi ip\ell)=\nu_{\sigma^{p\ell}(j)}(s)=\nu_j(s) ν j ( s + 2 π i p ℓ ) = ν σ p ℓ ( j ) ( s ) = ν j ( s ) なので f f f は一価の解析関数である。e p log t = t p e^{p\log t}=t^p e p l o g t = t p だから f ( t ) f(t) f ( t ) は r ( t p , ⋅ ) r(t^p,\cdot) r ( t p , ⋅ ) の根で、補題 4.1 から f f f は有界である。リーマンの除去可能特異点定理(命題 6.3[正則関数の強力な性質] )により f f f は ∣ t ∣ < ρ 1 / p |t|<\rho^{1/p} ∣ t ∣ < ρ 1/ p で解析的に延びる。κ = e s \kappa=e^s κ = e s と κ 1 / p = e s / p \kappa^{1/p}=e^{s/p} κ 1/ p = e s / p について f ( e 2 π i ℓ / p κ 1 / p ) = ν j ( s + 2 π i ℓ ) = ν σ ℓ ( j ) ( s ) f(e^{2\pi i\ell/p}\kappa^{1/p})=\nu_j(s+2\pi i\ell)=\nu_{\sigma^\ell(j)}(s) f ( e 2 π i ℓ / p κ 1/ p ) = ν j ( s + 2 π i ℓ ) = ν σ ℓ ( j ) ( s ) なので、この巡回の根は f ( ω κ 1 / p ) f(\omega\kappa^{1/p}) f ( ω κ 1/ p ) (ω p = 1 \omega^p=1 ω p = 1 )で尽くされる。すべての巡回についてこれを行えば主張を得る。
∎
p j = 1 p_j=1 p j = 1 の巡回は κ = 0 \kappa=0 κ = 0 で解析的な根です。p j ≥ 2 p_j\ge 2 p j ≥ 2 の巡回では、κ = 0 \kappa=0 κ = 0 のまわりを一周すると p j p_j p j 個の根が巡回的に入れ替わり、根は κ 1 / p j \kappa^{1/p_j} κ 1/ p j 程度の速さで分かれます。既出の例 5.6[作用素の収束と摂動] は p = 2 p=2 p = 2 の巡回です。p j p_j p j は κ = 0 \kappa=0 κ = 0 での根の重複度を超えません。
例 4.4 (例外点と PT 対称な 2 準位系 )
T ( κ ) = ( i κ κ − i ) T(\kappa)=\begin{pmatrix}i&\kappa\\\kappa&-i\end{pmatrix} T ( κ ) = ( i κ κ − i ) の固有値は λ = ± κ 2 − 1 \lambda=\pm\sqrt{\kappa^2-1} λ = ± κ 2 − 1 である。T ( κ ) T(\kappa) T ( κ ) は(複素)対称行列だが、Hermite 行列ではない。κ = ± 1 \kappa=\pm1 κ = ± 1 で二つの固有値は 0 0 0 に合流し、T ( ± 1 ) T(\pm1) T ( ± 1 ) は対角化できない(T ( 1 ) 2 = 0 T(1)^2=0 T ( 1 ) 2 = 0 、T ( 1 ) ≠ 0 T(1)\ne 0 T ( 1 ) = 0 )。κ = 1 \kappa=1 κ = 1 のまわりで κ − 1 = h \kappa-1=h κ − 1 = h とおくと λ = ± 2 h ( 1 + h / 2 ) 1 / 2 = ± 2 h 1 / 2 ( 1 + h 4 + O ( h 2 ) ) \lambda=\pm\sqrt{2h}\,(1+h/2)^{1/2}=\pm\sqrt2\,h^{1/2}\bigl(1+\tfrac{h}{4}+O(h^2)\bigr) λ = ± 2 h ( 1 + h /2 ) 1/2 = ± 2 h 1/2 ( 1 + 4 h + O ( h 2 ) ) で、長さ 2 2 2 の巡回の Puiseux 級数になっている。実の κ \kappa κ では、∣ κ ∣ < 1 |\kappa|<1 ∣ κ ∣ < 1 で固有値は純虚数、∣ κ ∣ > 1 |\kappa|>1 ∣ κ ∣ > 1 で実数で、例外点 κ = ± 1 \kappa=\pm1 κ = ± 1 がその境目である(PT 対称性の破れ)。κ = 0 \kappa=0 κ = 0 の近くでは固有値 ± i 1 − κ 2 \pm i\sqrt{1-\kappa^2} ± i 1 − κ 2 は解析的で、分岐は起こらない。
自己共役な族では、実の κ \kappa κ に対して固有値がすべて実数であることが、分岐を禁じます。以下の証明では、Puiseux 展開の各巡回について、実の κ \kappa κ で値がすべて実数という条件を課すと巡回の長さが 1 1 1 に限られる、という順で進めます。
定理 5.1 (Rellich の定理(行列版) ) ✓
T ( κ ) T(\kappa) T ( κ ) を ∣ κ ∣ < ρ 0 |\kappa|<\rho_0 ∣ κ ∣ < ρ 0 で解析的な m × m m\times m m × m 行列の族とし、実の κ \kappa κ では Hermite 行列であるとする。このとき ρ > 0 \rho>0 ρ > 0 と、∣ κ ∣ < ρ |\kappa|<\rho ∣ κ ∣ < ρ で解析的な関数 μ 1 , … , μ s \mu_1,\dots,\mu_s μ 1 , … , μ s (互いに異なる関数)、および解析的な行列値関数 P 1 , … , P s P_1,\dots,P_s P 1 , … , P s があって、次が成り立つ。
実の κ ∈ ( − ρ , ρ ) \kappa\in(-\rho,\rho) κ ∈ ( − ρ , ρ ) で μ j ( κ ) \mu_j(\kappa) μ j ( κ ) は実数、P j ( κ ) P_j(\kappa) P j ( κ ) は直交射影で、P j P k = δ j k P j P_jP_k=\delta_{jk}P_j P j P k = δ j k P j 、∑ j P j = 1 \sum_jP_j=1 ∑ j P j = 1 。
T ( κ ) = ∑ j μ j ( κ ) P j ( κ ) T(\kappa)=\sum_j\mu_j(\kappa)P_j(\kappa) T ( κ ) = ∑ j μ j ( κ ) P j ( κ ) がすべての ∣ κ ∣ < ρ |\kappa|<\rho ∣ κ ∣ < ρ で成り立つ。P j P_j P j の階数は κ \kappa κ によらない。
証明(定理 5.1) 固有値の解析性。 特性多項式 det ( λ − T ( κ ) ) \det(\lambda-T(\kappa)) det ( λ − T ( κ )) に 定理 4.3 を適用し、一つの巡回 ( p , f ) (p,f) ( p , f ) を取る。実の κ > 0 \kappa>0 κ > 0 に対して κ 1 / p \kappa^{1/p} κ 1/ p を正の実数 t t t に取ると、f ( ω t ) f(\omega t) f ( ω t ) (ω p = 1 \omega^p=1 ω p = 1 )はすべて Hermite 行列 T ( κ ) T(\kappa) T ( κ ) の固有値なので実数である。κ < 0 \kappa<0 κ < 0 に対しては κ 1 / p = e i π / p ∣ κ ∣ 1 / p \kappa^{1/p}=e^{i\pi/p}|\kappa|^{1/p} κ 1/ p = e iπ / p ∣ κ ∣ 1/ p と取れるので、f ( ω e i π / p t ) f(\omega e^{i\pi/p}t) f ( ω e iπ / p t ) も実数である。f ( t ) = ∑ k a k t k f(t)=\sum_ka_kt^k f ( t ) = ∑ k a k t k と展開する。t > 0 t>0 t > 0 で f ( t ) f(t) f ( t ) が実数なので、実変数の実数値関数としての Taylor 係数 a k a_k a k は実数である。t > 0 t>0 t > 0 で f ( e i π / p t ) = ∑ k a k e i π k / p t k f(e^{i\pi/p}t)=\sum_ka_ke^{i\pi k/p}t^k f ( e iπ / p t ) = ∑ k a k e iπ k / p t k も実数なので a k e i π k / p a_ke^{i\pi k/p} a k e iπ k / p も実数で、a k ≠ 0 a_k\ne 0 a k = 0 なら e i π k / p = ± 1 e^{i\pi k/p}=\pm1 e iπ k / p = ± 1 、すなわち p ∣ k p\mid k p ∣ k 。よって f ( t ) = h ( t p ) f(t)=h(t^p) f ( t ) = h ( t p ) (h h h は解析的)で、巡回の p p p 個の根 f ( ω κ 1 / p ) = h ( κ ) f(\omega\kappa^{1/p})=h(\kappa) f ( ω κ 1/ p ) = h ( κ ) はすべて等しい。巡回の根は互いに異なるので p = 1 p=1 p = 1 である。したがって 0 < ∣ κ ∣ < ρ 0<|\kappa|<\rho 0 < ∣ κ ∣ < ρ の根はすべて κ = 0 \kappa=0 κ = 0 で解析的な関数 μ 1 , … , μ s \mu_1,\dots,\mu_s μ 1 , … , μ s の値で、補題 4.2 から各 μ j \mu_j μ j の重複度は κ ≠ 0 \kappa\ne 0 κ = 0 で一定である。μ j \mu_j μ j は実軸上で実数である。
射影の解析性。 μ j − μ k \mu_j-\mu_k μ j − μ k (j ≠ k j\ne k j = k )は恒等的に 0 0 0 でない解析関数なので、ρ \rho ρ を小さくすれば 0 < ∣ κ ∣ < ρ 0<|\kappa|<\rho 0 < ∣ κ ∣ < ρ では零点をもたない。そこで
P j ( κ ) = ∏ k ≠ j T ( κ ) − μ k ( κ ) μ j ( κ ) − μ k ( κ ) P_j(\kappa)=\prod_{k\ne j}\frac{T(\kappa)-\mu_k(\kappa)}{\mu_j(\kappa)-\mu_k(\kappa)} P j ( κ ) = k = j ∏ μ j ( κ ) − μ k ( κ ) T ( κ ) − μ k ( κ ) は 0 < ∣ κ ∣ < ρ 0<|\kappa|<\rho 0 < ∣ κ ∣ < ρ で解析的である。実の κ ≠ 0 \kappa\ne 0 κ = 0 では T ( κ ) T(\kappa) T ( κ ) は相異なる固有値 μ 1 ( κ ) , … , μ s ( κ ) \mu_1(\kappa),\dots,\mu_s(\kappa) μ 1 ( κ ) , … , μ s ( κ ) をもつ Hermite 行列なので、固有空間分解して各固有ベクトルに作用させれば、P j ( κ ) P_j(\kappa) P j ( κ ) は固有値 μ j ( κ ) \mu_j(\kappa) μ j ( κ ) の固有空間への直交射影であり、特に ∥ P j ( κ ) ∥ ≤ 1 \|P_j(\kappa)\|\le 1 ∥ P j ( κ ) ∥ ≤ 1 。P j P_j P j の各成分は解析関数の商なので、κ = 0 \kappa=0 κ = 0 は高々極である。極であるとすると、Laurent 展開の最低次 κ − N C − N \kappa^{-N}C_{-N} κ − N C − N (N ≥ 1 N\ge 1 N ≥ 1 、C − N ≠ 0 C_{-N}\ne 0 C − N = 0 )から実軸上で ∥ κ N P j ( κ ) ∥ → ∥ C − N ∥ ≠ 0 \|\kappa^NP_j(\kappa)\|\to\|C_{-N}\|\ne 0 ∥ κ N P j ( κ ) ∥ → ∥ C − N ∥ = 0 となるが、左辺は ∣ κ ∣ N |\kappa|^N ∣ κ ∣ N 以下なので 0 0 0 に収束し、矛盾する。よって P j P_j P j は ∣ κ ∣ < ρ |\kappa|<\rho ∣ κ ∣ < ρ で解析的である。等式 P j P k = δ j k P j P_jP_k=\delta_{jk}P_j P j P k = δ j k P j 、∑ P j = 1 \sum P_j=1 ∑ P j = 1 、T = ∑ μ j P j T=\sum\mu_jP_j T = ∑ μ j P j は実の κ ≠ 0 \kappa\ne 0 κ = 0 で成り立つので、一致の定理により円板全体で成り立つ。実の κ = 0 \kappa=0 κ = 0 でも直交射影であることは連続性による。冪等行列の階数はトレースに等しく、トレースは κ \kappa κ の連続な整数値関数なので一定である。
∎
証明の要は、「f ( t ) f(t) f ( t ) が 2 p 2p 2 p 本の半直線の上で実数値をとる」ことから t p t^p t p の関数であることを導く部分です。非 Hermite な例 4.4 では、実の κ \kappa κ でも固有値が純虚数になりうるので、この議論が使えません。
定理 5.2 (Rellich の定理(作用素版) ) ✓
§2 の設定で、E 0 E_0 E 0 の多重度を m m m とする。κ 1 > 0 \kappa_1>0 κ 1 > 0 と、( − κ 1 , κ 1 ) (-\kappa_1,\kappa_1) ( − κ 1 , κ 1 ) 上の実解析関数 E 1 ( κ ) , … , E m ( κ ) E_1(\kappa),\dots,E_m(\kappa) E 1 ( κ ) , … , E m ( κ ) (E j ( 0 ) = E 0 E_j(0)=E_0 E j ( 0 ) = E 0 、重複を許す)、および各実 κ \kappa κ で正規直交な ψ 1 ( κ ) , … , ψ m ( κ ) ∈ D ( A ) \psi_1(\kappa),\dots,\psi_m(\kappa)\in D(A) ψ 1 ( κ ) , … , ψ m ( κ ) ∈ D ( A ) があり、κ ↦ ψ j ( κ ) \kappa\mapsto\psi_j(\kappa) κ ↦ ψ j ( κ ) は D ( A ) D(A) D ( A ) のグラフノルムについて実解析的で、
A κ ψ j ( κ ) = E j ( κ ) ψ j ( κ ) , σ ( A κ ) ∩ ( E 0 − r , E 0 + r ) = { E 1 ( κ ) , … , E m ( κ ) } A_\kappa\psi_j(\kappa)=E_j(\kappa)\psi_j(\kappa),\qquad \sigma(A_\kappa)\cap(E_0-r,E_0+r)=\{E_1(\kappa),\dots,E_m(\kappa)\} A κ ψ j ( κ ) = E j ( κ ) ψ j ( κ ) , σ ( A κ ) ∩ ( E 0 − r , E 0 + r ) = { E 1 ( κ ) , … , E m ( κ )} が成り立つ。
証明(定理 5.2) 系 3.3 の T ( κ ) T(\kappa) T ( κ ) に 定理 5.1 を適用し、μ j \mu_j μ j 、P j P_j P j を得る。P j ( κ ) P_j(\kappa) P j ( κ ) は M \mathcal{M} M 上の解析的な冪等族で実軸上で直交射影なので、定理 3.2 を M \mathcal{M} M の中で適用した変換関数 U j ( κ ) U_j(\kappa) U j ( κ ) は実軸上でユニタリで、P j ( κ ) U j ( κ ) = U j ( κ ) P j ( 0 ) P_j(\kappa)U_j(\kappa)=U_j(\kappa)P_j(0) P j ( κ ) U j ( κ ) = U j ( κ ) P j ( 0 ) 。R a n P j ( 0 ) \mathrm{Ran}\,P_j(0) Ran P j ( 0 ) の正規直交基底 e e e に対し U ( κ ) U j ( κ ) e U(\kappa)U_j(\kappa)e U ( κ ) U j ( κ ) e を考え、すべての j j j についてまとめて ψ 1 , … , ψ m \psi_1,\dots,\psi_m ψ 1 , … , ψ m と番号を付け直し、対応する μ j \mu_j μ j を E 1 , … , E m E_1,\dots,E_m E 1 , … , E m とする。U ( κ ) U(\kappa) U ( κ ) と U j ( κ ) U_j(\kappa) U j ( κ ) は実軸上でユニタリなので ψ k ( κ ) \psi_k(\kappa) ψ k ( κ ) は正規直交で、系 3.3 の最後の等式から A κ ψ k = U T U j e = μ j U U j e = E k ψ k A_\kappa\psi_k=U\,T\,U_je=\mu_jU\,U_je=E_k\psi_k A κ ψ k = U T U j e = μ j U U j e = E k ψ k 。グラフノルムでの解析性は、ψ k = P ( κ ) ψ k \psi_k=P(\kappa)\psi_k ψ k = P ( κ ) ψ k と 補題 2.1 による。κ 1 \kappa_1 κ 1 は 1 / M 1/M 1/ M 、κ 0 \kappa_0 κ 0 と 定理 5.1 の ρ \rho ρ の最小値に取ればよい。
∎
例 5.3 (交差する固有値と順序づけた固有値 )
T ( κ ) = ( κ κ κ − κ ) T(\kappa)=\begin{pmatrix}\kappa&\kappa\\\kappa&-\kappa\end{pmatrix} T ( κ ) = ( κ κ κ − κ ) の固有値は ± 2 κ \pm\sqrt2\,\kappa ± 2 κ で、κ = 0 \kappa=0 κ = 0 で交差する。Rellich の定理が与える解析的な枝は 2 κ \sqrt2\kappa 2 κ と − 2 κ -\sqrt2\kappa − 2 κ であり、大きさの順に並べた固有値 2 ∣ κ ∣ \sqrt2|\kappa| 2 ∣ κ ∣ 、− 2 ∣ κ ∣ -\sqrt2|\kappa| − 2 ∣ κ ∣ は κ = 0 \kappa=0 κ = 0 で微分可能でさえない。固有値を「下から n n n 番目」と番号づけると解析性が失われるので、Rellich の定理は番号づけを取り替えることを許した主張である。
κ λ Hermite:交差 κ Re λ 非 Hermite:分岐 左:Hermite な族では、固有値は交差しても解析的な枝(直線)として続く(Rellich)。右:非 Hermite な族 [[0,1],[κ,0]] では、実の κ に対する固有値 ±√κ が κ = 0 で分岐し、κ < 0 では実軸を離れる(破線は実部 0 を表す)。
この節では m = 1 m=1 m = 1 とし、φ 0 \varphi_0 φ 0 を E 0 E_0 E 0 の単位固有ベクトル、Q 0 = 1 − ⟨ φ 0 , ⋅ ⟩ φ 0 Q_0=1-\langle\varphi_0,\cdot\rangle\varphi_0 Q 0 = 1 − ⟨ φ 0 , ⋅ ⟩ φ 0 、S = ∫ σ ( A ) ∖ { E 0 } ( λ − E 0 ) − 1 d E A ( λ ) S=\int_{\sigma(A)\setminus\{E_0\}}(\lambda-E_0)^{-1}dE_A(\lambda) S = ∫ σ ( A ) ∖ { E 0 } ( λ − E 0 ) − 1 d E A ( λ ) を簡約レゾルベントとします(命題 5.5[作用素の収束と摂動] と同じ記号)。関数計算から ∥ S ∥ ≤ 1 / d \|S\|\le 1/d ∥ S ∥ ≤ 1/ d 、S φ 0 = 0 S\varphi_0=0 S φ 0 = 0 、ψ ∈ D ( A ) \psi\in D(A) ψ ∈ D ( A ) について S ( A − E 0 ) ψ = Q 0 ψ S(A-E_0)\psi=Q_0\psi S ( A − E 0 ) ψ = Q 0 ψ です。また r = d / 2 r=d/2 r = d /2 と取り、対応する M M M を使います(このとき δ = d / 2 \delta=d/2 δ = d /2 )。
定理 6.1 (Rayleigh–Schrödinger 級数(全次数) ) ✓
∣ κ ∣ < 1 / ( 2 M ) |\kappa|<1/(2M) ∣ κ ∣ < 1/ ( 2 M ) で n ( κ ) = ⟨ φ 0 , P ( κ ) φ 0 ⟩ ≠ 0 n(\kappa)=\langle\varphi_0,P(\kappa)\varphi_0\rangle\ne 0 n ( κ ) = ⟨ φ 0 , P ( κ ) φ 0 ⟩ = 0 であり、
χ ( κ ) = P ( κ ) φ 0 n ( κ ) , E ( κ ) = ⟨ φ 0 , A κ P ( κ ) φ 0 ⟩ n ( κ ) \chi(\kappa)=\frac{P(\kappa)\varphi_0}{n(\kappa)},\qquad E(\kappa)=\frac{\langle\varphi_0,A_\kappa P(\kappa)\varphi_0\rangle}{n(\kappa)} χ ( κ ) = n ( κ ) P ( κ ) φ 0 , E ( κ ) = n ( κ ) ⟨ φ 0 , A κ P ( κ ) φ 0 ⟩ はそれぞれ D ( A ) D(A) D ( A ) (グラフノルム)と C \mathbb{C} C に値をとる解析関数である。実の ∣ κ ∣ < min { κ 0 , 1 / ( 2 M ) } |\kappa|<\min\{\kappa_0,1/(2M)\} ∣ κ ∣ < min { κ 0 , 1/ ( 2 M )} では E ( κ ) E(\kappa) E ( κ ) は A κ A_\kappa A κ の ( E 0 − r , E 0 + r ) (E_0-r,E_0+r) ( E 0 − r , E 0 + r ) 内の唯一の固有値で、χ ( κ ) \chi(\kappa) χ ( κ ) は ⟨ φ 0 , χ ( κ ) ⟩ = 1 \langle\varphi_0,\chi(\kappa)\rangle=1 ⟨ φ 0 , χ ( κ )⟩ = 1 と正規化された固有ベクトルである。べき級数 E ( κ ) = ∑ n ≥ 0 E n κ n E(\kappa)=\sum_{n\ge 0}E_n\kappa^n E ( κ ) = ∑ n ≥ 0 E n κ n 、χ ( κ ) = ∑ n ≥ 0 χ n κ n \chi(\kappa)=\sum_{n\ge 0}\chi_n\kappa^n χ ( κ ) = ∑ n ≥ 0 χ n κ n の係数は、χ 0 = φ 0 \chi_0=\varphi_0 χ 0 = φ 0 と漸化式
E n = ⟨ B φ 0 , χ n − 1 ⟩ , χ n = S ( − B χ n − 1 + ∑ k = 1 n − 1 E k χ n − k ) ( n ≥ 1 ) E_n=\langle B\varphi_0,\chi_{n-1}\rangle,\qquad \chi_n=S\Bigl(-B\chi_{n-1}+\sum_{k=1}^{n-1}E_k\chi_{n-k}\Bigr)\qquad(n\ge 1) E n = ⟨ B φ 0 , χ n − 1 ⟩ , χ n = S ( − B χ n − 1 + k = 1 ∑ n − 1 E k χ n − k ) ( n ≥ 1 ) で定まり、これらの級数は ∣ κ ∣ < 1 / ( 2 M ) |\kappa|<1/(2M) ∣ κ ∣ < 1/ ( 2 M ) で収束する。
証明(定理 6.1) ∥ P k ∥ ≤ r M k / δ = M k \|P_k\|\le rM^k/\delta=M^k ∥ P k ∥ ≤ r M k / δ = M k から ∥ P ( κ ) − P ( 0 ) ∥ ≤ ∑ k ≥ 1 ( ∣ κ ∣ M ) k = ∣ κ ∣ M / ( 1 − ∣ κ ∣ M ) \|P(\kappa)-P(0)\|\le\sum_{k\ge1}(|\kappa|M)^k=|\kappa|M/(1-|\kappa|M) ∥ P ( κ ) − P ( 0 ) ∥ ≤ ∑ k ≥ 1 ( ∣ κ ∣ M ) k = ∣ κ ∣ M / ( 1 − ∣ κ ∣ M ) で、∣ κ ∣ M < 1 / 2 |\kappa|M<1/2 ∣ κ ∣ M < 1/2 ならこれは 1 1 1 未満である。P ( 0 ) φ 0 = φ 0 P(0)\varphi_0=\varphi_0 P ( 0 ) φ 0 = φ 0 なので ∣ n ( κ ) − 1 ∣ < 1 |n(\kappa)-1|<1 ∣ n ( κ ) − 1∣ < 1 、特に n ( κ ) ≠ 0 n(\kappa)\ne 0 n ( κ ) = 0 。補題 2.1 から χ \chi χ はグラフノルムで、E E E は複素数値として解析的である。
実の小さな κ \kappa κ では P ( κ ) P(\kappa) P ( κ ) は A κ A_\kappa A κ の階数 1 のスペクトル射影なので、A κ χ ( κ ) = E ′ χ ( κ ) A_\kappa\chi(\kappa)=E'\chi(\kappa) A κ χ ( κ ) = E ′ χ ( κ ) となる固有値 E ′ E' E ′ があり、⟨ φ 0 , χ ⟩ = 1 \langle\varphi_0,\chi\rangle=1 ⟨ φ 0 , χ ⟩ = 1 から E ′ = ⟨ φ 0 , A κ χ ⟩ = E ( κ ) E'=\langle\varphi_0,A_\kappa\chi\rangle=E(\kappa) E ′ = ⟨ φ 0 , A κ χ ⟩ = E ( κ ) である。等式 A κ χ ( κ ) = E ( κ ) χ ( κ ) A_\kappa\chi(\kappa)=E(\kappa)\chi(\kappa) A κ χ ( κ ) = E ( κ ) χ ( κ ) の両辺は κ \kappa κ について解析的(A χ A\chi A χ は 補題 2.1 により、B χ B\chi B χ は B B B がグラフノルムについて有界なので解析的)なので、一致の定理により ∣ κ ∣ < 1 / ( 2 M ) |\kappa|<1/(2M) ∣ κ ∣ < 1/ ( 2 M ) 全体で成り立つ。級数 χ = ∑ κ n χ n \chi=\sum\kappa^n\chi_n χ = ∑ κ n χ n はグラフノルムで収束するので、κ n \kappa^n κ n の係数を比べて
( A − E 0 ) χ n = − B χ n − 1 + ∑ k = 1 n E k χ n − k ( n ≥ 1 ) . (A-E_0)\chi_n=-B\chi_{n-1}+\sum_{k=1}^{n}E_k\chi_{n-k}\qquad(n\ge 1). ( A − E 0 ) χ n = − B χ n − 1 + k = 1 ∑ n E k χ n − k ( n ≥ 1 ) . ⟨ φ 0 , χ ( κ ) ⟩ = 1 \langle\varphi_0,\chi(\kappa)\rangle=1 ⟨ φ 0 , χ ( κ )⟩ = 1 から ⟨ φ 0 , χ n ⟩ = δ n 0 \langle\varphi_0,\chi_n\rangle=\delta_{n0} ⟨ φ 0 , χ n ⟩ = δ n 0 。上式と φ 0 \varphi_0 φ 0 の内積を取ると、左辺は ⟨ ( A − E 0 ) φ 0 , χ n ⟩ = 0 \langle(A-E_0)\varphi_0,\chi_n\rangle=0 ⟨( A − E 0 ) φ 0 , χ n ⟩ = 0 、右辺の和は k = n k=n k = n の項だけが残るので E n = ⟨ φ 0 , B χ n − 1 ⟩ = ⟨ B φ 0 , χ n − 1 ⟩ E_n=\langle\varphi_0,B\chi_{n-1}\rangle=\langle B\varphi_0,\chi_{n-1}\rangle E n = ⟨ φ 0 , B χ n − 1 ⟩ = ⟨ B φ 0 , χ n − 1 ⟩ (B B B の対称性)。次に上式に S S S を作用させると、S ( A − E 0 ) χ n = Q 0 χ n = χ n S(A-E_0)\chi_n=Q_0\chi_n=\chi_n S ( A − E 0 ) χ n = Q 0 χ n = χ n (n ≥ 1 n\ge1 n ≥ 1 )と S χ 0 = S φ 0 = 0 S\chi_0=S\varphi_0=0 S χ 0 = S φ 0 = 0 から χ n \chi_n χ n の漸化式を得る。
∎
低次を書き下すと E 1 = ⟨ φ 0 , B φ 0 ⟩ E_1=\langle\varphi_0,B\varphi_0\rangle E 1 = ⟨ φ 0 , B φ 0 ⟩ 、χ 1 = − S B φ 0 \chi_1=-SB\varphi_0 χ 1 = − S B φ 0 、E 2 = − ⟨ B φ 0 , S B φ 0 ⟩ E_2=-\langle B\varphi_0,SB\varphi_0\rangle E 2 = − ⟨ B φ 0 , S B φ 0 ⟩ 、
χ 2 = S B S B φ 0 − E 1 S 2 B φ 0 , E 3 = ⟨ B φ 0 , S B S B φ 0 ⟩ − E 1 ⟨ B φ 0 , S 2 B φ 0 ⟩ \chi_2=SBSB\varphi_0-E_1S^2B\varphi_0,\qquad E_3=\langle B\varphi_0,SBSB\varphi_0\rangle-E_1\langle B\varphi_0,S^2B\varphi_0\rangle χ 2 = S B S B φ 0 − E 1 S 2 B φ 0 , E 3 = ⟨ B φ 0 , S B S B φ 0 ⟩ − E 1 ⟨ B φ 0 , S 2 B φ 0 ⟩
です。2 次までは既出の命題と一致します。収束半径の下界 1 / ( 2 M ) 1/(2M) 1/ ( 2 M ) は、摂動が A A A のスペクトルの隙間 d d d に比べてどれだけ小さいかで決まり、B B B が有界なら M ≤ ∥ B ∥ / δ = 2 ∥ B ∥ / d M\le\|B\|/\delta=2\|B\|/d M ≤ ∥ B ∥/ δ = 2∥ B ∥/ d なので、∣ κ ∣ ∥ B ∥ < d / 4 |\kappa|\|B\|<d/4 ∣ κ ∣∥ B ∥ < d /4 で級数が収束します。
次の定理は、固有値の微分が摂動の期待値に等しいことを述べます。量子化学で原子核に働く力を計算するときに使われる式です。
定理 6.2 (Hellmann–Feynman の定理 ) ✓
I I I を実区間、E ( κ ) E(\kappa) E ( κ ) (κ ∈ I \kappa\in I κ ∈ I )を A κ A_\kappa A κ の固有値、ψ ( κ ) ∈ D ( A ) \psi(\kappa)\in D(A) ψ ( κ ) ∈ D ( A ) を対応する単位固有ベクトルとし、κ ↦ ψ ( κ ) \kappa\mapsto\psi(\kappa) κ ↦ ψ ( κ ) は H \mathcal{H} H のノルムで連続とする。このとき E E E は I I I で微分可能で
E ′ ( κ ) = ⟨ ψ ( κ ) , B ψ ( κ ) ⟩ . E'(\kappa)=\langle\psi(\kappa),B\psi(\kappa)\rangle . E ′ ( κ ) = ⟨ ψ ( κ ) , B ψ ( κ )⟩ .
証明(定理 6.2) κ , κ ′ ∈ I \kappa,\kappa'\in I κ , κ ′ ∈ I について、A κ ′ = A κ + ( κ ′ − κ ) B A_{\kappa'}=A_\kappa+(\kappa'-\kappa)B A κ ′ = A κ + ( κ ′ − κ ) B と A κ A_\kappa A κ 、B B B の対称性から
E ( κ ′ ) ⟨ ψ ( κ ) , ψ ( κ ′ ) ⟩ = ⟨ ψ ( κ ) , A κ ′ ψ ( κ ′ ) ⟩ = E ( κ ) ⟨ ψ ( κ ) , ψ ( κ ′ ) ⟩ + ( κ ′ − κ ) ⟨ B ψ ( κ ) , ψ ( κ ′ ) ⟩ . E(\kappa')\langle\psi(\kappa),\psi(\kappa')\rangle=\langle\psi(\kappa),A_{\kappa'}\psi(\kappa')\rangle=E(\kappa)\langle\psi(\kappa),\psi(\kappa')\rangle+(\kappa'-\kappa)\langle B\psi(\kappa),\psi(\kappa')\rangle . E ( κ ′ ) ⟨ ψ ( κ ) , ψ ( κ ′ )⟩ = ⟨ ψ ( κ ) , A κ ′ ψ ( κ ′ )⟩ = E ( κ ) ⟨ ψ ( κ ) , ψ ( κ ′ )⟩ + ( κ ′ − κ ) ⟨ B ψ ( κ ) , ψ ( κ ′ )⟩ . よって κ ′ ≠ κ \kappa'\ne\kappa κ ′ = κ で
E ( κ ′ ) − E ( κ ) κ ′ − κ ⟨ ψ ( κ ) , ψ ( κ ′ ) ⟩ = ⟨ B ψ ( κ ) , ψ ( κ ′ ) ⟩ . \frac{E(\kappa')-E(\kappa)}{\kappa'-\kappa}\,\langle\psi(\kappa),\psi(\kappa')\rangle=\langle B\psi(\kappa),\psi(\kappa')\rangle . κ ′ − κ E ( κ ′ ) − E ( κ ) ⟨ ψ ( κ ) , ψ ( κ ′ )⟩ = ⟨ B ψ ( κ ) , ψ ( κ ′ )⟩ . κ ′ → κ \kappa'\to\kappa κ ′ → κ で ⟨ ψ ( κ ) , ψ ( κ ′ ) ⟩ → 1 \langle\psi(\kappa),\psi(\kappa')\rangle\to 1 ⟨ ψ ( κ ) , ψ ( κ ′ )⟩ → 1 、右辺は ⟨ B ψ ( κ ) , ψ ( κ ) ⟩ \langle B\psi(\kappa),\psi(\kappa)\rangle ⟨ B ψ ( κ ) , ψ ( κ )⟩ に収束するので主張を得る。
∎
固有ベクトルの微分可能性も、B ψ ( κ ′ ) B\psi(\kappa') B ψ ( κ ′ ) の連続性も使っていないことに注意してください。B ψ ( κ ) B\psi(\kappa) B ψ ( κ ) は固定したベクトルで、動くのは ψ ( κ ′ ) \psi(\kappa') ψ ( κ ′ ) だけです。単純固有値なら ψ ( κ ) = P ( κ ) φ 0 / ∥ P ( κ ) φ 0 ∥ \psi(\kappa)=P(\kappa)\varphi_0/\|P(\kappa)\varphi_0\| ψ ( κ ) = P ( κ ) φ 0 /∥ P ( κ ) φ 0 ∥ は連続なので仮定は満たされます。縮退がある場合は、Rellich の枝に沿って同じことが言えます。
系 6.3 (縮退した固有値の 1 次摂動 ) ✓
定理 5.2 の記号で、E 1 ′ ( 0 ) , … , E m ′ ( 0 ) E_1'(0),\dots,E_m'(0) E 1 ′ ( 0 ) , … , E m ′ ( 0 ) は P ( 0 ) B P ( 0 ) P(0)BP(0) P ( 0 ) B P ( 0 ) を K e r ( A − E 0 ) \mathrm{Ker}(A-E_0) Ker ( A − E 0 ) に制限した Hermite 作用素の固有値(重複度込み)であり、ψ j ( 0 ) \psi_j(0) ψ j ( 0 ) はその固有ベクトルである。
証明(系 6.3) 実の κ ≠ 0 \kappa\ne 0 κ = 0 で ( A − E 0 ) ψ j ( κ ) + κ B ψ j ( κ ) − ( E j ( κ ) − E 0 ) ψ j ( κ ) = 0 (A-E_0)\psi_j(\kappa)+\kappa B\psi_j(\kappa)-(E_j(\kappa)-E_0)\psi_j(\kappa)=0 ( A − E 0 ) ψ j ( κ ) + κ B ψ j ( κ ) − ( E j ( κ ) − E 0 ) ψ j ( κ ) = 0 。P 0 = P ( 0 ) P_0=P(0) P 0 = P ( 0 ) は A A A のスペクトル射影なので ψ ∈ D ( A ) \psi\in D(A) ψ ∈ D ( A ) について P 0 ( A − E 0 ) ψ = ( A − E 0 ) P 0 ψ = 0 P_0(A-E_0)\psi=(A-E_0)P_0\psi=0 P 0 ( A − E 0 ) ψ = ( A − E 0 ) P 0 ψ = 0 。P 0 P_0 P 0 を作用させて κ \kappa κ で割ると
P 0 B ψ j ( κ ) = E j ( κ ) − E 0 κ P 0 ψ j ( κ ) . P_0B\psi_j(\kappa)=\frac{E_j(\kappa)-E_0}{\kappa}\,P_0\psi_j(\kappa). P 0 B ψ j ( κ ) = κ E j ( κ ) − E 0 P 0 ψ j ( κ ) . ψ j \psi_j ψ j はグラフノルムで連続なので B ψ j ( κ ) → B ψ j ( 0 ) B\psi_j(\kappa)\to B\psi_j(0) B ψ j ( κ ) → B ψ j ( 0 ) (相対有界性)。κ → 0 \kappa\to 0 κ → 0 として P 0 B ψ j ( 0 ) = E j ′ ( 0 ) ψ j ( 0 ) P_0B\psi_j(0)=E_j'(0)\psi_j(0) P 0 B ψ j ( 0 ) = E j ′ ( 0 ) ψ j ( 0 ) 。ψ 1 ( 0 ) , … , ψ m ( 0 ) \psi_1(0),\dots,\psi_m(0) ψ 1 ( 0 ) , … , ψ m ( 0 ) は m m m 次元空間 K e r ( A − E 0 ) \mathrm{Ker}(A-E_0) Ker ( A − E 0 ) の正規直交基底なので主張が従う。
∎
物理の教科書で「縮退のある摂動論では、まず摂動を縮退空間で対角化する」と教える手順の根拠がこの系です。P 0 B P 0 P_0BP_0 P 0 B P 0 の固有値が重なっている場合には、1 次では枝が分かれず、2 次以降で(χ \chi χ の漸化式に相当する縮約で)分かれます。
Rayleigh–Ritz の方法(系 6.4[作用素の収束と摂動] )は固有値の上界を与えます。下界を得るには、試行関数がどれだけ固有ベクトルに近いかを測る量 ε \varepsilon ε と、スペクトルの隙間についての事前情報が必要です。以下、ψ ∈ D ( A ) \psi\in D(A) ψ ∈ D ( A ) 、∥ ψ ∥ = 1 \|\psi\|=1 ∥ ψ ∥ = 1 とし、
η = ⟨ ψ , A ψ ⟩ , ε = ∥ ( A − η ) ψ ∥ , ε 2 = ∥ A ψ ∥ 2 − η 2 \eta=\langle\psi,A\psi\rangle,\qquad \varepsilon=\|(A-\eta)\psi\|,\qquad \varepsilon^2=\|A\psi\|^2-\eta^2 η = ⟨ ψ , A ψ ⟩ , ε = ∥ ( A − η ) ψ ∥ , ε 2 = ∥ A ψ ∥ 2 − η 2
とおきます。μ ψ = ⟨ ψ , E A ( ⋅ ) ψ ⟩ \mu_\psi=\langle\psi,E_A(\cdot)\psi\rangle μ ψ = ⟨ ψ , E A ( ⋅ ) ψ ⟩ は確率測度で、∫ t d μ ψ = η \int t\,d\mu_\psi=\eta ∫ t d μ ψ = η 、∫ ( t − η ) 2 d μ ψ = ε 2 \int(t-\eta)^2d\mu_\psi=\varepsilon^2 ∫ ( t − η ) 2 d μ ψ = ε 2 です。
補題 7.1 (Weinstein の評価 ) ✓
d i s t ( η , σ ( A ) ) ≤ ε \mathrm{dist}(\eta,\sigma(A))\le\varepsilon dist ( η , σ ( A )) ≤ ε である。
証明(補題 7.1) D = d i s t ( η , σ ( A ) ) D=\mathrm{dist}(\eta,\sigma(A)) D = dist ( η , σ ( A )) とすると、μ ψ \mu_\psi μ ψ は σ ( A ) \sigma(A) σ ( A ) に台をもつので ε 2 = ∫ ( t − η ) 2 d μ ψ ≥ D 2 \varepsilon^2=\int(t-\eta)^2d\mu_\psi\ge D^2 ε 2 = ∫ ( t − η ) 2 d μ ψ ≥ D 2 。
∎
これは誤差 ε \varepsilon ε について 1 次の評価です。スペクトルの隙間が分かれば 2 次に改良できます。次の定理が述べる内容は、数直線の上に、固有値 λ \lambda λ の左右にスペクトルの無い区間 ( α , λ ) (\alpha,\lambda) ( α , λ ) と ( λ , β ) (\lambda,\beta) ( λ , β ) があるとき、近似値 η \eta η の周りの区間に λ \lambda λ が入る、ということです。下の図はその配置の模式図です。
定理 7.2 (Kato–Temple の不等式 ) ✓
α < η < β \alpha<\eta<\beta α < η < β で、σ ( A ) ∩ ( α , β ) = { λ } \sigma(A)\cap(\alpha,\beta)=\{\lambda\} σ ( A ) ∩ ( α , β ) = { λ } (1 点)とする。このとき
η − ε 2 β − η ≤ λ ≤ η + ε 2 η − α . \eta-\frac{\varepsilon^2}{\beta-\eta}\ \le\ \lambda\ \le\ \eta+\frac{\varepsilon^2}{\eta-\alpha}. η − β − η ε 2 ≤ λ ≤ η + η − α ε 2 .
証明(定理 7.2) t ∈ σ ( A ) t\in\sigma(A) t ∈ σ ( A ) なら t ≤ α t\le\alpha t ≤ α 、t = λ t=\lambda t = λ 、t ≥ β t\ge\beta t ≥ β のいずれかである。g ( t ) = ( t − α ) ( t − λ ) g(t)=(t-\alpha)(t-\lambda) g ( t ) = ( t − α ) ( t − λ ) は、t ≤ α t\le\alpha t ≤ α では二つの因子がともに 0 0 0 以下(λ > α \lambda>\alpha λ > α )、t ≥ β t\ge\beta t ≥ β ではともに正(λ < β \lambda<\beta λ < β )、t = λ t=\lambda t = λ では 0 0 0 なので、σ ( A ) \sigma(A) σ ( A ) 上で g ≥ 0 g\ge 0 g ≥ 0 。t − α = ( t − η ) + ( η − α ) t-\alpha=(t-\eta)+(\eta-\alpha) t − α = ( t − η ) + ( η − α ) 、t − λ = ( t − η ) + ( η − λ ) t-\lambda=(t-\eta)+(\eta-\lambda) t − λ = ( t − η ) + ( η − λ ) と ∫ ( t − η ) d μ ψ = 0 \int(t-\eta)d\mu_\psi=0 ∫ ( t − η ) d μ ψ = 0 から
0 ≤ ∫ g d μ ψ = ε 2 + ( η − α ) ( η − λ ) , 0\le\int g\,d\mu_\psi=\varepsilon^2+(\eta-\alpha)(\eta-\lambda), 0 ≤ ∫ g d μ ψ = ε 2 + ( η − α ) ( η − λ ) , すなわち λ ≤ η + ε 2 / ( η − α ) \lambda\le\eta+\varepsilon^2/(\eta-\alpha) λ ≤ η + ε 2 / ( η − α ) 。同様に h ( t ) = ( t − β ) ( t − λ ) h(t)=(t-\beta)(t-\lambda) h ( t ) = ( t − β ) ( t − λ ) も σ ( A ) \sigma(A) σ ( A ) 上で非負で、0 ≤ ∫ h d μ ψ = ε 2 + ( η − β ) ( η − λ ) 0\le\int h\,d\mu_\psi=\varepsilon^2+(\eta-\beta)(\eta-\lambda) 0 ≤ ∫ h d μ ψ = ε 2 + ( η − β ) ( η − λ ) から λ ≥ η − ε 2 / ( β − η ) \lambda\ge\eta-\varepsilon^2/(\beta-\eta) λ ≥ η − ε 2 / ( β − η ) 。積分 ∫ t 2 d μ ψ = ∥ A ψ ∥ 2 \int t^2d\mu_\psi=\|A\psi\|^2 ∫ t 2 d μ ψ = ∥ A ψ ∥ 2 が有限なのは ψ ∈ D ( A ) \psi\in D(A) ψ ∈ D ( A ) による。
∎
α β η λ σ(A) σ(A)
Kato–Temple の不等式の配置(模式図)。スペクトルは α 以下、λ の 1 点、β 以上にあり、近似値 η と、λ を含む区間(太線)を示す。
系 7.3 (Temple の不等式 ) ✓
λ 0 = inf σ ( A ) \lambda_0=\inf\sigma(A) λ 0 = inf σ ( A ) が孤立固有値で、σ ( A ) ∖ { λ 0 } ⊂ [ β , ∞ ) \sigma(A)\setminus\{\lambda_0\}\subset[\beta,\infty) σ ( A ) ∖ { λ 0 } ⊂ [ β , ∞ ) 、η < β \eta<\beta η < β とする。このとき
η − ε 2 β − η ≤ λ 0 ≤ η . \eta-\frac{\varepsilon^2}{\beta-\eta}\ \le\ \lambda_0\ \le\ \eta . η − β − η ε 2 ≤ λ 0 ≤ η .
証明(系 7.3) 上界は λ 0 = inf σ ( A ) ≤ ∫ t d μ ψ = η \lambda_0=\inf\sigma(A)\le\int t\,d\mu_\psi=\eta λ 0 = inf σ ( A ) ≤ ∫ t d μ ψ = η による。下界は 定理 7.2 を α < min { λ 0 , η } \alpha<\min\{\lambda_0,\eta\} α < min { λ 0 , η } として適用すればよい(λ 0 ≤ η \lambda_0\le\eta λ 0 ≤ η なので α < λ 0 \alpha<\lambda_0 α < λ 0 と取れば十分)。
∎
例 7.4 (調和振動子の基底エネルギーの下界 )
A = − d 2 / d x 2 + x 2 A=-d^2/dx^2+x^2 A = − d 2 / d x 2 + x 2 の固有値は 1 , 3 , 5 , … 1,3,5,\dots 1 , 3 , 5 , … である。試行関数 ψ a ( x ) = ( a / π ) 1 / 4 e − a x 2 / 2 \psi_a(x)=(a/\pi)^{1/4}e^{-ax^2/2} ψ a ( x ) = ( a / π ) 1/4 e − a x 2 /2 (a > 0 a>0 a > 0 )について、ψ a ′ ′ = ( a 2 x 2 − a ) ψ a \psi_a''=(a^2x^2-a)\psi_a ψ a ′′ = ( a 2 x 2 − a ) ψ a から A ψ a = ( a + ( 1 − a 2 ) x 2 ) ψ a A\psi_a=(a+(1-a^2)x^2)\psi_a A ψ a = ( a + ( 1 − a 2 ) x 2 ) ψ a 。Gauss 積分 ⟨ x 2 ⟩ = 1 / ( 2 a ) \langle x^2\rangle=1/(2a) ⟨ x 2 ⟩ = 1/ ( 2 a ) 、⟨ x 4 ⟩ = 3 / ( 4 a 2 ) \langle x^4\rangle=3/(4a^2) ⟨ x 4 ⟩ = 3/ ( 4 a 2 ) を使うと
η = a 2 + 1 2 a , ∥ A ψ a ∥ 2 = 1 + 3 ( 1 − a 2 ) 2 4 a 2 , ε 2 = ( a 2 − 1 ) 2 2 a 2 . \eta=\frac{a^2+1}{2a},\qquad \|A\psi_a\|^2=1+\frac{3(1-a^2)^2}{4a^2},\qquad \varepsilon^2=\frac{(a^2-1)^2}{2a^2}. η = 2 a a 2 + 1 , ∥ A ψ a ∥ 2 = 1 + 4 a 2 3 ( 1 − a 2 ) 2 , ε 2 = 2 a 2 ( a 2 − 1 ) 2 . a = 2 a=2 a = 2 では η = 5 / 4 \eta=5/4 η = 5/4 、ε 2 = 9 / 8 \varepsilon^2=9/8 ε 2 = 9/8 で、β = 3 \beta=3 β = 3 とした 系 7.3 は 17 / 28 ≤ λ 0 ≤ 5 / 4 17/28\le\lambda_0\le 5/4 17/28 ≤ λ 0 ≤ 5/4 を与える。ところが ψ a \psi_a ψ a は偶関数なので、偶関数の部分空間(A A A を簡約する)で考えれば次の固有値は 5 5 5 で、β = 5 \beta=5 β = 5 として 19 / 20 ≤ λ 0 ≤ 5 / 4 19/20\le\lambda_0\le 5/4 19/20 ≤ λ 0 ≤ 5/4 に改良される。真の値は 1 1 1 である。対称性で隙間を広げることが、Temple の不等式を実用的にする典型的な工夫である。
有限の箱の中で計算した作用素は、全空間の作用素とは別の Hilbert 空間に住んでいます。部分空間の上の作用素で近似する状況を扱うために、強レゾルベント収束(定義 3.1[作用素の収束と摂動] )を拡張します。
定義 8.1 (一般化強レゾルベント収束 ) ↗
H n ⊂ H \mathcal{H}_n\subset\mathcal{H} H n ⊂ H を閉部分空間、P n P_n P n をその上への直交射影、H n H_n H n を H n \mathcal{H}_n H n の上の自己共役作用素、H H H を H \mathcal{H} H の上の自己共役作用素とする。すべての z ∈ C ∖ R z\in\mathbb{C}\setminus\mathbb{R} z ∈ C ∖ R について ( H n − z ) − 1 P n → ( H − z ) − 1 (H_n-z)^{-1}P_n\to(H-z)^{-1} ( H n − z ) − 1 P n → ( H − z ) − 1 が H \mathcal{H} H 上で強収束するとき、H n H_n H n は H H H に一般化強レゾルベント収束 するという。f ( H n ) P n f(H_n)P_n f ( H n ) P n は H n ⊥ \mathcal{H}_n^\perp H n ⊥ 上で 0 0 0 とおいた H \mathcal{H} H 上の作用素と考える。
補題 8.2 (芯による判定 ) ✓
D D D を H H H の芯とし、各 φ ∈ D \varphi\in D φ ∈ D について、ある n 0 n_0 n 0 以降 φ ∈ D ( H n ) \varphi\in D(H_n) φ ∈ D ( H n ) かつ H n φ → H φ H_n\varphi\to H\varphi H n φ → H φ となるとする。このとき H n H_n H n は H H H に一般化強レゾルベント収束し、P n → 1 P_n\to 1 P n → 1 は強収束する。
証明(補題 8.2) z ∈ C ∖ R z\in\mathbb{C}\setminus\mathbb{R} z ∈ C ∖ R を固定し、R n = ( H n − z ) − 1 P n R_n=(H_n-z)^{-1}P_n R n = ( H n − z ) − 1 P n 、R = ( H − z ) − 1 R=(H-z)^{-1} R = ( H − z ) − 1 とおく。∥ R n ∥ ≤ 1 / ∣ I m z ∣ \|R_n\|\le 1/|\mathrm{Im}\,z| ∥ R n ∥ ≤ 1/∣ Im z ∣ 。φ ∈ D \varphi\in D φ ∈ D 、n ≥ n 0 n\ge n_0 n ≥ n 0 なら ( H n − z ) φ ∈ H n (H_n-z)\varphi\in\mathcal{H}_n ( H n − z ) φ ∈ H n なので R n ( H n − z ) φ = φ R_n(H_n-z)\varphi=\varphi R n ( H n − z ) φ = φ 、したがって
R n ( H − z ) φ − R ( H − z ) φ = R n ( H − z ) φ − R n ( H n − z ) φ = R n ( H φ − H n φ ) → 0. R_n(H-z)\varphi-R(H-z)\varphi=R_n(H-z)\varphi-R_n(H_n-z)\varphi=R_n(H\varphi-H_n\varphi)\to 0 . R n ( H − z ) φ − R ( H − z ) φ = R n ( H − z ) φ − R n ( H n − z ) φ = R n ( H φ − H n φ ) → 0. D D D は H H H の芯なので ( H − z ) D (H-z)D ( H − z ) D は H \mathcal{H} H で稠密である(( H − z ) − 1 (H-z)^{-1} ( H − z ) − 1 は有界で、D D D はグラフノルムで D ( H ) D(H) D ( H ) に稠密)。R n R_n R n が一様有界なので、稠密集合上の収束から全体での強収束が従う。また φ ∈ D \varphi\in D φ ∈ D なら n ≥ n 0 n\ge n_0 n ≥ n 0 で P n φ = φ P_n\varphi=\varphi P n φ = φ で、D D D は稠密、∥ P n ∥ ≤ 1 \|P_n\|\le 1 ∥ P n ∥ ≤ 1 なので P n → 1 P_n\to 1 P n → 1 も強収束する。
∎
定理 8.3 (一般化強レゾルベント収束とスペクトル ) ✓
H n H_n H n が H H H に一般化強レゾルベント収束するとする。
f ∈ C 0 ( R ) f\in C_0(\mathbb{R}) f ∈ C 0 ( R ) なら f ( H n ) P n → f ( H ) f(H_n)P_n\to f(H) f ( H n ) P n → f ( H ) は強収束する。
λ ∈ σ ( H ) \lambda\in\sigma(H) λ ∈ σ ( H ) なら、任意の ε > 0 \varepsilon>0 ε > 0 に対して N N N があり、n ≥ N n\ge N n ≥ N で σ ( H n ) ∩ ( λ − ε , λ + ε ) ≠ ∅ \sigma(H_n)\cap(\lambda-\varepsilon,\lambda+\varepsilon)\ne\emptyset σ ( H n ) ∩ ( λ − ε , λ + ε ) = ∅ 。
証明(定理 8.3) (1) A \mathcal{A} A を ( t − i ) − 1 (t-i)^{-1} ( t − i ) − 1 、( t + i ) − 1 (t+i)^{-1} ( t + i ) − 1 が生成する多項式の代数とすると、命題 3.3[作用素の収束と摂動] の証明のとおり A \mathcal{A} A は C 0 ( R ) C_0(\mathbb{R}) C 0 ( R ) で一様ノルムについて稠密である。P n P_n P n は H n H_n H n のレゾルベントと可換で P n 2 = P n P_n^2=P_n P n 2 = P n なので、g ∈ A g\in\mathcal{A} g ∈ A に対し g ( H n ) P n g(H_n)P_n g ( H n ) P n は ( H n ∓ i ) − 1 P n (H_n\mp i)^{-1}P_n ( H n ∓ i ) − 1 P n の積の有限和であり、各因子はノルム 1 1 1 以下で強収束するから g ( H n ) P n → g ( H ) g(H_n)P_n\to g(H) g ( H n ) P n → g ( H ) は強収束する。∥ f ( H n ) P n − g ( H n ) P n ∥ ≤ ∥ f − g ∥ ∞ \|f(H_n)P_n-g(H_n)P_n\|\le\|f-g\|_\infty ∥ f ( H n ) P n − g ( H n ) P n ∥ ≤ ∥ f − g ∥ ∞ と ε / 3 \varepsilon/3 ε /3 論法で一般の f f f に及ぶ。
(2) 否定すると、ε > 0 \varepsilon>0 ε > 0 と部分列で σ ( H n k ) ∩ ( λ − ε , λ + ε ) = ∅ \sigma(H_{n_k})\cap(\lambda-\varepsilon,\lambda+\varepsilon)=\emptyset σ ( H n k ) ∩ ( λ − ε , λ + ε ) = ∅ 。0 ≤ f ≤ 1 0\le f\le 1 0 ≤ f ≤ 1 、台が ( λ − ε , λ + ε ) (\lambda-\varepsilon,\lambda+\varepsilon) ( λ − ε , λ + ε ) に含まれ、I = ( λ − ε / 2 , λ + ε / 2 ) I=(\lambda-\varepsilon/2,\lambda+\varepsilon/2) I = ( λ − ε /2 , λ + ε /2 ) で f = 1 f=1 f = 1 となる連続関数を取ると f ( H n k ) = 0 f(H_{n_k})=0 f ( H n k ) = 0 で、(1) から f ( H ) = 0 f(H)=0 f ( H ) = 0 。一方 f ( H ) E H ( I ) = E H ( I ) ≠ 0 f(H)E_H(I)=E_H(I)\ne 0 f ( H ) E H ( I ) = E H ( I ) = 0 (λ ∈ σ ( H ) \lambda\in\sigma(H) λ ∈ σ ( H ) )で矛盾する。
∎
例 8.4 (箱への制限 )
H = L 2 ( R ) \mathcal{H}=L^2(\mathbb{R}) H = L 2 ( R ) 、H n = L 2 ( − n , n ) \mathcal{H}_n=L^2(-n,n) H n = L 2 ( − n , n ) (外側を 0 0 0 として埋め込む)、H n H_n H n を ( − n , n ) (-n,n) ( − n , n ) 上の Dirichlet ラプラシアン(定義域 H 2 ( − n , n ) ∩ H 0 1 ( − n , n ) H^2(-n,n)\cap H^1_0(-n,n) H 2 ( − n , n ) ∩ H 0 1 ( − n , n ) )、H = − d 2 / d x 2 H=-d^2/dx^2 H = − d 2 / d x 2 (定義域 H 2 ( R ) H^2(\mathbb{R}) H 2 ( R ) )とする。D = C c ∞ ( R ) D=C_c^\infty(\mathbb{R}) D = C c ∞ ( R ) は H H H の芯で、φ ∈ D \varphi\in D φ ∈ D の台が ( − n , n ) (-n,n) ( − n , n ) に入れば φ ∈ D ( H n ) \varphi\in D(H_n) φ ∈ D ( H n ) 、H n φ = − φ ′ ′ = H φ H_n\varphi=-\varphi''=H\varphi H n φ = − φ ′′ = H φ 。補題 8.2 から H n → H H_n\to H H n → H は一般化強レゾルベント収束である。σ ( H n ) = { ( k π / 2 n ) 2 : k ≥ 1 } \sigma(H_n)=\{(k\pi/2n)^2:k\ge 1\} σ ( H n ) = {( k π /2 n ) 2 : k ≥ 1 } は n → ∞ n\to\infty n → ∞ で [ 0 , ∞ ) = σ ( H ) [0,\infty)=\sigma(H) [ 0 , ∞ ) = σ ( H ) を埋めつくし、定理 8.3 (2) と整合している。
演習 9.1 標準
A = − d 2 / d x 2 + x 2 A=-d^2/dx^2+x^2 A = − d 2 / d x 2 + x 2 、B = x B=x B = x とし、E 0 = 1 E_0=1 E 0 = 1 (固有ベクトル h 0 h_0 h 0 )について 定理 6.1 の漸化式で E 1 , E 2 , E 3 E_1,E_2,E_3 E 1 , E 2 , E 3 を計算せよ。また x 2 + κ x = ( x + κ / 2 ) 2 − κ 2 / 4 x^2+\kappa x=(x+\kappa/2)^2-\kappa^2/4 x 2 + κ x = ( x + κ /2 ) 2 − κ 2 /4 を使って厳密な E ( κ ) E(\kappa) E ( κ ) を求め、比較せよ。Hermite 関数について x h n = ( n + 1 ) / 2 h n + 1 + n / 2 h n − 1 xh_n=\sqrt{(n+1)/2}\,h_{n+1}+\sqrt{n/2}\,h_{n-1} x h n = ( n + 1 ) /2 h n + 1 + n /2 h n − 1 、A h n = ( 2 n + 1 ) h n Ah_n=(2n+1)h_n A h n = ( 2 n + 1 ) h n を使ってよい。
解答 S h n = h n / ( 2 n ) Sh_n=h_n/(2n) S h n = h n / ( 2 n ) (n ≥ 1 n\ge1 n ≥ 1 )、S h 0 = 0 Sh_0=0 S h 0 = 0 である。E 1 = ⟨ h 0 , x h 0 ⟩ = 0 E_1=\langle h_0,xh_0\rangle=0 E 1 = ⟨ h 0 , x h 0 ⟩ = 0 。B h 0 = h 1 / 2 Bh_0=h_1/\sqrt2 B h 0 = h 1 / 2 なので χ 1 = − S B h 0 = − h 1 / ( 2 2 ) \chi_1=-SBh_0=-h_1/(2\sqrt2) χ 1 = − S B h 0 = − h 1 / ( 2 2 ) 、E 2 = ⟨ B h 0 , χ 1 ⟩ = − 1 / 4 E_2=\langle Bh_0,\chi_1\rangle=-1/4 E 2 = ⟨ B h 0 , χ 1 ⟩ = − 1/4 。B χ 1 = − ( h 2 + h 0 / 2 ) / ( 2 2 ) B\chi_1=-(h_2+h_0/\sqrt2)/(2\sqrt2) B χ 1 = − ( h 2 + h 0 / 2 ) / ( 2 2 ) なので χ 2 = S ( − B χ 1 + E 1 χ 1 ) = h 2 / ( 8 2 ) \chi_2=S(-B\chi_1+E_1\chi_1)=h_2/(8\sqrt2) χ 2 = S ( − B χ 1 + E 1 χ 1 ) = h 2 / ( 8 2 ) 、E 3 = ⟨ h 1 / 2 , h 2 / ( 8 2 ) ⟩ = 0 E_3=\langle h_1/\sqrt2,h_2/(8\sqrt2)\rangle=0 E 3 = ⟨ h 1 / 2 , h 2 / ( 8 2 )⟩ = 0 。厳密には平行移動 x ↦ x + κ / 2 x\mapsto x+\kappa/2 x ↦ x + κ /2 がユニタリなので E ( κ ) = 1 − κ 2 / 4 E(\kappa)=1-\kappa^2/4 E ( κ ) = 1 − κ 2 /4 で、E 2 = − 1 / 4 E_2=-1/4 E 2 = − 1/4 、他の E n = 0 E_n=0 E n = 0 (n ≠ 0 , 2 n\ne 0,2 n = 0 , 2 )と一致する。B = x B=x B = x は A A A 有界(∥ x ψ ∥ 2 ≤ ⟨ ψ , A ψ ⟩ ≤ ∥ A ψ ∥ ∥ ψ ∥ \|x\psi\|^2\le\langle\psi,A\psi\rangle\le\|A\psi\|\|\psi\| ∥ x ψ ∥ 2 ≤ ⟨ ψ , A ψ ⟩ ≤ ∥ A ψ ∥∥ ψ ∥ から相対限界 0 0 0 )なので定理の仮定を満たす。
演習 9.2 標準
例 7.4 の A A A について、奇関数の試行関数 ψ ( x ) = c x e − x 2 \psi(x)=c\,xe^{-x^2} ψ ( x ) = c x e − x 2 (a = 2 a=2 a = 2 、c c c は正規化定数)を取る。η = 15 / 4 \eta=15/4 η = 15/4 、ε 2 = 27 / 8 \varepsilon^2=27/8 ε 2 = 27/8 を確かめ、(a) 全空間で α = 1 \alpha=1 α = 1 、β = 5 \beta=5 β = 5 とした Kato–Temple の不等式、(b) 奇関数の部分空間での Temple の不等式から、固有値 3 3 3 の評価を求めよ。
解答 a a a を一般にして ψ = c x e − a x 2 / 2 \psi=c\,xe^{-ax^2/2} ψ = c x e − a x 2 /2 とすると ψ ′ ′ = ( a 2 x 3 − 3 a x ) c e − a x 2 / 2 \psi''=(a^2x^3-3ax)ce^{-ax^2/2} ψ ′′ = ( a 2 x 3 − 3 a x ) c e − a x 2 /2 なので A ψ = ( 3 a + ( 1 − a 2 ) x 2 ) ψ A\psi=(3a+(1-a^2)x^2)\psi A ψ = ( 3 a + ( 1 − a 2 ) x 2 ) ψ 。重み x 2 e − a x 2 x^2e^{-ax^2} x 2 e − a x 2 に関する平均で ⟨ x 2 ⟩ = 3 / ( 2 a ) \langle x^2\rangle=3/(2a) ⟨ x 2 ⟩ = 3/ ( 2 a ) 、⟨ x 4 ⟩ = 15 / ( 4 a 2 ) \langle x^4\rangle=15/(4a^2) ⟨ x 4 ⟩ = 15/ ( 4 a 2 ) だから η = 3 a + ( 1 − a 2 ) ⋅ 3 / ( 2 a ) = 3 ( a 2 + 1 ) / ( 2 a ) \eta=3a+(1-a^2)\cdot 3/(2a)=3(a^2+1)/(2a) η = 3 a + ( 1 − a 2 ) ⋅ 3/ ( 2 a ) = 3 ( a 2 + 1 ) / ( 2 a ) 、∥ A ψ ∥ 2 = 9 a 2 + 9 ( 1 − a 2 ) + 15 ( 1 − a 2 ) 2 / ( 4 a 2 ) \|A\psi\|^2=9a^2+9(1-a^2)+15(1-a^2)^2/(4a^2) ∥ A ψ ∥ 2 = 9 a 2 + 9 ( 1 − a 2 ) + 15 ( 1 − a 2 ) 2 / ( 4 a 2 ) 、ε 2 = 3 ( a 2 − 1 ) 2 / ( 2 a 2 ) \varepsilon^2=3(a^2-1)^2/(2a^2) ε 2 = 3 ( a 2 − 1 ) 2 / ( 2 a 2 ) 。a = 2 a=2 a = 2 で η = 15 / 4 \eta=15/4 η = 15/4 、ε 2 = 27 / 8 \varepsilon^2=27/8 ε 2 = 27/8 。(a) 1.05 = 15 / 4 − ( 27 / 8 ) / ( 5 / 4 ) ≤ 3 ≤ 15 / 4 + ( 27 / 8 ) / ( 11 / 4 ) = 15 / 4 + 27 / 22 ≈ 4.98 1.05=15/4-(27/8)/(5/4)\le 3\le 15/4+(27/8)/(11/4)=15/4+27/22\approx 4.98 1.05 = 15/4 − ( 27/8 ) / ( 5/4 ) ≤ 3 ≤ 15/4 + ( 27/8 ) / ( 11/4 ) = 15/4 + 27/22 ≈ 4.98 。(b) 奇関数の部分空間(A A A を簡約する)では最小固有値が 3 3 3 、次が 7 7 7 なので、β = 7 \beta=7 β = 7 で 141 / 52 ≈ 2.71 ≤ 3 ≤ 15 / 4 141/52\approx 2.71\le 3\le 15/4 141/52 ≈ 2.71 ≤ 3 ≤ 15/4 。対称性で隙間を広げると評価が大きく改善する。
演習 9.3 難
V V V を R \mathbb{R} R 上の有界連続関数で lim ∣ x ∣ → ∞ V ( x ) = 0 \lim_{|x|\to\infty}V(x)=0 lim ∣ x ∣ → ∞ V ( x ) = 0 、H = − d 2 / d x 2 + V H=-d^2/dx^2+V H = − d 2 / d x 2 + V が負の固有値 λ 1 < 0 \lambda_1<0 λ 1 < 0 を最小固有値にもつとする。H n H_n H n を ( − n , n ) (-n,n) ( − n , n ) 上の Dirichlet 境界条件の − d 2 / d x 2 + V -d^2/dx^2+V − d 2 / d x 2 + V とする。λ 1 ( H n ) = inf σ ( H n ) \lambda_1(H_n)=\inf\sigma(H_n) λ 1 ( H n ) = inf σ ( H n ) について λ 1 ( H n ) ≥ λ 1 \lambda_1(H_n)\ge\lambda_1 λ 1 ( H n ) ≥ λ 1 かつ λ 1 ( H n ) → λ 1 \lambda_1(H_n)\to\lambda_1 λ 1 ( H n ) → λ 1 を示せ。(H H H の形式定義域は H 1 ( R ) H^1(\mathbb{R}) H 1 ( R ) 、H n H_n H n のそれは H 0 1 ( − n , n ) H^1_0(-n,n) H 0 1 ( − n , n ) で、後者は外側を 0 0 0 として前者に含まれることを使ってよい。)
解答 下に有界な自己共役作用素 T T T とその形式 q q q について inf σ ( T ) = inf { q [ ψ ] : ψ ∈ Q ( q ) , ∥ ψ ∥ = 1 } \inf\sigma(T)=\inf\{q[\psi]:\psi\in Q(q),\ \|\psi\|=1\} inf σ ( T ) = inf { q [ ψ ] : ψ ∈ Q ( q ) , ∥ ψ ∥ = 1 } である。実際、表現定理(定理 4.2[自己共役性の判定] (2))から q [ ψ ] = ∫ t d μ ψ ( t ) ≥ inf σ ( T ) ∥ ψ ∥ 2 q[\psi]=\int t\,d\mu_\psi(t)\ge\inf\sigma(T)\|\psi\|^2 q [ ψ ] = ∫ t d μ ψ ( t ) ≥ inf σ ( T ) ∥ ψ ∥ 2 で、逆向きは D ( T ) ⊂ Q ( q ) D(T)\subset Q(q) D ( T ) ⊂ Q ( q ) とミニマックス原理 定理 6.3[作用素の収束と摂動] の μ 1 ( T ) = inf σ ( T ) \mu_1(T)=\inf\sigma(T) μ 1 ( T ) = inf σ ( T ) による。H 0 1 ( − n , n ) ⊂ H 1 ( R ) H^1_0(-n,n)\subset H^1(\mathbb{R}) H 0 1 ( − n , n ) ⊂ H 1 ( R ) で、そこでは二つの形式が一致するので、下限の比較から λ 1 ( H n ) ≥ λ 1 \lambda_1(H_n)\ge\lambda_1 λ 1 ( H n ) ≥ λ 1 。また H 0 1 ( − n , n ) ⊂ H 0 1 ( − n − 1 , n + 1 ) H^1_0(-n,n)\subset H^1_0(-n-1,n+1) H 0 1 ( − n , n ) ⊂ H 0 1 ( − n − 1 , n + 1 ) から λ 1 ( H n ) \lambda_1(H_n) λ 1 ( H n ) は非増加である。逆向きには、H H H の正規化固有関数 u u u と ∥ u − φ ∥ H 1 \|u-\varphi\|_{H^1} ∥ u − φ ∥ H 1 が小さい φ ∈ C c ∞ \varphi\in C_c^\infty φ ∈ C c ∞ を取れば、形式が H 1 H^1 H 1 ノルムについて連続なので q [ φ ] / ∥ φ ∥ 2 q[\varphi]/\|\varphi\|^2 q [ φ ] /∥ φ ∥ 2 は λ 1 \lambda_1 λ 1 に近く、φ \varphi φ の台が ( − n , n ) (-n,n) ( − n , n ) に入る n n n では λ 1 ( H n ) ≤ q [ φ ] / ∥ φ ∥ 2 \lambda_1(H_n)\le q[\varphi]/\|\varphi\|^2 λ 1 ( H n ) ≤ q [ φ ] /∥ φ ∥ 2 。よって λ 1 ( H n ) → λ 1 \lambda_1(H_n)\to\lambda_1 λ 1 ( H n ) → λ 1 。(定理 8.3 だけからは、λ 1 \lambda_1 λ 1 の近くに σ ( H n ) \sigma(H_n) σ ( H n ) の点があることしか分からない。最小値であることの比較が、ミニマックス原理から来る部分である。)
T. Kato, Perturbation Theory for Linear Operators , 2nd ed., Springer, 1976 — 第 II 章(有限次元の摂動論:Puiseux 級数、変換関数、Rellich の定理)、第 VII 章(解析族)。
F. Rellich, Perturbation Theory of Eigenvalue Problems , Gordon and Breach, 1969.
M. Reed, B. Simon, Methods of Modern Mathematical Physics IV: Analysis of Operators , Academic Press, 1978 — 第 XII 章(解析摂動、Rayleigh–Schrödinger 級数)。
T. Kato, “On the upper and lower bounds of eigenvalues”, Journal of the Physical Society of Japan 4 (1949), 334–339.
G. Temple, “The theory of Rayleigh’s principle as applied to continuous systems”, Proceedings of the Royal Society of London A 119 (1928), 276–293.
J. Weidmann, Lineare Operatoren in Hilberträumen, Teil I: Grundlagen , Teubner, 2000 — レゾルベント収束(一般化強レゾルベント収束を含む)の章。
執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 2 回(Claude Code)・最終 2026-10-12