Abstract: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.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to accurately calculate the ground - state properties in strongly correlated electron systems. Specifically, the paper proposes a new and efficient numerical algorithm - the Path - Integral Renormalization Group method (PIRG) for studying the ground - state properties of strongly correlated fermion systems. This method aims to overcome the negative - sign problem encountered by existing methods such as Quantum Monte Carlo (QMC) when dealing with certain systems, and can be applied to systems of any dimension and lattice structure. ### Main problems 1. **Negative - sign problem**: QMC will encounter the negative - sign problem when dealing with certain strongly correlated electron systems (such as fermion systems and frustrated spin systems), resulting in difficult - to - control statistical and rounding errors in calculation results. 2. **Spatial dimension limitation**: Although the Density Matrix Renormalization Group (DMRG) is not affected by the negative - sign problem, its spatial renormalization process makes it mainly applicable to one - dimensional systems. 3. **Computational efficiency and accuracy**: Existing numerical methods have problems of low computational efficiency and insufficient accuracy when dealing with large - scale systems. ### Solutions The PIRG method solves the above problems through the following steps: 1. **Initial state selection**: Any type of initial state can be selected, and an initial state close to the ground state is usually selected. 2. **Projection operation**: The projection of the ground state is achieved by introducing projection operators $\exp(-\tau H)$ or $1-\tau H$. 3. **Truncation operation**: Since the projection operation will increase the number of ground states, it is necessary to select important states and maintain the dimension $L$ of the stored Hilbert subspace. 4. **Iterative process**: Through the use of appropriate projection operators and truncation operations, numerical renormalization is achieved, so that the stored $L$-dimensional subspace gradually approaches the ground state. This renormalization process needs to be repeated until the lowest energy of the truncated Hilbert subspace converges to the lowest energy. ### Advantages - **No negative - sign problem**: The PIRG method is not affected by the negative - sign problem, so it can be applied to Hamiltonians of any type and dimension. - **Wide range of application**: It can be applied to systems of any dimension and lattice structure. - **High accuracy**: Experimental results show that when dealing with the two - dimensional Hubbard model, the PIRG method can achieve an accuracy comparable to that of the Quantum Monte Carlo method with less than 400 selected ground states. ### Mathematical representation - **Ground - state wave function**: \[ |\psi_g\rangle=\lim_{\tau\rightarrow\infty}\exp[-\tau H]|\phi_0\rangle \] - **Projection process**: \[ |\psi_g\rangle=\lim_{\Delta\tau\rightarrow0}\lim_{n\rightarrow\infty}(\exp[-\Delta\tau H])^n|\phi_0\rangle \] - **Hamiltonian**: \[ H = H_k + H_U \] where, \[ H_k=-\sum_{<i,j>,\sigma}t_{ij}(c^\dagger_{i\sigma}c_{j\sigma}+h.c.) \] \[ H_U = U\sum_i(n_{i\uparrow}-\frac{1}{2})(n_{i\downarrow}-\frac{1}{2})=U\sum_i n_{i\uparrow}n_{i\downarrow}-U\left(\frac{M - N}{4}\right) \] ### Conclusion The PIRG method provides a new and efficient numerical algorithm that can accurately calculate the ground - state properties of strongly correlated electron systems, especially performing well in dealing with the negative - sign problem and multi - dimensional systems. This method

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

Path-Integral Renormalization Group Treatments for Many-Electron Systems with Long-Range Repulsive Interactions

Perturbational method for the numerical evaluation of the ground-state energy of the Heisenberg model

Partial Renormalization of Quasiparticle Interactions

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

The (extended) dynamical mean field theory combined with the two-particle irreducible functional renormalization-group approach as a tool to study strongly-correlated systems

Functional renormalization group without functional integrals: implementing Hilbert space projections for strongly correlated electrons via Hubbard X-operators

Quantum many-body simulations of the two-dimensional Fermi-Hubbard model in ultracold optical lattices

RENORMALIZATION GROUP APPROACHES TO STRONGLY CORRELATED ELECTRON SYSTEMS

A functional renormalization group approach to electronic structure calculations for systems without translational symmetry

Two dimensional quantum lattice models via mode optimized hybrid CPU-GPU density matrix renormalization group method

Quasi-particle functional Renormalisation Group calculations in the two-dimensional t-t'-Hubbard model

Efficient Simulation of 1D Long-Range Interacting Systems at Any Temperature

Covariant Bethe-Salpeter Approximation in Models of Strongly Correlated Electron Systems

Elective-momentum Ultrasize Quantum Monte Carlo Method

From Heisenberg to Hubbard: An initial state for the shallow quantum simulation of correlated electrons

Lower Bounds on Ground-State Energies of Local Hamiltonians through the Renormalization Group

Tensor Renormalization Group for fermions

High-efficiency quantum Monte Carlo algorithm for extracting entanglement entropy in interacting fermion systems

Functional Renormalization Group for fermions on a one dimensional lattice at arbitrary filling