#### Abstract

In this paper we propose a generalised holographic framework to describe superconductors. We first unify the description of s-, p-, and d-wave superconductors in a way that can be easily promoted to higher spin. Using a semianalytical procedure to compute the superconductor properties, we are able to further generalise the geometric description of the hologram beyond the AdS-Schwarzschild Black Hole paradigm and propose a set of higher-dimensional metrics which exhibit the same universal behaviour. We then apply this generalised description to study the properties of the condensate and the scaling of the critical temperature with the parameters of the higher-dimensional theory, which allows us to reproduce existing results in the literature and extend them to include a possible description of the newly observed f-wave superconducting systems.

#### 1. Introduction

In holographic approaches to superconductivity, the critical behaviour is encapsulated by the dynamics of a dual theory, where interacting fields propagate in the bulk of a fictitious compactified extra dimension and whose geometry features dictate the evolution of the system with temperature and its reaction to external electromagnetic fields. Different localisations into the extra dimensions correspond then to snapshots of the target theory at different energies. This picture is based on the AdS/CFT correspondence [1–3], the idea that certain strongly coupled theories can be described by an extra dimensional dual theory which contains weak gravitational interactions.

Based on this framework, models for s-, p-, and d-wave superconductivity have been built in the literature [4–10], which provide a computational framework to describe superconducting properties of certain materials.

Despite describing a universal phenomenon, criticality, all these models are based on different setups and link to different target theories.

The aim of this article is to provide a systematic approach to holographic superconductors, which would reveal universal attributes, as well as to put to test some of the features that are commonly introduced in the models. Moreover, this approach could provide insights into models for superconductivity that might be generalised by higher spin configurations. We will indeed propose a formulation of holographic superconductivity valid for f-wave, g-wave, or higher state superconductor. We will further generalise this framework to introduce extra dimensional duals which go beyond the AdS-Schwarzschild Black Hole paradigm. These generalisations will be possible in a semianalytical framework based on matching asymptotic solutions of the differential equations ruling the system.

For this matching approach to work, backreactions of the matter fields on the metric are neglected. As it is discussed in Ref. [5], this limit does not qualitatively change the behaviour of the condensate at .

#### 2. Building a Holographic Superconductor

##### 2.1. The Generalised Action

Holographic superconductivity can be realised in an extra dimensional setup by enforcing a field to condense about a Black Hole (BH) horizon in an AdS background.

To emulate the electromagnetic properties of the material, the matter field—in the following the rank tensor —is endowed with a charge under some gauge symmetry associated with electromagnetism. Then, both and the gauge boson are minimally coupled to gravity; see Figure 1 for a schematic depiction. Further, nonminimal couplings to the metric might as well be included for the purpose of stabilising the theory, in a way that will be shortly discussed.

Superconducting theories have been formulated in both 2 and 3 space dimensions; in the former case, the materials are layered. The generalised action that is introduced in this section allows to describe -dimensional superconductors for arbitrary . However, when singling out particular ansatz for the fields, the calculations are performed for layered 2-dimensional superconductors for convenience, where is then .

On describing superconductors in spacetime dimensions, the +1-dimensional dual gravitational theory will in general look different depending on the particular realisation. But even so, there are a number of common features in every model which this paper aims to point out. The central claim of this work is that the essential dynamical aspects giving rise to a critical behaviour in holographic superconductor theories are captured by a generic action of the form

where is the gravitational strength, is a length scale that coincides with the AdS radius in a pure AdS metric, and the metric signature is . The minimal contributions to the Lagrangian read where is a Lorentz index, , , and is the covariant derivative in the curved spacetime defined by the metric . Meanwhile, the extra piece is introduced as containing all the extra nonminimal coupling terms between the fields and gravity that might be present in each particular realisation. As it has been discussed [9, 10], though not at all essential for the formation of the condensate, the presence of such terms is in general important in safeguarding the stability of the higher-dimensional theory after the introduction of couplings between the matter field, gravity, and electromagnetism. More specifically, said stability would require the insertion of an extra set of properly weighted dimensional terms, a subject set beyond the scope of this paper; see Refs. [10, 11] for a deeper insight into the issue for the case .

A similar approach to Equation (2) has been pursued in Ref. [12] to describe colour superconductivity, i.e., superconductors in the colour pairing of quarks in dense matter QCD.

Upon singling out a specific model, the geometrical properties of the superconductivity carriers in the dual quantum field theory (QFT) are determined by the rank of the tensor field . In particular, this incorporates the s-, p-, and d-wave models present in the literature. together with mixtures of these configurations. Meanwhile, descriptions of higher geometrical complexities (e.g., f-wave with ) are still lacking. However, also for higher spins , the corresponding fields can be described by a rank- tensor with appropriate constraints as we will discuss in Section 4.

Though Equation (2) provides a convenient framework for the study of holographic superconductors in all generality, there are two unsettling aspects of such an action that should be pointed out. Their discussion will be postponed to subsequent sections, but a quick overview of these issues should be of use at this stage already: (1)The kinetic part of the minimal Lagrangian in Equation (2) contains not only the physical degrees of freedom but additional unphysical ones that should be removed from the theory for consistency. This is to say, the Lagrangian does not describe by itself the correct number of propagating degrees of freedom(2)In curved spacetimes, Equation (2) deviates for from known actions for the regarding higher spin fields because it does not account for curvature corrections. This issue will be explained thoroughly in Section 4, where corrected models are presented

As a result, the theory defined just by the Lagrangian in Equation (2) ought to be distinguished from a full higher-dimensional theory that does not suffer from the mentioned problems, for which both and terms in Equation (1) are considered. The latter kind of theories will thus from now on be referred to as the “physical” ones, whereas the theory that is defined by Equation (2) and bears the mentioned inconsistencies will receive the label of “naïve.” Note, though, that physicality here refers to the higher-dimensional theory we build as a hologram of the *real* superconductor system, the target theory. The lack of some features in this higher-dimensional space does not imply an inherent flaw in the target theory, the actual system we want to describe.

Nevertheless, two questions immediately arise from this distinction:
(1)Of what use is the *naïve* theory, then, if it does not comply with its consistency conditions?(2)Is there a feasible way to restore physicality?

In Section 4, the equations of motion (EOMs) for and 2 for the *naïve* theory will be given explicitly, and the inconsistencies in the gravity theory are discussed. Then, it will be shown that the EOMs for some known *physical* theories eventually contain no additional qualitatively relevant contributions when compared to those obtained by using the *naïve* action.

Therefore, an answer can be advanced for the first of these questions: it is expected that the *naïve* model be reliable at the phenomenological level even if there are inconsistencies from the point of view of the gravity theory. Thus, a study of the formation of the condensate by means of the *naïve* EOMs is justified and will be carried out in Section 3.

The second question is of a rather more technical nature, but as it will be explained in Section 4.2, the nonphysicalities in Equation (2) can be erased as long as the EOMs that are produced by said action are accompanied by appropriate constraints over the fields. This is equivalent to saying that a suitable expression of is given. However, while these constraints are easy to grasp in some limits—e.g., by antisymmetrising the kinetic term for and fixing to be Lorentz-symmetric, transversal, and traceless for in the free limit—the general case for higher spin and interacting fields is more involved and not iteratively generalisable. Already the formulation of the constraint equations for is well beyond the scope of this work. The correct formulation of such constraints is in fact an open problem in the field [11], for which solutions may only exist in a limited number of background metrics. Nevertheless, the issue is worthy of consideration and will be tackled to some extent in Section 4.

##### 2.2. Ansatz for the Fields

Once the theory has been established, particular ansatz must be implemented for the fields. These ansatz are required to be in accordance with the constraints imposed over the fields that were just discussed in the paragraph above. Typically, the gauge field only depends on the holographic coordinate and is parametrised as

Regarding the matter field, the different tensor components are fixed so that a single AdS-radial degree of freedom remains, while the required geometrical properties are reproduced. The particular ansatz adopted in this article are shown in Table 1. Note that in the s-wave () case, spatial isotropy is preserved in the QFT dual, while by picking the component in the p-wave () case generates a preferred direction. Finally, this particular ansatz for the d-wave superconductor () [9] generates a condensate in the plane with translation invariance that breaks the rotational symmetry down to , with the condensate flipping its sign under rotations on that plane.

##### 2.3. The Generalised Warped Metric

It has already been stated that the superconductivity phase transition can be triggered by an AdS BH background metric in the 4-dimensional gravitational theory. Here, this statement will be briefly clarified and, at last, generalisations to the AdS-Schwarzschild metric will be introduced.

The background metric under consideration is time-independent and fully symmetric in the QFT spatial dimensions. Therefore, it only depends on the holographic coordinate and can be parametrised as

In the following, the probe limit is assumed, meaning that the matter fields do not backreact on the metric. In this case, the matching approach that relies on asymptotic expansions of the fields provides analytical solutions for the condensates that will be computed further on.

Although the probe limit is known to be not numerically accurate for , the scaling of the condensate approaching the critical temperature is not affected. This was shown by comparing the probe limit to the numerically obtained fully backreacting solutions, e.g., in [5] for s-wave superconductors. The same applies to the condensate scaling and the dependence of the critical temperature on the charge density for p-wave geometries, as shown in [13–16].

Thus, the probe limit is expected to deliver a valid description of the superconductor’s properties close to the critical temperature. The main effects due to backreactions are expected to change the value of the condensate at , albeit with a lesser impact than the uncertainties arising from the matching point dependence of the semianalytical approach which is discussed in Appendix B.

In particular, this means that the metric components are fixed and asymptotically becomes the metric of anti-de Sitter space, i.e., as .

The holographic correspondence identifies QFTs at finite temperature as dual to an AdS space with an event horizon with radius , whose introduction renders a physical cut-off in the AdS coordinate. Then, the associated Hawking temperature of the horizon is interpreted as the temperature of the dual QFT. Customarily, for an AdS-Schwarzschild BH located at , the event horizon is set at , yielding the metric with which constitutes an exact solution of the gravity EOM inferred from Equation (1) when the presence of the remaining matter fields is neglected.

The holographic coordinate is conventionally redefined as which ranges in the interval to facilitate the calculations. This means that one usually works with

Notice also that the AdS-asymptotic region is located near , while the horizon lies at .

While the former is a well-motivated choice for the metric, it is not the only one that could be made. For this reason, in this paper, other possibilities will be explored with the motivation of continuing the study of holographic superconductors in full generality. This generalisation will be implemented via the ansatz that retains explicitly the essential feature of the AdS-BH metric , i.e., the first-order zero at the horizon with

Notice that recovering flat AdS metric in the asymptotic region constrains the function to fulfil as .

As it will be discussed in Appendix A, different choices for the power of have been studied, but all the cases under consideration failed to reproduce the characteristic temperature behaviour of the condensate.

Therefore, the compact form of the metric in Equation (8) captures the essential features needed to produce a possible holographic model of critical phenomena at finite temperature.

Besides the canonical example of the AdS-Schwarzschild metric, there are various instances of extra dimensional theories which would be good candidates for a holographic superconductor while still neglecting backreaction. In particular, the correct pole structure in Equation (8) is found on the Hirn-Sanz metrics [17, 18], which were employed to develop duals for QCD and for new physics at the electroweak scale, the AdS-dilaton metric [19], a specific type of metric designed to reproduce the Regge trajectories, and the Sakai-Sugimoto metric [20, 21], inspired by Witten’s string theory model of D8 branes [22]. Specific forms of the function for these proposals are summarised below: where is a polynomial of up to . Different choices of would lead to the same universal behaviour but slightly different details. This will be discussed further in Section 3.5.

##### 2.4. Hypothetical Higher-Order Superconductors

Apart from its compact form that summarises models for s-, p-, and d-wave superconductors, the *naïve* model presented in Equation (2) allows to speculate about higher-order superconductors, for instance about the f-wave superconductor with . Indeed, as it will be shown in Section 3, the condensates for the *naïve* models share the temperature behaviour with the physical models for . Therefore, it is expected that the *naïve* action can provide a valid approximation of the condensates for higher spin superconductors.

To give one example, for possible f-wave () superconductors, the matter field is taken to be a rank-3 tensor with the action

The dynamical field is singled out by the ansatz with all remaining components set to zero; compare also Table 1.

Note that the mass term in Equation (11) in principle would break gauge invariance. Nevertheless, note that this *naïve* action should be regarded as an effective description that considers only the relevant degrees of freedom. Some ideas to render the theory consistent will be addressed in Section 4.

##### 2.5. The Critical Temperature and Formation of Scalar Hair

After generalising both the field content and the geometry of the hologram, we briefly discuss the process of formation of a condensate.

In the presence of the background metric and the just introduced matter field with physical components , critical behaviour is triggered by the interplay between electrostatic forces and gravity. The BH attractiveness can be measured in terms of its surface temperature which is interpreted as the temperature of the dual field theory, too.

Once equipped with the notion of a finite temperature in the QFT, a critical temperature is determined by the dynamics of the gravity theory such that for electrostatic repulsion overcomes gravitational attraction and a superconducting layer can condense over the horizon. For it to actually occur, spacetime must be endowed with a negative curvature metric, as it happens in AdS geometry, while in a flat geometry scenario, the superconducting layer would be blown out to infinity.

Then, for , the free energy is minimised by a nontrivial configuration . This is to say, the gauge symmetry is spontaneously broken and develops a nonvanishing vacuum expectation value (VEV), leading to the formation of a condensate of the operator in the dual QFT, . This can be interpreted as the Black Hole developing scalar hair below . Meanwhile, for , the condensate gets swallowed and the system is out of the superconducting phase; the trivial solution is the only viable choice.

The existence of a critical temperature below which the condensate emerges is the first hallmark of the second order transition giving rise to superconductivity.

#### 3. The Condensate

##### 3.1. The Generalised EOMs

Eventually and under the ansatz in Table 1 and Equation (4), the theory contains two dynamical scalar fields, and , besides gravity. Under the assumption of a fixed metric—i.e., neglecting backreactions of these fields upon the metric itself—a set of two EOMs can therefore be produced from the action.

For the class of models defined by Equation (2), the ansatz in Table 1, and the metric given by Equation (5) with described by Equation (8), a closed expression for an arbitrary tensor of rank can be written down for the “minimal” part of the action density, . It is the simplest in terms of the AdS variable and reads where the factor denotes the multiplicity of the dynamical field within the tensor field itself, cf. Table 1.

A convenient reparametrisation, inspired by Ref. [10], consists in writing

This allows to write the EOMs in the compact form below: where a dynamical mass for the field has been defined as

It must be highlighted that these equations are identical to those in Refs. [5, 6, 9, 10], save for the explicitly -dependent term . Thus, the difference between the naïve models and those in the literature can be absorbed into the dynamical mass for .

Taking into account the rescaling of in Equation (13), the EOMs for p-wave superconductors also agree with those in [23].

At this point, the EOMs above are all that is needed to unravel the dynamics of the presented generalised holographic superconductor. Though their exact analytical solution is not known, in the remainder of this section, it will be shown that, by using the semianalytical matching procedure presented in Ref. [24], one can capture the correct behaviour of the condensate in every case. The reliability of this method for the analysis of the equations under considerations ascribes to the absence of important features between (horizon) and (AdS-asymptotic region) and involves the calculation of asymptotic solutions near both regions that are then matched together at an intermediate point.

##### 3.2. Solutions near the Horizon

The construction of semianalytical solutions requires the implementation of boundary conditions near the horizon that enforce regularity of Equations (14) and (15). That being so, the former equation imposes on the gauge field , while cancelling all divergences in the latter requires that

For convenience, the mass is expressed in units of the curvature from now on, with the dimensionless parameter .

Assuming that both and are regular for and keeping terms up to second order in , one can expand the fields around the horizon in the following manner:

for . The second derivatives at are recovered by plugging these expansions into the EOMs. Substituting the results back into the expansion yields the asymptotic values of the fields in terms of the two still unset boundary values and , where the following set of coefficients has been conveniently defined:

##### 3.3. Solutions in the Asymptotic Region

At the AdS boundary (), must vanish and be finite to ensure normalisability of the fields. Therefore, the equations for and asymptotically decouple and an independent expansion for each field can be constructed. In the case of the gauge field, one can easily infer that where and and are identified with the chemical potential and the charge density at the horizon, respectively. Meanwhile, since the metric becomes asymptotically AdS, the matter field behaves as where defines the condensate of the dual operator corresponding to the matter field, denoted by , and is the mass scaling dimension of said field.

The latter can be readily determined to be

Note that the definition of in this work might differ from the convention in other articles. For convenience, the matter field has been rescaled in Equation (13) to give easier matching expressions. But accordingly, the scaling of and with differs by a power of . Having this transformation in mind, Equation (28) also perfectly agrees with Equation (14) in [23] for .

##### 3.4. The Matching Procedure

At last, the matching is realised by imposing the Dirichlet and von Neumann boundary conditions at some intermediate point over each pair of asymptotic solutions. This gives rise to a set of four algebraic equations: whereby the four remaining free parameters , , , and can be uniquely determined.

As demonstrated in Ref. [24], the solutions depend only mildly on the precise choice of provided that it is far from its two limiting values; see Appendix B for more details. That being so, the convenient choice will be adopted in the following.

Defining the critical temperature allows to write the solution for for the generalised model under consideration as where one can write

##### 3.5. Discussion

Equations (31), (32), and (33) are the main results of this paper and can be used to describe phenomenological aspects of holographic superconductors, including the parametric dependence of the condensate and critical temperature.

To begin with, as the equations explicitly show, the scaling of the condensate with temperature is of the form

for , exactly as it is expected from a superconducting system. This is a result which was obtained in Ref. [24] for s-wave superconductors and in an AdS-BH metric, but in this paper, it has been shown that such a dependence emerges also for higher-order superconductors and for more general metrics.

The analytical derivation of this scaling relation with the temperature of the condensate, Equation (34), is worth describing in more detail: this scaling behaviour has been shown to be completely determined by the equations for , Equation (29), and is independent of the equations for , Equation (30). This is to say, to this level, the square root behaviour of the condensate appears solely because of the dynamics of and is seen as a consequence of the equilibrium between gravity and the electromagnetic force (cf. Section 2.5). Qualitatively, one can understand this behaviour by the following distinction: the field that carries the information about the material determines the scale of critical behaviour and the external EM field carries the information on the dependence of the condensate with this temperature. Schematically,

Additionally, Equations (32) and (33) retain their explicit dependence on the parameters of the gravity theory, namely, (i)the mass in units of the curvature radius : (ii)the free parameters of the generalised metric: and

So, they provide detailed analytical information about the formation of the condensate from the gravity side and allow for a study of this parameter dependence. Notice in fact that a well-defined condensate requires so any point in the parameter space that does not produce this required outcome will not give rise to a superconducting behaviour. In particular, notice that rendering and real and positive requires the following conditions to be fulfilled: (i) and , which are automatically satisfied for non-vanishing ; those conditions simply encode the necessity of a coupling between and the gauge field that is not negligible with respect to the gravitational pull of the event horizon(ii) which imposes a relation between and :

Notice also that the positive definiteness of Equation (28) prevents approaching the ill-defined limit [25]. Moreover, also from the mentioned equation, one can infer that the requirement implies a lower bound for the mass: which is not but the Breitenlohner-Freedman (BF) bound [25] for this theory.

Figure 2 shows the dependence of the different condensates with the temperature in units of as described by Equation (32) and puts them in comparison with the condensates obtained from some of the models in literature.

Choosing the masses for all models at their respective BF bound, Equation (38), the anomalous dimension is in all cases (cf. Equation (28)). Then, for near , the condensates in Equation (32) can be approximated by

After normalising by , the condensates have the same dimension for all values of and differ only by the coefficient .

The right panel focuses on the scaling of either curve near the critical temperature in a double-logarithmic form. From Equation (39), one reads of the slope in the double-logarithmic form being 1/2 while the offset is shifted by .

For further comparison between the *naïve* and literature models, see Section 4.3. Moreover, see Section 4 for an outline of the essential features of the latter models together with some of their technical aspects.

It remains to look at some of the properties and parameter dependencies of the critical temperature as defined in Equation (31). Figure 3 exhibits the scaling of with and for s-, p-, d-, and f-wave superconductors (from left to right), under the assumption that as in the Schwarzschild BH case . This is a simplifying assumption that is required in order to reduce the amount of parameters, but the result should not differ much from the general case.

For values of sufficiently far away from 0, demonstrates a simple scaling with and is well defined by its only pole at all the way down to the BF bound. However, for approaching zero, it displays a more interesting behaviour. In this regime, the leading contributions to the critical temperature in Equation (31) are given by the terms within and . Correspondingly, it simplifies to

Therefore, the critical temperature visibly diverges both when and . This behaviour can be spotted in Figure 3, where the lighter lines corresponding to smaller values of diverge long before reaching the BF-bound.

It is worth noticing the different parametric dependence of and for the typical AdS-BH case and other proposals in Equation (10); see Table 2. Note that is typically negative.

So far, all temperatures have been presented in units of the charge density . However, the distinctive quantity to relate the models to existing or hypothetical physical systems is the critical temperature in Kelvin. Therefore, the charge density is fixed such that the critical temperature of the s-wave superconductor gives 77 K at the BF bound and for the parameter choices and as in the Schwarzschild BH case.

The critical temperature in Kelvin for some choices of obtained from Equation (31) is shown in Figure 4 as a function of . Note that for s-, p-, and g-wave superconductors, a finite value of the critical temperature can be defined in the Schwarzschild BH limit , while it diverges before reaching this limit for d- and f-wave superconductors.

The semianalytical matching approach pursued in this article requires a mild dependence on the matching point . As it will be discussed in Appendix B, this is the case for close to the BF bound. However, the curves in Figure 4 imply that the AdS Schwarzschild BH background metric does not provide a critical temperature at the BF bound for d- and f-wave superconductors. Accordingly, this fact suggests to take into account generalised background metrics as it has been motivated in Section 2.3.

Apart from their condensate, holographic superconductors are typically characterised by the optical conductivity. However, the analytical matching approach presented here breaks down for the conductivity, see Appendix C for details.

#### 4. More on Higher Spin Fields

##### 4.1. Naïve Versus Physical Models

The *naïve* model introduced in Equations (1) and (2) and the thereupon derived EOM for , Equation (15), bear the two main shortcomings that were already mentioned in Section 2.1. First and foremost, the kinetic part of Equation (2) is inconsistent with the axioms of quantum field theory in curved spacetime as it leads to unphysical particle states and it does not describe the correct number of propagating degrees of freedom. Secondly, Equation (2) does not exactly reproduce the known models for p- and d-wave superconductors and only the s-wave case matches the models in the literature, despite producing a similar superconducting phenomenology.

It is the latter problem that can be demonstrated the easiest to tackle and will be the focus of attention in the following section.

It should then be illustrative to enumerate the differences between the generalised, *naïve* model defined by Equation (2) in the noninteracting limit and the well-established Lagrangians that describe free tensor fields of ranks and for different spacetimes. By doing so, the aims of the writers are two-fold: this discussion shall contribute to explaining where the differences originate from and will hopefully cast some light on some of the issues that could arise for higher spin fields.

Rank- tensors are suitable to describe the matter part of a p-wave holographic superconductor, and the right description of the dynamics of such fields in the free limit for Minkowski spacetimes is known to be given by the Proca Lagrangian [26].
whose mass term would break gauge invariance. Note again that this action should be regarded as an effective action, containing only the relevant degrees of freedom to describe superconductivity. However, for p-wave superconductors in flat spacetime and with , the *naïve* model reduces to the matter Lagrangian
that differs from Equation (41) by the lack of antisymmetrisation of the kinetic term and leads to the Klein-Gordon like EOM
where is the 4-dimensional D’Alembertian operator.

The imprint that the mentioned difference between Equation (41) and Equation (42) leaves on the dynamics of the theory depends crucially on the considered background spacetime. Note that in the so far considered Minkowski spacetime, the Proca equation inferred from the Lagrangian in Equation (41) reads
which simplifies in the Lorentz gauge to the EOM obtained from the *naïve* action, Equation (43). This is to say, in Minkowski spacetime, the lack of antisymmetrisation does not change the dynamics.

Next, one can see with little effort that turning on the minimal couplings and promoting the partial derivatives to gauge-covariant ones, , does not introduce any difference either. This means, minimally coupling the fields in Minkowski spacetime does not introduce any deviation in the dynamics derived from the naïve action with respect to the physical action.

In contrast, covariantising the partial derivatives in curved spacetimes, the second term in Equation (44) no longer vanishes in the Lorentz gauge and introduces additional terms with respect to the naïve EOM. To see this explicitly, one can write the Proca equation for in the presence of the generalised background metric and switch on the minimal couplings to the gauge field, obtaining

that coincides with the SU(2) Yang-Mills model for a p-wave superconductor [6]. Meanwhile, the EOM obtained from Equation (42) is (notice that trading for via Equation (13), Equation (46) becomes Equation (15) for )

which contains two further contributions to the 0^{th} order term in that originate from derivatives of the metric components and would not be present in flat space. However, it is noteworthy that these two terms enter the equation at the level of the mass term. Consequently, they only affect the definition of the critical temperature but not the shape of the condensate and do not crucially affect the dynamics of the theory after their presence is accounted for; see Section 4.3 for a proof of this fact. As promised, they hence do not alter the phenomenology of the model.

In d-wave holographic superconductors, the matter field is described by a rank-2 tensor whose exact dynamics in the free limit and for a Minkowski background are now given by the Fierz-Pauli Lagrangian [27] where the trace is defined using the Minkowski metric . The naïve model for vanishing couplings in Minkowski spacetime yields a rather more succinct expression for :

Now, this Lagrangian turns out to coincide exactly with Equation (47) in the traceless and transversal gauge given by

Accordingly, the naïve EOM equals the Fierz-Pauli equation in this gauge, in which they read

Thus, once more, both the physical and the naïve theories turn out to display the same dynamics.

Yet again, difficulties arise when considering curved spacetimes that are even more problematic to address than in the former case. Remarkably and as it is discussed in Ref. [11], adding gauge couplings can also introduce instabilities in the theory. However, although this is true on a technical level, it will be shown that the EOMs obtained from the naïve action and the minimally coupled extension of the Fierz-Pauli action still agree for the chosen ansatz for .

In curved spacetimes, the case is indeed more difficult than that of the p-wave superconductor. In particular, covariantising the Fierz-Pauli equation and introducing minimal couplings does not suffice to produce the required number of independent constraint equations [10]. Consequently, spurious propagating degrees of freedom are left in the theory that lead to, in principle, disastrous inconsistencies such as the emergence of ghosts or the loss of hyperbolicity.

Several approaches exist to cure this kind of difficulties, one of which is built in the BHRY model presented by Ref. [10]. As will be shown in Section 4.2, this model is constructed by including in the action all the possible terms of scaling dimension that are quadratic in the spin-2 field, weighted by coefficients that are treated as Lagrange multipliers. By fixing these coefficients, one ensures the existence of a correct number of independent constraint equations that retain the physical degrees of freedom only. This process leads to the BHRY Lagrangian that contains a number of nonminimal couplings between the spin-2 field, the gauge field and gravity, and a more involved kinetic sector:

However, in spite of its rather involved form, this Lagrangian simplifies in the traceless and transversal gauge for , in which it reads
i.e., it only differs from the *naïve* action by the curvature correction term involving the Riemann tensor . The corresponding EOM is

It is more illustrative to compare the equations arising from each model explicitly. Making use of the ansatz for in Table 1, the *naïve* EOM translates into
which coincides with that of the CKMWY model in Ref. [9], while the BHRY equation corrected by the curvature term delivers (the Riemann tensor for the generalised metric has been calculated with the subpackage *xCoba* of the Mathematica package *xAct* [28])

Hence once more, notwithstanding the formal issues of the spin-2 case, the differences between the models reduce to a series of curvature corrections that enter the EOMs at the level of the dynamical mass. Then, and as it happened in the p-wave case, the corrections only change the definition of the critical temperature and should not affect the model phenomenologically.

Notice also that, in terms of , the latter equation takes the simpler form in which the curvature corrections have been absorbed by the redefined field. One consequence is that after such a redefinition of the field, the EOMs of the naïve model are identical to those of the s-wave case; see for instance Ref. [4, 5].

After describing the procedure for p- and d-wave superconductors, it is time to address the problem for a field of rank . Even though the correct theory describing the dynamics of such higher spin fields in the general case is still to be formulated, there are some cases that are well known in the literature with which the model presented in this paper can be compared. In fact, the simplest possible situation of a massless higher spin field in Minkowski spacetime is well-understood, the Fronsdal equation [29] providing the right description of the propagating degrees of freedom in that case. It reads which reduces to in the traceless and transversal gauge

Trivially, this coincides with the EOM arising from Equation (2) in the massless limit .

But again, considering curved spacetimes seemingly puts one in dire straits. The main issue in this case is the noncommutativity of the covariant derivatives, which prevents the second term in Equation (57) from vanishing and introduces curvature correction terms that correspond to the higher spin generalised version of those that have explicitly been found in the p- and d-wave cases.

A generalised procedure exists to construct the required set of constraint equations to cure all instabilities in the rank- case. This method will be outlined in Section 4.2, but it should not be regarded as an universal remedy, as for larger , it becomes involved and there is no systematic generalisation of every of its steps.

Therefore, it would be useful to have an operative method to calculate these curvature corrections without having to explicitly consider constraint equations. In the AdS background, several worthwhile attempts to describe the dynamics of free higher spin fields have been made in the literature [9, 30–32]. One approach displayed in Ref. [33] is to regard the curvature correction terms as belonging to a dynamical mass , yielding an equation of the form in the traceless and transversal gauge, where the mass takes the form in spacetime dimensions.

For the presented superconducting theories, though, the correct curvature corrections are hard to write down in a general form for higher spins. The difficulty of this enterprise arises from two sources: (1)The nonpurely AdS metric(2)The couplings with the electromagnetic field

So as to visualise the first issue, one can consider the EOM for the matter fields in terms of the initial field , i.e., before performing the redefinition of the matter field defined by Equation (13). In those terms, the EOM for general takes the form where the dynamical mass is shifted by a pair of curvature correction terms

and the coefficients and weighting those terms are

Now, in a pure AdS background where , the , , and terms all scale with and can be collected into a single dynamical mass term proportional to

Hence, the curvature corrections become only a shift to an effective mass and are effectively erased from the EOMs. But unfortunately, this does not hold for arbitrary where the derivative of the metric is not proportional to , spoiling the convenient relation from AdS. Meanwhile, the coupling to the gauge field also complicates the argumentation and sabotages some ideas. In fact, one very promising path to present the EOMs in a neat form leads down to exploiting the invariance of an action of the form Equation (2) under gauge transformations for free fields in AdS-like theories, as shown for instance in Ref. [19]. In that paper, it is shown that by electing an appropriate gauge, the fields can be rescaled in such a way that the action only depends on partial derivatives of these fields and the EOMs can directly be read off. But even though this method could in general take care of general metrics by adapting the field rescaling, it is easy to see that it eventually breaks down when including gauge field interactions via gauge-covariant derivatives.

##### 4.2. Constructing Complete Models for Higher Spin Fields

To date, no complete, pathology-free formulations of higher spin theories have been found. However, several methods exist in the literature which suggest that doing so is formally possible, at least for a limited number of spacetime geometries that include AdS. Therefore, after having pointed out the differences between the naïve and physical models, it will be interesting to provide more detail about how it can be done in the general case by putting together a sufficient set of constraint equations.

The procedure is explained in Ref. [11] and it is summarised here for the sake of completeness. The main idea was already anticipated in the examples in the first part of this section. While transversality and tracelessness conditions are enough to eliminate unphysical degrees of freedom in free theories, further constraints are required when interactions are switched on in a higher spin field theory. For consistency, as a first requisite, one needs that the same number of degrees of freedom as in the flat theory is conserved when introducing interactions. Another crucial point is the preservation of causality, i.e., the absence of superluminal motion must also be enforced.

The general algorithm that allows to construct a theory fulfilling those requirements goes at follows: let be the number of degrees of freedom contained in the *naïve* theory that will be denoted as following Ref. [11]. By applying the principle of least action, a set of equations can be obtained from the theory that are not all independent of each other. In particular, the equations can be rewritten and split in two different categories:
(i)A set of equations which fix of the second derivatives of the degrees of freedom(ii)The remaining equations which do not contain accelerations. These are called *primary constraint equations*

As a consequence, the system of equations that directly arises from the action does not completely characterise the dynamics of the theory, and additional constraints must be brought in from elsewhere. One can recur, at this point, to conservation of the primary constraints. By doing so, a new set of *secondary constraint equations* is found, containing second derivatives which may or not fix the still undetermined degrees of freedom.

If they do not, one can then iteratively repeat this procedure until all second derivatives are fixed by the equations. Typically, for rank-*ℓ* tensor fields, a total of ℓ iterations are necessary, which generically require to rewrite the equations in an appropriate form. From this general insight, it becomes clear that the actual calculation of the constraints will be complicated for higher spin models and cannot be generalised automatically to arbitrary .

For an explicit example, the case of the d-wave superconductor will be considered. The problem was in fact solved for this case in Ref. [10] to set up the BHRY superconductor. The model aims to describe a massive spin-2 field in spacetime dimensions with a total of physical degrees of freedom. However, in the field being represented by a symmetric tensor with free elements, the theory contains degrees of freedom that are spurious. Accordingly, extra constraint equations are required, so as to fix both the values of those fields and their momenta.

The most general action containing only terms of scaling that are quadratic in the rank-2 tensor field can be written as where the coefficients are the Lagrange multipliers that were referred to in Section 4 when commenting on the BHRY model.

Let denote the EOM for the tensor component . Then, it is possible to divide the equations in the two categories defined before. As it is shown in Ref. [10], using the usual convention for Greek and Latin Lorentz indices: (i) are dynamical equations which determine the second derivative (ii) do not contain accelerations. Therefore, they are a set of primary constraints

In order to fix all the degrees of freedom, one applies the conservation of the dynamical equations. This amounts to write . Combined with the trace equation , it can be shown that doing so provides an extra set of secondary constraints. That way, only 2 of the required constraints are missing that can be obtained by iteration. In particular, the penultimate constraint originates from the second divergence , combined with the trace equation again. Finally, the last constraint is computed from the time derivative of the first divergence, .

Following these steps, the coefficients to in the most general Lagrangian Equation (65) can be fixed to give Equation (51). This Lagrangian defines a theory for the d-wave holographic superconductor that contains the right number of degrees of freedom. However, as it is discussed once again in Ref. [10], the theory obtained by this method has several limitations. For instance, the constraint equations only hold if the background is fixed, and so it cannot account for backreactions on the metric.

##### 4.3. Comparison of the Condensates between the Naïve and Physical Models

Now that the Abelian-Higgs s-wave, SU(2) Yang-Mills p-wave, and CKMWY and BHRY d-wave model have been introduced, everything is set up to resume the comparison between the *naïve* and physical models which was started in Section 3.5 including Figure 2.

Before comparing the condensates for the *naïve* and the mentioned physical models, it is useful to return for a moment to the case of general seeing as it provides some tools to relate the individual models. As it has been motivated in Section 2.2 and shown explicitly for , a spin- field can be represented by a tensor of rank accompanied by appropriate constraints. The general EOM for with coefficients and that incorporate the curvature corrections implemented in the physical models has already been given in Equations (61) and (62) in Section 4.1 in terms of .

However, in order to compare with the calculation in Section 3, it is convenient to work in terms of , cf. Equation (13). Then, the equation reads

When working in terms of the arbitrary curvature-correction coefficients and , one finds that, by following the same steps as in Sections 3.2 and 3.3, the only differences appear in the coefficients of the expansion near the horizon, which now read

For the *naïve* model, the coefficients are simply given by and . Plugging those values in allows to recover the expressions found in Section 3.2. It is crucial how these coefficients enter into Equations (31)–(33): the critical temperature and the prefactor of the condensate are modified, but the scaling with for that is distinctive of a superconducting condensate remains unaffected.

As it has been discussed in Section 4.1, the action and the EOMs differ between the *naïve* and physical models for p- and d-wave superconductors. Meanwhile, the s-wave models are in agreement with each other. Now, the expressions above allow to compare these differences both at the level of the EOMs and to eventually quantify them. Starting from Equation (32) and using the updated parameters, Equations (67)–(71), the modifications for the physical models can be regarded as and terms (cf. Equation (66)). Those terms enter at the same level as the mass parameter . Accordingly, and as it is apparent from the discussion in previous sections, near the AdS-asymptotic boundary, they can be absorbed into a dynamical mass:

Then, the anomalous dimension reads

Equation (72) implies that, after an appropriate case-wise choice of the value of , the values of and thus for the *naïve* and physical models can easily be matched. So, in order to better compare the different models, the asymptotic dynamical mass Equation (72) can be chosen at the BF bound, , resulting in . The corresponding values of together with the parameters and for each model are presented in Table 3.

Then, looking back at Figure 2, which shows the condensates as a function of the temperature for these values of , the common features as well as the differences are clearly visible. Notice that since the anomalous operator dimension agrees for all models, using the normalisation in the figure, all curves exhibit the same scaling with the temperature. The only differences originate from the coefficient of the condensate defined by Equation (33) that results in an overall multiplicative factor for the curves.

Naturally, due to the scalar nature of the matter field and the simplicity of its description, the *naïve* s-wave model agrees with the Abelian-Higgs s-wave model and so do the respective condensates. Meanwhile, the two curves for *naïve* and SU(2) Yang-Mills p-wave superconductors coincide because, by chance, the coefficient is the same for both models and the changes in due to the different choice of are cancelled equally by in Equation (68).

This is no longer the case for the d-wave superconductor, since and thus agree between the *naïve* and the CKMWY models (which coincide at the level of the action), but they both differ from the BHRY model. Even so, as it has already been stated, the differences appear only in the overall factor of the condensate. By taking the ratio of the coefficients for the *naïve* and BHRY p-wave model,

one finds that, from a phenomenological point of view, the models differ by the order of a 25% in their prediction of the magnitude of the superconducting condensate.

#### 5. Conclusions

In this paper, we presented a generalisation of the holographic approach towards superconductors.

First, a *naïve* model that describes the action of a matter field with spin of arbitrary value , coupled to an electromagnetic field, has been introduced in the bulk. Afterwards, the widely used AdS Schwarzschild Black Hole metric has been generalised to sets of metrics with a first-order zero at the horizon.

Eventually, equations of motion for the relevant degrees of freedom of the theory have been solved by means of the semianalytical matching approach presented in Ref. [24].

One great virtue of the *naïve* model consists in the fact that it can reproduce the expected scaling of the condensate with the temperature, , independently of . Notice that this was not at all straightforward, as the *naïve* model is not a self-consistent field theory in -dimensional curved spacetime for any value of : although the model agrees with the Fronsdal equations in Minkowski spacetime, it does neither account for correction terms in curved spacetime nor describe the correct amount of physical degrees of freedom. It hence needs to be endowed with suitable constraint equations, in order to be promoted to a consistent *physical* field theory.

In spite of this, the results for the condensate and the critical temperature retain the correct dependence on the free parameters of the theory and allow for an analytical study for any value of considered in the literature so far.

In particular, the mentioned model reproduces the Abelian-Higgs model for s-wave superconductors and the CKMWY model for d-wave superconductors studied in Ref. [5] and Ref. [9], respectively. On the other hand, the equations of motion derived in the *naïve* model differ slightly from the SU(2) Yang-Mills model for p-wave superconductors of Ref. [6] and the BHRY model for d-wave superconductors of Ref. [10]. The differences can be regarded as a change in the dynamical mass and only affect the critical temperature’s dependence on the extra dimensional parameters, whereas the scaling with the charge density is universal and leads to the same superconductor phenomenology. Moreover, a direct comparison of the numerical results shows that the coefficient of the condensate computed with the *naïve* d-wave model differs by approximately 25% from the *physical* BHRY d-wave model.

In order to address the important point of the difference in the equations of motion for the relevant degrees of freedom that arises when comparing the *naïve* model and the consistent field theory for a given spin , a technique to compute the constraint equations has been summarised. For the d-wave superconductors, the constraint equations were calculated explicitly to shed some light on how to obtain the *physical* equations of motions from the action in Equation (2). We found that in all cases studied in the literature, the equations of motion derived using the *naïve* model single out the essential features a holographic model for superconductivity must encapsulate.

Our results for s-, p-, and d-wave holographic superconductors encourage us to speculate about the behaviour of possible higher spin superconductors. In recent works [34], the existence of f-wave superconductivity was conjectured on systems like Sr_{2}RuO_{4}. More recently [35], f-wave superconductivity was observed in AV_{3}Sb_{5}, where , Rb, or Cs. Using our general description in the *naïve* model for any , we have been able to describe new types of possible superconductivity behaviour, including the f-wave case, which could help describing these new systems.

#### Appendix

#### A. Higher-Order Singularities at the Horizon

In Section 2.3, the background metric has been generalised to the form introduced by Equation (8) in an effort to fully generalise the description of holographic superconductors. However, the function in Equation (8) defines the singularity at the horizon to be of first order. One might then ask whether functions of the more general form could also give rise to the desired superconducting behaviour. In this appendix, it will be briefly argued that the case provides in fact the simplest scenario for the phenomenological description of superconductors.

To look at this issue, one may reproduce the steps in Section 3 that lead to the analytical calculation of the superconducting condensate. When examining the regularity conditions near the horizon, Equation (18) is modified into where is given by Equation (13) and is regular at .

Starting with , this implies that there are two cases: (i): then, in order to satisfy the EOM for at the horizon that is analogous to Equation (15) in the case , must be satisfied. As a consequence, one also finds using the Taylor expansion near the horizon. This means that the matter field disappears up to second order at the horizon and the condensate vanishes(ii): let be given by a power-law solution of the form near the horizon; then, it is . Therefore, scales with a larger power of , contradicting Eq. (A.2).

In both cases, does not contain a quadratic term which, as can be inferred from the results in Section 3, is one crucial ingredients for the emergence of a condensate. Consequently, there is no condensate in this case.

For , there are similarly two options: (i): if that is the case, the EOM for similar to Equation (15) enforces and, after the Taylor expansion Equation (21), also . Again, vanishes again up to second order at the horizon(ii): if had a power-law form , then . This implies by Equation (A.2). That is, for , solutions could exist, but no analytical solution could be found in this case

Therefore, no condensate emerges when , either, for the family of solutions under consideration. To summarise then, the condensate vanishes for , at least when considering the family of solutions obtained from power series by enforcing regularity conditions on the fields.

#### B. Dependence of the Analytical Results on the Matching Point

In Section 3, an analytical method was presented that allows to disentangle the dynamics of holographic superconductors. Said method consists in the construction of asymptotic solutions of the EOMs of a theory both near the event horizon and the AdS-asymptotic region, which are then glued together at an intermediate point. In this appendix, a technical aspect of this calculation is briefly clarified: the importance of the choice of the matching point to the coefficient of the condensate.

Keeping an explicit dependence on , the condensate reads where the dependent coefficient is given by and and are the coefficients defined by Equations (22) and (23), respectively. This expression coincides with Equation (33) when .

The only problematic point within the range of definition of (recall that ) is the pole which arises at whenever . Notably, as Equation (B.2) shows, the pole scales as . Therefore, the smaller the scaling dimension of the condensate is, the steeper and narrower is the growth of towards the pole at , making the overall dependence with milder. This behaviour, which can be spotted in Figure 5, justifies the choice of the value of that saturates the BF-bound for an optimal performance of the semianalytical calculation method that has been used in this paper.

#### C. The Electric Conductivity

Two different quantities typically characterise holographic superconductors. One is the condensate, whose emergence has been readily checked in this paper using the analytical matching method presented in Section 3. Additionally, the optical conductivity is worth being considered, since it generically exhibit a recognisable gap at some frequency in superconducting systems. With this as its outstanding feature, the conductivity can be used to classify systems as superconductors.

However, little attention has been brought to the conductivity throughout this paper. The reason for the lack of focus on that quantity is that, as it will be shown in this last appendix, the calculation method that has been made use in this paper breaks down when an analytical computation of the conductivity is attempted.

Nevertheless, we find it illustrative to take a look at the conductivity to show that the presented matching approach could provide a suitable starting point to derive semianalytical calculation methods for the conductivity. The modification of the analytical method to account for the conductivity is left for a future work.

In order to understand why this breaking happens, the results emerging from the analytical method will be put in comparison with the numerical outcome for the simple case of the s-wave superconductor, .

So as to obtain the conductivity, one starts by considering the fluctuations of the EM gauge field in a direction transversal to the AdS-radial dimension . Since the s-wave superconductor is rotationally invariant in its spatial dimensions, it is sufficient to consider fluctuations in the x-direction only:

Then, the EOM for this fluctuation in a generic BH background reads

As described in Refs. [4, 5], numerical solutions to this equation can be straightforwardly obtained once the EOMs for and have been solved. The resulting numerical conductivity then displays the expected behaviour, as it can be observed in Figure 6, where its real and imaginary parts are represented in red. In particular, its real part is negligible for smaller frequencies and becomes asymptotically unity as the frequency surpasses , shaping the aforementioned gap around that frequency.

Then, motivated by its successful description of the condensates (c.f. Section 3), one may attempt to apply the same analytical method to the calculation of the conductivity. For that purpose, the gauge field fluctuations near the horizon can be approximated by since near the horizon, with while the metric asymptotically becomes .

Afterwards, the coefficients , , and are fixed by expanding and solving the EOM Equation (C.2) up to .

On the other hand, in the asymptotic region, the fluctuations can be written as

As done in the case of the condensate, the expansion near the horizon, Equation (C.3), is matched to the asymptotic expansion near the AdS boundary, Equation (C.4), by imposing the Dirichlet and von Neumann boundary condition at the matching point , i.e.,

These two equations determine the coefficients and . According to, e.g., Ref. [5], the conductivity can be defined in terms of those coefficients as

However, when the analytical coefficients are plugged in this definition, the outcome is very far from resembling the reliable numerical solution described before. This is also to be seen in Figure 6, where the real and imaginary part of the analytically calculated conductivity are put in comparison with the numerical outcome. In fact, one can see that while some gap-like structure develops for , the analytical conductivity does not retain its proper normalisation and grows indefinitely as . This is a nonphysical behaviour in what respects the description of superconductors that indicates the breaking of this method of calculation.

To understand the reason for this failure, one may look at the expansion in Equation (C.3), whose prefactor forecasts an oscillating behaviour of the solution near the horizon. This oscillating solution is then matched to an a priori nonoscillating expansion near the asymptotic boundary, Equation (C.4). Therefore, a high dependence on the matching point is expected that spoils any hope of finding a proper expression for the conductivity using this method. Nevertheless, more refined matching procedures, beyond the one employed in this paper, may lead to a successful analytical treatment of the conductivity properties of SCs. We leave this study for a future work.

#### Data Availability

No data were used to support this study. Whenever reference is made to previous work, the corresponding paper is cited.

#### Disclosure

A preprint has previously been published on arXiv [36].

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

VS would like to thank Daniel Errandonea for pointing out relevant works in f-wave superconductivity and Fëanor Reuben Mitchell Ares for discussions on the analytical matching procedure. This work was supported by the grant PROMETEO 2019/08 and the Spanish FPA 2017-85985-P. FE is supported by the Generalitat Valenciana with the grant GRISOLIAP/2020/145. VS acknowledges support from the UK Science and Technology Facilities Council ST/L000504/1. This project has received funding support from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement 860881-HIDDeN.