Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0407184

Optimized Non-Orthogonal Localized Orbitals for Electronic Structure Calculations: Improved Linear Scaling Quantum Monte Carlo

Fernando Agustín Reboredo, Andrew Williamson

Published Jul 7, 2004• 0 citations
arXivHTMLPDF

Abstract

We derive an automatic procedure for generating a set of highly localized, non-orthogonal orbitals for linear scaling quantum Monte Carlo calculations. We demonstrate the advantage of these orbitals in calculations of the total energy of both semiconducting and metallic systems by studying bulk silicon and the homogeneous electron gas. For silicon, the improved localization of these orbitals reduces the computational time by a factor five and the memory by a factor of six compared to localized, orthogonal orbitals. For jellium, we demonstrate that the total energy is converged for orbitals truncated within spheres with radii 7-8 $r_s$, opening the possibility of linear scaling QMC calculations for realistic metallic systems.

DOI
10.48550/arxiv.cond-mat/0407184

Related papers

Optimized nonorthogonal localized orbitals for linear scaling quantum Monte Carlo calculations

Fernando Agustín Reboredo, Andrew Williamson

We derive an automatic procedure for generating a set of highly localized, nonorthogonal orbitals for linear scaling quantum Monte Carlo (QMC) calculations. We demonstrate the advantage of these orbitals for calculating the total energy of both semiconducting and metallic systems by studying bulk silicon and the homogeneous electron gas. For silicon, the improved localization of these orbitals reduces the computational time by a factor of 5 and the memory by a factor of 6 compared to localized, orthogonal orbitals. For jellium at typical metallic densities, we demonstrate that the total energy is converged to 3 meV per electron for orbitals truncated within spheres with radii $7{r}_{s}$, opening the possibility of linear scaling QMC calculations for realistic metallic systems.

Physical Review B•Mar 23, 2005•39 citations

Linear-Scaling Quantum Monte Carlo Calculations

Andrew Williamson, Randolph Q. Hood, Jeffrey C. Grossman

A method is presented for using truncated, maximally localized Wannier functions to introduce sparsity into the Slater determinant part of the trial wave function in quantum Monte Carlo calculations. When combined with an efficient numerical evaluation of these localized orbitals, the dominant cost in the calculation, namely, the evaluation of the Slater determinant, scales linearly with system size. This technique is applied to accurate total energy calculation of hydrogenated silicon clusters and carbon fullerenes containing 20-1000 valence electrons.

Physical Review Letters•Nov 26, 2001•157 citations

Linear-scaling quantum Monte Carlo technique with non-orthogonal localized orbitals

Dario Alfè, Michael J. Gillan

We have reformulated the quantum Monte Carlo (QMC) technique so that a large part of the calculation scales linearly with the number of atoms. The reformulation is related to a recent alternative proposal for achieving linear-scaling QMC, based on maximally localized Wannier orbitals (MLWO), but has the advantage of greater simplicity. The technique we propose draws on methods recently developed for linear-scaling density functional theory. We report tests of the new technique on the insulator MgO, and show that its linear-scaling performance is somewhat better than that achieved by the MLWO approach. Implications for the application of QMC to large complex systems are pointed out.

arXiv:cond-mat/0404578•Journal of Physics Condensed Matter•Jun 7, 2004•38 citations

Toward Linear Scaling Auxiliary-Field Quantum Monte Carlo with Local Natural Orbitals

Jo S. Kurian, Hong‐Zhou Ye, Ankit Mahajan +2 more

We develop a local correlation variant of auxiliary-field quantum Monte Carlo (AFQMC) 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 the 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.

arXiv:2308.12430•Journal of Chemical Theory and Computation•Dec 19, 2023•10 citations

Optimized Non-Orthogonal Localized Orbitals for Linear Scaling Quantum Monte Carlo calculations

Andrew Williamson, Fernando Agustín Reboredo, Giulia Galli

APS•Mar 1, 2004•0 citations

Linear-scaling quantum Monte Carlo calculations

Andrew Williamson, Randolph Q. Hood, Jeffrey C. Grossman

Abstract This study demonstrates linear scaling with system size in fixednode diffusion QMC calculations with pseudopotentials for carbon fullerenes and for hydrogenated silicon clusters with nearly 1000 valence electrons. The range covered is C20 to C18o and SiH4 to Si211H140, The calculations were carried out with standard pseudopotentials (leaving four electrons per C or Si atom) and Slater Jastrow trial functions with orbitals from LDA density functional calculations. The key elements leading to linear scaling were constructions of sparse determinants and simplification of orbital functions. The LDA orbitals were obtained with a plane-wave basis set which was transformed to yield Wannier functions with maximum localization. These functions were then fitted to cubic spline expressions with a limited number of grid points and cut-off distances matching those of the plane-wave orbitals. The computational effort for evaluating trial functions of this type was clearly seen to behave in a linear fashion. Exploration of the effects of variation of cut-off radius revealed that energies, sensitive only to node structures, increased rapidly for cut-offs less than 7 bohr but were unaffected for cutoffs greater than 7 bohr. The accuracies of the calculations were checked in other ways. These included comparison calculations for SiH3 and Si5H12 using Gaussian, original plane-wave, and Wannier function basis sets which gave close agreement for energies.

•Jun 18, 2007•0 citations
We derive an automatic procedure for generating a set of highly localized, non-orthogonal orbitals for linear scaling quantum Monte Carlo calculations. We demonstrate the advantage of these orbitals in calculations of the total energy of both semiconducting and metallic systems by studying bulk silicon and the homogeneous electron gas. For silicon, the improved localization of these orbitals reduces the computational time by a factor five and the memory by a factor of six compared to localized, orthogonal orbitals. For jellium, we demonstrate that the total energy is converged for orbitals truncated within spheres with radii 7-8 rsr_srs​, opening the possibility of linear scaling QMC calculations for realistic metallic systems.
We derive an automatic procedure for generating a set of highly localized, nonorthogonal orbitals for linear scaling quantum Monte Carlo (QMC) calculations. We demonstrate the advantage of these orbitals for calculating the total energy of both semiconducting and metallic systems by studying bulk silicon and the homogeneous electron gas. For silicon, the improved localization of these orbitals reduces the computational time by a factor of 5 and the memory by a factor of 6 compared to localized, orthogonal orbitals. For jellium at typical metallic densities, we demonstrate that the total energy is converged to 3 meV per electron for orbitals truncated within spheres with radii 7rs7{r}_{s}7rs​, opening the possibility of linear scaling QMC calculations for realistic metallic systems.
A method is presented for using truncated, maximally localized Wannier functions to introduce sparsity into the Slater determinant part of the trial wave function in quantum Monte Carlo calculations. When combined with an efficient numerical evaluation of these localized orbitals, the dominant cost in the calculation, namely, the evaluation of the Slater determinant, scales linearly with system size. This technique is applied to accurate total energy calculation of hydrogenated silicon clusters and carbon fullerenes containing 20-1000 valence electrons.
We have reformulated the quantum Monte Carlo (QMC) technique so that a large part of the calculation scales linearly with the number of atoms. The reformulation is related to a recent alternative proposal for achieving linear-scaling QMC, based on maximally localized Wannier orbitals (MLWO), but has the advantage of greater simplicity. The technique we propose draws on methods recently developed for linear-scaling density functional theory. We report tests of the new technique on the insulator MgO, and show that its linear-scaling performance is somewhat better than that achieved by the MLWO approach. Implications for the application of QMC to large complex systems are pointed out.
We develop a local correlation variant of auxiliary-field quantum Monte Carlo (AFQMC) 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 the 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.
Abstract This study demonstrates linear scaling with system size in fixednode diffusion QMC calculations with pseudopotentials for carbon fullerenes and for hydrogenated silicon clusters with nearly 1000 valence electrons. The range covered is C20 to C18o and SiH4 to Si211H140, The calculations were carried out with standard pseudopotentials (leaving four electrons per C or Si atom) and Slater Jastrow trial functions with orbitals from LDA density functional calculations. The key elements leading to linear scaling were constructions of sparse determinants and simplification of orbital functions. The LDA orbitals were obtained with a plane-wave basis set which was transformed to yield Wannier functions with maximum localization. These functions were then fitted to cubic spline expressions with a limited number of grid points and cut-off distances matching those of the plane-wave orbitals. The computational effort for evaluating trial functions of this type was clearly seen to behave in a linear fashion. Exploration of the effects of variation of cut-off radius revealed that energies, sensitive only to node structures, increased rapidly for cut-offs less than 7 bohr but were unaffected for cutoffs greater than 7 bohr. The accuracies of the calculations were checked in other ways. These included comparison calculations for SiH3 and Si5H12 using Gaussian, original plane-wave, and Wannier function basis sets which gave close agreement for energies.