結合クラスター法

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

最もよく使われるのは、量子化学計算化学)におけるポスト-ハートリー-フォック第一原理計算がある。CC法は、ハートリーフォック分子軌道法を基本にして、電子相関を考慮する指数関数クラスター演算子を使って多電子波動関数を構成する。CC法を用いて、小さい分子や中程度の大きさの分子について最も正確な計算を行うことができる[1][2][3]

波動関数

[編集]

系のハミルトニアンとすると、時間依存しないシュレーディンガー方程式(またはハミルトニアン固有値方程式)は以下で表される。

ここでエネルギー固有状態はエネルギー固有値である。多電子系についてはこの方程式は解けない。CC法ではこのエネルギー固有状態を既知の関数で表して、この方程式の近似解を求める。

最低エネルギー状態の波動関数とエネルギーは、それぞれEで表される。他のCC法(運動方程式結合クラスター法多参照結合クラスター法など)を用いれば、系の励起状態(と基底状態)の近似解も求めることができる[4] [5]

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

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

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

クラスター演算子

[編集]

クラスター演算子は以下のように表される。

ここではすべての1励起の演算子、は全ての2励起の演算子で、以下続いていく。1粒子励起演算子と2粒子励起演算子 はそれぞれ、ハートリーフォック法で求めた基底状態を1励起スレーター行列式の線形結合と2励起スレーター行列式の線形結合に変換する。

第二量子化を用いることで、この励起演算子を求める問題は、生成消滅演算子の係数を求める問題へと書き換えることができる。

ここで生成消滅演算子で、i, jは占有軌道を、a, bは非占有軌道を表す。近似解を得るためには未知の係数について解くことが必要である。

指数関数演算子テイラー級数に展開できる。例えばの項まで用いた場合、

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

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

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

結合クラスター方程式

[編集]

結合クラスターシュレーディンガー方程式は、

結合クラスター方程式の解は、上記の第二量子化の方法だと、係数tの組である。そのような方程式はいくらでも作れるが、普通は繰り返し解かれる方程式の組を打ち切る。

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

ここでより、適当な励起行列の組の全体がわかる。これらの方程式の関係を明らかにするため、より分かりやすい形に書き換える。を結合クラスターシュレーディンガー方程式の両辺に作用させる。に射影すると、

標準的なCCSD法では、

相似変換されたハミルトニアンは以下で定義され、BCH形式英語版で書くことができる。

この相似変換されたハミルトニアンはエルミート演算子ではない。一般の量子化学パッケージ(ACES IINWChem英語版など)では結合クラスター方程式を繰り返し解く。

CC法の種類

[編集]

CC法は、の定義での最大励起数で分類される。CC法の省略記号は"CC"(coupled cluster)の後ろに以下のような記号を付け加える。

  1. S - 1励起
  2. D - 2励起
  3. T - 3励起
  4. Q - 4励起

従って、CCSDTにおける演算子

丸括弧の中の記号は、その記号の部分については摂動論計算がされたことを意味する。たとえば、CCSD(T) ならば、

  1. 結合クラスター法である。
  2. 1励起と2励起は完全に含まれている。
  3. 3励起については摂動論で計算されている。

との内容を意味する。

脚注

[編集]
  1. ^ Kümmel, H. G. (2002). “A biography of the coupled cluster method”. In Xian, R. F.; Brandes, T.; Gernoth, K. A. et al.. Recent progress in many-body theories Proceedings of the 11th international conference. Singapore: World Scientific Publishing. pp. 334–348. ISBN 9789810248888 
  2. ^ Cramer, Christopher J. (2002). Essentials of Computational Chemistry. Chichester: John Wiley & Sons, Ltd.. pp. 191–232. ISBN 0-471-48552-7 
  3. ^ Shavitt, Isaiah; Bartlett, Rodney J. (2009). Many-Body Methods in Chemistry and Physics: MBPT and Coupled-Cluster Theory. Cambridge University Press. ISBN 978-0521818322 
  4. ^ Koch, Henrik; Jo̸rgensen, Poul (1990). “Coupled cluster response functions”. The Journal of Chemical Physics 93: 3333. Bibcode1990JChPh..93.3333K. doi:10.1063/1.458814. 
  5. ^ Stanton, John F.; Bartlett, Rodney J. (1993). “The equation of motion coupled-cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties”. The Journal of Chemical Physics 98: 7029. Bibcode1993JChPh..98.7029S. doi:10.1063/1.464746. 

関連項目

[編集]

外部リンク

[編集]