ハミルトニアンの選び方で固体を考える 第2部 第6章 §2.6.5 ★★ 状態密度を滑らかな関数として扱うために、結晶の大きさを無限にし、準位の幅を 0 とし、次元とホッピングの範囲を限っている。特異点はこれらで丸まる
削ったもの
この章の状態密度は、4つのものを削って得たものである。
1つめは結晶の有限の大きさである。N が有限なら準位は飛び飛びで、式 (2.6.6) の D(ε) はデルタ関数の和である。滑らかな D(ε) は N→∞ の極限でしか意味を持たず、しかもその極限は点ごとの値の極限ではない。極限が意味を持つのは積分した状態密度 N(ε) の側で、周期境界条件のもとでは固有値が格子点の上の εn(k) そのものなので、N(ε) は Brillouin 域の積分の Riemann 和になり、等エネルギー面の測度が 0 である限り収束する。次の例題で見るように、有限の N では van Hove 特異点の所にデルタ関数の大きな重みが集まることがあるが、その重みは N→∞ で 0 になる。
2つめは準位の幅である。この章の εn(k) は、寿命が無限に長い1電子の準位として扱った。実在の電子は、不純物・フォノン・他の電子に散乱されて有限の寿命 τ を持ち、準位はエネルギーの幅 ℏ/τ 程度にぼける。状態密度はこの幅で平均されるので、対数発散は有限の山(高さはおよそ 2π2t1lnℏ/τt)に、x の折れは丸い肩になる。不規則さが強くなると、帯の端にも状態が染み出す。
3つめは次元である。層状の結晶を2次元の正方格子として扱うと、層の間のホッピング t⊥ を削ったことになる。t⊥ を戻すとε=−2t(coskxa+coskya)−2t⊥coskzc(c は層の間隔)で、状態密度は式 (2.6.17) と同じ形の平均Dsite(ε)=π1∫0πdθDsitesq(ε+2t⊥cosθ) になる。t⊥≪t なら ε=εa で式 (2.6.15) を入れて
Dsite(εa)≈π1∫0π2π2tdθln2t⊥∣cosθ∣16t=2π2t1lnt⊥16t(2.6.18)
を得る(π1∫0πln∣cosθ∣dθ=−ln2 を使った)。2次元の対数発散は、t⊥ によって ln(t/t⊥) の高さで頭打ちになり、εa±2t⊥ の2つの3次元の折れに分かれる。
4つめはホッピングの範囲と軌道の数である。この章の数値例は、1サイトに軌道1個、最近接のホッピングだけを残した。例題 2.6.4 で見たように、第2近接を入れると特異点の位置と係数は変わる。ただし特異点の種類と最小個数は定理 2.6.6 で決まっていて、ホッピングの細部によらない。状態密度の形のうち、ホッピングを削っても残るのはこの位相的な部分であり、削ると変わるのは特異点の位置と係数である。
近似の裏付け
型 B(極限) | 小さいもの 1/N→0(単位胞の数 N が無限大。状態密度は積分した N(ε) の意味で比べる) | 証明 証明済み(周期境界条件の tight-binding では固有値が格子点上の εn(k) そのもので、N(ε) は Riemann 和として収束する。等エネルギー面の測度が 0 であることが前提。点ごとの D(ε) は収束しない。連続の周期 Schrödinger 作用素でも、積分した状態密度 N(ε) は境界条件によらず同じ極限になる(第V巻 第12章)) | 行き先 近似の数理 第12章
近似の裏付け
型 A(射影) | 小さいもの 各格子点に局在した原子軌道1個が張る部分空間(最近接のホッピングだけを残す) | 証明 条件付き(原子間距離が大きい深い井戸の極限で、トンネル積分の形まで証明されている〔Simon 1984、Helffer–Sjöstrand 1984〕。第2近接以降を落とす誤差の実在の結晶での大きさは要確認) | 行き先 近似の数理 第4章・近似の数理 第12章
例題 2.6.5 有限の正方格子で状態密度を数える ★★
オリジナルN×N 個の格子点からなる正方格子(N は偶数、周期境界条件、εa=0)の準位 ε=−2t(cosN2πm1+cosN2πm2)(m1,m2=0,…,N−1)を考える。
- (1) ちょうど ε=0 になる準位の数が 2N−2 であることを示せ。N=100 ではそれは全体の何割か。
- (2) N=100 で ∣ε∣≤0.1t の準位の割合を、式 (2.6.16) を積分した値 ∫−0.1t0.1tDsitedε と比べよ。
- (3) (1)(2) から、有限の N の状態密度と N→∞ の状態密度がどういう意味で近いかを述べよ。
解答
(1) 例題 2.6.2 により ε=0 は ∣kx∣+∣ky∣=π/a の菱形の上である。mi を −N/2<mi≤N/2 の範囲に取り直すと、条件は ∣m1∣+∣m2∣=N/2 になる。菱形の4辺の上の格子点は 4⋅2N=2N 個だが、頂点 (±N/2,0) は m1=N/2 と m1=−N/2 が同じ点(N だけずれた点)なので2個が1個に、(0,±N/2) も同様に1個になる。したがって 2N−2 個である。N=100 では 198 個で、全体 104 個の 1.98% である。割合はおよそ 2/N で、N→∞ で 0 になる。
(2) 数値で数えると、N=100 の ∣ε∣≤0.1t の準位の割合は 0.0590 で、このうち 0.0198 が (1) の ε=0 ちょうどの準位である。式 (2.6.16) を数値積分すると ∫−0.1t0.1tDsitedε=0.0616 である。幅 0.2t の窓で数えた割合は、4% 程度の差で積分の値に近い。
(3) 有限の N の状態密度はデルタ関数の和で、ε=0 には 2/N 程度の重みを持つデルタ関数がある。N→∞ の状態密度は ε=0 で対数発散するが、有限の幅の窓で積分した量(窓の中の準位の割合)どうしで比べれば、両者は近い。状態密度の N→∞ の極限は、点ごとの値でなく、積分した量 N(ε) の意味で取るものである。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com