ハリス汎関数のソースを表示

[[密度汎関数理論]] (DFT) において、'''ハリス汎関数'''（ハリスはんかんすう、{{lang-en-short|Harris functional}}）とは、[[コーン–シャム方程式|コーン–シャム密度汎関数理論]]の非[[自己無撞着]]な近似手法である<ref>{{cite journal | doi = 10.1103/PhysRevB.31.1770 | title = Simplified method for calculating the energy of weakly interacting fragments | year = 1985 | last1 = Harris | first1 = J. | issue=4| journal = Physical Review B | volume = 31 | pages = 1770–1779|bibcode = 1985PhRvB..31.1770H }}</ref>。ハリス汎関数は複合系の全エネルギーを、孤立した部分系の[[電子密度]]の関数として与える。電子密度が収束した状態から変動したとき、コーン-シャム法と比較して、ハリス汎関数により与えられるエネルギーの変動はずっと小さい{{要出典|date=2017-08}}。

厳密な電子密度<math> \rho_0(\vec r) </math>とは異なる、近似的な電子密度<math>\rho(\vec r)</math>が得られていると仮定する。[[電子相関|交換相関ポテンシャル]]<math> v_{\mathrm{xc}}(\vec r) </math>およびハートリーポテンシャル<math> v_{\mathrm H}(\vec r) </math>を近似的な電子密度<math>\rho(\vec r)</math>に基づき構築する。これらのポテンシャルによりコーン-シャム方程式を解くことができ、その固有値が得られる。この固有値の和はしばしばバンドエネルギー<math> E_{\text{band}}</math>と呼ばれる。ここで<math>i</math>は占有コーン-シャム軌道を表す添字である。

:<math>E_{\text{band}} = \sum_i \epsilon_i,</math>

ハリスエネルギー汎関数は次式で表される。

:<math>E_{\text{Harris}} = \sum_i \epsilon_i - \int dr^3 v_{\mathrm{xc}}(\vec r) \rho(\vec r) - \frac{1}{2} \int dr^3 v_{\mathrm H}(\vec r) \rho(\vec r) + E_{\mathrm{xc}}[\rho]</math>

ハリスエネルギー<math>E_{\text{Harris}}</math>と厳密な全エネルギーの差は近似的な電子密度の誤差の2乗<math>O((\rho-\rho_0)^2)</math>にスケールすることが、ハリスにより発見された。それゆえ、多くの系において、ハリス汎関数は十分に正確な値を与えうる。電子密度が変化するような自己無撞着な方法にも適用可能ではあるが、元々はこのようなエネルギーの評価を与えるために考案された。

DFTB+やFireball<ref>{{cite journal | doi = 10.1103/PhysRevB.64.195103 | title = Further developments in the local-orbital density-functional-theory tight-binding method | year = 2001 | last1 = Lewis | first1 = James P. | last2 = Glaesemann | first2 = Kurt R. | last3 = Voth | first3 = Gregory A. | last4 = Fritsch | first4 = Jürgen | last5 = Demkov | first5 = Alexander A. | last6 = Ortega | first6 = José | last7 = Sankey | first7 = Otto F.| journal = Physical Review B | volume = 64 | issue = 19 |pages = 195103|bibcode = 2001PhRvB..64s5103L }}</ref>、Hotbitといった多くの密度汎関数強束縛法(DFTB)はハリス汎関数に基づいている。これらの手法では、多くの場合、自己無撞着なコーン-シャムDFT計算を行わず、ハリスエネルギー汎関数を用いて全エネルギーを見積もる。自己無撞着にコーン-シャム方程式を解く従来のDFTコードと比較して、これらのコードはずっと高速であることが多い。

他にも、ハリス汎関数の考え方を取り入れた手法として、希ガス原子やアルカリ原子を取り扱うためのGordon-Kimポテンシャル関数や、[[分極|分極効果]]を取り入れたPGKRD (Pseudized Gordon-Kim Response Density) 法が挙げられる<ref>{{cite journal | title=分子シミュレーションにおける相互作用エネルギー関数の課題と動向 | year=2009 | | volume=11 | issue=3 | journal=アンサンブル | publisher=分子シミュレーション研究会 | author=三上益弘 | doi=10.11436/mssj.11.3_3 }}</ref>。

コーン–シャム法のエネルギーは[[変分法|変分的]]（基底状態エネルギーよりは決して低い値を取らない）であるのに対して、ハリス汎関数によるエネルギーはanti-varitational（反変分法的: 基底状態エネルギーよりは決して高い値を取らない）であると元々信じられていた<ref>{{cite journal | doi=10.1088/0953-8984/2/10/018 |issue = 10| title = Extremal properties of the Harris energy functional |year=1990 | last1=Zaremba | first1=E. | journal=Journal of Physics: Condensed Matter | volume=2 | pages=2479|bibcode = 1990JPCM....2.2479Z }}</ref>が、これは誤りであるという結論が出ている<ref>{{cite journal | doi=10.1103/PhysRevLett.66.3265 | issue = 25 | pmid=10043743 | title = Does the Harris energy functional possess a local maximum at the ground-state density? | year = 1991 | last1=Robertson | first1 = I. J. | last2 = Farid | first2 = B. | journal = Physical Review Letters | volume = 66 | pages=3265–3268 | bibcode=1991PhRvL..66.3265R}}</ref><ref>{{cite journal | doi=10.1103/PhysRevB.48.11602 | issue = 16 | title = Extremal properties of the Harris-Foulkes functional and an improved screening calculation for the electron gas | year = 1993 | last1=Farid | first1 = B. | last2 = Heine | first2 = V. | last3 =Engel | first3=G. E. | last4=Robertson | first4= I. J.  | journal = Physical Review B | volume = 48 | pages=11602–11621|bibcode = 1993PhRvB..4811602F }}</ref>。

== 出典 ==
{{Reflist|2}}

{{physics-stub}}
{{chem-stub}}
{{DEFAULTSORT:はりすはんかんすう}}
[[Category:密度汎関数理論]]