The one-dimensional 1/4-filled Hubbard model with alternating interactions

The one-dimensional Hubbard model with different on-site interaction on the even (U_a) and odd (U_b) sites is considered within the framework of the weak coupling approach. In the case of a 1/4-filled band the dynamical nonequivalence of sites leads to the appearance of Umklapp processes in the system and to the dynamical generation of a commensurability gap in the charge excitation spectrum for U_a ≠ U_b and U_a>0 or U_b>0. Depending on the relation between the bare coupling constants the system shows four different regimes of behaviour in the infrared limit corresponding to normal metal, nonmagnetic insulator, antiferromagnetic insulator and superconducting states. The extended model including interaction between particles on nearest and next-nearest neighbour sites is also considered.     

Resonance state in the 2D attractive Hubbard model

The systematic change of a resonance state with high momenta is studied with increasing particle density in the 2D attractive Hubbard model. Within the conserving self-consistent T-matrix approximation, we present the spectral functions for the one and two particle Green's functions as well as the self-energy. In the small density limit, the resonant state becomes stable and the result from the self-consistent calculations shows a good agreement with that from a simple analytical calculation. As particle density is increased, the resonance state acquires a short lifetime due to the increasing decay into two free particles.     

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

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

Decorated star-triangle relations and exact integrability of the one-dimensional Hubbard model

The exact integrability of the one-dimensional Hubbard model is demonstrated with the help of a novel set of triangle relations, the decorated star-triangle relations. The covering two-dimensional statistical mechanical model obeys the star-triangle or Yang-Baxter relation. A conjecture is presented for the eigenvalues of the transfer matrix.     

有限温度下一维Hubbard模型的化学势泛函理论研究

通过数值方法求解了有限温度下一维均匀Hubbard模型的热力学Bethe-ansatz方程组,得到了在给定温度和相互作用强度情况下,比热c、磁化率χ和压缩比κ随化学势μ的变化图像.基于有限温度下一维均匀Hubbard模型的精确解,利用化学势(μ)-泛函理论研究了一维谐振势下的非均匀Hubbard模型,给出了金属态和Mott绝缘态下不同温度情况时局域粒子密度n_i和局域压缩比_κi随格点的变化情况.     

Quasiparticle scattering interference in the renormalized Hubbard model

In this paper, we study the quasiparticle scattering interference phenomenon in the presence of a single impurity within the renormalized Hubbard model. By calculating the energy and momentum dependence of the Fourier-transformed local density of states in the full Brillouin zone, we can qualitatively describe the main features of the quasiparticle scattering interference phenomenon in cuprate superconductors using a single point-like impurity. In particular, we show that with increasing energy, the position of the peak along the nodal ([0, 0] → [π, π]) direction moves steadily to a large momentum region, while the position of the peak along the antinodal ([0, 0] → [π, 0]) direction moves toward the center of the Brillouin zone.     

On ground states of the one-dimensional ferromagnetic XXZ model

We obtain the complete list of pure infinite volume ground states for the one-dimensional ferromagnetic XXZ model.     

Thermoelectric power of high-T c superconductors using extended Hubbard model

The thermoelectric power (TEP) for a one dimensional lattice has been studied by using the extended Hubbard model in the limitU≠∞, whereU is the on-site Coulomb interaction. A new expression for TEP, derived in this study, has been found to successfully reproduce the experimental results of the insulating samples of different hole- and electron-doped high temperature superconducting systems.     

Flow equations and the strong-coupling expansion for the Hubbard model

Applying the method of continuous unitary transformations to a class of Hubbard models, we reexamine the derivation of thet/U expansion for the strong-coupling case. The flow equations for the coupling parameters of the higher order effective interactions can be solved exactly, resulting in a systematic expansion of the Hamiltonian in powers oft/U, valid for any lattice in arbitrary dimension and for general band filling. The expansion ensures a correct treatment of the operator products generated by the transformation, and only involves the explicit recursive calculation of numerical coefficients. This scheme provides a unifying framework to study the strong-coupling expansion for the Hubbard model, which clarifies and circumvents several difficulties inherent to earlier approaches. Our results are compared with those of other methods, and it is shown that the freedom in the choice of the unitary transformation that eliminates interactions between different Hubbard bands can affect the effective Hamiltonian only at ordert ³/U² or higher.     

Simulations of the Hubbard model

We discuss results of simulations of the Hubbard model of interacting electrons on a lattice. We start with a brief discussion of methodology and point out some of the outstanding problems. We then discuss results of simulations of the model in three, two, and one dimension, particularly in connection with its magnetic and superconducting properties. We conclude with a brief discussion of future directions.     

Analyticity in Hubbard Models

The Hubbard model describes a lattice system of quantum particles with local (on-site) interactions. Its free energy is analytic when t is small, or t ²/U is small; here, is the inverse temperature, U the on-site repulsion, and t the hopping coefficient. For more general models with Hamiltonian H=V+T where V involves local terms only, the free energy is analytic when T is small, irrespective of V. There exists a unique Gibbs state showing exponential decay of spatial correlations. These properties are rigorously established in this paper.     

Strong-coupling perturbation theory of the Hubbard model

The strong-coupling perturbation theory of the Hubbard model is presented and carried out to order (t/U)⁵ for the one-particle Green function in arbitrary dimension. The spectral weight is expressed as a Jacobi continued fraction and compared with new Monte-Carlo data of the one-dimensional, half-filled Hubbard model. Different regimes (insulator, conductor and short-range antiferromagnet) are identified in the temperature-hopping integral (T,t) plane. This work completes a first paper on the subject (Phys. Rev. Lett. 80, 5389 (1998)) by providing details on diagrammatic rules and higher-order results. In addition, the non half-filled case, infinite resummations of diagrams and the double occupancy are discussed. Various tests of the method are also presented. Received 25 October 1999     

半满Hubbard模型的新处理方法

从标准Hubbard模型出发,对半满狭带强关联系统,在讨论自旋为δ的电子运动时,忽略自旋为-δ的电子在格点间的跳跃,得到不对称哈密顿量,并运用平均场近似求得相应准粒子谱．在绝对零度时,与Green函数近似解作了比较;在有限温度时,讨论了从绝缘体到金属相变的可能．     

Zero finite-temperature charge stiffness within the half-filled 1D Hubbard model

Even though the one-dimensional (1D) Hubbard model is solvable by the Bethe ansatz, at half-filling its finite-temperature T>0$T>0$ transport properties remain poorly understood. In this paper we combine that solution with symmetry to show that within that prominent T=0$T=0$ 1D insulator the charge stiffness D(T)$D\left(T\right)$ vanishes for T>0$T>0$ and finite values of the on-site repulsion U$U$ in the thermodynamic limit. This result is exact and clarifies a long-standing open problem. It rules out that at half-filling the model is an ideal conductor in the thermodynamic limit. Whether at finite T$T$ and U>0$U>0$ it is an ideal insulator or a normal resistor remains an open question. That at half-filling the charge stiffness is finite at U=0$U=0$ and vanishes for U>0$U>0$ is found to result from a general transition from a conductor to an insulator or resistor occurring at U=U_c=0$U=U_c=0$ for all finite temperatures T>0$T>0$. (At T=0$T=0$ such a transition is the quantum metal to Mott-Hubbard-insulator transition.) The interplay of the η$\eta$-spin SU(2)$SU\left(2\right)$ symmetry with the hidden U(1)$U\left(1\right)$ symmetry beyond SO(4)$SO\left(4\right)$ is found to play a central role in the unusual finite-temperature charge transport properties of the 1D half-filled Hubbard model.     

Generalized Hartree-Fock theory and the Hubbard model

The familiar unrestricted Hartree-Fock variational principles is generalized to include quasi-free states. As we show, these are in one-to-one correspondence with the one-particle density matrices and these, in turn, provide a convenient formulation of a generalized Hartree-Fock variational principle, which includes the BCS theory as a special case. While this generalization is not new, it is not well known and we begin by elucidating it. The Hubbard model, with its particle-hole symmetry, is well suited to exploring this theory because BCS states for the attractive model turn into usual HF states for the repulsive model. We rigorously determine the true, unrestricted minimizers for zero and for nonzero temperature in several cases, notably the half-filled band. For the cases treated here, we can exactly determine all broken and unbroken spatial and gauge symmetries of the Hamiltonian.Dedicated to Philippe Choquard on his 65th birthday.     

Quantum phase diagram of the half filled Hubbard model with bond-charge interaction

Using quantum field theory and bosonization, we determine the quantum phase diagram of the one-dimensional Hubbard model with bond-charge interaction X in addition to the usual Coulomb repulsion U at half-filling, for small values of the interactions. We show that it is essential to take into account formally irrelevant terms of order X  . They generate relevant terms proportional to X²$X^2$ in the flow of the renormalization group (RG). These terms are calculated using operator product expansions. The model shows three phases separated by a charge transition at U=Uc$U=U_c$ and a spin transition at U=Us>Uc$U=U_s>U_c$. For U<Uc$U<U_c$ singlet superconducting correlations dominate, while for U>Us$U>U_s$, the system is in the spin-density wave phase as in the usual Hubbard model. For intermediate values Uc<U<Us$U_c<U<U_s$, the system is in a spontaneously dimerized bond-ordered wave phase, which is absent in the ordinary Hubbard model with X=0$X=0$. We obtain that the charge transition remains at Uc=0$U_c=0$ for X≠0$X\ne 0$. Solving the RG equations for the spin sector, we provide an analytical expression for Us(X)$U_s(X)$. The results, with only one adjustable parameter, are in excellent agreement with numerical ones for X<t/2$X<t/2$ where t is the hopping.     

Magnetism in the single-band Hubbard model

A self-consistent spectral density approach (SDA) is applied to the Hubbard model to investigate the possibility of spontaneous ferro- and antiferromagnetism. The starting point is a two-pole ansatz for the single-electron spectral density, the free parameter of which can be interpreted as energies and spectral weights of respective quasiparticle excitations. They are determined by fitting exactly calculated spectral moments. The resulting self-energy consists of a local and a non-local part. The higher correlation functions entering the spin-dependent local part can be expressed as functionals of the single-electron spectral density. Under certain conditions for the decisive model parameters (Coulomb interaction U, Bloch bandwidth W, band occupation n, temperature T) the local part of the self-energy gives rise to a spin-dependent band shift, thus allowing for spontaneous band magnetism. As a function of temperature, second-order phase transitions are found away from half-filling, but close to half-filling, the system exhibits a tendency towards first-order transitions. The non-local self-energy part is determined by use of proper two-particle spectral densities. Its main influence concerns a (possibly spin-dependent) narrowing of the quasiparticle bands with the tendency to stabilize magnetic solutions. The non-local self-energy part disappears in the limit of infinite dimensions. We present a full evaluation of the Hubbard model in terms of quasiparticle densities of states, quasiparticle dispersions, magnetic phase diagram, critical temperatures (T_c, T_N) as well as spin and particle correlation functions. Special attention is focused on the non-locality of the electronic self-energy, for which some rigorous limiting cases are worked out.     

Rigorous solution of a Hubbard model extended by nearest‐neighbour Coulomb and exchange interaction on a triangle and tetrahedron

The Hubbard model extended by both nearest‐neighbour (nn) Coulomb correlation and nearest‐neighbour Heisenberg exchange is solved rigorously for a triangle and tetrahedron. All eigenvalues and eigenvectors are given as functions of the model parameters in a closed analytical form. For fixed electron numbers we found a multitude of level crossings, both in the ground state and in the excited states in dependence on the various model parameters. By coupling an ensemble of clusters to an electron bath we get the cluster gas model or the cluster gas approximation, if an extended array of weak‐interacting clusters is considered. The grand‐canonical potential Ω (μ, T, h) and the electron occupation N (μ, T, h) of the related cluster gases were calculated for arbitrary values (attractive and repulsive) of the three interaction constants. For the cluster gases without the additional interactions we found various steps in N (μ, T = 0, h = 0) higher than one. The reason is the degeneration of ground states differing in their electron occupation by more than one electron. For the triangular cluster gas we have one such degeneration point. For the tetrahedral cluster gas two. As a consequence, we do not find areas with one electron in the μ‐U ground‐state phase diagram of the triangular cluster gas or with one, two and five electrons in the case of the tetrahedral cluster gas. The degeneration point of the triangular cluster gas can not be destroyed by an applied magnetic field. This holds also for the lower degeneration point of the tetrahedral cluster gas. Otherwise, the upper degeneration point breaks down at a critical magnetic field h_c. The dependence of h_c on U shows a maximum for strong on‐site correlation. The influence of nn‐exchange and nn‐Coulomb correlation on the ground‐state phase diagrams is calculated. Whereas antiferromagnetic nn‐exchange breaks the degeneration points of the tetrahedral cluster gas partially only, a repulsive nn‐Coulomb correlation lifts the underlying degeneracies completely. Otherwise both ferromagnetic nn‐exchange and attractive nn‐Coulomb interaction stabilise the degeneration points. The consequences of the cluster gas results for extended cluster arrays are discussed.