arXiv:2004.09767v4 [cond-mat.supr-con] 7 Jun 2020 odtosbtentecretadvlaeo h elec- the of voltage and current the between conditions unu hs lpjnto ( junction slip phase a quantum and for a junction other Josephson each nanowire-based to one-dimensional dual paper Hamiltonians previous two con- a introduced In we quark freedom. elucidate asymptotic and and finement chromodynamics QCD quantum for of environment desktop powerful phenomena attracted a have superinsulat- as superinsulators attention be Recently, to materials. observed ing experimentally been have inta sa h uljnto oJspsnjunction. Josephson TiN to studies, junction dual previous the as In is junc- superinsulator that / tion superconductor / superinsulator a uao ntecs fcnesdmte a rediscovered was al. matter et condensed superin- Diamantini of by The Hooft of case ’t confinement. the concept in quark by The sulator explaining 1978 theory in wire. a conceived as in first Achieved was InOx superinsulator by with performed loop was ing slip al. phase et quantum Astafiev coherent demonstra- of experimental tion an sub- and al. was wires, et slip superconducting phase Mooij the- quantum by A coherent (SC). mitted of superconductors study known to oretical phenomenon mainly dual is a which as superinsulator, phe- a called is other nomenon The system. superconducting dimensional n21,Ynd tal. Furthermore, et Yoneda environment. mass 2012, dissipative in finite a a in had traveled vortices and which in model nomenological nasl-ulapoiain hn n19,Kr¨amer and 1998, in Then, Doniach arrays approximation. junction self-dual a charge-vortex Josephson in planar describing to theory equivalent coupling gauge Abelian Simons ulpeoeo faJspsnjnto ( junction Josephson a a as of known mainly phenomenon is called dual which phenomenon slip, a have phase is quantum superconductivity, them coherent of in One states attention. dual attracted be to considered ult ewe oeetqatmpaesi n oeho j Josephson and slip phase quantum coherent between Duality nrcn eae,totpso hnmn,wihare which phenomena, of types two decades, recent In 7 lordsoee ueisltr sn phe- a using superinsulators rediscovered also 12,13 ASnmes 74.20.z B numbers: dual PACS the and superinsulators shown. also of is phenomenology theory, the is which erlal ihnteVlanapoiain h possibil The approximation. zer Villain at the dimension equ 1 within + is reliable 2 This be the in method. superinsulator Hamiltonian and conductivity dual the using investigated 4 h ult ewe oeetqatmpaesi n Josephs and slip phase quantum coherent between duality The hc meddi ag superconduct- large a in embedded which .INTRODUCTION 1. nodrt elz igeclrversion single-color a realize to order in 1 ih nvriyo ehooy 02UaoiNishihasama- Umanori 50-2 Technology, of University Aichi 2 colo niern,TkoDniUiest,5SnuAsah Senju 5 University, Denki Tokyo Engineering, of School 6 n19 saproi ie Chern- mixed periodic a as 1996 in 8 1 eicvrdsprnuaosfrom superinsulators rediscovered - 3 9 sflxqat hog nano through quanta flux As InO , eemndb h ulHmloinmethod Hamiltonian dual the by determined QPSJ 10 .Yoneda, M. n NbTiN and ,adapiddual applied and ), JJ naone- a in ) 1 Ωae:Jn ,2020) 9, June (ΩDated: 11 .Niwa, M. 2-5,Japan 120-855, films 14 5 , 2 anapoiaini h 1 + 2 the in approximation lain eueteda aitna ehdfrbtenthe between 2, for Section the method In and Hamiltonian follows: model dual as JJ the organized is use paper we this of rest a odce.Orrslsso hti h 1+1 the in that show results Hamiltonians Our two these conducted. was from temperature zero at system e oteter f2+1 + 2 of theory the to per anproeo hswr a oetn h hoyof theory the extend to was 1 work + this 1 of The purpose phenomena. main equivalent completely su- and are slip perinsulator phase quantum coherent temperature, zero at agdda anisotropic dual gauged h ili prxmto na2+1 + 2 a in approximation Villain the rvosscin eso htteei xc ult be- duality exact is the there between the in that tween and derived show we relationship space, section, amplitude lattice previous the compact and phase a of the variables in conjugate order canonical and infinite nonlinear with Lagrangians the Hamiltonians the on the based between the 3, transformation between Section Legendre relations In am- dual the a constants. and various and phase relations, the dual between plitude a derive to temperature ,teda izugLna ( Ginzburg-Landau dual the 6, otayt eto ,w eieada transformation dual a derive we the 4, from Section to contrary ae,dsuswehrtemnmmui fcag con- this charge of of conclusions unit minimum the the summarize discus- whether we discuss and 8, paper, summary Section the in In sion approach. energy effective ut ftepeiu eto ydrvn ultrans- dual a deriving re- anisotropic by the the of section validity from previous the formation the prove we of 4, sults Section In discussed. the between diagram phase the n tt)adthe and state) ing ueisltrtasto na1+1dmninl( 1 + (1 dimensional 1 + 1 a sys- in dual superconductor- transition the a superinsulator of proposed. discussion construct was general method, to a Hamiltonian Furthermore, theory dual the general called A tem, circuit. tric ntesprnuao a hw.I eto ,w cal- we the 7, of Section value In critical shown. the flux was electric culated of superinsulator confinement the gauged of in the possibility of the approximation and model, field mean the from n .Motohashi M. and d t fteda izugLna theory, Ginzburg-Landau dual the of ity ytm nnnwrssoni h rvospa- previous the in shown nanowires on systems eprtr.Ti ehdpoe to proved method This temperature. o DAXY vln otedaiybtensuper- between duality the to ivalent JJ Ster,wihi microscopic a is which theory, CS oe n the and model h,Aci4304,Japan 443-0047, Aichi cho, njnto nnnset was nanosheets in junction on oe othe to model -h,Aah-uTokyo Adachi-ku i-cho, QPS QPSJ 2 tt sprnuaigsae was state) (superinsulating state nto nananosheet a in unction XY oe nannseta zero at nanosheet a on model d ytm nnnset.The nanosheets. on systems ( QPSJ DAXY AXY DGL XY JJ d QPS ytm nScin5, Section In system. oe ihgueby gauge with model oe.Additionally, model. hoywsderived was theory ) ( tt (superconduct- state d oe yteVil- the by model ) AXY ytm nSection In system. mltd ythe by amplitude oe othe to model ) d system QPSJ d ) 2
finement in a superinsulator is 2e or e, and, finally, we other hand, when the superinsulator region acts strongly, describe the possibility of dual BCS theory. In Appendix the QPSJ state becomes strong, and the Hamiltonian in A, the numerical evaluation of the anisotropic massless this case is as follows: lattice Green’s function is shown at the origin of x = 0, M 2 M which is necessary for loop correction. In Appendix B, 2 H = E N˜˜ (x) + E 1 cos θ˜(x) , (3) QPS L θ S − ∇j the effective energy approach for the QPSJ model is ex- x=1 j=1 x=1 plained. X h i X X h i
where, N˜˜ (x) N˜˜ (x, τ) means the magnetic flux quan- θ ≡ θ tum numbers θ˜(x) θ˜(x, τ) of the magnetic flux quan- 2. DUAL HAMILTONIAN METHOD BETWEEN tum field. E ≡Φ 2/2L and E 2eV /2π are THE JOSEPHSON JUNCTION (JJ) MODEL AND L ≡ 0 S ≡ c THE QUANTUM PHASE SLIP the inductive energy per magnetic flux quantum and the JUNCTION(QPSJ) MODEL ON A NANOSHEET QPS amplitude, respectively, where Vc is the critical volt- AT ZERO TEMPERATURE age. From the Hamiltonian of Eq.(3), the dual Josephson equations are as follows:
~ ˜ ˜ ∂θ (x) 2Nθ˜ (x) V˜ (x)= i = EL, Φ0 ∂τ Φ0 2π I˜ (x)= E sin θ˜(x) , (4) j 2e S ∇j
where, V˜ (x) and I˜j (x) (j = 1, 2) are the dual voltage and the j components of dual current in the QPSJ, re- spectively. As the first step of the dual Hamiltonian method, two dual conditions8,14 between Eq.(2) and the FIG. 1: A checkered nanosheet consisting of superconductors dual equations of Eq.(4) are assumed as follows: (black squares) and superinsulators (white squares). V (x) I˜(x) , I (x) V˜ (x) ,. (5) ≡ ≡ Consider a checkered nanosheet consisting of super- conductors and superinsulators as shown in Figure 1. In where, I √I 2 + I 2 and I˜ I˜2 + I˜2 are the intensity ≡ 1 2 ≡ 1 2 such a 2d checkered nanosheet, when the superconductor of the current in the JJ andq the intensity of the dual region acts strongly, the JJ state becomes strong, and current in the QPSJ, respectively. When the condition the Hamiltonian in this case is as follows: of Eq.(5) is imposed between Eqs.(2) and (4), the follow- M M 2 ing two relational expressions between the phase and the H = E [N (x)]2 + E [1 cos θ (x)], (1) number of particles between dual systems are derived as JJ c θ J − ∇j x=1 x=1 j=1 X X X shown below. One of them is the relationship between x the phase θ of the Cooper pair and the magnetic flux where, Nθ(x) Nθ( , τ) means the number of the par- ˜ ≡ quantum numbers Nθ˜, and the other is the relationship ticles in the Cooper pair that is the canonical conju- ˜ gate to the phase θ of the Cooper pair, and the space between the phase θ of the magnetic flux quantum and difference14,15 of the phase θ (x) θ (x, τ) is defined by the number Nθ of the Cooper pair, as follows: ≡ j θ(x, τ) θ(x, τ) θ (x ae, τ); τ ~β (β 1/kBT is the re- verse∇ temperature),≡ − x−= (x , x ≡) and M≡ l2/a2 (l and 1 2 1 2 N (x)= sin2 θ˜(x), ≡ θ 2π v ∇j aare the space length and lattice spacing, respectively)are uj=1 uX the imaginary time, the lattice coordinate points, and the t 2 total number of lattices in the two-dimensional lattice 1 2 2 N˜˜ (x)= sin j θ (x), (6) space, respectively; and E (2e) /2C and E Φ I /2π θ 2π v ∇ c J 0 c uj=1 ≡ ≡ uX are the charging energy per Cooper pair and the Joseph- t son energy, respectively, where C,Ic and Φ0 = h/2e are In Eq.(6), the leftmost equation shows that the Cooper the capacitance, the critical current and the magnetic pair number is proportional to the strength of the flux flux-quantum, respectively. From the Hamiltonian of Eq. quantum current, and the next equation Is shown that (2), Josephson’s equations are as follows: the flux quantum number is proportional to the strength of the Cooper pair current. Since the Cooper pair cur- ~ ∂θ (x) 2Nθ (x) V (t)= i = Ec, rent and the flux quantum current form a closed loop, 2e ∂τ 2e 2π the number of Cooper pairs and the flux quantum num- I (x)= E sin θ (x) . (2) j −Φ J ∇j ber can be considered as the winding number of the cur- 0 rent loop formed from each dual current. If the relation- where, V (x) and Ij (x) (j =1, 2) are the voltage and the j ships described in Eq.(6) are satisfied, the relationship components of the current in the JJ, respectively. On the between the QPS amplitude and the charging energy per 3 single charge, and the relationship between the Joseph- function of the JJ state (superconducting state) at zero son energy and the inductive energy per magnetic flux- temperature is as follows: quantum, are as follows: Z exp E′0 + E′ MM Z′ , AXY ≡ − J J τ AXY Ec EL n o2 ES = 2 , EJ = 2 . (7) 2π 2π Z′ Dθ exp E′0 cos θ (x)+ E′ cos θ (x) . (15) AXY ≡ J ∇τ J ∇j x,τ j=1 Furthermore, the inductance and capacitance are related Z X X to the critical current and the critical voltage, respec- Mτ M Mτ M π dθ(x,τ) tively, as follows: where and Dθ 2π are sums x ≡ τ=1x=1 ≡ τ=1x=1−π Φ 2e and pathP integrals,P P respectively,R Q inQ 2R + 1d lattice space L = 0 , C = , (8) 2πIc 2πVc x = (x, τ). Similarly, from the Lagrangian in Eq.(14), the partition function of the QPSJ state (superinsulating The Lagrangians at zero temperature in Eqs.(1) and state) at zero temperature is as follows: Eq.(2) are as follows: ′0 ′ ′ ZDAXY exp E S + E S Mτ Mx Z DAXY , M 2 ≡ − E0 J 2 n o2 LJJ = [ τ θ (x)] +EJ [1 cos j θ (x)] , (9) − 2 ∇ − ∇ ′ ′0 ′ x=1 j=1 Z DAXY Dθ˜exp E cos τ θ˜(x)+ E S cos j θ˜(x) . (16) X X ≡ S ∇ ∇ x j=1 Z X X M 0 2 2 The partition functions ZAXY and ZDAXY in Eqs.(15) and ES ˜ ˜ LQPS = τ θ (x) +ES 1 cos j θ (x) , (10) (16) are the starting points for the analysis in the follow- − 2 ∇ − ∇ x=1 j=1 X h i X h i ing sections. 0 0 where EJ and ES can be considered as the imaginary time components of the JJ and the QPSJ, respectively, 3. DUAL TRANSFORMATION BETWEEN THE and are defined as follows: JJ MODEL AND THE QPSJ MODEL ON THE ~2 ~2 NANOSHEET 0 0 EJ 2 , ES 2 , (11) ≡ 2a0 Ec ≡ 2a0 EL This section describes the nonlinear Legendre trans- where a τ /M is an imaginary time spacing in the 0 ≡ max τ formation between the Hamiltonian and the Lagrangian time dimension; and τmax and Mτ are an imaginary time with canonical conjugate variables of the infinite order in length and an imaginary time division number, respec- a compact 2+1d lattice space on a nanosheet. Herein, we tively. Eqs.(7) and (11) can be summarized as a relation- show that there is exact duality between the JJ model ship between the Josephson energy and the QPSJ energy and the QPSJ model. For the partition functions ZJJ as follows: in the Lagrangian of Eq.(9), the auxiliary field N (x)
′0 1 ′0 1 by Hubbard-Stratonovich transformation is introduced E J 2 ′ , E S 2 ′ , (12) ≡ 4π E S ≡ 4π E J as follows: where E′ E a /~, E′0 E0 a /~, E′ E a /~ Z = Dθ DNexp iN(x) θ (x) E′ N(x)2 J J 0 J J 0 S S 0 JJ ∇τ − C ′0 ≡0 ~ ≡ ≡ x and E S ES a0/ represent the nondimension energies, Z Z Xn respectively.The≡ first terms of Eq.(9) and (10) are ex- 2 E′ [1 cos θ (x)] , (17) pressed in a quadratic form for the imaginary time dif- − J − ∇j j=1 ference of each phase, however, the rightmost terms in X o Eqs.(9) and (10) are expressed in a cosine form for the It can be seen that the auxiliary field N (x) is the same spatial difference of each phase. Here, the cosine form as the number Nθ (x) of the Cooper pairs introduced in also introduces approximately to the leftmost terms of Eq.(1).The canonical conjugate momentum pθ (x) with Eqs.(9) and (10) in consideration of the periodicity in respect to θ (x) in the Lagrangian of the lattice space of the lattice space as follows: Eq.(9) is defined as follows:.
M 2 ~ 0 0 ipθ(x) i N(x)= a0EJ τ θ (x). (18) LAXY = E 1 cos τ θ (x) +EJ 1 cos j θ (x) , (13) − J − ∇ − ∇ ≡ − ∇ x=1 j=1 X h i X h i On the other hand, for the partition functions ZAXY in Eq.(15), the auxiliary field N(x) by Hubbard- M 2 Stratonovich transformation is introduced as follows: 0 ˜ ˜ LDAXY = ES 1 cos τ θ (x) +ES 1 cos xθ (x) . (14) − − ∇ − ∇ 2 x=1 j=1 Z = Dθ DNexp i2N (x) sin [ θ (x)/2] E′ N(x) X h i X h i AXY ∇τ − C x Z Z Xn where, LAXY and LDAXY are equivalent to LJJ and LQPS, 2 ′ respectively, and the AXY and DAXY model, respec- E J [1 cos j θ(x)] , (19) − − ∇ j=1 tively. From the Lagrangian in Eq.(13) the partition X o 4
Again, it can be seen that the auxiliary field N (x) is Using Eqs.(7), (11) and (24), the following relationships the same as the number Nθ (x) of the Cooper pairs in- are obtained: troduced in Eq.(1). The canonical conjugate momentum 2 2 0 pθ(x) with respect to θ(x) in the Lagrangian of the com- Ec [Nθ (x)] = E 1 cos τ θ (x) =ES 1 cos j θ˜(x) , − J − ∇ − ∇ x x x j=1 pact lattice space of Eq.(13) is defined as follows: X Xh i X X h i (25) ip (x) i~N (x)= 2a E0 sin θ (x)/2 , (20) θ ≡ θ − 0 J ∇τ 2 2 h i 0 EL N˜˜ (x) = E 1 cos τ θ˜(x) =EJ 1 cos j θ (x) . In Eq.(20), if the linear approximation of sin( τ θ/2) θ − S − ∇ − ∇ x x x j=1 ∇ ≃ X h i Xh i X X h i τθ/2 holds, it matches the canonical conjugate mo- (26) ∇mentum in Eq.(18).Therefore, the contents of the curly bracket of the exponential function between Eqs.(15) and Eq.(25) means that the charging energy by the charge 2e (19) can be considered as a nonlinear Legendre transfor- in the JJ is equal to the QPSJ energy, that is, the con- mation introduced by the canonical conjugate momen- densation energy of the magnetic flux. Eq.(26) means tum of Eq.(20). Similarly, for the partition function that the charging energy by the flux quantum Φ in the ˜ 0 ZDAXY in Eq.(16), the auxiliary field N (x) by Hubbard- QPSJ is equal to the JJ energy, that is, the condensation Stratonovich transformation is introduced as follows: energy of the Cooper pair. By establishing the relation- ship between Eqs.(25) and (26), between the Hamiltonian Z = Dθ˜ DN˜exp i2N˜(x) sin θ˜(x)/2 E′ N˜(x)2 DAXY ∇τ − L Eqs.(1) and (3), and, between the Lagrangian Eqs.(13) x Z Z Xn h i and (14), it can be shown that each of them is self- 2 ′ dual. Thus, from the canonical conjugate momentum E S 1 cos j θ˜(x) , (21) − − ∇ of Eqs.(20) and (22), and the relationship between the j=1h io X phase and amplitude of Eq.(6) derived in the previous It can be seen that the auxiliary field N˜ (x) is the same section, an exact duality between the JJ and QPSJ mod- ˜ els has been proven. From the Josephson’s equations as the magnetic flux quantum numbers Nθ˜ (x) in Eq.(3). The canonical conjugate momentump ˜˜ (x) with respect Eq.(2) and the dual Josephson’s equations Eq.(4), the θ 8,14 to θ˜(x) in the Lagrangian of the compact lattice space electrical conductance G can be derived as follows : of Eq.(14) is defined as follows: ˜ dNθ (x) EJ Nθ˜ (x) G = GQ = GQ , (27) − dN˜˜ (x) ES Nθ (x) ~ ˜ 0 ˜ θ ip˜θ˜ (x) i Nθ˜ (x)= 2a0ES sin τ θ (x)/2 , (22) ≡ − ∇ 2 h i where GQ (2e) /h is the quantum conductance. From Therefore, the contents of the curly brackets of the ex- Eq.(27), the≡ following elliptical orbit can be drawn: ponential function between Eqs.(16) and (21) can also be considered as thenonlinear Legendre transformation −1 2 2 1 g Nθ + gN˜ = , introduced by the canonical conjugate momentum of θ˜ η Eq.(22). Since Eqs.(20) and (22) have a half-angle re- 1/2 g (EJ /ES) , η EJ ES/Γ, (28) lationship, a half-angle version of Eq.(6) is introduced as ≡ ≡ follows: where Γ is an arbitrary integral constantp with an energy dimension. From Eq.(27), one can draw an elliptical orbit 2 1 N (x)= sin2 θ˜(x)/2 , as shown in Figure 2. Figure 2 (a) shows a JJ junction θ 2π v ∇j uj=1 state (superconducting state) for EJ ES, Figure 2 (b) uX h i ≫ t shows a quantum Hall state (Bose semiconductor state) 2 1 2 for E = E , and Figure 2 (c) shows a QPS junction state N˜˜ (x)= sin θ (x)/2 , (23) J S θ 2π v ∇j uj=1 (superinsulating state) for EJ ES .The results shown uX h i ≪ 12 15 16 t in Figure 2 are very similar to those in references ’ ’ . Eq.(23) is equivalent to Eq.(6) within the range of linear Compare the result of Eq.(27) with the result of the quan- approximation. From Eqs.(20), (21), and (23), the rela- tum Hall effect shown below:. tionship between the phase θ (x) of the Cooper pair and the phase θ˜(x) of the magnetic flux quantum field is as G = GQν, (29) follows: where ν is the Landau level occupancy and is defined as: 0 2 2 2 4 a0EJ 2 1 2 N (x)= sin θ (x) /2 = sin θ˜(x)/2 , Ne θ − ~2 ∇τ π2 ∇j ν , (30) j=1 h i X h i ≡ Nφ 2 4 a E0 2 ˜ 2 0 S 2 ˜ 1 2 N (x)= sin τ θ (x)/2 = sin j θ (x) /2 , where Ne and Nθ are the electron number and the num- θ˜ − ~2 ∇ π2 ∇ j=1 h i X h i ber of magnetic fluxes (vortex number), respectively. In (24) our model, the correspondence is N N and N N˜˜. e → θ φ → θ 5
1/2 FIG. 3: Schematic phase diagram for g (EJ /ES ) -η (a) EJ ES (b) EJ = ES ≡ ≡ ≫ √EJ ES/Γ. Schematic phase diagram with g on the horizontal axis and η on the vertical axis.
Villain approximation.25-27.First, we apply the Villain ′ approximation to Z AXY introduced in Eq.(15) as follows:
(E′0 ) 2 Z R Dθ exp − J v θ 2πn QV ≡ QV 2 ∇τ − 0 x,τ Z X{n} Xn (E′ ) 2 2 +− J v θ 2πn , (32) 2 ∇j − j j=1 X o (c) EJ ES ≪ where ZQV is the Villain approximation of the parti- ′ 2 ′0 MMτ tion function Z AXY , and RQV [Rv(EJ ) Rv(E J )] is the Villains normalization≡ parameter,, where FIG. 2: Elliptical orbit with the number Nθ of the Cooper −1 Rv(E) 2π(E)vI0(E), and (E)v 2 ln[I0(E)/I1(E)] . pairs versus the number N˜˜ of the magnetic flux quantum ≡ ≡{− } θ I (E) and I (E) represents the modified Bessel func- 0 p 1 tions of order zero and order one, respectively. The 2 ∞ When the Hall conductivity of the quantum Hall effect in summation symbols are used for ≡ the bulk is derived using the Kubo formula, the Landau {n} x,τj=0 nj (x,τ)=−∞ P Q Q P level occupancy ν can also be expressed as the topological the integer fields nj (x, τ). The partition function in quantum number, also called the Chern number. From Eq.(32) is equivalent to the Euclidean version of the quantum vortex dynamics in a film of superfluid helium Eqs.(23), (27), (29) and (30), it can be seen that the 25,28 quantum Hall state is obtained if the following relational introduced by Kleinert . For Eq.(32) the following expression is satisfied: identities associated with the Jacobi theta function are used: ˜ ∞ ∞ EJ 1 cos θ (x) = ES 1 cos θ (x) , (31) E 1 b2 − − exp − (θ 2πn)2 = exp − + ibθ , (33) h i h i 2 − √ 2E n=−∞ l=−∞ 2πE Eq.(31) means that the energy of the Josephson junction X X is equal to the energy of the quantum phase slip. In the As a result, Eq.(32) can be rewritten as follows: vicinity of the region where the Landau level occupancy 2 ν is represented by the Chern number, it is expected that b2 (x) b2 (x) Z =C δ exp − 0 + − j , (34) not only the quantum Hall phase but also a topological QV QV ∇j bj ,0 2(E′0 ) 2(E′ ) x,τ J v j=1 J v insulator phase and a topological metal phase exist. If X{b} X X 1/2 η > 1 and g (EJ /ES) ≦ 1 , the region means Bose where CQV is a normalization parameter defined by 16 19≡ 1/2 2 0 MMτ insulator - , and if η > 1 and g (EJ /ES) ≧ 1, [I (E′ ) I (E′ )] , and b (x) represents auxiliary mag- ≡ 0 J 0 J i the region means Bose metal16,20-24. Figure 3 shows a netic fields with integer values. The integer dual vector 26 schematic phase diagram for g versus η. potentialsa ˜i(x) (i =0, 1, 2) is introduced as follows : b (x)= ε a˜ (x) = ( ˜a) (x), (35) i ijl∇j l ∇× i 4. DUAL TRANSFORMATION FROM THE AXY MODEL TO THE GAUGED DAXY MODEL where εijl is the LeviCivita symbol of three dimensions. BY VILLAIN APPROXIMATION IN A 2+1D By using the dual transformations of Eq.(35), the follow- SYSTEM ing Eq.(36) can be obtained:
˜a 2 2 ( ˜a) 2 ( )0 j In this section, we perform the dual transformation ZQV CQV δ∇ ε ∇ a˜ ,0exp − ∇× + − ∇× , (36) ≡ j ijl j l ′0 2(E′ ) x 2 E J v j=1 J v between the AXY model and the DAXY model from the X{a˜} X X 6
Substituting Poisson’s formula in Eq.(37) to (36): gab is defined as follows:
∞ ∞ ∞ 2 γ 0 0 E0 f (˜a )= Dα˜′ f(α˜′ ) δ exp i2π l˜ α˜′ gab , γ J v j j j ∇j ˜lj ,0 j j , (37) 0 1 0 , (43) ≡ ≡ (EJ ) ˜ x j=0 v (˜aj )=−∞ −∞Z (lj )=−∞ 0 0 1 X X X X then Eq..(38) is as follows: where γ is an anisotropic parameter in the JJ model. ′ Integrating over the continuous-valued gauge fieldsα ˜1 2 ′ α˜′ (x) andα ˜2 for Eq.(42) yields following the partition function: ′ − ∇× 0 ZQV =CQV Dα˜ i δ ˜ exp ∇j lj ,0 ′0 ′ 2 ′0 ′ ′ 2 E ZQV =C QV δ˜ exp 2π (E ) l˜j (x) 0 (x x ) ˜lj (x ) ˜ x J v li(x),0 J v Z li ′ − V − X X {l } x,x { } Xi Xh ˜ 2 2 ′ 2 2 ′ ′ ′ α (x) 2π (E J ) l˜0 (x) 0 (x x ) l˜0 (x ) , (44) − ∇× j v + + i2π ˜l (x) α˜′ (x) , (38) − V − 2(E′ ) j j i j=1 J v j=0 −1 −1 X X where C′ C det ηab ¯ 2 det ¯ 2 , and QV ≡ QV − ∇a∇b −∇0∇0 the anisotropic massless lattice potential (or lattice Poisson’s formula in Eq.(37) converts the integer-valued 25 Green’s function) 0 (x) is defined as: vector potentialsa ˜i to the continuous-valued vector po- V ′ 25,28 1 tentialsα ˜i. Following the quantum vortex dynamics , (x) − (x), (45) ′ V0 ≡ gab ¯ the Euclidean Lagrangian density ofα ˜i is as follows: ∇a∇b Now, from this lattice potential, the ”split lattice po- 2 2 2 2 β˜ (x) β˜ (x) ′ 0 i ′ (x) LQV(x)= + i2π ˜lj(x) α˜ j (x), (39) tential” 0 obtained by dividing the ”core lattice ′0 2(E′ ) − V ′ 2 E J v i=1 J v j=0 potential” (0) and ”split difference operator” are in- X X V0 ∇ i troduced as follows: where, β˜i(i = 1, 2) and β˜0 can be considered as a dual 1 ′ (x) − (x) , electric field and a dual magnetic field in a 2 + 1d dual V 0 ≡ gab ¯′ ′ ∇ a∇ b electromagnetic field, respectively, and are defined as fol- ′ i lows: i ∇ (46) ∇ ≡ 1 (0) gab ¯ −V0 − ∇a∇b β˜ (x) α˜′ (x) α˜′ (x), q 0 ≡ ∇1 2 − ∇2 1 Thus, we have the following identity: ′ ′ β˜1 (x) 2α˜ 0 (x) 0α˜ 2 (x), 2 ′0 ˜ ′ ′ ˜ ′ 2 ′ ˜ ′ ′ ˜ ′ ≡ ∇ − ∇ exp 2π E J lj (x) 0(x x )lj (x ) 2π (E J )vl0 (x) 0 (x x ) l0 (x ) − v V − − V − ˜ ˜′ ˜′ x,x′ β2 (x) 0α 1 (x) 1α 0 (x), (40) Xh i ≡ ∇ − ∇ 2 2 ˜ ′ ˜′ ′ α′ 2 α − ∇ × j =C′ Dα˜′ exp − ∇ × 0 + If the 1,2 componentse ˜ ande ˜ of the dual electric field i 0 ′ , (47) 1 2 2 E′ 2 (E J ) Z x J v j=1 v are set ase ˜1 β˜2, ande ˜2 β˜1, respectively, the dual X X ≡ ≡ − Maxwell’s equations from Lagrangian of Eq.(39) are as 1 1 ′ ¯˜′ ′ 2 ¯˜′ ′ 2 where C det ˜ l det ˜ . Substituting follows: ≡ −∇l∇ −∇0∇0 Eqs.(46) andh (47) in Eq.(44)i h gives:i 1 ˜ [ 1e˜1 (x)+ 2e˜2 (x)] = i2πl0 (x), 1 ′ Z =C′′ Dα˜′ exp f ab ′ α˜′ ′ α˜′ δ (E J ) ∇ ∇ QV QV i − 0 ˜li,0 v 2 E′ ∇ × a ∇ × b x " J v # ˜ Z li 1 1 X {X} ˜ ˜ ′0 2β0 (x) ′ 0e˜1 (x)= i2πl1 (x), 2 (E J )v ∇ − (E J )v ∇ 2 ′0 ˜2 2 ′ ˜2 ˜ ˜′ exp 2π 0 (0) E J lj 2π 0 (0) E J l0 + i2π lj α j , (48) × − V v − V v x j=0 1 1 ˜ ˜ X X ′ 0e˜2 (x) ′0 1β0 (x)= i2πl2 (x), (41) −(E J ) ∇ − (E ) ∇ 1 1 v J v ′′ ′ ¯˜′ ′ 2 ¯˜′ ′ 2 where C C det ˜ l det ˜ 0 , and the met- QV ≡ QV −∇l∇ −∇0∇ For Eq.(38), in order to integrate out for continuous- ric f ab is definedh as:: i h i ′ valued vector potentials α˜ i, the partition function when ˜′ 1 0 0 the axial gauge-fixing condition α 0 = 0 is imposed is as ab f 0 γ 0 . (49) follows: ≡ 0 0 γ ′ ′ 1 ZQV =CQV Dα˜ 1 Dα˜ 2 δ ˜ exp − ∇j lj ,0 ′0 ˜ 2 E J v Furthermore, from Eqs.(15), (32), and (34), the following Z Z li h {X} identity can be obtained: 2 ⊥ ⊥ ⊥ ˜′ ⊥ ′ ˜′ ′ ˜⊥ ˜′ α i (x) Dij (x, x ) α j (x )+i2π lj (x) α i (x) , 2 2 ′0 2 2 ′ 2 × ′ δ exp 2π (0) E ˜l 2π (0) E ˜l + i2π ˜l α˜′ j=0 ˜li,0 0 J j 0 J 0 j j x,x − V v − V v X X i ˜ x j=0 Xli X X ⊥ ′ ⊥ ab ¯ ¯ ⊥ ⊥ { } Dij (x, x ) δij g a b + i j , (42) 2 ≡− ∇ ∇ ∇ ∇ ˜ ˜′ ˜ ˜′ 1 cos j θ 2πα j 1 cos 0θ 2πα 0 j=1 − ∇ − Dθ˜exp( 1) − ∇ − + h i, (50) where the orthogonal symbol as a superscript indicates ≈ − h 2 ′ i P 2 ′0 Z x 4π 0 (0) E J 4π 0 (0) E J ′ ⊥ ′ X V v V v that only the componentsα ˜1 andα ˜2 exist, and the metric 7
As an analogy to the relational expression between the the AXY model by the Villain approximation. First, ap- ′ JJ energy and the QPSJ energy in Eq.(12),the following ply the Villain approximation to Z DAXY introduced in relational expression is obtained from Eq.(50). Eq.(16) as follows:
1 1 ′0 ′ 2 ′0 ′ E S ES E , (E S) , (51) ˜ − v ˜ 2 − v ˜ 2 S v 2 ′ v 2 ′0 ZQDV RQDV Dθ exp τ θ 2πn˜0 + j θ 2πn˜j , (56) ≡ 4π 0(0)(E J ) ≡ 4π 0(0)(E ) ≡ 2 ∇ − 2 ∇ − v J v x j=1 V V Z X{n˜} X X Eqs.(12) and (51) differ only in the coefficients of the core MMτ where R [R (E′ )2R (E′0 )] , and Z repre- lattice potential V0 (0), therefore, Eq.(51) is a reasonable DQV ≡ v S v S QDV result. Substituting Eqs.(50) and (51) in Eq.(48) gives: sents the Villain approximations of the partition func- ′ tion Z DAXY . Using the Jacobi theta function of Eq.(33), f ab ′′ ˜′ ′ ˜′ ′ ˜′ ˜ Eq.(56) can be rewritten as follows: ZQV =C QV Dα i exp − 0 α α Dθ exp( 1) ′ ∇ × a ∇ × b − x " 2 E J v # Z X Z 2 2 2 ˜2 ˜b (x) bj (x) ′0 ˜ ˜′ ′ ˜ ˜′ 0 (E S) 1 cos 0θ 2πα 0 +(E S) 1 cos j θ 2πα j , (52) ZQDV =CQDV δ ˜ exp − + − , (57) × v − ∇ − v − ∇ − ∇j bj ,0 ′0 ′ x j=1 2(E S)v 2(E S)v X h i X h i ˜b x j=1 X{ } X X Eq.(52) represents the gauge-coupled QPSJ partition ′ 2 ′0 MMτ ˜ ˜′ where CQDV [I0(E J ) I0(E J )] , and bi (x) represents function by the dual gauge field α i. For the coefficient ≡ ′0 auxiliary dual magnetic fields with integer values. Inte- E J v of the dual gauge field energy, use the relation- ships in Eqs.(51) and (7), and next, when these non- ger vector potentials ai (x)(i =0, 1, 2) are introduced as dimensional energy constants have been converted to the follows : original energy dimension, Eq. (and next, when these ˜bi(x)=εijl j al(x) = ( a) (x), (58) non-dimensional energy constants have been converted ∇ ∇× i to the original energy dimension, Eq. (52) will be as By using the dual transformations of Eq.(58), the follow- follows:) will be as follows: ing Eq.(59) is obtained:
2 ab a0 (2e) f 2 2 Z =C′′ Dα˜′ exp (0) ′ α˜′ ′ α˜′ Dθ˜ a 2 ( a) QV QV i −~ 0 ( )0 j x V 2C ∇ × a ∇ × b − ∇× Z τ, " #Z ZQDV =CQDV δ∇j εijl∇j al,0exp − ∇×0 + ′ , (59) X 2 E′ 2(E S) {a} x S v j=1 v a 2 X X X 0 0 ˜ ˜′ ˜ ˜′ exp − (E )v 1 cos 0θ 2πα 0 +(ES)v 1 cos j θ 2πα j , (53) × ~ S − ∇ − − ∇ − τ,x j=1 X h i Xh i Using Poisson’s formula in the following Eq.(37) for Furthermore, by scaling 2e to the non-dimensional dual Eq.(59): ′ gauge field α˜i, we introduce a new dual gauge fieldα ˜i as 2 1 2 Z =C Dα′ δ exp − ( α′) follows: QDV QDV i ∇j lj ,0 2(E′ ) ∇× j x S v j=1 Z X{li} Xh X ′ 2eα˜ i (x) α˜i (x), (54) 2 ≡ 1 ′ 2 ′ + − ( α ) + i2π lj α j , (60) 2 E′0 ∇× 0 By the transformation of Eq.(54), Eq.(53) is transformed S v j=0 X i as follows: The Euclidean Lagrangian density of α′ is as follows: ab j ′′ a0 f a0 ZQV =C QV Dα˜i exp − ( α˜) ( α˜) Dθ˜exp − ~ 2µ˜′′ ∇× a ∇× b ~ 2 2 τ,x 1 1 Z X Z 2 2 ′ LQDV(x)= ′ βi (x)+ ′0 β0 (x) i2π lj (x) α j (x), (61) 2 2(E ) 2(E ) − S v i=1 S v j=0 0 ˜ ˜ E 1 cos 0θ qmα˜0 +ES 1 cos j θ qmα˜j , (55) X X × S − ∇ − − ∇ − τ,x j=1 X h i X h i where,βi(i =1, 2) and β0 can be considered as an electric where µ˜′′ C/V (0), and q 2π/2e=Φ /~ represents a unit field and a magnetic field in a 2+1d electromagnetic field, ≡ 0 m ≡ 0 magnetic charge. Eq.(55) shows that the pure AXY respectively, and are defined as follows: model (JJ model without gauge coupling) has been dual β (x) α′ (x) α′ (x), transformed to the gauged DAXY model (gauged QPSJ 0 ≡ ∇1 2 − ∇2 1 β (x) α′ (x) α′ (x), model) by the Villain approximation in the 2+1d system. 1 ≡ ∇2 0 − ∇0 2 ′ ′ In other words, the gauge QPSJ model is a frozen lattice β2 (x) 0α 1 (x) 1α 0 (x), (62) ≡ ∇ − ∇ dual superconductor25,29-31, for the JJ model without gauge coupling. If the 1,2 components e1 and e2 of the electric field are set as e1 β2, and e2 β1, respectively, the Maxwell’s equations≡ from the Lagrangian≡− of Eq.(61) are as follows: 5. DUAL TRANSFORMATION FROM THE 1 DAXY MODEL TO THE GAUGED AXY MODEL [ 1e1 (x)+ 2e2 (x)] = i2πl0 (x), (E′ ) ∇ ∇ BY VILLAIN APPROXIMATION IN A 2 + 1 D S v SYSTEM 1 1 ′0 2β0 (x) ′ 0e1 (x)= i2πl1 (x), (E S)v ∇ − (E S)v ∇ In this section, contrary to the previous section, we 1 1 ′ 0e2 (x) ′0 1β0 (x)= i2πl2 (x), (63) show the dual transformation from the DAXY model to −(E S)v ∇ − (E S)v ∇ 8
′ Integrating over the continuous-valued gauge fields α 1 Furthermore, by scaling Φ0 to the non-dimensional gauge ′ ′ and α 2 for Eq.(60) yields the following partition function: field α i, we introduce a new gauge field αi as follows:
′ 2 ′0 ˜ ′ ′ ′ ZQDV =C QDV δli(x),0exp 2π E J lj (x) 0(x x )lj (x ) Φ0α i (x) αi (x), (70) ′ − v V − ≡ X{li} Xx,x h 2 ′ ′ ′ Eq.(69) is transformed as follows: 2π E J l0 (x) ˜0(x x )l0 (x ) , (64) − v V − − i − ′′ a0 1 ab ′ ′ a0 1 1 ZQDV =C QDV Dαi exp − f ( α) ( α) Dθexp − ′ ab 2 2 ~ 2µ ∇ × a ∇ × b ~ where C C det η ¯ det ¯ , τ,x QDV ≡ QDV − ∇a∇b −∇0∇0 Z X Z and the anisotropic massless lattice potential (or lattice 2 0 (E )v 1 cos( 0θ 2qα0) +(EJ )v 1 cos( j θ 2qαj) , (71) ˜ × J − ∇ − − ∇ − Green’s function) V0 (x) is defined as: τ,x j=1 X h i X h i 1 ˜ ˜ 0 (x) ab− (x) , where µ L/ 0 (0), 2q 2e/~=2π/Φ0 represents a unit V ≡ g˜ ¯ a b ∇ ∇ Cooper≡ pairV charge, which≡ is twice the unit charge 0 γ˜ 0 0 E q e/~=2π/2Φ0. Eq.(70) shows that the pure DAXY g˜ab , γ˜ S v , (65) ≡ ≡ 0 10 ≡ (E ) model (JJ model without gauge coupling) has been 0 01 S v dual transformed to the gauged AXY model (gauged JJ From this lattice potential, we introduce a ”split lattice model) by the Villain approximation in the 2 + 1d sys- ′ tem. In other words, the gauge JJ model is a frozen potential” ˜ 0 (x) ˜0 (x) ˜0 (0) δx,0 which is obtained by V ≡ V −V lattice dual superconductor for the QPS model without dividing the ”core lattice potential” ˜0 (0) and ”split dif- ference operator” ′ , as follows: V gauge coupling. ∇ j 1 ˜′ (x) − (x) , 0 ab ¯′ ′ V ≡ g˜ a b 6. MEAN FIELD ANALYSIS OF THE GAUGED ∇ ∇ QPSJ MODEL ON THE NANOSHEET
′ i i ∇ ∇ ≡ ab 1 ˜0 (0) g˜ ¯ a b In this section, we introduce the mean field approxima- − V − ∇ ∇ tion to the partition function of Eq.(55) and discuss its q Similar from Eq.(47) to (51) in the previous section, the phase transition. From Eq. (55), the partition function partition function ZDQV can be written as follows: excluding the constant part is newly defined as ZGQPJ of the gauged QPSJ model, and using unit vectors of two ′′ ′ 1 ab ′ ′ ′ ′ ˜ ZQDV =C QDV Dα i exp − f ( α ) ( α ) Dθ exp( 1) ˜ ˜ ˜ ′0 ∇ × a ∇ × b − real components Ui = [cos θ, sin θ], Eq. (55) is rewritten x "2 E S v # Z X Z 2 as follows: ′0 ˜ ′ ′ ˜ ′ (E )v 1 cos 0θ 2πα 0 +(E J )v 1 cos j θ 2πα j , (66) × J − ∇ − − ∇ − 1 x j=1 ab ′ ˜ ′ ˜ h i h i ZGQPJ = Dα˜i (x)exp − f ( α)a( α)b (x) X X 2µ˜′ ∇ × ∇ × Z x 1 1 X ′′ ′ ¯˜′ ˜′ 2 ¯˜′ ˜′ 2 2 where C QDV C QDV det l l det 0 0 , and the ′ ≡ −∇ ∇ −∇ ∇ Dθ˜exp E S d˜ U˜l (x) R˜α˜U˜l (x) (72) ˜ab × metric f is definedh as: i h i x ( x ) Z X X Xl=1 ˜′ ~ ˜′′ 1 0 0 where µ µ /a0, and the lattice difference operator ab ≡ f˜ 0γ ˜ 0 (67) R˜ is defined as: ≡ α˜ 0 0γ ˜ 2 ′0 1 ¯ ¯ E S R˜α˜ 1+ D˜iD˜ i +˜γD˜τ D˜τ , γ˜ , (73) ′0 ′ ≡ ˜ ≡ E′ where (E J )v and (E J )v are defined as follows: 2d i=1 ! S X ′0 1 ′ 1 (E )v , (E J )v (68) where d˜ 2+˜γ represents anisotropic dimensional con- J ≡ 4π2 ˜ (0)(E′ ) ≡ 4π2 ˜ (0) E′0 0 S v 0 S v ≡ ¯ V V stants of the gauged QPSJ; and D˜i and D˜ i are for- If 0 (0) ˜0 (0) holds, Eqs.(68) and (51) are completely ward and backward covariant lattice derivatives, respec- V ≡V equivalent. Eq.(66) represents the gauge-coupled JJ par- tively, and are defined, for example, for a complex field ′ tition function by the gauge field α j . For the coeffi- U˜ = U˜1 +iU˜2, as follows: (E′0 ) cient S v of the gauge field energy, use the relation- x D˜ U˜ (x, τ) U˜ (x + i, τ) e−iqmα˜i( ,τ) U˜ (x, τ) ,, ship between Eqs.(68) and (7), and then, when these i ≡ − ¯ iqmα˜i(x−i,τ) non-dimensional energy constants are converted to the D˜ iU˜ (x, τ) U˜ (x, τ) U˜ (x i, τ) e , (74) original energy dimension, Eq. (68) will be as follows: ≡ − − ˜ ˜¯ 2 the same applies to D and D . In Eq.(72), we introduce a (Φ ) τ τ Z =C′′ Dα′ exp − 0 ˜ (0)(0) 0 f ab( ′ α′) ( ′ α′) Dθ QDV QDV i ~ V0 2L ∇ × a ∇ × b ˜ τ,x " # two sets of real two-component fieldsu ˜l and ψl (l =1, 2) Z X Z 2 which satisfy the following identity: a0 0 ′ ′ exp − (EJ )v 1 cos( 0θ 2πα0) +(EJ )v 1 cos( j θ 2παj ) , × ~ − ∇ − − ∇ − ∞ ∞ τ,x j=1 ˜ ˜ X h i Xh i dψ1dψ2 ˜ ˜ d˜u1d˜u2 2 exp ψl u˜1 Ul =1, (75) (69) −∞ −∞ (2πi) − − Z Z n o 9
2 ∞ ˜ 1 ab ′ ′ dψl ZGQPJ = Dα˜iexp − f ( α˜) ( α˜) d˜ul ˜′ ∇ × a ∇ × b 2πi x 2µ x −∞ Z X YYl=1 Z 2 exp E′ d˜u˜ R˜ u˜ ψ˜ u˜ + ln I ( ψ˜ ) , (76) × S l α˜ l − l l 0 | l| x FIG. 4: Schematic diagram of a superinsulator on a nanosheet XXl=1h i at zero temperature. where we have defined the functional integrals of θ˜ as follows: Cooper pair and the anti-Cooper pair12, 13 is confined π dθ˜ 2 within the superinsulator by the electric flux-tubes. exp ψ˜ U˜ =exp ln I ( ψ˜ ) , (77) 2π l l 0 | l| x Z−π ( x l=1 ) x,l Figure 4 shows the confinement of electric flux in the Y X X X n o superinsulator between the Cooper pair and the anti- where I ( ψ˜ ) ( ψ˜ ψ˜2 + ψ˜2 ) represents the modified Cooper pair. Whether this confinement picture in which 0 | | | |≡ 1 2 the charge of 2e for a superinsulator is the smallest unit Bessel functions of aq zeroth-order integer. In Eq.(76), of charge is correct will be discussed again in Section 8 performing the integrals overu ˜ fields, we obtain the par- l (Summary and discussion). From Eq.(81), three DGL tition function by the complex field ψ˜ ψ˜1 + iψ˜2 and ∗ ≡ equations are derived as follows: ψ˜ ψ˜1 iψ˜2. ≡ − ~ 2 i 1 ′ia ′ ˜ Φ0 ∗ ∗ Φ0 2 ∞ ˜ ˜∗ j = η ∂ B = ψ˜ ∂ ψ˜ ∂ ψ˜ ψ˜ ψ˜ A˜ dψdψ ′ ∗ Φ i i i, (82) ZGQPJ = Dα˜i exp F ψ,˜ ψ˜ , α˜i µ˜ × a 2mΦ − − mΦ | | ′ ˜ − x −∞ 4πE Sd! n o Z Y Z n o 2 ′ ˜ ˜∗ 1 ab ′ ˜ ′ ˜ 1 ˜ˆ ˜ˆ F ψ, ψ , α˜i f ( α)a( α)b + ′ ψ ln I0( ψ ) 2 ≡ ˜′ ∇ × ∇ × 4E d| | − | | 1 0 γ~Φ Φ x 2µ S 0 ′ a ′ B˜ 0 ˜∗ ˜ ˜∗ ˜ 0 ˜ 2 ˜ X jΦ = η ∂ = ψ ∂0ψ ∂0ψ ψ ψ A0, (83) ˆ 1 µ˜ × a 2mΦ − − mΦ | | ψ˜ (x) R˜ 2 ψ˜ (x), (78) n o ≡ α˜
∗ 2 2 2 In Eq.(78), since ψ˜ and ψ˜ can be regarded as the or- 1 γ 2 i~∂i Φ0A˜i ψ˜+ i~∂0 Φ0A˜0 ψ˜+2β ψ˜ ψ˜= αεψ˜, 2m − − 2m − − | | − der parameters of the superinsulator (i.e., the disorder Φi=1 Φ X parameters of the superconductor), the non-dimensional (84) ′ ∗ free energy F (ψ,˜ ψ˜ , α˜i) can be Landau expanded for 4 2 2 25 where B˜ (x)= ε ∂′ A˜ (x) represents a dual magnetic terms up to ψ , Diψ and Dτ ψ as follows : i iab a b | | | | | | flux density. Eq.(82) and (83) represent the current
′ ∗ 1 ab ′ ′ density and Eq.(84) represents the nonlinear Schrdinger F DGL(ψ,˜ ψ˜ , α˜i) f ( α˜) ( α˜) ≡ ˜′ ∇ × a ∇ × b x 2µ equation for superinsulators. The thermodynamic crit- X 2 ˜ 1 2 2 1 1 2 1 4 ical dual magnetic field Hc for the dual magnetic field + Diψ˜ +˜γ Dτ ψ˜ + 1 ψ˜ + ψ˜ , (79) ˜ ′ ˜ H˜ = B˜ /µ˜ is as follows: 8d | | ! 4 E Sd − | | 64| | ) 0 0 i=1 X
′ 2 2 ′MF F DGLis the non-dimensional DGL energy of the superin- ˜ ′ α ε α E S Hc (E S )= = 1 ′ , (85) sulator or QPSJ model on the nanosheet in d˜ 2+˜γ s 2βµ˜ √2βµ˜ − E S ! dimension at zero temperature. Therefore, the≡ critical ′MF where H˜c is the dimension of the voltage. Both have a values E S according to the mean field approximation ′ dimension of length, the penetration depth λ˜(E′ ) and of QPS amplitude E S are as follows: S ′ the coherent length ξ˜(E S) in the super insulator, which 1 1 E′MF = , (80) both have the dimension of length, are as follows: S ≡ d˜ 2+˜γ 1 βm E′ 2 The continuous limit in Eq.(79), is as follows: λ˜ (E′ )= Φ S S 2 ′ ′MF , (86) ab 2 2 s2πµα˜ Φ0 E S E S ˜ ˜∗ ˜ 2 f ′ A˜ ′ A˜ 1 ~ ˜ ˜ FDGL(ψ, ψ , Ai) dx0 d x (∂ )a(∂ )b + i ∂i Φ0Ai ψ − ≡ 2˜µ × × 2mΦ − − Z Z ( i=1 2 X γ˜ ~ ˜ ˜ ˜ 2 ˜ ˜ 4 + i ∂0 Φ0A0 ψ +α ˜ε˜ ψ + β ψ 1 2mΦ − − | | | | ~ ′ 2 ′ E S ξ˜(E )= S ′ ′MF . (87) ~ ′MF ′ ~ √2mΦα E S E S E S E S ˜ − α˜ 2 , ε˜ ′− , β 2 , (81) ≡ 4a0a ≡ E S ≡ 64a0a From Eqs.(86) and (87), the dual GinzburgLandau pa- 2 rameter κ is as follows: whereµ ˜ C/a 0 (0), mΦ 4d˜~a0is a pseudo-mass of mag- ≡ V ≡ 2 netic flux having a dimension of J s , and A˜µ repre- ˜ ′ ′ λ (E S) mΦ β 2π ′ 2 ′ · κ˜ (E S)= = = √8πB˜c(E S) λ˜ (E S), (88) sents dual vector potentials having a dimension of [C/m]. ˜ ′ ~Φ πµ˜ 2e ξ (E S) 0 s Therefore, the order parameter ψ˜ of the superinsulator is ˜ gauge coupled to the U(1) dual gauge field Aµ by a unit where B˜c =µ ˜H˜c represents the thermodynamic criti- magnetic charge qm 2π/2e=Φ0/~, and when ψ˜ is in the cal dual magnetic flux density. From the analogy of the condensed state, As≡ shown in Figure 4, the pair of the classification of type I and type II superconductors, the 10 following classifications by type I and type II superinsu- lators are formed from the sign of the surface energyσ ˜SN and the value of the GinzburgLandau parameter κ at the superinsulatornormal insulator boundary. i) type-I superinsulator
1 κ<˜ , σ˜SN >0. √2 ii) intermediate of type I and type II superinsulator
1 κ˜ = , σ˜SN =0. √2 iii) typeII superinsulator
1 κ>˜ , σ˜SN <0. FIG. 5: The critical value and the tricritical point as a func- √2 tion of anisotropy parameter in mean field approximation with up to two loop corrections. Where,γ ˜ is the anisotropy From the analogy with the mixed state of the type-II 2 parameter, and (E′ ) loop and (E′ )tri are the critical value superconductor, the possibility of the existence of the S c S mixed state in the case of the type-II superinsulator is and the tricritical point of QPS amplitud, respectively. expected, and, from the analogy with the Abrikosov mag- netic flux lattice of the type-II superconductor, the exis- tence of an electric flux lattice is also expected in the case of the type-II superinsulator on a nanosheet.
7. ESTIMATING THE CRITICAL VALUE OF THE QPS AMPLITUDE BY THE EFFECTIVE ENERGY APPROACH
In the previous section, we derived the thermodynamic critical dual magnetic field H˜c, the penetration depth λ˜, and the coherent length ξ˜ from mean field analysis, all of which depended on the difference between the QPS am- ′ ′MF plitude E S and its mean field critical value E S . Since the mean field approximation is a very rough approxi- FIG. 6: The critical value of the charging energy for various mation, in this section, we show the calculation result of values of the anisotropy parameter γ in Josephson junction the mean energy approximation with the contribution of (JJ) models as a function ofγ ˜. fluctuations up to the two loop corrections by the effec- tive energy approach shows in Appendix B. The critical ′ 2loop ′ 8. SUMMARY AND DISCUSSION value (E S)c of the QPS amplitude E S in the mean field approximation with up to the two loop corrections is given to Eq.(B14) as a function of the anisotropy pa- The conclusions of this paper are summarized below. rameterγ ˜. The results are plotted in Figure 5. Similarly, First, using the dual Hamiltonian method, the phase and ′ tri Figure 5 shows the results for the tricritical point (E S) amplitude relationship between the JJ and QPSJ models of the QPS amplitude given in Eq.(B15). From Eqs. (51) without gauge coupling on 2+1d nanosheets at zero tem- ′ and (11), the QPS amplitude E S has the following rela- perature was determined, and the relationships between ′ tionship with the charging energy E c : various constants were derived. Furthermore, the exact duality between the JJ model and the QPSJ model based ′ 2 ′ E c =2π 0 (0)(E S)v, (89) on the nonlinear Legendre transformation between the V Lagrangian and the Hamiltonian using canonical conju- where 0 (0) is a massless lattice Green’s function having gate variables of infinite order in a compact 2+1d lattice an anisotropyV parameter γ in the JJ model defined by space was demonstrated. A dual transformation from the Eq.(43).Using Eq.(69), we have plotted in Figure 6 the AXY model to the gauged DAXY model by the Villain ′ ′ critical value (E c)c of the charging energy E c for various approximation in the 2+1d system was derived.there are values of the anisotropy parameter γ in JJ models as a two main differences between the dual transformation by function ofγ ˜ the dual Hamiltonian method and the dual transforma- 11
ab 2 2 tion by the Villain approximation. One is that Eqs.(12) ˜ ˜∗ ˜ 2 f ′ A˜ ′ A˜ 1 ~ ˜ ˜ FDGL(ψ, ψ , Ai) dx0 d x (∂ )a(∂ )b + i ∂i 2Φ0Ai ψ ≡ ( 2˜µ × × 2mΦ − − and (51) differ only in the core potential 0 (0). Another Z Z i=1 X V difference is that, in the case of the dual transformation 2 γ˜ ˜ ˜ 2 ˜ ˜ 4 + i~∂0 2Φ0A˜0 ψ +˜αε˜ ψ + β ψ (92) by the Villain approximation, there is a gauge coupling 2mΦ − − | | | | ˜′ by the dual gauge field α i, however, in the dual trans- formation by the dual Hamiltonian method, there is no where m2Φ is the effective mass of the flux pair (vortex gauge coupling by the dual gauge field. The gauge cou- pair in the same rotation direction = superinsulator) and ′ ˜ pling by the dual gauge field α˜ i is gauge coupled with ψ is the wave function of the flux pair. On the other hand, the U(1) dual gauge field by the unit magnetic charge the GinzburgLandau free energy for Eq.(71) is as follows: ˜ab 2 2 qm 2π/2e by the scaling of 2e introduced in Eq.(54), ∗ 2 f ′ A ′ A 1 ~ ˜ FGL(ψ, ψ , Ai) dx0 d x (∂ )a(∂ )b + ( i ∂i Φ0Ai) ψ ≡ ≡ ( 2µ × × 2mc − − and, as shown in Figure 4, it was shown that the electric Z Z i=1 2 X γ 2 4 + ( i~∂0 Φ0A0) ψ˜ +αε ˜ ψ + β˜ ψ flux of the Cooper pair and the anti-Cooper pair in units 2mc − − | | | | of 2e was confined in the superinsulator. In the following,
′MF ′ we consider whether the superinsulator’s picture of con- E J E J ε ′− , (93) finement in 2e units that is, confinement by a pair consist- ≡ E J ing of a Cooper pair and an anti-Cooper pair is correct where m is the effective mass of the Cooper pair and ψ or incorrect. In the case of the magnetic flux confine- c is the wave function of the Cooper pair. It is known that ment for the superconductor, the magnetic flux quantum the microscopic theory of superconductivity with respect Φ = h/2e, which is the minimum unit of magnetic flux, 0 to the GinzburgLandau theory can be described by the is confined. However, if Eq.(55) is correct, in the confine- following BCS Hamiltonian (HBCS ): ment of charges in the superinsulator, 2e = h/Φ0, which is twice the elementary charge e = h/(2Φ0),which is the 3 + 1 ~ 2 HBCS d xϕσ (x) [i ∂j eAj(x)] µ ϕσ(x) minimum unit of charge, is confined. This is clearly in- ≡ 2me − − consistent with the superconducting case. Therefore, in σX=↑,↓ΩZ g 3 + + ′ the case of a superinsulator, it should also be considered | | d xϕσ (x) ϕσ′ (x) ϕσ (x) ϕσ(x), (94) − 2 ′ correct to assume that confinement occurs in units of the Xσ,σ ΩZ elementary charge e = h/(2Φ0), which is the minimum unit of charge. In other words, since the scaling by 2e in where me is the effective mass of the electron, and ϕσ (x) Eq.(54) in Section 4 was completely artificial, one could is the electron field having a spin subscript σ and is a change Eq. (54) to scaling by e as follows: fermion satisfying the anti-commutation relation + ′ 3 ′ ϕ (x) , ϕ ′ (x ) = δ ′ δ (x x ) . eα˜′ (x) α˜ (x), (90) σ σ + σσ i ≡ i − From Eqs.(93) and (94), in microscopic theory, the mat- By the transformation of Eq.(90), Eq.(55) is transformed ter field ϕ is a fermion field which is coupled to the gauge as follows: σ field by an elementary charge e. On the other hand, in ˜ ′′ a0 1 ab a0 the GinzburgLandau theory, the material field ψ is a bo- ZQV =C QV Dα˜i exp − f ( α˜) ( α˜) Dθ˜exp − ~ ˜′′ ∇× a ∇× b ~ τ,x 2µ son field which is coupled to the gauge field by 2e. The Z X Z 2 magnetic flux quantum in the superconductor have a re- 0 ˜ ˜ E 1 cos 0θ 2qmα˜0 +ES 1 cos j θ 2qmα˜j , (91) × S − ∇ − − ∇ − pulsive force, form a vortex lattice, and never intersect τ,x j=1 X h i X h i with each other, so they can be regarded as elementary excitations of fermions.Therefore, from the analogy be- When Eq.(91) is compared with Eq.(55), the coupling tween the BCS theory and the GinzburgLandau theory magnetic charge is 2qm = 2Φ0/~ =2π/e, which is twice ~ of superconductivity, it is expected that the microscopic the unit magnetic charge qm = Φ0/ = 2π/2e. Fig- theory of superinsulators with respect to the DGL theory ure 7 shows the confinement of the electric flux be- in Eq.(92) can be described by the following dual BCS tween the positive and the negative elementary charges Hamiltonian: in the superinsulator, that is, it differs from Figure 2 4. To obtain the result of Eq.(91), Ec introduced 3 + 1 ~ ˜ HDBCS d xϕ˜σ˜ (x) i ∂j Φ0Aj (x) µ˜ ϕ˜σ˜ (x) ≡ 2mΦ − − σX=↑,↓ΩZ h i g˜ 3 + + ′ | | d xϕ˜σ˜ (x)ϕ ˜σ˜′ (x)ϕ ˜σ˜ (x)ϕ ˜σ˜ (x), (95) − 2 ′ σ,X˜ σ˜ ΩZ
where mΦ is the magnetic flux (vortex) quantum, andϕ ˜σ FIG. 7: Schematic diagram of a superinsulator with an ele- is the magnetic flux quantum field having the pseudo- mentary charge e on a nanosheet at zero temperature. spin subscriptσ ˜, and is a fermion satisfying the anti- commutation relationship in Eq.(1) is not (2e)2/2C but rather e2/2C. In this + ′ 3 ′ case, the DGL free energy for Eq.(91) is as follows: ϕ˜σ˜ (x) , ϕ˜ ′ (x ) = δσ˜σ˜′ δ (x x ) . σ˜ + − 12
Similar to the relationship between Eqs.(93) and (94), in the microscopic theory of superinsulators of Eq.(94), the TABLE I: Values of the anisotropic hopping coefficient hn for values of n up to 10 for values of the anisotropic parameter material fieldϕ ˜σ is a fermion field which is coupled to γ˜ = 0 to 1 in the case of 2 < d˜ 3, where d˜ 2 +γ ˜. the gauge field by the flux quantum Φ0. On the other ≤ ≡ hand, in the DGL theory of Eq.(92), the material field γ˜ h2 h4 h6 h8 h10 ψ˜ is a boson field and is coupled to the gauge field by 2Φ0. The construction of the dual BCS theory, which 0 4 36 400 4900 63504 is a microscopic theory of superinsulators on nanosheets, not only elucidates the microscopic mechanisms for su- 0.1 4.02 36.003 410.83602 5125.514241 67964.55133 perinsulators and quantum phase slips, but also the mi- croscopic mechanisms for quark confinement and asymp- 0.2 4.08 36.048 443.77728 5820.335539 81960.10908 totic freedom. This theory is expected to play an impor- 0.3 4.18 36.243 500.13058 7040.109553 107387.453 tant role as a powerful model of the QCD phenomenon 12,13 on desktop . 0.4 4.32 36.768 582.09792 8880.292915 147593.9992
0.5 4.5 37.875 692.8125 11480.27344 207665.9648 9. ACKNOWLEDGMENTS 0.6 4.72 39.888 836.38912 15029.23717 294849.9426 I would like to thank all the faculty and staff of Aichi University of Technology. 0.7 4.98 43.203 1017.98898 19773.88112 419126.4121 0.8 5.28 48.288 1243.89888 26028.09861 593959.5603 Appendix A: Anisotropic lattice Green’s function 0.9 5.62 55.683 1521.62482 34184.79254 837254.033 m (0) at the source x = 0 V 1.0 6 66 1860 44730 1172556 Perform a numerical evaluation of the anisotropic mas- sive lattice Green’s function (lattice potential) m (0) at the origin x = 0. V anisotropic massless lattice Green’s function 0 (0) at the ∞ origin x=0 is evaluated as the sum of the power-seriesV up 0 1 hn m ( ) µν 2 = n+1 , (A1) −g ∇¯ µ∇ν +m 2 ˜ to n =10, it becomes as shown in TABLEIII and Figure V ≡ n=0 [ m + 2d ] P 6: where, d˜ 2+˜γ is the anisotropic dimension, and hn is the anisotropic≡ hopping coefficient of the anisotropic massive lattice Green’s function m (0) at the origin x = 0, and is introduced as follows: V
n − γ˜n j hn n! − 2 Hj , (A2) ≡ [ (n j) ]! j! j=0,2,4 { 2 } P 25 where Hn represents the isotropic hopping coefficients , for example, in the two-dimensional case, H0 = 1, H2 = 4, H6 = 36, H8 = 400,, and in the three-dimensional case, H0 = 1, H2 = 6, H6 = 90, H8 = 1860,,. TABLE I lists examples ofγ ˜ = 0.1, 0.2, ...., 0.8, 0.9, and 1.0. The asymptotic behavior of the anisotropic massive lattice Green’s function m(0) at small values of m in the case of 2 < d 3 is asV follows: ≤ FIG. 8: Asymptotic massless lattice Green’s function 0 (0) at the origin x = 0 versus the anisotropic parameterγ ˜.V −1/2 ˜ 2 ˜ √ 2 ∞ 0 γ˜ 2d m +4d m ∆hn m ( )= − + n+1 V 4π 2 ˜ p p m +2d n=0,2,4.. m2 +2d˜ X p n+1 n ˜ γ˜n−j γ˜−1/2 2d˜ (2n)! 2d ∆hn n! Hj , (A3) ≡ 2 − 4π 2n 2 n +1 Appendix B: Appendix B. Effective energy j=0,2,4 [(n j)/2 ]! j! p 2 (n!) X { − } approach for the QPSJ model TABLE II shows the asymptotic anisotropic hopping co- efficient ∆hn for values of n up to 10 for values of the Consider the order parameter representation of the su- anisotropic parameterγ ˜ =0to1inthecaseof2 < d 3. per insulator by the partition function of Eq.(78). To From Eq.(A3) and TABLE II , when the value of≤ the simplify the problem, we deal with the case where there 13
below: TABLE II: Asymptotic anisotropic hopping coefficient ∆hn for values of n up to 10 for values of the anisotropic parameter −1 ψ˜ G ˜ (x1, x2)ab ˆ = V2 (x1, x2)ab ˆ γ˜ = 0 to 1 in the case of 2 < d˜ 3, where d˜ 2 +γ ˜. ψ ψ˜=α ψ˜=α ≤ ≡ δ2F ′ ˆ = E′ , ψ˜ (x) α˜ , (B1) γ˜ ∆h2 ∆h4 ∆h6 ∆h8 ∆h10 S ˆ ˆ l l δψ˜ (x ) δψ˜ (x ) ˆ ≡ a 1 b 2 ψ˜=α D E
0.1 -1.16603 -0.7561 -0.3789 27.6816 549.214 ˜ˆ where αl is an expectation of ψl (x), and, in general, all 0.2 -0.64231 0.10562 4.26659 59.7304 786.122 n-th order one-particle irreducible graphs involving the vertex functions can be computed as follows: 0.3 -0.4334 0.38867 5.27033 62.4803 769.857 n ′ ψ˜ ′ δ F Vn (x1, x2 , xn)ab =E S , (B2) 0.4 -0.32319 0.50921 5.42088 60.563 744.348 · ·· a1,a2···,an ˜ˆ ˜ˆ ˜ˆ δψa (x1) δψa (x2) δψa (xn) 1 2 ··· n ψ=α
0.5 -0.25823 0.56803 5.36908 59.2432 755.792 The 2 2 matrix G−1(x , x ) defined in Eq.(A1) is × ψ˜ 1 2 ab −1 L 0.6 -0.2182 0.60248 5.34728 60.225 818.965 divided into longitudinal Gψ (x, y) and transverse −1 T Gψ (x, y) parts, which are respectively parallel and or- 0.7 -0.19353 0.6296 5.46036 64.2966 944.418 thogonal to the expected value of the field αl, as follows:
0.8 -0.17904 0.65818 5.7662 72.0562 1145.97 −1 1 1 µν G (x1, x2) = δabδx ,x 1+ g ¯ µ ν (x1, x2) ηab ψ˜ ab ′ ˜ 1 2− ˜ ∇ ∇ 2E Sd 2d 0.9 -0.17169 0.69306 6.30443 84.1746 1443.87 = G−1(x , x )LP L + G−1(x , x )T P T ψ˜ 1 2 ab ψ˜ 1 2 ab,
1.0 -0.16955 0.73705 7.10848 101.516 1866.95 −1 L 1 ′ 1 µν G (x, y) δxy 2E Sdη˜ L 1+ g ¯ µ ν (x, y) , ψ ≡ ′ ˜ − ˜ ∇ ∇ 2E Sd 2d −1 T 1 ′ 1 µν G (x, y) δxy 2E Sdη˜ T 1+ g ¯ µ ν (x, y) , (B3) ψ ≡ ′ ˜ − ˜ ∇ ∇ 2E Sd 2d L T where ηab is a 2 2 matrix; and Pab and Pab are longitudi- TABLE III: Asymptotic massless lattice Green’s function nal and transverse× projection matrices, respectively: 0 (0) at the origin x=0 for values of n up to 10 for values of V the anisotropic parameterγ ˜ = 0 to 1 in the case of 2 1loop 1 1 0.4 0.331151594 F ′ = T r log − gµν ¯ 2MM ˜ ∇µ∇ν τ 2d 2 2 0.5 0.311864992 m m 1 µν +T r log 1 g ¯ µ ν , ˜ − − ˜ ˜ ∇ ∇ 2d 2d 2d 0.6 0.29608701 2 ′ 2 2 ˜ ˜ ′ I1 (E S) m 4d 2d E S 1 ′ , (B5) ≡ − ( − I0 (E S) ) 0.7 0.282786758 where I0 (E) and I1 (E) are order zero and order one the 0.8 0.271328982 modified Bessel functions, respectively. The first trace means zero mass fluctuations (Goldstone modes) and the 0.9 0.261294231 second trace represents massive fluctuations and can be 1.0 0.252390927 calculated by hopping expansion as follows: m2 m2 1 T r log 1 gµν ¯ ˜ − − ˜ ˜ ∇µ∇ν 2d 2d 2d ∞ 1 m2 n h = 1 n , (B6) − n − ˜ ˜ n n=2 2d (2d) is no coupling of gauge fields. First, to derive the one- X loop effective field theory, the free correlation function in where hn is the anisotropic hopping coefficient of the the presence of a non-vanishing background field is shown anisotropic massive lattice Green’s function Vm (0) at 14 ′ the source x = 0 shown in Appendix A. The one-loop For ˜b d˜, i.e., at the limit of small E S, the free energy free energy can be expressed by the hopping expansion due to≪ the mean field approximation with up to the two of Eq.(B6) as follows: loop corrections for the up to order ˜b4 is as follows: ∞ ˜ n ˜ n ′1loop 1 hn bηL bηT ψ˜2 ˜b2 1 F = − + , (B7) F ′ = log I (α) η2 + η2 1+ hˆ 2 n ˜ ˜ ˜ 0 ˜ L T ˜ 2 n=0,2,4,... " d ! d ! # 4b − − 2d 2d X ˜3 ˜4 ˆ b 2 2 b 4 4 1 h4 3(η ˙T ) + (η ˙L) 3η +3η 1 + ′ ˜2 ˜2 T L ˜ ˜2 where ˜b E Sd˜. Next, from the one-particle irreducible −6d − 2d ( − 2d 12d ! ≡ h i diagrams and , the free energy of the two loop 1 2 2 ˜5 + 3η ˙T η 4ηLηT η˙T +2ηLη˙LηT + η η¨L +0 b , ˜ T − L corrections is derived⊖ as follows, respectively: ψ ˆ 2 ˆ 4 2 h2 2 γ˜ γ˜ , h4 6˜γ + (24D 12)γ ˜ + (6 24D)˜γ, (B12) 2 2 ≡ − ≡ − − ′ 1 ... F = − 4Q¨( α˜ ) 3 TrG(0) +6 TrG(0) +8Q ( α˜ ) 4! | | | | h i h i ˜ ˜ ¯2 ¯ .... ¯ 2 The result of finding the minimum value ψ0 of ψ accord- 12Gˆ (0)+6 TrGˆ (0) Gˆ (0) + 16 Q ( α˜ )3 Gˆ (0) , (B8) × | | ing to Eq.(B12) is as follows: h i h i 2 2 ˜ 1 1 ... ¯ 2 ¯ 8 1 1/b + 32∆2 F ′ = 4Q¨ +8Q( α˜ ) 18Gˆ(x)Tr Gˆ(x) +36Gˆ(x) ˜ − ⊖ 2 3! | | ψ0 = v , x u 1 64∆ 1 X h h i i u − 3 t ... ¯2 ¯ ... 2 ¯ ˜4 ˜2 ˜ ˆ ˜2 ˜2 ˜ ˆ +36 4Q¨ 8Q Gˆ (x) Gˆ(x)+6 8Q Gˆ(x) b 580d 612d + 51h4 + 68b d 2d + h2 , (B9) − ∆1 , h i ≡ 4096d˜4 ˙ 38˜b4d˜2 60˜b3d˜2 2d˜+ hˆ 3˜b5 34d˜2 60d˜+5hˆ where the dotted accent for Q ( α˜ ) is defined by the mod- − 2 − − 4 | | ∆2 , (B13) ˙ ˆ ≡ ˜4˜ ified derivative i.e., Q (1/ 2˜α )dQ/d α˜ , and G(x)ab is 768d b defined as: ≡ | | | | Therefore, for the mean field approximation with up ′ to the two loop corrections of E S, the critical point 1 µν 2loop tri Gˆ(x1, x2) ψˆa (x1) ψˆb (x2) 1+ g ¯ µ ν Gψ(x1, x2) ′ ′ ab ≡ ≡ ˜ ∇ ∇ ab (E S)c and the triple critical point (E S) are respec- 2d D¯ E tively as follows: Gˆ (x) ψaGˆ(x) ψb, TrGˆ (0) Gˆ(0) , (B10) ≡ ab ≡ ll Xl 1 (E′ )2loop = 1 120D3 +8D2 15˜γ2 + 30˜γ +8 S c (D + γ) − where the trace refers only to the index of the 2 2 ma- n h −1 × 4 trix Gˆ(x) . To calculate the free energy of the two-loop +4D 60˜γ3 + 120˜γ2 58˜γ 45 +2γ 105˜γ3 58˜γ 45 /192(D +˜γ) , (B14) ab − − − − ˆ i o correction, introduce the hopping expansion of G(x)ab, as follows: n ∞ 2 ˜b ˜bη ˜ ˆ ˜2 ˜ ˆ ˜ ˆ ˆ n+1 1156 2d + h2 +64 580d 612d+51h4 34 2d+h2 G(x, y)ab = h(x, y) , ′ tri v − − ˜ ˜ (E S) =ur , (B15) d n=0 d ! u 580d˜2 612d˜+51hˆ ab u 4 X t − h (x, y) 2dδ˜ + gµν ¯ ,, (B11) ≡ xy ∇µ∇ν 1 J. E. Mooij and C. J. P. M. Harmans, New Journal of Novel Magnetism 51, 32-47 (2019). Physics 7, 219 (2005). 14 M. Yoneda, et al., Int. J. Mod. Phys. B Vol. 33, No. 25, 2 J. E. Mooij and Yu. V. Harmans, Nature Physics 2, 169172 1950291 (2019). (2006). 15 M. C. Diamantini, et al, Gauge Topological Nature of the 3 J. E. Mooij, New J. Phys 17, 033006 (2015). Superconductor-Insulator Transition. arXiv:1710.10575. 4 O. V. Astafiev, et al, Nature 484, 355358 (2012). 16 M. C. Diamantini, et al, Bosonic topological insulator in- 5 G.’ t. Hooft, Nucl. Phys B 138, 1-25 (1978). termediate state in the superconductor-insulator transi- 6 M. C. Diamantini, P. Sodano, and C. A. Trugenberger, tion. arXiv:1906.07969. Nucl. Phys. B 474, 641-677 (1996). 17 S. M Hollen, et al, Phys. Rev. B 90, 140506(R) (2014). 7 A. Kr¨amer and S. Doniach, Phys. Rev. Lett 81, 3523-3527 18 Brayden Ware, et al, Phys. Rev. B 92, 195105 (2015). (1998). 19 X. Zhang, and A. Schilling, Phys. Rev. B 97, 214524 8 M. Yoneda, M. Niwa, and M. Motohashi, Physica Scripta (2018). 2012, T151 (2012). 20 D. Das,and S. Doniach, Phys. Rev. B 60, 1261 (1999). 9 V. M. Vinokur, et al, Nature 452, 613-615 (2008). 21 D. Das,and S. Doniach, Phys. Rev. B 64, 134511 (2001). 10 M. Ovadia, et al, Scientific Reports 5, 13503 (2015). 22 P. Phillips,and D. Dalidovich, SCIENCE 302, 243-247 11 M. Mironov, et al, Scientific Reports 8, 4082 (2018). (2003). 12 M. C. Diamantini, et al, Nature Comm. Phys 1, 77 (2018). 23 J. Wu,and P. Phillips, Phys. Rev. B 73, 214507 (2006). 13 M. C. Diamantini, et al, Journal of Superconductivity and 24 P. Phillips, SCIENCE 366, 1450-1451 (2019). 15 25 H. Kleinert, Gauge fields in condensed matter. Vol. 1: Su- frozen lattice superconductor (Cambridge University Press, perflow and vortex lines. Disorder fields, phase transitions 2007). (World Scientific, Singapore, 1989). 30 T. Neuhaus, A. Rajantie, and K. Rummukainen, Nuclear 26 J. Villain, J. de Phys 36, 581 (1975). Physics B 119, 915-917 (2003). 27 W. Janke and H. Kleinert, Nuclear Physics B 270, 135-153 31 T. Neuhaus, A. Rajantie, and K. Rummukainen, Phys. (1986). Rev. B 67, 014525 (2003) 28 H. Kleinert, Int. J. Engng. Sci. 23, 927-935 (1985). 29 I. Herbut, A Modern Approach to Critical Phenomena. :7.1