Topological Nature of High Temperature Superconductivity

M. C. Diamantini,1 C. A. Trugenberger,2 and V.M. Vinokur3, 4 1NiPS Laboratory, INFN and Dipartimento di Fisica e Geologia, University of Perugia, via A. Pascoli, I-06100 Perugia, Italy 2SwissScientific Technologies SA, rue du Rhone 59, CH-1204 Geneva, Switzerland 3Terra Quantum AG, St. Gallerstrasse 16A, CH-9400 Rorschach, Switzerland 4University of Twente, Faculty of Science and Technology and MESA+ Institute of Nanotechnology, 7500 AE Enschede, The Netherlands The key to unraveling the nature of high-temperature superconductivity (HTS) lies in resolving the enigma of the pseudogap state. The pseudogap state in the underdoped region is a distinct thermodynamic phase character- ized by nematicity, temperature-quadratic resistive behavior, and magnetoelectric effects. Till present, a general description of the observed universal features of the pseudogap phase and their connection with HTS was lack- ing. The proposed work constructs a unifying effective field theory capturing all universal characteristics of HTS materials and explaining the observed phase diagram. The pseudogap state is established to be a phase where a charged magnetic monopole condensate confines Cooper pairs to form an oblique version of a superinsulator. The HTS phase diagram is dominated by a tricritical point (TCP) at which the first order transition between a fundamental Cooper pair condensate and a charged magnetic monopole condensate merges with the continuous superconductor-normal metal and superconductor-pseudogap state phase transitions. The universality of the HTS phase diagram reflects a unique topological mechanism of competition between the magnetic monopole condensate, inherent to antiferromagnetic-order-induced Mott insulators and the Cooper pair condensate. The obtained results establish the topological nature of the HTS and provide a platform for devising materials with the enhanced superconducting transition temperature.

INTRODUCTION emergence of the magnetoelectric (Kerr) effect, see [3,7–9] and references therein, evidencing that the PG state is a dis- The origin and underlying mechanism of high-temperature tinct thermodynamic phase. Recent shot noise measurements superconductivity (HTS) has been remaining a subject of deep in the PG region detected pairing of the charge carri- heated debate for more than three decades after its discovery. ers [4]. Finally, the evolution from the Mott insulator at p = 0 Even the metallic state above the superconducting transition to superconductivity upon increasing p suggests a quantum temperature Tc seem to differ strikingly from metals in which superconductor-insulator transition (SIT) [10] at the putative ∗ ∗ conventional superconductors transform [1–5]. In the under- transition point p where the T (p) line continued into the doped regime of HTSs the pseudogap state (PG) opens up dome hits the p-axis, see Fig.1. ∗ at the temperature T much above Tc, see [5] and references therein. The PG state exhibits striking properties including Building on the idea that it is quantum criticality associ- nematicity [3], temperature-quadratic resistive behavior, mag- ated with the SIT that dominates the physics of HTS and netoelectric effects [1], and carries a signature of pairing well neighboring phases, we adopt the machinery of the topologi- above Tc [4]. These exotic behaviors suggest that understand- cal gauge theory of the SIT [12, 25] for the description of the ing the nature of HTS requires gaining a fundamental insight HTS phase diagram. Deriving from the established Fermi- 2 into the nature of the PG state. liquid-like R ∝ T behavior of the resistance in the PG regime, combined with the fact that the dominant charge car- The behavior of HTS materials is graphically summarized riers in the pseudogap state carry the charge 2e, we demon- in the phase diagram in Fig.1. At very low doping, strong strate that the PG regime represents a novel state, the oblique electron interactions in HTS, responsible for strong Cooper superinsulator, in which the electric conductance is mediated pairing [5], localize the antiferromegnetic state into a Mott in- by symmetry-protected surface fermionic modes. This en- sulator [1,6,7]. Further doping, i.e. adding p holes per Cu ables us to construct a theory that encompasses all the ob- atom, turns the electrons itinerant. Thus, the superconduc- served features of the PG and HTS phases within a unified tivity dome forms in the p ≈ 0.05 ÷ 0.3 interval, followed universal description. We show that the key ingredient of the arXiv:2009.01763v2 [cond-mat.supr-con] 13 Feb 2021 by a Fermi liquid at higher doping, see phase diagram in HTS mechanism is the formation of a condensate of dyons, Fig.1. The elusive PG state emerges in the underdoped region. particles carrying both electric and magnetic charge, which Above the temperature T ∗, associated with the onset of the appear as topologically-charged magnetic monopoles at the

PG, a metallic sheet resistance linear in temperature, R ∝ T, end of open vortices playing the role of Dirac strings [23, 24]. is observed. At the same time, below temperature, T ∗∗ < T ∗, The superconducting dome of HTS appears as a coexistence 2 a distinct switch to a quadratic R ∝ T dependence occurs, phase of Cooper pair and dyon condensates. The PG state holding till superconducting fluctuations set in, see [1,6] and above Tc in the underdoped region is a pure dyon conden- references therein. The full switching on of the PG state at sate characterized by magnetoelectric and nematic effects, in T=T ∗∗ is accompanied by a nematic phase transition and the a perfect agreement with the experiment. In this state there 2 are symmetry-protected fermionic surface dyons of charge 2e The electromagnetic response of superinsulators is de- which live on a bubble-like texture and, as a result, cause the scribed [25] by the compact QED [21] in the phase where observed quadratic in T resistance. The high transition tem- magnetic monopoles [22] form a Bose condensate [23, 24]. perature of the HTS is explained by the topological obstruc- The inclusion of the θ-term in the compact QED theory re- tion of splitting Cooper pairs into single electrons in pres- sults in the Euclidean lattice gauge model [28] ence of magnetic monopoles of half the fundamental magnetic Z X −S charge, i.e. by the Dirac quantization condition. Z = DAµe ,

{nµ},{mµ} X 1  2  θ  S = F − 2πS + iA n + m , (2) THE MODEL AND FREE ENERGY OF THE DYONIC 4 f 2 µν µν µ µ 2π µ STATE x

where the integers nµ and mµ = (1/2)KµαβS αβ are the conserved charge and magnetic monopole currents, respec- We start constructing our model with stating its two key tively, Kµαβ is the lattice BF operator [12] (see Methods), and ingredients, the relevance of the SIT for the HTS mecha- the dimensionless coupling f is an effective strength of the nism and the established Fermi liquid-like behaviour of the Coulomb interaction in the material. Magnetic monopoles are charge carriers in the PG state [1], which implies that they are topological excitations arising in the electromagnetic response symmetry-protected surface states. Three possible states in- due to the compact nature of the electromagnetic gauge poten- volved in the SIT are superconductors, bosonic topological tial. Equation (2) implies that, in the presence of the topolog- insulators (BTI) [15] and superinsulators [12, 20, 25]. Super- ical θ-term, magnetic monopoles acquire an electric charge conductors do not carry surface states. In three dimensions, θ/2π and become dyons. This is known as Witten effect [29]. the BTI, in which neither time-reversal symmetry breaking Having in mind Cooper pairs, we let the electric charge unit nor topological order occur, also cannot provide the neces- be 2e. Dirac quantization qg/(4π~) = n/2, with q and g being sary fermionic surface degrees of freedom [19]. The PG state electric and magnetic charges, and an integer n [22], requires should thus be associated with superinsulators. Recalling the unit of magnetic charge be π/e. that the pseudogap state exhibits also the magnetoelectric ef- fect [1], we conclude that the Lagrangian describing the elec- Integrating out the gauge field Aµ, one obtains the model tromagnetic response of the system should include the so- formulated solely in terms of charge and monopole currents, called θ-term, Lθ. As we now demonstrate, the θ-term [27] X − transforms the standard superinsulator into a new state, which Z = e S , (3) we call oblique superinsulator, and which does indeed host {nµ},{mµ} symmetry-protected fermionic surface dyons of charge 2e. X  θ  1  θ  π2 X 1 S = f 2 n + m n + m + m m µ 2π µ −∇2 µ 2π µ f 2 µ −∇2 µ The θ-term is given by the expression x x 1 θ θ +iπ nµ Kµαβ Mαβ(4), L = F µναβF = E · B , (1) −∇2 θ 32π2 µν αβ 4π2 where the integers Mµν = (1/2)S µµναβS αβ are such that where Fµν = ∂µAν − ∂νAµ is the electromagnetic tensor, mµ = ∆ν Mµν (S µ is the lattice shift operator, see Methods). 2 Aν = (ϕ, −A) is the gauge potential, (θ/4π )δi j represents the The last term represents the topological interaction between pseudoscalar part of the magnetoelectric coupling, or polariz- charges and monopoles. When the Dirac quantization condi- ability αi j ≡ [∂M j/∂Ei]B=0 ≡ [∂Pi/∂B j]E=0, E = −∇ϕ − ∂tA tion is satisfied, it falls out of the partition function since it and B = ∇ × A are the electric and magnetic fields, respec- is always an integer multiple of 2π, as can be easily seen by tively, P and M are the respective polarization and magne- representing the conserved charge current as nµ = KµαβXαβ, tization. Here we use natural units c=1, with c being the with the gauge conditions ∆αXαβ = ∆βXαβ = 0, so that Xαβ light velocity in the material, ~=1, ε0=1, Greek letters de- contains four degrees of freedom as nµ, and using the relation note space-time indices and a sum is implied over repeating (16) in Methods. indices. Since the θ-term can be presented as a full deriva- µναβ tive, ∝ ∂µ( Aν∂αAβ), it is a pure surface term on bounded Keeping only self-interaction terms, one finds that the ac- spaces. On a compact Euclidean torus T4, the integral of the θ- tion associated with a particle is proportional to the length N term is a topological invariant contributing a factor exp(inθ) to of its Euclidean world-line, representing a lattice string [30]. the partition function, hence θ is 2π-periodic. The θ-term de- As usual in a statistical field theory, this Euclidean action is scribes the topological coupling of a vector and a pseudovec- equivalent to the “energy” of an equivalent statistical model tor field, therefore, it is odd under parity and time-reversal in one more spatial dimension, where the coupling constant symmetries. The values of θ compatible with the time-reversal plays the role of “temperature”. Since the entropy of a lattice symmetry are thus θ = 0 (mod 2π) and θ = π (mod 2π). string is also proportional to its length, the “free energy” of 3 a particle with the electric charge n and magnetic charge m Λ becomes / 7 )

350 �

( 6 "  2 2 # ∗ � θ π � ! Oblique F f 2G n m G m2 − µ N , � 5 = (0) + + 2 (0) (5) 300 SI 2π f 4 3 MC + CP CP + MC where G(0) is the value of the lattice Coulomb potential at 250 condensate condensate ) 2 µ ≈ � coinciding points and ln7, since at each step the non- ( 0.8 1.0 1.2 1.4 � �/�" backtracking string has to choose among 7 possible directions 200 �∗∗ to continue. if the factor in the brackets is positive, the free energy is minimal for N = 0. If, instead, it is negative, the 150 free energy is minimized by N = ∞. This means that par- PG ticles with quantum numbers resulting in the positive factor 100 in brackets in Eq. (5) are suppressed and exist only as short- living fluctuations, while particles with quantum numbers for 50 which the factor in brackets is negative form Bose conden- AF SC ∗ sates. If both Bose condensates are possible, the one with the � � 0 lowest free energy is stable, but the phases may coexist. The 0 0.05 0.1 0.15 0.2 0.25 condensation condition � (ℎ��� �������) f 2  θ 2 π η n m η m2 < , FIG. 1. A sketch of the temperature-doping phase diagram of + + 2 1 (6) π 2π f HTS. The phase diagram presents different phases, including anti- ferromagnetism (blue), superconductivity (purple to dark blue), and where η = πG(0)/µ, describes the interior of a tilted ellipse pseudogap (green) domains. The dashed gray line at T ∗ marks the p 2 p p 2 p with the semiaxes rn= π/ f 1/η and rm= f /π 1/η on onset of the pseudogap phase detected by deviation from the linear ∗∗ an integer lattice of electric and magnetic charges. The quan- R ∝ T resistance dependence. The gray dashed-point line, T , √ 2 tity 1/ η defines the overall scale of the ellipse, while the indicates the emergence of square R ∝ T resistance dependence, 2 nematicity, and magnetoelectric effects. The pc marks the tricritical coupling f /π is the ratio between its semiaxes. Varying the ∗ ∗ coupling f and angle θ leads to complex phase structures con- point at the dome maximum, and p , where T hits the x-axis, is the position of the putative quantum Mott-insulator-superconductor sisting of alternating sequences of three possible phases [28], phase transition at T = 0. Inset: The calculated phase transition a superconducting phase consisting of a pure charge conden- lines in the vicinity of the tricritical point. sate (m=0), a Coulomb phase with no condensate, and an oblique confinement phase, where condensed particles carry both electric and magnetic charges [31]. not necessarily unity. Only the shift θ → θ + 4π, accompanied by n → n − 2m, yields the same charge and statistics. Had charges been fermions, then the correct periodicity, indeed, THE PHASE DIAGRAM AND HTS MECHANISM would have been θ → θ + 2π, but fermions do not Bose con- dense. This statistical Witten effect [32] is the same reason The seminal reference [28] restricted its analysis to θ ∈ why strong bosonic topological insulators are also periodic [0, 2π[ because of the periodicity of the charge spectrum un- only under shifts θ → θ + 4π, contrary to fermionic ones [19]. der shifts θ→θ + 2π. This periodicity is evident in Eqs. [5,6], since the shift θ→θ + 2π can be compensated by the corre- sponding shift n→n − m. As a consequence of this restric- Since the correct periodicity in the angle θ is 4π, and not tion, the superconducting and oblique confinement phases can 2π, the case θ = 2π is generically different from that with never be adjacent, but are always separated by the Coulomb θ = 0. For θ = 2π the (n, m) lattice is not tilted, but all phases [28]. Importantly, however, the correct periodicity in θ m = const , 0 lines are shifted in their n-values with respect is 4π. While the charge of dyonic composites is indeed peri- to the x-axis m = 0, as shown in Fig.2. Cranking down the odic under shifts θ → θ + 2π, their statistics is not, if charges coupling π/ f 2 from high (corresponding to weak Coulomb are bosons, as has been pointed out in [32]. The statistics of interactions) to small (corresponding to strong Coulomb in- a composite of n bosonic electric charges and m magnetic teractions) values, first elongates the ellipse along the x-axis, monopoles with θ = 0 is (−1)nm, where +1 corresponds to then passes it through a circle, and, finally, elongates it along bosons and −1 to fermions [33, 34]. If we turn on the an- the y-axis. At η > 1, the nearly circles forming at the mid-way gle θ, the total electric charge changes to q = n + mθ/2π contain only the origin of the lattice, which corresponds to the but the statistics remains unchanged [35]. If we compensate Coulomb insulating phase. At η < 1, instead, the deformation the shift θ → θ + 2π by the corresponding transformation processes from horizontally elongated ellipses containing the n → n − m, to leave the total charge unchanged, we change lattice point (n = 1, m = 0) and describing the superconduc- the statistics of the composite by the factor (−1)m2 , which is tor, Fig.2a, through a circle containing both kinds of lattice 4

� = 0 � < 1 a � b � c �

(1,1) (2,1) (1,1) (2,1) (1,1) (2,1)

(−1,0) (0,0) (1,0) � (−1,0) (0,0) (1,0) � (−1,0) (0,0) (1,0) �

(−1, −1) (−1, −1) (−1, −1)

FIG. 2. The ellipse technique for calculating emerging phases at T = 0 and η < 1. The x-axis shows the integer dyonic charges q = m + n, for θ = 2π, while the y-axis lists integer magnetic charges m. a: Superconducting state, where only points n = ±1 (electric charge) and m = 0 (blue circles) fall inside the ellipse. b: First order phase transition between the superconductor and oblique superinsulator as the ellipse contains both superconducting and dyonic unit charge (q = ±1, m = ±1) integer points (red rim green circles). c: Oblique superinsulator phase.

� all scale of the ellipse shrinks, see Fig.3, by the scale factor S (T) shown in the inset in Fig.3, see Methods. Thus, lattice points that were within the ellipse at T = 0 fall out as the (1,1) (2,1) temperature grows and the system crosses from the under-the- dome HTS coexistence phase to the pseudogap- and then to a metallic phase as illustrated in Fig.3.

� (−1,0) (0,0) (1,0) As we have already mentioned, the phase harboring the 1.0 magnetic monopole condensate is a superinsulator [12, 20, 25, 0.9 0.8 26]. We, therefore, name the material hosting an oblique con-

) 0.7 finement phase an oblique superinsulator. The numerically (−1, −1) � 0.6 S( 0.5 computed phase diagram in the vicinity of the tricritical point 0.4 is shown in the inset in Fig.1, where we have identified the 0 2 4 6 8 10 coupling constant π/ f 2 with the doping normalized to unity at �!�/Λ the Tc maximum. The values of Tc = O(100) K are obtained

FIG. 3. Temperature evolution and deconfinement at p 1, the boundary of the superconducting coexist. The particles in both these coexisting condensates dome marks the region where the magnetic monopoles cease carry unit charge, but for the ‘superconducting’ particles the to exist. Above this phase boundary, the Cooper pairs are free charge is fundamental, while ‘oblique confinement’ carriers to split into single electrons to form a normal metal like in an acquire charge due to the Witten effect. Since G(0) = 0.155 ordinary superconductor. In the underdoped region p/pcr < 1, and µ ≈ ln7, we have η=0.25<1. It is exactly this coexis- the phase boundary marks the temperature at which the funda- tence regime that is realized, without an intervening normal mental charge condensate vanishes, but the charged monopole insulating phase. Upon increasing the temperature, the over- condensate survives at higher temperatures. This implies that 5 the oblique superinsulator forms above the dome. We iden- that this superconducting transition temperature Tc should co- tify this state as the pseudogap state that sets in between the incide with the temperature of the formation of the monopole ∗ dome and the temperature T [1]. Note that the charge con- condensate, TM. However, a precise self-consistent deriva- densate disappears at T > Tc not because of the collapse of the tion of the enhanced transition temperature Tc requires a mi- coexisting monopole-Cooper pair mixture to the true ground croscopic theory of Cooper pairing in the presence of half- state, as one would have expected in the case of the adiabatic magnetic monopoles, which goes beyond the scope of this over-cooling transition, but because Cooper pairs excitations work and is a subject of a forthcoming publication. Note that now exist in an out-of-condensate state acquire a large finally that, at variance to what occurs in 3D HTS systems, mass. In other words, Cooper pairs viewed as charge excita- magnetic monopoles in 2D are not real particles but rather tun- tions acquire an energy gap in their excitation spectrum above neling events, instantons, which do not impose any topolog- Tc, hence do not contribute to transport in precise analogy to ical obstruction on Cooper pair breaking. This explains why usual electrons in a usual insulator. In the proposed model, even a monolayer of HTS material has essentially the same the physics of the high-Tc cuprates is dominated by the tri- transition temperature Tc as the bulk sample [18] in stark con- critical point at p/pc = 1 and kBTc/Λ = 4.9, where the first- trast to behaviors of conventional 2D superconducting films order transition within the superconducting dome meets two which do not harbor real 3D monopoles. As a result, con- continuous transitions to a normal metal (in the overdoped re- ventional superconducting films, even hosting the SIT and the gion) and to the pseudogap state (in the underdoped region). superinsulating state, demonstrate a significantly reduced Tc This universal structure is an implication of the presence of compared to that of the bulk samples of the same materials. the magnetic monopole condensate in cuprates. Furthermore, monopoles ensure that in the underdoped region the charge THE NATURE OF THE PSEUDOGAP STATE carriers appear as paired well above Tc [4]: the out of con- densate Cooper pairs are protected from splitting into single electrons by the Dirac quantization condition. In a superconductor, the electromagnetic response is me- −2 µ diated by photons with mass (1/2)λL AµA , leading to the µ −2 The key ingredients of HTS and its elevated transition tem- London equations j = λL Aµ. The corresponding elec- perature emerging from our model are magnetic monopoles tromagnetic response of an oblique superinsulator has been and the magnetoelectric effect encoded in the θ-term. Both of derived in [36, 37]. As in superconductors, the condensate them are of the topological origin and both are prominently generates a massless mode that combines with the photon 3D effects. A word of caution is, however, needed at this to produce a massive excitation. Contrary to superconduc- point. When discussing dimensionality effects in HTS one tors, however, this massless mode is not a scalar represent- must be very careful to specify precisely which effect one is ing the phase of the charge order parameter but, rather, the addressing. It is well established, e.g. that the electric re- “dual magnetic phase”, embedded into an antisymmetric ten- sponse in the fluctuation region is described by the 2D be- sor gauge field Bµν, with the three-tensor field strength Hµνρ = havior of the resistivity and the Hall coefficient, see [41, 42] ∂µBνρ + ∂νBρµ + ∂ρBµν invariant under the gauge symmetries and references therein, reflecting the dominant quasi-2D na- of the second kind, Bµν → Bµν + ∂µλν − ∂νλµ. Since B0i are ture of the macroscopic electric response of cuprates. At Lagrange multipliers with no time derivatives appearing in the the same time, the topological effects are of the 3D nature. kinetic term Hµνρ, and the gauge function λµ is itself defined This is perfectly consistent with the fact that the unit cell only up to a derivative, the antisymmetric tensor contains in- of the HTS materials contains two neighboring CuO planes, deed a single degree of freedom. The full electromagnetic which makes it possible for the Cooper pairs not only to response of the oblique superinsulator is move along the ab-planes but also to tunnel between them 1    µν  in the c-direction. Therefore the HTS materials, even mono- L = − Bµν + Fµν B + Fµν + 4 f 2 layers containing a single unit cell, are essentially 3D with respect to monopole formation and host real 3D quantized θ   µναβ   1 µνα Dirac monopoles. These monopoles carry a half of the funda- + Bµν + Fµν  Bαβ + Fµν + HµναH . (7) mental magnetic charge since the corresponding unit charge 32π2 12Λ2 is the charge of the Cooper pairs. A single electron, in the Since fields enter only in the combination Bµν+Fµν, the gauge presence of these half monopoles, would have acquired an ad- invariance of the second kind of Bµν is ensured by the simul- ditional spin S = ~(eg/4π) = ~/4, which is impossible due taneous shift Aµ→Aµ−λµ, leading to three degrees of freedom to the non-Abelian nature of the rotation group in 3D. As a of mass, see Methods, result, the emergence of single electrons is prevented by the s !2 energy barrier associated with this topological restriction and f Λ 4π  θ 2 m = + . (8) is suppressed in the superconducting phase of HTS materi- θ 4π f 2 π als. This topological obstruction ties Tc self-consistently to the condensation temperature of magnetic monopoles, inde- Two of them belong in the original gauge field and one is pendently of the pairing gap. At this point, we can conjecture from the new field due to the monopole condensate. Alterna- 6 tively, the original electromagnetic field tensor Fµν can be re- coordinate and which one is the momentum). Setting the absorbed by the new field in a generalized Stuckelberg¨ mech- Lagrange multiplier λ0=0 (Weyl gauge) after recording the anism [36] dual to the familiar Anderson-Higgs mechanism. Gauss law constraint b=−Φ0 it implements, gives The Stuckelberg¨ mechanism is the Anderson-Higgs mecha- M nism without a Higgs field or, in other words, the limit of 1 2 2 2 H∂M = πξ + v (∂iξ) , (12) the Anderson-Higgs mechanism in which the Higgs field be- 2M 2 comes infinitely heavy and predates both the Anderson and which describes a massless degree of freedom as can be seen   the Higgs versions. That the phase of the order parameter and 2 2 P 2 from the Hamiltonian equations of motion ∂0 − v i ∂i ξ = ∂ ϕ−A the gauge field appear as a combination ( µ µ), means that 0. The gauge field λµ is not a dynamical field because there it is the Stuckelberg¨ mechanism, much anterior to Anderson’s are no “electric fields” appearing in the Lagrangian, and it realization, that works. Usually, it is the photon that “eats up” is expressed entirely in terms of the vortex configuration via the phase of the order parameter. The dual Stuckelberg¨ for- the Chern-Simons Gauss law constraint [38]. This constraint, 0 mulation, in which the scalar is encoded in the antisymmetric b=−Φ , implies that the momentum πξ conjugate to the phase tensor field, presents an opposite situation. To conclude here, ξ is not an electric charge, as it would be in superconductors, the electromagnetic response of the oblique superinsulator ex- but the vortex number itself. However, these boundary vor- presses fully in terms of the antisymmetric tensor field with tices carry also unit electric charge (2e), as can be seen by µ three degrees of freedom and mass given by Eq. (8). including the electromagnetic coupling (θ/2π)AµΦ . They are thus dyons themselves. Finally, it follows from the formula- The mass (8) is the sum of normal and topological contri- tion (10) of the boundary Lagrangian, that these dyons acquire butions. In the limit f 1, the topological contribution domi- fractional statistics via the boundary Chern-Simons term [39]. nates and the mass eventually diverges as f → ∞, Λ → ∞. In In particular, for the relevant case θ = 2π they become this limit, the bulk dynamics is frozen because only the topo- fermions. Thus, the topological limit of the oblique superin- logical contribution, sulators for θ = 2π contains only boundary states which are θ     massless fermionic dyons carrying unit magnetic and electric L = B + F µναβ B + F , (9) top 32π2 µν µν αβ µν charges. The canonical structure and the massless character of these boundary fermions rests on the U(1) combined gauge survives in the Lagrangian. The topological regime is realized symmetry λµ → λµ + ∂µχ and ξ → ξ − χ inherited from the deep in the underdoped region, where f 1. bulk. Any boundary perturbation that leaves this symmetry intact does not affect the dynamics of the boundary fermions, Let us now consider the oblique superinsulator on an open which are symmetry-protected surface states. Boundary per- manifold M with the boundary ∂M. The fields Bµν and Aµ turbations that break time-reversal symmetry can affect their in Eq. (7) contain transverse and longitudinal modes, the lat- statistics by effectively changing the value of θ. As a conse- ter being encoded in the gauge field λµ defined by Bµν = quence, these surface fermions behave as a Fermi liquid. ∂µλν − ∂νλµ and the scalar ξ defined by Aµ = ∂µξ. In the topological limit the mass (8) diverges and all bulk transverse modes get frozen. The longitudinal modes, however, con- DISCUSSION AND CONCLUSION tribute a total derivative into Eq. (9). When the model is de- fined on a bounded space, this total derivative for the longitu- The developed theory provides a consistent, unified in- dinal modes is all that remains and it gives rise to a boundary sight into the observed universal properties of the PG state and theory governed by the Lagrangian adjacent phases emerging near the tricritical point. First, the 2 symmetry-protected surface fermions in the dyon condensate θ µαν θ µ M 2 2 θ 2 2 L∂M = λµ ∂αλν + λµΦ − v (∂iξ) − (10)b form a Fermi liquid hence should exhibit resistivity ρ ∝ T , 8π2 2π 2 2π2 M which is exactly what is seen in the experiment [1]. As de- θ θ θ M θ2 λ µαν∂ λ λ 0 − bξ˙ − v2 (∂ ξ)2 − b2 , rived above, dyons are intimately associated with the pres- = 2 µ α ν + 0Φ i 2 (11) 8π 2π π 2 2π M ence of the θ-term in the electromagnetic response, presented µ µαβ where Φ = (1/2π) ∂α∂βξ represents the vortex current by Eq. (1), which naturally explains the onset of magnetoelec- i j and b = (1/2π) ∂iλ j plays the role of the “charge” canoni- tric effects in the pseudogap state, exactly as in strong topo- cally conjugate to the phase ξ. The kinetic terms for the sur- logical insulators, see [40]. The observation that supercon- face modes are added-in a posteriori and, correspondingly M ducting fluctuations and quantum correction contributions to is a non-universal mass scale and v is the non-universal prop- the resistive behavior near the dome are perfectly described agation speed of the surface modes. by the standard 2D formulas, see for example, [41, 42], indi- cates that the parameters of the boundary fermions are close The Hamiltonian of the boundary theory is derived by set- to those of the standard normal quasiparticle excitations in the 2 ting the canonical momenta πξ=−(θ/π)b and πλ2 =−(θ/4π )λ1 corresponding materials. The fact that the structurally sim- (as usual in pure Chern-Simons theories, one has a choice of ple model compound Hg1201 exhibits negligible residual re- deciding which of two components assumes the role of the sistivity [1] complies perfectly well with the notion that for 7 the symmetry-protected 2D boundary fermions localization These two can be considered as thermal and quantum contri- is absent. Importantly, while the cuprates exhibit quasi-2D butions, respectively. physics, the dyon condensate responsible for HTS proper- ties and the phase diagram structure is essentially of three- The resolution of the linear-T puzzle lies in the fact that in dimensional origin. This seeming contradiction is immedi- the presence of magnetic monopoles, the out-of-condensate ately resolved by the fact that the two adjacent CuO planes Cooper pairs behave as fermions [33, 34]. At low tempera- that present the unit cell of cuprates are enough to ensure the tures, these fermions scatter primarily by bogolons, the coher- existence of Dirac monopoles and the 3D-like interaction be- ent fluctuations of the residual dyon condensate domains [49]. tween monopoles and Cooper pairs. This explains why even At temperatures higher than the Bloch-Gruneisen¨ tempera- a few- and monolayer cuprates may exhibit practically the ture TBG ' 2skF/kB, where kF is the Fermi momentum of the same superconducting transition temperature Tc as bulk sam- fermionic Cooper pairs and s is the sound velocity of the ples [18, 42]. bogolons, the resistance due to this scattering mechanism is

linear in T, R∝T [49]. In HTS the zero-temperature super- Turning to the detailed Hall resistance experimental fluid density ns, and, accordingly, the effective ‘Fermi en- ∗ ∗ data [43], we recall that in optimally doped cuprates the ergy,’ F ' nsξ/4m , associated with fermionic Cooper pairs sign reversal of the Hall resistance is attributed to the con- are relatively small [50] (here ξ is the superconducting coher- tribution of vortices carrying in the core excessive electric ence length, and m∗ is the effective mass of the carriers). Ac- charge [42, 44]. It is natural to conjecture that, in the un- cordingly, the Bloch-Gruneisen¨ temperature TBG is low, espe- derdoped regime in the vicinity of the quantum critical point cially if the bogolon sound velocity s is also small. Thus, marking the onset of the superconducting dome, dyons that the universal scattering time τ ∼ ~/kBT (in physical units) carry both electric and magnetic charge take up on the role of for T > TBG arises similarly to the usual high temperature vortices with the excessive electric charge. Moreover, recall- electron-phonon scattering time. To understand the origin of ing further that the data on the charge Berezinskii-Kosterlitz- the different slopes at low and high temperatures, note that Thouless transition into the superinsulating state [45] indicate the linear dependence can be understood in the framework of that the boundary fermionic states are harbored by the per- the Anderson orthogonality catastrophe (AOC) approach. In- colation Chalker-Coddington structure, which is known to be deed, one expects that all the relevant microscopic times are a natural host for gapless boundary states [46], one naturally shorter than the scattering time τ. According to the AOC, τ is expects erratic behavior with noticeable random components, related to the heat, Q, generated by the scattering processes exactly as observed in the experiment [43]. However, further by the universal relation [51] ~/τ = f Q, where f 'O(1) theoretical study of the dyon-related Hall effect for a quanti- is some numerical function of the microscopic parameters. tative description of the experiment is needed. Adopting Q = kBTδS , where δS is the entropy change as- sociated with the scattering and weakly (logarithmically) de- An underdoped antiferromagnetic phase of cuprates hold- pending on temperature [51], one arrives at the universal rela- ing exactly a single hole per Cu site, realizes a 2D spin square tion 1/τ = αkBT/~, α = O(1), giving rise to the linear-in-T Heisenberg antiferromagnet. Thus, the symmetry of the lattice resistivity. Since the microscopic origin of scattering at high is C4. In such a lattice the adjacent spins are opposite while and low temperatures is different, one expects the different the diagonal ones are parallel. As a result, the combined elec- scattering-related entropy changes, hence different slopes. tric/magnetic symmetry of the lattice reduces to C2. Since due to the θ-term, the boundary fermions are dyons and carry both At the high magnetic field limit, the number of the available electric and magnetic charge, they feel the combined C2 sym- degenerate ground states occupied by bosons in the area A per- metry, hence nematicity with the nematic director going along pendicular to the magnetic field B is N = (1/2π)(eB/~) (see, the CuO lattice diagonals. This essential characteristic of the e.g. [52]). This is the number of available scattering centres PG state [3,8] goes hand-in-hand with the T 2 resistance. at T = 0 and the same reasoning that leads to the linear-in-T resistance leads to the observed linear-in-B behaviour of the

Now we discuss the linear, R∝T, resistance behavior in the intercept. In other words, at high fields the role of the temper- ∗ overdoped region, observed at high magnetic fields destroying ature T is taken by ~ωc, where ωc = eB/m c is the cyclotron the Cooper pair condensate at low temperatures and declared frequency [53]. a major puzzle of condensed-matter physics [47]. While its generic character in HTS is established experimentally and is Finally, as we have mentioned above, the boundary associated with a universal scattering rate, its origin remains a fermionic states are harbored by the percolation Chalker- mystery. A careful experimental analysis [48] reveals that the Coddington structure in accord with experimental findings slopes of the linear terms in the low-temperature regime and of [43]. It implies that one can identify T ∗ with the temper- in the high-temperature, inherently metallic phase, are notice- ature at which the ”bubbles” hosting massless boundary dyon able different. Moreover, at sufficiently low temperatures and states form finite clusters but do not extend across the whole high magnetic field B there is an intercept which is also linear system. Then the T ∗∗ line in the phase diagram is the line in B [48], so that the resistivity assumes the form ρ = aT +bB. of the percolation transition at which the Chalker-Coddington 8 structure spreads across the system and the pseudogap state This can be achieved only by introducing integers k ∈ [−b, b], with nematicity, the T 2 resistance and magnetoelectric effects called Matsubara frequencies, and restrict the allowed values is fully formed. to k0 = πk/b. Correspondingly, integrals over k0 have to be replaced by sums over Matsubara frequencies,

k=b ! APPENDIX Z π   X π πk dk0 f k0 → f . (19) β b −π k=−b Lattice BF term As a consequence, at finite temperatures G(0) becomes To formulate the gauge-invariant lattice BF-term, we fol- 1 kX=+b π Z π dk1dk2dk3 low [12] and introduce the lattice BF operators G(0, T) = . (2π)4 b  πk 2 P3  ki 2 k=−b −π 4 sin + 4 sin K ≡ S  , Kˆ ≡  ˆ Sˆ , 2b i=1 2 µνρ µ µανρ∆α µνρ µναρ∆α ρ (13) (20) where This affects primarily the parameter η which becomes a function of the temperature,√ η(T) = πG(0, T)/µ. Introducing ∆µ f (x) ≡ f (x + `µˆ) − f (x) , S µ f (x) ≡ f (x + `µˆ) , the scale function S (T) = G(0)/G(0, T) we obtain the result ˆ ˆ that, at finite temperatures the overall scale of the ellipse (6) ∆µ f (x) ≡ f (x) − f (x − `µˆ) , S µ f (x) ≡ f (x − `µˆ) , (14) p p shrinks by S (T), 1/η → S (T) 1/η. are the forward and backward lattice derivative and shift op- erators, respectively. Summation by parts on the lattice inter- changes both the two derivatives (with a minus sign) and the two shift operators; gauge transformations are defined using Topological two-form mass term in 3D the forward lattice derivative. The two lattice BF operators are interchanged (no minus sign) upon summation by parts on Let us consider the following model for an antisymmetric the lattice and are gauge invariant in the sense that tensor bµν in 3D [36] K ∆ = K ∆ = ∆ˆ K = 0 , 1 1 θ µνρ ν µνρ ρ µ µνρ L h hµνα − b bµν b µναβb , = 2 µνα 2 µν + 2 µν αβ (21) Kˆ µνρ∆ρ = ∆ˆ µKˆ µνρ = ∆ˆ νKˆ µνρ = 0 . (15) 12Λ 4 f 32π

They also satisfy the equations where hµνα = ∂µbνα + ∂νbαµ + ∂αbµν is the three-form field strength, Λ has dimension [mass] and f and θ are dimension-   2 Kˆ µνρKρλω = − δµλδνω − δµωδνλ ∇ less. The first two terms are the generalization to two forms     of the Proca Lagrangian for a massive vector field. The third + δ ∆ ∆ˆ − δ ∆ ∆ˆ + δ ∆ ∆ˆ − δ ∆ ∆ˆ , (16) µλ ν ω νλ µ ω νω µ λ µω ν λ term is topological, since it is metric-independent.  2  Kˆ µνρKρνω = KµνρKˆ ρνω = 2 δµω∇ − ∆µ∆ˆ ω , (17) The equations of motions of this model are

2 ˆ Λ2 Λ2θ where ∇ = ∆µ∆µ is the lattice Laplacian. ∂ hµαβ + bαβ − αβγδb = 0 . (22) µ f 2 8π2 γδ

Contracting with ∂α we obtain the condition Finite Temperature Deconfinement Transition f 2θ ∂ bµν + hν = 0 , (23) µ 4π2 In field theory, a finite temperature T is introduced by µ µναβ formulating the action on a Euclidean time of finite length where h = (1/6) hναβ is the dual field strength. Finally, γ β = 1/T, with periodic boundary conditions (we have reab- contracting (22) with νγαβ∂ and using the above condition sorbed the Boltzmann constant into the temperature). If the we get original field theory model is defined on a Euclidean lattice of s !2  2 spacing `, then β is quantized in integer multiples of 1/` as  2 2 µ f Λ 4π θ ∂ + m h = 0 , mθ = + . (24) β = b/`, with b an integer. θ 4π f 2 π The lattice Coulomb Green’s function G(0) at coinciding This shows that there is a topological contribution to the mass. points is 1 Z π 1 Acknowledgements G(0) = d4k . (18) 4  i 2 We are delighted to thank Prof. Ivan Bozoviˇ c´ for illumi- (2π) −π P3 k i=0 4 sin 2 nating discussions. The work at the University of Twente (V.M.V.) was partly supported by the NWO grant. M.C.D. 0 At finite temperatures, k is both defined on a Brillouin zone thanks CERN, where she completed this work, for kind hos- 0 0 of length 2π and also invariant under shifts k → k + 2π/`β. pitality. 9

DATA AVAILABILITY (2013). [20] V.M. Vinokur, et al. Superinsulator and quantum synchroniza- tion. Nature 452, 613 – 615 (2008). Data sharing not applicable to this article as no datasets [21] A. Polyakov, Compact gauge fields and the infrared catastro- were generated or analyzed during the current study. phe. Physics Letters B 59, 82-84 (1975). [22] P.Goddard and D. I. Olive, Magnetic monopoles in gauge fields theories, Rep. Progr. Phys. 41, 1357-1437 (1978). [23] C. A. Trugenberger, M. C. Diamantini, N. Poccia, F. S. Nogeuira, and V.M. Vinokur, Magnetic monopoles and superinsulation [1] N. Barisi˘ c,´ et al. Universal sheet resistance and revised phase in Josephson junction arrays. Quantum Reports 2, 388 – 399 diagram of the cuprate high-temperature superconductors. (2020). PNAS 110, 12235 – 12240 (2013). [24] M. C. Diamantini, C. A. Trugenberger, and V.M. Vinokur, [2] B. Keimer, S. A. Kivelson, M. R. Norman, S. Uchida, and Quantum magnetic monopole condensate. In press. J. Zaanen, From quantum matter to high-temperature supercon- [25] M. C. Diamantini, C. A. Trugenberger, and V.M. Vinokur, Con- ductivity in copper oxides. Nature 518, 179–186 (2015). finement and asymptotic freedom with Cooper pairs. Comm. [3] Y.Sato, et al. Thermodynamic evidence for a nematic phase Phys. 1, 77 (2018). transition at the onset of the pseudogap in YBa2Cu3Oy. Nature [26] M. C. Diamantini et al, Direct probe of the interior of an electric Physics 13, 1074 – 1078 (2017). pion in a Cooper pair superinsulator. Communications Physics [4] P.Zhou, et al. Electron pairing in the pseudogap state revealed 3, 142 (2020). by shot noise in copper oxide junctions. Nature 572, 493 (2019). [27] F. Wilczek, Two applications of axion electrodynamics. Phys. [5] C. Proust and L. Taillefer, The Remarkable Underlying Ground Rev. Lett. 58, 1799 – 1803 (1987). States of Cuprate Superconductors. Annu. Rev. Condens. Mat- [28] J. L. Cardy and E. Rabinovici, Phase structure of Z(p) models in ter Phys. 10, 409 – 429 (2019). presence of theta parameter. Nucl. Phys. B 205, 1-16 (1982). [6] C. Proust, B. Vignolle, J. Levallois, S. Adachi, and [29] E. Witten, Dyons of Charge 2θ/2π. Phys. Lett. 86, 283-287 N. E. Hussey, Fermi liquid behavior of the in-plane resis- (1979). tivity in the pseudogap state of YBa2Cu4O8. PNAS 113, 13564 [30] T. Banks, R. Myerson, and J. Kogut, Phase transitions in – 13659 (2016). Abelian lattice gauge theories. Nucl. Phys. B 129, 493-510 [7] S. Mukhopadhyay, et al. Evidence for a vestigial nematic state (1977). in the cuprate pseudogap phase. PNAS 116, 13249 – 13254 [31] G. ’t Hooft, Topology of the gauge condition and new confine- (2019). ment phases in non-Abelian gauge theories. Nucl. Phys. B 190, [8] S. A. Kivelson, E. Fradkin, and V.J. Emery, Electronic liquid- 455-478 (1981). crystal phases of a doped Mott insulator. Nature, 393, 550–553 [32] M. A. Metlitski, C. L. Kane, and M. P.A. Fisher, Bosonic topo- (1998). logical insulator in three dimensions and the statistical Witten [9] L. Nie, G. Tarjus, and S. A. Kivelson, Quenched disorder and effect. Phys. Rev. B 88, 035131 (2013). vestigial nematicity in the pseudogap regime of the cuprates. [33] R. Jackiw and C. Rebbi, Spin from isosopin in a gauge theory. PNAS, 111, 7980 – 7985 (2014). Phys. Rev. Lett. 36, 1116 – 1119 (1976). [10] G. Q. Zheng, P.L. Kuhns, A. P.Reyes, B. Liang, and C. T. Lin, [34] P.Hasendfratz and G. ’t Hooft, A fermion-boson puzzle in a Critical point and the nature of the pseudogap of single-layered gauge theory. Phys. Rev. Lett. 36, 1119 – 1122 (1976). copper-oxide Bi2Sr2−xLaxCuO6+δ superconductors. Phys. Rev. [35] A. S. Goldhaber, R. MacKenzie, and F. Wilczek,Field correc- Lett. 94, 047006 (2005). tions to induced statistics. Mod. Phys. Lett. 4, 21-31 (1989). [11] A. M. Goldman, Superconductor-insulator transitions. Int. J. [36] F. Quevedo and C. A. Trugenberger, Phases of antisymmetric Mod. Phys. B 24 4081– 4101 (2010). tensor field theories. Nucl. Phys. B 501, 143 – 172 (1997). [12] M. C. Diamantini, P.Sodano, and C. A. Trugenberger, Gauge [37] M. C. Diamantini, F. Quevedo, and C. A. Trugenberger, Confin- theories of Josephson junction arrays. Nuclear Physics B 474, ing Strings with Topological Term. Phys. Lett. B 396, 115 – 121 641 – 677 (1996). (1997). [13] D. Das and S. Doniach, Existence of a Bose metal at T = 0. [38] R. Jackiw, Topics in planar physics. in Physics, Geometry and Phys. Rev. B 60, 1261 – 1275 (1999). Topology, H. C. Lee ed. Plenum Press, New York (1989). [14] A. Kapitulnik, S. A. Kivelson, and B. Spivak, Anomalous met- [39] F. Wilczek, Fractional statistics and anyon superconductivity. als: failed superconductors. Rev. Mod. Phys. 91, 011002 World Scientific, Singapore (1990). (2019). [40] A. M. Essin, J. E. Moore, and D. Vanderbilt, Magnetoelectric [15] M. C. Diamantini, et al. Bose topological insulator intermedi- Polarizability and Axion Electrodynamics in Crystalline Insu- ate state in the superconductor-insulator transition. Phys. Lett. lators. Phys. Rev. Lett. 102, 146805 (2009). A 384 126570 (2020). [41] A. Pomar, M. V.Ramallo, J. Mosqueira, C. Torron,´ and F. Vidal, [16] Y.-M.Lu and A. Vishwanath, Theory and Classification of in- Fluctuation-induced in-plane conductivity, magnetoconductiv- teracting integer topological phases in two dimensions: a ity, and diamagnetism of Bi2Sr2CaCu2O8 single crystals in Chern-Simons approach. Phys. Rev B 86, 125119 (2012). weak magnetic fields. Phys. Rev. B 54, 7470 – 7480 (1996). [17] X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, Symmetry pro- [42] F. Zhao, et al. Sign-Reversing Hall Effect in Atomically tected topological orders and the group cohomology of their Thin High-Temperature Bi2.1Sr1.9CaCu2.0O8+δ Superconduc- symmetry group. Phys. Rev. B 87, 155114 (2013). tors. Phys. Rev. Lett. 122, 247001 (2019). [18] Y.Yu, et al. High-temperature superconductivity in monolayer [43] J. Wu, A. T. Bollinger, Y.Sun and I. Bozovic, Hall effect in Bi2Sr2CaCu2O8+δ. Nature 575, 156 – 163 (2019). quantum critical charge-cluster glass, PNAS 113 4284 – 4289 [19] A. Vishwanath and T. Senthil, Physics of Three-Dimensional (2016). Bosonic Topological Insulators: Surface-Deconfined Criticality [44] M. V.Feigelman, V.B. Geshkenbein, A. I. Larkin, and and Quantized Magnetoelectric Effect. Phys. Rev. X 3 011016 V.M. Vinokur, Sign change of the flux-flow Hall effect in 10

HTSC. JETP Lett. 62, 834 – 837 (1995). Bogolon-mediated electron scattering in graphene and hybrid [45] A. Yu.Mironov, et al. Charge Berezinskii-Kosterlitz-Thouless Bose Fermi systems. Phys. Rev. B 99 115408 (2019). transition in superconducting NbTiN films. Scientific Reports [50] V.J. Emery and S. A. Kivelson, Importance of phase fluctua- 8 4082 (2018). tions in superconductors with small superfluid density. Nature [46] T. Ochiai, Gapless surface states in a three- 374, 434 – 437 (1995). dimensional Chalker-Coddington type network model. [51] A. V.Lebedev and V.M. Vinokur, Heat generation due to An- http://arxiv.org/abs/1510.04033v1 (2015). derson catastrophe in mesoscopic devices. Phys. Rev. B, in [47] A. Legros, et al. Universal T-linear resistivity and Planckian press. dissipation in overdoped cuprates. Nature Physics 15, 142 – [52] A. A. Abrikosov, Fundamentals of the Theory of Metals. North 147 (2019). Holland, Amsterdam, 1988. [48] P.Giraldo-Gallo et. al. Scale-invariant magnetoresitance in a [53] A. Glatz, A. A. Varlamov, and V.M. Vinokur, Fluctuation spec- cuprate superconductor. Science 361, 479-481 (2018). troscopy of disordered two-dimensional superconductors. Phys. [49] M. Sun, K. H. S. Villegas, V.M. Kovalev, and G. Savenko, Rev. B 84, 104510 (2011).