Abstract
Using highprecision numerical analysis, we show that dimensional gauge theories holographically dual to dimensional EinsteinMaxwellChernSimons theory undergo a quantum phase transition in the presence of a finite charge density and magnetic field. The quantum critical theory has dynamical scaling exponent , and is reached by tuning a relevant operator of scaling dimension . For magnetic field above the critical value , the system behaves as a Fermi liquid. As the magnetic field approaches from the high field side, the specific heat coefficient diverges as , and nonFermi liquid behavior sets in. For the entropy density becomes nonvanishing at zero temperature, and scales according to . At , and for small nonzero temperature , a new scaling law sets in for which . Throughout a small region surrounding the quantum critical point, the ratio is given by a universal scaling function which depends only on the ratio .
The quantum phase transition involves nonanalytic behavior of the specific heat and magnetization but no change of symmetry. Above the critical field, our numerical results are consistent with those predicted by the Hertz/Millis theory applied to metamagnetic quantum phase transitions, which also describe nonanalytic changes in magnetization without change of symmetry. Such transitions have been the subject of much experimental investigation recently, especially in the compound SrRuO, and we comment on the connections.
July 4, 2019
Holographic Metamagnetism, Quantum Criticality,
and Crossover Behavior ^{1}^{1}1This work was supported in part by NSF grant PHY0757702.
Eric D’Hoker and Per Kraus
Department of Physics and Astronomy
University of California, Los Angeles, CA 90095, USA
dhoker@physics.ucla.edu; pkraus@ucla.edu
1 Introduction and summary of results
The AdS/CFT correspondence provides a precise and powerful tool for the study of thermodynamics, statistical mechanics, and transport properties in a variety of 4dimensional gauge theories at finite temperature, charge density, and magnetic field. In the large and large ‘t Hooft coupling limits these gauge theories are holographically dual to certain black brane solutions (with both electric and magnetic charges) to 5dimensional EinsteinMaxwell theory. This theory includes a ChernSimons term, whose coupling captures the strength of the chiral anomaly in the Maxwell current. For the special value (in our conventions) the bulk EinsteinMaxwell theory is a consistent supersymmetric truncation of Type IIB or Mtheory [1, 2, 3], while for not equal to this value, supersymmetry is lost.
The study of thermodynamics and transport properties in 4dimensional YangMills theory is of direct interest to the physics of heavy ion collisions at RHIC (and soon at the LHC), where quarks and gluons are subject to high temperatures, strong magnetic fields, and large currents. Of great interest as well is the possibility of studying novel phases of finite density matter at low temperatures. By varying the charge density and magnetic field one can search for signals of zero temperature quantum phase transitions [4], as is done experimentally in, for example, the heavy fermion compounds and high temperature superconductors.
The thermodynamics of the relevant solutions to 5dimensional EinsteinMaxwell theory can be understood analytically in two limiting cases. In the absence of magnetic fields, the relevant supergravity solution is the electrically charged (ReissnerNordstrom) black brane in . At high temperature, its entropy density scales as , as dictated by scale invariance in the UV, and in accord with weak coupling YangMills theory. As tends to 0, the black brane tends to its extremal form with nearhorizon geometry, and nonvanishing entropy density. On the one hand, this AdS factor is needed in “semiholographic” models of nonFermi liquid behavior [5, 6, 7, 8] (see also [9]). On the other hand the ground state entropy density is rather exotic from the Fermi surface perspective. It is exotic also from the point of view of the CFTs that arise in the AdS/CFT duality, which typically contain massless charged bosons that would be expected to condense and lead to a unique ground state.
The other limiting case consists of vanishing electric charge density and finite magnetic field . At high temperature the entropy density again scales as , but at low temperature the analysis of [10] showed that . The overall numerical coefficient in this linear scaling law was computed analytically by taking advantage of the near horizon AdS factor that emerges in this regime. On the CFT side, the low temperature physics was seen to be controlled by a gas of fermions arising from the lowest Landau level of the 4dimensional gauge theory in the presence of a magnetic field.
The general case of nonzero and was studied in [11], and the low temperature thermodynamics were found to depend crucially on the value of the ChernSimons coupling . This was characterized in terms of flows in parameter space as the temperature was lowered. These flows head towards three distinct fixed points, depending on the value of . For , the geometry near the horizon can be thought of as deformed by the presence of the magnetic field. The entropy density was found to be nonvanishing at . Precisely at , the solutions support a nearhorizon warped AdS factor (the warped solutions were studied in the context of topologically massive gravity in [12, 13]), and the entropy density is again nonvanishing at extremality. For (which includes the supersymmetric value ), the presence of a moderate strength (or larger) magnetic field was seen to lead to a precipitous drop in the entropy density as the temperature was lowered. However, our numerics were found to break down in the combined regions of very low temperatures and small magnetic fields, presumably as a result of one of our choices of gauge and, as a result, our analysis stopped short of obtaining reliable data for the entropy density at low over the full range of magnetic fields.
1.1 Summary of results
In the present paper, we shall apply a simple remedy to the gauge choice problem which plagued the low temperature numerical work of [11], and carry out highprecision numerical analyses down to ultralow temperatures, for a wide range of values of .
By virtue of scale invariance, only dimensionless combinations of quantities, such as , afford any intrinsic physical meaning. Thus, we shall introduce normalized, dimensionless, magnetic field , temperature , and entropy density via the following relations,
(1.1) 
The results are displayed schematically in Figure 1 and summarized below.

Our main result is that a continuous quantum phase transition occurs at
(1.2) 
For , the entropy density goes to zero linearly with temperature, , reflecting the same Fermi liquid physics that was seen in the large limit. Here, however, the linear behavior is occurring at finite charge density where the result does not follow from conformal invariance. Note that the specific heat , at constant , has the same behavior as the entropy density, since .

On approaching from the large side, the coefficient of the linear term in the versus relation is found to diverge according to,
(1.3) This signals a breakdown of Fermi liquid behavior.

Approaching the quantum critical point along the temperature axis at fixed , the entropy density exhibits a new power law scaling,
(1.4) 
For the zero temperature entropy density is found to be nonzero. Near it obeys the scaling,
(1.5) 
The entropy density in the vicinity of the fixed point can be expressed in terms of a scaling function as,
(1.6) This is to be contrasted with regions far from the critical point, where the entropy density is a nontrivial function of two dimensionless combinations of , and .
From our results, we can infer that the quantum critical theory lives in spacetime dimensions (that is, there are no long range correlations in the remaining 2 spatial directions along the boundary), has dynamical critical exponent (since the value found numerically for is consistent with ), and has a relevant operator with scaling dimension . The relevant operator corresponds to a change of away from .
1.2 Comparison with known quantum critical systems
It is illuminating to place these results within the context of known quantum critical systems. Although thermodynamic quantities such as the specific heat and magnetization behave in a nonanalytic fashion across the phase transition, we note that there is no change of symmetry associated with the transition. A finite temperature example of such behavior is the liquidgas transition in water. A zero temperature version is a metamagnetic quantum critical point [14, 15], whose behavior closely parallels our system (a finite temperature version of a holographic metamagnetic phase transition in the D4D8 system was studied in [16]). Metamagnetism refers to a sharp change in the magnetization of a material as an external magnetic field is tuned through some nonzero value. At finite temperature, metamagnetic transitions typically consist of first order lines terminating at a second order critical point, as in the liquidvapor case. If the critical point can be brought to zero temperature by adjusting some control parameter, one obtains a metamagnetic quantum critical point.
The standard approach to such a quantum critical point is based on the Hertz/Millis theory [17, 18, 19], for which the effective action in momentum space is,
(1.7) 
The real bosonic field represents the local magnetization, and the above action can be obtained by integrating out fermions at 1loop. Under scale transformations acting as , we see that should be assigned scale dimension 3, and hence . Similarly, has scale dimensions . These assignments match what we found for our system, which furthermore corresponds to , since the Landau level quantization only allows low energy modes to propagate parallel to .
The action (1.7) is only meant to be applied for . Indeed, the behavior of our system in the region , with its nonzero ground state entropy density, cannot be described by this action alone.
Metamagnetic quantum criticality in the compound SrRuO has been the subject of extensive experimental investigation in the past few years (see [15] and references therein), and we will comment on this connection in Section 3.8.
A large number of other AdS/CFT examples undergoing phase transitions have been studied in the literature, both continuous and discontinuous, and at finite and zero temperature; for example [20, 21, 22, 23, 24, 25, 26, 27, 28]. However, our setup seems particularly attractive and is nicely related to real experimental systems. In particular, unlike other examples of quantum phase transitions we do not have to add any extra ingredients in the way of scalar fields or probe branes. We employ only a metric and an Abelian gauge field, with an EinsteinMaxwellChernSimons action that is known to describe all supersymmetric AdS theories related by compactification of Type IIB or Mtheory [1, 2, 3]. Thus our framework is both simple and universal.
The remainder of this paper is organized as follows. In Section 2, we spell out the setup of the holographic calculations, including the specification of initial data at the horizon, asymptotic data at the boundary, and the construction of regular gauge choices. In Section 3, we present the results of highprecision numerical solutions to the reduced EinsteinMaxwellChernSimons equations, identify the quantum critical point, and critical behavior in its vicinity. We also compare our results with those from the effective Hertz/Millis theory, and comment on the connection with metamagnetic quantum criticality observed in real materials like SrRuO. A discussion of the results and open avenues for future research is given in Section 4.
2 Holographic calculations
In this section we spell out some of the technical details of our computations. Results of the numerical calculations will be presented in section 3, and the impatient reader may wish to jump there.
The starting point for our holographic calculations is 5dimensional EinsteinMaxwell theory with a ChernSimons term. Throughout this paper, the ChernSimons coefficient will be considered fixed at its supersymmetric value . A detailed discussion of the action, including boundary terms needed for holographic renormalization, and the construction of the boundary current and stress tensor may be found in [11]. Here, we shall limit our discussion to the field equations and the asymptotic behavior of the fields near the horizon, and at the asymptotic boundary. The EinsteinMaxwell field equations are given by and,
(2.1) 
Uniformity and constancy in time of the magnetic field and the charge density allow us to restrict to a spacetime translation invariant Ansatz, given by,
(2.2) 
The functions depend only on the radial coordinate , while the magnetic field is constant by the Bianchi identity. A gauge choice has been made here for the coordinate in order to put the Ansatz in canonical form with matching coefficients of its first two terms in . The reduced field equations were given in [11].
2.1 Data at the horizon
The reduced field equations are to be solved subject to regularity conditions at the horizon and at the asymptotic boundary . For the purpose of numerical analysis, it will be convenient to parametrize solutions in terms of data at the horizon which satisfy the regularity conditions at the horizon. Regularity of the full solution, including at the asymptotic boundary, must then be verified numerically for each set of data.
We begin by spelling out the data at the horizon. (This discussion will parallel the one presented in [11], but there will be important differences motivated by the need to remedy the gauge choice problems alluded to in the Introduction.) The (outer) horizon at , and the Hawking temperature are defined by,
(2.3) 
Our numerical analysis will always be carried out at , even though may become very small; thus, we are free to rescale , and set . By rescaling also , we may set . Invariance of the Ansatz under symmetry, (under which , , and ), allows us to set . With these choices, the fields at the horizon take the form,
(2.4) 
The reduced field equations of [11] for , combined with the requirement of regularity at the horizon, dictate certain relations amongst the remaining data at the horizon, namely , , and . They are given by (for ),
(2.5) 
Given , the data and are related to one another by the first equation of (2.1). Keeping either or as an independent free parameter, the remaining initial data and are uniquely determined by the last two equations of (2.1).
2.2 Gauge fixing and regularity
The parametrization of the horizon data presented above is more general than the one given in [11], since here the parameter is kept unspecified, while was set to 0 in [11]. The argument invoked to set was the covariance of the Ansatz of (2) under boosts in the direction. Under a boost by velocity , the spacetime coordinates transform as usual, and with , and where cannot exceed the speed of light, . The transformation under boosts of the functions , however, must be accompanied by a transformation of the holographic coordinate , which is required to restore the boosted Ansatz back to the canonical form of (2). As a result, the Maxwell fields transform as follows,
(2.6) 
where . The ratio transforms as . Clearly, provided , a boost by may be used to set to zero. The problem, however, is that the coordinate transformation required to accompany this boost may be singular on some of the functions . For example, the transformation law under a boost by velocity of the function is given by,
(2.7) 
We have verified numerically that, in the region of low where our numerics were found to break down in [11], the denominator in (2.7) indeed crosses zero as is increased away from the horizon at , thereby rendering this boost transformation singular.
The remedy to this problem is simple: the parameter should be left unspecified, thereby eliminating the need to perform boost transformations on the fields. It will be convenient to parametrize solutions by the values , so that the remaining initial data at the horizon, , and are uniquely determined by (2.1). Not every assignment of will produce a regular solution. Also, two regular solutions may be related to one another by a regular boost, and thus be physically equivalent. The parameter space of all regular solutions may be described as follows. To every pair we assign the maximal interval on the real line such that for every value of , the solution specified by the parameters is regular. The end points of the interval correspond to the limits where the velocity of the solution tends to the speed of light. When computing boost invariant physical quantities, such as the magnetic field , and the rest frame temperature , charge density , and entropy density , the datum may be chosen to be any value in the interval .
2.3 Data at the asymptotic boundary
For regular solutions the functions have the following asymptotic behavior as , (keeping only leading contributions),
(2.8) 
In these coordinates, the conformal boundary metric is , and the solution’s velocity is . Rescaling by , and by , while combining an transformation with a boost by a velocity , restores the coordinates to the standard Minkowski metric, and yields the following expressions for the physical magnetic field , temperature , charge density , and entropy density in the rest frame of the solution,^{2}^{2}2In the CFT, the rest frame corresponds to a statistical ensemble weighted by the Boltzmann factor . Boosting produces additional chemical potentials multiplying momenta and currents. (see [11] for derivations),
(2.9) 
Here, , and the normalized entropy density is defined as , where is the total entropy in volume Vol, and is the 5dimensional Newton constant.
By virtue of scale invariance, only dimensionless combinations of quantities, such as , afford any intrinsic physical meaning. Thus, the results on phase transitions and associated critical exponents at critical points, aimed for in this paper, will all be derived by evaluating the dependence between the dimensionless physical quantities and , for various (fixed) values of , which were all defined in (1.1).
2.4 Numerical fine points
The initial data for any regular solution is a pair and a value , where the interval was defined so that the data produce a regular solution. For certain values of , the interval may be empty. For , numerical analysis yields the following bound,
(2.10) 
The precise form of the critical curve in the plane inside of which is not known analytically, but may be obtained numerically, as was done in [11].
Away from the low temperature regime the interval is sufficiently large that it is possible, at least for , to make a single uniform choice for and still cover most of the parameter space. A convenient choice is . However, as the temperature is lowered shrinks, and one is forced to tune to greater precision. This effect becomes especially pronounced in the region of low magnetic field where a ground state entropy develops. However, by specifying more general tunings for as we vary we are able to fully cover the low temperature region.
The approach to ultralow temperatures, which is needed in various parts of our numerical work, requires a high degree of finetuning of the horizon initial data and the gauge choice . It also requires evaluating the asymptotic data, such as and at large values of , which we typically have taken to range from to . With such high degrees of finetuning, and extended ranges of integration, the issue of numerical accuracy and numerical stability of the calculations becomes of utmost importance. Our calculations were performed with 15 to 20 digits accuracy, and the ODEs were solved with absolute and relative error tolerances ranging from to for the lowest temperatures. Stability of the results was checked versus changing the asymptotic value of , the absolute and relative error tolerances, and the number of digits.
3 Numerical results
In this section, we present our numerical results, organized as a function of the magnitude of the magnetic field , starting at large .
3.1 Large regime
In [11] the case (quoted there as ) was studied at the supersymmetric value . As displayed in Figure 3 of [11] the low temperature entropy density was seen to drop well below its value, and appeared to be heading towards zero. However, the numerics broke down at temperature , and so did not allow for an exploration of ultralow temperatures. By adjusting as we lower the temperature, we can do much better, as shown in Figure 2 below. The entropy density clearly vanishes at zero temperature. A linear behavior, is evident in the approach to zero temperature, as shown in the right panel of Figure 2.
By repeating this analysis for other sufficiently large values of , we find similar behavior: the entropy density vanishes at zero temperature, and does so linearly at very low temperatures. This linear behavior is characteristic of Fermi liquids. The textbook intuitive explanation for the linear dependence is that at low temperature the smoothed out step function form of the FermiDirac distribution implies that only electrons within energy of the Fermi energy contribute. We indeed expect our system to be described by a theory of light fermions in this large , low regime, since the magnetic field raises up all energies except for the lowest fermion Landau level. Our numerical results are a pleasing confirmation of this intuition.
In the limit it was found in [10] that a near horizon AdS factor emerged at low energies, and the resulting dimensional conformal invariance could be used to explain the behavior. This is no longer true away from this limit, as the charge density introduces an additional scale into the problem, and from the numerics we can see that the near horizon geometry becomes deformed away from AdS.
3.2 Approaching
Experimentally, a breakdown of Fermi liquid behavior upon tuning an external magnetic field can often be seen in a divergence of the specific heat coefficient, defined as,
(3.1) 
Since , we can equivalently write at low temperatures. In our system, stabilizes at a constant value for large , but is seen to diverge at a critical value . Numerically, the critical value is found to be bounded as follows,
(3.2) 
which results in the the value quoted already in (1.2), namely . To characterize the divergence we display a plot of versus , as in Figure 3. The straight line shows that the low temperature entropy behaves as
(3.3) 
as for fixed near (but larger than) . Upon approaching we find that the linear regime with is confined to an ever smaller low temperature window.
3.3 Scaling at the critical magnetic field
Next we set and again study the low temperature behavior of the entropy density. We find a new scaling law, , which numerically extends over at least four orders of magnitude in temperature, , as shown in Figures 4. This nontrivial power law manifestly represents nonFermi liquid behavior, analogous to what is seen in real materials at a quantum critical point.
Towards ultralow temperatures, the numerical behavior of will ultimately turn over to be linear in (for ) or to be a nonzero constant for (see the next subsection). This deviation from scaling is caused by the fact that the value of is known only numerically from (1.2), and we are never able to sit at precisely the value of .
3.4 The quantum critical region and crossover
The scaling law extends to the vicinity of the critical magnetic field , but the range over which this scaling law holds shrinks as is increased. The corresponding behavior of is illustrated in Figure 5.
First, on the right panel of Figure 5, is plotted versus for . The scaling law exhibited in the preceding subsection, is clearly recovered here. At sufficiently high temperature, the behavior eventually crosses over to the dependence controlled by the UV theory. This is clearly shown on the right panel of Figure 5, where the crossover region may be identified with the temperature interval .
Second, on the left panel of Figure 5, the flows of as a function of at various fixed values of are shown. Curves a, b, c, d, e, and f clearly exhibit the turnover from scaling behavior to linear behavior at ultralow temperatures. Given this crossover behavior, we know that the corresponding values of must all be above the critical magnetic field , with value closest to critical corresponding to curve f with . Curve g behaves completely differently at ultralow temperatures, and is seen to tend towards a nonzero constant. Given this behavior, we know that the corresponding value of must be below critical. Combining both results gives , as announced in (3.2), and (1.2).
Finally, as is increased, the region over which the scaling law holds becomes smaller and smaller, being overtaken by the law at high joining the law at low , and ultimately will have disappeared altogether by the time is reached (this result is not shown in the figures).
3.5 Low region
Decreasing the magnetic field further, namely , we find a completely different low temperature behavior: the entropy now stabilizes to a nonzero value at . We plot this limiting value in Figure 6. In [11] we speculated that an infinitesimally small magnetic field might be sufficient to remove the ground state entropy of the solution. We now see that this speculation is incorrect – only magnetic fields accomplish this.
The extremal entropy turns on continuously below the critical magnetic field. A curve that passes through all the points to better than accuracy is given by
(3.4) 
The prefactor of is fixed by the value for the purely electrically charged black brane, which is of course known analytically. It is remarkable that this function fits the data so well, but we are unable to say whether it represents an exact analytical result.
To excellent accuracy, and as exhibited by (3.4), we find
(3.5) 
near the critical point. If we increase the temperature at fixed magnetic field in this regime, we again find a region obeying the same scaling as described in section 3.3.
3.6 Scaling region near the critical point
The fact that the entropy density at the critical magnetic field has a simple power law dependence on temperature is indicative of a quantum critical point. Recall that since we are working in terms of dimensionless quantities as defined by the scalings in the asymptotically AdS region, a priori the dimensionless entropy density is allowed to be an arbitrary function of the dimensionless temperature . This interpretation can be sharpened further by writing down a scaling form for the entropy in the vicinity of the critical point
(3.6) 
where, according to our results, the scaling function has asymptotic behavior
(3.7) 
for some constants . In general, the entropy is a function of two dimensionless combinations of , , and , whereas near the critical point the claim is that it can be written as a function of only one variable. All of this is of course standard from the general theory of classical and quantum critical phenomena.
It is further useful to recall that for a quantum critical theory in spatial dimensions, with dynamical critical exponent , and with a relevant coupling of scale dimension , the entropy density will take the scaling form
(3.8) 
Comparing to (3.7) we read off , , and .
We have verified that the scaling form (3.7) conforms to our numerical results near the critical point at and . From the numerics we can reconstruct the form of the scaling function, which is shown in Figure 7.
We can use this to determine the constants appearing in (3.7) to be,
(3.9) 
As noted in the Introduction, the standard approach to modelling a magnetically tuned quantum critical point of the type that we are seeing is based on the Hertz/Millis theory. We consider the action in (1.7) with . To the extent that this action captures the low energy degrees of freedom of our theory, it is natural to compute its finite temperature entropy density and compare to our results. This cannot be entirely correct for a number of reasons, not least that there is no way of explaining the ground state entropy density in this framework, but the comparison is instructive nonetheless. In the free field limit we just need the dispersion relation implied by (1.7), which is
(3.10) 
The partition function is
(3.11) 
from which we can extract the entropy density as
(3.12) 
We can write the result in a scaling form as
(3.13) 
and then compare with the scaling function coming from gravity. The comparison is shown in the right panel of Figure 7. We only compare in the region , since it is only in this region that the Hertz/Millis action applies. Nothing fixes the normalization of the magnetic field appearing in (1.7) compared to that in gravity, and so we have allowed ourselves to adjust this relationship to achieve a good fit. Similarly, we have introduced a parameter to adjust the overall normalizations. The functions match surprisingly well, although it is unclear how much significance to attach to this.
3.7 The quantum critical point
Initial data yielding our closest approach to the quantum critical point is given by
(3.14) 
This choice of along with the gauge choice leads to the following values for the asymptotic parameters of (2.3),
(3.15) 
from which we can compute the temperature and entropy density as
(3.16) 
From (3.14) we see that there is a large hierarchy between the size of the electric and magnetic charges at the horizon in our chosen coordinates. Nevertheless, measured at infinity in the rest frame the ratio of these quantities is of order unity, . This happens because there are large rescalings that occur in the process of integrating out from the horizon to the asymptotic regions, as illustrated by the magnitude of the parameters in (3.15). Using (2.9) these can convert a large hierarchy into one of order unity.
It is instructive to examine the metric functions for this near critical solution. In Figure 8 we provide plots of these functions in the near (and not so near) horizon region. Near the horizon, one can think of these metric functions as representing a deformation of the configuration
(3.17) 
This represents AdS in “boosted coordinates” with magnetic flux. We expect that, as we get nearer and nearer to the critical point, the metric functions will approach (3.17) to greater and greater accuracy. On the other hand, it is crucial that there always be some deviation away from (3.17), as this solution is purely magnetic with vanishing electric charge density. Presumably what happens is that there is a finetuned limit in which all of the charge is carried by the bulk supergravity fields outside the near horizon region; this is possible by virtue of the ChernSimons term. It would of course be extremely useful if this limiting solution could be found analytically.
We close this subsection with an issue which has not yet been conclusively resolved by our numerical studies. For large , all flows towards lower end in the purely magnetic fixed point solution of [10], whose nearhorizon geometry is . More precisely, it can be shown numerically that all flows for end up at the purely magnetic fixed point. Where do the flows for end? Numerically, we have established that flows for less than but close to come near the purely magnetic fixed point, but at temperatures so low (on the order of ) that further numerical analysis can be pursued only at the cost of calculation times and memory which exceed those of our computers. Thus, the logical possibility that the flows with end up on the critical curve, before the purely magnetic fixed point has been attained, cannot be excluded at this time. If this situation is in fact the one that occurs, then it would have a surprising similarity to the dynamics of the theory, derived in [11].
3.8 Metamagnetic quantum criticality in SrRuO
We close this section with a few words about a real system that parallels ours in a number of ways. Quantum criticality and novel phases in SrRuO have been the subject of much experimental and theoretical interest in the past few years; e.g. [15, 34]. SrRuO is a layered structure, which for a large magnetic field perpendicular to the layers exhibits a line of first order metamagnetic phase transitions at finite temperature, ending at a finite temperature critical point. By including a component of magnetic field in the plane of the layers, the critical point can be brought to zero temperature. As in our case, the transition occurs at finite magnetic field, and involves no change of symmetry. Away from the critical point the system behaves as a Fermi liquid, with entropy density linear in temperature. As the critical point is approached, the linear term diverges as , just as we found (in this case is approximately Tesla.) In recent work, the complete “entropy landscape” of SrRuO at finite temperature and magnetic field has been mapped out [15]. In very pure samples, as one tries to sit right on top of the critical point one finds instead that a new phase emerges, which is believed to be a spatially anisotropic nematic phase [33, 34]. This has been described as nature’s solution to the problem of avoiding a finite entropy density at zero temperature. The parallels with our system are evident (though an obvious difference is that our critical theory effectively sees only a single spatial direction, compared to the two inplane directions in SrRuO), and lead us to speculate whether in our case the extremal black hole phase is unstable and gives way to a new phase with reduced symmetry, as in [35].
4 Discussion
We have found a holographic description of a quantum critical point, reached by tuning a magnetic field at finite density, that nicely resembles examples seen in the real world. In particular, approached from the high field side, we have Fermi liquid behavior; as the magnetic field is lowered the specific heat coefficient diverges, and we enter a regime of nonFermi liquid behavior with nontrivial scaling properties.
The most exotic property as compared to known physical examples is the existence of a ground state entropy on the low field side of the transition. The presence of a nonvanishing entropy at for would again seem to contradict the third “law” of thermodynamics [11], as it did in the absence of magnetic fields. Actually, our results show that the lifting of the ground state entropy as a function of an external magnetic field is a subtle and dynamical issue. Also, we expect the ground state entropy to be lifted by instabilities driving the system towards inhomogeneity, or by turning on further fields besides an external magnetic field.
It is worth emphasizing again the universality of our results: they apply to all supersymmetric AdS examples arising from IIB or Mtheory, since all such theories admit a consistent truncation to the EinsteinMaxwellChernSimons theory used here [1, 2, 3]. We did not have to introduce any model building devices in the way of probe branes or scalar fields.
There are many open questions and avenues for further investigation. Many of our numerical results cry out for an analytical derivation. In particular, one would expect to be able to derive the value of the dynamical exponent, and the dimension of the relevant operator governing the critical theory. It may similarly be possible to understand these results microscopically on the gauge theory side.
All of our results were presented for the supersymmetric value of the ChernSimons coupling, , but it would be useful to consider other values as well. Our expectation is that as is increased the critical point will retain its character but with moving to a smaller value. An interesting question would then be whether we reach at finite , for if so there would no longer be an extremal black hole phase.
In this work we have only considered the thermodynamics, but the existence of the critical theory implies scaling behavior of correlation functions with respect to frequency and momentum. Computing these would be valuable in pinning down the precise connection to the Hertz/Millis theory.
Finally, to more closely model real systems it would be very interesting to construct a version of our system giving rise to critical behavior in two or three spatial dimensions.
Acknowledgments
It is a pleasure to thank David Berenstein, Sudip Chakravarty, Tom Faulkner, Gary Horowitz, Don Marolf, Eric Perlmutter, Joe Polchinski, and Matt Roberts for helpful discussions.
References
 [1] A. Buchel and J. T. Liu, “Gauged supergravity from type IIB string theory on Y(p,q) manifolds,” Nucl. Phys. B 771 (2007) 93 [arXiv:hepth/0608002].
 [2] J. P. Gauntlett, E. O Colgain and O. Varela, “Properties of some conformal field theories with Mtheory duals,” JHEP 0702 (2007) 049 [arXiv:hepth/0611219].
 [3] J. P. Gauntlett and O. Varela, “Consistent KaluzaKlein Reductions for General Supersymmetric AdS Solutions,” Phys. Rev. D 76 (2007) 126007 [arXiv:0707.2315 [hepth]].
 [4] S. Sachdev, “Quantum Phase Transitions”, Cambridge University Press (2001)
 [5] H. Liu, J. McGreevy and D. Vegh, “NonFermi liquids from holography,” arXiv:0903.2477 [hepth].
 [6] M. Cubrovic, J. Zaanen and K. Schalm, “Fermions and the AdS/CFT correspondence: quantum phase transitions and the emergent Fermiliquid,” arXiv:0904.1993 [hepth].
 [7] T. Faulkner, H. Liu, J. McGreevy and D. Vegh, “Emergent quantum criticality, Fermi surfaces, and AdS2,” arXiv:0907.2694 [hepth].
 [8] T. Faulkner and J. Polchinski, “SemiHolographic Fermi Liquids,” arXiv:1001.5049 [hepth].
 [9] S. J. Rey, “String Theory On Thin Semiconductors: Holographic Realization Of Fermi Points And Surfaces,” Prog. Theor. Phys. Suppl. 177, 128 (2009) [arXiv:0911.5295 [hepth]].
 [10] E. D’Hoker and P. Kraus, “Magnetic Brane Solutions in AdS,” JHEP 0910, 088 (2009) [arXiv:0908.3875 [hepth]].
 [11] E. D’Hoker and P. Kraus, “Charged Magnetic Brane Solutions in and the fate of the third law of thermodynamics,” arXiv:0911.4518 [hepth].
 [12] D. Anninos, W. Li, M. Padi, W. Song and A. Strominger, “Warped AdS Black Holes,” JHEP 0903, 130 (2009) [arXiv:0807.3040 [hepth]].
 [13] G. Compere and S. Detournay, “Boundary conditions for spacelike and timelike warped AdS spaces in topologically massive gravity,” JHEP 0908, 092 (2009) [arXiv:0906.1243 [hepth]].
 [14] A. J. Millis, A. J. Schofield, G. G. Lonzarich and S. A. Grigera, “Metamagnetic quantum criticality in metals,” Phys. Rev. Lett. 88, 217204 (2002).
 [15] A. W. Rost, R. S. Perry, J.F. Mercure, A. P. Mackenzie, and S. A. Grigera “Entropy Landscape of Phase Formation Associated with Quantum Criticality in ”, Science Vol. 325. no. 5946, pp. 1360  1363 (2009)
 [16] G. Lifschytz and M. Lippert, “Holographic Magnetic Phase Transition,” Phys. Rev. D 80, 066007 (2009) [arXiv:0906.3892 [hepth]].
 [17] J. A. Hertz, “Quantum critical phenomena,” Phys. Rev. B 14 1165 (1976).
 [18] A. J. Millis, “Effect of a nonzero temperature on quantum critical points in itinerant fermion systems”, Phys. Rev. B 48, 7183â7196 (1993)
 [19] H. v. Lohneysen, A. Rosch, M. Mojta, and P. Wolfle, “Fermiliquid instabilities at magnetic quantum phase transitions”, Rev. Mod. Phys. 79 , 1015â1075 (2007)
 [20] E. Witten, “Antide Sitter space, thermal phase transition, and confinement in gauge theories,” Adv. Theor. Math. Phys. 2, 505 (1998) [arXiv:hepth/9803131].
 [21] A. Chamblin, R. Emparan, C. V. Johnson and R. C. Myers, “Charged AdS black holes and catastrophic holography,” Phys. Rev. D 60, 064018 (1999) [arXiv:hepth/9902170].
 [22] A. Parnachev and D. A. Sahakyan, “Chiral phase transition from string theory,” Phys. Rev. Lett. 97, 111601 (2006) [arXiv:hepth/0604173].
 [23] D. Mateos, R. C. Myers and R. M. Thomson, “Holographic phase transitions with fundamental matter,” Phys. Rev. Lett. 97, 091601 (2006) [arXiv:hepth/0605046].
 [24] J. L. Davis, M. Gutperle, P. Kraus and I. Sachs, “Stringy NJL and GrossNeveu models at finite density and temperature,” JHEP 0710, 049 (2007) [arXiv:0708.0589 [hepth]].
 [25] S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, “Building a Holographic Superconductor,” Phys. Rev. Lett. 101, 031601 (2008) [arXiv:0803.3295 [hepth]].
 [26] J. L. Davis, P. Kraus and A. Shah, “Gravity Dual of a Quantum Hall Plateau Transition,” JHEP 0811, 020 (2008) [arXiv:0809.1876 [hepth]].
 [27] N. Evans, A. Gebauer, K. Y. Kim and M. Magou, “Holographic Description of the Phase Diagram of a Chiral Symmetry Breaking Gauge Theory,” arXiv:1002.1885 [hepth].
 [28] K. Jensen, A. Karch and E. G. Thompson, “A Holographic Quantum Critical Point at Finite Magnetic Field and Finite Density,” arXiv:1002.2447 [hepth].
 [29] G. T. Horowitz and M. M. Roberts, “Zero Temperature Limit of Holographic Superconductors,” JHEP 0911, 015 (2009) [arXiv:0908.3677 [hepth]].
 [30] K. Goldstein, S. Kachru, S. Prakash and S. P. Trivedi, “Holography of Charged Dilaton Black Holes,” arXiv:0911.3586 [hepth].
 [31] M. Cadoni, G. D’Appollonio and P. Pani, “Phase transitions between ReissnerNordstrom and dilatonic black holes in 4D AdS spacetime,” JHEP 1003, 100 (2010) [arXiv:0912.3520 [hepth]].
 [32] J. Gauntlett, J. Sonner and T. Wiseman, “Quantum Criticality and Holographic Superconductors in Mtheory,” JHEP 1002, 060 (2010) [arXiv:0912.0512 [hepth]].
 [33] V. Oganesyan, S. A. Kivelson and E. Fradkin, “Quantum theory of a nematic Fermi fuid,” Phys. Rev. B 64, 195109 (2001).
 [34] E. Fradkin, S.A. Kivelson, M. J. Lawler, J. P. Eisenstein, and A. P. Mackenzie, “Nematic Fermi Fluids in Condensed Matter Physics”, arXiv:0910.4166
 [35] S. Nakamura, H. Ooguri and C. S. Park, “Gravity Dual of Spatially Modulated Phase,” arXiv:0911.0679 [hepth].