Supersymmetric Lifshitzlike backgrounds from SYM with heavy quark density
Abstract:
We examine a class of gravity backgrounds obtained by considering the backreaction of a spatially uniform density of mutually BPS Wilson lines or heavy quarks in SUSY YangMills theory. The configurations preserve eight supercharges and an subgroup of the Rsymmetry. They are obtained by considering the BPS geometries associated to smeared string/D3brane (F1D3) intersections. We argue that for the (partially) localized intersection, the geometry exhibits a flow from in the UV to a novel IR scaling solution displaying anisotropic Lifshitzlike scaling with dynamical critical exponent , hyperscaling violation and a logarithmic running dilaton. We also obtain a twoparameter family of smeared BPS solutions on the Coulomb branch of SYM exhibiting Lifshitz scaling and hyperscaling violation. For a certain parametric range these yield IR geometries which are conformal to , and which have been argued to be relevant for fermionic physics.
1 Introduction
The finite density physics of strongly interacting quantum systems presents significant challenges for existing theoretical frameworks. The questions of physical interest range from the behaviour of quark/baryonic matter at high densities [1, 2, 3] to the physics of quantum critical points in condensed matter systems [4]. The AdS/CFT correspondence [5] provides an avenue for tackling such phenomena within toy models that can be constructed and analyzed using the framework of gauge/string duality [6, 7]. One of the directions that has attracted considerable attention in recent years is the possible emergence of Lifshitz [8, 9] and Lifshitzlike (with hyperscaling violation) [10, 11, 12] scaling symmetries as infrared (IR) descriptions of strongly coupled field theories with holographic duals. A particularly fascinating feature of such scaling is the possibility of observing logarithmic violation of the area law for entanglement entropies indicating the presence of hidden Fermi surfaces [10, 11], and other possible signatures of holographic fermionic physics [13].
The aim of this paper is to exhibit some supersymmetric backgrounds in string theory resulting from the backreaction of distributions of fundamental strings (F1strings) in the presence of a large number of D3branes, giving rise to (IR) geometries exhibiting a range of Lifshitzlike scalings with and without hyperscaling violation^{1}^{1}1String/Mtheory realizations of backgrounds with hyperscaling violation (e.g. from lightlike compactification of AdS plane waves) and their properties have been studied in [14, 15].. Our primary motivation for considering such configurations arises from the interpretation of macroscopic fundamental strings as heavy quarks or Wilson lines [16] in the holographic dictionary between strongly coupled SUSY YangMills theory in the ’t Hooft large limit, and IIB string theory in . Spatially uniform distributions of macroscopic strings correspond to a finite density of heavy quarks in the gauge theory^{2}^{2}2It can also be argued that such configurations in the bulk should control the IR geometry produced by the backreaction of massive and massless flavour D7branes with chemical potential for baryon number. We thank D. Mateos and J. Tarrio for stimulating discussions on this point.. It has been found in [17] that a specific nonsupersymmetric smearing of such string sources, preserving global symmetry, produces a flow from to a Lifshitz geometry with critical exponent and a logarithmically running dilaton (see also [18]).
Spatially homogeneous dimensional spacetime metrics exhibiting Lifshitzlike scaling with dynamical critical exponent and hyperscaling violation have the form [10, 11, 12]
(1) 
so that under the rescalings: , and , the line element transforms covariantly as .
We will encounter such scaling backgrounds with , preserving 8 supersymmetries, as IR descriptions of SYM at strong coupling with a smeared density of supersymmetric Wilson lines or heavy quarks. In the large limit, we take the number density of such heavy quarks to scale as , so that it becomes necessary to include their backreaction on the field theory vacuum. The holographic duals to such states are described by a class of BPS backgrounds in IIB supergravity whose general properties will be presented elsewhere [19]. The specific class of backgrounds that we consider in this article are those which have the form of certain well known intersecting brane configurations [20, 21, 22]. The configurations within this class can be thought of as pairs of quarks and antiquarks which are mutually supersymmetric (by having opposite internal orientations) and are distributed uniformly in the spatial directions of the gauge theory. A particular solution within this class has appeared in [15, 23] where its Lifshitzlike scaling properties have been pointed out.
We find that one category of the intersecting F1D3 configurations naturally leads to an exact scaling solution in the infrared, exhibiting Lifshitzscaling with which is mildly broken by a logarithmically running dilaton. The background preserves supersymmetry and an internal symmetry. This result is noteworthy since the Lifshitz scaling behaviour matches that of [17] which dealt with an preserving nonsupersymmetric setup. This is indicative that the scaling is probably a universal feature of SYM with finite quark density (at strong coupling and large). It points toward the possibility that when the quark flavours are made dynamical (by introducing D7branes for instance [24]) and their backreaction taken into account at finite baryon density, then the low temperature IR physics may well be controlled by a similar scaling solution (see [25] for closely related discussions). A concrete framework where this possibility can be further investigated is the unquenched smeared flavour brane setup of [26, 27, 28, 29].
Within the setting of the supersymmetric configurations in this paper, constructing the full flow from to the IR scaling solution remains a challenging task due to the reduced isometry of the dual gravity backgrounds. However, we analytically examine the departure of the flow induced by the supersymmetric smeared strings in the UV, and find close similarities to the symmetric case of [17] wherein it was possible to obtain the entire flow numerically.
The supersymmetric intersecting brane configurations present us with a further interesting route for obtaining nontrivial scaling solutions. This is because of the existence of supersymmetric moduli spaces at zero temperature. In particular, it is possible to move out to a generic point on the Coulomb branch moduli space (where the gauge group is generically Higgsed to a product of Abelian factors) of SYM and examine the effect of a finite quark density. In the dual gravity description this corresponds to a general distribution of D3branes at large. We consider such continuous distributions of the D3branes and F1strings that are compatible with the supersymmetries preserved by the intersecting brane configurations. Specifically, we introduce distributions for these sources with power law scalings in the IR and solve the supergravity equations for such power law density functions. We thus find a two parameter family of Lifshitzlike scaling solutions with hyperscaling violation coefficients. Imposing the weak energy condition on each of the source distributions restricts the allowed values of and .
Perhaps the most interesting result within the class of Coulomb branch configurations is the appearance of a solution with and fixed (satisfying ). These give rise to geometries that are conformal to with vanishing entropy at zero temperature. Such geometries have been argued to be relevant for holographic descriptions of fermionic physics [13].
All the solutions we find bear close resemblance to scaling solutions in EinsteinMaxwelldilaton theories [9, 30, 31, 32]. In particular the dilaton softly breaks the scaling invariance via logarithmic running. Furthermore since the dilaton runs to zero in the IR, it renders corrections important in the deep IR.
2 Smeared strings and Sym
SYM theory possesses a global Rsymmetry which corresponds to the isometry of the fivesphere of the dual gravity background at strong coupling. A supersymmetric Wilson line in SYM at strong coupling is computed by an infinite probe F1string stretching from the boundary of along the radial AdS coordinate, whilst remaining localized at a point on the internal [16]. Such a configuration preserves and more generally, is associated to the MaldacenaWilson line in some representation ,
(2) 
where is a unit vector in and , the six adjoint scalars in the multiplet. We focus our attention on straight (timelike) BPSWilson lines. There also exists a gravity dual realization of such supersymmetric or BPSWilson lines for generic representations in terms of D3 and D5brane embeddings in [33, 34, 35, 36]. In particular, for representations whose size scales as , the gravity duals of the Wilson lines necessarily require inclusion of backreaction of the corresponding D3/D5brane configurations on , leading to the BPS ‘bubbling’ geometries of [37].
Our goal in this paper is to model a particular state with finite heavy quark density by taking a very specific BPS Wilson line configuration and smearing it uniformly along all spatial directions of the theory. The unsmeared configuration consists of a mutually BPS quarkantiquark pair i.e. antiparallel lines placed at antipodal points of the internal [16, 38]. This preserves an global symmetry and, prior to taking the near horizon limit of the D3branes, can be viewed as a supersymmetric F1D3 intersection. A specific delocalized version of this intersection has been discussed in [15, 23] and shown to lead to Lifshitzlike scaling. The backreacted geometries resulting from the spatially smeared configurations are BPS, preserving 8 real supercharges.
We will consider two distinct classes of BPS configurations in this paper:

In the first category lies the socalled partially localized F1D3 intersection [22], where F1strings are smeared along the relative transverse directions (spatial coordinates) on the D3branes, and the D3branes are not smeared. We will argue that the corresponding gravity background represents a flow from to a Lifshitzlike geometry.

The second class involves delocalized D3brane and F1string distributions, allowing for a general smearing density for both sets of sources. The solution of [15, 23] is a special case within this category of solutions. We find a large class of zero temperature scaling solutions with a range of dynamical critical exponents and hyperscaling violation coefficients.
2.1 Metric ansatz in 10D
Given the symmetries and supersymmetries of such configurations, we show in an accompanying paper [19], that the type IIB fermionic variations vanish if one takes the following ansatz for the Einstein frame metric, the RamondRamond fluxes and NeveuSchwarz field strength :
(3) 
Note that we follow the conventions of [37] and our normalization of differs from the usual one by a factor of 4. Here is the metric of the unit radius foursphere with volume form , the warp factor and dilaton are functions of the and coordinates only. The gauge theory lives in the 4D spacetime spanned by the coordinates . In addition, the tendimensional complex spinor parametrising supersymmetry variations must satisfy the projection conditions
(4) 
leading to BPS solutions. Here is the ten dimensional complex conjugation matrix. The projections are those associated to fundamental strings and D3branes, respectively.
Within this ansatz, the geometry can be viewed locally as a fibration of over the  plane. Therefore, these are time independent backgrounds with an isometry and symmetry under translations and rotations in .
2.2 (Partially) Localized F1D3 intersection
In this paper we will only study solutions with . This precludes a charge density for D5branes that could be interpreted as baryon density as in [17]. It is also physically clear why this is the case: we are considering a uniform density of quark and antiquark pairs (with opposite orientations). With this choice, the backgrounds represent F1D3 BPS intersections of the kind analyzed in [20, 21, 22]. We can make contact with the standard form of the metric for such backgrounds with the identifications:
(5) 
where and are the harmonic functions associated to the F1strings and D3branes respectively. It is well known that in these types of solutions (e.g. [22]) there is a smearing along the relative transverse directions of the strings/branes i.e. those which are parallel to the world volume of some of the branes (strings), but perpendicular to others. They are constructed by the socalled Ôharmonic ruleÕ: to each type of brane (string) there is a corresponding function which is harmonic in the space transverse to their worldvolumes. In order to source two sets of branes in the same solution one simply multiplies the two harmonic functions together in the appropriate way.
vacuum:
The vacuum solution, namely , is recovered when
(6) 
This is indicated in Fig.1. The fivesphere is obtained by fibering the along a semicircle enclosing the origin.
Backreaction of string sources:
The presence of macroscopic string sources will lead generically to a nonvanishing NSNS threeform flux . The form of in our ansatz (3) suggests that the strings must lie parallel to the axis (see Fig.1). Now, both the functions and are nontrivial, and for sources localized in the  plane they are determined by the Bianchi identity for and equation of motion for , namely,
(7) 
These are Laplacelike equations for (partially) localized sources smeared along certain directions transverse to them (the transverse and the spatial ). Since , we also know that is only a function of the coordinate,
(8) 
This is what we would expect for a localized, infinite string source oriented along the axis. The additive integration constant has been set to unity by the requirement that the dilaton vanishes for large . Here is proportional to the number density of Fstrings , which also measures the density of heavy quark sources in the gauge theory:
(9) 
being the ’t Hooft coupling of the theory. The number density of heavy quarks must scale as in order to keep fixed in the large limit, and to consistently include their backreaction on the background.
The normalizations are fixed by demanding that the equation of motion for yield a number density of macroscopic Fstrings localized at ,
(10) 
Setting the AdSradius to unity, the ’t Hooft coupling of the gauge theory is related to the tension of a fundamental string as . In addition, Newton’s constant in ten dimensions is given by , with .
With this solution for , the full background is determined by , which solves eq.(7) subject to the requirement that the geometry is asymptotically, locally . It is not possible to solve this equation analytically, but we can analyze its UV and IR asymptotics.
2.3 Lifshitz IR
In the IR limit, which corresponds to small , we have that . The general solution for the harmonic function in this limit is (e.g. [22]):
(11) 
We may view this as a ‘near horizon’ limit for the Fstrings. In fact can take the more general form, similar to that for multicentred D3brane sources, , but this generalization is not of particular interest to us at this point. It is easy to see that the solution for follows from a scaling argument^{3}^{3}3Allowing for a source density on the right hand side of (7), we expect that under , , for a scaling solution we should have . Comparing the two sides we find that and for some function . The differential equation for , yields the general solution .. The second important remark is that such solutions cannot correspond to localized F1D3 intersections, since increases as a positive power of for large . The D3branes are in fact delocalized along the axis in this limit (). The reasons behind such ‘spontaneous delocalization’ were explored in [39], including the case which is precisely the Sdual of the setup we are studying, namely D1branes smeared along the spatial directions of D3branes. In this case the spontaneous delocalization was attributed to the analogue of the ColemanMerminWagner theorem in the 0+1 dimensional quantum mechanics of the intersection.
One can confirm that the D3branes are indeed not localized by noting that (in this limit and with for simplicity) depends only on ,
(12)  
Therefore, the D3brane charge can be computed by where is a path joining two points on the curves, say, and . For any two points on the axis, which sits at , the integral along a path joining them vanishes trivially, and therefore the D3brane charge is not localized at the origin, nor at any other finite point on the axis. It also follows, by allowing arbitary deformations of the curve, that there are no localized D3brane sources for any finite values of and . In fact the sources must be viewed as being placed at infinity. This is also related to the fact that the in the geometry does not shrink anywhere in the  plane in this IR limit. This should be contrasted with wherein the D3brane charge is obtained by choosing a path with endpoints on the axis and, importantly, enclosing the origin, as in Fig.1.
This picture therefore suggests that the Fstrings have pulled apart the D3branes so that for scales set by , they appear to be ‘localized’ at infinity. Within the full flow from the UV AdS asymptotics we expect the D3branes to have a nontrivial distribution (and to therefore be delocalized) along the axis. Since the bulk geometry depends on two holographic directions, and (or and ), the D3brane distribution will emerge from numerical analysis of the entire flow, a study which we postpone for the future.
Plugging in the expressions for and in this IR limit, defining a new radial coordinate and after appropriate coordinate rescalings, we find
This is an exact, supersymmetric solution to the type IIB equations. The constants and can only be determined upon embedding this IR solution into the complete flow, with AdS UV asymptotics, triggered by the string charge density. The IR background metric exhibits a scaling symmetry under
(14) 
which is an anisotropic Lifshitzlike scaling, with dynamical critical exponent , in the 4D gauge theory. Note that the scaling symmetry is not exact, as it is broken by the logarithmic running of the dilaton with energy scale.
Let us point out certain noteworthy aspects of the scaling solution (2.3). The IR Lifshitz scaling is realized in an unusual fashion due to the nontrivial dependence of the metric on two holographic directions, and , and this will further affect the scaling properties of observables such as the entanglement entropy and thermodynamic quantities. It is remarkable that the Lifshitz scaling, with exactly the same running of the dilaton, was also observed in [17, 18] within the nonsupersymmetric symmetric ansatz for smeared Wilson lines in SYM. This suggests that the emergent Lifshitz symmetry with may be a universal feature of a dense state of heavy quarks introduced in SYM, independent of the details of their internal orientation. A further difference between the symmetric setup of [17, 18] and the supersymmetric configuration above, is that whilst the former has corresponding to D5brane/baryon charge density, the latter has . However this detail does not appear to affect the dynamical critical exponent in the IR. A supersymmetric configuration with is explored in [19].
We can recover a more conventional form of the Lifshitzlike metric by taking a further limit of the IR solution, namely , assuming that generally. After some coordinate rescalings, this yields
(15) 
where is a constant. The foursphere has a constant size in this limit, and the background (15) is not, by itself, an exact solution to the supergravity equations of motion; it is a limiting form of the IR (small ) metric.
An interesting aspect of the IR physics that is manifest in this limit, is that it is effectively a sixdimensional geometry (after reduction on ). The emergence of an extra noncompact coordinate in the IR geometry suggests that the IR field theory of the smeared heavy quark impurities in SYM has an emergent dimension. Such an interpretation would be consistent with the presence of a delocalized distribution of Dbranes that deconstruct an extra dimension [40]. Given the complicated nature of the putative RG flow solution, and the dual field theory being at strong coupling, it is not easy to pinpoint a direct test of this proposal. Our study of smeared strings on the Coulomb branch in the next section, lends support to this picture. It is known that Coulomb branch configurations in SYM deconstruct a higher dimensional field theory [41]. Within such an interpretation the deconstructed dimension should be viewed as being compact, its size being determined by the inverse spacing between two adjacent Dbranes. The anisotropy between the deconstructed coordinate and the spatial coordinates is automatic due to the distribution of heavy quark sources.
Entanglement entropy:
An important probe of the IR physics of the gauge theory is the entanglement entropy, which can be computed using the RyuTakayanagi prescription [42]. For the metric in eq.(2.3), this is not a straightforward calculation since the components depend both on and coordinates and the solution for the extremal surface (with a prescribed boundary) could depend nontrivially on both these dimensions.
We can, however, make a simplifying approximation by assuming that the physics in the deep infrared should be governed by sufficiently small and , so that . Then we can use the approximate form (15) to calculate the entanglement entropy of a ‘strip’ in , specified by
(16) 
with at a UVslice of the geometry, . The effective area functional for the 3D surface with the strip as its boundary, can be defined as [42]
(17) 
Extremizing the functional and extracting its finite part in the usual way (e.g. [10, 12]) we find
(18) 
where we have omitted various constants of proportionality and traded Newton’s constant in ten dimensions for a factor of according to the holographic dictionary. The scaling of the entanglement entropy with is characteristic of hyperscaling violation with negative [10, 12]. If we view the (IR) gauge theory as having 3 spatial dimensions (e.g. with the coordinate compactified), then .
The hyperscaling violation could also be directly inferred by reducing the ten dimensional metric (15) to five dimensions, by assuming the coordinate to be compact. The reduction to five dimensions yields
(19) 
which exhibits Lifshitz scaling with and hyperscaling violation with . For the configurations on the Coulomb branch studied below, the emergence of the extra spatial coordinate can be interpreted via deconstruction and the size of this dimension is controlled by the inverse spacing between D3branes [41, 40].
We stress that the calculation of the entanglement entropy outlined above is only an approximation. It would be interesting to have a more precise computation in the exact scaling background (2.3), and to check whether these qualitative expectations are reproduced.
2.4 UV AdS asymptotics
We will now attempt to understand how the solution to eq.(7) will modify the UV AdS asymptotics. In order to find a solution which asymptotes to , and also includes string sources, we need to use the general form (8) for . Whilst the problem is not analytically tractable, we can make progress by noting that eq.(7) exhibits translational invariance in . Fourier transforming with respect to this variable [43] yields
(20) 
where primes denote derivatives with respect to . The equation has irregular singular points of order at and order at . Around either of these points, the equation can be solved as a power series using a socalled Thomé expansion [44]. Around , this leads to solutions of the type,
(21) 
where the coefficients can be determined via a recursion relation and can be either positive or negative. A slightly different expansion which has overlap with the large limit is a formal series expansion in powers of – the string or ‘heavy quark’ density:
(22) 
Substituting into (7), and solving the resulting equations order by order in , we find,
(23)  
for the first few terms in the expansion ( can be either positive or negative). Ei is the exponential integral function. The correct superposition of the Fouriertransformed solutions must reproduce in the limit . This picks out the required combination,
(24) 
In terms of the standard radial coordinate of AdS and the polar angle on the (, ) to linear order in , we have,
(25) 
The dilaton, on the other hand, is exactly determined by
(26) 
and vanishes at i.e. the north and south poles of the in the UV geometry, as expected for F1string sources localized at . It is instructive to look at the corrections to the metric components in the above expansion,
(27) 
with similar corrections appearing for . The corrections to AdS asymptotics are characteristic of the backreaction due to string sources, as also noted in [45, 46, 17]. A similar analysis for the internal directions suggests that upon inclusion of the backreaction from the strings, the does not shrink at , and that there are curvature singularities at these points. We should stress, however, that we cannot trust the expansion in in the vicinity of the string sources at . Away from these points (fixed generic and large ) the geometry is asymptotically, locally .
To summarize, we have argued in this section that the introduction of an density of heavy quarks (specifically, mutually BPS quarkantiquark pairs), preserving some supersymmetry in SYM at large and strong coupling, induces a flow from to a scaling Lifshitzlike IR background with . The anisotropic scaling solution in the IR is somewhat novel due to the presence of two noncompact holographic directions, which scale differently to ensure that the resulting metric is scale invariant.
3 Coulomb branch solutions and hyperscaling violation
A notable feature of the supersymmetric intersections discussed above is that they inherit the Coulomb branch moduli space of the theory (where is Higgsed to generically). In particular a probe D3brane (spanning ) placed at any point in the  plane, in the general background (3), experiences no force due to an exact cancellation between the DiracBornInfeld (DBI) and WessZumino terms in the probe action. A similar cancellation occurs (between the NambuGoto action and the coupling to ) for a probe Fstring oriented parallel to the axis at any value of the coordinate. This points towards the existence of more general smeared solutions where both D3branes and the Fstrings are smeared with some distribution functions on the  plane. We show that for certain choices of such distribution functions the backreacted (IR) geometries can exhibit a range of Lifshitzlike scalings with hyperscaling violation.
General solutions with D3branes on the Coulomb branch (within our symmetric ansatz) can be obtained by considering a general D3brane density, so that
(28) 
A general density function preserves all supersymmetries whilst also leading to a violation of the Bianchi identity for due to the extended source distribution.
Since we want to interpret our solutions below as IR limits of Coulomb branch distributions it is useful to illustrate this with a simple example which allows to make contact with the F1D3 intersection of [15, 23]. Let us first consider a uniform distribution of D3branes, in an interval along the axis with
(29) 
and no macroscopic string sources. For large the geometry asymptotes to , whilst in the IR, for small and , we obtain
(30) 
This is a scaling regime where all metric components in (3) only depend on powers of . Upon performing a reduction of (3) on and the coordinate, the resulting 5D geometry precisely matches the reduction on a torus of the solution in 11D SUGRA [41]. This symmetric Coulomb branch configuration can therefore be viewed as a flow at strong coupling and large, to an IR theory which appears to be (a subsector of) the 6D superconformal theory with supersymmetry realized on a stack of M5branes (compactified on a twotorus). The Coulomb branch configuration therefore deconstructs a higher dimensional field theory. In the deconstruction picture, the sizes of the extra dimensions are controlled by the inverse spacings between individual D3branes on the Coulomb branch. The spacings determine the masses of the lightest Wbosons and dyonic states on the Coulomb branch, which in turn are related to the KaluzaKlein harmonics of the deconstructed compact dimensions. We would now like to understand how the IR dynamics on the Coulomb branch is modified by smeared, macroscopic fundamental strings viewed as heavy quarks in the gauge theory.
3.1 Smeared F1D3 intersections
In the presence of macroscopic string sources, both and will be nontrivial and the Bianchi identity for modified by D3brane sources becomes
(31) 
We will focus attention on distributions that, at least in some limit as in the above example, only depend on , so that . Hence, for simplicity, we also assume that is a function of alone. Similarly the equation of motion for with a general smeared distribution of Fstrings leads to the equation
(32) 
As pointed out above, probe Fstrings experience no force when they are aligned along the axis and placed at any value of . Therefore a general dependent distribution should preserve supersymmetry. In the Appendix we show that such smearing of sources satisfies a calibration condition which ensures that the string/brane embeddings we study respect the supersymmetries. The equations of motion follow from the type IIB supergravity action coupled to the smeared D3brane and string sources,
(33) 
The actions for the sources are the smeared versions of the NambuGoto and DBI actions with a particular choice of smearing forms and , respectively
(34)  
(35) 
where we have defined . The sources alter the equations of motion for the dilaton, , and the metric. In particular they contribute to the stress tensor. Imposing the weak energy condition (WEC) on the source stress tensor leads to positivity of the source density functions
(36) 
where is some timelike vector. For more details on the equations of motion and their consistency with supersymmetry, we refer the reader to the Appendix.
To obtain scaling solutions, we simply choose power laws for the smearing densities
(37) 
Then the equations of motion for and are also solved by power laws
(38) 
Positivity of the source densities requires . Further imposing the null energy condition (NEC) on the complete stress tensor leads to a weaker condition on the parameters which is consistent with the WEC . The case gives us the homogeneous solution of [15, 23]. This is the situation wherein the D3branes and strings are all at , but the branes are uniformly smeared along the axis.
Inserting the general scaling solution into the metric (3) and reducing carefully over the and the coordinate (which we treat as a deconstructed compact dimension) to fivedimensional Einstein frame we obtain a family of Lifshitz geometries with hyperscaling violation:
(39) 
Here the radius is and we have changed the radial coordinate according to and rescaled the rest. The coefficients are related to the exponents of the sources as
(40) 
We stress that the emergence of hyperscaling violation along with Lifshitz scaling in the 4D spacetime is already evident in the ten dimensional metric (3), upon direct substitution of the scaling solutions (38). Reduction to 5D changes the hyperscaling violation parameter, but leaves the Lifshitz exponent unchanged.
The case needs to be treated separately, as we do below. Given that the weak energy conditions for the sources, considered separately, restrict the values of the parameters to lie in the range , the allowed values for and in this class of solutions are shown in Fig. 2. Note that if we also require for stability as argued in [12], only the solutions shown in the lower right corner of Fig. 2 with survive. These correspond to D3brane distributions with . Finally, there are no solutions with , which would lead to a logarithmic violation of the area law for entanglement entropy.
: Conformally :
As stated before, the case has to be treated separately, since the change of variables used for the radial coordinate is not well defined for that value of the parameter. Importantly as , both and diverge,
(41) 
In the allowed range for , we have . The limit is also interesting because it corresponds to a uniform distribution of D3branes on the  plane,
(42) 
so that the harmonic function is,
(43) 
and as before. Hence we obtain a oneparameter () family of solutions. Again, substituting this into the tendimensional metric and reducing to the fivedimensional Einstein frame with the change of variable we arrive at a metric conformal to AdS,
(44) 
where the radius is now . Geometries with an IR factor conformal to have been used to describe locally quantum critical theories and argued to have certain properties desirable for a holographic description of Fermi surfaces [13]. The zero temperature entropy density associated to the conformally background is vanishing since and the entanglement entropy has no extensive finite contribution. As pointed out in [13], the entanglement entropy of a strip in this case has a rather puzzling behaviour in that the extremal surface only exists for a very specific value of the strip width.
In all the geometries, including the conformally backgrounds, curvature singularities appear in the deep IR (see e.g.[47]). In particular, since the dilaton is given by , and runs to zero in the IR, the string frame metric leads to divergent curvatures and the solutions will be expected to receive large corrections. One can Sdualize the backgrounds, and then the problem becomes one of strong coupling in the IR. It is a compelling question to ask what low energy physics arises and resolves such curvature singularities [48, 49]. This may be an interesting avenue to explore since we have some knowledge of the microscopic field theory the geometries describe.
4 Discussion
The supersymmetric scaling geometries we have discussed in this paper were obtained from previously known intersecting brane configurations along with a new ingredient, namely, generic source distributions compatible with the supersymmetries. Our motivation was to understand whether any universal features emerge within a tractable holographic framework, when a state of finite (quark) density is introduced into a known field theory (with a large string/gravity dual). The physical interpretation of the setup explored in this paper is that it corresponds to a uniform density of mutually BPS quarks and antiquarks (with opposite orientations) in SYM. In particular, these configurations all have vanishing three form flux, . The presence of nonzero can be interpreted as a finite baryon density, as was the case in [17] (see also [50] for related discussions). In a forthcoming publication [19] we will derive the general form for the BPS configurations involving smeared strings with nonvanishing .
One of the larger aims of this study is to understand whether the simplistic picture of backreacting smeared heavy quarks can be embedded into holographic setups with dynamical flavours, most notably the smeared D3D7 system explored, for instance, in [27, 28, 29].
An important observation in this paper is the emergence of a IR scaling regime from the partially localized intersection described by homogeneous (sourceless) equations. Despite having different global symmetries and no supersymmetry, the same scaling was observed in [17, 18]. The value of can now be understood within the intersecting brane framework. Following the results of [22], it is easy to verify that, more generally, partially localized F1Dp intersections give rise to anisotropic scaling with , when .
Finding finite temperature, black brane generalizations of the solutions discussed in this paper appears difficult. For the partially localized intersection (with ), the solutions are functions of two coordinates, and , and this makes the problem of ‘blackening’ the solutions challenging. We leave this question for future investigation. When the sources are smeared, as for the Coulomb branch configurations, the zero temperature delocalization is possible only because the theory has a moduli space of vacua. Finite temperature lifts such moduli spaces so that the free energy is minimized at the origin of moduli space^{4}^{4}4Away from the origin, various degrees of freedom are rendered massive, and their contribution to the entropy thus decreases.. A possible approach towards the finite temperature physics of smeared Coulomb branch distributions is to introduce chemical potential for charges under the global symmetry group, i.e. to consider rotating brane configurations along the lines of [51, 52].
Within the context of the Coulomb branch intersctions, it is interesting to note that by engineering appropriate string/D3brane sources we should be able to obtain zero temperature flows between different scaling geometries.
Acknowledgments:
We would like to thank Paolo Benincasa, Jerome Gaillard, Tim Hollowood, David Mateos, Carlos Nuez, Andy O’Bannon, Alfonso Ramallo and Javier Tarrio for illuminating discussions and helpful comments on various aspects of this work. We acknowledge partial support from the STFC grant ST/G000506/1. The work of A. F. was supported by STFC grant ST/J00040X/1.
Appendix : Calibration conditions and equations of motion
We work in Einstein frame in the conventions of [37]. We restrict our discussion to the situation with , and dependent density distributions and harmonic functions. Smearing the fundamental strings along all the transverse directions, the Nambu–Goto action plus coupling to the NS form are schematically,
(45) 
with a particular choice of smearing form
(46) 
The function of the radial coordinate describes the string charge distribution. Associated to a string worldsheet embedding, one introduces a calibration form, given essentially by the induced metric on the string,
(47) 
In the presence of smeared sources the calibration condition, which ensures that the embedding of the branes/strings respects supersymmetry, has to be modified (see for instance [29]). In the case of fundamental strings the modified condition is
(48) 
Using that , one can easily check that this condition is verified in our ansatz, establishing that general backreacting string distributions are consistent with the supersymmetries of the setup. Similarly, we introduce a set of D3 branes extending along the 4D Minkowski directions and smeared on the transverse coordinates , with action
(49) 
In this case the smearing form is
(50) 
while the calibration form reads
(51) 
where again the function parametrizes the brane charge along the radial direction. The pertinent calibration condition
(52) 
is also straightforwardly satisfied, again in line with the expected supersymmetry of the backgrounds.
The presence of these smeared sources alters several equations of motion, that now read
(53) 
The equations for and are not modified and turn out to be automatically satisfied within our ansatz. Notice that couples electrically to the strings through the Nambu–Goto action, so the Bianchi remains intact. As is customary the Bianchi for coincides with its equation of motion. The stress tensors for the smeared sources are
(54) 
where for arbitrary forms and we have defined the operations
(55) 
Using the known ansatz for the forms, all the equations above boil down to the equations determining the harmonic functions sourced by the charge distributions
(56) 
References
 [1] K. Fukushima and T. Hatsuda, “The phase diagram of dense QCD,” Rept. Prog. Phys. 74, 014001 (2011) [arXiv:1005.4814 [hepph]].
 [2] K. Rajagopal and F. Wilczek, “The Condensed matter physics of QCD,” In *Shifman, M. (ed.): At the frontier of particle physics, vol. 3* 20612151 [hepph/0011333].
 [3] J. CasalderreySolana, H. Liu, D. Mateos, K. Rajagopal and U. A. Wiedemann, “Gauge/String Duality, Hot QCD and Heavy Ion Collisions,” arXiv:1101.0618 [hepth].
 [4] S. Sachdev, “The Quantum phases of matter,” arXiv:1203.4565 [hepth].
 [5] J. M. Maldacena, “The Large N limit of superconformal field theories and supergravity,” Adv. Theor. Math. Phys. 2, 231 (1998) [hepth/9711200].
 [6] S. A. Hartnoll, “Lectures on holographic methods for condensed matter physics,” Class. Quant. Grav. 26, 224002 (2009) [arXiv:0903.3246 [hepth]].
 [7] C. P. Herzog, “Lectures on Holographic Superfluidity and Superconductivity,” J. Phys. A 42, 343001 (2009) [arXiv:0904.1975 [hepth]].
 [8] S. Kachru, X. Liu and M. Mulligan, “Gravity Duals of Lifshitzlike Fixed Points,” Phys. Rev. D 78, 106005 (2008) [arXiv:0808.1725 [hepth]].
 [9] M. Taylor, “Nonrelativistic holography,” arXiv:0812.0530 [hepth].
 [10] N. Ogawa, T. Takayanagi and T. Ugajin, “Holographic Fermi Surfaces and Entanglement Entropy,” JHEP 1201, 125 (2012) [arXiv:1111.1023 [hepth]].
 [11] L. Huijse, S. Sachdev and B. Swingle, “Hidden Fermi surfaces in compressible states of gaugegravity duality,” Phys. Rev. B 85, 035121 (2012) [arXiv:1112.0573 [condmat.strel]].
 [12] X. Dong, S. Harrison, S. Kachru, G. Torroba and H. Wang, “Aspects of holography for theories with hyperscaling violation,” JHEP 1206, 041 (2012) [arXiv:1201.1905 [hepth]].
 [13] S. A. Hartnoll and E. Shaghoulian, “Spectral weight in holographic scaling geometries,” JHEP1207, 078 (2012) [arXiv:1203.4236 [hepth]].
 [14] K. Narayan, “On Lifshitz scaling and hyperscaling violation in string theory,” Phys. Rev. D 85, 106006 (2012) [arXiv:1202.5935 [hepth]]. K. Narayan, T. Takayanagi and S. P. Trivedi, “AdS plane waves and entanglement entropy,” JHEP 1304, 051 (2013) [arXiv:1212.4328 [hepth]]. K. Narayan, “Nonconformal brane plane waves and entanglement entropy,” Phys. Lett. B 726, 370 (2013) [arXiv:1304.6697 [hepth]].
 [15] H. Singh, “Lifshitz/Schródinger Dpbranes and dynamical exponents,” JHEP 1207, 082 (2012) [arXiv:1202.6533 [hepth]]; H. Singh, “Lifshitz to AdS flow with interpolating brane solutions,” JHEP 1308, 097 (2013) [arXiv:1305.3784 [hepth]].
 [16] J. M. Maldacena, “Wilson loops in large N field theories,” Phys. Rev. Lett. 80, 4859 (1998) [hepth/9803002]; S. J. Rey and J. T. Yee, “Macroscopic strings as heavy quarks in large N gauge theory and antide Sitter supergravity,” Eur. Phys. J. C 22, 379 (2001) [hepth/9803001].
 [17] S. P. Kumar, “Heavy quark density in N=4 SYM: from hedgehog to Lifshitz spacetimes,” JHEP 1208, 155 (2012) [arXiv:1206.5140 [hepth]].
 [18] T. Azeyanagi, W. Li and T. Takayanagi, “On String Theory Duals of Lifshitzlike Fixed Points,” JHEP 0906, 084 (2009) [arXiv:0905.0688 [hepth]].
 [19] A. F. Faedo, B. Fraser and S. P. Kumar (to appear)
 [20] A. A. Tseytlin, “Harmonic superpositions of Mbranes,” Nucl. Phys. B 475, 149 (1996) [hepth/9604035].
 [21] J. P. Gauntlett, “Intersecting branes,” In *Seoul/Sokcho 1997, Dualities in gauge and string theories* 146193 [hepth/9705011].
 [22] D. Youm, “Partially localized intersecting BPS branes,” hepth/9902208.
 [23] P. Dey and S. Roy, “Lifshitzlike spacetime from intersecting branes in string/M theory,” JHEP 1206, 129 (2012) [arXiv:1203.5381 [hepth]]; P. Dey and S. Roy, “Holographic entanglement entropy of the near horizon 1/4 BPS FD bound states,” Phys. Rev. D 87, no. 6, 066001 (2013) [arXiv:1208.1820 [hepth]].
 [24] A. Karch and E. Katz, JHEP 0206, 043 (2002) [hepth/0205236].
 [25] S. A. Hartnoll, J. Polchinski, E. Silverstein and D. Tong, “Towards strange metallic holography,” JHEP 1004, 120 (2010) [arXiv:0912.1061 [hepth]].
 [26] F. Benini, F. Canoura, S. Cremonesi, C. Nunez and A. V. Ramallo, JHEP 0702, 090 (2007) [hepth/0612118]; F. Bigazzi, A. L. Cotrone and A. Paredes, “KlebanovWitten theory with massive dynamical flavors,” JHEP 0809, 048 (2008) [arXiv:0807.0298 [hepth]].
 [27] F. Bigazzi, A. L. Cotrone, J. Mas, D. Mayerson and J. Tarrio, “D3D7 QuarkGluon Plasmas at Finite Baryon Density,” JHEP 1104, 060 (2011) [arXiv:1101.3560 [hepth]].
 [28] F. Bigazzi, A. L. Cotrone and J. Tarrio, “Charged D3D7 plasmas: novel solutions, extremality and stability issues,” arXiv:1304.4802 [hepth].
 [29] C. Nunez, A. Paredes and A. V. Ramallo, “Unquenched Flavor in the Gauge/Gravity Correspondence,” Adv. High Energy Phys. 2010, 196714 (2010) [arXiv:1002.1088 [hepth]].
 [30] K. Goldstein, S. Kachru, S. Prakash and S. P. Trivedi, “Holography of Charged Dilaton Black Holes,” JHEP 1008, 078 (2010) [arXiv:0911.3586 [hepth]].
 [31] K. Goldstein, N. Iizuka, S. Kachru, S. Prakash, S. P. Trivedi and A. Westphal, “Holography of Dyonic Dilaton Black Branes,” JHEP 1010, 027 (2010) [arXiv:1007.2490 [hepth]].
 [32] C. Charmousis, B. Gouteraux, B. S. Kim, E. Kiritsis and R. Meyer, “Effective Holographic Theories for lowtemperature condensed matter systems,” JHEP 1011, 151 (2010) [arXiv:1005.4690 [hepth]]; B. Gouteraux and E. Kiritsis, “Generalized Holographic Quantum Criticality at Finite Density,” JHEP 1112, 036 (2011) [arXiv:1107.2116 [hepth]].
 [33] N. Drukker and B. Fiol, “Allgenus calculation of Wilson loops using Dbranes,” JHEP 0502, 010 (2005) [hepth/0501109].
 [34] S. Yamaguchi, “Wilson loops of antisymmetric representation and D5branes,” JHEP 0605, 037 (2006) [hepth/0603208].
 [35] J. Gomis and F. Passerini, “Holographic Wilson Loops,” JHEP 0608, 074 (2006) [hepth/0604007].
 [36] S. A. Hartnoll and S. P. Kumar, “Higher rank Wilson loops from a matrix model,” JHEP 0608, 026 (2006) [hepth/0605027].
 [37] E. D’Hoker, J. Estes and M. Gutperle, “Gravity duals of halfBPS Wilson loops,” JHEP 0706, 063 (2007) [arXiv:0705.1004 [hepth]].
 [38] N. Drukker and V. Forini, “Generalized quarkantiquark potential at weak and strong coupling,” JHEP 1106, 131 (2011) [arXiv:1105.5144 [hepth]].
 [39] D. Marolf and A. W. Peet, “Brane baldness versus superselection sectors,” Phys. Rev. D 60, 105007 (1999) [hepth/9903213].
 [40] N. ArkaniHamed, A. G. Cohen, D. B. Kaplan, A. Karch and L. Motl, “Deconstructing (2,0) and little string theories,” JHEP 0301, 083 (2003) [hepth/0110146].
 [41] C. Hoyos, “Higher dimensional conformal field theories in the Coulomb branch,” Phys. Lett. B 696, 145 (2011) [arXiv:1010.4438 [hepth]].
 [42] S. Ryu and T. Takayanagi, “Holographic derivation of entanglement entropy from AdS/CFT,” Phys. Rev. Lett. 96, 181602 (2006) [hepth/0603001]; S. Ryu and T. Takayanagi, “Aspects of Holographic Entanglement Entropy,” JHEP 0608, 045 (2006) [hepth/0605073].
 [43] S. Arapoglu, N. S. Deger and A. Kaya, “A Note on supergravity solutions for partially localized intersecting branes,” Phys. Lett. B 578, 203 (2004) [hepth/0306040].
 [44] S. Slavyanov and W. Lay, “ Special Functions: A Unified Theory based on Singularities”, Oxford University Press, 2000.
 [45] E. I. Guendelman and A. Rabinowitz, “The Gravitational field of a hedgehog and the evolution of vacuum bubbles,” Phys. Rev. D 44, 3152 (1991).
 [46] M. Headrick, “Hedgehog black holes and the Polyakov loop at strong coupling,” Phys. Rev. D 77, 105017 (2008) [arXiv:0712.4155 [hepth]].
 [47] G. T. Horowitz and B. Way, “Lifshitz Singularities,” Phys. Rev. D 85, 046008 (2012) [arXiv:1111.1243 [hepth]]. K. Copsey and R. Mann, “Pathologies in Asymptotically Lifshitz Spacetimes,” JHEP 1103, 039 (2011) [arXiv:1011.3502 [hepth]]. K. Copsey and R. Mann, “Singularities in Hyperscaling Violating Spacetimes,” JHEP 1304, 079 (2013) [arXiv:1210.1231 [hepth]].
 [48] S. Harrison, S. Kachru and H. Wang, “Resolving Lifshitz Horizons,” arXiv:1202.6635 [hepth].
 [49] J. Bhattacharya, S. Cremonini and A. Sinkovics, JHEP 1302, 147 (2013) [arXiv:1208.1752 [hepth]].
 [50] H. Y. Chen, K. Hashimoto and S. Matsuura, JHEP 1002, 104 (2010) [arXiv:0909.1296 [hepth]].
 [51] P. Kraus, F. Larsen and S. P. Trivedi, “The Coulomb branch of gauge theory from rotating branes,” JHEP 9903, 003 (1999) [hepth/9811120].
 [52] A. Brandhuber and K. Sfetsos, “Wilson loops from multicenter and rotating branes, mass gaps and phase structure in gauge theories,” Adv. Theor. Math. Phys. 3, 851 (1999) [hepth/9906201].