Fermionic probes of local quantum criticality in one dimension
We study the spectral functions of fermionic operators in dimensional Super Yang-Mills theory with 16 supercharges at finite density using the holographically dual D1-brane geometry. This system exhibits quasi-particle peaks indicating the existence of Fermi points about which excitations have a finite relaxation time, in contrast with the Tomonaga-Luttinger model. The finite width may be attributed to the non-trivial interactions of the probe operators with the background density matrix, modeled holographically as a charged black hole. We show that the fermionic correlators can in fact be deduced from known results for fermion probes of the charged AdS black hole background, owing to some remarkable coincidences in supergravity truncations.
Interacting Fermi systems in one spatial dimension (1d) provide the textbook example of departure from conventional Fermi liquid theory. For reasons associated with the restricted phase space in 1d, particle-hole pairs are indistinguishable from collective excitations. Whilst in higher dimensional fermionic systems have low energy fermionic quasiparticles, in 1d we have bosonic particles – spin and charge density waves.
The prototype model for interacting 1d fermions is the Tomonaga-Luttinger (TL) model which involves linearising the single particle dispersion about the Fermi points. This has the advantage of being exactly soluble via bosonisation Haldane (1981) and is a conformal field theory. Retaining nonlinearity in the spectrum has recently shown to be possible and the spectral function for nonlinear Luttinger (NLL) liquids has recently been computed in perturbation theory Khodas et al. (2007); Imambekov and Glazman (2009).
In this paper we are going to explore the dynamics of fermions in a non-conformal field theory, Super Yang-Mills with 16 supercharges in dimensions (SYM) at finite density. We will extract the spectral function for this system in the planar limit by exploiting the fact that it admits a holographic dual Itzhaki et al. (1998) for intermediate energy scales, , , in terms of the near-horizon geometry of a stack of parallel D1-branes in Type IIB supergravity. Here we introduced the ’t Hooft coupling whose engineering dimension is .
As SYM is a non-conformal theory, one anticipates that the spectral function will deviate significantly from the TL form.
While the example we focus on also has a local quantum critical region, we shall in fact find that the effective holographic master field equations that govern the spectral function are identical despite the seemingly different microscopic nature of the field theories. To be precise we find that the spectral functions for fermionic operators in SYM are identical to those obtained in the dimensional field theory of M2-branes. This is made possible by some remarkable coincidences in supergravity as we shall see in due course.
Ii The gravitational set-up
SYM for is described holographically in terms of the supergravity geometry obtained by taking the near-horizon limit of a stack of D1-branes. As we are interested in the dynamics of the theory in the grand canonical ensemble with non-zero temperature and charge density, we consider configurations where the D1-branes are non-extremal and spinning in the transverse . This turns on -charges along the Cartan and for simplicity we focus on the equally-spinning case when these are all equal.
Whilst the dimensional (D) Type IIB solution for this background configuration is easily written down Harmark and Obers (2000), it is convenient to work with a KK-reduced solution in 3d which is obtained by reducing Type IIB equations of motion on a with internal angular momentum David et al. (2010). This Lagrangian comprises of a dilaton , a gauge potential, and the metric, with the effective action
where . The spinning D1-brane is described as a charged, dilatonic black hole solution of (1). Employing a prescient radial parametrisation given by the dilaton, , it is given by
The horizon is located at and the black hole and dual field theory temperature is given by . Alert readers will recognise the close similarity between this solution and the four-dimensional Reissner-Nordstrom AdS solution, which will have an important role to play in the following.
ii.1 Fermionic probes
The field theory SYM has fermionic operators residing in super-multiplets whose lowest component comprises of scalar chiral primaries, , transforming symmetric traceless representations of the R-symmetry group. We recall that the Lagrangian for SYM can be obtained by dimensionally reducing the 10d SYM Lagrangian on a and hence apart from the gluons it contains adjoint scalars which transform as vectors under .
The dimensions of the operators can be computed by examining the two-point function in the vacuum. In the energy regime where the dual supergravity description is valid these take on integral values with . Equivalently, this spectrum can be described in terms of fermionic supergravity fields propagating in effective 3d D1-brane geometry by studying linearised KK fluctuations about the vacuum D1-brane solution. A comprehensive analysis of this reduction was carried out in Morales and Samtleben (2002) and it was found that the masses are related to their conformal dimension via
To compute the spectral function of using supergravity it suffices to consider a quadratic fermion effective action. For a fermionic mode with mass , we use
where , is the spin-connection, and . Underlined indices refer to the local tangent space. Note that this action is part of a larger interacting theory resulting from the reduction of IIB supergravity, where the mass term appears with a dilatonic coupling as a consequence of this reduction.
Given this set-up it is a straightforward to analyze the Dirac equation resulting from (6) in the background of the spinning D1-brane geometry. The fermionic correlators are obtained by a suitable modification of the prescription given in Iqbal and Liu (2009) to account for the fact that we have a geometry that is not asymptotically AdS. However, a simple trick allows us to reduce the computation to a previously solved problem and moreover provides a useful perspective. Therefore we will now take a detour into the mysterious world of supergravity truncations.
Iii Supergravity truncations
It has long been known that the dilatonic branes geometries are conformal to AdS Boonstra et al. (1999). It is for instance clear that (2) is asymptotically conformal to AdS. What is more surprising is that this geometry has a hidden AdS origin. For one, it can immediately be oxidised to a Reissner-Nordstrom black hole in AdS using the ansatz (the gauge field lifts trivially)
where denotes lift direction. This metric, in turn, solves field equations of the Einstein-Maxwell system with negative cosmological constant,
Connections between these two systems have been observed previously Peet and Polchinski (1999); Mateos et al. (2007); Kanitscheider and Skenderis (2009); David et al. (2009, 2010). Here we go further and illuminate the origin of these connections as part of a more general story. The basic picture is presented in Fig. 1.
Starting with D supergravity we may reduce on and perform a consistent truncation to the well known gauged supergravity de Wit and Nicolai (1987). The massless modes one retains are the graviton transforming in the singlet of , gravitino in the , while the fermions transform as and the scalars fill out and .
Likewise for Type IIA/B supergravity in D where only a dilaton and 3-form field strength are retained, there remarkably is also a consistent truncation when reduced on Cvetic et al. (2000) to gauged supergravity in D. The spectrum comprises of the graviton, a gauge field, 36 real scalars, transforming as and a singlet of , gravitinos transforming as and fermions in the .
We have now at hand two different gauged supergravity theories with gauge group. For 4D solutions where no pseudo-scalars are sourced (the ), the resulting effective bosonic theories are related through an reduction at the massless, nonlinear level. Whilst this is not a consistent truncation, it does allow us to relate a large class of solutions without pseudo-scalars, including the situation of current interest.
Given this observation it follows abstractly that the Dirac equation arising from (6) maps directly into the fermion wave equation on Reissner-Nordstrom AdS background studied in Liu et al. (2009); Faulkner et al. (2009). The map arises from simply KK-reducing on a spatial direction of this AdS, which amounts to restricting to the zero-momentum sector in the -direction.
Iv Spectral functions
iv.1 Numerical results at , ,
Fermion Green’s functions in theories with asymptotically AdS dual geometries are readily obtained using standard holographic techniques Iqbal and Liu (2009). In Appendix A we present the method adapted to the D1 geometry (2).
iv.2 Characterising the Fermi points
Semi-analytic expressions for the correlators can be obtained for small frequencies, , using a matched asymptotic expansion Faulkner et al. (2009). For the M2-brane geometry this relies on the fact that the geometry (7) interpolates between a near horizon and an asymptotically geometry. These results may be dimensionally reduced to yield expressions for Green’s functions of the D1 system. Indeed, the analogous matching calculation for the D1 system yields identical results, as its quantum critical behaviour in the IR is determined by an geometry consistent with the semi-holographic picture of Faulkner and Polchinski (2010).
The location of each peak is determined by the existence of a normalisable, in-falling fermion mode at . These peaks are indicated by green crosses in the example of Fig. 2. In the vicinity of each such peak we obtain,
where indicates the momentum of the peak in question, . , and are real and can be determined numerically. is the retarded Green’s function computed in the near horizon :
iv.3 Peak broadening
The type of quasi-particle behaviour depends on and there is no stable quasi-particle associated with peaks satisfying Faulkner et al. (2011). When , the magnitude of the AdS contribution (10) is subdominant to and in the vicinity of the pole at we obtain,
We have retained the leading imaginary contribution in the denominator, , which gives a width to the peak away from and a finite quasi-particle relaxation time.
The Green’s function (11) is the same as a slice of the M2 theory, in particular, the quasi-particle width is the same. In contrast, the restricted phase space in the TL model leads to zero-width quasi-particle excitations. However, the 1d NLL models – which modify TL by an correction to the single-particle dispersion relation and density-density interaction – exhibit Green’s functions of the form (11) near the pole, with width Khodas et al. (2007). The TL and NLL models also displays non-trivial power-law behaviour in the real-part of the Green’s function, while (11) resembles a free-fermion-like linear dispersion near . Thus the NLL models demonstrate that the differences between physics of Fermi points and higher dimensional Fermi surfaces need not be drastic, as indeed we find here.
Of course, neither the TL nor the NLL models approximate the D1 system.
We studied a holographic model of a fermionic quantum liquid with a concrete dual QFT: SYM with 16 supercharges. Thus our holographic results provide in-roads to a system which may be tackled in other ways, such as lattice field theory. In the interest of phenomenology we focused on generic probes of this system; we leave supergravity specific probes for future study. It would also be interesting to study condensation of the adjoint scalars.
Where possible we confirm exact agreement between the direct holographic calculations in the D1 geometry and dimensional reduction of existing results with AdS asymptotics. The spectral functions for SYM indicate the existence of quasi-particle peaks indicative of a Fermi surface, whose shape is controlled by the IR geometry and exhibit a range of behaviours depending on .
The D1 spectral functions presented here differ from those of the TL model. The spectral support has a finite width about the quasi-particle peak (see e.g., Fig. 2) indicative of nonlinear effects broadening the Fermi peak and is further consistent with lack of conformal invariance of SYM.
Of course, the TL model is idealised since the single particle spectrum is linearised about the Fermi points. Incorporating non-linearity in the dispersion enables the redistribution of quasi-particle momenta, leading to finite quasi-particle lifetimes Khodas et al. (2007); Imambekov and Glazman (2009). In the holographic model, rather than having explicit modification of the fermionic spectrum, we find that the effects of the charged black hole background provide a non-linearisation of the spectrum.
A curious feature of our analysis is the connection we have drawn between the SYM and the M2-brane theories. Whether this is an artefact of supergravity or holds deeper significance remains to be better understood. For the current discussion we note that it is rather amusing to realise that the analysis of fermionic correllators in the M2-brane theories can be related to a model in 1d, reminiscent of old ideas attributing strange metal phase of high superconductors to Luttinger liquid-like dynamics.
Acknowledgements.We thank J. Gauntlett, S. Gentle, V. Hubeny, H. Liu, J. McGreevy and S. Minwalla for helpful discussions. MR is supported in part by an STFC rolling grant. BW is supported by a Royal Commission for the Exhibition of 1851 Research Fellowship.
Appendix A Green’s functions in 3d
While we have emphasised the connection between the 3d and 4d effective actions in the main text, retarded Green’s functions for may be obtained without invoking this using standard holographic techniques. We introduce and pick a plane wave basis, . We further define projections under as We supplement (6) with the boundary term, where is the induced metric on the boundary. The variation of the total fermion action , when on-shell is a boundary integral by virtue of the equations of motion:
Counterterms are needed if , and for there is a logarithm in the boundary expansion. For convenience we will focus here on . Retarded Green’s functions are given by imposing in-falling boundary conditions at the horizon, and reading off data at the boundary, where The integration in the bulk is most convenient for the ratio (dropping subscripts), since its equation of motion resulting from (6) is first order,
where . Hence the Green’s function is easily obtained by integrating (13) from to the boundary where .
- preprint: DCPT-11/27
- Luttinger liquids have been modelled holographically in terms of fermionic fields probing a BTZ black hole Balasubramanian et al. (2010).
- Recently fermionic correlators were examined in Einstein-Maxwell-dilaton gravity in Iizuka et al. (2011) albeit without the dilaton coupling in the mass term.
- In addition one also projects the four-component spinor equation at into a pair of independent equations for two-component spinors. One equation matches precisely that which arises from (6), whilst the other matches with .
- It is however amusing to note that excitations about a Fermi point in the D1 system resemble the qualitative behaviour of relaxation via three particle scattering, the decay mode responsible for .
- F. D. M. Haldane, J. Phys. C 14, 2585 (1981).
- M. Khodas, M. Pustilnik, A. Kamenev, and L. Glazman, Physical Review B 76, 155402 (2007).
- A. Imambekov and L. Glazman, Science 323, 228 (2009), ISSN 0036-8075.
- N. Itzhaki, J. M. Maldacena, J. Sonnenschein, and S. Yankielowicz, Phys. Rev. D58, 046004 (1998), eprint hep-th/9802042.
- V. Balasubramanian, I. Garcia-Etxebarria, F. Larsen, and J. Simon (2010), eprint 1012.4363.
- H. Liu, J. McGreevy, and D. Vegh (2009), eprint 0903.2477.
- T. Faulkner, H. Liu, J. McGreevy, and D. Vegh (2009), eprint 0907.2694.
- T. Faulkner, N. Iqbal, H. Liu, J. McGreevy, and D. Vegh (2011), eprint 1101.0597.
- S.-S. Lee, Phys.Rev. D79, 086006 (2009), eprint 0809.3402.
- M. Cubrovic, J. Zaanen, and K. Schalm, Science 325, 439 (2009), eprint 0904.1993.
- T. Faulkner and J. Polchinski (2010), eprint 1001.5049.
- T. Harmark and N. A. Obers, JHEP 01, 008 (2000), eprint hep-th/9910036.
- J. R. David, M. Mahato, S. Thakur, and S. R. Wadia (2010), eprint 1008.4350.
- J. F. Morales and H. Samtleben, JHEP 08, 042 (2002), eprint hep-th/0206247.
- N. Iizuka, N. Kundu, P. Narayan, and S. P. Trivedi (2011), eprint 1105.1162.
- N. Iqbal and H. Liu, Fortsch. Phys. 57, 367 (2009), eprint 0903.2596.
- H. Boonstra, K. Skenderis, and P. Townsend, JHEP 9901, 003 (1999), eprint hep-th/9807137.
- A. W. Peet and J. Polchinski, Phys. Rev. D59, 065011 (1999), eprint hep-th/9809022.
- D. Mateos, R. C. Myers, and R. M. Thomson, JHEP 05, 067 (2007), eprint hep-th/0701132.
- I. Kanitscheider and K. Skenderis (2009), eprint 0901.1487.
- J. R. David, M. Mahato, and S. R. Wadia, JHEP 04, 042 (2009), eprint 0901.2013.
- B. de Wit and H. Nicolai, Nucl.Phys. B281, 211 (1987).
- M. Cvetic, H. Lu, and C. N. Pope, Phys. Rev. D62, 064028 (2000), eprint hep-th/0003286.