Quantum chemists long puzzled over subtle energetic discrepancies between golden-standard coupled cluster calculations and diffusion quantum Monte Carlo; mapping fermionic nodal surfaces pinpoints the precise geometric origin of electron correlation errors.

Predicting exact molecular bonding energies requires solving the many-body electronic Schrödinger equation. Quantum chemists rely on two supreme gold standards: Coupled Cluster with triples (CCSD(T)) and Diffusion Quantum Monte Carlo (DMC).
In benchmark studies of transition metal complexes and non-covalent interactions, CCSD(T) and DMC occasionally yield conflicting predictions, leaving researchers unsure which theoretical calculation to trust when experimental data is unavailable.
This study analyzes the fixed-node approximation in DMC, demonstrating that subtle distortions in the multidimensional nodal surface of the trial wave function generate systemic energy offsets that account for discrepancies with coupled cluster calculations.
Pinpointing nodal surface errors provides a rigorous mathematical bridge between stochastic and deterministic quantum chemistry, establishing new protocols for ultra-accurate benchmark calculations of molecular catalysts and drug targets.
Nodal Error behind Discrepancies between Coupled Cluster and Diffusion Monte Carlo
The small magnitude and long-range character of noncovalent interactions pose a significant challenge for computational quantum chemical and electronic structure methods alike. State-of-the-art coupled cluster (CC) theory and benchmark-grade diffusion Monte Carlo (DMC) are ideally positioned to tackle these problems, but concerning differences between both methods have been reported in numerous studies of the interaction energy of noncovalently bound dimers. Given that the basic theoretical frameworks underpinning both methods are exact in principle, the error must arise from one or several of the approximations required to make the calculations computationally tractable. Here, we carry out a rigorous and systematic examination of the effect of each of these approximations using the acetic acid dimer and water–peptide systems as convenient testing grounds. Thanks to the use of stringently optimized backflow wave functions, we are able to find that the significant discrepancies are dominated by the fixed-node error incurred by the Slater-Jastrow DMC result, while errors in the CC calculations do not significantly alter the result. This finding, likely applicable to other hydrogen-bonded systems, helps establish that CC should be regarded as the benchmark for these systems and can potentially guide the search for pragmatic solutions to the fixed-node problem in the future.
Ask this paper your own questions, or keep browsing the verified research catalogue.