A novel variational approach is outlined for treating electron correlation in many-particle systems and is applied to the one-band Hubbard model. The authors' treatment is based on the concept of fermion pairs and their representation in terms of ideal-boson operators. In this representation approximate eigen-states of the Hamiltonian are obtained in the form of an antisymmetrised product of pair orbitals. For the ground state of the Hubbard model these orbitals are then determined by the use of the variational principle. Because of the complexity of the resulting equation a simplified ansatz is used for the orbitals, which leads to the same energy expression as is obtained by the traditional Bogolubov method. The trial states used in both approaches are, however, quite different. They find reasonable qualitative agreement with well established results for the Hubbard model in one dimension even in those cases where the Bogolubov or related approaches fail. They are also able to show that in the atomic limit the Bogolubov energy bound is the best upper bound that can be obtained by using a single determinant of pair orbitals. Therefore, the only way to improve the ground-state energy in the atomic limit is by means of configuration interaction.
A novel variational approach is outlined for treating electron correlation in many-particle systems and is applied to the one-band Hubbard model. The authors' treatment is based on the concept of fermion pairs and their representation in terms of ideal-boson operators. In this representation approximate eigen-states of the Hamiltonian are obtained in the form of an antisymmetrised product of pair orbitals. For the ground state of the Hubbard model these orbitals are then determined by the use of the variational principle. Because of the complexity of the resulting equation a simplified ansatz is used for the orbitals, which leads to the same energy expression as is obtained by the traditional Bogolubov method. The trial states used in both approaches are, however, quite different. They find reasonable qualitative agreement with well established results for the Hubbard model in one dimension even in those cases where the Bogolubov or related approaches fail. They are also able to show that in the atomic limit the Bogolubov energy bound is the best upper bound that can be obtained by using a single determinant of pair orbitals. Therefore, the only way to improve the ground-state energy in the atomic limit is by means of configuration interaction.
We study the one-dimensional Hubbard model for interacting electrons in a crystal lattice, using two methods: solution of the Lieb-Wu equations in the atomic limit of no intersite hopping, and variational techniques to obtain a lower bound for the ground state energy. From the Lieb-Wu equations we derive an exact expression for the ground state energy in the atomic limit. We obtain our result by exploiting the spin symmetry of the Hubbard Hamiltonian. Our energy is lower than that obtained by previous published results, and we are able to explain why our results differ. Our method also gives information about the magnetic ordering of the ground state as a function of the number of electrons. Our results are valid for arbitrary electron density values. The variational techniques we use are applied to the two-particle reduced-density matrix (2-RDM), because this quantity suffices to obtain expectation values of all one- and two-particle operators, and because it is less computationally intensive than obtaining the complete N -particle wave function. Because the 2-RDM must be positive semidefinite, we use an efficient programming algorithm known as the Semidefinite Programming Algorithm. A naive search through the space of two-particle matrices encounters a difficult problem, known as the " N -representability problem", which is to assure that the trial density matrices actually come from a wave function that has the required exchange symmetry. This problem is ameliorated by imposing constraints on the variational search that eliminate trial matrices with the incorrect exchange symmetry. These constraints express known physical properties of the ground state (for example its total spin and its translational symmetry) and the so-called D, Q , and G two-positivity conditions that express the required exchange symmetry. To calculate matrix elements, we use the two-particle eigenstates of the reduced Hubbard Hamiltonian; this choice of basis significantly decreases the number of matrix elements needed for energy calculations. After an optimal two-particle density matrix has been obtained, we calculate additional properties of the ground state, including spin correlation functions, chemical potentials, and the momentum distribution. We compare all of these properties to values obtained by exact numerical solution of the Schrödinger equation.
The many-boson system with repulsive interactions is treated by a variational method based on a variational trial state of an exponential pair-excitation type obtained by a generalization of that of Bogoliubov; the treatment is closely related to an intermediate-coupling approximation with respect to pairs. The nonlinear integral equation which determines the variational ground state is derived, and various properties of this ground state and the associated excited states are examined. The resultant low-lying spectrum lies below that of Bogoliubov by an amount proportional to the total number of particles. The variational principle is shown to produce rigorous energy eigenvalue differences for the pair part of the Hamiltonian. The variational states, however, still exhibit unphysical features characteristic of pair-excitation states: The pair correlation function does not go to zero at zero particle separation and the phonon spectrum exhibits a gap above the ground state. It is suggested that these features can be removed by using states which take into account excitation of momentum-conserving groups of more than two particles.
A detailed account is given of a recently developed pair theory of the Hubbard Hamiltonian for the case of a half-filled band. The resulting pair Hamiltonian is an exact reformulation of the original Hubbard Hamiltonian in terms of pair operators which satisfy elementary boson commutation relations. The Hamiltonian takes the familiar form of a sum of a quadratic part representing independent-pair energies and a quartic part representing pair interactions. The matrix elements entering the pair Hamiltonian are the pair energies and products of pair wave functions, which are formally given as the solutions of a certain eigenvalue equation. By the solving of the eigenvalue problem, explicit analytical expressions are obtained for both pair energies and wave functions for an arbitrary number of dimensions. For the case M=1, where M is the number of spin-down electrons, pair interactions do not enter since there is only one electron-hole pair, and the exact energies of the Hubbard Hamiltonian are obtained. This will be explicitly verified for the one-dimensional case; here we demonstrate that the pair energies agree with the exact energies found by other authors for this case.
Various methods and approximation schemes are used to study many-electron interacting systems. Two important many-particle models, the Anderson model and the Hubbard model, and their electromagnetic properties have been investigated in many parameter regimes, and applied to physical systems. An Anderson single-impurity model Hamiltonian based calculation of the magnetic susceptibility is performed for YbN in the presence of crystal fields using an alteration of the Non-Crossing Approximation proposed by Zwicknagl et.al., incorporating parameters obtained from ab initio band structure calculations. It yields good agreement with experimental data. For the Anderson lattice model, a variational scheme which uses specific many-electron wavefunctions as basis is applied to both one- and two-dimensional systems represented by symmetric Anderson lattice Hamiltonians. Without much computational effort, the ground state energy is well approximated, especially in strong-coupling limit. Some electronic properties are examined using the variational ground state wavefunction. The one-dimensional Hubbard model has been solved exactly for small-size clusters by diagonalizing the Hamiltonian in the basis of many-electron Bloch states. The results for the energy spectrum and eigenfunctions of the ground state and low-lying excited states are presented. Also, mean field calculations of the two-dimensional single-band Hubbard model and Cu-O lattice model (three-band Hubbard model) are carried out for various physical quantities including the energy, occupation probability, staggered magnetization, momentum distribution Fermi surface and density of states, by using a projection operator formalism. To develop a systematic approach to solving many-electron problems, the many-particle partition function for the free electron gas system is explored using a cumulant expansion scheme. Starting from the ground state, the partition function can be approximated to any order in terms of excitation energy. Its application to interacting systems such as the Anderson model and the Hubbard model is briefly discussed.
Abstract The study of ground-state properties of the Fermi-Hubbard model is a long-lasting task in the research of strongly correlated systems, and is crucial for the understanding of notable quantum phenomena including superconductivity and magnetism. 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.
The single-band Hubbard Hamiltonian study faces a serious limitation and difficulty as we move away from finite - size lattices to larger N - dimensional lattices. Thus there is the needto develop the means of overcoming the finite - size lattice defects as we pass on to a higher dimension.In this work, a quantitative approximation to the one-band Hubbard model is presented using a variational analytic approach. The goal of this work, therefore, is to explore quantitatively the lowest ground-state energy and the pairing correlations in 3D N x N x N lattices of the Hubbard model. We developed the unit step model as an approximate solution to the single-band Hubbard Hamiltonian to solve variationallythe correlation of two interacting elections on a three-dimensional cubic lattice. We also showed primarily how to derive possible electronic states available for several even and odd3D lattices, although, this work places more emphasis on a 3D 5 x 5 x 5 lattice. The results emerging from our present study compared favourablywith the results of Gutzwillervariational approach (GVA) and correlated variational approach (CVA), at thelarge limit of the Coulomb interaction strength (U/4t). It is revealed in this study, that the repulsive Coulomb interaction which in part leads to the strong electronic correlations, would indicate that the two electron system prefer not to condense into s-wave superconducting singlet state (s = 0), at high positive values of the interaction strength.