固体物理学では、固体の電子バンド構造(または単にバンド構造)は、電子がその中で持つことができるエネルギーレベルの範囲と、電子が持つことができないエネルギーの範囲(バンドギャップまたは禁制帯と呼ばれる)を表します。
バンド理論は、原子または分子の大きな周期格子中の電子に許容される量子力学的波動関数を調べることで、これらのバンドとバンドギャップを導出します。バンド理論は、電気抵抗や光吸収など、固体の多くの物理的特性を説明するために効果的に用いられており、あらゆる固体デバイス(トランジスタ、太陽電池など) の理解の基礎となっています。

電子バンドとバンドギャップの形成は、固体中の電子に関する2つの相補的なモデルで説明できる。[1] : 161 1つ目は、電子が物質内をほぼ自由に運動すると仮定する、ほぼ自由電子モデルである。このモデルでは、電子状態は自由電子平面波に似ており、結晶格子による摂動はわずかである。このモデルは電子分散関係の起源を説明するが、バンドギャップの説明は微妙である。[2] : 121
2つ目のモデルは、電子が個々の原子に強く束縛されている逆の極限から始まります。孤立した単一の原子の電子は、離散的なエネルギー準位を持つ原子軌道を占有します。2つの原子が十分に接近し、原子軌道が重なり合うと、電子は原子間をトンネルすることができます。このトンネル効果により、原子軌道は異なるエネルギーを持つ分子軌道に分裂(混成)します。 [2] : 117–122
同様に、多数のN個の同一原子が集まって結晶格子などの固体を形成する場合、原子の原子軌道は近くの軌道と重なり合う。[3]各離散エネルギー準位はN個の準位に分裂し、それぞれが異なるエネルギーを持つ。マクロな固体中の原子数は非常に大きい(N ≈ 10 22)ため、互いに混成する軌道の数も非常に多い。このため、隣接する準位はエネルギー的に非常に接近している(10 −22 eV)、[4] [5] [6]であり、連続体、エネルギーバンドを形成すると考えられる。
このバンド形成は、主に原子の最外殻電子(価電子)の特徴であり、化学結合と電気伝導に関与しています。内殻電子軌道は大きく重なり合わないため、そのバンドは非常に狭くなっています。
バンドギャップとは、エネルギーバンドの幅が有限であるために、どのバンドにもカバーされない残りのエネルギー領域です。バンドの幅はそれぞれ異なり、その幅はバンドの元となる原子軌道の重なり具合によって決まります。隣接する2つのバンドが、エネルギー領域を完全にカバーするには幅が足りない場合があります。例えば、コア軌道(1s電子など)に関連するバンドは、隣接する原子間の重なり合いが小さいため、非常に狭くなっています。その結果、コアバンド間には大きなバンドギャップが生じる傾向があります。高次のバンドは比較的大きな軌道と重なり合いが大きく、高エネルギーでは徐々に幅が広くなるため、高エネルギーではバンドギャップは存在しません。
バンド理論は固体の量子状態に対する近似に過ぎず、多数の同一原子または分子が結合した固体に適用される。バンド理論が妥当であるために必要な仮定は以下の通りである。
上記の仮定は多くの重要な実際的な状況で破れており、バンド構造を使用する場合にはバンド理論の限界を注意深く確認する必要があります。


バンド構造の計算では、結晶格子の周期的性質を利用し、その対称性を活用します。格子周期ポテンシャル内の電子について単電子シュレーディンガー方程式を解くと、解として ブロッホ電子が得られます。 ここで、 kは波数ベクトルと呼ばれます。kの各値に対して、シュレーディンガー方程式には複数の解があり、バンドインデックスnでラベル付けされます。n はエネルギーバンドに番号を付ける単純なものです。これらのエネルギー準位はそれぞれkの変化に伴って滑らかに変化し、滑らかな状態バンドを形成します。各バンドに対して、そのバンド内の電子の分散関係を表す関数E n ( k )を定義できます。
ブリルアンゾーン(結晶格子と相関する波動ベクトル空間(逆格子空間)内の多面体)内では、波動ベクトルは任意の値をとります。ブリルアンゾーン外の波動ベクトルは、ブリルアンゾーン内の状態と物理的に同一の状態に対応します。ブリルアンゾーン内の特別な高対称性点/線には、Γ、Δ、Λ、Σなどのラベルが付けられます(図1参照)。
バンドの形状を波数ベクトルの関数として視覚化するのは困難で、4次元空間でE対k x、k y、k zのプロットが必要になるからです。科学文献では、対称点 (多くの場合、それぞれ Δ、Λ、Σ、または[100]、[111]、[110] とラベル付けされている) を結ぶ直線に沿ってkの値に対するE n ( k )の値を示すバンド構造のプロットがよく見られます。[7] [8]バンド構造を視覚化する別の方法は、波数ベクトル空間で一定エネルギー等値面をプロットし、特定の値に等しいエネルギーを持つすべての状態を表示することです。フェルミ準位に等しいエネルギーを持つ状態の等値面は、フェルミ面として知られています。
エネルギーバンドギャップは、バンドギャップを囲む状態の波数ベクトルを使用して分類できます。
電子バンド構造は通常、結晶性物質と関連付けられますが、準結晶や非晶質固体もバンドギャップを示すことがあります。これらは結晶のような単純な対称性を欠き、正確な分散関係を決定することが通常不可能であるため、理論的に研究するのはやや困難です。その結果、固体の電子バンド構造に関する既存の理論的研究は、ほぼすべて結晶性物質に焦点を当ててきました。
状態密度関数g ( E )は、 E付近の電子エネルギーにおける、単位体積あたり、単位エネルギーあたりの電子状態数として定義されます
状態密度関数は、バンド理論に基づく効果の計算において重要です。フェルミの黄金律(光吸収率の計算)においては、励起可能な電子の数と電子の終状態の数の両方を与えます。また、電気伝導率の計算では移動可能な状態の数を与え、電子散乱率の計算では散乱後の終状態の数を与えます。[要出典]
バンドギャップ内のエネルギーの場合、g ( E )=0です。

熱力学的平衡において、エネルギーEの状態が電子で満たされている尤度は、パウリの排他原理を考慮した熱力学的分布であるフェルミ–ディラック分布によって与えられる。 ここで、
物質中の電子密度は、フェルミ・ディラック分布と状態密度の積分で表されます。
バンドの数は無限であり、したがって状態の数も無限であるにもかかわらず、これらのバンドに配置できる電子の数は有限である。電子数の好ましい値は静電気力学の結果である。物質の表面は帯電していても、物質の内部は電荷が中性であることを好む。電荷中性の条件とは、 N / Vが物質内の陽子密度と一致しなければならないことを意味する。これを実現するために、物質は静電気力によってバンド構造をエネルギー的に上下にシフトさせ(それによってg ( E )をシフトさせ)、フェルミ準位に関して正しい平衡状態に達するまで調整する。
固体には、原子が無限のエネルギー準位を持つのと同様に、無限の数のバンドが許容されます。しかし、ほとんどのバンドはエネルギーが高すぎるため、通常は無視されます。[9]逆に、コア軌道に関連する非常に低いエネルギーバンド( 1s電子 など)も存在します。これらの低エネルギーのコアバンドも、常に電子で満たされているため不活性であり、通常は無視されます。[10] 同様に、物質はバンド構造全体に複数のバンドギャップを持ちます。
エレクトロニクスやオプトエレクトロニクスに関係する最も重要なバンドとバンドギャップは、フェルミ準位付近のエネルギーを持つものです。フェルミ準位付近のバンドとバンドギャップには、材料に応じて特別な名前が付けられています。
仮説は、回折の力学理論で一般的に扱われているブロッホの定理を用いた、周期結晶格子における電子波の特殊なケースです。すべての結晶はブラヴェ格子によって特徴付けられる周期構造であり、それぞれのブラヴェ格子に対して、3つの逆格子ベクトル( b⁻⁻、b⁻⁻、b⁻⁻ )の集合に周期性を包含する逆格子を決定できます。ここで、直接格子と同じ周期性を持つ任意の周期ポテンシャルV ( r )は、逆格子ベクトルに関連する成分のみがゼロでない フーリエ級数として展開できます。したがって、展開は次のように表すことができます。 ここで、任意の整数集合( m⁻⁻ 、m⁻⁻ 、m⁻⁻ )に対して、K = m⁻⁻ b⁻ + m⁻ b⁻ + m⁻ b⁻です
この理論から、特定の物質のバンド構造を予測する試みは可能ですが、電子構造計算のためのほとんどの ab initio 法では、観測されるバンドギャップを予測できません。
ほぼ自由電子近似では、電子間の相互作用は完全に無視されます。この近似により、周期ポテンシャル内の電子は、隣接する逆格子ベクトル間の一定の位相シフトを除けば、波動関数とエネルギーが波動ベクトルに関して周期的であるというブロッホの定理が適用できます。周期性の帰結は、ブロッホの定理によって数学的に記述されます。ブロッホの定理によれば、固有状態の波動関数は、 結晶格子上で ブロッホ関数が周期的であるような形、すなわち、
ここで指数nはn番目のエネルギーバンドを指し、波数ベクトルkは電子の運動方向に関連し、rは結晶内の位置、Rは原子サイトの位置である。[12] :179
NFEモデルは、隣接する原子間の距離が狭い金属などの材料で特に有効です。このような材料では、原子軌道と隣接する原子上のポテンシャルの重なりが比較的大きくなります。この場合、電子の波動関数は(修正された)平面波で近似できます。アルミニウムのような金属のバンド構造は、空格子近似に非常に近くなります。
ほぼ自由電子の近似の反対の極端な考え方では、結晶内の電子が構成原子の集合体のように振舞うと仮定する。この強結合モデルでは、時間に依存しない単一電子シュレーディンガー方程式の解が 原子軌道の線形結合で十分に近似されると仮定する。[12] : 245–248 ここで、係数は、この形式の最善の近似解を与えるように選択される。添え字n は原子エネルギーレベルを指し、R は原子サイトを指す。この考え方を使用するより正確なアプローチは、次のように定義されるワニエ関数を使用する。 [12] : 式 42 p. 267 [13] では、 はブロッホの定理の周期部分であり、積分はブリルアンゾーン上で行われる。ここで添え字nは結晶内のn番目のエネルギーバンドを指す。ワニエ関数は原子軌道のように原子サイトの近くに局在するが、ブロッホ関数で定義されているため、結晶ポテンシャルに基づく解と正確に関連している。異なる原子サイトR上のワニエ関数は直交する。ワニエ関数は、 n番目のエネルギーバンド のシュレーディンガー解を次のように表すために用いられる。
TB モデルは、原子軌道と隣接原子上のポテンシャルの重なりが限られている材料でうまく機能します。たとえば、 Si、GaAs、 SiO 2、ダイヤモンドなどの材料のバンド構造は、原子 sp 3軌道に基づいた TB ハミルトニアンでうまく説明できます。遷移金属では、幅広い NFE伝導帯と狭い埋め込み TB d バンドを説明するために、 TB-NFE 混合モデルが使用されます。ワニエ関数の原子軌道部分のラジアル関数は、擬ポテンシャル法を使用することで最も簡単に計算されます。NFE、TB、または組み合わせた NFE-TB バンド構造の計算は、[14] 擬ポテンシャル法に基づく波動関数近似で拡張されることもあり、以降の計算の経済的な出発点としてよく使用されます。
KKR法は、「多重散乱理論」またはグリーン関数法とも呼ばれ、ハミルトニアンではなく逆遷移行列Tの定常値を求める。変分法による実装はコリンガ、コーン、ロストカーによって提案され、コリンガ・コーン・ロストカー法と呼ばれることが多い。[15] [16] KKR法またはグリーン関数法の最も重要な特徴は、(1)問題の2つの側面、すなわち構造(原子の位置)と散乱(原子の化学的性質)を分離すること、(2)グリーン関数が合金やその他の無秩序系に適応可能な電子特性の局所的記述への自然なアプローチを提供することである。この近似の最も単純な形は、重なり合わない球(マフィン型と呼ばれる)を原子の位置を中心に配置する。これらの領域内では、電子が受けるポテンシャルは、与えられた原子核を中心に球対称となるように近似される。残りの格子間領域では、遮蔽されたポテンシャルは定数として近似されます。原子中心球と格子間領域間のポテンシャルの連続性が確保されます。
近年の物理学文献では、電子構造とバンドプロットの大部分は密度汎関数理論(DFT)を用いて計算されています。これはモデルではなく理論、すなわち凝縮物質物理学の微視的な第一原理理論であり、電子密度の汎関数に交換相関項を導入することで電子-電子多体問題に対処しようとします。DFTで計算されたバンドは、多くの場合、角度分解光電子分光法(ARPES)などによって実験的に測定されたバンドと一致することが分かっています。特に、バンド形状は通常、DFTによってよく再現されます。しかし、実験結果と比較すると、DFTバンドには系統的な誤差もあります。特に、DFTは絶縁体と半導体のバンドギャップを約30~40%系統的に過小評価するようです。[17]
DFT は系の基底状態の特性(全エネルギー、原子構造など)のみを予測する理論であり、励起状態の特性は DFT では決定できないと一般的に信じられています。これは誤解です。原理的には、基底状態密度をその特性にマッピングする関数が与えられれば、DFT は系のあらゆる特性(基底状態または励起状態)を決定できます。これがホーエンベルク・コーンの定理の本質です。[18]ただし、実際には、基底状態密度を物質内の電子の励起エネルギーにマッピングする既知の関数は存在しません。したがって、文献で DFT バンド プロットとして引用されているものは、DFTコーン・シャム エネルギー、つまり物理的な解釈がまったくない架空の非相互作用系、コーン・シャム 系のエネルギーの表現です。コーン=シャムの電子構造を、システムの真の準粒子電子構造と混同してはならない。また、コーン=シャムのエネルギーにはクープマンスの定理は成立しない。これは、準粒子エネルギーの近似として真に考えられるハートリー=フォックのエネルギーには成立する。したがって、原理的には、コーン=シャムに基づくDFTはバンド理論ではなく、すなわち、バンドやバンドプロットの計算に適した理論ではない。原理的には、時間依存DFTを用いて真のバンド構造を計算できるが、実際にはこれはしばしば困難である。一般的なアプローチは、ハートリー=フォックの正確な交換反応の一部を組み込んだハイブリッド汎関数の使用である。これは半導体のバンドギャップ予測を大幅に改善するが、金属やワイドバンドギャップ材料に対しては信頼性が低い。[19]
電子間相互作用の多体効果を含むバンドを計算するには、いわゆるグリーン関数法に頼ることができます。実際、システムのグリーン関数がわかれば、システムの基底状態(全エネルギー)と励起状態の観測量の両方が得られます。グリーン関数の極は準粒子エネルギー、つまり固体のバンドです。システムの自己エネルギーがわかれば、ダイソン方程式を解いてグリーン関数を計算できます。固体などの実際のシステムでは、自己エネルギーは非常に複雑な量であり、問題を解くには近似値が必要となるのが通常です。そのような近似値の 1 つがGW 近似です。これは、自己エネルギーがグリーン関数Gと動的遮蔽相互作用Wの積 Σ = GWとして数式化されることからこう呼ばれています。このアプローチは、バンドプロット(およびスペクトル関数などのそれ以外の量)の計算に適しており、完全に第一原理的な方法で定式化することも可能です。重力近似は、絶縁体および半導体のバンドギャップを実験と一致するように与え、それによって系統的なDFTの過小評価を修正するようです。
ほぼ自由電子の近似は電子バンド構造の多くの特性を記述できるが、この理論の1つの結果として、各単位セルの電子数は同じであると予測される。電子数が奇数であれば、各単位セルに不対電子が存在し、価電子帯が完全に占有されていないことが予想され、物質は導体になる。しかし、単位セルあたりの電子数が奇数であるCoOなどの物質は絶縁体であり、この結果と直接矛盾する。この種の物質はモット絶縁体として知られており、この矛盾を説明するには詳細な電子間相互作用(バンド理論では結晶ポテンシャルに対する平均効果としてのみ扱われる)を含める必要がある。ハバード模型はこれらの相互作用を含めることができる近似理論である。これは、ほぼ自由電子近似と原子極限とのギャップを埋めようとする、いわゆる動的平均場理論の中で非摂動論的に扱うことができる。しかし、正式には、この場合の状態は相互作用しないわけではなく、バンド構造の概念はこれらのケースを説明するのに適切ではありません。
バンド構造の計算は、理論固体物理学における重要なテーマです。上記のモデルに加えて、以下のようなモデルもあります。
バンド構造は複素数の波数ベクトルに一般化されており、その結果、表面や界面において興味深い 複素バンド構造と呼ばれるものが生じます
各モデルは、固体の種類によっては極めてよく記述できる一方、そうでないものもあります。ほぼ自由電子モデルは金属には有効ですが、非金属には有効ではありません。強結合モデルは、金属ハロゲン化物塩(例:NaCl)などのイオン性絶縁体に対して非常に正確な記述が可能です。
実空間におけるフェルミ準位に対するバンド構造の変化を理解するために、バンド構造のプロットは、多くの場合、まずバンド図の形で簡略化されます。バンド図では、縦軸はエネルギー、横軸は実空間を表します。水平線はエネルギー準位を表し、ブロックはエネルギーバンドを表します。これらの図の水平線が斜めになっている場合、準位またはバンドのエネルギーは距離とともに変化します。図的に言えば、これは結晶系内に電場が存在することを示しています。バンド図は、異なる物質を接触させた場合の一般的なバンド構造特性を相互に関連付けるのに役立ちます