Charge-4e superconductivity from nematic superconductors in 2D and 3D Shao-Kai Jian,1 Yingyi Huang,2 and Hong Yao2, 3, ∗ 1Condensed Matter Theory Center, Department of Physics, University of Maryland, College Park, Maryland 20742, USA 2Institute for Advanced Study, Tsinghua University, Beijing 100084, China 3State Key Laboratory of Low Dimensional Quantum Physics, Tsinghua University, Beijing 100084, China (Dated: February 8, 2021) Charge-4e superconductivity as a novel phase of matter remains elusive so far. Here we show that charge-4e phase can arise as a vestigial order above the nematic superconducting transition temperature in time-reversal-invariant nematic superconductors. On the one hand, the nontrivial topological defect|nematic vortex|is energetically favored over the superconducting phase vortex when the nematic stiffness is less than the superfluid stiffness; consequently the charge-4e phase emerges by proliferation of nematic vortices upon increasing temperatures. On the other hand, the Ginzburg-Landau theory of the nematic superconductors has two distinct decoupling channels to either charge-4e orders or nematic orders; by analyzing the competition between the effective mass of the charge-4e order and the cubic potential of the nematic order, we find a sizable regime where the charge-4e order is favored. These two analysis consistently show that nematic superconductors can provide a promising route to realize charge-4e phases, which may apply to candidate nematic superconductors such as PbTaSe2 and twisted bilayer graphene. Introduction.|Featuring the condensation of quartets 4e superconductivity in these high-temperature super- with four times the fundamental electron charges, charge- conductors has not been observed so far. Apart from 4e superconductivity [1{4] is intrinsically distinct from PDW superconductors, the nematic superconductor that the conventional charge-2e superconductivity discovered breaks both charge conservation and lattice rotational more than a century ago. One hallmark of charge-4e symmetry also hosts an order parameter that carries superconductors is the magnetic flux quantization with multiple quantum numbers. The additional topolog- period hc=4e, which is half that of the usual super- ical defect in nematic superconductors is created by conducting flux quantum. Unlike the conventional su- the intersection of domain walls separating different perconductivity that is well described by the seminal ground states. More importantly, many experimental Bardeen{Cooper{Schrieffer (BCS) theory, many proper- progresses have been made in achieving the nematic su- ties of charge-4e superconductivity remain less well un- perconductivity in various systems. For instance, re- derstood. The mystery of this phase is not only be- cent experiments found that the spin susceptibility be- cause it defies any BCS-type analysis, but also due to the low the superconducting temperature breaks the three- lack of experimental realization so far. Previous studies fold lattice rotational symmetry in doped topological in- in Ref. [3] suggested that the charge-4e superconductiv- sulator CuxBi2Se3 [26{29], where the breaking of rota- ity can arise as a vestigial long-range order above the tional symmetry suggests that the order parameter in transition temperature of a pair-density-wave (PDW) su- the superconducting phase is a two-dimensional Eu rep- perconductor [5{17], whose order parameter varies peri- resentation of point group D3d [30{32]. Providing the odically in the real space. The additional breaking of growing experimental evidence of the nematic super- translational symmetry in the PDW state is essential in conductivity in various systems, including doped topo- realizing the charge-4e superconductivity as it provides logical insulators CuxBi2Se3 [26{29], SrxBi2Se3 [33, 34] nontrivial topological defects|dislocations. As a conse- and NbxBi2Se3 [35, 36], superconducting topological quence, if the dislocation proliferates as the temperature semimetals PbTaSe2 [37], and more recently twisted bi- is raised, the equilibrium state will restore the transla- layer graphene [38, 39], as well as increasing interest in tional symmetry by doubling the original charge-2e con- novel charge-4e ordering, it is thus interesting to ask arXiv:2102.02820v1 [cond-mat.str-el] 4 Feb 2021 densate, and therefore lead to the charge-4e phase. whether the multiple-component superconducting order parameters in nematic superconductors is able to realize The crucial ingredient of the underlying condensate charge-4e phases above their charge-2e transition tem- having more than one quantum number [18{20] provides perature. a general guide to search for charge-4e superconductors. It was proposed that in certain two-dimensional high- In this paper, we answer this question in the affirma- temperature superconductors (e.g. La2−xBaxCuO4 and tive; namely, we provide two analysis, suitable for 2D and La1:6−xNd0:4SrxCuO4 [21{23]) exhibiting stripe super- 3D, respectively, to support the possibility of a vestigial conducting orders [24, 25], the charge-4e superconduc- charge-4e superconducting phase from nematic supercon- tivity may occur above the transition temperature [3]. ductors. In 2D, because a vortex is a point-like object, we Nevertheless, the experimental signature of the charge- analyze the fate of various topological defects to deter- 2 mine the phase diagram. After identifying three distinct topological defects that are responsible for three differ- ent orders out of the nematic superconducting phase, we find that the competition between the superfluid stiffness and the elastic constant leads to a rich phase diagram, as shown in Fig.1. In particular, if the elastic constant of the material is less than one third of the superfluid stiff- ness, the energetically favored nematic vortices are pro- liferated, resulting in the novel charge-4e phase, when the temperature gets raised above the transition tem- FIG. 1. Schematic phase diagram of the defect theory. κ, perature. In 3D, we employ a Ginzburg-Landau theory ρ and T denote the nematic stiffness, the superfluid stiffness near the transition point where the nematic supercon- and temperature, respectively. The solid lines refer to phase boundaries. Blow the dashed line, the three-fold anisotropy is ducting order parameter is well described by a simple relevant from the lowest order calculation. Nevertheless, the field theory up to quartic terms. We are able to show true transition to the charge-4e phase belongs to the three- analytically that the effective mass of the charge-4e or- state Potts model universality class [40] where the three-fold der is less than that of the nematic order by treating the anisotropy is relevant at the transition point. three-fold anisotropy perturbatively, which indicates the charge-4e order is favored. We further confirm our re- sults by numerically solving the saddle-point equation. The single-valueness of nematic superconducting or- Our result in 3D suggest that materials whose dispersion der parameter ∆± uniquely determines that the topolog- along the third direction is weaker than the in-plane dis- ical defects are given by (δθ; δφ) = (2π; 0); (0; 2π); (π; π), persion such as in SrxBi2Se3 are promising in realizing where δθ, δφ denote the winding of the phase around a the charge-4e phase. defect. Physically, they correspond to the superconduct- Charge-4e phase from proliferating topological defect.| ing vortex, the nematic double vortex and the supercon- We start by analyzing the symmetries in the nematic su- ducting half-vortex binding with a single nematic vortex, perconducting phase providing an appropriate language respectively. In the nematic superconducting phase, the for clarifying the topological defects, and then turn to effective Hamiltonian characterizing the phase modes is renormalization group analysis of the defect theory. As ρ κ H = (@θ)2 + (@φ)2 − g cos 6πφ, (3) mentioned above, we consider the order parameter of the 2 2 6 nematic superconducting phase being a two-component T where ρ is the superfluid stiffness (superfluid density) complex boson ∆ = (∆x; ∆y) carrying Eu represen- and κ is the nematic stiffness (elastic constant). No- tation in D3d group, and each complex component ad- ditionally hosts the U(1) quantum number correspond- tice that in writing this Hamiltonian, we have changed ing to the charge conservation, i.e., it is the conven- the compactification of the phase mode from 2π to 1 for tional charge-2e condensate. In the phase respecting convenience. Note that cos 6πφ is allowed owing to the time-reversal symmetry, the relative phase between two three-fold anisotropy, i.e., the action is invariant under components ∆ ; ∆ is pined at 0 or π. As we focus on φ!φ+1=3. Thermal proliferation of topological defects x y ~ ~ the time-reversal-invariant nematic superconductor, we can be then described by the dual bosons φ and θ [3], hereafter assume the relative phase between ∆x and ∆y T T is pined at 0 or π in the analysis of low-energy physics. H = (@θ~)2 + (@φ~)2 − g cos 6πφ − g cos 2πθ~ 2ρ 2κ 6 2;0 In terms of the basis ∆ = ∆ ±i∆ , one can bring the ± x y ~ ~ ~ the four-dimensional field configuration into three phase −g0;2 cos 2πφ − g1;1 cos πθ cos πφ, (4) modes and one amplitude mode. Moreover, the time re- where g , g and g are couplings characterizing the versal invariance makes one of the phase modes heavier, 2;0 0;2 1;1 strength of creating (annihilating) each kind of topolog- thus we have two phase modes and one amplitude mode ical defects, respectively, and T is the temperature. in the low-energy sector, The standard renormalization group flow for Eq. (4) i(θ+φ) i(θ−φ) ∆+ = j∆je ; ∆− = j∆je : (1) to the lowest-order reads where θ and φ denote the two phase modes. The φ field dg πρ 2;0 = 2 − g ; (5) describes the U(1) rotation between amplitudes of two dl T 2;0 components (i.e. the spatial rotation), whereas the θ field dg πκ 0;2 = 2 − g ; (6) is the U(1) phase conjugated to the global charge.