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

Projector quantum Monte Carlo with matrix product states

Sebastian Wouters, Brecht Verstichel, Dimitri Van Neck, Garnet Kin‐Lic Chan

Published Jul 3, 2014• 29 citations
arXivHTMLPDF

Abstract

We marry tensor network states (TNS) and projector quantum Monte Carlo (PMC) to overcome the high computational scaling of TNS and the sign problem of PMC. Using TNS as trial wave functions provides a route to systematically improve the sign structure and to eliminate the bias in fixed-node and constrained-path PMC. As a specific example, we describe phaseless auxiliary-field quantum Monte Carlo with matrix product states (MPS-AFQMC). MPS-AFQMC improves significantly on the density-matrix renormalization group (DMRG) ground-state energy. For the ${J}_{1}\ensuremath{-}{J}_{2}$ model on two-dimensional square lattices, we observe with MPS-AFQMC an order of magnitude reduction in the error for all couplings, compared to DMRG. The improvement is independent of walker bond dimension, and we therefore use bond dimension 1 for the walkers. The computational cost of MPS-AFQMC is then quadratic in the bond dimension of the trial wave function, which is lower than the cubic scaling of DMRG. The error due to the constrained-path bias is proportional to the variational error of the trial wave function. We show that for the ${J}_{1}\ensuremath{-}{J}_{2}$ model on two-dimensional square lattices, a linear extrapolation of the MPS-AFQMC energy with the discarded weight from the DMRG calculation allows us to remove the constrained-path bias. Extensions to other tensor networks are briefly discussed.

DOI
10.1103/physrevb.90.045104

Related papers

Unbiasing fermionic auxiliary-field quantum Monte Carlo with matrix product state trial wavefunctions

Tong Jiang, Bryan O’Gorman, Ankit Mahajan +1 more

In this work, we report, for the first time, an implementation of fermionic auxiliary-field quantum Monte Carlo (AFQMC) using matrix product state (MPS) trial wavefunctions, dubbed MPS-AFQMC. Calculating overlaps between an MPS trial and arbitrary Slater determinants up to a multiplicative error, a crucial subroutine in MPS-AFQMC, is proven to be #P-hard. Nonetheless, we tested several promising heuristics in successfully improving fermionic phaseless AFQMC energies. We also proposed a way to evaluate local energy and force bias evaluations free of matrix product operators. This allows for larger basis set calculations without significant overhead. We showcase the utility of our approach on one- and two-dimensional hydrogen lattices, even when the MPS trial itself struggles to obtain high accuracy. Our work offers a new set of tools that can solve currently challenging electronic structure problems with future improvements.

arXiv:2405.05440•Physical Review Research•Jan 10, 2025•18 citations

Combination of tensor network states and Green's function Monte Carlo

Mingpu Qin

We propose an approach to study the ground state of quantum many-body systems in which tensor network states, specifically projected entangled pair states (PEPSs), and Green's function Monte Carlo (GFMC) are combined. PEPSs, by design, encode the area law which governs the scaling of entanglement entropy in quantum systems with short-range interactions but are hindered by the high computational complexity scaling with bond dimension $D$. GFMC is a highly efficient method, but it usually suffers from the infamous negative sign problem which can be avoided by the fixed-node approximation in which a guiding wave function is utilized to modify the sampling process. The trade-off for the absence of the negative sign problem is the introduction of systematic error by the guiding wave function. In this work, we combine these two methods, PEPS and GFMC, to take advantage of both of them. PEPSs are very accurate variational wave functions, while at the same time, only contractions of the single-layer tensor network are necessary in GFMC, which reduces the cost substantially. Moreover, energy obtained in GFMC is guaranteed to be variational and lower than the variational energy of the guiding PEPS wave function. Benchmark results of the ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model on the square lattice are provided.

arXiv:2006.15608•Physical review. B./Physical review. B•Sep 24, 2020•6 citations

Implementing advanced trial wave functions in fermion quantum Monte Carlo via stochastic sampling

Zhi-Yu Xiao, Tao Xiang, Zixiang Lu +2 more

We introduce an efficient approach to implement correlated many-body trial wave functions in auxiliary-field quantum Monte Carlo (AFQMC). To control the sign/phase problem in AFQMC, a constraint is derived from an exact gauge condition but is typically imposed approximately through a trial wave function or trial density matrix, whose quality can affect the accuracy of the method. Furthermore, the trial wave function can also affect the efficiency through importance sampling. The most natural form of the trial wave function has been single Slater determinants or their linear combinations. More sophisticated forms, for example, with the inclusion of a Jastrow factor or other explicit correlations, have been challenging to use, and their implementation is often assumed to require a quantum computer. In this work, we demonstrate that a large class of correlated wave functions, written in the general form of multi-dimensional integrals over hidden or auxiliary variables times Slater determinants, can be implemented as a trial wave function by coupling the random walkers to a generalized Metropolis sampling. We discuss the fidelity of AFQMC with stochastically sampled trial wave functions, which are relevant to both quantum and classical algorithms. We illustrate the method and show that an efficient implementation can be achieved, which preserves the computational scaling of AFQMC. We test our method in molecules under bond stretching and in transition metal diatomics. Significant improvements are seen in both accuracy and efficiency over typical trial wave functions, and the method yields total ground-state energies systematically within chemical accuracy. The method can be useful for incorporating other advanced wave functions, for example, neural quantum state wave functions optimized from machine learning techniques, or for other forms of fermion quantum Monte Carlo.

arXiv:2505.18519•The Journal of Chemical Physics•Oct 22, 2025•5 citations

NNQS-AFQMC: Neural Network Quantum States Enhanced Fermionic Quantum Monte Carlo

Zhi-Yu Xiao, Bowen Kan, Huan Ma +2 more

We introduce an efficient approach to implement neural network quantum states (NNQS) as trial wave functions in auxiliary-field quantum Monte Carlo (AFQMC). NNQS are a recently developed class of variational ansätze capable of flexibly representing many-body wave functions, though they often incur a high computational cost during optimization. AFQMC, on the other hand, is a powerful stochastic projector approach for ground-state calculations, but it normally requires an approximate constraint via a trial wave function or trial density matrix, whose quality affects the accuracy. Recently, it has been shown (Xiao et al., arXiv2505.18519) that a broad class of highly correlated wave functions can be integrated into AFQMC through stochastic sampling techniques. In this work, we apply this approach and present a direct integration of NNQS with AFQMC, allowing NNQS to serve as high-quality trial wave functions for AFQMC with manageable computational cost. We test the NNQS-AFQMC method on the challenging nitrogen molecule (N 2 ) at stretched geometries. Our results demonstrate that AFQMC with an NNQS trial wave function can attain near-exact total energies, highlighting the potential of AFQMC with NNQS to overcome longstanding challenges in strongly correlated electronic structure calculations. We also outline future research directions for improving this promising methodology.

arXiv:2507.07540•Journal of Chemical Theory and Computation•Oct 1, 2025•2 citations

Green's function Monte Carlo combined with projected entangled pair state approach to the frustrated J1-J2 Heisenberg model

Lin, He-Yu, Yibin Guo, Rong-Qiang He +2 more

The tensor network algorithm, a family of prevalent numerical methods for quantum many-body problems, aptly captures the entanglement properties intrinsic to quantum systems, enabling precise representation of quantum states. However, its computational cost is notably high, particularly in calculating physical observables like correlation functions. To surmount the computational challenge and enhance efficiency, we propose integrating the Green's function Monte Carlo (GFMC) method with the projected entangled pair state (PEPS) ansatz. This approach combines the high-efficiency characteristics of Monte Carlo with the sign-free nature of tensor network states and proves effective in addressing the computational bottleneck. To showcase its prowess, we apply this hybrid approach to investigate the antiferromagnetic J1-J2 Heisenberg model on the square lattice, a model notorious for its sign problem in quantum Monte Carlo simulations. Our results reveal a substantial improvement in the accuracy of ground-state energy when utilizing a preliminary PEPS as the guiding wave function for GFMC. By calculating the structure factor and spin-spin correlation functions, we further characterize the phase diagram, identifying a possible columnar valence-bond state phase within the intermediate parameter range of 0.52<J2/J1<0.58. This comprehensive study underscores the efficacy of our combined approach, demonstrating its ability to accurately simulate frustrated quantum spin systems while ensuring computational efficiency.

•Jan 1, 2024•1 citation

Matrix Product States for Lattice Field Theories

Mari Carmen Bañuls, Krzysztof Cichy, J. I. Cirac +2 more

The term Tensor Network States (TNS) refers to a number of families of states that represent different ansätze for the efficient description of the state of a quantum many-body system.Matrix Product States (MPS) are one particular case of TNS, and have become the most precise tool for the numerical study of one dimensional quantum many-body systems, as the basis of the Density Matrix Renormalization Group method.Lattice Gauge Theories (LGT), in their Hamiltonian version, offer a challenging scenario for these techniques.While the dimensions and sizes of the systems amenable to TNS studies are still far from those achievable by 4-dimensional LGT tools, Tensor Networks can be readily used for problems which more standard techniques, such as Markov chain Monte Carlo simulations, cannot easily tackle.Examples of such problems are the presence of a chemical potential or out-of-equilibrium dynamics.We have explored the performance of Matrix Product States in the case of the Schwinger model, as a widely used testbench for lattice techniques.Using finite-size, open boundary MPS, we are able to determine the low energy states of the model in a fully non-perturbative manner.The precision achieved by the method allows for accurate finite size and continuum limit extrapolations of the ground state energy, but also of the chiral condensate and the mass gaps, thus showing the feasibility of these techniques for gauge theory problems.

arXiv:1310.4118•International Symposium on Lattice Field Theory (LATTICE)•Apr 28, 2014•31 citations
We marry tensor network states (TNS) and projector quantum Monte Carlo (PMC) to overcome the high computational scaling of TNS and the sign problem of PMC. Using TNS as trial wave functions provides a route to systematically improve the sign structure and to eliminate the bias in fixed-node and constrained-path PMC. As a specific example, we describe phaseless auxiliary-field quantum Monte Carlo with matrix product states (MPS-AFQMC). MPS-AFQMC improves significantly on the density-matrix renormalization group (DMRG) ground-state energy. For the J1\ensuremath−J2{J}_{1}\ensuremath{-}{J}_{2}J1​\ensuremath−J2​ model on two-dimensional square lattices, we observe with MPS-AFQMC an order of magnitude reduction in the error for all couplings, compared to DMRG. The improvement is independent of walker bond dimension, and we therefore use bond dimension 1 for the walkers. The computational cost of MPS-AFQMC is then quadratic in the bond dimension of the trial wave function, which is lower than the cubic scaling of DMRG. The error due to the constrained-path bias is proportional to the variational error of the trial wave function. We show that for the J1\ensuremath−J2{J}_{1}\ensuremath{-}{J}_{2}J1​\ensuremath−J2​ model on two-dimensional square lattices, a linear extrapolation of the MPS-AFQMC energy with the discarded weight from the DMRG calculation allows us to remove the constrained-path bias. Extensions to other tensor networks are briefly discussed.
In this work, we report, for the first time, an implementation of fermionic auxiliary-field quantum Monte Carlo (AFQMC) using matrix product state (MPS) trial wavefunctions, dubbed MPS-AFQMC. Calculating overlaps between an MPS trial and arbitrary Slater determinants up to a multiplicative error, a crucial subroutine in MPS-AFQMC, is proven to be #P-hard. Nonetheless, we tested several promising heuristics in successfully improving fermionic phaseless AFQMC energies. We also proposed a way to evaluate local energy and force bias evaluations free of matrix product operators. This allows for larger basis set calculations without significant overhead. We showcase the utility of our approach on one- and two-dimensional hydrogen lattices, even when the MPS trial itself struggles to obtain high accuracy. Our work offers a new set of tools that can solve currently challenging electronic structure problems with future improvements.
We propose an approach to study the ground state of quantum many-body systems in which tensor network states, specifically projected entangled pair states (PEPSs), and Green's function Monte Carlo (GFMC) are combined. PEPSs, by design, encode the area law which governs the scaling of entanglement entropy in quantum systems with short-range interactions but are hindered by the high computational complexity scaling with bond dimension DDD. GFMC is a highly efficient method, but it usually suffers from the infamous negative sign problem which can be avoided by the fixed-node approximation in which a guiding wave function is utilized to modify the sampling process. The trade-off for the absence of the negative sign problem is the introduction of systematic error by the guiding wave function. In this work, we combine these two methods, PEPS and GFMC, to take advantage of both of them. PEPSs are very accurate variational wave functions, while at the same time, only contractions of the single-layer tensor network are necessary in GFMC, which reduces the cost substantially. Moreover, energy obtained in GFMC is guaranteed to be variational and lower than the variational energy of the guiding PEPS wave function. Benchmark results of the ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model on the square lattice are provided.
We introduce an efficient approach to implement correlated many-body trial wave functions in auxiliary-field quantum Monte Carlo (AFQMC). To control the sign/phase problem in AFQMC, a constraint is derived from an exact gauge condition but is typically imposed approximately through a trial wave function or trial density matrix, whose quality can affect the accuracy of the method. Furthermore, the trial wave function can also affect the efficiency through importance sampling. The most natural form of the trial wave function has been single Slater determinants or their linear combinations. More sophisticated forms, for example, with the inclusion of a Jastrow factor or other explicit correlations, have been challenging to use, and their implementation is often assumed to require a quantum computer. In this work, we demonstrate that a large class of correlated wave functions, written in the general form of multi-dimensional integrals over hidden or auxiliary variables times Slater determinants, can be implemented as a trial wave function by coupling the random walkers to a generalized Metropolis sampling. We discuss the fidelity of AFQMC with stochastically sampled trial wave functions, which are relevant to both quantum and classical algorithms. We illustrate the method and show that an efficient implementation can be achieved, which preserves the computational scaling of AFQMC. We test our method in molecules under bond stretching and in transition metal diatomics. Significant improvements are seen in both accuracy and efficiency over typical trial wave functions, and the method yields total ground-state energies systematically within chemical accuracy. The method can be useful for incorporating other advanced wave functions, for example, neural quantum state wave functions optimized from machine learning techniques, or for other forms of fermion quantum Monte Carlo.
We introduce an efficient approach to implement neural network quantum states (NNQS) as trial wave functions in auxiliary-field quantum Monte Carlo (AFQMC). NNQS are a recently developed class of variational ansätze capable of flexibly representing many-body wave functions, though they often incur a high computational cost during optimization. AFQMC, on the other hand, is a powerful stochastic projector approach for ground-state calculations, but it normally requires an approximate constraint via a trial wave function or trial density matrix, whose quality affects the accuracy. Recently, it has been shown (Xiao et al., arXiv2505.18519) that a broad class of highly correlated wave functions can be integrated into AFQMC through stochastic sampling techniques. In this work, we apply this approach and present a direct integration of NNQS with AFQMC, allowing NNQS to serve as high-quality trial wave functions for AFQMC with manageable computational cost. We test the NNQS-AFQMC method on the challenging nitrogen molecule (N 2 ) at stretched geometries. Our results demonstrate that AFQMC with an NNQS trial wave function can attain near-exact total energies, highlighting the potential of AFQMC with NNQS to overcome longstanding challenges in strongly correlated electronic structure calculations. We also outline future research directions for improving this promising methodology.
The tensor network algorithm, a family of prevalent numerical methods for quantum many-body problems, aptly captures the entanglement properties intrinsic to quantum systems, enabling precise representation of quantum states. However, its computational cost is notably high, particularly in calculating physical observables like correlation functions. To surmount the computational challenge and enhance efficiency, we propose integrating the Green's function Monte Carlo (GFMC) method with the projected entangled pair state (PEPS) ansatz. This approach combines the high-efficiency characteristics of Monte Carlo with the sign-free nature of tensor network states and proves effective in addressing the computational bottleneck. To showcase its prowess, we apply this hybrid approach to investigate the antiferromagnetic J1-J2 Heisenberg model on the square lattice, a model notorious for its sign problem in quantum Monte Carlo simulations. Our results reveal a substantial improvement in the accuracy of ground-state energy when utilizing a preliminary PEPS as the guiding wave function for GFMC. By calculating the structure factor and spin-spin correlation functions, we further characterize the phase diagram, identifying a possible columnar valence-bond state phase within the intermediate parameter range of 0.52<J2/J1<0.58. This comprehensive study underscores the efficacy of our combined approach, demonstrating its ability to accurately simulate frustrated quantum spin systems while ensuring computational efficiency.
The term Tensor Network States (TNS) refers to a number of families of states that represent different ansätze for the efficient description of the state of a quantum many-body system.Matrix Product States (MPS) are one particular case of TNS, and have become the most precise tool for the numerical study of one dimensional quantum many-body systems, as the basis of the Density Matrix Renormalization Group method.Lattice Gauge Theories (LGT), in their Hamiltonian version, offer a challenging scenario for these techniques.While the dimensions and sizes of the systems amenable to TNS studies are still far from those achievable by 4-dimensional LGT tools, Tensor Networks can be readily used for problems which more standard techniques, such as Markov chain Monte Carlo simulations, cannot easily tackle.Examples of such problems are the presence of a chemical potential or out-of-equilibrium dynamics.We have explored the performance of Matrix Product States in the case of the Schwinger model, as a widely used testbench for lattice techniques.Using finite-size, open boundary MPS, we are able to determine the low energy states of the model in a fully non-perturbative manner.The precision achieved by the method allows for accurate finite size and continuum limit extrapolations of the ground state energy, but also of the chiral condensate and the mass gaps, thus showing the feasibility of these techniques for gauge theory problems.