Arxiv:2105.03639V2 [Cond-Mat.Str-El] 17 May 2021 Share Similar Features in Their Phase Diagram
Total Page:16
File Type:pdf, Size:1020Kb
Pseudogap and superconductivity emerging from quantum magnetic fluctuations: a Monte Carlo study Weilun Jiang,1, 2 Yuzhi Liu,1, 2 Avraham Klein,3 Yuxuan Wang,4 Kai Sun,5 Andrey V. Chubukov,6 and Zi Yang Meng7, 1, ∗ 1Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 2School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100190, China 3Department of Physics, Faculty of Natural Sciences, Ariel University, Ariel, Israel 4Department of Physics, University of Florida, Gainesville, FL 32601 5Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA 6School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455, USA 7Department of Physics and HKU-UCAS Joint Institute of Theoretical and Computational Physics, The University of Hong Kong, Pokfulam Road, Hong Kong SAR, China (Dated: May 18, 2021) The origin of the pseudogap behavior, found in behavior around the putative quantum critical points many high-Tc superconductors, remains one of the (QCP). These features have led to the proposal that greatest puzzles in condensed matter physics. One soft quantum-critical fluctuations of the order parame- possible mechanism is fermionic incoherence, which ter serve as the source for the universal behavior and near a quantum critical point allows pair formation mediate singular interaction that gives rise to super- but suppresses superconductivity. Employing quan- conductivity with nontrivial pairing symmetry, strange tum Monte Carlo simulations of a model of itiner- metal behavior, and intertwined orders. ant fermions coupled to ferromagnetic spin fluctu- In many unconventional superconductors, most no- ations, represented by a quantum rotor, we report tably the cuprates, there is a third feature: the “pseu- numerical evidence of pseudogap behavior, emerg- dogap”, a gap-like feature in the fermionic spectrum ing from pairing fluctuations in a quantum-critical above the superconducting phase. Despite decades of non-Fermi liquid. Specifically, we observe enhanced investigation, the origin (or origins) of the pseudo- pairing fluctuations and a partial gap opening in gap remain intensely debated. One class of proposals the fermionic spectrum. However, the system re- names exotic, possibly topological order in the particle- mains non-superconducting until reaching a much hole channel as the origin [1–3], while another points to lower temperature. In the pseudogap regime the pairing fluctuations in the strong coupling regime [4– system displays a “gap-filling” rather than “gap- 7]. closing” behavior, consistent with experimental ob- It is clear, that understanding of the coupling be- servations. Our results provide the first unambigu- tween fermionic excitations near the Fermi surface ous lattice model realization of a pseudogap state in (FS) and bosonic quantum critical fluctuations [9–13] a strongly correlated system, driven by supercon- is needed to describe these three features. The develop- ducting fluctuations. ment of quantum Monte Carlo (QMC) algorithms for a class of models of this type has created a feasible Introduction way to study this physics in an unbiased manner [14– Even though unconventional and high-Tc superconduc- 16] . In QMC models, FS fermions couple to bosonic tivity arises in a diverse set of materials, many of them fluctuations, representing certain collective modes of arXiv:2105.03639v2 [cond-mat.str-el] 17 May 2021 share similar features in their phase diagram. One low-energy fermions [17–23]. The bosonic part is prominent feature is a superconducting (SC) dome, bestowed with independent (non-fermionic) dynamics which emerges near the termination point of a non-SC and can be tuned to criticality to mimic the situation in phase with either spin or charge order. The second fea- real materials. Crucially, these models are free of the ture is anomalous transport and non Fermi-liquid (nFL) sign-problem plaguing most fermionic QMC, allowing for a realistic test of theory and comparison with exper- iment. In this work, we consider a variant of a quantum crit- ∗ [email protected] ical model, in which the bosons represent critical fer- 2 FIG. 1. Model and Phase diagram. (a) Sketch of the model in Eq. (1). Deep yellow dots and grids of the top and bottom layers represent fermion degrees of freedom with nearest hopping strength t1 and next-nearest hopping t2. Blue arrows and grid in the middle layers denote bosonic parts with an unit vector representing θ 2 [0; 2π) of the rotor on each site. The on-site coupling K between fermions and bosons is shown by the vertical dashed lines. The system size is L × L. (b) U − T phase diagram of the model obtained from QMC simulation. The inset zooms in to the vicinity of the pseudogap (PG, yellow), ferromagnetic (FM, blue) and superconducting (SC, orange) regions. The blue points on the FM phase boundary are determined by finite size scaling with fixed T or U. Notably, for U = 5:9, as temperature gets lower, the system first enters into the FM phase at T ≈ 0:13, then exits it at T ≈ 0:08. The yellow points of the PG boundary are determined from the onset of a pseudogap in the single-particle spectrum, as shown in Fig.2. The red points denoting an onset of s-wave superconductivity are determined from the onset of a full gap in the spectrum as well as Kosterlitz-Thouless scaling of the pairing susceptibility. The maximum of superconducting phase transition temperature Tc is approximately 0.05. The green dot and dashed line, are the phase boundary of the (uncoupled) quantum rotor model [8]. romagnetic spin fluctuations (a “spin-fermion” model). a gap-like feature in the density of states (DOS). Unlike Compared to earlier works on ferromagnetism, the cou- in a BCS superconductor, the size of the gap remains pling strength of our model is stronger in two aspects. roughly independent on temperature, while the DOS First, the spin system is an XY quantum rotor model becomes progressively depleted (filled) upon lowering that is inherently more strongly fluctuating than an (raising) temperature. Importantly, the scaling behav- Ising model, analyzed earlier [24, 25]. Second, the cou- ior of the pairing susceptibility clearly shows that the pling constant K between the fermionic and bosonic system is not in a superconducting state. We thus iden- sectors is set to larger values than in previous works. tify the spectral gap in such a state as a pseudogap. As we see below, these modifications allow us to re- We note that the “gap-filling” behavior observed veal new behavior, including a pseudogap. in our numerical results has also been observed In the normal state, at low enough temperatures we in tunneling and photoemission experiments on the found in the bosonic sector near the QCP an over- cuprates [26], and has been obtained in a class of damped dynamics with linear frequency response (z = γ−models of quantum-critical pairing [6]. Our results, 2 scaling). This is different from the z = 3 behav- obtained from unbiased large-scale QMC simulations, ior, found in Ising systems, and is a result of a non- confirm the existence of a pseudogap behavior from conservation of the order parameter in our model. pairing fluctuations in a quantum-critical system with In the temperature range, where the bosonic sus- itinerant fermions. ceptibility is linear in frequency, we observed several The quantum-critical spin dynamics and normal remarkable features. The uniform susceptibility devi- state fermionic properties that we found are consis- ates from Curie-Weiss behavior and actually becomes tent with recent theoretical predictions for nFLs at fi- weaker at smaller T . In the fermionic sector, we found nite temperature, obtained within the modified Eliash- 3 berg theory [25, 27, 28]. This allows us to benchmark fluctuations driven pseudogap phenomena. our simulations and extract relevant parameters from the observables (see Supplementary Information (SI) Model for details). The onset temperature for pseudogap be- We consider a model of itinerant fermions coupled to havior, TPG, and superconducting Tc, extracted from SO(2) quantum rotors, as shown in Fig.1(a) (rotors QMC, are consistent with theoretical predictions (see are in the middle layer). The model is described by Ref.6 and Methods). Our results therefore provide the ^ ^ ^ ^ first unambiguous numerical realization of the transi- H = Hqr + Hf + Hqr-f; (1) tion from nFL to pseudogap and eventually to super- conductivity, lending support to the scenario of pairing where U X X H^ = L^2 − t cos θ^ − θ^ qr 2 i b i j i hi;ji ^ X + X + X Hf = −t1 c^iσλc^jσλ − t2 c^iσλc^jσλ − µ n^iσλ hi;jiσλ hhi;jiiσλ iσλ K X H^ = = − c^+ σxc^ · cos θ^ +c ^+ σyc^ · sin θ^ : (2) qr-f 2 iλ iλ i iλ iλ i iλ ^ The first term Hqr describes a quantum rotor model on This setup then allows us to analyze the universal be- a square lattice. Here L^i is the angular momentum of havior near a QCP with high numerical accuracy and ^ ^ large system sizes. 2D rotor θi at site i. The second term Hf describes two identical copies of spin-1=2 fermions on a square lattice, with layer index λ = 1 and 2 representing the Results top and the bottom layers. Fermions in each layer can We express all quantities in units of tb. In the simu- hop between nearest-neighbor (next-nearest-neighbor) lations we set K = 4, t1 = 1, t2 = 0:2 and µ = 0. sites with hopping amplitudes t1 (t2), and the chemical We varied U and the temperature T and constructed potential µ controls the fermion density. The last term the phase diagram of the model, Fig.1 (b), which fea- ^ Hqr-f couples quantum rotors and fermions via an on- tures a paramagnetic-ferromagnetic transition and sev- site ferromagnetic interaction that tends to align XY eral other transitions/phases.