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

Technische Samenvatting: Tailored Coupled Cluster-methode met Sample-Based Quantum Diagonalization

Probleemstelling
De nauwkeurige voorspelling van selectiviteit in titanium-gebaseerde metallocene katalytische reacties voor de productie van 1-hexene vereist het berekenen van vrije-energieverschillen tussen concurrerende overgangstoestanden (ethyleeninsertie versus bèta-hydride-overdracht) met een nauwkeurigheid binnen 1 kcal/mol. Huidige elektronische structuurmethoden kampen met een dichotomie: Density Functional Theory (DFT) is computationeel efficiënt maar mist vaak de benodigde nauwkeurigheid door functional-afhankelijkheid, terwijl hoogwaardige golffunctiemethoden (bijv. CCSD(T)) voor realistische transitiemetaalsystemen computationeel onhaalbaar zijn. Bovendien zijn quantumcomputing-benaderingen zoals Sample-Based Quantum Diagonalization (SQD) momenteel beperkt door hardwarebeperkingen tot actieve ruimtes die te klein zijn om de volledige dynamische correlatie te vangen die vereist is voor chemisch accurate resultaten in katalytische systemen. Bijgevolg blijft het gelijktijdig behandelen van zowel statische correlatie (binnen de actieve ruimte) als dynamische correlatie (van externe orbitalen) een significante uitdaging.

Methodologie
Deze studie stelt een hybride quantum-klassiek kader voor dat Sample-Based Quantum Diagonalization (SQD) integreert met Tailored Coupled Cluster (TCC) theorie om de beperkingen van de omvang van de actieve ruimte en de behandeling van correlatie aan te pakken.

1. Doelsysteem: Een minimaal moleculair model van een kationische ansa-cyclopentadienyl-areen titanium katalysator werd gebruikt om de relatieve energieën van twee overgangstoestanden te onderzoeken, T3 (ethyleeninsertie) en T4 (bèta-hydride-overdracht). Om een consistente orbitaalstructuur te waarborgen, werd de T3-toestand uitgebreid met een expliciet ethyleenmolecuul (T3(+C2H4)), waarbij het aantal atomen en elektronen overeenkomt met T4.
2. SQD voor Statische Correlatie: SQD werd ingezet om de elektronische structuur van de actieve ruimte te behandelen. Met behulp van de LUCJ (Local-Unitary Cluster Jastrow) ansatz op IBM's Heron quantum processor, samplet SQD stochastisch bitstrings om een geprojecteerde Hamiltoniaan binnen een subruimte te construeren. Dit maakt de behandeling van actieve ruimtes (tot 30 ruimtelijke orbitalen) mogelijk die combinatorisch te groot zijn voor klassieke Complete Active Space Configuration Interaction (CASCI). De methode maakt gebruik van self-consistent configuratieherstel om hardwareruis te mitigeren en correcte elektronenaantallen en spin te waarborgen.
3. TCC voor Dynamische Correlatie: De door SQD afgeleide golffuncties dienden als referentietoestanden voor Tailored Coupled Cluster (TCC) berekeningen. In dit kader worden de clusteramplituden voor de actieve ruimte ($\hat{T}_{AS}$) "tailored" (vastgezet) op basis van de SQD-resultaten, waardoor statische correlatie en hogere-orde excitaties (triples, quadruples) binnen de actieve ruimte worden gevangen. De resterende amplitudes voor de externe ruimte ($\hat{T}_{Ext}$) worden geoptimaliseerd via standaard Coupled Cluster Singles and Doubles (CCSD) vergelijkingen om dynamische correlatie te vangen.
4. Extensies: De studie heeft verder een perturbatieve triples-extensie, TCC(T), toegepast om de nauwkeurigheid te verbeteren. Er werd een energiecorrectieschema gebruikt om ervoor te zorgen dat de actieve ruimte-energie de SQD-niveaus reflecteert in plaats van de lagere CCSD-benadering.
5. Computationele Setup: Alle post-DFT single-point berekeningen werden uitgevoerd in de gasfase om elektronische correlatie-effecten te isoleren, hoewel geometrie-optimalisaties rekening hielden met oplosmiddeleffecten (tolueen) via het Polarizable Continuum Model (PCM).

Belangrijkste Resultaten

• Convergentieproblemen met SQD alleen: Wanneer SQD geïsoleerd werd toegepast, konden de relatieve energieën van de overgangstoestanden niet worden geconvergeerd binnen de momenteel haalbare actieve ruimtes. De resultaten waren gevoelig voor de omvang van de actieve ruimte, wat aangeeft dat statische correlatie alleen onvoldoende is voor dit systeem.
• Succes van SQD-TCC: De hybride SQD-TCC benadering leverde redelijk geconvergeerde resultaten voor de relatieve energieën op. Door dynamische correlatie te integreren via het TCC-kader, slaagde de methode erin de behandeling van elektronische correlatie te balanceren.
• Impact van Statische Correlatie: Er werd een significante discrepantie (groter dan 1 kcal/mol) waargenomen tussen de standaard CCSD(T) en de SQD-TCC(T) resultaten. Dit benadrukt dat statische correlatie-effecten cruciaal zijn in dit titanium-gebaseerde systeem en dat het negeren ervan (zoals bij standaard single-reference CCSD(T)) leidt tot significante fouten in de voorspelde selectiviteit.
• Schaalbaarheid van de Actieve Ruimte: De studie demonstreerde dat SQD de behandeling van actieve ruimtes (bijv. 30 orbitalen) mogelijk maakt die praktisch ontoegankelijk zijn voor klassieke CASCI-methoden gezien de beschikbare computationele middelen (tot 192 CPU-cores en 4 TiB geheugen).

Betekenis en Claims
Het artikel beweert dat de integratie van quantum-afgeleide multireferentie golffuncties met klassieke coupled cluster theorie een levensvatbaar pad biedt naar chemisch accurate simulaties voor complexe katalytische systemen. Specifiek toont het werk aan dat:

1. Hybride Noodzaak: Zowel statische als dynamische elektronische correlaties zijn essentieel voor een betrouwbare beschrijving van de onderliggende chemie in titanium-metallocene katalyse.
2. Methodologische Vooruitgang: Het SQD-TCC framework overbrugt succesvol de kloof tussen de huidige mogelijkheden van quantum computing (beperkte actieve ruimtes) en de vereisten voor realistische katalytische simulaties (grote actieve ruimtes + dynamische correlatie).
3. Selectievoorspelling: De behandeling van elektronische correlatie beïnvloedt de voorspelde selectiviteit van de reactie aanzienlijk, waarbij de hybride methode een robuustere benadering biedt dan standaard klassieke high-level methoden voor systemen met sterke statische correlatie.

De auteurs concluderen dat hoewel de huidige studie zich richt op potentiele energie-oppervlakken in de gasfase, de methodologie een fundament legt voor toekomstig werk dat vrije-energiebijdragen en oplosmiddeleffecten op het niveau van gecorreleerde golffuncties incorporeert.