اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدOn the dangers of partial diagrammatic summations: Benchmarks for the two-dimensional Hubbard model in the weak-coupling regime
109
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We study the two-dimensional Hubbard model in the weak-coupling regime and compare the self-energy obtained from various approximate diagrammatic schemes to the result of diagrammatic Monte Carlo simulations, which sum up all weak-coupling diagrams up to a given order. While dynamical mean-field theory provides a good approximation for the local part of the self-energy, including its frequency dependence, the partial summation of bubble and/or ladder diagrams typically yields worse results than second order perturbation theory. Even widely used self-consistent schemes such as GW or the fluctuation-exchange approximation (FLEX) are found to be unreliable. Combining the dynamical mean-field self-energy with the nonlocal component of GW in GW+DMFT yields improved results for the local self-energy and nonlocal self-energies of the correct order of magnitude, but here, too, a more reliable scheme is obtained by restricting the nonlocal contribution to the second order diagram. FLEX+DMFT is found to give accurate results in the low-density regime, but even worse results than FLEX near half-filling.
قيم البحث
اقرأ أيضاً
Partial bosonisation of the two-dimensional Hubbard model focuses the functional renormalisation flow on channels in which interactions become strong and local order sets in. We compare the momentum structure of the four-fermion vertex, obtained on t
he basis of a patching approximation, to an effective bosonic description. For parameters in the antiferromagnetic phase near the onset of local antiferromagnetic order, the interaction of the electrons is indeed well described by the exchange of collective bosonic degrees of freedom. The residual four-fermion vertex after the subtraction of the bosonic exchange contribution is small. We propose that similar partial bosonisation techniques can improve the accuracy of renormalisation flow studies also for the case of competing order.
We find that imposing the crossing symmetry in the iteration process considerably extends the range of convergence for solutions of the parquet equations for the Hubbard model. When the crossing symmetry is not imposed, the convergence of both simple
iteration and more complicated continuous loading (homotopy) methods are limited to high temperatures and weak interactions. We modify the algorithm to impose the crossing symmetry without increasing the computational complexity. We also imposed time reversal and a subset of the point group symmetries, but they did not further improve the convergence. We elaborate the details of the latency hiding scheme which can significantly improve the performance in the computational implementation. With these modifications, stable solutions for the parquet equations can be obtained by iteration more quickly even for values of the interaction that are a significant fraction of the bandwidth and for temperatures that are much smaller than the bandwidth. This may represent a crucial step towards the solution of two-particle field theories for correlated electron models.
We study a ferromagnetic instability in a doped single-band Hubbard model by means of dynamical mean-field theory with the continuous-time quantum Monte Carlo simulations. Examining the effect of the strong correlations in the system on the hypercubi
c and Bethe lattice, we find that the ferromagnetically ordered state appears in the former, while it does not in the latter. We also reveal that the ferromagnetic order is more stable in the case that the noninteracting DOS exhibits a slower decay in the high-energy region. The present results suggest that, in the strong-coupling regime, the high-energy part of DOS plays an essential role for the emergence of the ferromagnetically ordered state, in contrast to the Stoner criterion justified in the weak interaction limit.
Cooperation and competition between the antiferromagnetic, d-wave superconducting and Mott-insulating states are explored for the two-dimensional Hubbard model including nearest and next-nearest-neighbor hoppings at zero temperature. Using the variat
ional cluster approach with clusters of different shapes and sizes up to 10 sites, it is found that the doping-driven transition from a phase with microscopic coexistence of antiferromagnetism and superconductivity to a purely superconducting phase is discontinuous for strong interaction and accompanied by phase separation. At half-filling the system is in an antiferromagnetic Mott-insulating state with vanishing charge compressibility. Upon decreasing the interaction strength U below a certain critical value of roughly U=4 (in units of the nearest-neighbor hopping), however, the filling-dependent magnetic transition changes its character and becomes continuous. Phase separation or, more carefully, the tendency towards the formation of inhomogeneous states disappears. This critical value is in contrast to previous studies, where a much larger value was obtained. Moreover, we find that the system at half-filling undergoes the Mott transition from an insulator to a state with a finite charge compressibility at essentially the same value. The weakly correlated state at half-filling exhibits superconductivity microscopically admixed to the antiferromagnetic order. This scenario suggests a close relation between phase separation and the Mott-insulator physics.
We show that the numerically exact bold-line diagrammatic theory for the $2d$ Hubbard model exhibits a non-Fermi-liquid (NFL) strange metal state, which is connected to the SYK NFL in the strong-interaction limit. The solution for the doped system fe
atures the expected phenomenology with the NFL near half-filling at strong couplings and in a wide temperature range enclosed by the atomic state at high temperatures and a Fermi liquid at low temperatures. We demonstrate, however, that this behavior in the weakly doped regime is due to the unphysical branch of the Luttinger-Ward functional. On the other hand, our analysis shows that the NFL physics is realized at larger doping.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
