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/0007304 · Journal of the Physical Society of Japan

Path-Integral Renormalization Group Method for Numerical Study of Strongly Correlated Electron Systems

Masatoshi Imada, Tsuyoshi Kashima

Published Sep 1, 2000• 78 citations
arXivHTMLPDF

Abstract

A numerical algorithm for studying strongly correlated electron systems is proposed. The groundstate wavefunction is projected out after numerical renormalization procedure in the path integral formalism. The wavefunction is expressed from the optimized linear combination of retained states in the truncated Hilbert space with numerically chosen basis. This algorithm does not suffer from the negative sign problem and can be applied to any type of Hamiltonian in any dimension. The efficiency is tested in examples of the Hubbard model where the basis of Slater determinants is numerically optimized. We show results on fast convergence and accuracy achieved with small number of retained states.

DOI
10.1143/jpsj.69.2723

Related papers

Path-Integral Renormalization Group Method for Numerical Study on Ground States of Strongly Correlated Electronic Systems

Tsuyoshi Kashima, Masatoshi Imada

A new efficient numerical algorithm for interacting fermion systems is proposed and examined in detail. The ground state is expressed approximately by a linear combination of numerically chosen basis states in a truncated Hilbert space. Two procedures lead to a better approximation. The first is a numerical renormalization, which optimizes the chosen basis and projects onto the ground state within the fixed dimension, $L$, of the Hilbert space. The second is an increase of the dimension of the truncated Hilbert space, which enables the linear combination to converge to a better approximation. The extrapolation $L\to\infty$ after the convergence removes the approximation error systematically. This algorithm does not suffer from the negative sign problem and can be applied to systems in any spatial dimension and arbitrary lattice structure. The efficiency is tested and the implementation explained for two-dimensional Hubbard models where Slater determinants are employed as chosen basis. Our results with less than 400 chosen basis indicate good accuracy within the errorbar of the best available results as those of the quantum Monte Carlo for energy and other physical quantities.

arXiv:cond-mat/0104140•Journal of the Physical Society of Japan•Aug 15, 2001•89 citations

First-principles path-integral renormalization-group method for Coulombic many-body systems

Masashi Kojo, Kikuji Hirose

An approach for obtaining the ground state of Coulombic many-body systems is presented. This approach is based on the path-integral renormalization-group method with nonorthogonal Slater determinants, is free of the negative sign problem, and can handle higher dimensional systems with consideration of the correlation effect. Furthermore, it can be easily extended to the multicomponent quantum systems that contain more than two kinds of quantum particles. According to our results obtained with the present approach, it achieves the same accuracy as the variational Monte Carlo method with a few Slater determinants and enables us to study the entire ground state consisting of electrons and nuclei without the need to use the Born-Oppenheimer approximation.

Physical Review A•Oct 23, 2009•4 citations

Finite-Temperature Variational Monte Carlo Method for Strongly Correlated Electron Systems

Kensaku Takai, Kota Ido, Takahiro Misawa +2 more

A new computational method for finite-temperature properties of strongly correlated electrons is proposed by extending the variational Monte Carlo method originally developed for the ground state. The method is based on the path integral in the imaginary-time formulation, starting from the infinite-temperature state that is well approximated by a small number of certain random initial states. Lower temperatures are progressively reached by the imaginary-time evolution. The algorithm follows the framework of the quantum transfer matrix and finite-temperature Lanczos methods, but we extends them to treat much larger system sizes without the negative sign problem by optimizing the truncated Hilbert space on the basis of the time-dependent variational principle (TDVP). This optimization algorithm is equivalent to the stochastic reconfiguration (SR) method that has been frequently used for the ground state to optimally truncate the Hilbert space. The obtained finite-temperature states allow an interpretation based on the thermal pure quantum (TPQ) state instead of the conventional canonical-ensemble average. Our method is tested for the one- and two-dimensional Hubbard models and its accuracy and efficiency are demonstrated.

arXiv:1510.05352•Journal of the Physical Society of Japan•Feb 18, 2016•21 citations

Calculating ground-state properties of correlated fermionic systems with BCS trial wave functions in Slater determinant path-integral approaches

Ettore Vitali, Peter Rosenberg, Shiwei Zhang

We introduce an efficient and numerically stable technique to make use of a BCS trial wave function in the computation of correlation functions of strongly correlated quantum fermion systems. The technique is applicable to any projection approach involving paths of independent-fermion propagators, for example, in mean-field or auxiliary-field quantum Monte Carlo (AFQMC) calculations. Within AFQMC, in the absence of the sign problem, the methodology allows the use of a BCS reference state which can greatly reduce the required imaginary time of projection and improves Monte Carlo sampling efficiency and statistical accuracy for systems where pairing correlations are important. When the sign problem is present, the approach provides a powerful generalization of the constrained-path AFQMC technique which usually uses Slater determinant trial wave functions. As a demonstration of the capability of the methodology, we present benchmark results for the attractive Hubbard model, both spin balanced (no sign problem) and with a finite spin polarization (with sign problem).

arXiv:1905.05012•Physical Review A•Aug 23, 2019•14 citations

Path-Integral Renormalization Group Method

Masatoshi Imada

Path‐integral renormalization group (PIRG) method has been developed for studying strongly correlated electron systems. In a wide range of systems including Hubbard‐type models, this method overcomes a number of difficulties known in the Monte‐Carlo‐type methods such as the negative sign problem. This method has been combined with a procedure of quantum number projection and grand canonical ensemble method as well, which contribute to a wider applicability. The quantum‐number projected PIRG enables calculations of excited spectra with a specified momentum or spin. We review recent numerical results on strongly correlated electron systems studied by this method. By using the methods, we determine the phase diagram of the two‐dimensional Hubbard model in the parameter space of the onsite interaction U, strength of the geometrical frustration effects defined by the ratio between next‐nearest to nearest neighbor transfer integrals, t′/t, and the chemical potential. It reveals severe competitions of various phases. The phase diagram contains a remarkable nonmagnetic Mott insulator phase with gapless and dispersionless spin excitations, sandwitched by the first‐order Mott transition and the antiferromagnetic transition. The first‐order character becomes more continuous one with increasing the frustration effect.

AIP conference proceedings•Jan 1, 2003•1 citation

Quantum-number projection in the path-integral renormalization group method

T. Mizusaki, Masatoshi Imada

We present a quantum-number projection technique which enables us to exactly treat spin, momentum, and other symmetries embedded in the Hubbard model. By combining this projection technique, we extend the path-integral renormalization-group method to improve the efficiency of numerical computations. By taking numerical calculations for the standard Hubbard model and the Hubbard model with next-nearest-neighbor transfer, we show that the present extended method can extremely enhance numerical accuracy and that it can handle excited states, in addition to the ground state.

arXiv:cond-mat/0311005•Physical Review B•Mar 23, 2004•69 citations
A numerical algorithm for studying strongly correlated electron systems is proposed. The groundstate wavefunction is projected out after numerical renormalization procedure in the path integral formalism. The wavefunction is expressed from the optimized linear combination of retained states in the truncated Hilbert space with numerically chosen basis. This algorithm does not suffer from the negative sign problem and can be applied to any type of Hamiltonian in any dimension. The efficiency is tested in examples of the Hubbard model where the basis of Slater determinants is numerically optimized. We show results on fast convergence and accuracy achieved with small number of retained states.
A new efficient numerical algorithm for interacting fermion systems is proposed and examined in detail. The ground state is expressed approximately by a linear combination of numerically chosen basis states in a truncated Hilbert space. Two procedures lead to a better approximation. The first is a numerical renormalization, which optimizes the chosen basis and projects onto the ground state within the fixed dimension, LLL, of the Hilbert space. The second is an increase of the dimension of the truncated Hilbert space, which enables the linear combination to converge to a better approximation. The extrapolation L→∞L\to\inftyL→∞ after the convergence removes the approximation error systematically. This algorithm does not suffer from the negative sign problem and can be applied to systems in any spatial dimension and arbitrary lattice structure. The efficiency is tested and the implementation explained for two-dimensional Hubbard models where Slater determinants are employed as chosen basis. Our results with less than 400 chosen basis indicate good accuracy within the errorbar of the best available results as those of the quantum Monte Carlo for energy and other physical quantities.
An approach for obtaining the ground state of Coulombic many-body systems is presented. This approach is based on the path-integral renormalization-group method with nonorthogonal Slater determinants, is free of the negative sign problem, and can handle higher dimensional systems with consideration of the correlation effect. Furthermore, it can be easily extended to the multicomponent quantum systems that contain more than two kinds of quantum particles. According to our results obtained with the present approach, it achieves the same accuracy as the variational Monte Carlo method with a few Slater determinants and enables us to study the entire ground state consisting of electrons and nuclei without the need to use the Born-Oppenheimer approximation.
A new computational method for finite-temperature properties of strongly correlated electrons is proposed by extending the variational Monte Carlo method originally developed for the ground state. The method is based on the path integral in the imaginary-time formulation, starting from the infinite-temperature state that is well approximated by a small number of certain random initial states. Lower temperatures are progressively reached by the imaginary-time evolution. The algorithm follows the framework of the quantum transfer matrix and finite-temperature Lanczos methods, but we extends them to treat much larger system sizes without the negative sign problem by optimizing the truncated Hilbert space on the basis of the time-dependent variational principle (TDVP). This optimization algorithm is equivalent to the stochastic reconfiguration (SR) method that has been frequently used for the ground state to optimally truncate the Hilbert space. The obtained finite-temperature states allow an interpretation based on the thermal pure quantum (TPQ) state instead of the conventional canonical-ensemble average. Our method is tested for the one- and two-dimensional Hubbard models and its accuracy and efficiency are demonstrated.
We introduce an efficient and numerically stable technique to make use of a BCS trial wave function in the computation of correlation functions of strongly correlated quantum fermion systems. The technique is applicable to any projection approach involving paths of independent-fermion propagators, for example, in mean-field or auxiliary-field quantum Monte Carlo (AFQMC) calculations. Within AFQMC, in the absence of the sign problem, the methodology allows the use of a BCS reference state which can greatly reduce the required imaginary time of projection and improves Monte Carlo sampling efficiency and statistical accuracy for systems where pairing correlations are important. When the sign problem is present, the approach provides a powerful generalization of the constrained-path AFQMC technique which usually uses Slater determinant trial wave functions. As a demonstration of the capability of the methodology, we present benchmark results for the attractive Hubbard model, both spin balanced (no sign problem) and with a finite spin polarization (with sign problem).
Path‐integral renormalization group (PIRG) method has been developed for studying strongly correlated electron systems. In a wide range of systems including Hubbard‐type models, this method overcomes a number of difficulties known in the Monte‐Carlo‐type methods such as the negative sign problem. This method has been combined with a procedure of quantum number projection and grand canonical ensemble method as well, which contribute to a wider applicability. The quantum‐number projected PIRG enables calculations of excited spectra with a specified momentum or spin. We review recent numerical results on strongly correlated electron systems studied by this method. By using the methods, we determine the phase diagram of the two‐dimensional Hubbard model in the parameter space of the onsite interaction U, strength of the geometrical frustration effects defined by the ratio between next‐nearest to nearest neighbor transfer integrals, t′/t, and the chemical potential. It reveals severe competitions of various phases. The phase diagram contains a remarkable nonmagnetic Mott insulator phase with gapless and dispersionless spin excitations, sandwitched by the first‐order Mott transition and the antiferromagnetic transition. The first‐order character becomes more continuous one with increasing the frustration effect.
We present a quantum-number projection technique which enables us to exactly treat spin, momentum, and other symmetries embedded in the Hubbard model. By combining this projection technique, we extend the path-integral renormalization-group method to improve the efficiency of numerical computations. By taking numerical calculations for the standard Hubbard model and the Hubbard model with next-nearest-neighbor transfer, we show that the present extended method can extremely enhance numerical accuracy and that it can handle excited states, in addition to the ground state.