Towards A Non-Relativistic Holographic Superfluid
Towards A Non-Relativistic Holographic Superfluid
Towards A Non-Relativistic Holographic Superfluid
,
hence it is often called the Number operator we prefer Mass to disambiguate the various meanings of number.
1
and reentrance of the normal phase at low temperatures for suciently large background density.
It is tempting to interpret this re-entrance as signaling a zero temperature quantum phase tran-
sition as the background mass density is tuned. However, the re-entrant normal state is again the
simple 1-parameter deformation which we do not expect to be the true equilibrium groundstate, so
we do not expect this probe analysis to be the end of the story. Meanwhile, it remains possible that
the system is in fact reentrant for all values of the background Mass density as T 0, where our
probe approximation becomes unreliable. Resolving these puzzles, however, requires going beyond
the truncated probe approximation discussed in this paper; we leave them to future study.
The plan of the paper is as follows. In Section 2 we quickly describe the basic strategy and
computational setup, with various details elaborated in Appendices. In Section 3 we explore the
phenomenology and phase stucture of holographic superuids outside a Schrodinger black hole (an
analogous study in the background of a Schrodinger soliton [15] is performed in Appendix A while
this is not in the same ensemble as the black hole, it provides an alternate example with surprising
physics of its own). We close in Section 5 with a summary and list of next steps.
2 The Setup
Our basic strategy involves studying an Abelian Higgs system,
L
probe
=
1
e
2
_
1
4
F
2
[T[
2
m
2
[[
2
_
, (2.1)
as a perturbation around the planar Schrodinger black hole background,
ds
2
=
_
f +
(f 1)
2
4(K 1)
_
dt
2
Kr
4
+
1 +f
r
2
K
dt d +
K 1
K
d
2
+
dx
2
r
2
+
dr
2
f r
2
. (2.2)
in the probe limit, e
2
. Here, f = 1 r
4
(T)
4/3
, K = 1 + r
2
2
and the metric is given in
string frame. One can think of this as a rather extreme truncation of the charged Schrodinger black
hole system [16, 17] where we drop the coupling of the vector to the scalar and massive vector of
the black hole background, or simply as a holographic toy model. The geometry is controlled by
two physical parameters, the background mass density, , and the temperature, T,with the horizon
located at the radial coordinate r
H
= (T)
1/3
.
For spatially homogeneous solutions, we can without loss of generality set
A = 0 and take
= (r) and A = A
t
(r)dt +A
f(4 f)r
_
f
_
q
2
A
2
+ 2q
2
r
2
A
A
t
+m
2
K
1/3
_
q
2
(f 1)
2
4(K 1)
_
A
2r
4
H
2
A
t
_
2
_
= 0 (2.3)
fr
2
A
t
_
2
f
3
(7K 4)
_
r
K
A
_
2 +f(f 1) +
(f 1)
2
K 1
_
1
Kr
A
2q
2
K
1/3
2
A
t
= 0 (2.4)
fr
2
A
_
4K 2
2 +K
3
f
_
r
K
A
4(K 1)
r
3
K
A
t
2q
2
K
1/3
2
A
= 0 (2.5)
Note that A
,= A
t
.
2
2.1 Asymptotic Behavior and the Holographic Dictionary
Near the boundary at r = 0, the vector components behave as,
A
t
=
Q
+
Q
r
2
+. . . , A
= M
o
+
M
r
2
+. . . (2.6)
where the various . . . represent various (possibly non-normalizable) terms whose coecients are
entirely xed by the equations of motion and the values of these integration constants,
Q
,
Q
,
M
o
and
M
.
3
As usual,
Q
represents the chemical potential per unit charge, which eectively sets
the zero of energy in the boundary theory the gauge-invariant bulk quantity that becomes the
boundary hamiltonian acting on the operator dual to the bulk matter eld of charge q is (i
t
+q A
t
);
at the boundary, for plane waves e
it
, this becomes (+q
Q
). Thus, one insertion of the charged
operator Oe
it
costs E = ( +q
Q
). As usual,
Q
computes the induced charge density.
It might be tempting to think of M
o
as a chemical potential for the Mass operator,
M. However,
this is not quite right it is a superselection parameter. Recall that, holographically,
M
P
= i(
iqA
)[
, (2.7)
i.e.
M is the boundary value of the gauge-invariant -momentum in the bulk. The mass eigenvalue
of a boundary operator dual to a bulk eld with -momentum and charge q is thus M = (q M
o
),
where M
o
= A
correspond to distinct
NRCFTs, not to a xed theory with dierent background elds turned on. In particular, as we
will see momentarily, the dimensions of various boundary operators depend on A
, an unfamiliar
eect.
M
computes the Mass density coupled to A
.
As for our charged scalar, near the boundary at r 0 it behaves as
1
r
+
2
r
+
+. . . , (2.8)
where
= 2
_
4 +m
2
+q
2
M
2
o
. (2.9)
(Note that we will occasionally write
1
and
2
for
and
+
, respectively.) In the window
1<
1
as the source with
2
O) the response. We will mostly focus on the alternate
quantization with
2
and O)
1
for reasons which will become clear in the next section.
3
In particular, the leading term for At runs as 2M log(r). While formally the dominant term, it is determined
by the equations of motion and M and thus does not represent an independent mode of the system. Importantly, due
to factors of the inverse metric, this log running does not lead any components of the bulk stress tensor to diverge.
A complete holographic renormalization of this system would settle the dictionary, but is beyond the scope of the
present paper; for the moment we simply take the above dictionary as a provisional interpretation which is supported
by the consistency of the results below. Interestingly, while At has no log in fully backreacted charged-black-hole
solutions [16, 17], linearizing the Maxwell equation around these solutions does generate a log without changing any
other of the asymptotics of the vector, so this log is likely a simple consequence of an extreme truncation of the full
charged black hole system. It would be interesting to study the full system and see what, if anything, changes.
3
Importantly, the dimensions,
iqA
)[
, this quantity is nothing but the Mass eigenvalue of the dual operator,
M = qA
= q M
o
.
Our expression for the dimensions above then becomes,
= 2
4 +m
2
+M
2
, which is the
expected form [7, 8], including the quadratic dependence on M inside the radical.
Now, as discussed in [9, 10, 11], the free energy of the full system takes the form, F E+
M
M,
where
M
=
1
2
2
r
4
H
is determined by the background spacetime. The total free energy per insertion
of an operator dual to a bulk eld with -momentum , frequency , and coupled with charge q to
our gauge eld is thus F = ( +q
Q
) +
M
M, where M = ( q M
o
)
2.2 Near-Horizon Behavior and Setting Up the Calculation
In the bulk, we are thus left with a six-parameter family of solutions labeled by sources (
Q
, M
o
, )
and responses (
Q
,
M
, O)). Holographically, we expect boundary conditions at the horizon, where
the radial equations of motion degenerate, to impose three additional constraints. Together with
the two parameters T and of the background geometry, this should leave us with a ve-parameter
phase space. To verify this, we need to study the behaviour of our solutions near the horizon.
The equations of motion degenerate at the black hole horizon, so we must impose boundary
conditions to pick the appropriate solutions. As usual, it suces to impose regularity at the
horizon, which is in any case necessary for the validity of the probe approximation. Assuming
regularity, the equations of motion as presented in (2.3) (2.5) degenerate into three algebraic
equations relating the six horizon values of the elds and their derivatives, as expected,
_
A
(r
H
) 2r
4
H
2
A
t
(r
H
)
_
2
(r
H
) = 0 (2.10)
2r
2
H
A
t
(r
H
) +
_
2K
H
1
K
H
1
_
A
(r
H
) + 2q
2
r
H
K
4/3
H
2
(r
H
)A
t
(r
H
) = 0 (2.11)
4r
H
(r
H
) +
_
4K
H
(K
H
1)r
4
H
q
2
A
2
t
(r
H
) +m
2
K
1/3
_
(r
H
) = 0 . (2.12)
This suggests a simple numerical strategy for constructing superuid states of our holographic
NRCFT. To specify a solution to the full equations of motion, we x any three of
Q
, M
o
, ,
Q
,
M
and O) at the boundary and impose the above regularity conditions at the horizon. Since we are
interested in spontaneously generated condensates, we will generally set = 0. The resulting two-
point boundary value problem can be solved numerically in various ways. The most straightforward
is a brute-force shooting method, as typically employed in the relativistic case.
In sweeping out parameter space, however, we must be careful to vary the parameters of the
NRCFT while holding the NRCFT itself xed i.e., while holding the spectrum of quantum
numbers xed. This is straightforward in AdS, where xing the set of dimensions reduces to xing
the bulk mass m
2
of the bulk scalar. Here, however, the dimension and Mass M of the boundary
scalar operator depend on the asymptotic value of A
as M = q A
and
= 2
4 +m
2
+M
2
.
Before sweeping out parameter space, then, we must x A
=M
o
. As weve already set = 0,
xing the system thus leaves us with a three-parameter phase space labeled by
Q
, and T.
4
This peculiar behavior that the denition of the boundary CFT depends on the boundary
behavior of the bulk elds is a very general phenomenon in Schrodinger holography. Indeed,
renormalizing the boundary stress tensor, say, or other operators in the boundary NRCFT, requires
counterterms which are local in time and space, but which depend explicitly on dimensions, ,
and thus on the asymptotic values of A
_
4 f
2 + 7K
3K
_
ra
x
+
_
2
r
4
(K 1)
f
2q
2
K
1/3
2
_
a
x
= 0 . (2.13)
Near the horizon, this reduces to,
(
2
+ (
4
r
3
H
d
d
)
2
)a
x
() 0 (2.14)
where r = r
H
. The infalling solutions thus takes the form,
A
x
= a
0
e
it
(r r
H
)
i
4T
(1 +a
1
(r r
H
) +. . . ) (2.15)
Near the boundary,
A
x
= A
0
+A
2
r
2
2
+. . .
A short computation then veries that the conductivity is given by,
() =
J
x
)
E
x
)
= i
J
x
)
A
x
)
= i
A
2
A
0
Note that, since we are solving a linear equation but only care about this ratio, the overall scale of
A
x
is immaterial. We can use this freedom to set a
0
= 1, which simplies the numerical problem.
Notably, we can analytically determine the -dependence of for large and small via standard
power-series analysis. Importantly, the scaling in the superuid phase to be independent of and
M. At small frequency, we nd,
Im[( 1)]
1
(2.16)
while for large we have,
Re[( 1)]
1/3
Im[( 1)]
1/3
(2.17)
This last result unsurprisingly diers from the AdS case, where Re[( 1)] = 1. Reassuringly,
they both match numerical results presented below, a nice sanity check.
5
3 Phases of a Schr odinger Superuid
Thus armed, we now get down to the business of nding a superuid state in our non-relativistic
holographic CFT and exploring its phase diagram. A priori, the phase space is fairly high-
dimensional specifying a point involves xing and M to x the theory, then tuning
Q
, T
and to sweep out the phase diagram. For simplicity, we will begin by picking convenient values
= 6/5, M = 1/2 and
Q
= 1/8, then dial the background Mass density, . This will reveal a
zero-temperature quantum phase transition at a critical value
For very small , the geometry remains essentially AdS until very close to the boundary, so we
expect most low-energy physics such as superuid condensation to very closely track familiar
AdS results. This turns out to be almost correct, modulo a surprise well explore shortly.
0.0 0.2 0.4 0.6 0.8 1.0 1.2
0.00
0.02
0.04
0.06
0.08
0.10
TT
c
O
1
0 2 4 6 8 10
0.00
0.05
0.10
0.15
Tc
Re
1.26
1.05
1.
0.96
0.88
0.65
0.37
0.29
0.24
0.19
0.16
0.08
0.05
0.01
TTc
0 2 4 6 8 10
0.0
0.2
0.4
0.6
0.8
1.0
Tc
Im
1.26
1.05
1.
0.96
0.88
0.65
0.37
0.29
0.24
0.19
0.16
0.08
0.05
0.01
TTc
Figure 1: At small , the behavior of the superuid is essentially the same as in AdS, with a 2
nd
order mean-eld phase transition at the onset of superconductivity at T
c
, including the familiar
gap-and-pole form in the AC conductivity, leaving us in a happy superuid state at T = 0. Here,
(
Q
, ) = (1/8, 1/16), with T
C
= 0.505.
Figure 1a shows the condensate as a function of temperature for the alternate quantization
(O)
1
) for = 1/16. As is clear by eye and can be checked precisely from the numerics,
the resulting condensate turns on at T = T
c
with classic mean-eld behavior (
c
=
1
2
) and grows
as the temperature is lowered. Figures 1b and 1c then show the real and imaginary parts of
the AC conductivity for various temperatures indicated by color, from high (violet) to low (red).
These demonstrate the appearance of a superconducting state at T
c
, with the gap growing as the
temperature is lowered. Note, too, that the conductivity in the superuid phase has Im[(
0)] 1/, while Re[( )] Im[( )]
1/3
. This scaling is expected on general
grounds, so gives us condence in our numerical results.
6
As we increase the background number density, the story changes dramatically. Figure 2 shows
the same plots as Figure 1 but with = 1 rather than = 1/16. The most obvious dierence
is that the order parameter vanishes at suciently low temperature, T T
L
, doing so again with
mean-eld behavior. As is clear form the nite value of Re[(0)], the extreme low-temperature
phase is again metallic.
0.0 0.2 0.4 0.6 0.8 1.0 1.2
0.00
0.01
0.02
0.03
0.04
0.05
0.06
0.07
TT
c
O
1
0 2 4 6 8 10
0.0
0.5
1.0
1.5
Tc
Re
1.28
1.
0.95
0.72
0.36
0.24
0.15
0.1
0.07
0.06
0.01
TTc
0 2 4 6 8 10
0.0
0.2
0.4
0.6
0.8
1.0
Tc
Im
1.28
1.
0.95
0.72
0.36
0.24
0.15
0.1
0.07
0.06
0.01
TTc
Figure 2: At large , in addition to the original transition to a superconducting state at T
c
, the
system now exhibits reentrance of the normal phase at a new low-temperature 2
nd
order transition
at T
L
, again with mean-eld exponents. Below T
L
, behaves like the normal gas. Here (
Q
, ) =
(1/8, 1) with T
C
= 0.149 and T
L
= 0.009
Consider now the behavior of the system at zero temperature as a function of the background
number density, . As 0, the system is superconducting. As 1, the metallic phase is reen-
trant. At some critical
. Figure 3
shows the same system at slightly above
0 2 4 6 8 10
0.0
0.5
1.0
1.5
T
c
Re
1.28
1.01
0.95
0.6
0.33
0.2
0.06
0.02
0.01
0.
TT
c
0 2 4 6 8 10
2
1
0
1
2
3
T
c
Im
1.28
1.01
0.95
0.6
0.33
0.2
0.06
0.02
0.01
0.
TT
c
Figure 3: At intermediate , we again have a 2
nd
order mean-eld transition into a superuid
state at T
c
. At low temperatures, however, the system undergoes a non-mean-eld transition to an
apparently insulating state. Here (
Q
, ) = (1/8, 3/8) with T
C
= 0.123.
standard mean-eld behavior. The zero temperature behavior, however, diers dramatically from
mean-eld expectations; rather, at low temperature, the condensate decays exponentially, as does
the superuid density, while the normal density remains vanishing and the conductivity heavily
suppressed at small but non-vanishing , suggesting that the T = 0 state is not metallic. It is
tempting to read this as indicating a translationally-invariant insulating phase.
7
However, numerical results in this region should be taken with a sizeable grain of salt. Indeed,
at suciently low temperature, the numerics simply fail to converge. More physically, in this
regime, the probe approximation is becoming dangerously unreliable the matter eld proles
which generate the required boundary values grow rapidly deep in the bulk (and in particular
near the horizon) as we approach T = 0 or
, so we can be
quite condent that the system is denitely re-entrant at suciently large .
3.2 High Temperature Condensates and the Free Energy
The surprise alluded to above involves the high-temperature limit. Figure 4 shows the same system
but now extending to higher temperatures. The surprise is the appearance of a high tempera-
ture condensate at T T
H
. Troublingly, the condensate appears to grow without bound as the
temperature increases.
4
0 1 2 3 4 5
0.00
0.02
0.04
0.06
0.08
0.10
0.12
T
O
1
2
1
2
_
M
(
M
Q
+
Q
M
o
) Moreover, if we hold xed the asymptotic values of A
i
(corresponding to
xing the values of the chemical potential
Q
=A
t
[
2
1
We can thus compute the relative free energy density (T
A
T
B
)
between any two states A and B connected by such a ow by integrating S along the ow,
T
B
T
A
= T
_
B
A
S
E
V
D
= T(
1
2
)
_
B
A
2
d
1
(3.1)
where V
D
is the volume of the boundary theory and the integral is performed along the ow specied
above. By construction, this agrees with what we would get by evaluating the fully holographically
renormalized free energy for each solution and subtracting. Happily, this allows us to compute the
correct free energy without having to worry about the full holographic renormalization of the theory
(for further comparison between holographic renormalization and our method, see [14],[21],[22]).
Now consider the case of our holographic superuid in alternate quantization, where
2
= is the
source and
1
=O) is the response. In this case, the curve
1
(
2
) is multi-valued over
2
=0, with
one solution corresponding to the trivial vacuum, O)=0, and one to the nontrivial condensate,
O), =0. As outlined above, these two solutions are connected by a very specic ow in the (
1
,
2
)
plane. To compute the properly renormalized relative free energy, then, all we must do is nd this
ow and integrate along it,
T
C
T
N
= T(
1
2
)
_
C
N
2
d
1
(3.2)
where the integration is again along the ow dened above. If this dierence is negative, the
condensate is thermodynamically favored.
Figure 5 plots two such ows. On the left we have a ow connecting the trivial vacuum (
1
=
2
=0)
to a symmetry-breaking vacuum (
1
,=0,
2
=0) in the alternate quantization in the low-temperature
regime, with the ow indicated by the solid line and the direction of ow dening the direction of
integration. The area under the curve, corresponding to the free energy of the condensed state, is
negative. On the right is the analogous ow in the high-temperature regime here the free energy
is positive. We thus deduce that the low-temperature condensate is thermodynamically stable,
while the high-temperature condensate is unstable, in the alternate quantization.
What about the standard quantization? Figure 6a plots the condensate in standard quantization
as a function of temperature. Note that there is no separate low vs high temperature condensate,
just a single continuous instability whose prole grows with temperature. Figure 6b then shows a
typical ow at typical temperature. Importantly, the enclosed area is negative for every temper-
ature, indicating a thermodynamic instability even at arbitrarily high temperature. This is why
we quietly chose the alternate quantization in Section 2. It would be interesting to understand the
meaning of this instability in detail.
9
0.00 0.05 0.10 0.15
0.6
0.4
0.2
0.0
2
0.00 0.05 0.10 0.15
0.0
0.1
0.2
0.3
2
Figure 5: (a) Low T condensate has smaller free energy than non-condensed phase, T
C
T
N
< 0.
(b) High T condensate has larger free energy than non-condensed phase, T
C
T
N
> 0. One can
determine the sign of T
C
T
N
from the orientation of the curve. Here (
Q
, ) = (3/8, 1/16).
0 50 100 150
0
200
400
600
800
1000
1200
1400
T
O
2
2
135
0 2 4 6 8
0
200
400
600
800
1000
1200
1400
2
135
Figure 6: (a) Condensate as a function of T for the standard quantization, with
2
= 13/5. (b)
Typical ow at generic temperature, indicating a thermodynamic instability at every temperature.
3.3 Varying
Q
and a Multicritical Point
The thermodynamic instability of the high-temperature condensate leads to an important physical
eect as we vary
Q
. Figure 7a plots T
c
(
Q
) and T
H
(
Q
), the critical temperatures for the low-
and high-temperature condensates as a function of the chemical potential
Q
. As we crank up
Q
,
holding all other parameters xed, T
c
increases while T
H
decreases. At a critical value,
, the
two critical points merge; above
.
However, we have already checked that the condensed phase is thermodynamically disfavored at
high temperatures. For
Q
>
T
c
2
nd
order phase transition
1
st
order phase transition
0.0 0.1 0.2 0.3 0.4
0
1
2
3
4
5
Q
T
0 1 2 3 4 5
0.00
0.05
0.10
0.15
0.20
T
O1
TH Q
2
nd
Tc Q
1
st
T Q
716
38
516
14
732
1364
0.195
0.192
0.192
0.191
316
18
0
Q
O1Tc
2
nd
order phase transition
O1T 1
st
order phase transition
Figure 7: At
Q
>
2
0.00 0.05 0.10 0.15 0.20 0.25
4
2
0
2
4
2
Figure 8: (a) Flow lines in the neighborhood of the critical temperature, for
Q
= 3/8. (b) The ow
switches direction discontinuously (lower solid to upper dashed curves) at a critical temperature T
indicated by the red dashed curve in Figure 7, leading to a rst order phase transition at T
. Here,
Q
=3/8 and =1/16.
temperature. Figure 8a shows the ows associated to points to the left and right of the critical
temperature where the low- and high-temperature instabilities meet (indicated by the red dashed
curve in Figure 7b). Figure 8b focuses in on the immediate neighborhood of the transition tem-
perature for
Q
>
. For all temperatures below the critical temperature, the ows go below the
horizontal axis, corresponding to a negative free energy and a thermodynamically stable conden-
sate. For all temperatures above the critical temperature, the ows go above the horizontal axis,
so the condensate is thermodynamically disfavored at high temperatures. Indeed, while the value
of the condensate is non-vanishing and in fact completely smooth as we ow through T
, the path
which carries us from the trivial state to the condensate changes discontinuously as we pass through
T
. As a result, the integrated area and thus the free energy also changes discontinuously at
T
. Moreover, as we take
Q
)) 0; this ensures that the latent heat of the transition goes to zero at the multicritical
point where the transition switches from 1
st
to 2
nd
order, as expected on general grounds.
The upshot of all of the above is that as we raise
Q
, the phase transition from high-temperature
11
metal to low-temperature superuid switches from 2
nd
order to 1
st
, with the transition occurring
at a multicritical point where the low- and high-temperature superuid phases collide. Near the
phase transition boundaries, including the multicritical point, the order parameter scales with
simple mean-eld exponents. More precisely, near the nite temperature 2
nd
order phase transition,
O) (T
c
T)
1/2
, while near the 1
st
order phase transition boundary when
Q
>
the condensate
jumps discontinuously at T
, with O(T
)) (
Q
)
1/2
. This can be succinctly encoded in a
simple mean-eld free energy, F() =
1
2
c
2
(T T
c
(
Q
))
2
+
1
4
c
4
(
Q
)
4
+
1
6
c
6
6
, with O)
and with coecients c
2
, c
4
, c
6
> 0.
3.4 Setting
Q
= 0 and the Persistence of Condensates
Playing with
Q
raises another interesting point. Fundamental to our construction is that the
scalar operator carries a charge q under a global symmetry of the boundary theory.
Q
tells us the
energy cost for adding a unit of this charge to the system. In AdS, superuid condensation is often
induced by tuning
Q
beyond a threshold. Is this also necessary in the non-relativistic case?
0.000 0.002 0.004 0.006 0.008 0.010
0.000
0.002
0.004
0.006
0.008
T
O
1
2
and f
s
= 1 r
4
/r
4
s
, with r
s
controlling both the gap and the radius of the
-direction, L
=
rs
. Here, the radial direction is cut o smoothly by a spacelike circle shrinking
rather than by a black hole horizon and indeed this solution was obtained by double Wick rotating
the Schrodinger black hole with compact -direction (which is spacelike near the horizon).
6
As in the black hole case, a probe superuid in the soliton geometry is characterized by ve
parameters: two dene the theory (the dimension and Mass M), two are properties of the
background which x thermodynamics quantities (the mass density and mass gap of the soliton,
m
G
=
1
L
Q
>
c
, by contrast, Im[( 0)]
1
Q
O
1
0 5 10 15 20
0.5
0.0
0.5
Im
0 5 10 15 20
0.2
0.0
0.2
0.4
Im
Figure 10: (a) We nd a 2
nd
order phase transition at
c
= 1.17. (b) When
Q
<
c
, Im[( 0)]
is nite, indicating an insulating phase. (c) For
Q
>
c
, we nd a superconducting
1
pole, as well
as two gapped poles at nite .
A.2 Varying
In the black hole case, tuning drove us through a superuid-conductor phase transition at zero
temperature. The (zero-temperature) soliton shows the same eect, with the spontaneous conden-
sate disappearing in a 2
nd
-order transition as passes through a critical value,
c
0.163, as
shown in Figure 11a. Here, however, the normal phase is an insulator (cf Figure 11b) with the
gap controlled by (
c
) and a double-pole structure as seen above. Interestingly, this double-
pole structure persists into the superconducting phase, with the gapped poles merging with the
zero-frequency poles at a nite value of
1
2
c
.
(a)
0.00 0.05 0.10 0.15
0.0
0.1
0.2
0.3
0.4
0.5
O
1
(b)
0 10 20 30 40 50
1.5
1.0
0.5
0.0
0.5
1.0
1.5
Im
1.09
0.997
0.701
0.319
0.106
c
Figure 11: (a) A 2
nd
order superuid-insulator quantum phase transition at
c
= 0.163. (b) For
>
c
, we nd an insulating gap. For
1
2
c
< <
c
, we nd a superconducting pole, plus two
isolated poles. For <
1
2
c
, these poles merge with the zero-frequency superconducting pole.
A.3 Varying m
G
and the Gap
In the above we have held m
G
xed. As it turns out, the only eect of varying m
G
is to rescale the
gaps in all of the above (see Figure 12). This ts nicely with naive intuition for the eect of the
compactication radius of the -direction.
15
0 20 40 60 80 100 120 140
1.0
0.5
0.0
0.5
1.0
Im
0 20 40 60 80 100 120 140
1.0
0.5
0.0
0.5
1.0
Im
Figure 12: Varying m
G
rescales the insulating gap. Here, m
G
=
1
4
(left),
1
16
(right).
B Scaling Symmetries
The system described by (2.1) and (2.2) enjoys three distinct scaling symmetries which we can
use to x various parameters to convenient values. Re-introducing G
N
and R
A
in the action and
metric, the three scaling symmetries act as:
scaling symmetry t x
i
r T r
H
ds
2
A
t
A
G
N
q m
1
0 0 0 0 0 0 0 4 0 2 2 0 2 2
2
0 0 0 0 0 0 0 0 2 2 2 4 2 0
3
2 0 1 1 2 1 1 0 0 2 0 0 0 0
(B.1)
The rst two symmetries can be used to x
1
16G
N
= 1 and R
A
= 1. The third, which is the basic
scaling symmetry of the Schrodinger system with dynamical exponent z = 2, can be used to x r
H
to a convenient reference value, r
o
. Given that T =
1
r
3
H
, this xes a relation between T and . To
access more general values of these parameters, corresponding to (r
H
, T
H
,
2
T
,
1