Monte Carlo simulations are used to study the charge-density-wave to spin-density-wave transition in the one-dimensional Hubbard model with on-site (
U) and nearest-neighbor (
V) repulsion. The transition is found to be continuous for small
U and first order for large
U, with the crossover occurring around $U\ensuremath{\sim}3$. The transition line deviates somewhat from the theoretical prediction
U=2V towards the charge-density-wave phase for intermediate coupling. A qualitative explanation for the results is given by use of strong-coupling arguments.
We present a generalization of the recently proposed variational cluster perturbation theory to extended Hubbard models at half-filling with repulsive nearest neighbor interaction. The method takes into account short-range correlations correctly by the exact diagonalization of clusters of finite size, whereas long-range order beyond the size of the clusters is treated on a mean-field level. For one dimension, we show that quantum Monte Carlo and density-matrix renormalization-group results can be reproduced with very good accuracy. Moreover we apply the method to the two-dimensional extended Hubbard model on a square lattice. In contrast to the one-dimensional case, a first order phase transition between spin density wave phase and charge density wave phase is found as function of the nearest-neighbor interaction at onsite interactions $U\ensuremath{\geqslant}3t$. The single-particle spectral function is calculated for both the one-dimensional and the two-dimensional system.
This thesis addresses the question of a strongly interacting many-body problem: the Hubbard model on one- and two-dimensional lattices. The technique employed is the small-cluster method wherein the full many-body Hamiltonian for finite clusters with periodic boundary conditions (PBC's) is solved exactly. We first apply this technique to an extended one-dimensional chain and a two-dimensional square-lattice Hubbard model. We then construct a two-dimensional overlayer/substrate Hubbard model with first- and second-nearest-neighbor hopping, on-site and nearest-neighbor interactions, and overlayer-site orbital energy. The ground-state properties are investigated and discussed in a many-body picture.
The recently introduced concept of exact mean-field theories is applied to the extended simplified Hubbard model, which is a one-band extended Hubbard model (including nearest-neighbor interaction) where only the \ensuremath{\downarrow} electrons can hop. An exact mean-field Hamiltonian is derived for the model on a Bethe lattice in the limit of a large coordination number (Z\ensuremath{\rightarrow}\ensuremath{\infty}). The mean-field Hamiltonian is used to calculate the Green functions and the free energy at half-filling. This exact solution is analyzed analytically at small and large values of the on-site (U) and nearest-neighbor (V) interaction strengths. One finds that a phase transition occurs at sufficiently low temperatures for all U,V>0. The low-temperature phase is a charge-density wave for V>1/2U and a spin-density wave for V1/2U. Special attention is paid to (i) the critical temperature as a function of U and V, (ii) the order parameter as a function of temperature, and (iii) the density of states. The density of states shows that the low-temperature phase is an insulator for all U,V provided that V/U\ensuremath{\ne}1/2, while the high-temperature phase is a paramagnetic metal for all V if U2, and a paramagnetic insulator if U>2. Many detailed results are given.
We consider a generalized Hubbard model with on-site interaction U, nearest-neighbor repulsion V, and nearest-neighbor hopping for spin \ensuremath{\sigma}, which depends on the sum of particles ${\mathrm{m}}_{\mathrm{\ensuremath{\sigma}}\ensuremath{-}|\mathrm{A}\mathrm{m}}$ with opposite spin in the two sites involved. The hopping matrix elements are denoted by
tAA,
tAB,
tBB for ${\mathrm{m}}_{\mathrm{\ensuremath{\sigma}}\ensuremath{-}|\mathrm{A}\mathrm{m}}$=0,1,2, respectively. For 0
tABtAA=
tBB, we have determined the regions of parameters for which the ground-state (GS) of the one-dimensional system is a charge-density wave (CDW), a spin-density wave (SDW), or a metal (M), using Hartree-Fock, exact diagonalization of finite chains and quantum Monte Carlo. The results agree qualitatively with the exactly solvable limit of
tAB=0. For 0
tABtAA=
tBB, the GS is a M for sufficiently low values of U and V. In contrast, when
tAA+
tBB-
2tAB=0, our results suggest that the GS is either a CDW or a SDW, with the boundary between them lying near the line U=2V.
We use the density-matrix renormalization group to calculate the phase diagram of the one-dimensional extended Hubbard chain at half filling. We find that for one specific
U, the critical value of
Vc for the charge-density-wave versus spin-density-wave phase transition is smaller than that predicted from the Monte Carlo simulation. Our results are consistent with the strong- and weak-coupling limit results, whereas the previous results from the Monte Carlo simulations deviate from the weak-coupling limit result. As
U increases, we find that the first-order phase transition develops gradually. The tricritical point is estimated to occur at
Uc=4−−6, not
Uc=3 as previously predicted from the Monte Carlo simulation. Due to the topological difference, we observe that with the periodic boundary condition, the phase transition is more abrupt for finite chains of length
N=4M than for chains of
N=4M+2, where
M is an integer, for the same
U and
V. The difference will diminish in the thermodynamic limit.