Coupled Cluster

Coupled Cluster (CC) bezeichnet eine Methode zur Lösung der Schrödinger-Gleichung, die besonders in der Quantenchemie als eine von vielen post-Hartee-Fock ab initio Methoden verwendet wird. Es handelt sich um ein der Configuration Interaction ähnliches Verfahren. Die Vielteilchen-Wellenfunktion wird dabei in eine Basis aus Slater-Determinanten entwickelt, wodurch die Schrödinger-Gleichung auf ein Matrix-Eigenwertproblem reduziert wird. Die (teilweise) Diagonalisierung dieser Matrix liefert dann die Eigenzustände des quantenmechanischen Systems.<ref name=":0">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. </ref>

Da die Genauigkeit z. B. der Hartree-Fock (HF) Lösungen in der Regel nicht ausreichend ist, muss im Anschluss an z. B. eine HF-Rechnung eine korrelierte Rechnung durchgeführt werden, wobei die vorher berechneten unbesetzten Orbitale zum Einsatz kommen.<ref name=":0" />

Die Methode stammte ursprünglich aus der Kernphysik und wurde dort in den 1950er Jahren von Fritz Coester und Hermann Kümmel eingeführt. In den 1960er Jahren wurde sie von Jiří Čížek (später zusammen mit Josef Paldus) für Elektronenkorrelationen in Atomen und Molekülen formuliert und fand weite Verbreitung in der Quantenchemie. Für die Coupled Cluster Methode (CCM) erhielten Kümmel und Raymond Bishop 2005 die Feenberg-Medaille.

Grundlagen

Die CC-Theorie verspricht eine exakte Lösung der zeitunabhängigen Schrödinger-Gleichung:<ref name=":1">Jensen, Frank: Introduction to computational chemistry. Third edition Auflage. Chichester West Sussex, UK, ISBN 978-1-118-82599-0, S. 169–174. </ref>

<math>H\vert{\Psi}\rangle = E \vert{\Psi}\rangle </math>

Die Wellenfunktion wird in der CC-Theorie als Exponentialansatz beschrieben:<ref name=":1" />

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

Wobei <math>\vert{\Phi_0}\rangle</math>, 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). <math>T</math> ist der Cluster-Operator, der angewendet auf <math>\vert{\Phi_0}\rangle</math> 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.<ref>Trygve Helgaker, Jeppe Olsen, Poul Jorgensen: Molecular Electronic-Structure Theory. Reprint Auflage. Wiley-Blackwell, Chichester 2013, ISBN 978-1-118-53147-1. </ref>

Der Cluster-Operator wird in folgender Form geschrieben:<ref name=":1" />

<math> T=T_1 + T_2 + T_3 + \cdots </math>,

wobei <math>T_1</math> der Operator aller Einzelanregungen (singles), <math>T_2</math> der Operator aller Doppelanregungen (doubles),... ist.

Diese Anregungsoperatoren werden ausgedrückt als:<ref name=":1" />

<math>

T_1=\sum_{i}\sum_{a} t_{a}^{i} \hat{a}^{a}\hat{a}_{i}, </math>

<math>

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}, </math>

und für den allgemeinen <math>n</math>-fachen Clusteroperator:<ref name=":1" />

<math>

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}. </math>

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

<math>T_1</math>(Einteilchen-Cluster-Operator) und <math>T_2</math> (Zweiteilchen-Cluster-Operator) konvertieren die Referenzfunktion <math>\vert{\Phi_0}\rangle</math> 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 <math>T_1</math> und <math>T_2</math>, die in den resultierenden Ausdrücken erscheinen, mehr als doppelt angeregte Determinanten erzeugen. Das Bestimmen der unbekannten Koeffizienten <math>t_{a}^{i}</math> und <math>t_{ab}^{ij}</math> ist notwendig, um die approximierte Lösung <math>\vert{\Psi}\rangle</math> zu finden.

Der Exponentialoperator <math>e^{T}</math> kann in einer Taylor-Reihe ausgedrückt werden (hier werden nur <math>T_1</math> und <math>T_2</math> berücksichtigt):<ref name=":1" />

<math> e^{T} = 1 + T + \frac{1}{2!}T^2 + \cdots = 1 + T_1 + T_2 + \frac{1}{2}T_1^2 + T_1T_2 + \frac{1}{2}T_2^2 + \cdots </math>

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.<ref name=":1" />

Nomenklatur

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

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

So hat der <math>T</math>- Operator in CCSDT die Form:

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

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

Anwendung

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)<ref>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]). </ref> 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,<ref>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]). </ref> viel akkurater, als dies mit DFT-D Methoden möglich ist<ref>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]). </ref>. CCSD(T) wird als "goldener Standard" der Computerchemie bezeichnet.<ref>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]). </ref>

Einzelnachweise

Vorlage:Hinweisbaustein