Superconductivity, Superfluidity and Holography
This is a concise review of holographic superconductors and superfluids. We highlight some predictions of the holographic models and the emphasis is given to physical aspects rather than to the technical details, although some references to understand the latter are systematically provided. We include gapped systems in the discussion, motivated by the physics of high-temperature superconductivity. In order to do so we consider a compactified extra dimension (with radius ), or, alternatively, a dilatonic field. The first setup can also be used to model cylindrical superconductors; when these are probed by an axial magnetic field a universal property of holography emerges: while for large (compared to the other scales in the problem) non-local operators are suppressed, leading to the so called Little-Parks periodicity, the opposite limit shows non-local effects, e.g. the uplifting of the Little-Parks periodicity. This difference corresponds in the gravity side to a Hawking-Page phase transition.
Scuola Normale Superiore and INFN, Piazza dei Cavalieri 7, 56126 Pisa, Italy
Departamento de Física Teórica, Universidad Autónoma de Madrid and
Instituto de Física Teórica IFT-UAM/CSIC, Cantoblanco, 28049 Madrid, Spain
1 Introduction: motivations, background material and conclusions
Holographic dualities have become useful tools to explore strongly coupled theories; in particular, in the last few years systems at finite temperature and density of a U(1) charge have been extensively studied. The main motivation comes from condensed matter physics (for a review see ): if the U(1) charge is the electric charge and the U(1) symmetry is spontaneously broken one identifies these systems with superconductors (SCs); when the U(1) symmetry is global, instead, the spontaneously broken phase corresponds to superfluidity.
One of the main ideas behind the application of holography to superconductivity is the hope to shed light on the physics of unconventional SCs (such as high- SCs), which are not completely described by the weakly coupled theory: the Bardeen-Cooper-Schrieffer theory. Cuprate high- SCs, the prototype of unconventional materials, are obtained by doping a Mott insulator and as such they typically show in the ( versus “dopant concentration”) phase diagram an insulator phase close to the region where the U(1) symmetry is spontaneously broken . This calls for a way to describe both conductor/SC and insulator/SC transitions. Since any insulating material eventually conducts if it is probed by a strong enough external current, this in turn requires a method to study gapped phases.
In this article we will briefly review holographic SCs [3, 4] and superfluids (SFs) [5, 4]. The emphasis will be given to physical aspects rather than to the technical details, for which we will refer to the original papers. Given that all SF systems can be considered as SCs in the limit of non-dynamical electromagnetic (EM) gauge fields, for the sake of definiteness, we will use the terminology adopted in the literature of superconductivity, rather than that of superfluidity (unless otherwise stated). We will focus on the best understood version of holographic duality, the anti de Sitter/conformal field theory (AdS/CFT) correspondence. The discussion of the previous paragraph then motivates us to introduce a conformal symmetry breaking in the infrared (IR), which is needed to have a gap. Moreover, since much of the interesting physics of cuprate SCs is layered, we shall restrict our analysis to 2+1 dimensional systems, dual to gravitational theories on four dimensional spacetimes. Although beyond the scope of this review, it is useful to mention that other regions of cuprate phase diagrams require a departure from the theory of Fermi liquids, a challenging property which can be realized in holography .
In the rest of this section we shall introduce the concepts and general properties of superconductivity which are needed to understand the holographic results discussed in the following sections, providing at the same time the reader with the organization of the paper and the conclusions. We assume the minimal field content to describe superconductivity: a charged scalar , responsible for the breaking of the U(1) symmetry and a gauge field . The quantum effective action is a generic functional of gauge invariant operators
where we have introduced the field strength and . The first term in the expression above represents the contribution of local gauge invariant terms. depends instead on non-local gauge invariant objects, which may be required if the material has a non-trivial topology. Although the action in (1) is very general, its form simplifies in some limits. This allows us to extract model independent properties and, at the same time, identify the predictions of specific models, such as the holographic ones that will be discussed here.
One of these limits is the small field limit, in which is approximated by the Ginzburg-Landau (GL) action:
where and are real and positive parameters and we have rescaled the scalar field, , to have a canonically normalized kinetic term for . Also . This theory holds if the system is close enough to the symmetry breaking transition, where is small. Another case in which the action simplifies is the limit of slowly varying fields, such that can be approximated by a two-derivative functional
where and are unspecified functions of . When (1) is a good approximation the SF limit can be taken as : this corresponds to taking the gauge field non-dynamical and thus the corresponding U(1) symmetry global. Deep inside a uniform SC, where is close to the minimum of the potential and the massive gauge field arg is suppressed, (1) reduces to a functional that is quadratic in the fields; this allows us to define the inverse masses of and , respectively and :
Note that and cannot be computed in a model-independent way as their values strongly depend on the form of and . The distinctive property of the superconductive phase is , as opposed to the normal phase in which .
Starting from the general action in (1) it is possible to show in very general terms some of the most famous properties of SCs, such as the Meissner effect, the infinite DC conductivity and the Josephson effect , which have been obtained in holography in Refs. [4, 8, 9, 3, 10].
In this article (unless otherwise stated) we will focus on static configurations, in which everything is constant in time and the temporal component of the gauge field is identified with the chemical potential of the U(1) charge: . Also, in order to study the thermodynamics of these systems we introduce the free energy .
In section 3 we will holographically study the normal phase and explain how to describe conductors and insulators. The homogeneous phase, in which , will be studied in section 4 through the AdS/CFT correspondence. As we will see, this provides a rationale for the fact that superconductivity is suppressed at large rather than at small .
However, inhomogeneous configurations are also of great importance in SCs; e.g. they generically occur when the system is probed by an external EM field. Among the most famous properties of SCs we can certainly include the existence of vortex solutions, which are indeed strongly inhomogeneous configurations. The simplest example of such configurations are straight vortex lines, which can be described by an ansatz of the form and and the other components of the vector potential set to zero. Here are the usual polar coordinates which parametrize the Euclidean two dimensional space. By using the action in (1) one obtains the following large field behavior 
for SCs, while for SFs, where the magnetic field is frozen to an external constant value, , the condensate approaches the homogeneous value with a power law . In Eq. (1) and are unspecified constants, while the action in (1) implies that and . However, the exponential behavior in Eq. (1) tells us that higher derivative terms cannot be neglected as they are generically of the same order of the two-derivative terms. Fortunately, the only effect of the higher derivatives is to modify the values of and , such that we generically have and , as discussed in . and are respectively called penetration depth and coherence length and are important to characterize SCs. Although the field behavior in (1) is model independent, and , like and , can only be fixed once the model is specified. Thus their actual values are predictions of the particular model one considers; in section 5 we will discuss how to extract them from holography.
While vortex solutions exist for any SC, it is not always the case that there is a range of the external magnetic field such that the vortex phase is energetically favorable; when this is true the SC is of Type II and such range is denoted with (if such range does not exist we have instead a Type I SC). is the value of the external magnetic field above which the system is always in the normal phase. The vortex configuration for Type II SCs and slightly smaller than is known to be a triangular lattice of vortices independently on the specific model one considers: for those high magnetic fields the condensate is small and the GL theory can be applied to predict that configuration . is instead the value of below which the system is in the homogeneous superconducting phase. It can be computed through the model independent formula\@footnotemark\@footnotetextHere we use the normalization of such that the external current that generates through is coupled to the gauge field by means of the interaction term in the Lagrangian. (see e.g. )
where and are the free energies for the vortex and superconducting phase respectively. We note that the actual value of is a prediction of the specific model one considers because , and are model dependent. However, in the SF limit, , we always have ; in other words is non-trivial only if the magnetic field is dynamical.
In section 5 we will study the vortex phase in holography and illustrate how to compute the critical magnetic fields and show that the holographic SC is of type II. We will emphasize in particular the differences between the case in which there are no sources of conformal symmetry breaking, other than and , and the case in which conformal symmetry is broken.
Up to now we have not discussed the role of the non-local terms in Eq. (1). For reasons which will be clear soon, we would like to spend the rest of this introductory section to describe the simplest setup in which is relevant: cylindrical SCs threaded by an external magnetic field along the symmetry axis of the cylinder. To render the discussion even simpler we will take the deep Type II limit, , in which the (total) magnetic field coincides with the external one. Since there is a non-contractible loop on this geometry we can construct non-local gauge invariant objects:
where is the charge of the fundamental charge carriers (say the electrons for real world SCs) and the integrals are performed along the non-contractible loop. We will refer to and as the Wilson line and the fluxoid respectively and the dots represent other non-local objects. While the classical action does not depend on , quantum effects, such as the Aharonov-Bohm one (or Sagnac one [12, 13] in the SF literature), could introduce a dependence of on these quantities. A simple (but still general) way to visualize a dependence of this sort is to think that the coefficients of the local part of vary with : for example, in the domain of validity of the GL theory in (1), this corresponds to thinking of and as functions of . Such dependence ought to be suppressed in the classical limit and/or when the typical scale of the non-contractible loops, , is large compared to the other scales in the problem (such as and ) and therefore we should stay away from these limits to see any interest effect. Moreover, since quantum corrections are small in a weakly coupled theory, the biggest effects of the non-local quantities are expected in strongly coupled theories. This justifies the use of holography in this setup (see section 6). Before moving to holography, however, let us identify model-independent effects.
We shall consider the simplest case in which the magnetic field is constant so that it can be represented by a constant vector potential along , which parametrizes the compact spatial dimension, . Since everything is static and homogeneous an appropriate ansatz is
Inserting this ansatz in the quantum effective action, Eq. (1), we obtain
in other words, modulo non-local terms, the system must be periodic in the magnetic flux\@footnotemark\@footnotetextThe factor has been introduced in order to have a canonically normalized kinetic term for in Eq. (1). with period , because can be compensated by a unit shift of the integer . This is known as the Little-Parks effect . This phenomenon has been observed in experiments, which give , and thus it is considered as an evidence for Cooper pairing. If the non-local terms in (1) are non-negligible the Little-Parks does not generically occur; indeed nothing forces to be a function of the local combination only, but it may depend on and separately. The conclusion is that if the system is strongly coupled, far away from the classical limit and is small enough we could see an uplifting of the Little-Parks periodicity to an enhanced periodicity, set by the fundamental charge: . This would corresponds to resolving the internal structure of the composite condensing operator. Whether this really occurs is not model-independent.
2 The holographic model
The first step to define a holographic model based on the standard AdS/CFT correspondence is to introduce fields living on an asymptotically AdS space. The minimal field content required to describe SCs holographically is a charged scalar field dual to a condensing operator () and a gauge field dual to the electric current operator . As we have stated in the introduction, another ingredient we are interested in is a mechanism for conformal symmetry breaking in the IR; this can be achieved by introducing a real scalar (a dilaton) which acquires a non-trivial VEV. The class of actions we consider is 
where is the Ricci scalar, the Newton constant, the dilaton potential, the U(1) gauge coupling and is the radius of AdS. Also, and are generic functions of the dilaton, which do not fulfill particular properties, besides the fact that they never vanish (in order for the semiclassical approximation in the bulk to be justified). Also, .
We are interested in describing a static system with two dimensional rotation and translation invariance, for which the metric and the dilaton have the form
where is the two dimensional Euclidean metric. The temperature can be introduced by requiring the presence of a black hole; so there is a value of the holographic coordinate, , such that . Then , where the prime denotes the derivative with respect to .
Another physical situation that can be described with this setup is the case of multiply connected SCs, more precisely of cylindrical geometry: this can be simply achieved by compactifying one of the two -coordinates: . Then there is necessarily another set of solutions : those obtained from the black holes by exchanging the Euclidean time with ; these configurations have no event horizons and are commonly called AdS solitons. One reason the compactification of is interesting is because this procedure allows us to introduce a scale (other than the and ) and thus breaks conformal invariance even in the absence of the dilaton. For this reason from now on the case in which is compact will be discussed with a simplified field content, where is removed and . The AdS soliton (black hole) geometry is energetically favorable with respect to the black hole (AdS soliton) at sufficiently small (large) temperatures, (). Another reason why we are interested in compactifying a spatial coordinate is to study what happens when cylindrical SCs are threaded by magnetic fields: in these physical setups indeed a universal prediction of holography emerges, as it will be discussed in section 6 and anticipated in the introduction.
Although we include the possibility of conformal symmetry breaking in the IR, in this article we will always assume that the ultraviolet (UV) is conformally invariant, or in other words that the fields are approximated by an AdS configuration close to a value of , say : for close to zero, , and . Thus we can use the standard AdS/CFT dictionary, which relates the properties of a gravitational theory with those of a CFT. In particular the values of and are the sources of and in the CFT. Therefore, if one solves the bulk field equations with boundary conditions
the Green’s function for and are given by differentiating the on-shell action with respect to and ; e.g. the vacuum expectation values are given by
In order to keep the discussion as simple as possible we assume that and do not backreact on the geometry and the dilaton; this can be consistently achieved by taking the limit . A generalization to the case can be found in [17, 18]. A general form of asymptotically AdS dilatonic black hole solutions with the assumed symmetries was derived in . One important aspect of these black holes is that the dilaton is not generically constant but runs with the energy:
where is a real parameter of the potential such that and recovers the Schwarzschild black hole (S-BH) of Einstein’s gravity. Eq. (2) tells us that although the theory has a UV fixed point, , conformal invariance is broken and maximally in the IR, . In order for the configurations of Ref.  to be a solution, has to be appropriately chosen. In particular the requirement of an asymptotically AdS configuration space with cosmological constant implies that at the conformal point (conventionally ) where we have . As a result, for each value of , even , the S-BH with cosmological constant is always a solution. One finds  that the S-BH is favorable at high temperatures, while the dilaton-BH, Eq. (2), dominates in the low temperature region.
Having described the solutions of the dilatonic gravity systems we now want to understand the type of materials they correspond to. Some information can be gained by studying the conductivity, . This will also shed light on the nature of the SF phase transitions, which will be discussed in the next section.
To compute let us consider, on top of these geometries, a small plane wave along a spatial coordinate ,
which is induced by a small electromagnetic field, . Here and are real functions of . The system responds creating a current which is linear in the electric field : that is . Using the AdS/CFT dictionary, the first equation in (2), we have
When an event horizon is present regularity of the solution implies that the plane wave should be ingoing rather than outgoing from the horizon. It can be shown that this results in a non-vanishing Re and DC conductivity 
For an AdS soliton, which has no event horizon, regularity allows for a vanishing DC conductivity, which has been identified with an insulating behavior [20, 18]. Therefore, compactifying a spatial dimension allows us to realize a conductor/insulator transition as the temperature is lowered. Eq. (3), however, tells us that there is another way to suppress the DC conductivity and obtain such transition if a dilaton is present in the spectrum: if one chooses sufficiently small for large values of , the running in Eq. (2) implies that Re is small, especially in the low temperature limit where is large. In this setup the transition (discussed in the previous section) between the S-BH and dilaton-BH, which is obtained by lowering the temperature, represents a conductor/insulator transition. It is worth noting that insulating systems corresponds to solids in the fluid mechanical interpretation .
4 Superfluid phase transition
We now move on and study the simplest example of superconducting state: a static, homogeneous and isotropic material in which the U(1) symmetry is spontaneously broken. The ansatz is
The reader may wonder why we have introduced the temporal component of the gauge field, because having a non-vanishing seems already enough to describe the state we are interested in. The reason is that no regular solution is found when in (4), as shown in , and we want to exclude singular solutions. Physically this corresponds to the fact that there is only one parameter breaking conformal symmetry ( for the black hole and for the AdS soliton) in the absence of the dilaton\@footnotemark\@footnotetextThe dilaton cannot change this conclusion as it is regular for any and the function , Eq. (2), is assumed to be regular and non-zero., thus no phase transition can occur in this case. Therefore we set the UV boundary conditions
The first condition guarantees that the U(1) symmetry is spontaneously broken, while the second one, with a non-vanishing chemical potential , keeps the profile of different from zero and allows us to find a regular solution.
An immediate consequence of this argument is that when the temperature is large compared to the chemical potential a black hole suppresses the superconductivity. At small temperatures instead the system turns into a SC . In AdS/CFT there is therefore a reason for the fact that a metal becomes a SC at small rather than at large temperatures. For the AdS soliton the same conclusion can be reached, but with substituted by .
According to the results of the previous section these SF phase transitions can be either of the conductor/SC or insulator/SC type, depending on the behavior of the conductivity in the normal phase: the two cases correspond to a non-vanishing or vanishing DC conductivity respectively. In the fluid mechanical interpretation, the SF phase of a holographic system with a solid normal phase has been identified with a supersolid .
5 Dynamical gauge fields in AdS/CFT and superconductivity
In the holographic results we have discussed so far the dynamics of the EM field is not important. Therefore, according to the discussion presented in the introduction, they can be applied equally well to SFs and SCs. In this section instead we discuss some important effects of superconductivity, which crucially rely on the dynamics of the EM field.
However, we note that imposing the first boundary condition in (2) treats as an external source, which does not participate in the dynamics of the system. In order to have a dynamical we should integrate over all possible field configurations: (working in the Euclidean space)
where is the Euclidean bulk action computed on a solution of the field equations with boundary condition in (2). Also, we have introduced for generality a kinetic term for and an external current . Then defines as usual the generating functional for the Green functions of . Of course, if there are other operators in the theory, besides the gauge field, in Eq. (5) will depend on the corresponding sources as well; for example in the case of superconductivity discussed in the previous sections also depends on , which we introduced in (2).
In the semiclassical limit this procedure reduces to solving the Maxwell equations,
where we have used . The semiclassical limit corresponds to taking a large external current, and a limit on the parameters of the bulk theory such that becomes large; for example, for the theory defined in (2) this limit is and .
The lesson is therefore that switching from the Dirichlet boundary condition (2) to the condition above promotes to a dynamical field.
As far as the holographic SC model of section 2 is concerned, this procedure has been applied to find the Meissner effect and genuine SC vortices [4, 9], at least in the simplest case of a straight vortex line:
are the usual polar coordinate parametrizing the two dimensional Euclidean space and so here we assume that there are at least two non-compact dimensions in the CFT. This is not the case for the four dimensional AdS soliton; however, vortex solutions have been found on top of the five dimensional AdS soliton in . Using the ansatz in (5) the Neumann-like boundary condition in (5) becomes
while the requirement of spontaneous symmetry breaking and the presence of a finite charge density again fixes the other boundary conditions, Eq. (4). Regularity of the solutions instead fixes the conditions at and at the center of the vortex and physical conditions on the behavior at infinity, , imposes constraints on the remaining boundary.
This results in genuine SC profiles for the total magnetic field : for example for one recovers the Meissner effect, while for one observes the exponential damping of far away from the center of the vortex [4, 8, 9], Eq. (1), allowing for a holographic prediction for and . Such properties should be contrasted with the Dirichlet boundary condition for the magnetic field, which in polar coordinates is (with constant), corresponding to SF systems (see e.g. [21, 4, 8, 9, 22]).
Once the vortex solutions are obtained we can compute through Eq. (1). Indeed selecting the vortex and homogeneous superconducting phase we can compute while, integrating over the bulk à la Kaluza-Klein gives :
The fact that the dilaton appears in this formula has an important consequence on the behavior of at low temperatures , as we now explain. In the absence of Eq. (5) implies that as (in this limit the horizon is removed, ) and Eq. (1) forces to vanish. The physical reason why this emerges is because the theory without the dilaton is scale invariant and so should go to zero as . The presence of the dilaton breaks scale invariance and can avoid this conclusion, providing a non-vanishing : remains non-zero all the way down to if goes to zero fast enough as , a limit that always occurs close to the horizon as , Eq. (2). Interestingly, when remains non-zero at zero temperature, Eq. (3) tells us that the corresponding normal phase must have a suppressed DC conductivity in the low temperature region, resembling an insulator.
6 Multiply connected holographic superconductors and superfluids
As explained in section 1 an interesting setup to apply the holographic techniques are multiply connected SCs probed by magnetic fields. We shall consider the simplest case described at the end of section 1. In the presence of a compactified dimension scale invariance is broken and we will assume for simplicity that the bulk field content does not include the dilaton. Also, the ansatz in (1) corresponds holographically to
with UV boundary conditions . As already discussed, the system is in the S-BH phase at large radii and in the AdS soliton otherwise and, as shown in , such transition corresponds to an uplifting of the Little-Parks periodicity. This is a universal prediction of holography as it uniquely relies on the geometrical properties of the S-BH and the AdS soliton. The simplest way to understand this point is to look at the IR boundary condition (at ) for and , which regularity of the bulk solutions requires:
for the S-BH and
for the AdS soliton. While the conditions in (6) only depend on the local combination of and , that is , the conditions in (6) depend separately on the fluxoid and the Wilson line. This means that the non-local terms in the quantum effective action are suppressed (unsuppressed) in the S-BH (AdS soliton) phase. Correspondingly, according to the model-independent discussion of section 1, the system should be characterized by the Little-Parks periodicity for the S BH and an uplifting of such periodicity for the AdS soliton as shown in [23, 8].
We would like to thank Oriol Domènech, Marc Montull, Alex Pomarol, Oriol Pujolàs and Pedro J. Silva for collaborations and Daniele Dorigoni and Massimo Mannarelli for useful discussions. This work was partly supported by the EU ITN “Unification in the LHC Era”, contract PITN-GA-2009-237920 (UNILHC) and by MIUR under contract 2006022501.
-  Hartnoll S A 2009 Class. Quant. Grav. 26 224002 (Preprint arXiv:0903.3246 [hep-th]). Herzog C P 2009 J. Phys. A 42 343001 (Preprint arXiv:0904.1975 [hep-th]). McGreevy J 2010 Adv. High Energy Phys 723105 (Preprint arXiv:0909.0518 [hep-th]). Hartnoll S A 2009 Preprint arXiv:0909.3553 [cond-mat.str-el]. Sachdev S 2010 Preprint arXiv:1002.2947 [hep-th]. Pires A S T 2010 Preprint arXiv:1006.5838 [cond-mat.str-el]. Horowitz G T 2011 Class. Quant. Grav. 28 114008 (Preprint arXiv:1010.2784 [gr-qc]). Hartnoll S A 2011 Preprint arXiv:1106.4324 [hep-th].
-  Lee P A, Nagaosa N and Wen X 2006 Rev. Mod. Phys. 78 1785.
-  Hartnoll S A, Herzog C P and Horowitz G T 2008 Phys. Rev. Lett. 101 031601 (Preprint arXiv:0803.3295 [hep-th]).
-  Domènech O, Montull M, Pomarol A, Salvio A and Silva P J 2010 JHEP 1008 033 (Preprint arXiv:1005.1776 [hep-th]).
-  Herzog C P, Kovtun P K and Son D T 2009 Phys. Rev. D 79 066002 (Preprint arXiv:0809.4870 [hep-th]).
-  Liu H, McGreevy J and Vegh D 2011 Phys. Rev. D 83 065029 (Preprint arXiv:0903.2477 [hep-th]). Hartnoll S A, Polchinski J, Silverstein E and Tong D 2010 JHEP 1004 120 (Preprint [arXiv:0912.1061 [hep-th]). Pal S S 2012 Preprint arXiv:1209.3559 [hep-th].
-  Weinberg S 1986 Prog. Theor. Phys. Suppl. 86 43.
-  Montull M, Pujolàs O, Salvio A and Silva P J 2012 JHEP 1204 135 (Preprint arXiv:1202.0006 [hep-th]).
-  Salvio A 2012 JHEP 1209 134 (Preprint arXiv:1207.3800 [hep-th]).
-  Horowitz G T, Santos J E and Way B 2011 Phys. Rev. Lett. 106 221601 (Preprint arXiv:1101.3326 [hep-th]).
-  Kleiner W H, Roth L M and Autler S H 1964 Phys. Rev. A 133 1226.
-  Anandan J 1981 Phys. Rev. Lett. 47, 7.
-  Sato Y and Packard R E 2012 Rep. Prog. Phys. 75 016401.
-  Little W A, Parks R D 1962 Phys. Rev. Lett. 9, 9.
-  Hawking S W, Page D N 1983 Commun. Math. Phys. 87, 577.
-  Witten E 1998 Adv. Theor. Math. Phys. 2 505-532 (Preprint hep-th/9803131). Horowitz G T and Myers R C 1998 Phys. Rev. D59 026005 (Preprint hep-th/9808079).
-  Hartnoll S A, Herzog C P and Horowitz G T 2008 JHEP 0812 015 (Preprint arXiv:0810.1563 [hep-th]).
-  Horowitz G T and Way B 2010 JHEP 1011 011 (Preprint arXiv:1007.3714 [hep-th]).
-  Anabalon A 2012 JHEP 1206 127 (Preprint arXiv:1204.2720 [hep-th]).
-  Nishioka T, Ryu S and Takayanagi T 2010 JHEP 1003 131 (Preprint arXiv:0911.0962 [hep-th]).
-  Montull M, Pomarol A and Silva P J 2009 Phys. Rev. Lett. 103 091601 (Preprint arXiv:0906.2396 [hep-th]). Keranen V, Keski-Vakkuri E, Nowling S and Yogendran K P 2010 Phys. Rev. D 81 126012 (Preprint arXiv:0912.4280 [hep-th]).
-  Roychowdhury D 2012 Phys. Rev. D 86 106009 (Preprint arXiv:1211.0904 [hep-th]). Roychowdhury D 2013 Phys. Lett. B 718 1089 (Preprint arXiv:1211.1612 [hep-th]).
-  Montull M, Pujolàs O, Salvio A and Silva P J 2011 Phys. Rev. Lett. 107 181601 (Preprint arXiv:1105.5392 [hep-th]).