近似の数理 第3章 §3.3 ★★★ 原子・分子のハミルトニアンが自己共役であることの証明。クーロン特異性が「弱い」ことの定量化
ハミルトニアン (3.1) は、自己共役な H0=∑i(−21Δi) に、原点(原子核・電子の衝突点)で発散するクーロンポテンシャル V を足したものである。V は有界作用素でないので、和が自己共役である保証は自明でない。鍵は「V の特異性は運動エネルギーで抑え込める」という定量的な不等式で、その出発点が Hardy の不等式である。
Hardy の不等式
定理 3.12(Hardy の不等式(3次元) ★★★)
ψ∈H1(R3) について
∫R3∣x∣2∣ψ(x)∣2dx ≤ 4∫R3∣∇ψ(x)∣2dx.(3.5)
定数 4 はこれ以上小さくできない。
証明
滑らかで台がコンパクトな ψ について示せば、H1 の稠密性と両辺の連続性から一般の場合が従う。ベクトル場 F(x)=x/∣x∣2 を考える。r=∣x∣ として∇⋅F=∇⋅(x/r2)=3/r2−2x⋅x/r4=1/r2(3次元で一般に (d−2)/r2)。そこで
∫r2∣ψ∣2dx=∫∣ψ∣2∇⋅Fdx=−∫F⋅∇∣ψ∣2dx=−2Re∫r2x⋅ψˉ∇ψdx.
部分積分の境界項は、原点を半径 ε の球でくり抜いて評価すると ∣ψ∣2∣F∣×(面積)∼ε−1⋅ε2→0 である。右辺を Cauchy–Schwarz の不等式で評価すると、∣x/r2∣=1/r より
∫r2∣ψ∣2 ≤ 2(∫r2∣ψ∣2)1/2(∫∣∇ψ∣2)1/2.
両辺を (∫∣ψ∣2/r2)1/2 で割って2乗すれば (3.5) を得る。定数が最良であることは、r−1/2 に近い関数で比が 1 に近づくこと(例題 3.2 と表 3.1)で分かる。
Hardy の不等式が言っていること
∣ψ∣2/r2 の積分が大きいのは、ψ が原点の近くに集中しているときである。ところが原点の近くの狭い領域に閉じ込めると、不確定性によって勾配が大きくなる。Hardy の不等式は「原点への集中の度合い」が「運動エネルギー」で上から抑えられることを、定数まで込めて述べた不確定性原理である。ΔxΔp≥21 という分散の不等式とは違い、平均からのずれではなく原点からの距離そのものを測っている点が、ポテンシャルの評価に使える理由である。
相対有界性:クーロンポテンシャルの相対限界は0
定義 3.13(相対有界)
自己共役作用素 A と作用素 B について、D(A)⊂D(B) で、ある a,b≥0 について
∥Bψ∥≤a∥Aψ∥+b∥ψ∥(ψ∈D(A))
が成り立つとき、B は A に相対有界であるという。このような a の下限を相対限界という。
相対限界は「B の大きさのうち、A に比例する部分の係数」である。b はいくら大きくてもよく、問題は a のほうである。a が小さいことは、高エネルギー(∥Aψ∥ が大きい)の状態で B が A に比べて無視できることを意味する。
命題 3.14(クーロンポテンシャルの相対限界は0 ★★★)
L2(R3) で V(x)=Z/∣x∣ とする。任意の a>0 について、全ての ψ∈H2(R3) で
∥Vψ∥≤a∥Δψ∥+aZ2∥ψ∥.(3.6)
証明
Hardy の不等式より ∥Vψ∥=Z∥ψ/r∥≤2Z∥∇ψ∥。部分積分と Cauchy–Schwarz により∥∇ψ∥2=⟨ψ,−Δψ⟩≤∥ψ∥∥Δψ∥。相加相乗平均 uv≤ηu+v/(4η)(η>0)をu=∥Δψ∥、v=∥ψ∥ に使うと ∥∇ψ∥≤η∥Δψ∥+∥ψ∥/(4η)。合わせて∥Vψ∥≤2Zη∥Δψ∥+2ηZ∥ψ∥。a=2Zη とおけば Z/(2η)=Z2/a となる(検算スクリプト)。
a をいくら小さくしても成り立つ代わりに、b=Z2/a が大きくなる。物理的には、クーロンポテンシャルは「低エネルギーでは強いが、高エネルギー(短い波長)では運動エネルギーに比べて取るに足らない」ということである。1/r は長さの −1 乗、運動エネルギーは −2 乗なので、短距離では運動エネルギーが必ず勝つ。次元の数え方としてはこれだけのことだが、それを定数つきの不等式にしたのが (3.6) である。
多電子系への拡張は次のとおりである。L2(R3N) で xi についての項 Z/∣xi−R∣ は、他の変数を止めて xi について (3.6) を使い、残りの変数で積分すれば∥(Z/∣xi−R∣)ψ∥≤a∥Δiψ∥+(Z2/a)∥ψ∥ を得る。フーリエ変換すると∣ki∣2≤∑j∣kj∣2 なので ∥Δiψ∥≤∥Δψ∥(Δ=∑jΔj)。電子間反発 1/∣xi−xj∣ は、相対座標 u=xi−xj と重心座標 w=(xi+xj)/2 に変数変換するとu についての1体のクーロンポテンシャルになり、∇u=(∇i−∇j)/2 のフーリエ表示∣(ki−kj)/2∣2≤(∣ki∣2+∣kj∣2)/2 から同じ形の評価が出る。有限個の項の和なので、(3.1) のポテンシャル部分全体 V も H0 に対して相対限界0をもつ。
Kato–Rellich の定理と加藤の結果
定理 3.15(Kato–Rellich ★★★)
A を自己共役、B を対称作用素とし、D(A)⊂D(B) かつある a<1, b≥0 で∥Bψ∥≤a∥Aψ∥+b∥ψ∥(ψ∈D(A))とする。このとき A+B は D(A) 上で自己共役である。さらに A が下に有界なら A+B も下に有界である。
証明
A+B が D(A) 上で対称であることは明らか。判定法(定理 3.8、±i を ±iμ に置き換えても同じ証明が通る)により、μ>0 を大きくとって Ran(A+B∓iμ)=H を示せばよい。A は自己共役なので A−iμ は全単射で、(3.4) より ∥A(A−iμ)−1∥≤1、∥(A−iμ)−1∥≤1/μ。よって任意の ϕ について
∥B(A−iμ)−1ϕ∥≤a∥A(A−iμ)−1ϕ∥+b∥(A−iμ)−1ϕ∥≤(a+μb)∥ϕ∥.
μ>b/(1−a) なら K=B(A−iμ)−1 は ∥K∥<1 を満たし、1+K は Neumann 級数 ∑n(−K)n で逆をもつ。
A+B−iμ=(1+K)(A−iμ)
の右辺は二つの全射の合成なので全射である。+iμ も同様。下に有界性は第 3.4 節の後の注意で扱う。
定理 3.16(加藤敏夫(1951) ★★★)
原子・分子のハミルトニアン (3.1) は D(H)=H2(R3N) 上で自己共役であり、下に有界である。反対称な関数の部分空間に制限しても同じことが成り立つ。
証明は、H0 が H2 上で自己共役(フーリエ変換で実数値関数 ∑i∣ki∣2/2 の掛け算になる)であることと、命題 3.14 の多電子版で a<1 を選べることを、Kato–Rellich の定理に入れるだけである。反対称化は H と交換するので部分空間に制限しても自己共役性は保たれる。この結果が、量子化学の全ての計算に「計算している対象が数学的に存在する」という保証を与えている。しかも定義域は自由粒子と同じ H2 であり、クーロン特異性は波動関数の正則性を2階微分の2乗可積分性の水準では落とさない。波動関数が原子核の位置で尖る(カスプをもつ)ことは、この H2 の範囲に収まる現象である。カスプは第4章でガウス基底の収束が遅くなる理由として再び現れる。
例題 3.2 Hardy の不等式の具体的な関数での確認 ★★
オリジナル[標準的な題材を自作]
R3 の球対称な関数 ψ(x)=f(r) について、比Q[f]=4∫∣∇ψ∣2dx∫∣ψ∣2/r2dx を考える(Hardy の不等式は Q≤1)。
- (1) f=e−r のとき Q を求めよ。
- (2) f=e−r2 のとき Q を求めよ。
- (3) Q が 1 に近づくのは f がどのような形のときか。不等式の証明の等号条件から答えよ。
解答
球対称なので ∫∣ψ∣2/r2dx=4π∫0∞f2dr、∫∣∇ψ∣2dx=4π∫0∞f′2r2dr。
(1) ∫0∞e−2rdr=21、∫0∞e−2rr2dr=232=41。よって Q=4⋅1/41/2=21。
(2) ∫0∞e−2r2dr=21π/2。f′=−2re−r2 より∫0∞4r4e−2r2dr=4⋅8⋅223π/2=83π/2(∫0∞r4e−αr2dr=8α23π/α で α=2)。よってQ=4⋅83π/221π/2=31。
(3) 証明で使った不等式は Cauchy–Schwarz だけなので、等号は ∇ψ と (x/r2)ψ が比例するとき、すなわち f′=−rcf のときである。積分すると f=r−c で、Cauchy–Schwarz の両辺を一致させると c=1/2 を得る。f=r−1/2 は両辺の積分がどちらも対数発散して不等式の外にあるが、これを原点近傍と遠方で切り落とした関数、あるいは f=(1+r2)−s で s↓41 とした列で Q→1 になる(表 3.1。検算スクリプト)。
表 3.1 f=(1+r2)−s での Hardy の比 Q(数値積分)。s↓1/4 で最良定数 Q=1 に近づく。| s | 1 | 0.75 | 0.5 | 0.4 | 0.3 | 0.27 | 0.26 |
| Q | 0.500 | 0.556 | 0.667 | 0.750 | 0.889 | 0.951 | 0.974 |
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com