Abstract:The study of ground-state properties of the Fermi-Hubbard model is a long-lasting task in the research of strongly correlated systems. Owing to the exponentially growing complexity of the system, a quantitative analysis usually demands high computational cost and is restricted to small samples, especially in two or higher dimensions. Here, we introduce a variational method in the frame of fermionic Gaussian states, and obtain the ground states of one- and two-dimensional attractive Hubbard models via imaginary-time evolution. We calculate the total energy and benchmark the results in a wide range of interaction strength and filling factor with those obtained via exact two-body results, the density matrix renormalization group based on matrix product states (MPS), and projector Quantum Monte Carlo (QMC) method. For both 1D and 2D cases, the Gaussian variational method presents accurate results for total energy with a maximum systematic error ~4% in the intermediate interaction region. The accuracy of these results has negligible dependence on the system size. We further calculate the double occupancy and find excellent agreement with MPS and QMC, as well as the experimental results of cold quantum gases in optical lattices. The results suggest that the Gaussian pairing state is a good approximation to the ground states of attractive Hubbard model, in particular in the strong and weak coupling limits. Moreover, we generalize the method to the attractive Hubbard model with a finite spin-polarization, which can be mapped to the repulsive interaction case via particle-hole transformation, and obtain accurate results of ground state energy and double occupancy. Our work demonstrates the ability of the Gaussian variational method to extract ground state properties of strongly correlated many-body systems with negligible computational cost, especially of large size and in higher dimensions.

Generalized Hartree Fock Theory for Dispersion Relations of Interacting Fermionic Lattice Systems

Gaussian variational method to Fermi Hubbard model in one and two dimensions

General Hubbard Model for Strongly Interacting Fermions in an Optical Lattice and Its Phase Detection

Density functional theory in one-dimension for contact-interacting fermions

Effective Single-Band Models for Strongly Interacting Fermions in an Optical Lattice

New results for a photon-photon collider

Operator Valued Flow Equation Approach to the Bosonic Lattice Polaron: Dispersion Renormalization Beyond the Fröhlich Paradigm

Superfluidity of Fermions with Repulsive On-Site Interaction in an Anisotropic Optical Lattice Near A Feshbach Resonance

Extendibility of fermionic states and rigorous ground state approximations of interacting fermionic systems

Interacting chiral fermions on the lattice with matrix product operator norms

A generalized effective spin-chain formalism for strongly interacting spinor gases in optical lattice

A Variational Improvement of the Hartree-Fock Approach to the 2D Hubbard Model

Hartree-Fock-Bogoliubov Theory of Dipolar Fermi Gases

The Fractal-Lattice Hubbard Model

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

Functional renormalization group for fermions on a one dimensional lattice at arbitrary filling

Grassmann phase space dynamics of strongly-correlated fermion

A multi-layer multi-configurational time-dependent Hartree approach to lattice models beyond one dimension

Hartree-Fock Treatment of Fermi Polarons Using the Lee-Low-Pine Transformation

What can lattice DFT teach us about real-space DFT?