電子相関のソースを表示

[[Image:Electron correlation.svg|thumb|right|250px|電子相関エネルギーは厳密な非相対論的エネルギーとハートリー‐フォック近似によって得られたエネルギーとの差で定義される。]]
'''電子相関'''（でんしそうかん、{{lang-en-short|Electron correlation}}）とは、多電子系の電子構造における[[電子]]間の相互作用である。また'''電子相関エネルギー'''''E''<sub>corr</sub> とは、多電子系における正確な非相対論的エネルギー''E''<sub>exact</sub> と[[ハートリー‐フォック近似]]によって計算したエネルギー''E''<sub>HF</sub> との差として定義される。
:<math>E_\mathrm{corr} = E_\mathrm{exact} - E_\mathrm{HF}</math>

多電子系における電子間の相互作用をハートリー-フォック法で扱った場合、電子相関の一部しか取り込めていない。

電子相関を考慮するために[[ポスト-ハートリー-フォック法]]が使用される。また、[[コーン–シャム方程式|コーン–シャム]][[密度汎関数理論]]（DFT）は電子相関を電子密度から計算するために相関汎関数を使用する。

== 概要 ==
多電子系において1個の電子の位置を定めると、他の1個の電子の存在確率は第一の電子の位置によって異なる。このことを'''電子相関'''という。[[バンド理論]]を基にした、電子のバンドによる描像では、電子間の相互作用を[[平均場近似]]（例: [[ハートリー‐フォック近似]]）から求める。つまり電子間の位置の相関を平均的なものに置き換えてしまっている。この平均化による近似のため真の電子間相互作用は求まらない。電子相関は、電子間相互作用において、この平均場近似を越える部分のことを指している。

現実の系では、[[d軌道|3d]][[遷移金属]]や、[[電子配置|4f電子]]を持つ[[希土類元素]]からなる[[化合物]]系の電子相関が強く（[[強相関電子系]]）、通常のバンド理論（或いは[[バンド計算]]）では扱えない場合がある。

電子の相関エネルギーは、多電子系での真の[[基底状態]]の（全）[[エネルギー]]（以下、非相対論の場合として考える）から、ハートリー‐フォック近似（[[平均場近似]]を使用）によって得られる基底状態の（全）エネルギーを差し引いたものである。

「相関」は確率理論に由来する概念である<ref name="kutzelnigg">{{cite book|author=Kutzelnigg W. |year=2003|chapter= Theory of Electron Correlation|editor= Rychlewski J. (eds) |title=Explicitly Correlated Wave Functions in Chemistry and Physics|series= Progress in Theoretical Chemistry and Physics, vol 13|publisher= Springer|location= Dordrecht|doi=10.1007/978-94-017-0313-0_1}}</ref>。量子力学の確率解釈という点で電子密度と対密度を定義した時、対密度は個々の電子密度の積と一致しない。つまり、2つの電子は相関していると見なされなければならない。電子相関は以下の2つの主な起源を持つ。
# 電子は「フェルミ統計」に従う。すなわち、[[パウリの排他原理]]を満たす。パウリの排他原理によれば、系の波動関数は粒子の交換に関して反対称化（波動関数の符号が変わる）されていなければならない。
# 全ての電子の対の間にはクーロン反発が存在する。
これら2つの要素を分離できるとすると、それぞれ「フェルミ相関」と「クーロン相関」と呼べるだろう<ref name="kutzelnigg"/>。ここで留意すべき点は、「フェルミ相関」と一般的に呼ばれているもの（平行スピンを有する電子間の積分が0となる）は反対称性とはほとんど関係なく、粒子統計とは独立している（フェルミ粒子とボース粒子の両方に存在する）ということである。この意味でのフェルミ相関は[[スレイター行列式]]で表現された波動関数（フェルミ粒子）と[[ハートリー近似|ハートリー積]]で表現された波動関数（ボース粒子）の両方で起こる。ハートリー-フォック近似で無視されているのは「クーロン相関」である。

== 原子系および分子系 ==
[[量子化学]]の[[ハートリー＝フォック方程式|ハートリー＝フォック法]]内では、反対称[[波動関数]]は単一の[[スレイター行列式]]によって近似される。しかしながら、正確な波動関数は一般に単一の行列式を用いたのでは表現できない。単一の行列式による近似ではクーロン相関が考慮に入らないので、全電子エネルギーは[[ボルン–オッペンハイマー近似]]内の非相対論的[[シュレーディンガー方程式]]の厳密解とは異なる。つまり[[変分原理]]により、ハートリー＝フォック限界は常にこの正確なエネルギーよりも上側にある。この差が「相関エネルギー」と呼ばれる（{{仮リンク|ペル＝オロフ・レフディン|en|Per-Olov Löwdin}}による造語<ref>{{cite journal | first = Per-Olov | last = Löwdin | title = Quantum Theory of Many-Particle Systems. III. Extension of the Hartree–Fock Scheme to Include Degenerate Systems and Correlation Effects  | journal = [[Physical Review]] | volume = 97 | issue = 6 | pages = 1509–1520 | date = March 1955 | doi = 10.1103/PhysRev.97.1509 |bibcode = 1955PhRv...97.1509L }}</ref>）。相関エネルギーの概念はそれ以前に[[ユージン・ウィグナー|ウィグナー]]によって研究されていた<ref>{{Cite journal|last=Wigner|first=E.|date=1934-12-01|title=On the Interaction of Electrons in Metals|journal=Physical Review|volume=46|issue=11|pages=1002–1011|doi=10.1103/PhysRev.46.1002|bibcode=1934PhRv...46.1002W}}</ref>。

既に述べたように通常フェルミ相関と呼ばれる電子相関はHF近似においては既に考慮されており、2つの平行スピン電子は空間中の同じ点には存在できない。一方、クーロン相関はクーロン反発による電子の空間的位置間の相関を記述し、[[ロンドン分散力]]といった化学的に重要な効果の原因である。考慮している系全体の対称性あるいは全スピンに関連する相関も存在する。

相関エネルギーという単語は注意して使用されなければならない。まず、これは普通は相関を考慮した手法とハートリー＝フォックエネルギーとの差として定義される。しかし、HFにも既にある程度の相関が含まれているので、この差は完全な相関エネルギーではない。次に、相関エネルギーは近似に使用する[[基底関数系 (化学)|基底関数系]]への依存性が高い。「正確な」エネルギーは完全な相関と完全な基底関数系によって与えられる。

=== 動的電子相関と静的電子相関 ===
電子相関は動的相関と非動的（静的）相関に分けられることがある。

動的相関は電子の運動の相関（電子–電子の衝突散乱）である。ハートリー＝フォック法ではハミルトニアン演算子の特異性を除去するために必要な波動関数の相関カスプ条件を満たさない。つまり、近距離の電子–電子間相互作用が過剰に取り込まれている。動的電子相関は電子相関ダイナミクス<ref>J.H. McGuire, "Electron Correlation Dynamics in Atomic Collisions", Cambridge University Press, 1997</ref>の下で、[[配置間相互作用]]（CI）法を使って2電子励起配置により取り込むことができる。動的電子相関は正確なエネルギーを議論しなければならない場合に必要となる<ref name="wasada">{{cite journal|和書|title=分子軌道法計算プログラムGaussian 03　―その10―|author-和佐田（筒井） 祐子、和佐田 裕昭|journal=名古屋大学情報連携基盤センターニュース|volume=7|issue=3|year=2008|333–350|url=https://hdl.handle.net/2237/10600}}</ref>。

静的相関は、系の基底状態が2つ以上の（ほぼ）縮退した行列式を用いた場合にだけうまく記述される（本質的に波動関数を多配置にしなければならない<ref name="wasada"/>）分子について重要である。この場合、ハートリー＝フォック波動関数（ただ1つの行列式）は定性的に誤っている。[[多配置自己無撞着場]]（MCSCF）法はこの静的相関が考慮に入るが、動的相関は入らない。

ただし、動的電子相関と静的電子相関をはっきりと区分することはできない<ref name="wasada"/>。例えば、H<sub>2</sub>の結合開裂では、最安定構造での各配置関数の重みと開裂後の各配置関数の重みは連続的に変化するため、最安定構造での動的電子相関と開裂後の静的電子相関の境界を決めることはできない<ref name="wasada"/>。

もし励起エネルギー（基底状態と[[励起状態]]との間のエネルギー差）を計算したいならば、両方の状態が均等にバランスが取られていることを注意しなければならない（例: [[多参照配置間相互作用法]]）。

=== 手法 ===
簡単に言えば、ハートリー＝フォック法の分子軌道は、電子間の瞬間的な反発を含めるのではなく、その他全ての電子の平均場中を運動する個々の分子軌道中の電子のエネルギーを評価することによって最適化される。つまり平均場近似の一種である。

電子相関を考慮するためには、多くの[[ポスト-ハートリー＝フォック法]]が存在する。

==== 配置間相互作用（CI）====
欠けている相関を補正するための最も重要な手法の1つが[[配置間相互作用]]（CI）法である。基底行列式としてハートリー＝フォック波動関数から始め、補正された波動関数として基底行列式と励起行列式の線形結合 <math> \Phi_I</math> を取り、[[変分法 (解析力学)|変分原理]]に従って重み因子 <math>c_I</math> を最適化する。全ての可能な励起行列式を考慮した時は、Full-CIと呼ばれる。Full-CI波動関数中では、全ての電子が完全に相関している。小さくない分子に対しては、Full-CIは計算的にコストが高過ぎる。CI展開を打ち切ると、打ち切りのレベルに応じてよく相関した波動関数とよく相関したエネルギーが得られる。ハートリー＝フォック軌道の基底状態の軌道間で2つの電子の励起に相当する配位までを取り込むCI展開を二電子CI、4つの電子の励起相当の配位までを取り込むCI展開を四電子CI、
6つの電子の励起相当の配位までを取り込むCI展開を六電子CIなどという。しかし取り込む励起の数を増やせば近似の程度は向上するが、電子の多い系ではやはり組合せの数が爆発してしまう。

==== メラー＝プレセット摂動理論（MP2、MP3、MP4等） ====
摂動理論では相関したエネルギーが求まるが、通常はそれに伴う波動関数は求まらない。また摂動理論は変分的ではない。つまり計算されたエネルギーは正確な系のエネルギーの[[有界|上界]]を与えないことを意味する。相互作用量子原子（Interacting Quantum Atoms、IQA）分配を介して[[メラー＝プレセット法|メラー＝プレセット摂動理論エネルギー]]を分配することは可能である（しかし通常は相関エネルギーは分配されない）<ref>{{cite journal|last1=McDonagh|first1=James L.|last2=Vincent|first2=Mark A.|last3=Popelier|first3=Paul L.A.|title=Partitioning dynamic electron correlation energy: Viewing Møller-Plesset correlation energies through Interacting Quantum Atom (IQA) energy partitioning|journal=Chemical Physics Letters|date=October 2016|volume=662|pages=228–234|doi=10.1016/j.cplett.2016.09.019|bibcode=2016CPL...662..228M|doi-access=free}}</ref>。これは、{{仮リンク|Atoms in Molecules|en|Atoms in Molecules}}理論の拡張である。IQAエネルギー分配により、相関エネルギー寄与を個別の原子および原子相互作用から詳細に見ることが可能となる。IQA相関エネルギー分配は、[[結合クラスター法]]でも可能であることが示されている<ref>{{cite journal|last1=Holguín-Gallego|first1=Fernando José|last2=Chávez-Calvillo|first2=Rodrigo|last3=García-Revilla|first3=Marco|last4=Francisco|first4=Evelio|last5=Pendás|first5=Ángel Martín|last6=Rocha-Rinza|first6=Tomás|title=Electron correlation in the interacting quantum atoms partition via coupled-cluster lagrangian densities|journal=Journal of Computational Chemistry|date=15 July 2016|volume=37|issue=19|pages=1753–1765|doi=10.1002/jcc.24372|pmid=27237084|language=en|issn=1096-987X}}</ref><ref>{{cite journal|last1=McDonagh|first1=James L.|last2=Silva|first2=Arnaldo F.|last3=Vincent|first3=Mark A.|last4=Popelier|first4=Paul L. A.|title=Quantifying Electron Correlation of the Chemical Bond|journal=The Journal of Physical Chemistry Letters|volume=8|issue=9|date=12 April 2017|pages=1937–1942|doi=10.1021/acs.jpclett.7b00535|pmid=28402120|issn=1948-7185|url=https://www.research.manchester.ac.uk/portal/en/publications/quantifying-electron-correlation-of-the-chemical-bond(efa8241c-e4a6-48ac-92be-66e1eae59fb5).html}}</ref>。

==== 多配置自己無撞着場（MCSCF）====
複数の手法の組み合わせも可能である。例えば、静的相関を考慮するために[[多配置自己無撞着場]]法のためのいくつかのほぼ縮退した行列式を、動的相関のためにいくつかの打ち切りCI法を、小さな摂動（重要でない）行列式のためにいつかの摂動的[[仮設 (数学)|アンザッツ]]を持つことができる。これらの組み合わせの例には、[[完全活性空間摂動理論|CASPT2]]やSORCIがある。

==== 密度汎関数法 (DF) ====
密度汎関数理論（DFT）に基づく手法である。汎関数をどのようにとるかで色々な変種がある。

==== 結合クラスター展開 (CC) ====

==== 量子モンテカルロ法（QMC） ====

==== テンソルネットワーク (TN) ====

==== ニューラルネットワーク (NN) による機械学習(ML) ====

== 出典 ==
{{脚注ヘルプ}}
{{reflist}}

== 学習用文献 ==
* 川端有郷：「電子相関」、丸善（パリティ物理学コース）、ISBN 4-621-03679-3 (1992年2月29日)。
* 永長直人：「電子相関における場の量子論」、ISBN 978-4-00-005701-1 (1998年10月）。
* 山田耕作：「電子相関」、岩波書店（現代物理学叢書）、ISBN 978-4-00-006749-2 (2000年11月15日）。
* 斯波弘行：「電子相関の物理」、岩波書店、ISBN 978-4-00-007417-9 (2001年6月1日)。
* 藤原毅夫：「固体電子構造論 : 密度汎関数理論から電子相関まで」、内田老鶴圃、ISBN 978-4-7536-2302-0 (2015年4月20日)。

== 関連項目 ==
* [[電子状態]]
* [[完全活性空間]]

{{DEFAULTSORT:てんしそうかん}}
[[Category:固体物理学]]
[[Category:量子化学]]