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

Variational Monte Carlo simulations using tensor-product projected states

Olga Sikora, Hsueh-Wen Chang, Chung-Pin Chou, Frank Pollmann, Ying-Jer Kao

Published Apr 8, 2015• 18 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.91.165113

Related papers

Matrix-product-based projected wave functions ansatz for quantum many-body ground states

Chung-Pin Chou, Frank Pollmann, Ting-Kuo Lee

We introduce a projected wave function approach based on projection operators in the form of matrix-product operators (MPOs). Our approach allows us to variationally improve the short-range entanglement of a given trial wave function by optimizing the matrix elements of the MPOs while the long-range entanglement is contained in the initial guess of the wave function. The optimization is performed using standard variational Monte Carlo techniques. We demonstrate the efficiency of our approach by considering a one-dimensional model of interacting spinless fermions. In addition, we indicate how to generalize this approach to higher dimensions using projection operators which are based on tensor products.

arXiv:1201.6121•Physical Review B•Jul 19, 2012•14 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

Matrix product states with backflow correlations

Guglielmo Lami, Giuseppe Carleo, Mario Collura

By taking inspiration from the backflow transformation for correlated systems, we introduce a tensor network Ansatz which extends the well-established matrix product state representation of a quantum many-body wave function. This structure provides enough resources to ensure that states in dimensions larger than or equal to one obey an area law for entanglement. It can be efficiently manipulated to address the ground-state search problem by means of an optimization scheme which mixes tensor-network and variational Monte Carlo algorithms. We benchmark the Ansatz against spin models both in one and two dimensions, demonstrating high accuracy and precision. We finally employ our approach to study the challenging $S=1/2$ two-dimensional (2D) ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model, demonstrating that it is competitive with the state-of-the-art methods in 2D.

arXiv:2201.00810•Physical review. B./Physical review. B•Aug 10, 2022•14 citations

Tensor Network Representations of Parton Wave Functions

Ying-Hai Wu, Lei Wang, Hong-Hao Tu

Tensor network states and parton wave functions are two pivotal methods for studying quantum many-body systems. This work connects these two subjects as we demonstrate that a variety of parton wave functions, such as projected Fermi sea and projected fermionic or bosonic paired states, can be represented exactly as tensor networks. The results can be compressed into matrix product states with moderate bond dimensions so various physical quantities can be computed efficiently. For the projected Fermi sea, we develop an excellent compression scheme with high fidelity using maximally localized Wannier orbitals. Numerical calculations on two parton wave functions demonstrate that our method exceeds commonly adopted Monte Carlo methods in some aspects. It produces energy and correlation function with very high accuracy that is difficult to achieve using Monte Carlo method. The entanglement measures that were almost impossible to compute before can also be obtained easily using our method.

arXiv:1910.11011•Physical Review Letters•Jun 19, 2020•39 citations

Variational Monte Carlo method for fermionic models combined with tensor networks and applications to the hole-doped two-dimensional Hubbard model

Hui‐Hai Zhao, Kota Ido, Satoshi Morita +1 more

The conventional tensor-network states employ real-space product states as reference wave functions. Here, we propose a many-variable variational Monte Carlo (mVMC) method combined with tensor networks by taking advantages of both to study fermionic models. The variational wave function is composed of a pair product wave function operated by real-space correlation factors and tensor networks. Moreover, we can apply quantum number projections, such as spin, momentum, and lattice symmetry projections, to recover the symmetry of the wave function to further improve the accuracy. We benchmark our method for one- and two-dimensional Hubbard models, which show significant improvement over the results obtained individually either by mVMC or by tensor network. We have applied the present method to a hole-doped Hubbard model on the square lattice, which indicates the stripe charge/spin order coexisting with a weak $d$-wave superconducting order in the ground state for the doping concentration of less than 0.3, where the stripe oscillation period gets longer with increasing hole concentration. The charge homogeneous and highly superconducting state also exists as a metastable excited state for the doping concentration less than 0.25.

arXiv:1703.03537•Physical review. B./Physical review. B•Aug 1, 2017•42 citations

Monte Carlo simulations of two-dimensional fermion systems with string-bond states

Jeong-Pil Song, Rudolf Torsten Clay

We describe an application of variational Monte Carlo to two-dimensional fermionic systems within the recently developed tensor-network string-bond state ansatz. We use a combination of variational Monte Carlo and stochastic optimization to optimize the matrix-product state matrices representing the ground state. We present results for a two-dimensional spinless fermion model including nearest-neighbor Coulomb interactions and determine using finite-size scaling the phase boundary between charge-ordered insulating and metallic phases. This approach can treat frustrated systems and be easily extended to fermion models with spin.

arXiv:1305.4151•Physical Review B•Feb 3, 2014•8 citations
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.
We introduce a projected wave function approach based on projection operators in the form of matrix-product operators (MPOs). Our approach allows us to variationally improve the short-range entanglement of a given trial wave function by optimizing the matrix elements of the MPOs while the long-range entanglement is contained in the initial guess of the wave function. The optimization is performed using standard variational Monte Carlo techniques. We demonstrate the efficiency of our approach by considering a one-dimensional model of interacting spinless fermions. In addition, we indicate how to generalize this approach to higher dimensions using projection operators which are based on tensor products.
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.
By taking inspiration from the backflow transformation for correlated systems, we introduce a tensor network Ansatz which extends the well-established matrix product state representation of a quantum many-body wave function. This structure provides enough resources to ensure that states in dimensions larger than or equal to one obey an area law for entanglement. It can be efficiently manipulated to address the ground-state search problem by means of an optimization scheme which mixes tensor-network and variational Monte Carlo algorithms. We benchmark the Ansatz against spin models both in one and two dimensions, demonstrating high accuracy and precision. We finally employ our approach to study the challenging S=1/2S=1/2S=1/2 two-dimensional (2D) ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model, demonstrating that it is competitive with the state-of-the-art methods in 2D.
Tensor network states and parton wave functions are two pivotal methods for studying quantum many-body systems. This work connects these two subjects as we demonstrate that a variety of parton wave functions, such as projected Fermi sea and projected fermionic or bosonic paired states, can be represented exactly as tensor networks. The results can be compressed into matrix product states with moderate bond dimensions so various physical quantities can be computed efficiently. For the projected Fermi sea, we develop an excellent compression scheme with high fidelity using maximally localized Wannier orbitals. Numerical calculations on two parton wave functions demonstrate that our method exceeds commonly adopted Monte Carlo methods in some aspects. It produces energy and correlation function with very high accuracy that is difficult to achieve using Monte Carlo method. The entanglement measures that were almost impossible to compute before can also be obtained easily using our method.
The conventional tensor-network states employ real-space product states as reference wave functions. Here, we propose a many-variable variational Monte Carlo (mVMC) method combined with tensor networks by taking advantages of both to study fermionic models. The variational wave function is composed of a pair product wave function operated by real-space correlation factors and tensor networks. Moreover, we can apply quantum number projections, such as spin, momentum, and lattice symmetry projections, to recover the symmetry of the wave function to further improve the accuracy. We benchmark our method for one- and two-dimensional Hubbard models, which show significant improvement over the results obtained individually either by mVMC or by tensor network. We have applied the present method to a hole-doped Hubbard model on the square lattice, which indicates the stripe charge/spin order coexisting with a weak ddd-wave superconducting order in the ground state for the doping concentration of less than 0.3, where the stripe oscillation period gets longer with increasing hole concentration. The charge homogeneous and highly superconducting state also exists as a metastable excited state for the doping concentration less than 0.25.
We describe an application of variational Monte Carlo to two-dimensional fermionic systems within the recently developed tensor-network string-bond state ansatz. We use a combination of variational Monte Carlo and stochastic optimization to optimize the matrix-product state matrices representing the ground state. We present results for a two-dimensional spinless fermion model including nearest-neighbor Coulomb interactions and determine using finite-size scaling the phase boundary between charge-ordered insulating and metallic phases. This approach can treat frustrated systems and be easily extended to fermion models with spin.