Topological Gauge Theory of Josephson Junction Arrays: The Discovery of Superinsulation

Preprint

Review

Topological Gauge Theory of Josephson Junction Arrays: The Discovery of Superinsulation

Altmetrics

Downloads

Views

Comments

A peer-reviewed article of this preprint also exists.

Maria Cristina Diamantini^*

This version is not peer-reviewed

Submitted:

25 September 2023

Posted:

26 September 2023

You are already at the latest version

Alerts

Abstract

We review the topological gauge theory description of Josephson junction arrays (JJA), fabricated systems which exhibit the superconductor to insulator transition (SIT). This description revealed the topological nature of the phases appearing at the SIT and lead to the discovery of a new state of matter, the superinsulator, characterized by an infinite resistance, even at finite temperature, due to linear confinement of electric charges. This discovery is particularly relevant for the physics of superconducting films with emergent granularity, which are modeled by JJA, and share the same phase diagram.

Keywords:

Subject: Physical Sciences - Condensed Matter Physics

1. Introduction

Josephson junction arrays (JJA) [1] are fabricated materials which exhibit the superconductor to insulator transition (SIT) [2,3], a paradigmatic example of a quantum phase transition. Planar quadratic JJA are square lattices formed by superconducting island (typical size O(

μ

m)), fabricated on a substrate, with lattice spacing l (O(100nm)), Josephson coupling

E_{J}

and capacitances C between nearest neighbour islands. Each islands is also characterized by the phase of the local order parameter, and by a ground capacitance

C_{0}

. When

C ≫ C_{0}

, JJA have only two relevant energy scales,

E_{J}

, the energy scale associated with the tunneling of Cooper pairs between the islands, and

E_{C} = e^{2} / 2 C

, the charging energy, with e the electron charge. The parameters

E_{J}

and

E_{C}

can be traded for a massive parameter

ω_{P} = \sqrt{2 E_{C} E_{J}}

, which represents the plasma frequency of the array, and a dimensionless parameter

g = \sqrt{\frac{π^{2} E_{J}}{2 E_{C}}}

, which, as we will see, plays a crucial role in the phase structure of the theory. In experimental realizations of JJA, the transition between the different phases is achieved by varying

E_{J}

and taking

E_{C}

essentially fixed. Varying

E_{C}

is much more difficult.

In the classical limit,

E_{J} ≫ E_{C}

, global phase coherence is realized and the array behaves as a superconductor. Above a critical temperature the arrays undergoes a Bere- zinskii-Kosterlitz-Thouless (BKT) transition [4] and superconductivity is destroyed. When, instead,

E_{c} ≫ E_{J}

, the tunneling between the islands is suppressed and the arrays shows an insulating behaviour. As it has been shown in [5,6], however, this is not an insutaltors, but a new state of matter called superinsulator.

Superinsulators are emergent condensed matter states dual to superconductors, and exhibit an infinite resistance at finite temperatures. They were first theoretically predicted in [5], while the final form of the theory describing this new state of matter was established in [7]. Superinsulators have, than, been experimentally found in [8,9,10,11,12,13,14] in thin superconducting films which, close to the SIT, have an emergent granularity and behave as a self organized JJA.

The infinite resistance which characterize superinsulators is due to linear charge confinement [7] of both Cooper pairs and electrons in a magnetic monopole plasma, which squeezes electric field lines into electric flux tubes connecting charge-anticharge pairs, in analogy with the Meissner effects in superconductors [13]. Superinsulators realize, thus, an Abelian version of the dual superconductivity mechanism advocated by ‘t Hooft to explain quark confinement [15]. In this picture chromo-electric strings will form with quarks at their ends and, when quarks are pulled apart, it is energetically favourable to pull out of the vacuum additional quark-antiquark pairs and to form several short strings. Free colour charge can never be observed on distances larger than

1 / Λ_{Q C D}

and quarks are, thus, confined.

In [5] we derived a topological gauge theory description of JJA. For planar JJA the relevant term, dominant at large distances is the mixed Chern-Simons term [16], which is nothing else than the (2+1)-dimensional version of the BF term [17], which is the relevant term at long distances in higher dimensions. This description has made it possible to derive the quantum phase diagram of JJA and to understand the nature of the various phases. This result is particularly important since, as we pointed out before, JJA are a model for thin superconducting films near the SIT, and the gauge theory which describes JJA is also an effective theory for these materials at the SIT.

We have identified three possible phases: when a parameter

η

, which depends on the details of the array, is

η < 1

, there is a direct transition between the superinsulating phase and the superconducting phase. The superconducting phase is a new type of superconductivity, which we called type III superconductivity [18], characterized by the fact that the gap is opened trough a topological mechanism and not trough the Higgs menhanism and is not described by the standard Landau-Ginzurg theory. Superconductivity is destroyed by proliferation of vortices, implying a BKT transition in two dimensions and a Vogel-Fulcher-Tamman transition in three dimensions, and not by the breaking of Cooper pairs [18,19]. This superconductivity model can be relevant for high-

T_{C}

superconductivity [20]. When

η > 1

an intermediate phase opens up between the superinsulating and the superconducting phase [21,22]. This metallic phase in two spatial dimensions, called Bose metal (BM) since charge carriers are Cooper pairs, arising at the SIT and experimentally found in superconducting films [23,24,25,26,27,28,29], has for long time challenged the understanding of electronic fluids [30,31]. In fact it is believed that a metallic phase in two-dimensions cannot exist due to localization. Moreover it was not clear why certain films had this phase and others not. In [21] we have been able to identify this phase as a bosonic topological insulators and to establish, trough the parameter

η

, which depends on the film characteristics, when this phase appears.

In section 2 we will review the topological gauge theory formulation of JJA. We will, in particular, discuss the role of the kinetic term for vortices in the description of JJA. In section 3 we will derive the quantum phase strucure of JJA and, in section 4, we will describe the properties of the new phase of matter, theoretically predicted thanks to the topological gauge theory description of JJA, the superinsulator.

2. Topological gauge theory of JJA

We will consider JJA defined on a square lattice, so we start by defining our lattice notation. We will consider a 3-dimensional Euclidean lattice with sites denoted by

{x}

and directions, indicated by Greek letters, going from 0 to 2. The lattice spacing is

l_{0}

in the 0 direction (time) and l in the 1 and 2 directions (space).The forward and backward finite difference and shift operators are defined as (we denoted the lattice spacing as

l_{μ}

only to distinguish between the 0 and 1 and 2 directions)

\begin{matrix} Δ_{μ} f (x) = f (x + l_{μ} \hat{μ}) - f (x), S_{μ} f (x) = f (x + l_{μ} \hat{μ}), \\ {\hat{Δ}}_{μ} f (x) = f (x) - f (x - l_{μ} \hat{μ}), {\hat{S}}_{μ} f (x) = f (x - l_{μ} \hat{μ}), \end{matrix}

(1)

where

\hat{μ}

denotes a unit vector in direction

μ

. Summation by parts on the lattice interchanges both the two finite differences (with a minus sign) and the two shift operators. Using these definition we introduce the operators:

K_{μ ν} = S_{μ} ϵ_{μ α ν} Δ_{α}, {\hat{K}}_{μ ν} = ϵ_{μ α ν} {\hat{Δ}}_{α} {\hat{S}}_{ν},

(2)

where no summation is implied over the equal indices

μ

and

ν

. These two operators are interchanged by summation by parts with no minus sign. They allow us to define a gauge invariant version of the Chern-Simons operator

ϵ_{μ α ν} \partial_{α}

on the lattice [5]. Their products gives the lattice Maxwell operator:

K_{μ α} {\hat{K}}_{α ν} = {\hat{K}}_{μ α} K_{α ν} = - δ_{μ ν} Δ + Δ_{μ} {\hat{Δ}}_{ν},

(3)

where

Δ = {\hat{Δ}}_{μ} Δ_{μ}

is the 3D Laplace operator.

Our starting point is the Hamiltonian for a planar JJA with nearest neighbours Josephson couplings

E_{J}

, ground capacitances

C_{0}

and nearest neighbours capacitances C [1,32]:

\begin{matrix} H = \sum_{x} \frac{C_{0}}{2} V_{x}^{2} + \sum_{< x y >} \frac{C}{2} {(V_{y} - V_{x})}^{2} + E_{J} (1 - \cos (φ_{y} - φ_{x})), \\ H = \sum_{x} \frac{1}{2} V_{x} (C_{0} - C Δ) V_{x} + \sum_{x, i} E_{J} (1 - \cos (Δ_{i} φ_{x})), \end{matrix}

(4)

where

< x y >

indicates nearest neighbours and

Δ \equiv {\hat{Δ}}_{i} Δ_{i}

is the two-dimensional finite difference Laplacian. In what follows we will use natural units

c = 1

ℏ = 1

ε_{0} = 1

Each island is characterized by an electric potential

V_{x}

and by the phase of the order parameter

φ_{x}

, which are quantum-mechanically conjugated to the charges

E_{x}

on the islands. The charges are quantized in integer multiples of 2e (Cooper pairs),

E_{x} = 2 e q_{x}, q_{x} \in Z

Using the discrete version of Poisson’s equation:

(C_{0} - C Δ) V_{x} = E_{x},

(5)

and introducing the charging energy

E_{C} \equiv e^{2} / 2 C

, we can rewrite the Hamiltonian eq. (4) as:

H = \sum_{x} 4 E_{C} q_{x} \frac{1}{C_{0} / C - Δ} q_{x} + \sum_{x, i} E_{J} (1 - \cos (Δ_{i} φ_{x})) .

(6)

The integer charges

q_{x}

in eq. (6) interact via a two-dimensional Yukawa potential with mass

\sqrt{C_{0} / C} / ℓ

which, in the limit

C ≫ C_{0}

, becomes a 2-dimensional Coulomb potential.

The next step in our derivation of the gauge theory description is the construction of the phase-space path-integral representation [5] of the JJA. To this end we introduce a fictitious temperature

β = 1 / T

and write the partition function of the JJA as:

\begin{matrix} Z = \sum_{{q}} \int_{- π}^{+ π} D φ \exp (- S), \\ S = \int_{0}^{β} d t \sum_{x} i q_{x} {\dot{φ}}_{x} + 4 E_{C} q_{x} \frac{1}{C_{0} / C - Δ} q_{x} + \sum_{x, i} E_{J} (1 - \cos (Δ_{i} φ_{x})), \end{matrix}

(7)

which, in the limit

C_{0} ≫ C

, experimentally accessible, becomes

\begin{matrix} Z = \sum_{{q}} \int_{- π}^{+ π} D φ \exp (- S), \\ S = \int_{0}^{β} d t \sum_{x} i q_{x} {\dot{φ}}_{x} + 4 E_{C} q_{x} \frac{1}{- Δ} q_{x} \\ + \sum_{x, i} E_{J} (1 - \cos (Δ_{i} φ_{x})) . \end{matrix}

(8)

Continuous time must be treated for a discrete time with interval steps

l_{0}

, since our degrees of freedom changes only by integer steps, introducing, thus, also a forward and backward finite time differences

Δ_{0}

and

{\hat{Δ}}_{0}

. Using the Villain representation [33], we can get rid of the cosine at the price of introducing a set of integer link variables

a_{i}

. By introducing real charge currents

j_{i}

, we can express the quadratic term

{(Δ_{i} φ + 2 π a_{i})}^{2}

, which comes from the Villain approximation, as a Gaussian integrals over these variables, thereby obtaining:

\begin{matrix} Z = \sum_{{a_{i}}, {j_{0}}} \int D j_{i} \int_{- π}^{+ π} D φ \exp (- S), \\ S = \sum_{x} i j_{0} Δ_{0} φ + i j_{i} (Δ_{i} φ + 2 π a_{i}) + 4 ℓ_{0} E_{C} j_{0} \frac{1}{- Δ} j_{0} + \frac{1}{2 ℓ_{0} E_{J}} {j_{i}}^{2}, \end{matrix}

(9)

where we renamed

j_{0}

the integer charges

q_{x}

The longitudinal part of

a_{i}

is not physical, and can be reabsorbed in a redefinition of

φ

. As we will see, the transverse degrees of freedom encode the vortex degrees of freedom. In fact, it is important to notice that eq. (9) contains a kinetic term only for the charges, but not for vortices. While for overdamped junction this omission has no consequences, in the general case, this is not correct. As shown in [32], the kinetic term for vortices is generated by integration over charge fluctuations and must, thus, be included at tree level. Moreover in our case the Coulomb interaction is long-range, and the dissipation is reduced [34], making the vortex kinetic term relevant. The kinetic term for the vortices represents tunneling events between adjacent plaquettes of the lattice. These events are the generalization of quantum phase slips, which plays a crucial role in Josephson chain [35], to two dimensions. They can be seen as half-lines of simultaneous phase slips of opposite chirality, which end in the island between the two adjacent plaquettes. Ballistic vortex motion, represented by phase slips corresponding to one vortex tunneling from one plaquette to the other, has actually been exeperimentally observed in [36]. The kinetic term for vortices will involve the time derivative of

a_{i}

{(Δ_{0} a_{i})}^{2}

. Let us start by introducing a fictitious electric field

φ_{i} = K_{i μ} a_{μ}

and a real Lagrange multiplier

a_{0}

. Taking the coefficient of the kinetic term as

π^{2} / 4 ℓ_{0} E_{C}

, we write the vortex kinetic term and the charge Coulomb interaction as:

\begin{matrix} Z = \sum_{{a_{i}}, {j_{0}}} \int D a_{0} D j_{i} \int_{- π}^{+ π} D φ \exp (- S), \\ S = \sum_{x} i j_{0} (Δ_{0} φ + 2 π a_{0}) + i j_{i} (Δ_{i} φ + 2 π a_{i}) + \frac{1}{2 ℓ_{0} E_{J}} {j_{i}}^{2} + \frac{π^{2}}{4 ℓ_{0} E_{C}} φ_{i}^{2}, \end{matrix}

(10)

where now the Coulomb interaction between the charges follows from the Gauss law constraint associated with the Lagrange multiplier

a_{0}

. In this case we have chosen a particular value of the vortex mass for which the JJA is dual under interchanges of charges and vortices and of

π^{2} E_{J} \leftrightarrow 2 E_{C}

, or alternatively,

g \leftrightarrow 1 / g

. This is the self-dual approximation introduced in [5]. It is important to notice, however, that a different vortex mass will, in general, renormalize the value of

E_{C}

, and as consequence change g, the parameter which, as we will see drive the quantum phase transition. A bigger or smaller vortex mobility will, thus, influence the possible phases which are observed.

The next step is to notice that the integration over

φ

in eq. (10) gives the constraint

Δ_{μ} j_{μ} = 0

, the current

j_{μ}

is, thus, conserved and can be represented in terms of a fictitious gauge field

b_{μ}

j_{0} = K_{0 i} b_{i}, j_{i} = K_{i 0} b_{0} + K_{i j} b_{j},

(11)

In eq. (11)

b_{0}

is a real variable, and

b_{i}

and

a_{i}

are integers, which can be made real using Poisson’s formula:

\sum_{n_{μ}} f (n_{μ}) = \sum_{k_{μ}} \int d n_{μ} f (n_{μ}) e^{i 2 π n_{μ} k_{μ}},

(12)

so that all components of the gauge fields

a_{μ}

and

b_{μ}

are real, at the price of introducing integer link variables

Q_{i}

and

M_{i}

\begin{matrix} Z = \sum_{{Q_{i}}} \sum_{{M_{i}}} \int D a_{μ} D b_{μ} \int_{- π}^{+ π} D φ \exp (- S), \\ S = \sum_{x} i 2 π a_{μ} K_{μ ν} b_{ν} + \frac{1}{2 ℓ_{0} E_{J}} {j_{i}}^{2} + \frac{π^{2}}{4 ℓ_{0} E_{C}} ϕ_{i}^{2} + i 2 π a_{i} Q_{i} + i 2 π b_{i} M_{i} \\ + b_{i} ({\hat{K}}_{i 0} Δ_{0} φ + {\hat{K}}_{i j} Δ_{j} φ) + b_{0} {\hat{K}}_{0 i} Δ_{i} φ . \end{matrix}

(13)

The last step of our derivation consist in reabsorbing the quantity

({\hat{K}}_{i 0} Δ_{0} φ + {\hat{K}}_{i j} Δ_{j} φ)

in a redefinition of the integers

M_{i}

, and in introducing a set of integers

M_{0}

trough the definition

{\hat{K}}_{0 i} Δ_{i} φ = 2 π M_{0}

. This step is justified by noticing that

{\hat{K}}_{μ ν} Δ_{ν} φ

are the circulations of the array phases around the plaquettes orthogonal to the direction

μ

in 3D Euclidean space-time and are thus quantized as

2 π integers

. At this point, using again Poisson’s formula, we can express the integral over

φ

as a sum over the integer

M_{0}

, giving the final expression for the topological gauge theory describing JJA:

\begin{matrix} Z = \sum_{{Q_{i}}} \sum_{{M_{μ}}} \int D a_{μ} D b_{μ} \exp (- S), \\ S = \sum_{x} i 2 π a_{μ} K_{μ ν} b_{ν} + \frac{1}{2 ℓ_{0} E_{J}} {j_{i}}^{2} + \frac{π^{2}}{4 ℓ_{0} E_{C}} ϕ_{i}^{2} + i 2 π a_{i} Q_{i} + i 2 π b_{μ} M_{μ} . \end{matrix}

(14)

The first term in the action is the lattice version of the mixed Chern-Simons term, which, being linear in derivatives, is dominant at large distances. The second and third terms are the electric fields for the fictitious gauge fields

a_{μ}

and

b_{μ}

respectively. The effect of the CS term is to give a mass to these gauge fields without the Higgs mechanism [16]. The mass that these fields acquires is nothing else than the plasma frequency of the array:

m_{t o p} = ω_{p} = \sqrt{8 E_{C} E_{J}}

. The two kinetic terms have a coupling constant with dimensions of mass and are thus naively irrelevant. They cannot, however, be neglected since the correct topological limit

m_{t o p} \to \infty

has to be derived from the full theory including the kinetic terms in order to describe physical systems, otherwise the states will not be normalizable. The topological limit

E_{J} \to \infty, E_{C} \to \infty

is not well defined without specifying the value of their ratio g in this limit, and, as we will show, the phase diagram depends crucially on this.

The dual field strengths of the fictitious gauge fields,

j_{μ} = K_{μ ν} b_{ν}

and

φ_{μ} = K_{μ ν} a_{ν}

, represent charge and vortex fluctuations, respectively. The integers fields

Q_{i}

and

M_{i}

are the electric and magnetic topological excitations respectively. Together with the vortex number

M_{0}

, the latter form a 3-current

M_{μ}

which is conserved due to gauge invariance in the

b_{μ}

gauge sector,

{\hat{Δ}}_{μ} M_{μ} = 0

. Perfect duality is broken by the absence of the integer variable

Q_{0}

. As we will see, for the

M_{μ}

variables, the quantity

Δ_{0} M_{0} = \pm 1

corresponds to vortices that appear and disappears on the arrays, possible because topological charges are not conserved. Noether charges coming from symmetry are, instead, conserved and we do not have the variable

Q_{0}

. Contrary to charges, vortices are, thus, topological excitations, characterized by a topological quantum number. The configuration space of the theory of vortices decomposes into so-called superselection sectors, characterized by the integer total vortex number, which are connected via instantons, non-perturbative configurations representing quantum tunneling events between topological vacua [37]. As a consequence, charges are conserved but vortices are not and can “appear" and“disappear" via quantum tunneling events forming the instantons.

By rescaling the emergent gauge fields, expressed in their canonical dimensions, by

1 / 2 π

, and using lattice derivatives instead of finite differences, we can rewrite the action in eq. (13) as

S = \sum_{x} i \frac{l_{0} l^{2}}{2 π} a_{μ} K_{μ ν} b_{ν} + \frac{l_{0} l^{2}}{8 π^{2} E_{J}} {j_{i}}^{2} + \frac{l_{0} l^{2}}{16 E_{C}} ϕ_{i}^{2} + i l a_{i} Q_{i} + i l_{0} b_{0} M_{0} + l b_{i} M_{i} .

(15)

This equation is the limit for

ϵ = 1, μ \to \infty

and

v = 1 / \sqrt{ϵ μ} = 1 / \sqrt{μ} ≪ 1

of the action:

\begin{matrix} S = \sum_{x} & i \frac{ℓ^{2} l_{0}}{2 π} a_{μ} k_{μ ν} b_{ν} + \frac{ℓ^{2} l_{0} v^{2}}{8 π^{2} E_{J}} j_{0} j_{0} + \frac{l_{0} l^{2}}{8 π^{2} E_{J}} {j_{i}}^{2} + \\ \frac{ℓ^{2} l_{0} v^{2}}{16 E_{C}} φ_{0} φ_{0} + + \frac{l_{0} l^{2}}{16 E_{C}} ϕ_{i}^{2} + i l a_{i} Q_{i} + i l_{0} b_{0} M_{0} + l b_{i} M_{i} . \end{matrix}

(16)

We want now to compute the induced action for the

Q_{i}

and

M_{μ}

obtained by integrating over the emergent gauge fields. To this end we will add the term

i l_{0} a_{0} Q_{0}

, and the sum over

Q_{0}

in the partition function. This term will make the action completely self-dual with respect to electric and magnetic degrees of freedom. This will allow us to rewrite eq. (16) as:

\begin{matrix} S = \sum_{x} & i \frac{ℓ^{2} l_{0} v}{2 π} a_{μ} k_{μ ν} b_{ν} + \frac{ℓ^{2} l_{0} v}{8 π^{2} E_{J}} j_{μ} j_{μ} + \frac{ℓ^{2} l_{0} v}{16 E_{C}} φ_{μ} φ_{μ} + \\ + i l a_{μ} Q_{μ} + i l b_{μ} M_{μ}, \end{matrix}

(17)

where we rescaled

l_{0} = l / v, d_{0} \to d_{0}, j_{i} (φ_{i}) \to \frac{1}{v} j_{i} (φ_{i})

. In eq. (17) gauge invariance imposes the constraint

\hat{d_{μ}} Q_{μ} = \hat{d_{μ}} M_{μ} = 0

, but, since

\hat{d_{i}} Q_{i} = 0

we also have

\hat{d_{0}} Q_{0}

, implying that there are no “electric instantons”. To recover the exact result for JJA we should put

Q_{0} = 0

at the end of the calculation. Making the theory perfectly self-dual, however, will not change the nature of the possible phases, we thus use this approximation.

Integrating out the emergent gauge fields we get [5,6]:

\begin{matrix} S_{TOP} & = \sum_{x} v \frac{8 E_{C}}{ℓ} Q_{μ} \frac{1}{v^{4} m^{2} - d_{0} {\hat{d}}_{0} - v^{2} \nabla_{2}^{2}} Q_{μ} + \\ + v \frac{4 π^{2} E_{J}}{ℓ} M_{μ} \frac{1}{v^{4} {\tilde{m}}^{2} - d_{0} {\hat{d}}_{0} - v^{2} \nabla_{2}^{2}} M_{μ} + \\ + i \frac{2 π v^{6} m^{2}}{ℓ} Q_{μ} \frac{K_{μ ν}}{(d_{0} {\hat{d}}_{0} + v^{2} \nabla_{2}^{2}) (v^{4} m^{2} - d_{0} {\hat{d}}_{0} - v^{2} \nabla_{2}^{2})} M_{ν} . \end{matrix}

(18)

The last imaginary term in the action is a lattice version of the topological linking of electric and magnetic strings, which, due to the Dirac condition, in the limit

v m \to \infty

becomes an integer. This term is an integer also if we set

Q_{0} = 0

. We can, thus, drop it [5].

3. Quantum Phase Structure

The

T = 0

quantum phase structure of JJA is determined the behaviour of the integer fields

Q_{μ}

and

M_{μ}

. These can be decomposed into a transverse component,

{\hat{d}}_{μ} Q_{μ}^{T} = 0

{\hat{d}}_{μ} M_{μ}^{T} = 0

, representing closed electric and magnetic loops or infinitely long strings, and a longitudinal component,

{\hat{k}}_{μ ν} Q_{ν}^{L} = 0

{\hat{k}}_{μ ν} M_{ν}^{L} = 0

representing open electric and magnetic strings ending on electric and magnetic monopoles

q = {\hat{d}}_{μ} Q_{μ}^{L}

and

m = {\hat{d}}_{μ} M_{μ}^{L}

. Notice, however, that, in eq. (15) only the integers

Q_{i}

appear implying that

{\hat{d}}_{i} Q_{i} = 0

To determine when infinitely long electric or magnetic strings proliferate, one has to look at the energy-entropy balance determined by the parameters of the model. Let us start from eq. (18), where we ignore the imaginary term, as explained. Near the transition, we expect to have very long strings and large loops. These configuration are, in general, very random and we expect that forces between links in the same loop and in other ones cancel out. We, thus, retain only the self-interaction terms in eq. (18) [38]. We assign to a closed string made of N links and integer quantum numbers Q and

M_{M}

on all the lattice links forming the string and zero elsewhere an energy (equivalent to Euclidean action in statistical field theory)

S_{top} = π m ℓ G (m ℓ) [\frac{2 E_{C}}{π^{2} E_{J}} Q^{2} + \frac{π^{2} E_{J}}{2 E_{C}} M^{2}] N,

(19)

where

G (m ℓ)

is the diagonal element of the lattice kernel

G (m ℓ, x - y)

representing the inverse of the operator

ℓ^{2} (m^{2} - \nabla^{2})

. The kernel

G (m ℓ, x - y)

is defined by the equation

ℓ^{2} (m^{2} - \nabla^{2}) G (m l, x - y) = δ_{x - y, 0} .

(20)

Defining the Fourier transform

G (m ℓ, x - y) = \int_{- π}^{π} d^{3} k G (m ℓ, k) \exp (i k \cdot x)

we obtain

\int_{- π}^{π} d^{3} k G (m ℓ, k) ℓ^{2} (m^{2} - \nabla^{2}) e^{i k \cdot x} = \frac{1}{{(2 π)}^{3}} \int_{- π}^{π} d^{3} k e^{i k \cdot x} .

(21)

Applying, finally, the finite difference operator

ℓ^{2} (m^{2} - \nabla^{2})

to the exponential in the Fourier transform gives the final result

G (m l) \equiv G (m l, 0) = \frac{1}{{(2 π)}^{3}} \int_{- π}^{π} d^{3} k \frac{1}{{(m ℓ)}^{2} + \sum_{i = 0}^{2} 4 \sin {(\frac{k^{i}}{2})}^{2}} .

(22)

The string entropy, however, is also proportional to their length, being given by

μ N

with

μ \approx ln (5)

since at each step the non-backtracking strings can choose among 5 possible directions on how to continue. One can thus assign the free energy

F = π m ℓ G (m ℓ) [\frac{1}{g} Q^{2} + g M^{2} - \frac{1}{η}] N,

(23)

to a string of length

L = ℓ N

carrying electric and magnetic quantum numbers Q and M, respectively. Here we have introduced the dimensionless parameter

η = \frac{π m ℓ G (m ℓ)}{μ},

(24)

which, together with the ratio

g = \sqrt{\frac{π^{2} E_{J}}{2 E_{C}}}

fully determines the quantum phase structure, as we now show.

The ground state of the quantum model is found by minimizing its free energy as a function of N. When the energy term in eq. (23) dominates, the free energy is positive and consequently minimized by short closed loop configurations. When, instead, the entropy dominates, the free energy is negative and minimized by large strings, long closed loops and instantons that break the original

R

gauge symmetry down to

Z

. The condition for condensation of long strings with integer quantum numbers Q and M is thus given by

η \frac{1}{g} Q^{2} + η g M^{2} < 1 .

(25)

If two or more condensations are allowed, one has to choose the one with the lowest free energy.

This condition describes the interior of an ellipse with semiaxes

\begin{matrix} r_{Q} = \sqrt{g \frac{1}{η}}, \\ r_{M} = \sqrt{\frac{1}{g} \frac{1}{η}}, \end{matrix}

(26)

on a square lattice of integer electric and magnetic charges. The phase diagram is consequently found by simply recording which integer charges lie within the ellipse when the semi-axes are varied,

\begin{matrix} η & < 1 \to \{\begin{matrix} \sqrt{\frac{2 E_{C}}{π^{2} E_{J}}} < 1, electric condensation = type III superconductor, \\ \sqrt{\frac{2 E_{C}}{π^{2} E_{J}}} > 1, magnetic condensation = superinsulator, \end{matrix} \\ η & > 1 \to \{\begin{matrix} \sqrt{\frac{2 E_{C}}{π^{2} E_{J}}} < \frac{1}{η}, electric condensation = type III superconductor, \\ \frac{1}{η} < \sqrt{\frac{2 E_{C}}{π^{2} E_{J}}} < η, no condensation = bosonic topological insulator / Bose metal, \\ \sqrt{\frac{2 E_{C}}{π^{2} E_{J}}} > η, magnetic condensation = superinsulator . \end{matrix} \end{matrix}

For

η > 1

the SIT in JJA occurs via an intermediate phase in which both

Q_{μ}

and

M_{μ}

are diluted and the action for JJA reduces to the action of a bosonic topological insulator [39]:

S = \sum_{x} i 2 π a_{μ} K_{μ ν} b_{ν} .

(27)

The existence of this phase was first predicted in [5]. The SIT point [40] is known to be an universal quantum phase transition point with resistance R equal to the quantum of resistance

R_{Q}

. This topological state is often referred to as a Bose metal because it hosts symmetry-protected metallic edge states [30]. It forms due to the same competition of two quantum orders, the charge condensate and the vortex condensate that leads to the universal SIT point. This intermediate phase is separated by two quantum BKT transition from the superinsulating phase for

g < 1

and from the superconducting phase for

g > 1

[21]. For

g \geq 1

R \leq R_{Q}

, the saturation of the resistance has been experimentally observed in [41] in Al/InAs JJA using a gate voltage, which suppresses charges tunneling, thereby confirming the validity of our model. In terms of the array parameters, this means varying

E_{J}

keeping

E_{C}

fixed to increase the value of g.

The superconducting phase is the phase in which the electric topological defects proliferate, and it is this new type of superconductivity that we called type III [18], characterized by a topological gap and superconductivity lost trough a BKT transition both at

T = 0

and at finite T, as discussed in the introduction.

For

η < 1

, instead, when

\sqrt{2 E_{C} / π^{2} E_{J}}

lies between

η

and

1 / η

there is coexistence of electric and magnetic excitations indicating a direct first order phase transition between the superinsulator and the superconductor .

3.1. Superinsulators

The superinsulating phase is the phase in which the magnetic monopoles defects

M_{μ}

are dense. To understand the nature of this phase we will compute the induced effective action for the electromagnetic gauge potential

A_{μ}

minimally coupled to the electric

S \to S + i \sum_{x} l_{0} l^{2} A_{μ} j_{μ} = S + i \sum_{x, i} (l^{2} b_{0} F_{0} + l_{0} l b_{i} F_{i}),

(28)

where

F_{μ} = {\hat{k}}_{μ ν} A_{ν}

. The induced action is:

e^{- S_{eff} (A_{μ})} = \sum_{Q_{μ}, M_{μ}} \int_{a_{μ}, b_{μ}} D a_{μ} D b_{μ} e^{- S (a_{μ}, b_{μ}, Q_{μ}, M_{μ}, A_{μ})},

(29)

and the corresponding induced current is:

j^{μ} = \frac{1}{ℓ^{3}} \frac{δ}{δ A_{μ}} S_{eff} (A_{μ}) .

(30)

By noticing that the coupling with the external electromagnetic field can be taken into account by shifting

M_{0}

and

M_{i}

\sum_{x} (l_{0} b_{0} M_{0} + i b_{i} M_{i})

to:

M_{0} \to M_{0}^{'} = M_{0} + \frac{l^{2}}{2 π} F_{0}, M_{i} \to M_{i}^{'} = M_{i} + \frac{l l_{0}}{2 π} F_{i},

(31)

we obtain from eq. (18):

\begin{matrix} S_{top} (Q_{μ}, M_{μ}, A_{μ}) & = \sum_{x, i} ℓ_{0} ℓ^{2} g \frac{μ η ℓ}{4 π^{2}} [v {(F_{0} + \frac{2 π}{ℓ^{2}} M_{0})}^{2} + \frac{1}{v} {(F_{i} + \frac{2 π}{ℓ_{0} ℓ} M_{i})}^{2}] \\ + \sum_{x} \frac{1}{g} μ η Q_{μ} Q_{μ} + \sum_{x, i} i \frac{μ η m v ℓ}{2 π} (ℓ_{0} A_{0} Q_{0} + ℓ A_{i} Q_{i}), \end{matrix}

(32)

where, also in this case, we retained only self-interactions. This corresponds to the limit

ℓ_{0} ω_{P} = l_{0} m ≫ 1

In the superinsulating phase, the electric topological excitations

Q_{i}

are suppressed because of their large energy, so that

Q_{i} = 0

and we restore

Q_{0} = 0

. We obtain, thus, from eq. (32):

\begin{matrix} \sum_{{M_{i}}} \int_{- π}^{+ π} D A_{μ} e^{S (A_{μ}, M_{i})}, \\ S (A_{μ}, M_{i}) = \frac{g}{4 π ℓ_{0} ω_{P}} \sum_{x} {(F_{i} - 2 π M_{i})}^{2}, \end{matrix}

(33)

where we used

m = ω_{P}

and finite lattice differences. In eq. (33) the gauge fields

F_{i}

are periodic under the shift

F_{i} \to F_{i} + 2 π M_{i}

and are compact, they are angular variables defined on the interval

[- π, + π]

. What we obtain is the deep non-relativistic limit of Polyakov’s compact QED action [42,43], in which only electric fields survives. It is the compacteness of the gauge fields that allows the presence of magnetic monopoles (instantons). As in the relativistic case, as we will show, the presence of monopoles gives linear confinement of probe charges, which become bound by electric flux tubes. There is, however, a crucial difference with the relativistic case. To see this, we decompose

M_{i}

into its transverse and longitudinal components [42,43]:

\begin{matrix} M_{i} = {M^{T}}_{i} + {M^{L}}_{i}, {M^{T}}_{i} = ϵ_{i j} Δ_{j} n + ϵ_{i j} Δ_{j} ξ, n \in Z \\ {M^{L}}_{i} = Δ_{i} λ, Δ λ = {\hat{Δ}}_{i} Δ_{i} λ = m, \end{matrix}

(34)

where m are integer magnetic monopoles. The sum over

{M_{i}}

can be traded for the sum over

{n}

and

{q}

. The integers

{n}

are used to shift the integration domain for the gauge field

A_{μ}

[- \infty, + \infty]

, and the real variables

{ξ}

can be absorbed into the gauge field, giving an integral over the non-compact gauge field

A_{μ}

and a sum over the monopoles degrees of freedom in the partition function:

Z = Z_{0} Z_{inst} = Z_{0} \sum_{{m}} e^{- \frac{π g}{ℓ_{0} ω_{P}} \sum_{x} m_{x} \frac{1}{- Δ} m_{x}}

(35)

where

Z_{0}

is the Gaussian integral over

A_{μ}

. While in the relativistic case monopoles interact with a potential

1 / x

and are, thus, always in a plasma phase and always confine charges, in JJA, monopoles interact with the inverse of the spatial Laplacian, which gives a

(e_{eff}^{2} / 2 π) \log | x |

potential in 2d with

e_{eff}^{2} = π^{2} l_{0} ω_{P} / g

. Monopoles undergo, thus, a quantum BKT transition with g playing the role of a inverse temperature: for low values of g instantons are free and confine probes charges; instantons undergoes a confinng transition, which is actually the SIT, where probe charges are liberated.

To see how instantons modify the Coulomb potential and cause linear confinement of probe charges, we compute the expectation value of the Wilson loop operator

W (C)

, whith C a closed loop in 3D Euclidean space-time restricted to the plane formed by the Euclidean time and one of the space coordinates, which gives the interaction potential between to external probe charges of strength

\pm q_{ext}

〈 W (C) 〉 = \frac{1}{Z_{A_{μ}, M_{i}}} \sum_{{M_{i}}} \int_{- π}^{+ π} D A_{μ} e^{- \frac{g}{4 π ℓ_{0} ω_{P}} \sum_{x} {(F_{i} - 2 π M_{i})}^{2}} e^{i q_{ext} \sum_{C} A_{μ}},

(36)

where we absorbed a factor l in

A_{μ}

. A linear interaction between the probe charges will give rise to an area-law: [37,42]

〈 W (C) 〉 = e^{- σ A},

(37)

where A is the area of the surface S enclosed by the loop C and

σ

is called the string tension [37,42], which we will compute in what follows. Using the lattice Stoke theorem, one rewrites Eq. (36) as

〈 W (C) 〉 = \frac{1}{Z_{A_{μ}, M_{i}}} \sum_{{M_{i}}} \int_{- π}^{+ π} D A_{μ} e^{- \frac{g}{4 π ℓ_{0} ω_{P}} \sum_{x} {(F_{i} - 2 π M_{i})}^{2}} e^{i q_{ext} \sum_{S} S_{i} (F_{i} - 2 π M_{i})},

(38)

where the quantities

S_{i}

are unit vectors perpendicular to the plaquettes forming the surface S encircled by the loop C and vanish on all other plaquettes. We have also multiplied the Wilson loop operator by 1 in the form

\exp (- i 2 π q_{ext} M_{i})

on all plaquettes forming S.

Repeating the steps that lead to eq. (35) we obtain:

〈 W (C) 〉 = \frac{1}{Z_{m}} \sum_{{m}} e^{- \frac{π g}{ℓ_{0} ω_{P}} \sum_{x} m_{x} \frac{1}{- Δ} m_{x}} e^{i 2 π q_{ext} \sum_{S} {\hat{Δ}}_{i} S_{i} \frac{1}{- Δ} m_{x}},

(39)

where we negleted

Z_{0}

, since the integral over non-compact gauge field

A_{μ}

gives then the Gaussian fluctuations around the instantons m, which do not contribute to confinement [6].

We will now perform the sum over the instantons and introduce an auxiliary scalar field

χ

to write the quadratic term in the instantons in eq. (39) as a Gaussian integral over

χ

. At low g, in the deep superinsulationg regime we will use the dilute instanton approximation, summing only over single monopoles configuration

m = \pm 1

, obtainin thus [42]:

〈 W (C) 〉 = \frac{1}{Z_{χ}} \int D χ e^{- \frac{ℓ_{0} ω_{P}}{4 π g} \sum_{x} Δ_{i} χ Δ_{i} χ + \frac{8 π g}{ℓ_{0} ω_{P}} z (1 - \cos (χ + q_{ext} η))},

(40)

where the angle

η = 2 π {\hat{Δ}}_{i} S_{i} / (- Δ_{2})

represents a dipole sheet on the Wilson surface S and the monopole fugacity z is determined by the self-interaction as

z = e^{- \frac{π g}{ℓ_{0} ω_{P}} G (0)},

(41)

with

G (0)

being the inverse of the 2D Laplacian at coinciding arguments. Eq. (40) can be rewritten as:

〈 W (C) 〉 = \frac{1}{Z_{χ}} \int D χ e^{- \frac{ℓ_{0} ω_{P}}{2 π g} \sum_{x} \frac{1}{2} Δ_{i} (χ - q_{ext} η) Δ_{i} (χ - q_{ext} η) + μ^{2} (1 - \cos (χ))},

(42)

where we wrote:

μ^{2} = 4 π g z / ℓ_{0} ω_{P}

and shift

χ \to χ - q_{ext}

To compute the string tension

σ

we have to evaluate the integral eq. (42). We will do it both for Cooper pairs, with

q_{ext} = 1

and single electrons, corresponding to

q_{ext} = 1 / 2

in our case. We will evaluate eq. (42) using the saddle point approximation:

Δ_{2} χ_{cl} = q_{ext} Δ_{2} η + μ^{2} \sin χ_{cl},

(43)

valid for small g where the integral is dominated by the classical solution to the equation of motion and which reduces to a one-dimensional equation when we consider time and one spatial coordinate, e.g. in the

(t, x_{2})

plane far from the boundaries of S:

{\hat{Δ}}_{x 1} Δ_{x 1} χ_{cl} = - 2 π q_{ext} {\hat{Δ}}_{x 1} S_{1} + μ^{2} \sin χ_{cl} .

(44)

To solve this equation we will use the continuum limit [42] with boundary conditions

χ_{cl} \to 0

for

| x_{1} | \to \infty

. For Cooper pairs (

q_{ext} = 1

) we obtain:

χ_{cl} = sign (x_{1}) 4 \arctan e^{- μ | x_{1} |},

(45)

which gives eq. (37) with

σ

σ = \frac{ℏ ω_{P}}{ℓ} \sqrt{\frac{16}{π g ℓ_{0} ω_{P}}} \sqrt{z} = \frac{ℏ ω_{P}}{ℓ} \sqrt{\frac{16}{π g ℓ_{0} ω_{P}}} e^{- \frac{π g}{2 ℓ_{0} ω_{P}} G (0)} .

(46)

The string binds together charges, prevents charge transport on arrays of a sufficient size and is the origin of the infinite resistance characterizing superinsulation. If we consider, instead, a single electron probes,

q_{ext} = 1 / 2

, with the same boundaries conditions we obtain, for the string tensions

σ_{electrons} = \frac{1}{2} σ,

(47)

which implies that single electrons are also confined. This explains why charge transport mediated by thermally excited normal quasiparticles is not present in superinsulators.

From eq. (46) we can estimate the typical string size

ℓ_{string} = \sqrt{c ℏ / σ}

. Taking the following typical values for experimental JJA,

ℓ = 100

nm and

ω_{P} = 10

GHz, and taking for

ℓ_{0} ω_{P} = O (1000)

we arrive at

ℓ_{string} / ℓ \approx 150

lattice spacing, which represents the distance between the superconducting islands. This set a minimum dimension for the arrays to be able to accommodate an electric pion and to show superinsulations. To be able to see superinsulation in JJA, moreover, it is necessary to lower the value of g, that, in terms of the parameter of the arrays, implies increasing the vortex tunneling parameter

E_{C}

. While making

g > 1

has been experimentally achieved in [41], the opposite limit seems to be more difficult. It will require a "more insulating" substrate to be able to increase vortex mobility and govern

E_{C}

. Our predictions is that in this case it will be possible to observe superinsulation in JJA.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

Fazio, R. and van der Zant, H. Quantum phase transitions and vortex dynamics in superconducting networks. Physics Reports 355, 235 – 334 (2001). [CrossRef]
Fazio, R. et al., Helv. Phys. Acta 65, 228 (1992).
Penttila, J.S. et al., Phys. Rev. Lett. 82, 1004 (1999).
J.M. Kosterlitz, D.J. Thouless, Long range order and metastability in two dimen- sional solids and superfluids. Application of dislocation theory, Journal of Physics C: Solid State Physics 5 (1972) L124. [CrossRef]
Diamantini, M.C., Sodano, P. and Trugenberger, C. Gauge theories of josephson junction arrays.Nuclear Physics B 474, 641 – 677 (1996). [CrossRef]
Trugenberger, C., Diamantini, M. C., Poccia, N., Nogueira, F.S. and Vinokour, V.M., Magnetic monopoles and suprinsulationin Josephson junction arrays, Quant. Rep. 2, 388-399, (2020). [CrossRef]
Diamantini, M. C., Trugenberger, C. A. and Vinokur, V. M. Confinement and asymptotic freedom with cooper pairs. Communications Physics 1, 77 (2018). [CrossRef]
Krämer, A. and Doniach, S. Superinsulator phase of two-dimensional superconductors. Phys. Rev. Lett. 81, 3523–3527 (1998). [CrossRef]
Sambandamurthy, G., Engel, L. W., Johansson, A., Peled, E. and Shahar, D. Experimental evidence for a collective insulating state in two-dimensional superconductors. Phys. Rev. Lett., 94, 017003 (2005). [CrossRef]
Baturina, T. I. et al. Localized superconductivity in the quantum-critical region of the disorder-driven superconductor-insulator transition in TiN thin films. Phys. Rev. Lett. 99, 257003 (2007). [CrossRef]
Vinokur, V. M. et al. Superinsulator and quantum synchronization. Nature 452, 613–615 (2008). [CrossRef]
Baturina, T. I. and Vinokur, V. M. Superinsulator–superconductor duality in two dimensions. Annals of Physics 331, 236 – 257 (2013). [CrossRef]
Diamantini, M. et al. Direct probe of the interior of an electric pion in a Cooper pair superinsulator, Nature Comm. Phys.3:142, 2020. [CrossRef]
Mironov, A.,Diamantini, M. C., Trugenberger, C. A. and Vinokur, V. M., Relaxation eloctrodynamics of superinsulators, Sci. Rep. 12, 6733 (2022). [CrossRef]
‘t Hooft, G. On the phase transition towards permanent quark confinement. Nuclear Physics B 138, 1 – 25 (1978). [CrossRef]
Deser, S., Jackiw, R. & Templeton, S. Three-dimensional massive gauge theories. Phys. Rev. Lett. 48 975 (1982). [CrossRef]
M. Kalb and P. Ramond, Phys. Rev. D9 (1974) 2273.
Diamantini, M. C., Trugenberger, C. A. and Vinokur, V. M., Type III superconductivity, Advanced Science 2023, 2206523. [CrossRef]
Diamantini, M. C., Trugenberger, C. A. and Vinokur, V. M., How planar superconductors cure their infrared divergences, JHEP 10, 100 (2022). [CrossRef]
Diamantini, M. C., Trugenberger, C. A. and Vinokur, V. M.,Topological Nature of High Temperature Superconductivity, Advanced Quantum Technologies 2000135 (2021). [CrossRef]
Diamantini, M. et al., Bosonic topological insulator intermediate state in the superconductor-insulator transition, Phys. Lett. A 384 (2020) 126570. [CrossRef]
Diamantini, M. C., Trugenberger, C. A. and Vinokur, V. M., The superconductor-insulator transition in absence of disorder, Phys. Rev. B 103, 174516 (2021). [CrossRef]
D. Haviland, Y. Liu, A. Goldman, Onset of Superconductivity in the Two- Dimensional Limit, Physical Review Letters 62 (1989) 2180 – 2183. [CrossRef]
H.M. Jaeger, D.B. Haviland, B.G. Orr, A.M. Goldman, Onset of superconductiv- ity in ultrathin granular metal films, Physical Review B 40 (1989) 183 – 196. [CrossRef]
A. Hebard, M.A. Paalanen, Magnetic-Field-Tuned Superconductor-Insulator Transition in Two-Dimensional Films, Physical Review Letters 65 (1990) 927 – 930. [CrossRef]
N. Mason, A. Kapitulnik, Dissipation effects on the superconductor– insulator transition in 2D superconductors, Physical Review Letters 82 (1999) 5341–5344 (1999). [CrossRef]
N. Markovic , C. Christiansen, A. M. Mack, W. H. Huber, A. M. Goldman, Super- conductor–insulator transition in two dimensions, Physical Review B 60 (1999) 4320–4328. [CrossRef]
M.A. Steiner, N.P. Breznay, A. Kapitulnik, Approach to a superconductor- to-Bose-insulator transition in disordered films, Physical Review B77 (2008) 212501. [CrossRef]
A.T. Bollinger, et al, Superconductor-insulator transition in La2-xSrxCuO4 at the pair quantum resistance, Nature 472 (2011) 458–460. [CrossRef]
Das, D. & Doniach, S. Existence of a Bose metal at T = 0. Phys. Rev. B 60, 1261–1275 (1999). [CrossRef]
Bottcher, C. G. L. et al. Superconducting, insulating and anomalous metallic regimes in a gated two-dimensional semiconductor-superconductor array. Nature Physics 14, 1138-1144 (2018). Milton, K. A. Theoretical and experimental status of magnetic monopoles. Rep. Prog. Phys. 2006, 69, 1637-1712.
Fazio, R. and Schön, G. Charge and vortex dynamics in arrays of tunnel junctions. Phys. Rev. B 43, 5307–5320 (1991). [CrossRef]
Kleinert, H. Gauge Fields in Condensed Matter. Vol. 1: Superflow and Vortex Lines (Disorder Fields, Phase Transitions); Vol. 2: Stresses and Defects (Differential Geometry, Crystal Melting) World Scientific, 1989. [CrossRef]
van Otterlo A., Fazio, R. and Schön G. Quantum vortex dynamics in Josephson junction arrays. Physica B 203, 504–512 (1994). [CrossRef]
Zaikin, A. D., Golubev, D. S., van Otterlo, A. and Zimányi, G. T. Quantum phase slips and transport in ultrathin superconducting wires. Phys. Rev. Lett. 7, 1552–1555 (1997). [CrossRef]
van der Zant, H., Fritschy, F. C., Orlando, T. P. and Mooji, J. E. Ballistic motion of vortices in Josephson jucntion arrays. Europhys. Lett. 18, 343–512 (1992).
Coleman, S. Aspects of symmetry. Cambridge University Press, 1985.
Banks, T., Myerson, R. and Kogut, J., Nucl. Phys. B 129, 493 (1977).
G. Y. Cho and J. E. Moore, Ann. Phys. 326 (2011) 1515.
M.P.A.Fisher,G.Grinstein,S.M.Girvin,PresenceofQuantumDiffusioninTwo Dimensions: Universal Resistance at the Superconductor-Insulator Transition, Physical Review Letters 64 (1990) 587 – 590. [CrossRef]
Bøttcher, C. G. L. et al. Superconducting, insulating and anomalous metallic regimes in a gated two-dimensional semiconductor–superconductor array. Nature Physics, 14, 1138–1144 (2018). [CrossRef]
Polyakov, A. M. Gauge fields and strings. Contemp. Concepts Phys. 3, 1–301 (1987).
Polyakov, A. Compact gauge fields and the infrared catastrophe. Physics Letters B 59, 82 – 84 (1975).

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2020 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Topological Gauge Theory of Josephson Junction Arrays: The Discovery of Superinsulation

Abstract

1. Introduction

2. Topological gauge theory of JJA

3. Quantum Phase Structure

3.1. Superinsulators

Funding

Conflicts of Interest

References

MDPI Initiatives

Important Links

Subscribe