結合クラスター法

結合クラスター法（けつごうクラスターほう、クラスター展開法、CC法：Coupled Cluster）は多体系を記述するために使われる数値手法である。最もよく使われるのは、量子化学（計算化学）におけるポスト-ハートリー-フォック第一原理計算がある。CC法は、ハートリーフォック分子軌道法を基本にして、電子相関を考慮する指数関数クラスター演算子を使って多電子波動関数を構成する。CC法を用いて、小さい分子や中程度の大きさの分子について最も正確な計算を行うことができる。
波動関数.
系のハミルトニアンをformula_1とすると、時間依存しないシュレーディンガー方程式（またはハミルトニアンの固有値方程式）は以下で表される。
ここでformula_3はエネルギー固有状態、formula_4はエネルギー固有値である。多電子系についてはこの方程式は解けない。CC法ではこのエネルギー固有状態を既知の関数で表して、この方程式の近似解を求める。
最低エネルギー状態の波動関数とエネルギーは、それぞれformula_3と"E"で表される。他のCC法（運動方程式結合クラスター法や多参照結合クラスター法など）を用いれば、系の励起状態（と基底状態）の近似解も求めることができる
CC法では多電子系の波動関数を以下のように近似して、励起演算子を求める問題へと変換される。
ここでformula_7は通常はハートリー-フォック分子軌道から構成されたスレーター行列式である。formula_8は励起演算子で、formula_7に作用した場合、様々な励起状態を表すスレーター行列式の線形結合が作られる。詳しくは以下を参照。
配置間相互作用などとは違って、解の示量性を保証するため、この指数関数を用いる方法は適切である。しかしCC法の大きさについての無矛盾性は参照波動関数の大きさの無矛盾性に依存する。CC法の欠点は、変分原理を用いないところである。
クラスター演算子.
クラスター演算子は以下のように表される。
ここでformula_11はすべての1励起の演算子、formula_12は全ての2励起の演算子で、以下続いていく。1粒子励起演算子formula_11と2粒子励起演算子 formula_12 はそれぞれ、ハートリーフォック法で求めた基底状態formula_7を1励起スレーター行列式の線形結合と2励起スレーター行列式の線形結合に変換する。
第二量子化を用いることで、この励起演算子を求める問題は、生成消滅演算子の係数を求める問題へと書き換えることができる。
ここでformula_17とformula_18は生成消滅演算子で、"i, j"は占有軌道を、"a, b"は非占有軌道を表す。近似解formula_3を得るためには未知の係数formula_20とformula_21について解くことが必要である。
指数関数演算子formula_22はテイラー級数に展開できる。例えばformula_8をformula_12の項まで用いた場合、
となる。式には…とあるが、占有軌道の数は有限なので、可能な励起回数も有限であり、この級数は有限である。
"t"を求めるための計算量を少なくするために、formula_8の個々の励起演算子への展開は、3励起ぐらいまでで打ち切ることが多い。このアプローチは、たとえ4励起以上が許されたとしても、演算子へのformula_27, formula_28などの影響は小さいだろうという事実によって保証されている。さらに演算子formula_8の最高励起が"n"である場合、つまり
の場合でも、指数関数演算子のテイラー展開に非線形結合(formula_31など)が含まれているため、"n"回以上励起のスレイター行列式も波動関数formula_3に寄与する。よってformula_33で打ち切られたCC法は、最大"n"励起の配置間相互作用よりも多くの電子相関エネルギーを取り込む。
結合クラスター方程式.
結合クラスターシュレーディンガー方程式は、
結合クラスター方程式の解は、上記の第二量子化の方法だと、係数"t"の組である。そのような方程式はいくらでも作れるが、普通は繰り返し解かれる方程式の組を打ち切る。
未知の"q"個の係数"t"で波動関数を表した場合、"q"個の方程式が必要である。よって係数"t"は、特定の励起行列式に相当することが予想される。formula_35は、占有軌道"i,j,k..."を非占有軌道"a,b,c..." で置き換えることでformula_7から得られる行列式に相当する。よって"q"個の方程式が得られる。
ここでformula_38より、適当な励起行列の組の全体がわかる。これらの方程式の関係を明らかにするため、より分かりやすい形に書き換える。formula_39を結合クラスターシュレーディンガー方程式の両辺に作用させる。formula_40とformula_41に射影すると、
標準的なCCSD法では、
相似変換されたハミルトニアンformula_44は以下で定義され、で書くことができる。
この相似変換されたハミルトニアンはエルミート演算子ではない。一般の量子化学パッケージ（[[ACES (計算化学)|ACES II]]、など）では結合クラスター方程式を繰り返し解く。
CC法の種類.
CC法は、formula_8の定義での最大励起数で分類される。CC法の省略記号は"CC"(coupled cluster)の後ろに以下のような記号を付け加える。
よってCCSDTにおける演算子formula_8は
丸括弧の中の記号は、その記号の部分については[[摂動論]]計算がされたことを意味する。たとえば、CCSD(T) ならば、
という内容を意味する。
外部リンク.
[[Category:量子化学]]