DMRG study of ferromagnetism in a one-dimensional Hubbard model

The one dimensional Hubbard model with nearest and (negative) next-nearest neighbour hopping has been studied with the density-matrix renormalization group (DMRG) method. A large region of ferromagnetism has been found for finite density and finite on-site interaction.     

Thermodynamic properties of the one-dimensional Hubbard model

We discuss the thermodynamic behaviour of the one-dimensional Hubbard model in the narrow-band regime, where the intra-atomic Coulomb-repulsion is large compared to the bandwidth. An approximation scheme on a perturbative basis is developed which applies for all temperatures. First order perturbation theory is performed for arbitrary electron densities; second order perturbation theory is discussed in the case of the half-filled band. Also the one-particle Green's function is calculated. Our approximation agrees excellently with numerical calculations. By comparison with exact results, which are available for some special limits, the range of validity of our approximation is estimated.     

Ladder operator for the one-dimensional Hubbard model

The one-dimensional Hubbard model is integrable in the sense that it has an infinite family of conserved currents. We explicitly construct a ladder operator which can be used to iteratively generate all of the conserved current operators. This construction is different from that used for Lorentz invariant systems such as the Heisenberg model. The Hubbard model is not Lorentz invariant, due to the separation of spin and charge excitations. The ladder operator is obtained by a very general formalism which is applicable to any model that can be derived from a solution of the Yang-Baxter equation.     

Entanglement scaling in the one-dimensional Hubbard model at criticality

We derive exact expressions for the local entanglement entropy epsilon in the ground state of the one-dimensional Hubbard model at a quantum phase transition driven by a change in magnetic field h or chemical potential mu. The leading divergences of delta epsilon/delta h and delta epsilon/delta mu are shown to be directly related to those of the zero-temperature spin and charge susceptibilities. Logarithmic corrections to scaling signal a change in the number of local states accessible to the system as it undergoes the transition.     

On the one-dimensional Hubbard model away from half-filling

We study the spectrum of the one-dimensional Hubbard model near the half-filled band, i.e. for electronic densityn=1–, 1, and also in the low-density limit,n1. Ground state and excited states are described, by integral equations which are derived from the Bethe Ansatz. These equations, which can be solved analytically only in the case of half-filling (=0), are investigated systematically in both limits 1,n1. Appropriate expansions of the momenta and energies of the ground state and the excited states as well as correlation exponents are determined.Work performed within the research program of the Sonderforschungsbereich 341, Köln-Aachen-JülichSimilar equations have been derived by Matveenko [3] using a different method     

Mott-insulator-superfluid-liquid transition in a one-dimensional bosonic Hubbard model: Quantum Monte Carlo method

The values of the insulator gap Δ in one-dimensional systems of interacting bosons described by the Hubbard Hamiltonian are calculated at low temperatures by the quantum world-line Monte Carlo algorithm. The dependence of Δ on the size of the system, the temperature, and the parameters of the model is investigated. It is shown that a chain with N _a=50 sites is already sufficient to estimate the thermodynamic value of the critical quantity (t/U)_c for which a transition from the insulator into the superfluid state occurs in a commensurate system. To within the computational error, this value, (t/U)_c=0.300±0.005, agrees with the value (t/U)_c=0.304±0.002 obtained previously by the combined “exact diagonalization + renormalization-group analysis” method. The characteristic Kosterlitz-Thouless behavior of the insulator gap is demonstrated near the critical region: Δ∼exp[−b(1−t/t _c)^−1/2]. Pis’ma Zh. éksp. Teor. Fiz. 64, No. 2, 92–96 (25 July 1996)     

Correlations in the ground state of the one-dimensional Hubbard model

With eigenfunctional theory and a rigorous expression of exchange-correlation energy of a general interacting electron system, we study the ground state properties of the one-dimensional Hubbard model, and calculate the ground-state energy as well as the charge gap at half-filling for arbitrary coupling strength u=U/(4t) and electron density nc. We find that the simple linear approximation of the phase field works well in weak coupling case, but it becomes inappropriate as the on-site Coulomb interaction becomes strong where the fluctuations of the bosonic auxiliary field are strong. Then we propose a new scheme by adding Gutzwiller projection which suppresses the density fluctuations and the new results are quite close to the exact ones up to considerably strong coupling strength u=3.0 and for arbitrary electron density nc. Our calculation scheme is proved to be effective for strongly correlated electron systems in one dimension, and its extension to higher dimensions is straightforward.     

一维扩展离子Hubbard模型的相图研究

应用密度矩阵重整化群方法, 研究了存在交错离子势Δ时一维半满扩展Hubbard模型的相图. 通过计算关联函数、结构因子、位置算符等方法, 描绘了从Mott绝缘体-键有序绝缘体-Band 绝缘体的特性并给出了精确的相边界. 研究发现: 中间的键有序绝缘体相在相图中占据了很小的一部分区域, 当存在离子势Δ的情况下, 这个区域将会有所增大; 而当相互作用足够强时, 这个中间相消失. 给出了离子Hubbard模型(最近邻电子-电子相互作用V=0)的相图.     

Phase diagram of the one-dimensional half-filled extended Hubbard model

We determine the ground-state phase diagram of the one-dimensional half-filled Hubbard model with on-site (nearest-neighbor) repulsive interaction U (V) and nearest-neighbor hopping t using the density-matrix renormalization group technique. Based on the results of the excitation gaps, Luttinger-liquid exponents, and bond-order-wave (BOW) order parameter, we confirm that the BOW phase appears in a substantial region between the charge-density-wave (CDW) and spin-density-wave phases. Each phase boundary is determined by multiple means and it allows us to make a cross-check on the validity of our estimations. We also find that the BOW-CDW transition changes from continuous to first order at the tricritical point (U(t),V(t)) approximately (5.89 t,3.10 t) and the BOW phase shrinks to zero at the critical end point (U(c),V(c)) approximately (9.25 t,4.76 t).