Toward linear scaling auxiliary field quantum Monte Carlo with local natural orbitals

Jo S. Kurian,Hong-Zhou Ye,Ankit Mahajan,Timothy C. Berkelbach,Sandeep Sharma
2023-08-24
Abstract:We develop a local correlation variant of auxiliary field quantum Monte Carlo (AFQMC) that is based on local natural orbitals (LNO-AFQMC). In LNO-AFQMC, independent AFQMC calculations are performed for each localized occupied orbital using a truncated set of tailored orbitals. Because the size of this space does not grow with system size for a target accuracy, the method has linear scaling. Applying LNO AFQMC to molecular problems containing a few hundred to a thousand orbitals, we demonstrate convergence of total energies with significantly reduced costs. The savings are more significant for larger systems and larger basis sets. However, even for our smallest system studied, we find that LNO-AFQMC is cheaper than canonical AFQMC, in contrast with many other reduced-scaling methods. Perhaps most significantly, we show that energy differences converge much more quickly than total energies, making the method ideal for applications in chemistry and material science. Our work paves the way for linear scaling AFQMC calculations of strongly correlated systems, which would have a transformative effect on ab initio quantum chemistry.
Chemical Physics,Strongly Correlated Electrons
What problem does this paper attempt to address?
The paper aims to address the issue of high computational cost in the Auxiliary Field Quantum Monte Carlo (AFQMC) method when dealing with large-scale systems. Specifically, the researchers developed an AFQMC method based on Local Natural Orbitals (LNO), referred to as LNO-AFQMC. This method reduces computational complexity by performing AFQMC calculations independently for each localized occupied orbital and using a truncated set of customized orbitals. Since the spatial size does not increase with the system size, this method exhibits linear scaling properties. The study demonstrated the application of the LNO-AFQMC method on molecular problems containing hundreds to thousands of orbitals, proving the convergence of total energy and significantly reducing computational costs. The cost savings are more pronounced for larger systems and basis sets. Additionally, the research showed that the convergence rate of energy differences is much faster than that of total energy, making this method highly suitable for applications in the fields of chemistry and materials science. This work paves the way for linear-scaling AFQMC calculations of strongly correlated systems, which has a transformative impact on ab initio quantum chemistry.