Coupled Cluster

Die CC-Theorie ist die perturbative Variante der Vielelektronentheorie (MET) von Oktay Sinanoğlu. Da die MET schwierig zu berechnen ist, wird in der heutigen Computerchemie CC verwendet.[2][3][4]

Die CC-Theorie verspricht eine exakte Lösung der zeitunabhängigen Schrödinger-Gleichung:[5]

${\displaystyle H\vert {\Psi }\rangle =E\vert {\Psi }\rangle }$

Die Wellenfunktion wird in der CC-Theorie als Exponentialansatz beschrieben:[5]

${\displaystyle \vert {\Psi }\rangle =e^{T}\vert {\Phi _{0}}\rangle }$,

Wobei ${\displaystyle \vert {\Phi _{0}}\rangle }$, die Referenzwellenfunktion (bzw. das Referenzorbital) ist, welche normalerweise als Slaterdeterminante aus HF-Molekülorbitalen erhalten wird (aber auch andere Wellenfunktionen können benutzt werden, z. B.: aus MCSCF-Rechnungen). ${\displaystyle T}$ ist der Cluster-Operator, der angewendet auf ${\displaystyle \vert {\Phi _{0}}\rangle }$ eine Linearkombination aus angeregten Determinanten erzeugt.

Im Gegensatz zu anderen Ansätzen, wie z. B. der Konfigurationsinteraktion (CI) ist der Exponentialansatz größenkonsistent und -extensiv.[6]

Der Cluster-Operator wird in folgender Form geschrieben:[5]

${\displaystyle T=T_{1}+T_{2}+T_{3}+\cdots }$,

wobei ${\displaystyle T_{1}}$ der Operator aller Einzelanregungen (singles), ${\displaystyle T_{2}}$ der Operator aller Doppelanregungen (doubles),... ist.

Diese Anregungsoperatoren werden ausgedrückt als:[5]

${\displaystyle T_{1}=\sum _{i}\sum _{a}t_{a}^{i}{\hat {a}}^{a}{\hat {a}}_{i},}$

${\displaystyle T_{2}={\frac {1}{4}}\sum _{i,j}\sum _{a,b}t_{ab}^{ij}{\hat {a}}^{a}{\hat {a}}^{b}{\hat {a}}_{j}{\hat {a}}_{i},}$

und für den allgemeinen ${\displaystyle n}$-fachen Clusteroperator:[5]

${\displaystyle T_{n}={\frac {1}{(n!)^{2}}}\sum _{i_{1},i_{2},\ldots ,i_{n}}\sum _{a_{1},a_{2},\ldots ,a_{n}}t_{a_{1},a_{2},\ldots ,a_{n}}^{i_{1},i_{2},\ldots ,i_{n}}{\hat {a}}^{a_{1}}{\hat {a}}^{a_{2}}\ldots {\hat {a}}^{a_{n}}{\hat {a}}_{i_{n}}\ldots {\hat {a}}_{i_{2}}{\hat {a}}_{i_{1}}.}$

Im obigen Ausdruck werden mit ${\displaystyle ({\hat {a}}_{a}^{\dagger }=){\hat {a}}^{a}}$ and ${\displaystyle {\hat {a}}_{i}}$ die Erzeugungs- bzw. Vernichtungsoperatoren bezeichnet. Dabei stehen ${\displaystyle i}$ und ${\displaystyle j}$ für besetzte und ${\displaystyle a}$ sowie ${\displaystyle b}$ für unbesetzte (Teilchen-)Orbitale.

${\displaystyle T_{1}}$(Einteilchen-Cluster-Operator) und ${\displaystyle T_{2}}$ (Zweiteilchen-Cluster-Operator) konvertieren die Referenzfunktion ${\displaystyle \vert {\Phi _{0}}\rangle }$ in eine lineare Kombination der einfach und doppelt angeregten Slater-Determinanten, wenn sie ohne Exponentialfunktion angewendet werden. Wenn man den Exponentialclusteroperator auf die Wellenfunktion anwendet, kann man aufgrund der verschiedenen Potenzen von ${\displaystyle T_{1}}$ und ${\displaystyle T_{2}}$, die in den resultierenden Ausdrücken erscheinen, mehr als doppelt angeregte Determinanten erzeugen. Das Bestimmen der unbekannten Koeffizienten ${\displaystyle t_{a}^{i}}$ und ${\displaystyle t_{ab}^{ij}}$ ist notwendig, um die approximierte Lösung ${\displaystyle \vert {\Psi }\rangle }$ zu finden.

Der Exponentialoperator ${\displaystyle e^{T}}$ kann in einer Taylor-Reihe ausgedrückt werden (hier werden nur ${\displaystyle T_{1}}$ und ${\displaystyle T_{2}}$ berücksichtigt):[5]

${\displaystyle e^{T}=1+T+{\frac {1}{2!}}T^{2}+\cdots =1+T_{1}+T_{2}+{\frac {1}{2}}T_{1}^{2}+T_{1}T_{2}+{\frac {1}{2}}T_{2}^{2}+\cdots }$

Obwohl diese Reihe in der Praxis endlich ist, weil die Anzahl besetzter Molekülorbitale endlich ist, wie auch die Anzahl der Anregungen, ist sie immer noch sehr groß, so dass selbst moderne Hochleistungscomputer Schwierigkeiten mit der Berechnung haben. Die Berechnung wird daher nach einer bestimmten Anzahl an höhere Anregungstypen abgebrochen. So werden z. B. bei CCSD höhere Anregungstypen durch Produkte der singles- und doubles-Koeffizienten approximiert. Die Entwicklung läuft bis unendlich in singles- und doubles-Kombinationen, hat aber keine "reine" triples-Anregung, bzw. quadruples etc.[5]

Die Klassifizierung der traditionellen CC-Methoden beruht auf der höchsten Anzahl von Anregungen, die in der Definition von ${\displaystyle T}$ erlaubt sind. Die Abkürzungen für Coupled Cluster Methoden beginnen gewöhnlich mit den Buchstaben "CC", gefolgt von:

1. S – für eine einzelne Anregung (in der Terminologie der CC-Methode abgekürzt zu singles)
2. D – für doppelte Anregung (doubles)
3. T – für dreifache Anregung (triples)
4. Q – für vierfache Anregung (quadruples)

So hat der ${\displaystyle T}$- Operator in CCSDT die Form:

${\displaystyle T=T_{1}+T_{2}+T_{3}.}$

Buchstaben in runden Klammern weisen darauf hin, dass diese Begriffe auf der Grundlage der Störungstheorie berechnet werden. Zum Beispiel bedeutet CCSD(T):

1. CC mit einer vollständigen Behandlung von singles und doubles.
2. Eine Schätzung des verbundenen Tripelbeitrags wird nicht-iterativ unter der Verwendung der Mehrkörper-Perturbationstheorie berechnet.

Für hinreichend genaue Ergebnisse ist eine CCSD-Berechnung oft ausreichend. Bessere Ergebnisse werden aber z. B. mit CCSDT oder CCSDTQ erzielt, deren Nachteil aber die stark steigenden Kosten sind. In neuster computerchemischen Softwaretechnologien wird eine besondere Art von CCSD(T) angewandt: DLPNO-CCSD(T). DLPNO-CCSD(T)[7] steht für "domain based local pair-natural orbital" CCSD(T). Diese Methode beruht auf einer Annäherung an die originale CCSD Methode, aber durch die semi-kanonische Annäherung ist diese Methode deutlich kosteneffizienter, wodurch es Wissenschaftlern möglich ist, die Berechnung mittelgroßer bis großer Molekülsysteme durchzuführen. Dies war auf Basis der mathematisch akkurateren herkömmlichen CCSD(T) Methode unmöglich und eröffnet ganz neue Anwendungsbereiche in der Forschung, wie zum Beispiel akkurate Scans von Molekülpaaren,[8] viel akkurater, als dies mit DFT-D Methoden möglich ist[9]. CCSD(T) wird als "goldener Standard" der Computerchemie bezeichnet.[10]

1. Young, David C.: Computational chemistry : a practical guide for applying techniques to real world problems. Wiley, New York 2001, ISBN 0-471-33368-9, S. 25.
2. Jiří Čížek: On the Correlation Problem in Atomic and Molecular Systems. Calculation of Wavefunction Components in Ursell-Type Expansion Using Quantum-Field Theoretical Methods. In: The Journal of Chemical Physics. 45. Jahrgang, Nr. 11, 1966, S. 4256, doi:10.1063/1.1727484, bibcode:1966JChPh..45.4256C.
3. O. Sinanoğlu, K. Brueckner: Three approaches to electron correlation in atoms. Yale Univ. Press, 1971, ISBN 0-300-01147-4. and references therein
4. Oktay Si̇nanoğlu: Many-Electron Theory of Atoms and Molecules. I. Shells, Electron Pairs vs Many-Electron Correlations. In: The Journal of Chemical Physics. 36. Jahrgang, Nr. 3, 1962, S. 706, doi:10.1063/1.1732596, bibcode:1962JChPh..36..706S.
5. Jensen, Frank: Introduction to computational chemistry. Third edition Auflage. Chichester West Sussex, UK, ISBN 978-1-118-82599-0, S. 169–174.
6. Trygve Helgaker, Jeppe Olsen, Poul Jorgensen: Molecular Electronic-Structure Theory. Reprint Auflage. Wiley-Blackwell, Chichester 2013, ISBN 978-1-118-53147-1.
7. Yang Guo, Christoph Riplinger, Ute Becker, Dimitrios G. Liakos, Yury Minenkov: Communication: An improved linear scaling perturbative triples correction for the domain based local pair-natural orbital based singles and doubles coupled cluster method [DLPNO-CCSD(T)]. In: The Journal of Chemical Physics. Band 148, Nr. 1, 4. Januar 2018, ISSN 0021-9606, S. 011101, doi:10.1063/1.5011798 (scitation.org [abgerufen am 20. Mai 2021]).
8. Yury Minenkov, Edrisse Chermak, Luigi Cavallo: Accuracy of DLPNO–CCSD(T) Method for Noncovalent Bond Dissociation Enthalpies from Coinage Metal Cation Complexes. In: Journal of Chemical Theory and Computation. Band 11, Nr. 10, 4. September 2015, ISSN 1549-9618, S. 4664–4676, doi:10.1021/acs.jctc.5b00584 (acs.org [PDF; abgerufen am 20. Mai 2021]).
9. Yang Guo, Christoph Riplinger, Ute Becker, Dimitrios G. Liakos, Yury Minenkov: Communication: An improved linear scaling perturbative triples correction for the domain based local pair-natural orbital based singles and doubles coupled cluster method [DLPNO-CCSD(T)]. In: The Journal of Chemical Physics. Band 148, Nr. 1, 4. Januar 2018, ISSN 0021-9606, S. 011101, doi:10.1063/1.5011798 (scitation.org [abgerufen am 20. Mai 2021]).
10. Yang Guo, Christoph Riplinger, Ute Becker, Dimitrios G. Liakos, Yury Minenkov: Communication: An improved linear scaling perturbative triples correction for the domain based local pair-natural orbital based singles and doubles coupled cluster method [DLPNO-CCSD(T)]. In: The Journal of Chemical Physics. Band 148, Nr. 1, 4. Januar 2018, ISSN 0021-9606, S. 011101, doi:10.1063/1.5011798 (scitation.org [abgerufen am 20. Mai 2021]).