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:1008.4945 · Physical Review B

Nonstochastic algorithms for Jastrow-Slater and correlator product state wave functions

Eric Neuscamman, Hitesh J. Changlani, Jesse Kinder, Garnet Kin‐Lic Chan

Published Nov 17, 2011• 32 citations
arXivHTMLPDF

Abstract

Jastrow-Slater and correlator product state wave functions, two classes of quantum many-body wave functions, are commonly studied using Monte Carlo methods with the associated drawbacks of stochastic error. Here we show that efficient nonstochastic algorithms for these wave functions exist, both for observable evaluation and for optimization. The algorithms rely on the structure of these states as a product of local, commuting, invertible operators acting on a simple reference wave function. We describe the nonstochastic energy evaluation and optimization algorithms, and demonstrate them with applications to the Heisenberg and spinless and full Hubbard models. Our results demonstrate that the nonstochastic algorithms yield optimized wave functions and energies very close to those obtained with the variational Monte Carlo algorithm. Such algorithms provide new criteria for identifying new classes of wave functions for efficient computational simulation.

DOI
10.1103/physrevb.84.205132

Related papers

A Nonstochastic Optimization Algorithm for Neural-Network Quantum States

Xiang Li, Jia-Cheng Huang, Guang-Ze Zhang +4 more

Neural-network quantum states (NQS) employ artificial neural networks to encode many-body wave functions in a second quantization through variational Monte Carlo (VMC). They have recently been applied to accurately describe electronic wave functions of molecules and have shown the challenges in efficiency compared with traditional quantum chemistry methods. Here, we introduce a general nonstochastic optimization algorithm for NQS in chemical systems, which deterministically generates a selected set of important configurations simultaneously with energy evaluation of NQS. This method bypasses the need for Markov-chain Monte Carlo within the VMC framework, thereby accelerating the entire optimization process. Furthermore, this newly developed nonstochastic optimization algorithm for NQS offers comparable or superior accuracy compared to its stochastic counterpart and ensures more stable convergence. The application of this model to test molecules exhibiting strong electron correlations provides further insight into the performance of NQS in chemical systems and opens avenues for future enhancements.

arXiv:2305.12776•Journal of Chemical Theory and Computation•Nov 14, 2023•16 citations

Variational Monte Carlo simulations using tensor-product projected states

Olga Sikora, Hsueh-Wen Chang, Chung-Pin Chou +2 more

We propose an efficient numerical method, which combines the advantages of recently developed tensor-network based methods and standard trial wave functions, to study the ground-state properties of quantum many-body systems. In this approach, we apply a projector in the form of a tensor-product operator to an input wave function, such as a Jastrow-type or Hartree-Fock wave function, and optimize the tensor elements via variational Monte Carlo. The entanglement already contained in the input wave function can considerably reduce the bond dimensions compared to the regular tensor-product state representation. In particular, this allows us to also represent states that do not obey the area law of entanglement entropy. In addition, for fermionic systems, the fermion sign structure can be encoded in the input wave function. We show that the optimized states provide good approximations of the ground-state energy and correlation functions in the cases of two-dimensional bosonic and fermonic systems.

arXiv:1407.4107•Physical Review B•Apr 8, 2015•18 citations

Improved wave function for strongly correlated electronic systems

Kenji Kobayashi, Kaoru Iguchi

A scheme of optimization of a trial wave function in the variational Monte Carlo calculation is developed. It provides a more precise description of the ground state for strongly correlated electronic systems with highly improved trial wave function. A general type of variational wave function constructed from a linear combination of Slater determinants multiplied by Jastrow correlation factors is proposed. They can explicitly contain low-energy fluctuations and long-range correlations. Using this optimization method, we examine the properties of these wave functions by applying them to a one-dimensional Hubbard model for which the exact solution is available for comparison. It is found that these wave functions prove to be substantial improvements compared to the conventional Jastrow-type wave functions.

Physical review. B, Condensed matter•Jan 15, 1993•14 citations

Computing the energy of a water molecule using MultiDeterminants: A simple, efficient algorithm

Clark, B K, Morales, M A, McMinis, J +2 more

Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave functions are critical to ascertaining new physics. One such wave function is the multi-Slater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally, we implement this method and use it to compute the ground state energy of a water molecule.

•May 3, 2011•60 citations

Computing the energy of a water molecule using multideterminants: A simple, efficient algorithm

Bryan K. Clark, Miguel A. Morales, Jeremy McMinis +2 more

Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave functions are critical to ascertaining new physics. One such wave function is the multi-Slater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally, we implement this method and use it to compute the ground state energy of a water molecule.

arXiv:1106.2456•The Journal of Chemical Physics•Dec 27, 2011•71 citations

Generalized Lanczos algorithm for variational quantum Monte Carlo

Sandro Sorella

We show that the standard Lanczos algorithm can be efficiently implemented statistically and self-consistently improved, using the stochastic reconfiguration method, which has been recently introduced to stabilize the Monte Carlo sign problem instability. With this scheme a few Lanczos steps over a given variational wave function are possible even for large size as a particular case of a more general and more accurate technique that allows to obtain lower variational energies. This method has been tested extensively for a strongly correlated model like the $t\ensuremath{-}J$ model. With the standard Lanczos technique it is possible to compute any kind of correlation functions, with no particular computational effort. By using the fact that the variance $〈{H}^{2}〉\ensuremath{-}〈H{〉}^{2}$ is zero for an exact eigenstate, we show that the approach to the exact solution with few Lanczos iterations is indeed possible even for $\ensuremath{\sim}100$ electrons for reasonably good initial wave functions. The variational stochastic reconfiguration technique presented here allows in general a many-parameter energy optimization of any computable many-body wave function, including for instance generic long-range Jastrow factors and arbitrary site-dependent orbital determinants. This scheme improves further the accuracy of the calculation, especially for long-distance correlation functions.

arXiv:cond-mat/0009149•Physical review. B, Condensed matter•Jun 19, 2001•368 citations
Jastrow-Slater and correlator product state wave functions, two classes of quantum many-body wave functions, are commonly studied using Monte Carlo methods with the associated drawbacks of stochastic error. Here we show that efficient nonstochastic algorithms for these wave functions exist, both for observable evaluation and for optimization. The algorithms rely on the structure of these states as a product of local, commuting, invertible operators acting on a simple reference wave function. We describe the nonstochastic energy evaluation and optimization algorithms, and demonstrate them with applications to the Heisenberg and spinless and full Hubbard models. Our results demonstrate that the nonstochastic algorithms yield optimized wave functions and energies very close to those obtained with the variational Monte Carlo algorithm. Such algorithms provide new criteria for identifying new classes of wave functions for efficient computational simulation.
Neural-network quantum states (NQS) employ artificial neural networks to encode many-body wave functions in a second quantization through variational Monte Carlo (VMC). They have recently been applied to accurately describe electronic wave functions of molecules and have shown the challenges in efficiency compared with traditional quantum chemistry methods. Here, we introduce a general nonstochastic optimization algorithm for NQS in chemical systems, which deterministically generates a selected set of important configurations simultaneously with energy evaluation of NQS. This method bypasses the need for Markov-chain Monte Carlo within the VMC framework, thereby accelerating the entire optimization process. Furthermore, this newly developed nonstochastic optimization algorithm for NQS offers comparable or superior accuracy compared to its stochastic counterpart and ensures more stable convergence. The application of this model to test molecules exhibiting strong electron correlations provides further insight into the performance of NQS in chemical systems and opens avenues for future enhancements.
We propose an efficient numerical method, which combines the advantages of recently developed tensor-network based methods and standard trial wave functions, to study the ground-state properties of quantum many-body systems. In this approach, we apply a projector in the form of a tensor-product operator to an input wave function, such as a Jastrow-type or Hartree-Fock wave function, and optimize the tensor elements via variational Monte Carlo. The entanglement already contained in the input wave function can considerably reduce the bond dimensions compared to the regular tensor-product state representation. In particular, this allows us to also represent states that do not obey the area law of entanglement entropy. In addition, for fermionic systems, the fermion sign structure can be encoded in the input wave function. We show that the optimized states provide good approximations of the ground-state energy and correlation functions in the cases of two-dimensional bosonic and fermonic systems.
A scheme of optimization of a trial wave function in the variational Monte Carlo calculation is developed. It provides a more precise description of the ground state for strongly correlated electronic systems with highly improved trial wave function. A general type of variational wave function constructed from a linear combination of Slater determinants multiplied by Jastrow correlation factors is proposed. They can explicitly contain low-energy fluctuations and long-range correlations. Using this optimization method, we examine the properties of these wave functions by applying them to a one-dimensional Hubbard model for which the exact solution is available for comparison. It is found that these wave functions prove to be substantial improvements compared to the conventional Jastrow-type wave functions.
Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave functions are critical to ascertaining new physics. One such wave function is the multi-Slater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally, we implement this method and use it to compute the ground state energy of a water molecule.
Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave functions are critical to ascertaining new physics. One such wave function is the multi-Slater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally, we implement this method and use it to compute the ground state energy of a water molecule.
We show that the standard Lanczos algorithm can be efficiently implemented statistically and self-consistently improved, using the stochastic reconfiguration method, which has been recently introduced to stabilize the Monte Carlo sign problem instability. With this scheme a few Lanczos steps over a given variational wave function are possible even for large size as a particular case of a more general and more accurate technique that allows to obtain lower variational energies. This method has been tested extensively for a strongly correlated model like the $t\ensuremath{-}J$ model. With the standard Lanczos technique it is possible to compute any kind of correlation functions, with no particular computational effort. By using the fact that the variance $〈{H}^{2}〉\ensuremath{-}〈H{〉}^{2}$ is zero for an exact eigenstate, we show that the approach to the exact solution with few Lanczos iterations is indeed possible even for $\ensuremath{\sim}100$ electrons for reasonably good initial wave functions. The variational stochastic reconfiguration technique presented here allows in general a many-parameter energy optimization of any computable many-body wave function, including for instance generic long-range Jastrow factors and arbitrary site-dependent orbital determinants. This scheme improves further the accuracy of the calculation, especially for long-distance correlation functions.