Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:0808.3167 · Physical Review B

Low-energy theory of the t−t′−t″−U Hubbard model at half-filling: Interaction strengths in cuprate superconductors and an effective spin-only description of La2CuO4

Jean-Yves P. Delannoy, Michel J. P. Gingras, P. C. W. Holdsworth, A.–M. S. Tremblay

Published Jun 24, 2009• 57 citations
arXivHTMLPDF

Abstract

Spin-only descriptions of the half-filled one-band Hubbard model are relevant for a wide range of Mott insulators. In addition to the usual Heisenberg exchange, many other types of interactions, including ring exchange, appear in the effective Hamiltonian in the intermediate coupling regime. In order to improve on the quantitative description of magnetic excitations in the insulating antiferromagnetic phase of copper-oxide (cuprate) materials, and to be consistent with band-structure calculations and photoemission experiments on these systems, we include second- and third-neighbor hopping parameters, ${t}^{\ensuremath{'}}$ and ${t}^{\ensuremath{''}}$, into the Hubbard Hamiltonian. A unitary transformation method is used to find systematically the effective Hamiltonian and any operator in the spin-only representation. The results include all closed four-hop electronic pathways in the canonical transformation. The method generates many ring exchange terms that play an important role in the comparison with experiments on ${\text{La}}_{2}{\text{CuO}}_{4}$. Performing a spin-wave analysis, we calculate the magnon dispersion as a function of $U$, $t$, ${t}^{\ensuremath{'}}$, and ${t}^{\ensuremath{''}}$. The four parameters are estimated by fitting the magnon dispersion to the experimental results of Coldea et al. [Phys. Rev. Lett. 86, 5377 (2001)] for ${\text{La}}_{2}{\text{CuO}}_{4}$. The ring exchange terms are found essential, in particular to determine the relative sign of ${t}^{\ensuremath{'}}$ and ${t}^{\ensuremath{''}}$, with the values found in good agreement with independent theoretical and experimental estimates for other members of the cuprate family. The zero-temperature sublattice magnetization is calculated using these parameters and also found to be in good agreement with the experimental value estimated by Lee et al. [Phys. Rev. B 60, 3643 (1999)]. We find a value of the interaction strength $U\ensuremath{\simeq}8t$ consistent with Mott insulating behavior.

DOI
10.1103/physrevb.79.235130

Related papers

Néel order, ring exchange, and charge fluctuations in the half-filled Hubbard model

Jean-Yves P. Delannoy, Michel J. P. Gingras, P. C. W. Holdsworth +1 more

We investigate the ground state properties of the two-dimensional half-filled one band Hubbard model in the strong (large-$U$) to intermediate coupling limit (i.e., away from the strict Heisenberg limit) using an effective spin-only low-energy theory that includes nearest-neighbor exchange, ring exchange, and all other spin interactions to order $t{(t∕U)}^{3}$. We show that the operator for the staggered magnetization, transformed for use in the effective theory, differs from that for the order parameter of the spin model by a renormalization factor accounting for the increased charge fluctuations as $t∕U$ is increased from the $t∕U\ensuremath{\rightarrow}0$ Heisenberg limit. These charge fluctuations lead to an increase of the quantum fluctuations over and above those for an $S=1∕2$ antiferromagnet. The renormalization factor ensures that the zero temperature staggered moment for the Hubbard model is a monotonously decreasing function of $t∕U$, despite the fact that the moment of the spin Hamiltonian, which depends on transverse spin fluctuations only, in an increasing function of $t∕U$. We also comment on quantitative aspects of the $t∕U$ and $1∕S$ expansions.

arXiv:cond-mat/0412033•Physical Review B•Sep 16, 2005•53 citations

Magnetism in the Hubbard model: An effective spin Hamiltonian approach

Michael A. Tusch, Yolande H. Szczech, David E. Logan

We present an approach to the magnetic properties of the half-filled Hubbard model, based on an approximate mapping of its low-energy transverse spin excitations on to those of an effective underlying Heisenberg model, but with effective spin interactions which are self-consistently determined and not confined solely to nearest-neighbor couplings. The mapping is exact in strong-coupling and is found to be accurate over a very wide range of interaction strengths, down to weak coupling. At zero temperature, it permits ready evaluation at finite U of the one-loop effects of zero-point spin fluctuations on, e.g., the sublattice magnetization. At finite temperatures, thermodynamic properties of the system in the thermal paramagnet are studied via a physically transparent Onsager reaction field approach, which amounts to a self-consistent treatment of paramagnetic spin correlations. This is central not only in recovering the correct dimensionality dependence of antiferromagnetic long-ranged order, but also for the d=3 case of primary interest here yields a N\'eel temperature in close agreement with known strong- and weak-coupling limits. Spin correlation functions and magnetic susceptibilities also show very good agreement with quantum Monte Carlo calculations over an appreciable temperature range in which the low-lying transverse spin excitations are thermally dominant. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Mar 1, 1996•16 citations

Magnon spectra of cuprates beyond spin wave theory

Jiahui Bao, Matthias Gohlke, Jeffrey G. Rau +1 more

The usual starting point for understanding magnons in cuprate antiferromagnets such as La 2 CuO 4 is a spin model incorporating cyclic exchange, which descends from a one-band Hubbard model, and has parameters taken from fits based on non-interacting spin wave theory. Here we explore whether this provides a reliable description of experiment, using matrix product states (MPS) to calculate magnon spectra beyond spin wave theory. We find that analysis based on low orders of spin wave theory leads to systematic overestimates of exchange parameters, with corresponding errors in estimates of Hubbard t / U . Once these are corrected, the “standard” model provides a good account of magnon dispersion and lineshape in La 2 CuO 4 , but fails to fully capture the continuum observed at high energies. The extension of this analysis to CaCuO 2 and Sr 2 IrO 4 is also discussed.

arXiv:2407.20501•Physical Review Research•Mar 4, 2025•5 citations

Hubbard-model description of the high-energy spin-weight distribution in La 2 CuO 4

J. M. P. Carmelo, M. A. N. Araújo, Steven R. White +1 more

The spectral-weight distribution in recent neutron scattering experiments on the parent compound La${}_{2}$CuO${}_{4}$ (LCO), which are limited in energy range to about 450 meV, is studied in the framework of the Hubbard model on the square lattice with effective nearest-neighbor transfer integral $t$ and on-site repulsion $U$. Our study combines a number of numerical and theoretical approaches, including, in addition to standard treatments, density matrix renormalization group calculations for Hubbard cylinders and a suitable spinon approach for the spin excitations. The latter spin-$\frac{1}{2}$ spinons are the spins of the rotated electrons that singly occupy sites. These rotated electrons are mapped from the electrons by a uniquely defined unitary transformation, in which rotated-electron single and double occupancy are good quantum numbers for finite interaction values. Our results confirm that the $U/8t$ magnitude suitable to LCO corresponds to intermediate $U$ values smaller than the bandwidth $8t$, which we estimate to be $8t\ensuremath{\approx}2.36$ eV for $U/8t\ensuremath{\approx}0.76$. This confirms the unsuitability of the conventional linear spin-wave theory. Our theoretical studies provide evidence for the occurrence of ground-state $d$-wave spinon pairing in the half-filled Hubbard model on the square lattice. This pairing applies only to the rotated-electron spin degrees of freedom, but it could play a role in a possible electron $d$-wave pairing formation upon hole doping. We find that the higher-energy spin spectral weight extends to about 566 meV and is located at and near the momentum $[\ensuremath{\pi},\ensuremath{\pi}]$. The continuum weight energy-integrated intensity vanishes or is extremely small at momentum $[\ensuremath{\pi},0]$. This behavior of this intensity is consistent with that of the spin waves observed in recent high-energy neutron scattering experiments, which are damped at the momentum $[\ensuremath{\pi},0]$. We suggest that future LCO neutron scattering experiments scan the energies between 450 and 566 meV and momenta around $[\ensuremath{\pi},\ensuremath{\pi}]$.

arXiv:1209.1010•Physical Review B•Aug 20, 2012•8 citations

Effective Hamiltonian for cuprate superconductors

C. D. Batista, Armando A. Aligia

We compare numerical results of O and Cu occupations and optical and magnetic excitations for the three-band Hubbard model in a ${\mathrm{Cu}}_{4}$${\mathrm{O}}_{8}$ cluster, with those obtained with an effective Hamiltonian ${\mathit{H}}_{\mathit{e}\mathit{f}\mathit{f}}$ in which the Cu-O hopping ${\mathit{t}}_{\mathit{p}\mathit{d}}$ is eliminated by means of a canonical transformation. This transformation up to order ${\mathit{t}}_{\mathit{p}\mathit{d}}^{4}$ is not quantitatively valid when the relevant energy differences are of the order of \ensuremath{\surd}8 ${\mathit{t}}_{\mathit{p}\mathit{d}}$ or less. However, an ${\mathit{H}}_{\mathit{e}\mathit{f}\mathit{f}}$ with the form of the second-order one, but which retains all higher-order terms contained in a ${\mathrm{CuO}}_{4}$ (${\mathrm{Cu}}_{2}$${\mathrm{O}}_{7}$) cluster for zero or one added hole (one added electron) to the undoped system, leads to a very good agreement with the results using the full Hilbert space. This ${\mathit{H}}_{\mathit{e}\mathit{f}\mathit{f}}$ makes possible the numerical study of the three-band Hubbard model for realistic parameters, using larger clusters.

Physical review. B, Condensed matter•Apr 1, 1993•43 citations

Self-energy corrections in an antiferromagnet: Interplay of classical and quantum effects on quasiparticle dispersion

Pooja Srivastava, Avinash Singh

Self-energy corrections due to to fermion-magnon interaction are studied in the antiferromagnetic state of the $t\ensuremath{-}{t}^{\ensuremath{'}}\ensuremath{-}{t}^{\ensuremath{''}}$ Hubbard model within the noncrossing approximation in the full $U$ range from weak to strong coupling. The role of classical (mean-field) features of fermion and magnon dispersion, associated with finite $U,{t}^{\ensuremath{'}},{t}^{\ensuremath{''}}$, are examined on quantum corrections to quasiparticle energy, weight, one-particle density of states, etc. A finite-$U$ induced classical dispersion term, absent in the $t\ensuremath{-}J$ model, is found to play an important role in suppressing the quasiparticle weight for states near $\mathbf{k}=(0,0)$, as seen in cuprates. For intermediate $U$, the renormalized antiferromagnetic (AF) band gap is found to be nearly half of the classical value, and the weak-coupling limit is quite nontrivial due to strongly suppressed magnon amplitude. For finite ${t}^{\ensuremath{'}}$, the renormalized AF band gap is shown to vanish at a critical interaction strength ${U}_{c}$, yielding a spin fluctuation driven first-order AF insulator--paramagnetic metal transition. Quasiparticle dispersion evaluated with the same set of Hubbard model cuprate parameters, as obtained from a recent magnon spectrum fit, provides excellent agreement with angle-resolved photoemission spectroscopy data for ${\mathrm{Sr}}_{2}\mathrm{Cu}{\mathrm{O}}_{2}{\mathrm{Cl}}_{2}$.

arXiv:cond-mat/0402487•Physical Review B•Sep 13, 2004•9 citations
Spin-only descriptions of the half-filled one-band Hubbard model are relevant for a wide range of Mott insulators. In addition to the usual Heisenberg exchange, many other types of interactions, including ring exchange, appear in the effective Hamiltonian in the intermediate coupling regime. In order to improve on the quantitative description of magnetic excitations in the insulating antiferromagnetic phase of copper-oxide (cuprate) materials, and to be consistent with band-structure calculations and photoemission experiments on these systems, we include second- and third-neighbor hopping parameters, t\ensuremath′{t}^{\ensuremath{'}}t\ensuremath′ and t\ensuremath′′{t}^{\ensuremath{''}}t\ensuremath′′, into the Hubbard Hamiltonian. A unitary transformation method is used to find systematically the effective Hamiltonian and any operator in the spin-only representation. The results include all closed four-hop electronic pathways in the canonical transformation. The method generates many ring exchange terms that play an important role in the comparison with experiments on La2CuO4{\text{La}}_{2}{\text{CuO}}_{4}La2​CuO4​. Performing a spin-wave analysis, we calculate the magnon dispersion as a function of UUU, ttt, t\ensuremath′{t}^{\ensuremath{'}}t\ensuremath′, and t\ensuremath′′{t}^{\ensuremath{''}}t\ensuremath′′. The four parameters are estimated by fitting the magnon dispersion to the experimental results of Coldea et al. [Phys. Rev. Lett. 86, 5377 (2001)] for La2CuO4{\text{La}}_{2}{\text{CuO}}_{4}La2​CuO4​. The ring exchange terms are found essential, in particular to determine the relative sign of t\ensuremath′{t}^{\ensuremath{'}}t\ensuremath′ and t\ensuremath′′{t}^{\ensuremath{''}}t\ensuremath′′, with the values found in good agreement with independent theoretical and experimental estimates for other members of the cuprate family. The zero-temperature sublattice magnetization is calculated using these parameters and also found to be in good agreement with the experimental value estimated by Lee et al. [Phys. Rev. B 60, 3643 (1999)]. We find a value of the interaction strength U\ensuremath≃8tU\ensuremath{\simeq}8tU\ensuremath≃8t consistent with Mott insulating behavior.
We investigate the ground state properties of the two-dimensional half-filled one band Hubbard model in the strong (large-UUU) to intermediate coupling limit (i.e., away from the strict Heisenberg limit) using an effective spin-only low-energy theory that includes nearest-neighbor exchange, ring exchange, and all other spin interactions to order t(t∕U)3t{(t∕U)}^{3}t(t∕U)3. We show that the operator for the staggered magnetization, transformed for use in the effective theory, differs from that for the order parameter of the spin model by a renormalization factor accounting for the increased charge fluctuations as t∕Ut∕Ut∕U is increased from the $t∕U\ensuremath{\rightarrow}0$ Heisenberg limit. These charge fluctuations lead to an increase of the quantum fluctuations over and above those for an S=1∕2S=1∕2S=1∕2 antiferromagnet. The renormalization factor ensures that the zero temperature staggered moment for the Hubbard model is a monotonously decreasing function of t∕Ut∕Ut∕U, despite the fact that the moment of the spin Hamiltonian, which depends on transverse spin fluctuations only, in an increasing function of t∕Ut∕Ut∕U. We also comment on quantitative aspects of the t∕Ut∕Ut∕U and 1∕S1∕S1∕S expansions.
We present an approach to the magnetic properties of the half-filled Hubbard model, based on an approximate mapping of its low-energy transverse spin excitations on to those of an effective underlying Heisenberg model, but with effective spin interactions which are self-consistently determined and not confined solely to nearest-neighbor couplings. The mapping is exact in strong-coupling and is found to be accurate over a very wide range of interaction strengths, down to weak coupling. At zero temperature, it permits ready evaluation at finite U of the one-loop effects of zero-point spin fluctuations on, e.g., the sublattice magnetization. At finite temperatures, thermodynamic properties of the system in the thermal paramagnet are studied via a physically transparent Onsager reaction field approach, which amounts to a self-consistent treatment of paramagnetic spin correlations. This is central not only in recovering the correct dimensionality dependence of antiferromagnetic long-ranged order, but also for the d=3 case of primary interest here yields a N\'eel temperature in close agreement with known strong- and weak-coupling limits. Spin correlation functions and magnetic susceptibilities also show very good agreement with quantum Monte Carlo calculations over an appreciable temperature range in which the low-lying transverse spin excitations are thermally dominant. \textcopyright{} 1996 The American Physical Society.
The usual starting point for understanding magnons in cuprate antiferromagnets such as La 2 CuO 4 is a spin model incorporating cyclic exchange, which descends from a one-band Hubbard model, and has parameters taken from fits based on non-interacting spin wave theory. Here we explore whether this provides a reliable description of experiment, using matrix product states (MPS) to calculate magnon spectra beyond spin wave theory. We find that analysis based on low orders of spin wave theory leads to systematic overestimates of exchange parameters, with corresponding errors in estimates of Hubbard t / U . Once these are corrected, the “standard” model provides a good account of magnon dispersion and lineshape in La 2 CuO 4 , but fails to fully capture the continuum observed at high energies. The extension of this analysis to CaCuO 2 and Sr 2 IrO 4 is also discussed.
The spectral-weight distribution in recent neutron scattering experiments on the parent compound La2{}_{2}2​CuO4{}_{4}4​ (LCO), which are limited in energy range to about 450 meV, is studied in the framework of the Hubbard model on the square lattice with effective nearest-neighbor transfer integral ttt and on-site repulsion UUU. Our study combines a number of numerical and theoretical approaches, including, in addition to standard treatments, density matrix renormalization group calculations for Hubbard cylinders and a suitable spinon approach for the spin excitations. The latter spin-12\frac{1}{2}21​ spinons are the spins of the rotated electrons that singly occupy sites. These rotated electrons are mapped from the electrons by a uniquely defined unitary transformation, in which rotated-electron single and double occupancy are good quantum numbers for finite interaction values. Our results confirm that the U/8tU/8tU/8t magnitude suitable to LCO corresponds to intermediate UUU values smaller than the bandwidth 8t8t8t, which we estimate to be $8t\ensuremath{\approx}2.36$ eV for $U/8t\ensuremath{\approx}0.76$. This confirms the unsuitability of the conventional linear spin-wave theory. Our theoretical studies provide evidence for the occurrence of ground-state ddd-wave spinon pairing in the half-filled Hubbard model on the square lattice. This pairing applies only to the rotated-electron spin degrees of freedom, but it could play a role in a possible electron ddd-wave pairing formation upon hole doping. We find that the higher-energy spin spectral weight extends to about 566 meV and is located at and near the momentum $[\ensuremath{\pi},\ensuremath{\pi}]$. The continuum weight energy-integrated intensity vanishes or is extremely small at momentum $[\ensuremath{\pi},0]$. This behavior of this intensity is consistent with that of the spin waves observed in recent high-energy neutron scattering experiments, which are damped at the momentum $[\ensuremath{\pi},0]$. We suggest that future LCO neutron scattering experiments scan the energies between 450 and 566 meV and momenta around $[\ensuremath{\pi},\ensuremath{\pi}]$.
We compare numerical results of O and Cu occupations and optical and magnetic excitations for the three-band Hubbard model in a Cu4{\mathrm{Cu}}_{4}Cu4​O8{\mathrm{O}}_{8}O8​ cluster, with those obtained with an effective Hamiltonian Heff{\mathit{H}}_{\mathit{e}\mathit{f}\mathit{f}}Heff​ in which the Cu-O hopping tpd{\mathit{t}}_{\mathit{p}\mathit{d}}tpd​ is eliminated by means of a canonical transformation. This transformation up to order tpd4{\mathit{t}}_{\mathit{p}\mathit{d}}^{4}tpd4​ is not quantitatively valid when the relevant energy differences are of the order of \ensuremath{\surd}8 tpd{\mathit{t}}_{\mathit{p}\mathit{d}}tpd​ or less. However, an Heff{\mathit{H}}_{\mathit{e}\mathit{f}\mathit{f}}Heff​ with the form of the second-order one, but which retains all higher-order terms contained in a CuO4{\mathrm{CuO}}_{4}CuO4​ (Cu2{\mathrm{Cu}}_{2}Cu2​O7{\mathrm{O}}_{7}O7​) cluster for zero or one added hole (one added electron) to the undoped system, leads to a very good agreement with the results using the full Hilbert space. This Heff{\mathit{H}}_{\mathit{e}\mathit{f}\mathit{f}}Heff​ makes possible the numerical study of the three-band Hubbard model for realistic parameters, using larger clusters.
Self-energy corrections due to to fermion-magnon interaction are studied in the antiferromagnetic state of the $t\ensuremath{-}{t}^{\ensuremath{'}}\ensuremath{-}{t}^{\ensuremath{''}}$ Hubbard model within the noncrossing approximation in the full UUU range from weak to strong coupling. The role of classical (mean-field) features of fermion and magnon dispersion, associated with finite $U,{t}^{\ensuremath{'}},{t}^{\ensuremath{''}}$, are examined on quantum corrections to quasiparticle energy, weight, one-particle density of states, etc. A finite-UUU induced classical dispersion term, absent in the $t\ensuremath{-}J$ model, is found to play an important role in suppressing the quasiparticle weight for states near k=(0,0)\mathbf{k}=(0,0)k=(0,0), as seen in cuprates. For intermediate UUU, the renormalized antiferromagnetic (AF) band gap is found to be nearly half of the classical value, and the weak-coupling limit is quite nontrivial due to strongly suppressed magnon amplitude. For finite ${t}^{\ensuremath{'}}$, the renormalized AF band gap is shown to vanish at a critical interaction strength Uc{U}_{c}Uc​, yielding a spin fluctuation driven first-order AF insulator--paramagnetic metal transition. Quasiparticle dispersion evaluated with the same set of Hubbard model cuprate parameters, as obtained from a recent magnon spectrum fit, provides excellent agreement with angle-resolved photoemission spectroscopy data for Sr2CuO2Cl2{\mathrm{Sr}}_{2}\mathrm{Cu}{\mathrm{O}}_{2}{\mathrm{Cl}}_{2}Sr2​CuO2​Cl2​.