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.

A gauge constrained algorithm of VDAT at $\mathcal{N}=3$ for the multi-orbital Hubbard model

Precise ground state of multi-orbital Mott systems via the variational discrete action theory

Benchmark Study of the Two-Dimensional Hubbard Model with Auxiliary-Field Quantum Monte Carlo Method

Machine learning Hubbard parameters with equivariant neural networks

Improved Multi-Variable Variational Monte Carlo Method Examined by High-Precision Calculations of One-Dimensional Hubbard Model

Tangent Space Approach for Thermal Tensor Network Simulations of the 2D Hubbard Model

Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

The Hubbard model on the Bethe lattice via variational uniform tree states: metal-insulator transition and a Fermi liquid

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

Variational Monte Carlo Study of the Spiral State in the Hubbard Model.

Determinant Quantum Monte Carlo for the Half-Filled Hubbard Model with Nonlocal Density-Density Interactions

New Approach to Strong Correlation: Twisting Hubbard into the Orbital Hatsugai-Kohmoto Model

Anisotropic multi-orbital Hubbard model simulated with impurity approximation

Boosting Determinant Quantum Monte Carlo with Submatrix Updates: Unveiling the Phase Diagram of the 3D Hubbard Model

Complete phase diagram for three-band Hubbard model with orbital degeneracy lifted by crystal field splitting

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

Two-dimensional extended Hubbard model: doping, next-nearest neighbor hopping and phase diagrams

Momentum Distribution and Tunneling Density of States of One-Dimensional Fermionic SU(N) Hubbard Model

Extended Hubbard model on a C$_{20}$ molecule