Electrons in a molecule are like dancers trying to keep their distance from each other. But when two get too close, their dance hits a sharp snag — a so-called Coulomb hole (a deep dip in the interaction energy). Standard quantum calculations try to smooth over this snag by adding more and more dance moves, basis functions, and converge to the exact answer at a depressingly slow pace. The research group proposed not to smooth the hole but to redesign the stage itself. They applied a transformation that mathematically 'imprints' the sharp electron repulsion directly into the rules of the dance — into the Hamiltonian (the description of the system's total energy). This approach is known as the transcorrelated method, but until now it was impractical for large molecules due to the enormous computational cost.
The breakthrough came from two ideas. The first is data compression: the bulky integrals describing interactions were replaced with compact low-rank matrices using interpolative separable density fitting. Simply put, instead of a detailed map of every dancer's step, only key points are stored, from which the entire motion is reconstructed with almost no loss. The second idea is smart handling of triple interactions involving three electrons at once. These were reduced to effective pairwise interactions, drastically lowering complexity. All of this was implemented on several graphics processing units with automatic differentiation, making model construction a routine task even for very large basis sets.
The method was tested on two problems. For a linear chain of hydrogen atoms, calculations with the record basis set cc-pV5Z reached the thermodynamic limit and the complete basis set limit simultaneously — the error was about 1 millihartree per atom. For the benzene molecule (C₆H₆) with 1200 orbitals, accuracy matched the best benchmarks, and extrapolation to the complete basis set became much more reliable than with standard coupled-cluster methods. This means the method is ready to tackle realistic chemical systems where such accuracy was previously only a dream.
as written by the authors 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.
An accuracy of 1 millihartree per atom is like measuring the distance from Earth to the Moon and being off by just the thickness of a human hair.