Abstract This paper describes a group theoretical classification of the unrestricted Hartree–Fock (UHF) solution of the extended Hubbard model including nearest‐neighbor repulsion and exchange interactions. We derived systematically all types of the broken symmetry solution with a single‐order parameter for the ordering vector Q = (π,π). It is shown that there are four nonmagnetic states: a charge density wave (CDW), a charge current wave (CCW), and two bond‐order waves (BOW), and five magnetic states: an ordinary antiferromagnetic state (AF), a spin current wave (SCW), two axial spin bond‐order waves (ASBOW), and an helical spin bond‐order wave (HSBOW). We showed that these phases can be derived naturally from the irreducible decomposition of the Fourier transformed interaction and the self‐consistent field (SCF) condition. Some numerical results are given for the half‐filled case.
The ground states of the two-dimensional repulsive Hubbard model are studied within the unrestricted Hartree-Fock (UHF) theory. Magnetic and charge properties are determined by systematic, large-scale, exact numerical calculations, and quantified as a function of electron doping, h. In the solution of the self-consistent UHF equations, multiple initial configurations and simulated annealing are used to facilitate convergence to the global minimum. New approaches are employed to minimize finite-size effects in order to reach the thermodynamic limit. At low to moderate interacting strengths and low doping, the UHF ground state is a linear spin-density wave (l-SDW), with antiferromagnetic order and a modulating wave. The wavelength of the modulating wave is 2/h. Corresponding charge order exists but is substantially weaker than the spin order, hence holes are mobile. As the interaction is increased, the l-SDW states evolve into several different phases, with the holes eventually becoming localized. A simple pairing model is presented with analytic calculations for low interaction strength and small doping, to help understand the numerical results and provide a physical picture for the properties of the SDW ground state. By comparison with recent many-body calculations, it is shown that, for intermediate interactions, the UHF solution provides a good description of the magnetic correlations in the true ground state of the Hubbard model.
Nonlinear doping states in Peierls–Hubbard models are investigated through response functions in magnetic, optical and vibronic channels. A random phase approximation is used for their analysis on the basis of totally unrestricted Hartree–Fock solutions, which in general have broken symmetries. The validity of the present method is first discussed through a comparison with known exact results for a two-site two-electron system. Then this method is applied to a two-dimensional three-band system with one additional hole relative to stoichiometry, for modeling the lightly doped copper oxides. New excitation modes are induced by doping and observed in the different channels of response functions. These doping-induced modes distinguish between various doping states and helpful for the identification of doping states in specific materials.
We study the two-dimensional one-band Hubbard model on the square lattice through the spin rotation invariant Kotliar-Ruckenstein slave-boson method, with an emphasis on the electron-doped regime. The phase diagram determined within mean field approximation reveals three different magnetically ordered phases separated by first-order transitions, and a wide instability region as indicated by a negative compressibility. We show that this conundrum of the unstable mean field domain is resolved by spatially mixed phases of neighboring stable states. We also find charge density wave type instabilities in the antiferromagnetically ordered (AFM) state at wave vectors commensurate with the doping. The corresponding charge and spin density response functions are calculated from a Gaussian fluctuations formalism about a magnetic saddle point, for which we provide a detailed derivation. Further, we present density and spin excitation spectra in the AFM phase, allowing for a quantitative calculation of the spin wave parameters, the gap in the longitudinal spin excitation spectrum, and the density excitations in the upper Hubbard band.
We present an unbiased mean-field analysis of magnetic and charge orders in the two-dimensional Hubbard model on a square lattice, at both zero and finite temperatures. Unrestricted Hartree-Fock calculations on large finite lattices are complemented by solutions restricted to N\'eel and circular spiral order in the thermodynamic limit. The magnetic states are classified by a systematic scheme based on the dominant Fourier components of the spin texture. On finite lattices a whole zoo of ordering patterns appears. We show that many of these states are finite-size artifacts related to the limited choice of ordering wave vectors on a finite lattice. In the thermodynamic limit only three classes of states with a relatively simple structure survive: N\'eel, circular spiral, and stripe states. Stripes involve also charge order and can be unidirectional or bidirectional, with horizontal and/or vertical orientation. We present complete phase diagrams in the plane spanned by electron density and temperature, for a moderate Hubbard interaction and various choices of the next-nearest-neighbor hopping amplitude.
The two-dimensional Hubbard model is one of the most fundamental models in condensed matter physics. It captures many properties of cuprate high-temperature superconductors. Despite its apparent simplicity, solving the Hubbard model is a challenging task. In this work, we employ the Hartree-Fock approximation, as well as functional renormalization group methods, to construct phase diagrams of the Hubbard model and to shed light on the zoo of magnetic, charge ordered, and superconducting states. The Hartree-Fock approximation can self-consistently determine the expectation values of local charge and magnetic order parameters. Mean-field theories like this can lead to a deeper understanding of the competition, coexistence and transition of different orders. We start by employing the Hartree-Fock approximation in a completely unrestricted form on finite-size lattices and complement the results with calculations in the thermodynamic limit. We construct mean-field phase diagrams of the two-dimensional Hubbard model in the thermodynamic limit for different values of the next-nearest neighbor hopping as a function of filling and temperature. We find that the dominant magnetic orders at moderate coupling strength are Néel order around half-filling, followed by spiral order, and stripe order as the hole doping is increased. We also clarify the nature of the transition from magnetic spiral to stripe order. We then extend the Hartree-Fock scheme by a functional renormalization group flow, which generates new effective interactions and allows us to also access the superconducting phase of the Hubbard model. We present an improved phase diagram, featuring Néel, spiral, and stripe order coexisting with d-wave superconductivity at low temperatures. We analyze the coexistence and competition between these orders. Finally, we perform unrestricted Hartree-Fock calculations of a Hubbard model on the graphene lattice and the Bernal bilayer graphene lattice. We also construct magnetic phase diagrams for these geometries, which feature Néel order around half-filling and various stripe orders in a secondary magnetic dome.