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

Résumé Technique : Méthode de Coupled Cluster Adaptée avec Diagonalisation Quantique Basée sur l'Échantillonnage

Énoncé du Problème
La prédiction précise de la sélectivité des réactions catalytiques des métallocènes à base de titane pour la production de 1-hexène nécessite le calcul des différences d'énergie libre entre les états de transition concurrents (insertion d'éthylène vs transfert de $\beta$-hydrure) avec une précision de l'ordre de 1 kcal/mol. Les méthodes actuelles de structure électronique sont confrontées à une dichotomie : la théorie de la fonctionnelle de la densité (DFT) est efficace en termes de calcul mais manque souvent de la précision nécessaire en raison de la dépendance aux fonctionnelles, tandis que les méthodes de fonction d'onde de haut niveau (par exemple, CCSD(T)) sont trop coûteuses pour les systèmes de métaux de transition réalistes. De plus, les approches de calcul quantique telles que la Diagonalisation Quantique Basée sur l'Échantillonnage (SQD) sont actuellement limitées par des contraintes matérielles à des espaces actifs trop petits pour capturer toute la corrélation dynamique requise pour des résultats chimiquement précis dans les systèmes catalytiques. Par conséquent, traiter simultanément la corrélation statique (au sein de l'espace actif) et la corrélation dynamique (provenant des orbitales externes) reste un défi important.

Méthodologie
Cette étude propose un cadre hybride quantique-classique qui intègre la Diagonalisation Quantique Basée sur l'Échantillonnage (SQD) avec la théorie du Coupled Cluster Adapté (TCC) pour remédier aux limitations de la taille de l'espace actif et du traitement de la corrélation.

1. Système Cible : Un modèle moléculaire minimal d'un catalyseur de titane ansa-cyclopentadiényl-arène cationique a été utilisé pour étudier les énergies relatives de deux états de transition, T3 (insertion d'éthylène) et T4 (transfert de $\beta$-hydrure). Pour assurer un cadre orbitalaire cohérent, l'état T3 a été augmenté d'une molécule d'éthylène explicite (T3(+C2H4)), correspondant au nombre d'atomes et d'électrons de T4.
2. SQD pour la Corrélation Statique : La SQD a été employée pour traiter la structure électronique de l'espace actif. En utilisant l'ansatz LUCJ (Local-Unitary Cluster Jastrow) sur le processeur quantique IBM Heron, la SQD échantillonne stochastiquement des bitstrings pour construire un Hamiltonien projeté dans un sous-espace. Cela permet le traitement d'espaces actifs (jusqu'à 30 orbitales spatiales) qui sont combinatoirement trop grands pour la méthode classique de Configuration Interaction à Espace Actif Complet (CASCI). La méthode utilise une récupération de configuration auto-cohérente pour atténuer le bruit matériel et garantir le nombre correct d'électrons et de spin.
3. TCC pour la Corrélation Dynamique : Les fonctions d'onde dérivées de la SQD ont servi d'états de référence pour les calculs de Tailored Coupled Cluster (TCC). Dans ce cadre, les amplitudes de cluster pour l'espace actif ($\hat{T}_{AS}$) sont « adaptées » (fixées) sur la base des résultats de la SQD, capturant la corrélation statique et les excitations d'ordre supérieur (triples, quadruples) au sein de l'espace actif. Les amplitudes restantes pour l'espace externe ($\hat{T}_{Ext}$) sont optimisées via les équations standards de Coupled Cluster Singles and Doubles (CCSD) pour capturer la corrélation dynamique.
4. Extensions : L'étude a également appliqué une extension perturbative aux triples, TCC(T), pour améliorer la précision. Un schéma de correction d'énergie a été utilisé pour s'assurer que l'énergie de l'espace actif reflète le niveau SQD plutôt que l'approximation CCSD de niveau inférieur.
5. Configuration de Calcul : Tous les calculs de point unique post-DFT ont été effectués en phase gazeuse pour isoler les effets de corrélation électronique, bien que les optimisations de géométrie incluent des effets de solvant (toluène) via le Modèle de Continuum Polarisable (PCM).

Résultats Clés

• Problèmes de Convergence avec la SQD Seule : Lorsque la SQD était appliquée de manière isolée, les énergies relatives des états de transition ne pouvaient pas converger dans les espaces actifs actuellement réalisables. Les résultats étaient sensibles à la taille de l'espace actif, indiquant que la corrélation statique seule est insuffisante pour ce système.
• Succès de la SQD-TCC : L'approche hybride SQD-TCC a produit des résultats raisonnablement convergés pour les énergies relatives. En incorporant la corrélation dynamique via le cadre TCC, la méthode a réussi à équilibrer le traitement de la corrélation électronique.
• Impact de la Corrélation Statique : Un écart significatif (dépassant 1 kcal/mol) a été observé entre les résultats CCSD(T) standard et les résultats SQD-TCC(T). Cela souligne que les effets de corrélation statique sont critiques dans ce système à base de titane et que négliger ces effets (comme dans le cas du CCSD(T) standard à référence unique) conduit à des erreurs significatives dans la prédiction de la sélectivité.
• Scalabilité de l'Espace Actif : L'étude a démontré que la SQD permet le traitement d'espaces actifs (par exemple, 30 orbitales) qui sont pratiquement inaccessibles aux méthodes CASCI classiques compte tenu des ressources de calcul disponibles (jusqu'à 192 cœurs CPU et 4 TiB de mémoire).

Signification et Revendications
L'article affirme que l'intégration de fonctions d'onde multiréférences dérivées du quantique avec la théorie classique du coupled cluster offre une voie viable vers des simulations chimiquement précises pour des systèmes catalytiques complexes. Plus précisément, ce travail démontre que :

1. Nécessité Hybride : Les corrélations électroniques statiques et dynamiques sont toutes deux essentielles pour des descriptions fiables de la chimie sous-jacente dans la catalyse des métallocènes de titane.
2. Avancement Méthodologique : Le cadre SQD-TCC comble avec succès le fossé entre les capacités actuelles du calcul quantique (espaces actifs limités) et les exigences des simulations catalytiques réalistes (grands espaces actifs + corrélation dynamique).
3. Prédiction de la Sélectivité : Le traitement de la corrélation électronique influence de manière significative la sélectivité prédite de la réaction, la méthode hybride offrant une approche plus robuste que les méthodes classiques de haut niveau pour les systèmes présentant une forte corrélation statique.

Les auteurs concluent que, bien que l'étude actuelle se concentre sur les surfaces d'énergie potentielle en phase gazeuse, la méthodologie établit une base pour des travaux futurs intégrant les contributions d'énergie libre et les effets de solvant au niveau de la fonction d'onde corrélée.