Les électrons dans une molécule ressemblent à des danseurs qui cherchent à rester à distance les uns des autres. Mais lorsque deux d'entre eux se rapprochent de très près, leur danse connaît un raté brutal — apparaît ce qu'on appelle le puits de Coulomb (une profonde dépression dans l'énergie d'interaction). Les calculs quantiques ordinaires tentent de lisser ce raté en ajoutant toujours plus de figures de danse — les fonctions de base — et convergent vers la réponse exacte d'une manière désespérément lente. Un groupe de chercheurs a proposé de ne pas lisser le puits, mais de réaménager la scène elle-même. Ils ont appliqué une transformation qui « imprime » mathématiquement la répulsion brutale des électrons directement dans les règles de la danse — dans l'hamiltonien (la description de l'énergie totale du système). Cette approche est connue sous le nom de méthode transcorrélée, mais jusqu'ici elle était peu pratique pour les grandes molécules en raison de l'énormité des calculs.
La percée est venue de deux idées. La première est la compression des données : les volumineuses intégrales décrivant les interactions ont été remplacées par des matrices compactes de rang faible grâce à l'approximation séparable par interpolation de la densité. En termes simples, au lieu d'une carte détaillée de chaque pas du danseur, on ne stocke que les points clés, à partir desquels tout le mouvement se reconstitue presque sans perte. La deuxième idée est le traitement intelligent des interactions triples, où trois électrons interviennent simultanément. Elles ont été réduites à des paires effectives, ce qui a fortement diminué la complexité. Tout cela a été implémenté sur plusieurs processeurs graphiques avec différenciation automatique, et la construction du modèle est devenue une tâche routinière, même pour de très grandes bases.
La méthode a été testée sur deux problèmes. Pour la chaîne linéaire d'atomes de hydrogène, les calculs avec la base record cc-pV5Z ont atteint simultanément la limite thermodynamique et la limite de base complète — l'erreur était d'environ 1 millihartree par atome. Pour la molécule de benzène (C₆H₆) avec 1200 orbitales, la précision s'est révélée au niveau des meilleures références, et l'extrapolation vers la base complète est devenue beaucoup plus fiable que celle des méthodes standard de clusters couplés. Cela signifie que la méthode est prête à affronter des systèmes chimiques réalistes, où une telle précision relevait auparavant du rêve.
tel que rédigé par les auteurs, en anglais arXiv:2607.17314
The transcorrelated (TC) method dramatically accelerates the convergence of correlated calculations toward the complete-basis-set (CBS) limit by folding a Jastrow correlator into the Hamiltonian via a similarity transformation, incorporating the electron--electron cusp into the effective interaction. We make the TC framework practical for large systems and flexible, multi-center correlators by compressing the grid-evaluated TC integrals with the interpolative separable density-fitting (ISDF) approximation, combined with the effective two-body (xTC) treatment of the three-body operator. This low-rank representation reduces storage and integration costs by orders of magnitude, and a multi-GPU implementation with automatic differentiation of the correlator makes the construction routine for large basis sets. We demonstrate the resulting ISDF-xTC-CCSD method on the linear hydrogen chain, reaching the joint thermodynamic and CBS limits with basis sets up to cc-pV5Z in agreement with state-of-the-art many-body references to within about 1~mHa/atom, and on the benzene ground-state energy with up to 1200 orbitals (cc-pCV5Z), where the method attains state-of-the-art accuracy at the coupled cluster singles and doubles level and its CBS extrapolation is markedly more robust than that of conventional coupled-cluster methods.
Une précision de 1 millihartree par atome, c'est comme si vous mesuriez la distance de la Terre à la Lune et que vous vous trompiez de l'épaisseur d'un cheveu humain.