Ferromagnetism in the Hubbard model: instability of the Nagaoka state on the triangular,honeycomb and kagome lattices

In order to analyse the lattice dependence of ferromagnetism in the two-dimensional Hubbard model we investigate the instability of the fully polarised ferromagnetic ground state (Nagaoka state) on the triangular, honeycomb and kagome lattices. We mainly focus on the local instability, applying single spin flip variational wave functions which include majority spin correlation effects. The question of global instability and phase separation is addressed in the framework of Hartree-Fock theory. We find a strong tendency towards Nagaoka ferromagnetism on the non-bipartite lattices (triangular, kagome) for more than half filling. For the triangular lattice we find the Nagaoka state to be unstable above a critical density of n = 1.887 at U = ∞, thereby significantly improving former variational results. For the kagome lattice the region where ferromagnetism prevails in the phase diagram widely exceeds the flat band regime. Our results even allow the stability of the Nagaoka state in a small region below half filling. In the case of the bipartite honeycomb lattice several disconnected regions are left for a possible Nagaoka ground state.     

Ferromagnetism in the Hubbard model: Spin waves and instability of the Nagaoka state

We discuss two single spin flip variational wave functions describing spin wave excitations which were proposed earlier by Shastry, Krishnamurthy and Anderson (SKA) and by Basile and Elser (BE), respectively, in order to investigate the instability of the fully polarized ferromagnetic state (Nagaoka state) in the infinite U Hubbard model. We calculate the energy of these variational states for the square lattice and for multiple chains. At the zone boundary in the vicinity of the point (0, π) the spin wave energy is reduced substantially by the binding of the spin up hole to the flipped down spin. For the square lattice this leads to a critical hole density of δ_cr = 0.407 for the SKA spin wave and of δ_cr = 0.322 for the BE spin wave which implies remarkable improvements in comparison to the corresponding scattering states investigated previously.     

Ferromagnetism in the Hubbard model on the square lattice: Improved instability criterion for the Nagaoka state

The instability of the fully polarized ferromagnetic state (Nagaoka state) with respect to single spin flips is re-examined for the Hubbard model on the square lattice with a large family of variational wave functions which include correlation effects of the majority spins in the vicinity of the flipped spin. We find a critical hole density of δ_cr = 0.251 for U = ∞ and a critical coupling of U_cr = 77.7t. Both values improve previous variational results considerably.     

Ferromagnetism in the Hubbard model: instability of the Nagaoka state on the square lattice

The instability of the fully polarized ferromagnetic ground state (Nagaoka state) of the Hubbard model on the square lattice is investigated. We use single spin flip variational wave functions including majority spin correlation effects and calculate spin flip energies in the thermodynamic limit. With very local wave functions and with moderate numbers of variational parameters we reproduce the best known estimate for the critical hole density δ_cr = 0.29 and we obtain an estimate of U_cr = 63 t for the critical coupling which is considerably better than the best estimate of U_cr = 42 t previously known. The simplicity of our wave functions makes the physical origin of the various aspects of the instability particularly transparent.     

Absence of saturated ferromagnetism in the two-dimensional Hubbard model with two holes for U=∞

For a square Hubbard lattice with infinite repulsion energy U the following exact result has been obtained: the ferromagnetic state with maximum spin is not the ground state of the system if the number of holes is equal to two. Zh. éksp. Teor. Fiz. 113, 1000–1008 (March 1998)     

Spin-wave spectrum of a two-dimensional itinerant electron system: Analytic results for the incommensurate spiral phase in the strong-coupling limit

We study the zero-temperature spin fluctuations of a two-dimensional itinerant-electron system with an incommensurate magnetic ground state described by a single-band Hubbard Hamiltonian. We introduce the (broken-symmetry) magnetic phase at the mean-field (Hartree-Fock) level through a spiral spin configuration with characteristic wave vector Q different in general from the antiferromagnetic wave vector Q _AF, and consider spin fluctuations over and above it within the electronic random-phase (RPA) approximation. We obtain a closed system of equations for the generalized wave vector and frequency dependent susceptibilities, which are equivalent to the ones reported recently by Brenig. We obtain, in addition, analytic results for the spin-wave dispersion relation in the strong-coupling limit of the Hubbard Hamiltonian and find that at finite doping the spin-wave dispersion relation has a hybrid form between that associated with the (localized) Heisenberg model and that associated with the (long-range) RKKY exchange interaction. We also find an instability of the spin-wave spectrum in a finite region about the center of the Brillouin zone, which signals a physical instability toward a different spin- or, possibly, charge-ordered phase, as, for example, the stripe structures observed in the high-T _c materials. We expect, however, on physical grounds that for wave vectors external to this region the spin-wave spectrum that we have determined should survive consideration of more sophisticated mean-field solutions. Received 15 September 2000     

Coexistence of antiferromagnetic and ferromagnetic phase for ferromagnetic Kondo lattice model

To study the proposed phase separations in doped manganites, we performed Monte-Carlo calculations for the ferromagnetic Kondo lattice model with strong Hund's coupling between conduction electrons and localized spins. For the practical calculations, we adopted a one dimensional lattice and treated the spins of the localized t_2g electrons semi-classically. A direct evidence of the phase separation is observed from a snapshot of the spatial dependence of localized spins. No indication of the canted or spiral phases is found in the results of simulations. Further, the calculated results of the spin structure factor in the phase separation region are well compared with recent experiments. Received: 1st September 1998 / Revised: 30 October 1998 / Accepted: 27 November 1998     

Effective Hubbard model for Zn-doped CuO2 plane

In the framework of the cell-perturbation method for the original p-d model an effective two-band Hubbard model for the CuO₂ plane with Zn impurities is derived. Zn impurities are modelled by Wannir oxygen one-hole states at vacant Cu sites. The model is based on the results of band structure calculations carried out within the local-density approximation. Further reduction to an extended t-J model shows a large ferromagnetic superexchange interaction between the Cu spin with the nearest virtual oxygen spin in the Zn cell. Received 17 November 1998     

Vortices in quantum spin systems

We examine spin vortices in ferromagnetic quantum Heisenberg models with planar anisotropy on two-dimensional lattices. The symmetry properties and the time evolution of vortices built up from spin-coherent states are studied in detail. Although these states show a dispersion typical for wave packets, important features of classical vortices are conserved. Moreover, the results on symmetry properties provide a construction scheme for vortex-like excitations from exact eigenstates, which have a well-controlled time evolution. Our approach works for arbitrary spin length both on triangular and square lattices. Received 2 October 1998     

Variational approach to spin fluctuations in the Hubbard model

We study a one parameter variational wave function to improve the spin density wave ground state of the Hubbard model by inclusion of quantum spin fluctuations. Using a perturbative approach and novel lattice summation techniques we present analytical as well as numerical results for the correlation energies and the staggered magnetizations in one and two dimensions. We find ground state energies which are satisfyingly close to known exact results and are significantly lower than those of existing Gutzwiller and numerical treatments.     

Density-matrix renormalization study of the Hubbard model [4]on a Bethe lattice

The half-filled Hubbard model on the Bethe lattice with coordination number z=3 is studied using the density-matrix renormalization group (DMRG) method. Ground-state properties such as the energy per site E, average local magnetization , its fluctuations and various spin correlation functions are determined as a function of the Coulomb interaction strength U/t. The local magnetic moments increase monotonically with increasing Coulomb repulsion U/t showing antiferromagnetic order between nearest neighbors []. At large U/t, is strongly reduced with respect to the saturation value 1/2 due to exchange fluctuations between nearest neighbors (NN) spins [ for ]. shows a maximum for U/t=2.4-2.9 that results from the interplay between the usual increase of with increasing U/t and the formation of important permanent moments at large U/t. While NN sites show antiferromagnetic spin correlations that increase with increasing Coulomb repulsion, the next NN sites are very weakly correlated over the whole range of U/t. The DMRG results are discussed and compared with tight-binding calculations for U=0, independent DMRG studies for the Heisenberg model and simple first-order perturbation estimates. Received 8 February 1999 and Received in final form 14 June 1999     

On evaluation of transverse spin fluctuations in quantum magnets

A numerical method is described for evaluating transverse spin correlations in the random phase approximation. Quantum spin-fluctuation corrections to sublattice magnetization are evaluated for the antiferromagnetic ground state of the half-filled Hubbard model in two and three dimensions in the whole U/t range. Extension to the case of defects in the AF is also discussed for spin vacancies and low-U impurities. In the limit, the vacancy-induced enhancement in the spin fluctuation correction is obtained for the spin-vacancy problem in two dimensions, for vacancy concentration up to the percolation threshold. For low-U impurities, the overall spin fluctuation correction is found to be strongly suppressed, although surprisingly spin fluctuations are locally enhanced at the low-U sites. Received 27 April 1998 and Received in final form 13 August 1998     

Ferromagnetic ground state in the Kondo lattice model

We present a series of rigorous examples of the Kondo lattice model that exhibit full ferromagnetism in the ground state. The models are defined in one-, two- and three-dimensional lattices, and are characterized by a range of hopping terms, specific electron filling, and large ferromagnetic coupling. Our examples show that a sufficient strong but finite exchange coupling between conduction electrons and localized spins could overcome the competition from mobility of a finite density of electrons and drive the system from a paramagnetic phase to a ferromagnetic phase. We also establish a relation of ferromagnetism between the Hubbard model and Kondo lattice model. Meanwhile some rigorous results on ferromagnetism in the corresponding Hubbard model are presented. Received: 10 September 1997 / Revised: 15 October 1997 / Accepted: 17 October 1997     

Effective electron Hamiltonian for a quasi one-dimensional system. The (TMTSF)2X and the (TMTTF)2X systems

In this paper we have introduced a variational approach to investigate the ground state of a model which includes both the Holstein electron-phonon interaction and the extended Hubbard electron-electron interaction. We have considered a variational state for the phonon subsystem which generalizes the previous used forms. This state allows to take into account the possibility of extended phonon mediated correlations. The effective electron Hamiltonian, which we have obtained, includes first and second neighbor electron-electron interaction terms. We have treated exactly, through a Lanczos method, this effective model in the one-dimensional case. We have applied our method to two Bechgaard salts and in these cases we have estimated the correlation parameters. We have shown that the introduction of electron-phonon interaction allows an estimate of the on site U and nearest-neighbor V Coulomb repulsion, which are in agreement with the experimental optical spectra of the above mentioned two compounds. Received: 30 October 1997 / Revised: 28 January 1998 / Accepted: 10 April 1998     

Spin-triplet superconductivity induced by off-site Coulomb interaction in a triangular lattice

In order to study how superconductivity emerges when the charge fluctuation coexists with the spin fluctuation in a triangular lattice, we obtain a phase diagram against the nearest-neighbor Coulomb repulsion V and band filling n on the extended Hubbard model using the fluctuation exchange (FLEX) approximation. We find that a charge density wave (CDW) phase exists in the region where the value of V is large, and the f-wave spin-triplet pairing mediated by a strong charge fluctuation is dominant near the CDW phase.     

Heavy electron YBNi2B2C and giant exchange YbNiBC: 170Yb Mössbauer spectroscopy and magnetization studies

We present magnetic properties of the three-band Hubbard model in the para- and antiferromagnetic phase on a hypercubic lattice calculated with the Dynamical Mean-Field Theory (DMFT). To allow for solutions with broken spin-symmetry we extended the approach to lattices with AB-like structure. Above a critical sublattice magnetization one can observe rich structures in the spectral-functions similar to the t-J model which can be related to the well known bound states for one hole in the Neél-background. In addition to the one-particle properties we discuss the static spin-susceptibility in the paramagnetic state at the points and for different dopings . The -T-phase-diagram exhibits an enhanced stability of the antiferromagnetic state for electron-doped systems in comparison to hole-doped. This asymmetry in the phase diagram is in qualitative agreement with experiments for high-T_c materials. Received: 28 May 1998 / Revised and Accepted: 14 September 1998     

Low temperature thermodynamics of inverse square spin models in one dimension

We present a field-theoretic renormalization group calculation in two loop order for classical O(N)-models with an inverse square interaction in the vicinity of their lower critical dimensionality one. The magnetic susceptibility at low temperatures is shown to diverge like with a=(N-2)/(N-1) and . From a comparison with the exactly solvable Haldane-Shastry model we find that the same temperature dependence applies also to ferromagnetic quantum spin chains. Received: 20 February 1998 / Revised: 27 April 1998 / Accepted: 30 April 1998     

Magnetic susceptibility, heat capacity, and optical conductivity of LiPc and LiPcI

The magnetic susceptibility, using dc and electron spin resonance (ESR) methods, the specific heat, and the infrared properties of the one-dimensional molecular semiconductors lithium phthalocyanine (LiPc) and the iodinated compound LiPcI have been investigated for temperatures K. LiPc has a half-filled conduction band and is expected to be an organic metal. However, due to the strong Coulomb repulsion the system is a one-dimensional Mott-Hubbard insulator with a Hubbard gap of 0.75 eV as inferred from optical measurements. The localized electrons along the molecular stacks behave like a S = 1/2 antiferromagnetic spin chain. The spin susceptibility, as determined by ESR experiments, and the magnetic contribution to the heat capacity show a Bonner-Fisher type of behavior with an exchange constant K. LiPcI is an intrinsic narrow-gap semiconductor with an optical gap of 0.43 eV. In ESR experiments it is silent, indicating that all the unpaired electrons have been removed from the macrocycle via doping with iodine. Received: 16 June 1998 / Accepted: 14 July 1998     

An investigation of the 2D attractive Hubbard model

We present an investigation of the 2D attractive Hubbard model, considered as an effective model relevant to superconductivity in strongly interacting electron systems. We use both hybrid Monte-Carlo simulations and existing hopping parameter expansions to explore the low temperature domain. The increase of the static S-wave pair correlation with decreasing temperature, which depends weakly on the band filling in the explored temperature range, is analyzed in terms of an expected Kosterlitz-Thouless superconducting transition. Using both our data and previously published results, we show that the evidence for this transition is weak: If it exists, its temperature is very low. The number of unpaired electrons remains nearly constant with temperature at fixed attractive potential strength. In contrast, the static magnetic susceptibility decreases fast with temperature, and cannot be related only to pair formation. We introduce a method by which the Padé approximants of the existing series for the susceptibility give sensible results down to rather low temperature region, as shown by comparison with our numerical data. Received: 30 October 1996 / Revised: 23 October 1997 / Accepted: 29 January 1998     

Local environment and electron correlation effects on the magnetic properties of clusters

The electronic and magnetic properties of clusters are investigated in the framework of the Hubbard model by treating electron correlations effects in a saddle-point slave-boson approximation. The size dependent single-particle spectrum is calculated using a third moment real-space expansion of the local density of states. Results for the magnetic moments, magnetic order, average number of double occupations and hopping renormalizations are given as a function of the local coordination number z, for different representative values of the Coulomb interaction strength U/t and band filling n. Several transitions between paramagnetic, ferromagnetic and antiferromagnetic behaviors are obtained as a function of z. The environment dependence of the magnetic behavior and of the degree of electron delocalization is analyzed. Advantages and limitations of the present approach are discussed. Received: 8 January 1998 / Revised: 22 June 1998 / Accepted: 6 August 1998