Abstract:The two-particle self-consistent approach (TPSC) is a method for the one-band Hubbard model that can be both numerically efficient and reliable. However, TPSC fails to yield physical results deep in the renormalized classical regime of the bidimensional Hubbard model where the spin correlation length becomes exponentially large. We address the limitations of TPSC with improved approaches that we call TPSC+ and TPSC+SFM. In this work, we show that these improved methods satisfy the Mermin-Wagner theorem and the Pauli principle. We also show that they are valid in the renormalized classical regime of the 2D Hubbard model, where they recover a generalized Stoner criterion at zero temperature in the antiferromagnetic phase. We discuss some limitations of the TPSC+ approach with regards to the violation of the f-sum rule and conservation laws, which are solved within the TPSC+SFM framework. Finally, we benchmark the TPSC+ and TPSC+SFM approaches for the one-band Hubbard model in two dimensions and show how they have an overall better agreement with available diagrammatic Monte Carlo results than the original TPSC approach.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is that in the renormalized classical region of the two - dimensional Hubbard model, the two - particle self - consistent method (TPSC) fails to produce physical results. Specifically, when the spin - correlation length becomes very large, the TPSC method performs poorly in this region. To overcome this limitation, the authors propose two improved methods: TPSC+ and TPSC+SFM. ### Main problems 1. **Limitations of the TPSC method**: - The TPSC method fails in the renormalized classical region (i.e., the region where the spin - correlation length increases exponentially). - In particular, in the two - dimensional Hubbard model, the TPSC method cannot give correct physical results at low temperatures. 2. **Goals of the improved methods**: - **TPSC+**: Improve the TPSC method by introducing a partially dressed susceptibility so that it remains valid in the renormalized classical region. - **TPSC+SFM**: Further improve the TPSC+ method to solve the problems of violating the spin and charge conservation laws and the f - sum rule. ### Specific goals - **Satisfy the Mermin - Wagner theorem and the Pauli principle**: Ensure that the improved methods are physically reasonable. - **Restore the generalized Stoner criterion at zero temperature**: In the antiferromagnetic phase, the improved method should be able to restore this criterion. - **Improve the consistency with the existing diagrammatic Monte Carlo results**: In the weak - to - moderate coupling region, the improved method should show better consistency. ### Method overview - **TPSC+**: - Use the partially dressed susceptibility \(\chi^{(2)}\) to replace the non - interacting susceptibility \(\chi^{(1)}\). - Calculate the self - energy \(\Sigma^{(2)}\) and the Green's function \(G^{(2)}\) self - consistently to ensure the self - consistency of the method. - **TPSC+SFM**: - Further improve the partially dressed susceptibility to solve the problems of violating the spin and charge conservation laws and the f - sum rule in the TPSC+ method. - By choosing a specific form of the partially dressed susceptibility, ensure that the method is more physically reasonable. ### Mathematical expressions - **Definitions of self - energy and susceptibility**: - The expression of the self - energy \(\Sigma^{(2)}_{\sigma}(k, i\omega_n)\) is: \[ \Sigma^{(2)}_{\sigma}(k, i\omega_n) = U n_{-\sigma} + \frac{T}{N} \frac{U}{8} \sum_{q, i\nu_n} \left[ 3 U_{\text{sp}} \chi_{\text{sp}}(q, i\nu_n) + U_{\text{ch}} \chi_{\text{ch}}(q, i\nu_n) \right] G^{(1)}_{\sigma}(k + q, i\omega_n + i\nu_n). \] - The expression of the partially dressed spin susceptibility \(\chi^{(2)}_{\text{sp}}(q, i\nu_n)\) is: \[ \chi^{(2)}_{\text{sp}}(q, i\nu_n) = \frac{\chi^{(2)}(q, i\nu_n)}{1 - \frac{U_{\text{sp}}}{2} \chi^{(2)}(q, i\nu_n)}. \] - **Definition of spin - correlation length**: - In the renormalized classical region, the expression of the spin - correlation length \(\xi_{\text{sp}}\) is: \[ \xi_{\text{sp}} = \xi_0 \sqrt{\frac{U_{\text{sp}}}{2 \chi^{

An Improved Two-Particle Self-Consistent Approach

Improved Silbey-Harris Polaron Ansatz for the Spin-Boson Model

Improved effective vertices in the multi-orbital Two-Particle Self-Consistent method from Dynamical Mean-Field Theory

The Two-Particle Self-Consistent Approach for Multiorbital models: application to the Emery model

Multi-orbital two-particle self-consistent approach -- strengths and limitations

Two-particle self-consistent approach for broken symmetry phases

Magnetic correlations of a doped and frustrated Hubbard model: benchmarking the two-particle self-consistent theory against a quantum simulator

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Tangent Space Approach for Thermal Tensor Network Simulations of the 2D Hubbard Model

Benchmark Study of the Two-Dimensional Hubbard Model with Auxiliary-Field Quantum Monte Carlo Method

The (extended) dynamical mean field theory combined with the two-particle irreducible functional renormalization-group approach as a tool to study strongly-correlated systems

Beyond the instanton gas approach: dominant thimbles approximation for the Hubbard model

Covariant Bethe-Salpeter Approximation in Models of Strongly Correlated Electron Systems

Fermi surface evolution of the 2D Hubbard model within a novel four-pole approximation

A regularized second-order correlation method from Green's function theory

COM(3p) Solution of the 2D Hubbard Model: Momentum-Resolved Quantities

Testing the Monte Carlo–mean Field Approximation in the One-Band Hubbard Model

Benchmark of a modified Iterated Perturbation Theory approach on the 3d FCC lattice at strong coupling

Quasi-particle functional Renormalisation Group calculations in the two-dimensional t-t'-Hubbard model

Magnetic properties and superconductivity in the two-dimensional repulsive Hubbard model