結合クラスター法 (けつごうクラスターほう、クラスター展開法 、CC法 :Coupled Cluster)は多体系 を記述するために使われる数値手法である。最もよく使われるのは、量子化学 (計算化学 )におけるポスト-ハートリー-フォック 第一原理計算 がある。CC法は、ハートリーフォック 分子軌道法 を基本にして、電子相関 を考慮する指数関数クラスター演算子を使って多電子波動関数を構成する。CC法を用いて、小さい分子や中程度の大きさの分子について最も正確な計算を行うことができる[1] [2] [3] 。
波動関数 [ 編集 ] 系のハミルトニアン を H ^ {\displaystyle {\hat {H}}} とすると、時間依存しないシュレーディンガー方程式 (またはハミルトニアン の固有値方程式 )は以下で表される。
H ^ | Ψ ⟩ = E | Ψ ⟩ {\displaystyle {\hat {H}}\vert {\Psi }\rangle =E\vert {\Psi }\rangle } ここで | Ψ ⟩ {\displaystyle \vert {\Psi }\rangle } はエネルギー固有状態 、 E {\displaystyle E\ } はエネルギー固有値である。多電子系についてはこの方程式は解けない。CC法ではこのエネルギー固有状態を既知の関数で表して、この方程式の近似解を求める。
最低エネルギー状態の波動関数とエネルギーは、それぞれ | Ψ ⟩ {\displaystyle \vert {\Psi }\rangle } とE で表される。他のCC法(運動方程式結合クラスター法 や多参照結合クラスター法 など)を用いれば、系の励起状態 (と基底状態 )の近似解も求めることができる[4] [5] 。
CC法では多電子系の波動関数を以下のように近似して、励起演算子を求める問題へと変換される。
| Ψ ⟩ = e T ^ | Φ 0 ⟩ {\displaystyle \vert {\Psi }\rangle =e^{\hat {T}}\vert {\Phi _{0}}\rangle } ここで | Φ 0 ⟩ {\displaystyle \vert {\Phi _{0}}\rangle } は通常はハートリー-フォック分子軌道 から構成されたスレーター行列式 である。 T ^ {\displaystyle {\hat {T}}} は励起演算子で、 | Φ 0 ⟩ {\displaystyle \vert {\Phi _{0}}\rangle } に作用した場合、様々な励起状態を表すスレーター行列式の線形結合が作られる。詳しくは以下を参照。
配置間相互作用 などとは違って、解の示量性 を保証するため、この指数関数を用いる方法は適切である。しかしCC法の大きさについての無矛盾性 は参照波動関数の大きさの無矛盾性に依存する。CC法の欠点は、変分原理 を用いないところである。
クラスター演算子 [ 編集 ] クラスター演算子は以下のように表される。
T ^ = T ^ 1 + T ^ 2 + T ^ 3 + ⋯ {\displaystyle {\hat {T}}={\hat {T}}_{1}+{\hat {T}}_{2}+{\hat {T}}_{3}+\dotsb } ここで T ^ 1 {\displaystyle {\hat {T}}_{1}} はすべての1励起の演算子、 T ^ 2 {\displaystyle {\hat {T}}_{2}} は全ての2励起の演算子で、以下続いていく。1粒子励起演算子 T ^ 1 {\displaystyle {\hat {T}}_{1}} と2粒子励起演算子 T ^ 2 {\displaystyle {\hat {T}}_{2}} はそれぞれ、ハートリーフォック法で求めた基底状態 | Φ 0 ⟩ {\displaystyle \vert {\Phi _{0}}\rangle } を1励起スレーター行列式の線形結合と2励起スレーター行列式の線形結合に変換する。
第二量子化 を用いることで、この励起演算子を求める問題は、生成消滅演算子の係数を求める問題へと書き換えることができる。
T ^ 1 = ∑ i ∑ a t i a a ^ i a ^ a † T ^ 2 = 1 4 ∑ i , j ∑ a , b t i j a b a ^ i a ^ j a ^ a † a ^ b † , … {\displaystyle {\begin{aligned}{\hat {T}}_{1}&=\sum _{i}\sum _{a}t_{i}^{a}{\hat {a}}_{i}{\hat {a}}_{a}^{\dagger }\\{\hat {T}}_{2}&={\frac {1}{4}}\sum _{i,j}\sum _{a,b}t_{ij}^{ab}{\hat {a}}_{i}{\hat {a}}_{j}{\hat {a}}_{a}^{\dagger }{\hat {a}}_{b}^{\dagger },\dotsc \\\end{aligned}}} ここで a ^ † {\displaystyle {\hat {a}}^{\dagger }} と a ^ {\displaystyle {\hat {a}}} は生成消滅演算子 で、i, j は占有軌道を、a, b は非占有軌道を表す。近似解 | Ψ ⟩ {\displaystyle \vert {\Psi }\rangle } を得るためには未知の係数 t i a {\displaystyle t_{i}^{a}} と t i j a b {\displaystyle t_{ij}^{ab}} について解くことが必要である。
指数関数演算子 e T ^ {\displaystyle e^{\hat {T}}} はテイラー級数 に展開できる。例えば T ^ {\displaystyle {\hat {T}}} を T ^ 2 {\displaystyle {\hat {T}}_{2}} の項まで用いた場合、
e T ^ = 1 + T ^ + T ^ 2 2 ! + ⋯ = 1 + T ^ 1 + T ^ 2 + T ^ 1 2 2 + T ^ 1 T ^ 2 + T ^ 2 2 2 + ⋯ {\displaystyle 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 } となる。式には…とあるが、占有軌道の数は有限なので、可能な励起回数も有限であり、この級数は有限である。
t を求めるための計算量を少なくするために、 T ^ {\displaystyle {\hat {T}}} の個々の励起演算子への展開は、3励起ぐらいまでで打ち切ることが多い。このアプローチは、たとえ4励起以上が許されたとしても、演算子への T ^ 5 {\displaystyle {\hat {T}}_{5}} , T ^ 6 {\displaystyle {\hat {T}}_{6}} などの影響は小さいだろうという事実によって保証されている。さらに演算子 T ^ {\displaystyle {\hat {T}}} の最高励起がn である場合、つまり
T ^ = 1 + T ^ 1 + ⋯ + T ^ n {\displaystyle {\hat {T}}=1+{\hat {T}}_{1}+\dotsb +{\hat {T}}_{n}} の場合でも、指数関数演算子のテイラー展開に非線形 結合( T ^ 1 2 / 2 , T ^ 1 T ^ 2 , T ^ 2 2 / 2 {\displaystyle {\hat {T}}_{1}^{2}/2,{\hat {T}}_{1}{\hat {T}}_{2},{\hat {T}}_{2}^{2}/2} など)が含まれているため、n 回以上励起のスレイター行列式も波動関数 | Ψ ⟩ {\displaystyle \vert {\Psi }\rangle } に寄与する。よって T ^ n {\displaystyle {\hat {T}}_{n}} で打ち切られたCC法は、最大n 励起の配置間相互作用よりも多くの電子相関エネルギーを取り込む。
結合クラスター方程式 [ 編集 ] 結合クラスターシュレーディンガー方程式は、
H ^ e T ^ | Ψ 0 ⟩ = E e T ^ | Ψ 0 ⟩ {\displaystyle {\hat {H}}e^{\hat {T}}\vert {\Psi _{0}}\rangle =Ee^{\hat {T}}\vert {\Psi _{0}}\rangle } 結合クラスター方程式の解は、上記の第二量子化の方法だと、係数t の組である。そのような方程式はいくらでも作れるが、普通は繰り返し解かれる方程式の組を打ち切る。
未知のq 個の係数t で波動関数を表した場合、q 個の方程式が必要である。よって係数t は、特定の励起行列式に相当することが予想される。 t i j k . . . a b c . . . {\displaystyle t_{ijk...}^{abc...}} は、占有軌道i,j,k,... を非占有軌道a,b,c,... で置き換えることで | Φ 0 ⟩ {\displaystyle \vert {\Phi _{0}}\rangle } から得られる行列式に相当する。よってq 個の方程式が得られる。
⟨ Ψ ∗ | H ^ e T ^ | Ψ 0 ⟩ = E ⟨ Ψ ∗ | e T ^ | Ψ 0 ⟩ {\displaystyle \langle {\Psi ^{*}}\vert {\hat {H}}e^{\hat {T}}\vert {\Psi _{0}}\rangle =E\langle {\Psi ^{*}}\vert e^{\hat {T}}\vert {\Psi _{0}}\rangle } ここで | Ψ ∗ ⟩ {\displaystyle \vert {\Psi ^{*}}\rangle } より、適当な励起行列の組の全体がわかる。これらの方程式の関係を明らかにするため、より分かりやすい形に書き換える。 e − T ^ {\displaystyle e^{-{\hat {T}}}} を結合クラスターシュレーディンガー方程式の両辺に作用させる。 Ψ 0 {\displaystyle \Psi _{0}} と Ψ ∗ {\displaystyle \Psi ^{*}} に射影すると、
⟨ Ψ 0 | e − T ^ H ^ e T ^ | Ψ 0 ⟩ = E ⟨ Ψ ∗ | e − T ^ H ^ e T ^ | Ψ 0 ⟩ = E ⟨ Ψ ∗ | e − T ^ e T ^ | Ψ 0 ⟩ = 0 {\displaystyle {\begin{aligned}\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{aligned}}} 標準的なCCSD法では、
⟨ Ψ 0 | e − ( T ^ 1 + T ^ 2 ) H ^ e ( T ^ 1 + T ^ 2 ) | Ψ 0 ⟩ = E ⟨ Ψ S | e − ( T ^ 1 + T ^ 2 ) H ^ e ( T ^ 1 + T ^ 2 ) | Ψ 0 ⟩ = 0 ⟨ Ψ D | e − ( T ^ 1 + T ^ 2 ) H ^ e ( T ^ 1 + T ^ 2 ) | Ψ 0 ⟩ = 0 {\displaystyle {\begin{aligned}\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{aligned}}} 相似変換されたハミルトニアン H ¯ {\displaystyle {\bar {H}}} は以下で定義され、BCH形式 (英語版 ) で書くことができる。
H ¯ = e − T ^ H ^ e T ^ = H + [ H , T ] + ( 1 / 2 ) [ [ H , T ] , T ] + ⋯ {\displaystyle {\bar {H}}=e^{-{\hat {T}}}{\hat {H}}e^{\hat {T}}=H+[H,T]+(1/2)[[H,T],T]+\dotsb } この相似変換されたハミルトニアンはエルミート演算子ではない。一般の量子化学パッケージ(ACES II 、NWChem (英語版 ) など)では結合クラスター方程式を繰り返し解く。
CC法の種類 [ 編集 ] CC法は、 T ^ {\displaystyle {\hat {T}}} の定義での最大励起数で分類される。CC法の省略記号は"CC"(coupled cluster)の後ろに以下のような記号を付け加える。
S - 1励起 D - 2励起 T - 3励起 Q - 4励起 よってCCSDTにおける演算子 T ^ {\displaystyle {\hat {T}}} は
T ^ = T ^ 1 + T ^ 2 + T ^ 3 . {\displaystyle {\hat {T}}={\hat {T}}_{1}+{\hat {T}}_{2}+{\hat {T}}_{3}.} 丸括弧の中の記号は、その記号の部分については摂動論 計算がされたことを意味する。たとえば、CCSD(T) ならば、
結合クラスター法である。 1励起と2励起は完全に含まれている。 3励起については摂動論で計算されている。 という内容を意味する。
^ 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 ^ Cramer, Christopher J. (2002). Essentials of Computational Chemistry . Chichester: John Wiley & Sons, Ltd.. pp. 191–232. ISBN 0-471-48552-7 ^ Shavitt, Isaiah; Bartlett, Rodney J. (2009). Many-Body Methods in Chemistry and Physics: MBPT and Coupled-Cluster Theory . Cambridge University Press. ISBN 978-0521818322 ^ Koch, Henrik; Jo̸rgensen, Poul (1990). “Coupled cluster response functions”. The Journal of Chemical Physics 93 : 3333. Bibcode : 1990JChPh..93.3333K . doi :10.1063/1.458814 . ^ 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. Bibcode : 1993JChPh..98.7029S . doi :10.1063/1.464746 . 関連項目 [ 編集 ] 外部リンク [ 編集 ]