結合クラスター法のソースを表示

{{電子構造論}}
'''結合クラスター法'''（けつごうクラスターほう）とは、[[多体系]]を記述するために使われる数値手法である。日本語では'''クラスター展開法'''、'''CC法'''(Coupled Cluster)とも称す。

最もよく使われるのは、[[量子化学]]（[[計算化学]]）における[[ポスト-ハートリー-フォック]][[第一原理計算]]がある。CC法は、[[ハートリーフォック]][[分子軌道法]]を基本にして、[[電子相関]]を考慮する指数関数クラスター演算子を使って多電子波動関数を構成する。CC法を用いて、小さい分子や中程度の大きさの分子について最も正確な計算を行うことができる<ref>
{{ cite book | first1 = H. G. |last1 =  Kümmel |chapter = A biography of the coupled cluster method | editor1-first=R. F. | editor1-last = Bishop  | editor2-first= T.  | editor2-last = Brandes
 | editor3-first =  K. A. | editor3-last =  Gernoth | editor4-first =  N. R. |editor4-last = Walet |editor5-first = Y.
| title = Recent progress in many-body theories Proceedings of the 11th international conference| publisher =  World Scientific Publishing | location =  Singapore | year = 2002
| pages =334–348 | isbn= 9789810248888 | editor-last =Xian }}
</ref><ref>{{cite book
  | last = Cramer  | first = Christopher J.
  | title = Essentials of Computational Chemistry
  | publisher = John Wiley & Sons, Ltd.  | year = 2002  | location = Chichester
  | pages = 191–232  | isbn = 0-471-48552-7}}</ref><ref>{{Cite book|isbn=978-0521818322|publisher=Cambridge University Press |title=Many-Body Methods in Chemistry and Physics: MBPT and Coupled-Cluster Theory|year=2009| first1=Isaiah |last1=Shavitt | first2=Rodney J. | last2=Bartlett}}</ref>。

== 波動関数 ==
系の[[ハミルトニアン]]を<math>\hat{H}</math>とすると、時間依存しない[[シュレーディンガー方程式]]（または[[ハミルトニアン]]の[[固有値方程式]]）は以下で表される。
:<math>\hat{H} \vert{\Psi}\rangle = E \vert{\Psi}\rangle </math>
ここで<math>\vert{\Psi}\rangle</math>は[[エネルギー固有状態]]、<math>E \ </math>はエネルギー固有値である。多電子系についてはこの方程式は解けない。CC法ではこのエネルギー固有状態を既知の関数で表して、この方程式の近似解を求める。

最低エネルギー状態の波動関数とエネルギーは、それぞれ<math>\vert{\Psi}\rangle</math>と''E''で表される。他のCC法（[[運動方程式結合クラスター法]]や[[多参照結合クラスター法]]など）を用いれば、系の[[励起状態]]（と[[基底状態]]）の近似解も求めることができる<ref>{{cite journal | last1 =Koch | first1 =Henrik | last2 =Jo̸rgensen | first2 =Poul | title =Coupled cluster response functions | journal =The Journal of Chemical Physics | volume =93 | pages =3333 | year =1990 | doi =10.1063/1.458814|bibcode = 1990JChPh..93.3333K }}</ref>
<ref>{{cite journal |last1 = Stanton |first1 = John F. |last2 = Bartlett |first2 = Rodney J. |title = The equation of motion coupled-cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties |journal = The Journal of Chemical Physics |volume = 98 |pages = 7029 |year = 1993 |doi = 10.1063/1.464746|bibcode = 1993JChPh..98.7029S }}</ref>。

CC法では多電子系の波動関数を以下のように近似して、励起演算子を求める問題へと変換される。

:<math>\vert{\Psi}\rangle = e^{\hat{T}} \vert{\Phi_0}\rangle</math>

ここで<math>\vert{\Phi_0}\rangle</math>は通常はハートリー-フォック[[分子軌道]]から構成された[[スレーター行列式]]である。<math>\hat{T}</math>は励起演算子で、<math>\vert{\Phi_0}\rangle</math>に作用した場合、様々な励起状態を表すスレーター行列式の線形結合が作られる。詳しくは以下を参照。

[[配置間相互作用]]などとは違って、解の[[示量性]]を保証するため、この指数関数を用いる方法は適切である。しかしCC法の[[大きさについての無矛盾性]]は参照波動関数の大きさの無矛盾性に依存する。CC法の欠点は、[[変分原理]]を用いないところである。

== クラスター演算子 ==
クラスター演算子は以下のように表される。
:<math>\hat{T}=\hat{T}_1 + \hat{T}_2 + \hat{T}_3 + \dotsb</math>
ここで<math>\hat{T}_1</math>はすべての1励起の演算子、<math>\hat{T}_2</math>は全ての2励起の演算子で、以下続いていく。1粒子励起演算子<math>\hat{T}_1</math>と2粒子励起演算子 <math>\hat{T}_2</math> はそれぞれ、ハートリーフォック法で求めた基底状態<math>\vert{\Phi_0}\rangle</math>を1励起スレーター行列式の線形結合と2励起スレーター行列式の線形結合に変換する。

[[第二量子化]]を用いることで、この励起演算子を求める問題は、生成消滅演算子の係数を求める問題へと書き換えることができる。
:<math>\begin{align}
\hat{T}_1 &= \sum_{i}\sum_{a} t_{i}^{a} \hat{a}_{i}\hat{a}^{\dagger}_{a} \\
\hat{T}_2 &= \frac{1}{4}\sum_{i,j}\sum_{a,b} t_{ij}^{ab} \hat{a}_{i}\hat{a}_j\hat{a}^{\dagger}_{a}\hat{a}^{\dagger}_{b}, \dotsc \\
\end{align}</math>

ここで<math>\hat{a}^{\dagger}</math>と<math>\hat{a}</math>は[[生成消滅演算子]]で、''i, j''は占有軌道を、''a, b''は非占有軌道を表す。近似解<math>\vert{\Psi}\rangle</math>を得るためには未知の係数<math>t_{i}^{a}</math>と<math>t_{ij}^{ab}</math>について解くことが必要である。

指数関数演算子<math>e^{\hat{T}}</math>は[[テイラー級数]]に展開できる。例えば<math>\hat{T}</math>を<math>\hat{T}_2</math>の項まで用いた場合、

:<math>e^{\hat{T}}
= 1 + \hat{T} + \frac{\hat{T}^2}{2!} + \dotsb
= 1 + \hat{T}_1 + \hat{T}_2 + \frac{\hat{T}_1^2}{2} + \hat{T}_1\hat{T}_2 + \frac{\hat{T}_2^2}{2} + \dotsb</math>

となる。式には…とあるが、占有軌道の数は有限なので、可能な励起回数も有限であり、この級数は有限である。

''t''を求めるための計算量を少なくするために、<math>\hat{T}</math>の個々の励起演算子への展開は、3励起ぐらいまでで打ち切ることが多い。このアプローチは、たとえ4励起以上が許されたとしても、演算子への<math>\hat{T}_5</math>, <math>\hat{T}_6</math>などの影響は小さいだろうという事実によって保証されている。さらに演算子<math>\hat{T}</math>の最高励起が''n''である場合、つまり

:<math>\hat{T} = 1 + \hat{T}_1 + \dotsb + \hat{T}_n </math>

の場合でも、指数関数演算子のテイラー展開に[[非線形]]結合(<math>\hat{T}_1^2 / 2, \hat{T}_1\hat{T}_2, \hat{T}_2^2 / 2</math>など)が含まれているため、''n''回以上励起のスレイター行列式も波動関数<math>\vert{\Psi}\rangle</math>に寄与する。よって<math>\hat{T}_n</math>で打ち切られたCC法は、最大''n''励起の配置間相互作用よりも多くの電子相関エネルギーを取り込む。

== 結合クラスター方程式 ==
結合クラスターシュレーディンガー方程式は、
:<math>\hat{H} e^{\hat{T}} \vert{\Psi_0}\rangle = E e^{\hat{T}} \vert {\Psi_0}\rangle</math>
結合クラスター方程式の解は、上記の第二量子化の方法だと、係数''t''の組である。そのような方程式はいくらでも作れるが、普通は繰り返し解かれる方程式の組を打ち切る。

未知の''q''個の係数''t''で波動関数を表した場合、''q''個の方程式が必要である。よって係数''t''は、特定の励起行列式に相当することが予想される。<math>t_{ijk...}^{abc...}</math>は、占有軌道''i,j,k,...''を非占有軌道''a,b,c,...'' で置き換えることで<math>\vert{\Phi_0}\rangle</math>から得られる行列式に相当する。よって''q''個の方程式が得られる。

:<math>\langle {\Psi^{*}}\vert \hat{H} e^{\hat{T}} \vert{\Psi_0}\rangle = E \langle {\Psi^{*}} \vert e^{\hat{T}} \vert {\Psi_0}\rangle</math>

ここで<math>\vert{\Psi^{*}}\rangle</math>より、適当な励起行列の組の全体がわかる。これらの方程式の関係を明らかにするため、より分かりやすい形に書き換える。<math> e^{-\hat{T}} </math>を結合クラスターシュレーディンガー方程式の両辺に作用させる。<math> \Psi_0 </math>と<math> \Psi^* </math>に射影すると、

:<math>\begin{align}
\langle {\Psi_0}\vert e^{-\hat{T}} \hat{H} e^{\hat{T}} \vert{\Psi_0}\rangle &= E \\
\langle {\Psi^{*}}\vert e^{-\hat{T}} \hat{H} e^{\hat{T}} \vert{\Psi_0}\rangle &= E \langle {\Psi^{*}}\vert e^{-\hat{T}} e^{\hat{T}} \vert{\Psi_0}\rangle = 0
\end{align}</math>
標準的なCCSD法では、

:<math>\begin{align}
\langle {\Psi_0}\vert e^{-(\hat{T}_1+\hat{T}_2)} \hat{H} e^{(\hat{T}_1+\hat{T}_2)} \vert{\Psi_0}\rangle &= E \\
\langle {\Psi_{S}}\vert e^{-(\hat{T}_1+\hat{T}_2)} \hat{H} e^{(\hat{T}_1+\hat{T}_2)} \vert{\Psi_0}\rangle &=0 \\
\langle {\Psi_{D}}\vert e^{-(\hat{T}_1+\hat{T}_2)} \hat{H} e^{(\hat{T}_1+\hat{T}_2)} \vert{\Psi_0}\rangle &=0
\end{align}</math>

相似変換されたハミルトニアン<math>\bar{H}</math>は以下で定義され、{{ill|BCH形式|en|Baker–Campbell–Hausdorff formula}}で書くことができる。
:<math>\bar{H} = e^{-\hat{T}} \hat{H} e^{\hat{T}} = H + [H,T] + (1/2)[[H,T],T] + \dotsb</math>
この相似変換されたハミルトニアンはエルミート演算子ではない。一般の量子化学パッケージ（[[ACES (計算化学)|ACES II]]、{{ill|NWChem|en|NWChem}}など）では結合クラスター方程式を繰り返し解く。

== CC法の種類 ==
CC法は、<math>\hat{T}</math>の定義での最大励起数で分類される。CC法の省略記号は"CC"(coupled cluster)の後ろに以下のような記号を付け加える。
# S - 1励起
# D - 2励起
# T - 3励起
# Q - 4励起
従って、CCSDTにおける演算子<math>\hat{T}</math>は
:<math>\hat{T} = \hat{T}_1 + \hat{T}_2 + \hat{T}_3.</math>
丸括弧の中の記号は、その記号の部分については[[摂動論]]計算がされたことを意味する。たとえば、CCSD(T) ならば、
# 結合クラスター法である。
# 1励起と2励起は完全に含まれている。
# 3励起については摂動論で計算されている。
との内容を意味する。

== 脚注 ==
{{reflist}}

==関連項目==
* [[量子化学]]
* [[配置間相互作用]]

== 外部リンク ==
* [http://www.ccc.uga.edu/lec_top/cc/html/review.html A theoretical review and introduction to coupled-cluster theory]
* [http://www.ua.es/cuantica/docencia/otros/cc/cc.html An Introduction to Coupled-Cluster Theory]
* [http://student.ccbcmd.edu/~bhoffm30/instructional/ci/node14.html The Coupled Cluster (CC) Approach]

{{DEFAULTSORT:けつこうくらすたあほう}}
[[Category:量子化学]]