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:cond-mat/0309407 · Physical Review B

Variational cluster approach to spontaneous symmetry breaking: The itinerant antiferromagnet in two dimensions

Christopher Dahnken, Markus Aichhorn, Werner R. Hanke, Enrico Arrigoni, Michael Potthoff

Published Dec 10, 2004• 170 citations
arXivHTMLPDF

Abstract

Based on the self-energy-functional approach proposed recently [M. Potthoff, Eur. Phys. J. B 32, 429 (2003)], we present an extension of the cluster-perturbation theory to systems with spontaneously broken symmetry. Our method applies to models with local interactions and accounts for both short-range correlations and long-range order. Short-range correlations are accurately taken into account via exact diagonalization of finite clusters. Long-range order is described by variational optimization of a ficticious symmetry-breaking field. In comparison with related cluster methods, our approach is more flexible and, for a given cluster size, less demanding numerically, especially at zero temperature. An application of the method to the antiferromagnetic phase of the Hubbard model at half-filling shows good agreement with results from quantum Monte Carlo calculations. We demonstrate that the variational extension of the cluster-perturbation theory is crucial to reproduce salient features of the single-particle spectrum.

DOI
10.1103/physrevb.70.245110

Related papers

Charge ordering in extended Hubbard models: Variational cluster approach

Markus Aichhorn, Hans Gerd Evertz, Wolfgang von der Linden +1 more

We present a generalization of the recently proposed variational cluster perturbation theory to extended Hubbard models at half-filling with repulsive nearest neighbor interaction. The method takes into account short-range correlations correctly by the exact diagonalization of clusters of finite size, whereas long-range order beyond the size of the clusters is treated on a mean-field level. For one dimension, we show that quantum Monte Carlo and density-matrix renormalization-group results can be reproduced with very good accuracy. Moreover we apply the method to the two-dimensional extended Hubbard model on a square lattice. In contrast to the one-dimensional case, a first order phase transition between spin density wave phase and charge density wave phase is found as function of the nearest-neighbor interaction at onsite interactions $U\ensuremath{\geqslant}3t$. The single-particle spectral function is calculated for both the one-dimensional and the two-dimensional system.

arXiv:cond-mat/0402580•Physical Review B•Dec 3, 2004•99 citations

Cluster perturbation theory in Hubbard model exactly taking into account the short-range magnetic order in 2 × 2 cluster

Sergey V. Nikolaev, Sergei G. Ovchinnikov

The cluster perturbation theory is presented in the 2D Hubbard model constructed using X operators in the Hubbard-I approximation. The short-range magnetic order is taken into account by dividing the entire lattice into individual 2 × 2 clusters and solving the eigenvalue problem in an individual cluster using exact diagonalization taking into account all excited levels. The case of half-filling taking into account jumps between nearest neighbors is considered. As a result of numerical solution, a shadow zone is discovered in the quasiparticle spectrum. It is also found that a gap in the density of states in the quasiparticle spectrum at zero temperature exists for indefinitely small values of Coulomb repulsion parameter U and increases with this parameter. It is found that the presence of this gap in the spectrum is due to the formation of a short-range antiferromagnetic order. An analysis of the temperature evolution of the density of states shows that the metal-insulator transition occurs continuously. The existence of two characteristic energy scales at finite temperatures is demonstrated, the larger scale is associated with the formation of a pseudogap in the vicinity of the Fermi level, and the smaller scale is associated with the metal-insulator transition temperature. A peak in the density of states at the Fermi level, which is predicted in the dynamic mean field theory in the vicinity of the metal-insulator transition, is not observed.

Journal of Experimental and Theoretical Physics•Oct 1, 2010•15 citations

Two-particle self-consistent approach for broken symmetry phases

Lorenzo Del Re

Spontaneous symmetry breaking of interacting fermion systems constitutes a major challenge for many-body theory due to the proliferation of new independent scattering channels once absent or degenerate in the symmetric phase. One example is given by the ferro/antiferromagnetic broken symmetry phase (BSP) of the Hubbard model, where vertices in the spin-transverse and spin-longitudinal channels become independent with a consequent increase in the computational power for their calculation. We generalise the non-perturbative Two-Particle-Self-Consistent method (TPSC) to address broken SU(2) magnetic phases in the Hubbard model, offering an efficient approach that incorporates strong correlations. We show that in the BSP, the sum-rule enforcement of susceptibilities must be accompanied by a modified gap equation resulting in a renormalisation of the order parameter, vertex corrections and the preservation of the gap-less feature of the Goldstone modes. We then apply the theory to the antiferromagnetic phase of the Hubbard model in the cubic lattice at half-filling. We compare our results of double occupancies and staggered magnetisation to the ones obtained using Diagrammatic Monte Carlo showing excellent quantitative agreement. We demonstrate how vertex corrections play a central role in lowering the Higgs resonance with respect to the quasi-particle excitation gap in the spin-longitudinal susceptibility, yielding a well visible Higgs-mode.

arXiv:2312.16280•SciPost Physics•Mar 3, 2025•5 citations

Time-reversal symmetry breaking phase in the Hubbard model: A variational cluster approach study

Xiancong Lu, L. Chioncel, Enrico Arrigoni

We study the stability of the time-reversal symmetry breaking staggered-flux phase of a single band Hubbard model, within the variational cluster approach. For intermediate and small values of the interaction $U$, we find metastable solutions for the staggered-flux phase, with a maximum current per bond at $U\ensuremath{\approx}3.2$. However, allowing for antiferromagnetic and superconducting long-range order it turns out that in the region at and close to half filling the antiferromagnetic phase is the most favorable energetically. The effect of nearest-neighbor interaction is also considered. Our results show that a negative nearest-neighbor interaction and finite doping favors the stability of the staggered-flux phase. We also present preliminary results for the three-band Hubbard model obtained with a restricted set of variational parameters. For this case, no spontaneous time-reversal symmetry breaking phase is found in our calculations.

arXiv:1112.4562•Physical Review B•Mar 19, 2012•8 citations

Variational cluster approach to the Hubbard model: Phase-separation tendency and finite-size effects

Markus Aichhorn, Enrico Arrigoni, Michael Potthoff +1 more

Using the variational cluster approach, we study the transition from the antiferromagnetic to the superconducting phase of the two-dimensional Hubbard model at zero temperature. Our calculations are based on a method to evaluate the VCA grand potential which employs a modified Lanczos algorithm and avoids integrations over the real or imaginary frequency axis. Thereby, very accurate results are possible for cluster sizes not accessible to full diagonalization. This is important for an improved treatment of short-range correlations, including correlations between Cooper pairs in particular. We apply this improved method in order to investigate the cluster-size dependence of the phase-separation tendency that has been proposed recently on the basis of calculations for smaller clusters. While the energy barrier associated with phase separation rapidly decreases with increasing cluster size for both hole and electron doping, the extension of the phase-separation region behaves differently in the two cases. More specifically, our results suggest that phase separation remains persistent in the hole-doped case and disappears in the electron-doped case. We also study the evolution of the single-particle spectrum as a function of doping and point out the relevance of our results for experimental findings in electron and hole-doped materials.

arXiv:cond-mat/0607271•Physical Review B•Dec 22, 2006•89 citations

Variational Monte-Carlo Studies of Hubbard Model. II

Hisatoshi Yokoyama, Hiroyuki Shiba

As a continuation of a previous paper [J. Phys. Soc. Jpn. 56 (1987) 1490], the variational Monte-Carlo method is extended to include the antiferromagnetic long-range order. The theory is based on the Gutzwiller-type correlation factor and its effect is exactly taken into account by the Monte-Carlo procedure. An application is made to the half-filled-band case of one-dimensional lattice (50 sites), two-dimensional square lattice (up to 20×20 sites) and three-dimensional simple cubic lattice (6×6×6 sites). The result is qualitatively different from previous studies relying on the random-phase-type "Gutzwiller approximation." The variational energy for two and three dimensions is favorably compared with Hirsch's quantum Monte-Carlo data.

Journal of the Physical Society of Japan•Oct 15, 1987•153 citations
Based on the self-energy-functional approach proposed recently [M. Potthoff, Eur. Phys. J. B 32, 429 (2003)], we present an extension of the cluster-perturbation theory to systems with spontaneously broken symmetry. Our method applies to models with local interactions and accounts for both short-range correlations and long-range order. Short-range correlations are accurately taken into account via exact diagonalization of finite clusters. Long-range order is described by variational optimization of a ficticious symmetry-breaking field. In comparison with related cluster methods, our approach is more flexible and, for a given cluster size, less demanding numerically, especially at zero temperature. An application of the method to the antiferromagnetic phase of the Hubbard model at half-filling shows good agreement with results from quantum Monte Carlo calculations. We demonstrate that the variational extension of the cluster-perturbation theory is crucial to reproduce salient features of the single-particle spectrum.
We present a generalization of the recently proposed variational cluster perturbation theory to extended Hubbard models at half-filling with repulsive nearest neighbor interaction. The method takes into account short-range correlations correctly by the exact diagonalization of clusters of finite size, whereas long-range order beyond the size of the clusters is treated on a mean-field level. For one dimension, we show that quantum Monte Carlo and density-matrix renormalization-group results can be reproduced with very good accuracy. Moreover we apply the method to the two-dimensional extended Hubbard model on a square lattice. In contrast to the one-dimensional case, a first order phase transition between spin density wave phase and charge density wave phase is found as function of the nearest-neighbor interaction at onsite interactions $U\ensuremath{\geqslant}3t$. The single-particle spectral function is calculated for both the one-dimensional and the two-dimensional system.
The cluster perturbation theory is presented in the 2D Hubbard model constructed using X operators in the Hubbard-I approximation. The short-range magnetic order is taken into account by dividing the entire lattice into individual 2 × 2 clusters and solving the eigenvalue problem in an individual cluster using exact diagonalization taking into account all excited levels. The case of half-filling taking into account jumps between nearest neighbors is considered. As a result of numerical solution, a shadow zone is discovered in the quasiparticle spectrum. It is also found that a gap in the density of states in the quasiparticle spectrum at zero temperature exists for indefinitely small values of Coulomb repulsion parameter U and increases with this parameter. It is found that the presence of this gap in the spectrum is due to the formation of a short-range antiferromagnetic order. An analysis of the temperature evolution of the density of states shows that the metal-insulator transition occurs continuously. The existence of two characteristic energy scales at finite temperatures is demonstrated, the larger scale is associated with the formation of a pseudogap in the vicinity of the Fermi level, and the smaller scale is associated with the metal-insulator transition temperature. A peak in the density of states at the Fermi level, which is predicted in the dynamic mean field theory in the vicinity of the metal-insulator transition, is not observed.
Spontaneous symmetry breaking of interacting fermion systems constitutes a major challenge for many-body theory due to the proliferation of new independent scattering channels once absent or degenerate in the symmetric phase. One example is given by the ferro/antiferromagnetic broken symmetry phase (BSP) of the Hubbard model, where vertices in the spin-transverse and spin-longitudinal channels become independent with a consequent increase in the computational power for their calculation. We generalise the non-perturbative Two-Particle-Self-Consistent method (TPSC) to address broken SU(2) magnetic phases in the Hubbard model, offering an efficient approach that incorporates strong correlations. We show that in the BSP, the sum-rule enforcement of susceptibilities must be accompanied by a modified gap equation resulting in a renormalisation of the order parameter, vertex corrections and the preservation of the gap-less feature of the Goldstone modes. We then apply the theory to the antiferromagnetic phase of the Hubbard model in the cubic lattice at half-filling. We compare our results of double occupancies and staggered magnetisation to the ones obtained using Diagrammatic Monte Carlo showing excellent quantitative agreement. We demonstrate how vertex corrections play a central role in lowering the Higgs resonance with respect to the quasi-particle excitation gap in the spin-longitudinal susceptibility, yielding a well visible Higgs-mode.
We study the stability of the time-reversal symmetry breaking staggered-flux phase of a single band Hubbard model, within the variational cluster approach. For intermediate and small values of the interaction UUU, we find metastable solutions for the staggered-flux phase, with a maximum current per bond at $U\ensuremath{\approx}3.2$. However, allowing for antiferromagnetic and superconducting long-range order it turns out that in the region at and close to half filling the antiferromagnetic phase is the most favorable energetically. The effect of nearest-neighbor interaction is also considered. Our results show that a negative nearest-neighbor interaction and finite doping favors the stability of the staggered-flux phase. We also present preliminary results for the three-band Hubbard model obtained with a restricted set of variational parameters. For this case, no spontaneous time-reversal symmetry breaking phase is found in our calculations.
Using the variational cluster approach, we study the transition from the antiferromagnetic to the superconducting phase of the two-dimensional Hubbard model at zero temperature. Our calculations are based on a method to evaluate the VCA grand potential which employs a modified Lanczos algorithm and avoids integrations over the real or imaginary frequency axis. Thereby, very accurate results are possible for cluster sizes not accessible to full diagonalization. This is important for an improved treatment of short-range correlations, including correlations between Cooper pairs in particular. We apply this improved method in order to investigate the cluster-size dependence of the phase-separation tendency that has been proposed recently on the basis of calculations for smaller clusters. While the energy barrier associated with phase separation rapidly decreases with increasing cluster size for both hole and electron doping, the extension of the phase-separation region behaves differently in the two cases. More specifically, our results suggest that phase separation remains persistent in the hole-doped case and disappears in the electron-doped case. We also study the evolution of the single-particle spectrum as a function of doping and point out the relevance of our results for experimental findings in electron and hole-doped materials.
As a continuation of a previous paper [J. Phys. Soc. Jpn. 56 (1987) 1490], the variational Monte-Carlo method is extended to include the antiferromagnetic long-range order. The theory is based on the Gutzwiller-type correlation factor and its effect is exactly taken into account by the Monte-Carlo procedure. An application is made to the half-filled-band case of one-dimensional lattice (50 sites), two-dimensional square lattice (up to 20×20 sites) and three-dimensional simple cubic lattice (6×6×6 sites). The result is qualitatively different from previous studies relying on the random-phase-type "Gutzwiller approximation." The variational energy for two and three dimensions is favorably compared with Hirsch's quantum Monte-Carlo data.