Superconductor-insulator transition and topological nature of the Bose metal

M. C. Diamantini,1 C. A. Trugenberger,2 I. Lukyanchuk,3 and V.M. Vinokur4, 5 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 3University of Picardie, Laboratory of Condensed Matter Physics, Amiens, 80039, France 4Materials Science Division, Argonne National Laboratory, 9700 S. Cass Ave, Argonne, IL 60439, USA 5Computation Institute, University of Chicago, 5735 S. Ellis Avenue, Chicago, IL 60637, USA It has long been believed that at absolute zero can form only one quantum coherent state, a super- conductor. Yet, several two dimensional superconducting systems were found to harbor the superinsulating state with infinite resistance, a mirror image of , and a metallic state often referred to as Bose metal, characterized by finite longitudinal and vanishing Hall resistances. The nature of these novel and mysterious quantum coherent states is the subject of intense study. Here, we propose a topological gauge description of the superconductor-insulator transition (SIT) that enables us to identify the underlying mechanism of superinsula- tion as Polyakov’s linear confinement of Cooper pairs via instantons. We find a criterion defining conditions for either a direct SIT or for the SIT via the intermediate Bose metal and demonstrate that this Bose metal phase is a Mott topological insulator in which the -vortex liquid is frozen by Aharonov-Bohm interactions.

The dual quantum Aharonov–Bohm and Aharonov–Casher the first-order Lagrangian [28] L = (ϕN˙ − Nϕ˙)/2 − H(N, ϕ), effects [1,2] and a paradigm of topological phase transi- having as canonically conjugate pair of variables an ampli- tions by Berezinskii-Kosterlitz-Thouless (BKT) [3–6] brought tude N and the 2π-periodic phase ϕ. Depending on the pa- topology to prominence in condensed matter physics. The rameters of the Hamiltonian H(N, ϕ), one finds three possi- superconductor-insulator transition (SIT) in strongly disor- ble outcomes of the uncertainty principle following from the dered superconducting films and Josephson junction arrays commutation relation, [N, ϕ] = i~: (i) sharp ϕ = x and plane (JJA) [7–18], gripping both in a whole, hosts fascinating quan- waves ∼ e−ixN for H ≡ H(ϕ); (ii) sharp N = y and plane tum phases emerging from the intertwined charges and vor- waves ∼ e−iyϕ for H ≡ H(N); and (iii) eigenstates with fixed tices and their topological interactions. The nature of the SIT quantized uncertainties ∆N and ∆ϕ in the general case. Go- critical region and the phases it harbors still remain mysterious ing over to an infinite number of degrees of freedom, en- and are the subject of intense research and debate. Here we coded in a field Ψ = N exp(iϕ), one finds that all three re- develop a long-distance topological gauge theory of the SIT alizations of the uncertainty principle appear manifestly as and find, in its critical vicinity, three competing quantum or- distinct quantum states around the quantum tri-critical point ders: the superconductor, the superinsulator [14–17], and the of the SIT: a superconductor, comprising N-plane waves and quantum Bose metal (QM) [19, 20], which, as we will show, fixed ϕ (Cooper pair condensate); a superinsulator, ϕ-plane is a Mott topological insulator (TI) [21–23]. We find that, if waves with sharp N (vortex condensate) [14–17]; and a phase quantum fluctuations are weak, the SIT occurs as a first-order harboring both frozen charge and vortex fluctuations (neither direct superconductor-superinsulator transition with granular Cooper pairs nor vortices condense), which is the topological superconducting-insulating texture forming around the phase insulator [21–23]. Accordingly, the SIT is a quantum phase boundary. Strong quantum fluctuations, instead drive the SIT transition between the generic dual superconducting and su- via the intermediate TI phase. The three quantum phases near perinsulating states, which can be either a direct transition or the SIT epitomize the possible mechanisms for a gauge field and indirect one, going through the intermediate topological mass without Higgs fields, with the TI and the superconduc- insulator state. tor realizing topologically massive gauge models [24] and the We consider a lateral JJA, which sets the perfect stage for superinsulator implementing [25] Polyakov’s instanton-driven a gauge theory and adequately models superconducting disor- linear confinement with neutral mesons as excitations [26]. dered films [12, 29]. The properties of the JJA are controlled

Our findings pave the way to a direct probing of gauge topo- by the competition between the Josephson coupling energy EJ 2 logical models by desktop experiments on easily accessible and the charging energy EC = (2e )/2C of a single Josephson arXiv:1710.10575v2 [cond-mat.supr-con] 13 Jun 2018 superconducting systems. junction, C being the junction capacitance [7]. At EJ > EC, su- perconducting correlations win and the system is a supercon- In the framework of the gauge theory, the superconductor- ductor. At E < E , the Coulomb blockade turns the system insulator transition (SIT) [7–13] appears as a material realiza- J C insulating. tion of the field-theoretical S-duality, which, in its simplest form, expresses the invariance of Maxwell’s equations under The Hamiltonian for a lateral JJA, modeled as a 2D grid of the interchange of electric and magnetic fields in the presence superconducting granules coupled by weak links is of magnetic monopoles [27]. To gain an insight into the role 1 X X of S-duality in the physics of the SIT, let us consider its sim- H = V(C − C`2∆)V + E (1 − cos(2`∇ φ)), (1) 2 0 J l plest quantum mechanical example, a system endowed with x x,l 2

Superconductor Superinsulator a M b M 1 1 e Q Q e -1 0 1 -1 0 1 Topological insulator -1 -1 Topological insula

Superinsulator/Superconductor Topological insulator coexistence c M d M Superinsulator Superinsulator Superconductor 1 1 Superconductor Q Q -1 0 1 -1 0 1 -1 -1

FIG. 1. The configurations that minimize the string gas free energy are determined by the integer quantum numbers that fall in the interior of an ellipse with semi-axes given by the model parameters. (a) Superconductor: strings with electric quantum numbers condense. (b) Superinsula- tor: Strings with magnetic quantum numbers condense. (c) Coexistence of long electric and magnetic strings: this is an unstable configuration c near the first-order direct transition from a superinsulator to a superconductor, the strings with the quantum numbers that minimize their self-energy are the stable configurations. (d) Topological insulator/quantum metal: all strings are suppressed by their high self-energy. The resulting phase diagram is shown in panel (e).

where the sum is taken over all the points x of the grid, ` is Thus, the first three terms express the doubled topologically the lattice spacing, l is the vector connecting adjacent super- massive gauge theory [24] for a 2D gauge invariant massive conducting granules, V is the electric potential in the granule photon with no Higgs field. Note that, contrary to the pure x, ∇lφ is the superconducting phase difference between ad- Chern-Simons term used to model topological states in strong jacent granules, C0 is the granule capacitance to the ground, magnetic fields [30], this mixed, or doubled Chern-Simons and we assume C  C0. To construct the gauge theory of term does not violate parity and time-reversal (in Minkowski the SIT, we follow [12] and express the JJA partition func- space) since it involves a vector and a pseudovector gauge R tion Z = DVDϕ exp(−H/T) via an Euclidean action of a field [14]. The last two terms, comprising integer-valued fields coupled Coulomb gas of charges and vortices. Then we intro- Qµ and Mµ, describe the topological excitations – electric duce two gauge fields aµ and bµ which mediate these Coulomb and magnetic strings [26] – arising from the compactness of interactions [14] and, after a few transformations obtain (see the U(1) gauge group and representing (Cooper pair) charge Supplementary Material (SM)) the Euclidean action of a topo- and vortex excitations, respectively. Strings can be closed, logical gauge field theory: in which case they are the “world-lines” of charge-anticharge X  `3 `3 `3 and vortex-antivortex quantum fluctuations over the ground S = fµ fµ + i aµkµνbν + gµgµ + state, or else, infinitely long, representing then the “world- 8π2E π 16E x J c lines” of point charges and vortices. Infinitely long strings √    can also end on electric or magnetic (monopole) instantons, i` 2 aµQµ + bµ Mµ , (2) describing tunneling events [26]. where the sum runs over the Euclidean 3D discrete lattice, The gauge theory defined by the action (2) enables a com- Greek indices denote the coordinate axes and repeated indices prehensive description of the interplay of the distinct orders mean summation. We use natural units, c = 1, ~ = 1, but emerging in the critical vicinity of the SIT. Our first step is we will restore the physical units when necessary. The first to find the free energy associated with the charge and vor- and the third terms in this action describe Josephson cou- tex topological excitations. To that end, we integrate out the pling and Coulomb energies in the JJA expressed through the gauge fields in (2) and obtain (see SM) the free energy of gauge fields, fµ = kµνbν, gµ = kµνaν, and kµν = S µµανdα is a topological string of length L = N` carrying electric and the lattice Chern-Simons operator, formulated in terms of the magnetic quantum numbers Q and M, respectively, assigned lattice derivative dα and shift operator S µ f (x) = f (x + `µˆ) to every bond belonging to the string in a form guaranteeing gauge invariance (see SM). The sec- ! ond term, the so called mixed Chern-Simons term, describes 1 1 F = Q2 + gM2 − µηN , (3) the Aharonov-Bohm coupling between charges and vortices. g η 3

Q f f

FIG. 2. Phase transitions induced by an external magnetic field. The frustration parameter 0 6 f < 1 displaces the original ellipse along the magnetic axis: (a) Direct transition from a superconductor to a superinsulator for η < 1. (b) Transition from a superconductor to a topological insulator for η > 1. FIG. 3. Representative plot of η as a function of the parameter m` for µ = 0.25 ln(5). The critical value of m` that separates the domains of where we introduced the dimensionless coupling parameter weak and strong quantum fluctuations is (m`)c = 0.65. At m` > 0.65, p 2 g = π EJ/(2EC), with g=gc=1 corresponding to the SIT, and quantum fluctuations are strong enough so that η > 1 and the SIT at the dimensionless quantity η = πm`G(m`)/µ, reflecting the T = 0 occurs via the intermediate topological insulator phase. strength of quantum fluctuations, and µ ' ln(5) is the string entropy per bond. Here G(m`) is the diagonal element of the duce the frustration factor f = B/BΦ, where BΦ is the mag- lattice Green function G(x √− y) representing the inverse of the netic field corresponding to one flux quantum Φ = π~c/e 2 2 −∇2 0 operator ` (m ), m = 8EJ EC is the Chern-Simons mass, piercing one plaquette of the JJA. An external magnetic field corresponding to the JJA plasma frequency, which is the in- corresponds to a special case of condensed magnetic strings verse screening length associated with the Aharonov-Bohm with non-integer quantum number and results thus in a shift interactions. of the magnetic quantum number in the free energy and, ac- If F is negative, i.e. if cordingly in the condensation condition (4): M → M+ f . This (1/g)Q2 + gM2 < 1/η , (4) shift along the magnetic axis modifies the condensation condi- tions, see Fig.2. Let us assume g = 1+, where   1, so that the proliferation of loops of an arbitrary size, infinitely long the system is in a superconducting state but close to the SIT. 2 strings and instantons, i.e. the Bose condensation of charges Then for a direct SIT at η < 1 one finds fc = (1/2)(g −1) ≈ . and/or vortices, becomes energetically advantageous. The At η > 1, but still close to the tri-critical point, we set g = η+ condensation condition (4) implies that a particular conden- and find that the superconductor√ transforms into a topological 3/2 1/2 3/2 sate forms if the pair {Q, M} on a square lattice of integer insulator at fc = /η = (g − η) /η . electric and magnetic charges falls within the interior of an We now extend the gauge theory of the SIT and the analysis 1/2 1/2 ellipse with semi-axes rQ = (g/η) and rM = 1/(gη) . of the corresponding emergent phases onto disordered super- Figure1a shows the conditions for the formation of the conducting continuous media. To do so, we transcribe the

Cooper pair condensate comprising quantum fluctuations car- above results obtained in terms of the JJA parameters EC and rying ‘unit’ charges ±2e. In the superconducting state a charge EJ into the language of material characteristics of supercon- current does not induce any voltage since the condensate of ducting films. The role of the tuning parameter driving the electric strings prevents the formation of vortex loops that film across the SIT at zero magnetic field is taken by the re- would cause a finite resistance. In the dual, superinsulating sistance per square, R (or, equivalently, by the dimensionless state, see Fig.1b, M = ±1 and vortex-antivortex quantum conductance g = 4RQ/R). The disorder-driven SIT in films fluctuations form a Bose condensate which blocks the propa- is expected to occur at g = gc = 1 which corresponds to the p 2 gation of Cooper pairs. The conditions shown in Fig.1c, refer condition 2EC/π EJ = 1. A further connection between the to coexisting Cooper pair and vortex condensates, Q, M = ±1. material characteristics of the superconducting films and those 2 Finally, Fig.1d illustrates the situation where none of the con- of JJA is established by the relation [29] λ⊥ = cΦ0/(8π Ic), densates can form, Q = M = 0. This is the intermediate state 2 where λ⊥ = λL /d is the Pearl screening length, λL is the often referred to as a quantum (Bose) metal [19, 20]. Under bulk London penetration depth, d is the film thickness, and the condition of exact duality and charge-hole symmetry, the Ic ≡ (2eEJ/~) is the critical current of a single Josephson junc- Bose metal is characterized by the universal quantum resis- tion. Then, our dimensionless parameter m` acquires the form tance R = h/(4e2) at T = 0 [11]. In real materials the resis- Q 1 1 tance can significantly deviate from the quantum value. m` = , (5) To include the magnetic field-driven SIT in systems that are 8α κ already on the brink of the SIT i.e. have g ≈ 1, let us intro- where α = e2/(~c) is the fine structure constant, ` now plays 4 the role of the characteristic microscopic cutoff length of order (Bose) metal phase. of the superconducting coherence length ξ and thus κ = λ⊥/` To unravel the nature of the intermediate quantum Bose is the (2D) Landau parameter of the superconducting film. Fi- metal, whose mysterious nature is the subject of intense nally we have to generalize (2) to its non-relativistic version scrutiny [32–38], we return to the gauge action (2). When containing material parameters such a magnetic permeabil- neither Cooper pairs nor vortices condense, they experience ity µ and a dielectric constant  with light velocity in the the mutual statistical repulsion [39] due to the topological P √ P medium given by v = 1/ µPP (see SM). We note here that Aharonov-Bohm/Aharonov-Casher effects and get thus frozen JJA correspond formally to the limit µP → ∞ with P = 1 in a “doubled version” of an incompressible quantum fluid. As (see SM). For the superconducting films of interest in this pa- a consequence, the bulk dynamics is completely suppressed per, on the contrary, only the dielectric constant is a relevant and the flow of charges and vortices is supported only by parameter while we can safely set µP = 1. In this case the edge excitations. To see this, one has to calculate the induced µ 3 eff parameter η decomposes into three factors (see SM) charge current jind = (1/` )[δS (Aµ)/δAµ], where Aµ is the electromagnetic potential coupled to the probe electric charge. 1 1 πv2G˜(ml) η = . (6) In the phase where topological excitations do not condense, α κ 8µ the kinetic terms in the action (2) become irrelevant and the The first, universal factor indicates that the parameter η de- whole action is dominated by the topological doubled Chern- scribes the strength of quantum fluctuations: an intermediate Simons term [40], leading to an effective action in which the quantum metal/topological insulator phase occurs when these topological Chern-Simons mass suppresses all bulk conduc- are strong enough, i.e. η > 1. The second and third factors tances: embody all relevant material parameters, namely the Landau 4 2 eff   X ` e gηµ  2  parameter of the superconducting substrate and the dielectric S Aµ = Aµ −δµν∇ + dµdˆν Aν . (7) π2 constant. This latter appears only in the third term, where all x non-relativistic effects are bundled. In disordered supercon- This is the effective action of a bulk insulator. The functional ductors, where the coherence length ξ exceeds the mean free derivative with respect to the field Aµ yields the current in- 2 2 2 ~ path, λL = [4πnse /(mc )](τ∆/ ), where ns is the density of duced by external electromagnetic fields Fµν = ∂µAν − ∂νAµ superconducting electrons, which at T = 0 is equal to the to- as tal density n and τ is the transport scattering time. e 2 Accordingly, m` ' (α/4)(`dn2/3)(τ∆/~). If we associate the µ 2e gηµ µν e j = ` ∂νF , (8) ultraviolet cutoff ` with the minimal superconducting scale ξ, ind π2 2/3 the dependence m` ∼ ne reflects correctly the trend of cross- showing that only variations of external electric and magnetic ing over from the direct SIT to SMIT upon an increase of the field, but not homogenous field themselves, can cause an in- electron density in films. duced current, which is a typical electromagnetic response of Taking µ = ln(5) and an example for which r = an insulator. Therefore, the bulk conductances, both the lon- v2G˜(ml)/G(ml) = 4, see Fig.3, the magnitude of the parame- gitudinal and the Hall ones, vanish. ter (m`)c where all three phases meet, i.e. its magnitude cor- This is not the whole story though. The Chern-Simons responding to the tri-critical point, is (m`)c ≈ 0.65. However, effective action is not invariant under gauge transformations when describing real superconducting films one has to take ai = ∂iλ and bi = ∂iχ at the edges. Two chiral bosons [41] into account that ε and thus also v and the above parameter r λ = ξ + η and χ = ξ − η have to be introduced to restore are not fixed but, rather depend on the proximity of the system the full gauge invariance (see SM), exactly as it is done in the to the SIT and that ε diverges as g → 1 [17]. Since η decreases quantum Hall effect framework [42]. The full gauge invari- with increasing ε, the scenario of the SIT, direct vs. via the in- ance is restored by adding the edge action (we use continuum termediate quantum metal phase, depends on the magnitude of notation for simplicity’s sake) η(ε ) where ε is the maximal dielectric constant achieved max max Z Z  √  at the value of g = 1, where the correlation length associated 1 2 2  2  S edge = d x (∂0ξ∂sξ − ∂0η∂sη)+2e d x A0  ∂sχ , with the SIT compares to the lateral dimension of the film. π  2π  4 Taking an estimate of εmax ' 10 as a characteristic value for (9) the NbTiN film, where the divergent ε was observed on ap- where we have included the√ electromagnetic coupling of the proach to the magnetic field-driven SIT [31], and using a cor- edge charge density ρ = ( 2/2π)∂sχ in the As = 0 gauge. responding estimate of m` ≈ 0.26 for this material one obtains As in the case of the quantum Hall effect, the non-universal η < 1 and hence one expects this NbTiN film to exhibit a di- dynamics of the edge modes is generated by boundary effects rect SIT, as indeed is observed experimentally. Analogously, [42], which result in the Hamiltonian one can observe that near the SIT η < 1 for TiN, for which 1 Z h i m` ≈ 2.63, which thus also follows the direct SIT scenario, H = ds −v (∂ ξ)2 − v (∂ η)2 , (10) π s s as observed. Instead, for NbSi, one does not expect a ε di- vergence and, correspondingly, the estimate m` ≈ 21.6 gives where v is the velocity of propagation of the edge modes along η > 1 confirming a transition via an intermediate quantum the boundary. Upon adding this term, the total edge action 5

c a b

NM NM TI TI TI TI TI Temperature Temperature SI SC SI SC SC TI VL SI SC Temperature

gc gc SIT tuning parameter, g SIT tuning parameter, g Magnetic field

p 2 FIG. 4. Phase diagram of the vicinity of the SIT. Tuning the parameter g = π EJ/2EC, one drives the system across the SIT. The quantity η characterizes the strength of quantum fluctuations. Panels (a) and (b) show the SIT in zero magnetic field. (a) Weak quantum fluctuations, η < 1, where the system experiences the direct SIT at g = gc and T = 0. (b) If quantum fluctuations are strong, η > 1, an intermediate quantum state, Bose metal, which is demonstrated to be a topological insulator, opens up between the superconductor and the superinsulator at T = 0. The solid lines separating phases are deconfinement lines coinciding with the BKT transitions lines. The dashed lines are thawing lines where the frozen Chern-Simons topological insulator melts. At B = 0 thawing lines coincide with the superconducting, Tc and insulating, TI transition temperatures, respectively. At finite magnetic field Chern-Simons phase melts into either conventional insulator or vortex liquid (weak metal). (c) Sketch of the phase diagram in the T–B coordinates at the superconducting side and η > 1. becomes come strong enough so that the frozen state of intertwined, 1 Z non-condensed vortices and Cooper pairs would thaw and the S = d2 x (∂ − v∂ ) ξ∂ ξ − (∂ + v∂ ) η∂ η edge π 0 s s 0 s s TI transforms into an ordinary bulk insulator or into a vortex Z  √  liquid (bad metal). Taking into account non-relativistic effects 2  2  +2e d x A0  ∂sχ . (11) due to material parameters, the energy gap associated with the  2π  Chern-Simons√ mass in the critical vicinity of the SIT is (see ν/2 The equation of motion generated by this action is SM) m/  ∝ (g − gc) , where ν is the scaling exponent de- scribing the divergence of the dielectric constant of the film 2e 2e v∂sρ = E = ∂sA0 , (12) on approach to the SIT [17]. Then, by the the duality prin- 2π 2π ciple [10, 11], one would expect the Chern-Simons melting ν/2 i.e. the charge conduction with the quantum resistance RQ = temperature to scale as Tth ∝ |g − gc| , as shown in Fig.4 by h/(4e2). Thus, while all bulk conductances are suppressed dashed lines. by the large Chern-Simons mass, there remains the ballis- A qualitative description of the transformation of the su- tic edge conductance with the quantum resistance, which is perconducting/superinsulating phases into the TI and further characteristic of a topological insulator [21, 40], in this case a into a bulk insulator/vortex liquid is achieved by including Mott topological insulator. The edge modes in this topologi- the magnetic field driving the transition into our consider- cal insulator carry the charge 2e since the elemental carriers ation. To that end, let us first introduce the dimensionless mediating the conductivity are Cooper pairs. This is remi- temperature τ = T`/(~vc), with vc the light velocity in the niscent of recent observations of doubly charged egde exci- material. As derived in a previous section, the critical√ field 3/2 tations in the integer quantum Hall effect regime [43]. Then, for the superconductor-to-TI transition is fc = /η = a straightforward calculation (see SM) yields the TI sheet re- (g − η)1/2/η3/2. The important point is then that, at finite tem- 2 sistance R = h/(4e ) ≡ RQ, which corresponds to perfect peratures, the parameter η scales as η(T) = η0S (T), with a duality, where exactly one fluxon for each Cooper pair tra- scale factor S (T) > 1 and η0 = η(T = 0), as derived in verses the system. This coincides with the result by M. P.A. [25]. Using this scaling we obtain the equation that defines

Fisher [10] obtained by an elegant qualitative consideration. the SC↔TI transition temperature τTI as a function of the ap- In experiments, the resistance can significantly differ from RQ plied magnetic field, signaling deviations from strict duality [37]. Remarkably, our p prediction that in the Bose metal the Hall resistance should η0(1 − S (τTI)) + 0 f (τTI) = , (13) disappear is in full accord with the recent observation by [38]. 3/2 3/2 η0 S (τTI) We further expect that, at the temperature T = Tth associ- ated with the topological mass gap, thermal fluctuations be- where 0 = g − η0. The critical temperature for the SC↔TI 6

ξ was conjectured [16, 44] that it is this emergent granularity, SC SC see Fig.5, that serves as a material platform for the superinsu- SC lating state. Yet, while being generally accepted, the concept SC SC of emergent granularity is not thoroughly justified. Our gauge approach to the SIT enables us to put this concept on a firm SC SC field-theoretical foundation. SC To see this, note that there is a finite width strip embrac- ing the line of the direct SIT at g = g = 1 and η < 1, SC SC SC c where Cooper pairs and vortex condensates coexist. Calcu- SC SC lating from (2) at g > 1 and η < 1 the energy E of a supercon- ducting droplet of perimeter L immersed into a superinsulat- SC SC SC ing matrix of area A, we find (see SM) ! SC 1 1 SC E ∝ A − 1 − L2 + σL , (16) 4π g2 SC SC where the boundary contribution arises due to fluctuation- induced charge and vortex excitations within the opposite con- densates, respectively, and σ is a numerical coefficient. The FIG. 5. Self-induced electronic ‘granular’ structure in the crit- ∝ 2 ical vicinity of the SIT. At the first order direct superconductor- energy has a global minimum Eglobal A/g when the su- superinsulator transition (η < 1) the phase separation occurs and perconductor droplet fills the whole area A, but also a local the system breaks into the superconducting droplets immersed into minimum Elocal ∝ A at L = 0. These two minima are sepa- an insulating matrix and connected by weak Josephson links. The rated by a maximum at L ∝ g2/(g2 − 1). If close enough to the characteristic size of a droplet is & ξ. SIT, this maximum spreads over a scale exceeding the radius of the droplet. Hence, it becomes energetically advantageous to fragment a superconducting droplet into smaller ones. The transition at f = 0 is then defined by the equation fragmentation stops at the minimal dimension of order ξ set- ting the scale of the self-induced granularity. S τ  /η . ( TI) = 1 + 0 0 (14) There is a deep connection between the structure of the superinsulating state and field theory models of confine- The dimensionless thawing temperature, instead, is deter- ment [25]. To illustrate this, we consider again the effective mined by the topological mass splitting caused by the external action for electromagnetic gauge fields. In the superinsulating magnetic field (see SM), phase it is given by (see SM) m`   gµη X   τ = (1 − f /2) , (15) S eff A = 1 − cos 2e`2F  , (17) th ~ µ π2 µ vc 2 x,µ and the intermediate phase appears at zero magnetic field if where Fµ = kµνAν is the dual field strength. This action co- m`/(~vc) > τTI. The sketch of the resulting SC↔TI transi- incides with Polyakov’s compact quantum electrodynamics tion phase diagram in T-B coordinates is shown in Fig.4c (QED) model [26], in which a non-perturbative photon mass and reproduces the corresponding experimental phase dia- and linear confinement appear due to instanton effects. One gram of [38]. can thus view the dipole Cooper pair-anti-Cooper pair neu- Being an exemplary laboratory for theoretical and exper- tral excitations in the superinsulator as neutral “U(1) meson” imental study of the SIT and emerging phases in its critical excitations [25]. The charge BKT transition of the SIT cor- vicinity, Josephson junction arrays offer a perfect model for responds then to the deconfinement transition of the (2+1)- strongly disordered superconducting films [12, 29]. More- dimensional compact gauge theory [47]. This analogy has far over, a striking quantitative agreement between the experi- reaching implications. In particular, depending on their spin mentally measured characteristics of the superinsulating state quantum number, these neutral bound states of Cooper pairs in disordered TiN films [44] and its description in terms of a could decay into, or mix with photons via the one-loop vac- regular JJA [18] suggests that the close connection between uum polarization diagram (see SM). In this case, applying an these two systems goes well beyond the similarity of coarse- external electric field along the film would cause this to emit grained JJA and disordered film. The idea of emergent elec- light and/or flicker. This effect calls for further investigation. tronic granularity, i.e. that even in structurally homogeneous films the electronic texture in the critical vicinity of the SIT is an array of superconducting puddles connected by weak links Acknowledgements and immersed in an insulating matrix was pioneered in [45]. Over two past decades, this hypothesis evolved into a paradig- We are delighted to thank M. Vasin and Ya. Kopelevich matic attribute of the SIT, see [46] and references therein. It for illuminating discussions. M. C. D. thanks CERN, where 7 she completed this work, for kind hospitality. The work of [20] Das, D. & Doniach, S., Bose metal: gauge field fluctuations and V.M. V. was supported by the U.S. Department of Energy, Of- scaling for field-tuned quantum phase transitions. Phys. Rev. fice of Science, Materials Sciences and Engineering Division. B 64 134511 (2001). [21] Kane, C. L. & Mele, E. J., Z2 Topological Order and the Quan- tum Spin Hall Effect. Phys. Rev. B 95 146802 (2005). [22] Fu, L., Kane,C. L. & Mele, E. J., Topological insulators in three dimensions. Phys. Rev. Lett. 98 106803 (2007). [23] Moore, J. & Balents, L., Topological invariants of time-reversal- [1] Aharonov, Y., Bohm, D., Significance of Electromagnetic Po- invariant band structures. Phys. Rev. B 75 121306(R) (2007). tentials in the Quantum Theory. Phys. Rev. 115, 485–491 [24] Deser, S., Jackiw, R. & Templeton,§., Three-Dimensional Mas- (1959). sive Gauge Theories. Phys. Rev. Lett. 48 975 (1982). [2] Aharonov, Y., Casher, A., Topological Quantum Effects for [25] Diamantini, M. C., Trugenberger, C. A. & Vinokur, V.M., Con- Neutral Particles. Phys. Rev. Lett. 53, 319-321 (1984). finement and Asymptotic Freedom with Cooper pairs, under [3] Berezinskii, V.L., Destruction of long-range order in one- consideration for publication. dimensional and two-dimensional systems having a continuous [26] Polyakov, A. M., Gauge Fields and Strings, Harwood Aca- Sov. Phys.–JETP 32 symmetry group I. Classical systems. 493– demic Publisher, Chur (Switzerland) (1987). 500 (1970). [27] Goddard, P. & Olive, D. I., Magnetic monopoles in gauge field [4] Berezinskii, V.L., Destruction of Long-range Order in One- theories. Rep. Prog. Phys. 41 1357 (1978). dimensional and Two-dimensional Systems Possessing a Con- [28] Dunne,G. Jackiw,R., & Trugenberger, C. A., “Topological” Zh. Eksp. tinuous Symmetry Group. II. Quantum Systems. (Chern-Simons) quantum mechanics. Phys. Rev. D41 661-666 Theor. Fiz 61 Sov. Phys.– JETP 34 . , 1144 (1971). ( , , 610-616 (1989). (1971)). [29] Tinkham, M., Introduction to Superconductivity. McGraw-Hill, [5] Kosterlitz, J. M., Thouless, D. J., Long range order and metasta- Inc. (1996). bility in two dimensional solids and superfluids. (Application of [30] Wen, X.-G., & Zee,A., Classification of Abelian Quantum Hall Journal of Physics C: Solid State Physics, dislocation theory). States and Matrix Formulation of Topological Fluids, Phys. 5 , L124 (1972). Rev.B46 2290-2301 (1992). [6] Kosterlitz, J. M., Thouless, D. J., Ordering, metastability and [31] Mironov, A. Yu., Silevitch, D. M., Proslier, T., Postolova, S. V., phase transitions in two-dimensioal systems. J. Phys. C: Solid Burdastyh, M. V. Gutakovskii, A. K., Rosenbaum, T. F., State Phys. 6 1181–1203 (1973). Vinokur, V.M. & Baturina, T. I. Charge Berezinskii- [7] Efetov, K. B., Phase transition in granulated superconductors. Kosterlitz-Thouless transition in superconducting NbTiN Sov. Phys. JETP 51, 1015 – 1022 (1980). films. arXiv:1707.09679 (2017). [8] Haviland, D., Liu, Y., & Goldman, A., Onset of superconduc- [32] Dalidovich, D., & Phillips, P., Phase Glass is a Bose Metal: A tivity in the two-dimensional limit. Phys. Rev. Lett. 62, 2180 – New Conducting State in Two Dimensions. Phys. Rev. Lett. 89, 2183 (1989). 27001 (2002). [9] A. Hebard, A. & Paalanen, M. A., Magnetic-field-tuned [33] Phillips, P. & Dalidovich, D., The Elusive Bose Metal. Science superconductor-insulator transition in two-dimensional films. 302 243 – 247 (2003). Phys. Rev. Lett. 65, 927 – 930 (1990). [34] Breznay, N. P., Steiner, M. A., Kivelson, S. A., & Kapitulnik, A., [10] Fisher, M. P.A., Grinstein, G. & Girvin, S. M., Presence of Self-duality and a Hall-insulator phase near the superconductor- quantum diffusion in two dimensions: Universal resistance at to-insulator transition in indium-oxide films. Proc. Natl. Acad. the superconductor-insulator transition. Phys. Rev. Lett. 64, 587 Sci. 113, 280–285 (2016). – 590 (1990). [35] May-Mann, J. & Phillips, P.W., Vanishing hall conductance [11] Fisher, M. P.A., Quantum Phase Transitions in Disordered in the phase glass bose metal at zero temperature. arXiv Phys. Rev. Lett. 65 Two-Dimensional Superconductors. , 923 – 1708.09016, (2017). 926 (1990). [36] Ying, T., Dalmonte, M., Angelone, A., Mezzacapo, F., Zoller, P., [12] Fazio, R. & Schon,¨ G., Charge and Vortex Dynamics in Arrays & Pupillo, G., Cluster bose metals, arXiv 1606.04267, (2017). of Tunnel Junctions. Physical Review B 43, 5307 – 5320 (1991). [37] Kapitulnik, A., Kivelson, S. A., & Spivak, B., Anomalous met- Int. J. [13] Goldman, A., M., Superconductor-Insulator Transitions. als – failed superconductors. arXiv:1712.07215v1 [cond- Mod. Phys. 24 B , 4081 – 4101 (2010). mat.supr-con] (2017) [14] Diamantini, M. C., Sodano,P. & Trugenberger, C. A., Gauge [38] Breznay, N. P. & Kapitulnik, A., Particle-hole symmetry re- theories of Josephson junction arrays. Nuclear Physics B 474, veals failed superconductivity in the metallic phase of two- 641 – 677 (1996). dimensional superconducting films. Sci. Adv. 3, e1700612 ¨ [15] A. Kramer, A. & Doniach, S., Superinsulator phase of two- (2017). Phys. Rev. Lett. 81 dimensional superconductors. , 3523 – 3527 [39] Wilczek, F., Disassembling Anyons. Phys. Rev. Lett. 69 132- (1998). 135 (1992). [16] Vinokur, V.M. et al., Superinsulator and quantum synchroniza- [40] Cho, G. Y., & Moore,J. E., Topological BF Field Theory De- Nature 452 tion. , 613 – 615 (2008). scription of Topological Insulators. Ann. Phys. 326, 1515-1535 [17] Baturina, T. I. & Vinokur, V.M., Superinsulator-superconductor (2011). duality in two dimensions. Ann. Phys. 331 236 – 257 (2013). [41] Floreanini, R. & Jackiw, R., Self-dual fields as charge-density [18] Fistul, M. V., Vinokur, V.M. & Baturina, T. I., Collective solitons. Phys. Rev. Lett. 59 1873 – 1876 (1987). Cooper-Pair Transport in the Insulating State of Josephson- [42] Wen, X.-G., Theory of the Edge States in Fractional Quantum Phys. Rev. Lett 100 Junction Arrays. . , 086805 (2008). Hall Effects. Int. J. Mod. Phys. B6 1711-1762 (1992). T [19] Das, D. & Doniach, S., Existence of a Bose metal at = 0. [43] Choi, H. K., Sivan, I., Rosenblatt, A. Heiblum, M., Umansky, V. Phys. Rev. B 60, 1261 – 1275 (1999) & Mahalu, D., Robust Electron Pairing in the Integer Quan- tum Hall Effect Regime. Nature Communications 6:7435 1-7 8

(2015). (1994). [44] Baturina, T. I., Mironov, A., Yu., Vinokur, V.M., Bak- [46] Ganguly, R., Roy, I., Banerjee, A., Singh, H., Ghosa, A. & Ray- lanov, M. R. & Strunk, C., Localized superconductivity in the chaudhuri, P., Magnetic field induced emergent inhomogeneity quantum-critical region of the disorder-driven superconductor- in a superconducting film with weak and homogeneous disor- insulator transition in TiN thin films. Phys. Rev. Lett. 99, der. Phys. Rev. B96, 054509 (2017). 257003 (2007). [47] Svetitsky, B. & Yaffe, L. G., Critical behavior at finite tempera- [45] Kowal, D. & Ovadyahu, Z. Disorder induced granularity in an ture confinement transitions. Nucl. Phys. B210 423-447 (1982). amorphous superconductor. Solid State Commun. 90, 783–786