Tailored coupled cluster method with sample-based quantum diagonalization: Application to titanium-based metallocene catalytic reactions for 1-hexene production

Technisches Resümee: Tailored Coupled Cluster Methode mit Sample-Based Quantum Diagonalization

Problemstellung
Die präzise Vorhersage der Selektivität bei Titan-basierten Metallocen-Katalysereaktionen für die 1-Hexen-Produktion erfordert die Berechnung von Freie-Energie-Differenzen zwischen konkurrierenden Übergangszuständen (Ethylen-Insertion vs. Beta-Hydrid-Transfer) mit einer Genauigkeit innerhalb von 1 kcal/mol. Aktuelle elektronische Strukturmethoden stehen vor einem Dilemma: Die Dichtefunktionaltheorie (DFT) ist recheneffizient, mangelt jedoch oft an der notwendigen Genauigkeit aufgrund der Funktionalabhängigkeit, während hochgradige Wellenfunktionsmethoden (z. B. CCSD(T)) für realistische Übergangsmetallsysteme rechnerisch prohibitiv aufwendig sind. Zudem sind Quantencomputing-Ansätze wie die Sample-Based Quantum Diagonalization (SQD) derzeit durch Hardware-Beschränkungen auf aktive Räume begrenzt, die zu klein sind, um die volle dynamische Korrelation zu erfassen, die für chemisch genaue Ergebnisse in katalytischen Systemen erforderlich ist. Folglich bleibt die gleichzeitige Behandlung sowohl der statischen Korrelation (innerhalb des aktiven Raums) als auch der dynamischen Korrelation (aus externen Orbitalen) eine signifikante Herausforderung.

Methodik
Diese Studie schlägt ein hybrides Quanten-Klassik-Framework vor, das die Sample-Based Quantum Diagonalization (SQD) mit der Tailored Coupled Cluster (TCC)-Theorie integriert, um die Einschränkungen der aktiven Raumgröße und der Korrelationsbehandlung zu adressieren.

1. Zielsystem: Ein minimales molekulares Modell eines kationischen Ansa-Cyclopentadienyl-Aren-Titan-Katalysators wurde verwendet, um die relativen Energien zweier Übergangszustände, T3 (Ethylen-Insertion) und T4 (Beta-Hydrid-Transfer), zu untersuchen. Um einen konsistenten Orbitalrahmen zu gewährleisten, wurde der T3-Zustand durch ein explizites Ethylen-Molekül (T3(+C2H4)) ergänzt, welches die Atom- und Elektronenanzahl von T4 angleicht.
2. SQD für statische Korrelation: SQD wurde eingesetzt, um die Elektronenstruktur des aktiven Raums zu behandeln. Unter Verwendung des LUCJ (Local-Unitary Cluster Jastrow)-Ansatzes auf dem IBM Heron Quantenprozessor wurde SQD genutzt, um Bitstrings stochastisch zu sampeln, um einen projizierten Hamiltonian innerhalb eines Subraums zu konstruieren. Dies ermöglicht die Behandlung von aktiven Räumen (bis zu 30 räumliche Orbitale), die kombinatorisch zu groß für die klassische Complete Active Space Configuration Interaction (CASCI) sind. Die Methode nutzt eine selbstkonsistente Konfigurationswiederherstellung, um Hardware-Rauschen zu mildern und korrekte Elektronenanzahlen sowie Spins sicherzustellen.
3. TCC für dynamische Korrelation: Die aus SQD abgeleiteten Wellenfunktionen dienten als Referenzzustände für Tailored Coupled Cluster (TCC)-Berechnungen. In diesem Framework werden die Cluster-Amplituden des aktiven Raums ($\hat{T}_{AS}$) „maßgeschneidert“ (fixiert) basierend auf den SQD-Ergebnissen, wodurch die statische Korrelation und höhere Anregungen (Triples, Quadruples) innerhalb des aktiven Raums erfasst werden. Die verbleibenden Amplituden für den externen Raum ($\hat{T}_{Ext}$) werden über Standard-Coupled-Cluster-Singles-und-Doubles (CCSD)-Gleichungen optimiert, um die dynamische Korrelation zu erfassen.
4. Erweiterungen: Die Studie wandte zudem eine perturbatieve Triples-Erweiterung, TCC(T), an, um die Genauigkeit zu verbessern. Ein Energiekorrekturschema wurde verwendet, um sicherzustellen, dass die Energie des aktiven Raums die SQD-Ebene widerspiegelt und nicht die niedrigere CCSD-Approximation.
5. Rechenaufbau: Alle Post-DFT-Einpunktberechnungen wurden in der Gasphase durchgeführt, um elektronische Korrelationseffekte zu isolieren, obwohl die Geometrieoptimierungen unter Berücksichtigung von Solvenzeffekten (Toluol) mittels des Polarizable Continuum Model (PCM) erfolgten.

Wesentliche Ergebnisse

• Konvergenzprobleme bei alleiniger Anwendung von SQD: Wenn SQD isoliert angewendet wurde, konnten die relativen Energien der Übergangszustände innerhalb der derzeit praktikablen aktiven Räume nicht konvergiert werden. Die Ergebnisse waren sensitiv gegenüber der Größe des aktiven Raums, was darauf hindeutet, dass statische Korrelation allein für dieses System unzureichend ist.
• Erfolg von SQD-TCC: Der hybride SQD-TCC-Ansatz lieferte vernünftig konvergierte Ergebnisse für die relativen Energien. Durch die Einbeziehung der dynamischen Korrelation über das TCC-Framework gelang es der Methode, die Elektronenkorrelation erfolgreich auszubalancieren.
• Einfluss der statischen Korrelation: Eine signifikante Diskrepanz (über 1 kcal/mol) wurde zwischen den Standard-CCSD(T)- und den SQD-TCC(T)-Ergebnissen beobachtet. Dies verdeutlicht, dass statische Korrelationseffekte in diesem Titan-basierten System entscheidend sind und dass die Vernachlässigung dieser Effekte (wie bei Standard-Single-Reference-CCSD(T)) zu signifikanten Fehlern in der vorhergesagten Selektivität führt.
• Skalierbarkeit des aktiven Raums: Die Studie demonstrierte, dass SQD die Behandlung von aktiven Räumen (z. B. 30 Orbitale) ermöglicht, die klassische CASCI-Methoden angesichts der verfügbaren Rechenressourcen (bis zu 192 CPU-Kerne und 4 TiB Speicher) praktisch nicht zugänglich sind.

Bedeutung und Ansprüche
Die Arbeit behauptet, dass die Integration von quantengestützten Multireferenz-Wellenfunktionen mit der klassischen Coupled-Cluster-Theorie einen lebensfähigen Weg zur chemisch genauen Simulation komplexer katalytischer Systeme darstellt. Insbesondere zeigt die Arbeit, dass:

1. Hybride Notwendigkeit: Sowohl die statische als auch die dynamische Elektronenkorrelation sind essenziell für eine zuverlässige Beschreibung der zugrunde liegenden Chemie in der Titan-Metallocen-Katalyse.
2. Methodische Weiterentwicklung: Das SQD-TCC-Framework überbrückt erfolgreich die Lücke zwischen aktuellen Quantencomputing-Kapazitäten (begrenzte aktive Räume) und den Anforderungen für realistische katalytische Simulationen (große aktive Räume + dynamische Korrelation).
3. Selektivitätsvorhersage: Die Behandlung der Elektronenkorrelation beeinflusst die vorhergesagte Selektivität der Reaktion signifikant, wobei die hybride Methode einen robusteren Ansatz als Standard-High-Level-Klassik-Methoden für Systeme mit starker statischer Korrelation bietet.

Die Autoren kommen zu dem Schluss, dass sich die aktuelle Studie zwar auf Potenzialflächen in der Gasphase konzentriert, die Methodik jedoch die Grundlage für zukünftige Arbeiten bildet, die Freie-Energie-Beiträge und Solvenzeffekte auf der Ebene der korrelierten Wellenfunktion einbeziehen.