Origin and fate of the pseudogap in the doped Hubbard model

The relationship between the pseudogap and underlying ground-state phases has not yet been rigorously established. We investigated the doped two-dimensional Hubbard model at finite temperature using controlled diagrammatic Monte Carlo calculations, allowing for the computation of spectral properties in the infinite-size limit and with arbitrary momentum resolution. We found three distinct regimes as a function of doping and interaction strength: a weakly correlated metal, a correlated metal with strong interaction effects, and a pseudogap regime at low doping. We show that the pseudogap forms both at weak coupling, when the magnetic correlation length is large, and at strong coupling, when it is shorter. As the temperature goes to zero, the pseudogap regime extrapolates precisely to the ordered stripe phase found by ground-state methods.

Symmetry-Broken Perturbation Theory to Large Orders in Antiferromagnetic Phases

We introduce a spin-symmetry-broken extension of the connected determinant algorithm [Riccardo Rossi, Determinant diagrammatic Monte Carlo algorithm in the thermodynamic limit, Phys. Rev. Lett. 119, 045701 (2017).PRLTAO0031-900710.1103/PhysRevLett.119.045701]. The resulting systematic perturbative expansions around an antiferromagnetic state allow for numerically exact calculations directly inside a magnetically ordered phase. We show new precise results for the magnetic phase diagram and thermodynamics of the three-dimensional cubic Hubbard model at half-filling. With detailed computations of the order parameter in the low to intermediate-coupling regime, we establish the Néel phase boundary. The critical behavior in its vicinity is shown to be compatible with the O(3) Heisenberg universality class. By determining the evolution of the entropy with decreasing temperature through the phase transition we identify the different physical regimes at U/t=4. We provide quantitative results for several thermodynamic quantities deep inside the antiferromagnetic dome up to large interaction strengths and investigate the crossover between the Slater and Heisenberg regimes.

Optimizing Superconductivity: From Cuprates via Nickelates to Palladates

Motivated by cuprate and nickelate superconductors, we perform a comprehensive study of the superconducting instability in the single-band Hubbard model. We calculate the spectrum and superconducting transition temperature T_{c} as a function of filling and Coulomb interaction for a range of hopping parameters, using the dynamical vertex approximation. We find the sweet spot for high T_{c} to be at intermediate coupling, moderate Fermi surface warping, and low hole doping. Combining these results with first principles calculations, neither nickelates nor cuprates are close to this optimum within the single-band description. Instead, we identify some palladates, notably RbSr_{2}PdO_{3} and A_{2}^{'}PdO_{2}Cl_{2} (A^{'}=Ba_{0.5}La_{0.5}), to be virtually optimal, while others, such as NdPdO_{2}, are too weakly correlated.

Evaluating Second-Order Phase Transitions with Diagrammatic Monte Carlo: Néel Transition in the Doped Three-Dimensional Hubbard Model

Diagrammatic Monte Carlo-the technique for the numerically exact summation of all Feynman diagrams to high orders-offers a unique unbiased probe of continuous phase transitions. Being formulated directly in the thermodynamic limit, the diagrammatic series is bound to diverge and is not resummable at the transition due to the nonanalyticity of physical observables. This enables the detection of the transition with controlled error bars from an analysis of the series coefficients alone, avoiding the challenge of evaluating physical observables near the transition. We demonstrate this technique by the example of the Néel transition in the 3D Hubbard model. At half filling and higher temperatures, the method matches the accuracy of state-of-the-art finite-size techniques, but surpasses it at low temperatures and allows us to map the phase diagram in the doped regime, where finite-size techniques struggle from the fermion sign problem. At low temperatures and sufficient doping, the transition to an incommensurate spin density wave state is observed.

Homotopic Action: A Pathway to Convergent Diagrammatic Theories

The major obstacle preventing Feynman diagrammatic expansions from accurately solving many-fermion systems in strongly correlated regimes is the series slow convergence or divergence problem. Several techniques have been proposed to address this issue: series resummation by conformal mapping, changing the nature of the starting point of the expansion by shifted action tools, and applying the homotopy analysis method to the Dyson-Schwinger equation. They emerge as dissimilar mathematical procedures aimed at different aspects of the problem. The proposed homotopic action offers a universal and systematic framework for unifying the existing-and generating new-methods and ideas to formulate a physical system in terms of a convergent diagrammatic series. It eliminates the need for resummation, allows one to introduce effective interactions, enables a controlled ultraviolet regularization of continuous-space theories, and reduces the intrinsic polynomial complexity of the diagrammatic Monte Carlo method. We illustrate this approach by an application to the Hubbard model.