擬ポテンシャル

擬ポテンシャルには、次の2種類がある。

• 経験的に作られるもの（これは第一原理ではない）
• 第一原理計算の結果を利用して作られるもの

現在は精度の上からも、後者が使われることが多い。前述のように第一原理計算は計算コストが高いため、擬ポテンシャルの作製は、もっぱら球対称問題である原子に対して行われる。こうして作製した擬ポテンシャルが、化学結合をふくむ固体中の原子へどの程度利用できるかは、それ自体が複雑な問題である。問題の性質上、明確な答えは存在せず、概ね計算コストとの兼ね合いになる。様々な固体へ適用しても問題を起こさない擬ポテンシャルは、「トランスフェラビリティーが高い」と表現される。

擬ポテンシャルを利用する上での問題点は、内殻電子の寄与を無視するため、内殻電子が関与する物性（内殻励起、コアレベルシフトなど）には擬ポテンシャルを使った手法は事実上無力になることである（コアレベルシフトを擬ポテンシャル手法で扱おうとする試みは存在する）。また、非常に高い圧力下で内殻電子の寄与が物性に影響する（内殻電子の価電子化）ような状況でも、擬ポテンシャルによるバンド計算は対応できなくなる。光学応答の計算についても、擬波動関数が内核領域における波動関数を適切に表現していないため、問題が生じる。

ただし、どこまでを内殻電子とし、どこからを価電子として取り扱うかには任意性が残る。たとえばガリウムの3d軌道など浅い内殻電子を価電子として擬ポテンシャルの手法を利用することは可能である。そうした取り扱いでは、内殻軌道からの効果も部分的に計算に取り込むことが出来る。

エンリコ・フェルミは、原子核による自由中性子の散乱を記述するために擬ポテンシャル${\displaystyle V}$を導入した。^[1] 散乱体から遠く離れた中性子の波動関数は、球面波で表されるs 波の散乱波と入射平面波との和で表されると仮定する。よってポテンシャルは動径${\displaystyle r}$の関数で与えられる。

${\displaystyle V(r)={\frac {2\pi \hbar ^{2}}{m}}b\,\delta (r)}$

ここで${\displaystyle \hbar }$はプランク定数を${\displaystyle 2\pi }$で割ったもの、${\displaystyle m}$ は質量、${\displaystyle \delta (r)}$ は ディラックのデルタ関数、${\displaystyle b}$ は中性子散乱長、${\displaystyle r=0}$ は原子核の重心である。^[2] この${\displaystyle \delta }$関数のフーリエ変換によって中性子の形状因子が得られる。

以上は１つの原子核による中性子の散乱についてである。散乱体が多体系である場合のフェルミ擬ポテンシャルは次のように書ける。

${\displaystyle V_{n}=\sum _{n}{\frac {2\pi \hbar ^{2}}{m}}b_{n}\delta (r-R_{n}(t))}$

• アシュクロフトの擬ポテンシャル
• ハイネ-アバレンコフの擬ポテンシャル

• ノルム保存型擬ポテンシャル
• ウルトラソフト擬ポテンシャル

• 第一原理バンド計算
• Kleinman-Bylander近似
• ゴーストバンド
• 部分内殻補正
• トランスフェラビリティー
• 直交化された平面波 (OPW)
• フローズンコア近似