A Quantum Critical Point from Flavours on a Compact Space.
Abstract
We analyse a dimensional defect field theory on a two sphere in an external magnetic field. The theory is holographically dual to probe D5branes in global AdS background. At any finite magnetic field only the confined phase of the theory is realised. There is a first order quantum phase transition, within the confined phase of theory, ending on a quantum critical point of a second order phase transition. We analyse the condensate and magnetisation of theory and construct its phase diagram. We study the critical exponents near the quantum critical point and find that the second derivatives of the free energy, with respect to the bare mass and the magnetic field, diverge with a critical exponent of . Next, we analyse the meson spectrum of the theory and identify a massless mode at the critical point signalling a diverging correlation length of the quantum fluctuations. We find that the derivative of the meson mass with respect to the bare mass also diverges with a critical exponent of . Finally, our studies of the magnetisation uncover a persistent diamagnetic response similar to that in mesoscopic systems, such as quantum dots and nano tubes.
Keywords:
DIASSTP1408
^{†}^{†}institutetext: School of Theoretical Physics, Dublin Institute for Advanced Studies
10 Burlington Road, Dublin 4, Ireland.
\subheader
1 Introduction
Quantum phase transitions are phase transitions which take place at zero temperature and are thus driven by quantum as opposed to thermal fluctuations. Such phase transitions are realised by tuning a nontemperature control parameter such as pressure, chemical potential or magnetic field. Of particular interest are second order quantum phase transitions which, similarly to their thermal analogue, have a quantum critical point where the quantum fluctuations driving the transition diverge and become scale invariant. The properties of these systems near the phase transition are characterised by critical exponents, which are independent on the microscopic details of the systems – a property called Universality. At the critical point fluctuations of all size scales become important demanding a scale invariant description, naturally provided by conformal field theories. Universality can then be seen as a byproduct of the fact that there are relatively few scale invariant theories.
Although quantum phase transitions take place at zero temperature, at finite temperature quantum fluctuations compete with thermal fluctuations and the existence of a “quantum critical” regime is expected, bearing remnants of the quantum critical point even at finite temperature. Experimentally the quantum critical regime becomes manifest in an exotic physical behaviour often signalling the existence of a novel phase.
Quantum critical points are an active area of research in condensed matter community, where they appear in the phase diagram of high temperature superconductors Sachdev (). Clearly the existence of a quantum critical regime accessible at finite temperature is crucial for such applications, however away from the critical point one needs more than just an effective field theory and the analysis often relies on physical arguments and approximations rather than rigorous calculations. It is therefore, of a particular importance to come up with clean theoretical models, which not only allow full analytic (or quasi numerical) control, but are complex enough to release the vicinity of a quantum critical point. The holographic quantum critical point that we propose in this paper is an example of such a model and is potentially applicable to the qualitative description of the quantum critical regions in realistic condensed matter systems such as high temperature superconductors.
The original AdS/CFT correspondence establishes a duality between an super YangMills theory (SYM) in dimensions, which is a conformal field theory, and superstring theory on an AdS spacetime. There is no known quantum critical point described by SYM. In fact, although a vast amount of novel supergravity backgrounds dual to conformal field theories have been constructed, usually their construction is not associated with the critical regime of a particular theory. As a result relatively few holographic quantum critical points have been realised, usually in the context of holographic flavour dynamics.
A system with continuous quantum phase transition has been constructed in ref. Karch:2009ph (), however this construction involved a third order phase transition. A second order quantum phase transition has been realised in ref. Karch:2007br (), where the holographic gauge theory dual to the D3/D7 intersection at zero temperature and finite chemical potential has been explored. A quantum critical point has also been constructed in refs Jensen:2010vd (); Evans:2010iy (), where the holographic gauge theory dual to the D3/D5 intersection at finite magnetic field and chemical potential has been analysed. In this paper we porpoise a novel quantum critical point realised in a 1+2 dimensional defect field theory living on a maximal two sphere and subjected to an external magnetic field.
Our holographic set up involves probe D5branes in a global AdS spacetime. The D5branes are extended along the radial direction and wrap a maximal two sphere inside the AdS part of the geometry. The D5branes also wrap a two sphere inside the part of the geometry. The corresponding gauge theory is an SYM theory on a three sphere coupled to an hypermultiplet living on a maximal two sphere. The system is supersymmetric only at vanishing bare mass, when the D5brane wraps a maximal two sphere inside the . At finite bare mass supersymmetry is broken and the theory develops a fundamental condensate.
The properties of the theory at zero magnetic field can be further analysed by considering a pair of external probe quarks. A key observation is that at strong coupling the screening length is inversely proportional to the bare mass Karch:2009ph (). This suggests that for large bare masses the screening length is small compared to the radius of the two sphere and the theory is in a confined phase (separating the quarks would result in a pair production). For small bare masses the screening length is larger than the radius of the sphere and a pair of quarks can be separated without pair production – the theory is in a deconfined phase. The quantum phase transition between the two regimes is of a first order Karch:2009ph () and is triggered by the Casimir energy of the theory. In the holographic set up the phase transition is realised as a topology change transition for D5brane embedding. The confined phase of the theory corresponds to embeddings, which close above the origin of AdS by having the inside the part of the geometry shrink to zero size, referred to as “Minkowski” embeddings, while the deconfined phase of the theory corresponds to embeddings which reach all the way to the origin of AdS and have the inside the AdS part of the geometry shrinking, we will refer to them as the “Ball” embeddings.
At first glance subjecting the theory to an external magnetic field seems like an innocent exercise. Indeed, the second betti number of the two sphere is nonzero and turning on a magnetic field on the two sphere does not require additional electric current. In the holographic set up this can be achieved by either introducing a pure gauge field or by fixing the gauge field of the D5brane. However, although such solutions can be constructed relatively easy, for ball embeddings one encounters a problem. At the origin of AdS the radius of the two sphere shrinks to zero size, but the magnetic flux throughout the two sphere remains finite, as a result the norm of the field strength diverges at the origin, signalling the presence of a magnetic monopole. One can show that a magnetic monopole on the world volume of the D5brane would violate charge conservation of the RRflux unless some number of D3–branes are attached at the position of the monopole. On the other hand, doing so would deform the D5brane embedding due to the tension of the D3–branes.
In general constructing a balanced D3D5brane configuration can be a challenging task, but in our case a small miracle happens: the “Ball” D5brane embeddings are mimicked by Minkowski embeddings, which close at very small but finite distance above the origin of AdS thus avoiding divergent norm of the field strength and the need of a magnetic monopole. Furthermore, these embeddings fold along the part of the geometry, while still close to the origin, developing a throat with very small radius of the , mimicking the D3–branes needed to source the magnetic monopole in the unbalanced D3D5 configuration (see figure 3).
One can show that Minkowski embeddings cover the whole parametric space, completely replacing the “‘Ball” embeddings, this picture persist even for an infinitesimal magnetic field. The first order confinement/deconfiment phase transition also has its analogue within the confined phase described by the Minkowski embeddings.^{1}^{1}1This situation is the magnetic analogue of the finite density study of ref. Kobayashi:2006sb (), where the radial electric field of the probe D7brane demands the presence of an electric monopole for Minkowski embeddings. The fundamental strings required to source the monopole are realised as spikes of black hole embeddings mimicking Minkowski embeddings. However, unlike the phase transition at vanishing magnetic field, there is no associated topology change for the D5brane embeddings (only an approximate one between embeddings with and without an throat) as a result we find that at sufficiently strong magnetic field the first order phase transition ends on a critical point of a second order phase transition. Increasing further the magnetic field leads to the replacement of the phase transition by a smooth cross over. Focusing on the properties of the theory near the quantum critical point we find that all second derivatives of the free energy diverge as functions of the bare mass and the external magnetic field with critical exponent . In addition, our studies of the meson spectrum show that at the critical point there is a massless mode which we associate to the divergent correlation length of the quantum fluctuations. The derivative of the mass of this mode with respect to the bare mass also diverges at the phase transition with a critical exponent .
Another interesting property of the theory is its magnetic response. Our studies of the magnetisation and magnetic susceptibility show that the theory is a diamagnetic, which is not unusual for holographic gauge theories (see for example ref. Filev:2013vka ()). However, our system has peculiar properties in the limit of vanishing magnetic field. The analogue of the deconfined phase has a persistent diamagnetic response independent on the magnitude of the magnetic field. This behaviour is similar to the persistent diamagnetic current observed in realistic mesoscopic systems such as nano tubes and quantum dots Zelyak (). The fact that this phase is realised for small radius of the two sphere, when the size of the system is small and finite size quantum effects (such as the existence of Casimir energy) are important, makes the comparison to mesoscopic systems plausible. Furthermore, the fact that the observed effect is in the strongly coupled regime of the system suggests potential applications in the description of strongly coupled condensed matter systems, exhibiting persistent diamagnetic response.
Finally, for any finite value of the magnetic field the sable phase of the theory at vanishing bare mass has a negative condensate. This suggests that the global symmetry corresponding to rotations in the transverse space to the D5branes is spontaneously broken. In this way our system realises the effect of magnetic catalysis of chiral symmetry breaking confirming the universal nature of this phenomena magcat (). ^{2}^{2}2For recent holographic studies of this phenomena in dimensional systems see refs. Evans:2013jma (); Kristjansen:2013hma (); Filev:2014bna ()
The paper is organised as follows: In section 2 we review the properties of the theory at zero magnetic field studied in refs. Karch:2009ph (), Erdmenger:2010zm (). In section 3 we explore the properties of the theory at finite magnetic field. We show that at any finite magnetic field only Minkowski embeddings are physical and only the confined phase of theory is realised. We study the fundamental condensate of the theory as a function of the bare mass and show that for sufficiently small magnetic fields there is a first order phase transition within the confined phase, which ends on a critical point of a second order phase transition. We also calculate the magnetisation of the theory and its magnetic susceptibility. Finally, we study the critical behaviour of the theory near the quantum critical point and calculate various critical exponents. In section 4 we analyse the meson spectrum of the theory. We obtain the general quadratic action for all four scalar and six vector meson modes. We then focus on the spectrum of one of the scalar modes, which becomes massless at the quantum critical point, corresponding to the divergent correlation length of the quantum fluctuations. Finally, we focus on the ground state of the spectrum near the phase transition and calculate the corresponding critical exponent. Section 5 contains our conclusion.
2 The theory at zero magnetic field
The supergravity dual of SYM theory on a three sphere is an AdS spacetime in global coordinates with a metric given by:
(1) 
To introduce fundamental flavours in the quenched approximation to the dual field theory we will consider a probe D5brane extended along the radial coordinate and wrapping two spheres in both the and the AdS parts of the geometry. More explicitly if we parameterise the AdS and one of the by:
(2)  
then the D5brane probe is extended along the and directions. Note that the D5–brane respects the symmetry of the dial defect field theory by wrapping completely the boundary in the AdS part of the geometry. One can check that placing the D5brane at and allowing the D5brane to have a nontrivial profile along is consistent with the equations of motion. The most symmetric ansatz (which is also consistent with the general equations of motion) is to consider . The corresponding field theory is coupled to an fundamental hypermultiplet living on the two sphere. The theory is at most supersymmetric, and one can show that a nontrivial profile completely breaks supersymmetry. When supersymmetry is broken the theory develops a fundamental condensate, with nontrivial dependence on the bare mass of the theory. Furthermore, the theory has a first order phase transition corresponding to a topology change transition in the holographic set up. Indeed, generally the possible embeddings of the D5brane split into two classes  embeddings which close at some radial distance above the origin of AdS by having the inside the part of the geometry shrink to zero size (Minkowski embeddings), and embeddings which reach all the way to the origin of AdS and close there, which we call “Ball” embeddings. This construction has been investigated in ref. Karch:2009ph (), where it has been shown that the phase transition is of a first order. It has also been argued that this is a confinement/deconfinement phase transition driven by the Casimir energy of the theory on . The properties of this theory have been further investigated in ref. Erdmenger:2010zm (), where the meson spectrum of the theory has been studied in details.
In what follows, we will briefly review the properties of the theory at zero temperature and external magnetic field. For illustrative purposes it is instructive to introduce a new radial coordinate given by:^{3}^{3}3Note that in these coordinates the origin of AdS is at .
(3) 
which enables us to write part of the geometry as a conformal :
(4) 
The conformally part of the geometry can be split into a conformally part wrapped by the D5brane and a conformally part transverse to the D5brane:
(5) 
where and . The possible D5brane embeddings can then be visualised in the versus plane. While we will eventually present our results in these coordinates, technically it is more convenient to work in the original (,) coordinates. For the radial part of the corresponding DBI action (after integrating over the and the AdS) one obtains:
(6) 
The solution to the equation of motion for has the following asymptotic behaviour at large :
(7) 
and according to the AdS/CFT dictionary, the parameter is related to the bare mass of the fundamental multiplet via , while the parameter is proportional to the fundamental condensate of the theory . Note that to arrive at the last result one has to introduce appropriate counter terms Karch:2009ph () (see also Karch:2005ms ()). However, with this definition one obtains a constant nonzero value of the condensate at infinite bare mass, suggesting that an additional counter term is needed.
In general the equation of motion for derived from (6) has to be solved numerically. One can check that for both Minkowski and “Ball” embeddings, the equation of motion has predetermined boundary conditions at the closing point, in the sense that the solution depends only on one parameter. This enables us to obtain a perturbative analytic solution, which we feed into a numerical shooting technique. A set of possible D5brane embeddings is presented in figure 1. The red curves represent Ball embeddings reaching the origin (the black segment), while the blue curves represent Minkowski embeddings. One can see that the two classes are separated by a critical embedding (the dashed black curve), which has both the and the AdS shrinking at the origin.
By extracting the coefficients and from the asymptotic behaviour of the embeddings at large radial distances , one can generate the plot of the condensate versus bare mass presented in figure 2, where we have introduced the dimensionless parameters and . Note, that the physical meaning of the parameter follows from the relation , where is the radius of the defect (which is equal to radius of the field theory ) and is the t’Hooft coupling. At large the radius of the two sphere is large (at fixed bare mass), the Casimir energy is small and the theory is in a confined phase (see refs. Karch:2009ph (),Karch:2006bv ()). As decreases the Casimir energy increases and at there is a first order confinement/deconfinement phase transition triggered by the Casimir energy. One can also see from the plot that at large the condensate approaches a constant value. In fact one can show analytically that the parameter approaches . Since large corresponds to large bare masses this is a surprising result. One would expect that as the flavours become infinitely massive they will decouple from the theory and hence the fundamental condensate will vanish. This suggests that a finite counter term (that was not prescribed in refs. Karch:2009ph (),Karch:2005ms ()) is needed.
Given that the condensate approaches a constant value at large bare mass, and hence the free energy grows linearly with the bare mass, the appropriate counter term should be linear in the supergravity field . In fact, in order to preserve the symmetry of the Lagrangian (6), it should be proportional to . Indeed, one can show that a counter term proportional to , where is the induced metric on the asymptotic boundary spanned by the field theory directions of the defect, and is the corresponding scalar curvature, produces a term linear in ^{4}^{4}4Note that naively this choice of a counter term breaks the scaling symmetry of the boundary field theory, However, the flavoured gauge theory is conformal only for vanishing bare mass, when and the counter term vanishes. One can show that for both conformal symmetry and supersymmetry are explicitly broken and hence the proposed counter term does not break any further symmetries. However, this counter term induces a discontinuity of the condensate of the theory because . Since we are mainly interested in the properties of the theory near the phase transition, we will take the same approach as in ref. Karch:2009ph () and use the present definition of the condensate (being proportional to ).
3 The theory at finite magnetic field
3.1 General set up
We are interested in the properties of the theory at finite magnetic field, which (as advertised in the introduction) realises a quantum critical point. Because the second betti number of the two sphere is nonvanishing a magnetic field can be turned on without the introduction of external electric currents. A natural choice in this case is to consider a magnetic field proportional to the volume form of the twosphere. In order to turn on such a magnetic field, we turn on a pure gauge Bfield:
(8) 
One can check that for . Thus equivalently, we could introduce magnetic field by turning on the component of the U(1) gauge field. In both cases one should check^{5}^{5}5Note that as pointed out in ref. Chunlen:2014zpa () the analogous ansatz in the case of D7brane probe considered in ref. Filev:2012ch () does not satisfy the equations of motion for the gauge field. The physical reason for that is that the second betti number of the three sphere is zero and one has to introduce external currents to support such a magnetic field. The study considered in ref.Filev:2012ch () has been reinterpreted as a bottom up study, where one introduces an external flux (serving as external currents) to support the magnetic field on the D7brane. that the equations of motion for the gauge field are satisfied. The DBI action of the D5brane is then given by:
(9) 
defining for the equation of motion for the gauge field we obtain:
(10) 
One can easily check that for the ansatz , and given by equation (8), the equation of motion for the gauge field (10) is satisfied. The equation of motion for is derived by varying the DBI action (9):
(11) 
Naively, one expects that for moderate magnetic field the behaviour of the system should be similar to the zero magnetic field case. One expects that the embeddings would again split into “Ball” and “Minkowski” classes and there would be a first order phase transition between the two, at some critical value of the parameter . Furthermore, because of the universal nature of the effect of magnetic catalysis, one would expect that at strong magnetic fields the phase transition disappears and only the confined phase, with chiral symmetry breaking vacuum at vanishing bare mass, is stable. However, as we are going to see, something much more intricate happens and only part of our naive expectations are met.
3.2 Branes and magnetic monopoles
Let us calculate the norm of the Bfield , we obtain: . Therefore, for “Ball” embeddings the norm of the Bfield diverges at the origin of AdS. This simply reflects the fact that for our ansatz the magnetic flux through the of the defect remains constant, and since at the origin of AdS the twosphere shrinks to zero size, there should be a magnetic monopole^{6}^{6}6Note that by magnetic monopole we mean a charged object Sdual to an electric monopole. sitting there. The world volume of the D5brane is six dimensional and in six dimensions a magnetic monopole is a two dimensional object. In our case the monopole should be wrapping the , which remains finite at the origin. It has been known for quite a while that a Dp2 brane ending on a Dpbrane sources a magnetic monopole in the world volume of the Dpbrane Strominger:1995ac (). We can easily see this. The WessZummino term of the D5brane contains a coupling of the form:
(12) 
where is a four form RamondRamond potential. Now gauge invariance of the WessZummino action under the transformation suggests that the term:
(13) 
should vanish, where we have used that and is spanned by the time and directions at the origin of AdS. Therefore, we see that the magnetic monopole induces a three dimensional “boundary” of the D5brane and to preserve gauge invariance we have to attach some number () of D3branes to the D5brane. The variation of the WessZummino action of the D3branes with respect to the gauge transformation of would result in a term , which can be used to cancel the term in equation (13). For the number of D3branes required to source the magnetic monopole we obtain:
(14) 
To visualise the D3branes, that we need to attach to the D5brane, let us plot a typical “Ball” D5brane embedding. The resulting plot (for ) is presented in figure 3. The solid red curve represents a “Ball” D5brane embedding ending at the origin at . The green dashed segment represent a D3brane sitting near the origin of AdS, wrapping an and extended along and , The D3brane ends on the D5brane at , and closes smoothly at thus having a topology . Of course, this configuration is not stable. The tension of the D3brane will pull the D5brane and the configuration will change. Remarkably, the stable configuration is found among the Minkowski class of embeddings! The blue curve represents a Minkowski embedding mimicking the D3D5 configuration. The D3brane is realised by the Minkowski embedding as a throat of topology , where the AdS has a very small radius along the segment near the origin. Note that because Minkowski embeddings close above the origin, the AdS never shrinks to zero size, avoiding the need of a localised magnetic monopole.^{7}^{7}7Note that this situation is a magnetic analogue of the finite density study of ref. Kobayashi:2006sb (), where the radial electric field of the probe D7brane requires the presence of an electric monopole for Mikowski embedding, the fundamental strings required to source the monopole are realised as spikes of black hole embeddings mimicking Minkowski embeddings.
In fact, we can present further evidence that the throat of the D5brane realises D3branes. Using that it is located in a region , for the part of the DBI action corresponding to the throat (to leading order) we obtain:
(15) 
where:
(16) 
and we have used equation (14). Note that the integrant in the last equation in (15) is the volume form associated to the induced metric on a D3brane wrapping the and extended along , therefore the total expression is precisely the DBI action of a stack of D3branes wrapping the and extended along and for . In the next subsection we will use this set up to study the fundamental condensate of the theory.
3.3 Fundamental condensate
Our numerical studies confirm that Minkowski embeddings cover the entire parametric space and are the thermodynamically stable phase for any arbitrary small nonzero magnetic field. What is more, we find that the first order confinement/deconfinement phase transition, at vanishing magnetic field, has an analogue and there is a first order phase transition, within the confined phase, ending on a critical point of a second order quantum phase transition.
To obtain our numerical results we first introduce the dimensionless radial coordinate and study numerically the function , in these coordinates the Lagrangian corresponding to (11) becomes:
(17) 
As in the previous section we will present our results in terms of the radial coordinate (defined in equation (3)), whose dimensionless analogue is given by . The presence of external magnetic field does not introduce new counter terms and using the same definition for the condensate as in the previous section, we read the bare mass and the condensate of the theory from the asymptotic expansion of the field at large :
(18) 
where and is proportional to the fundamental condensate of the theory . We also use that the nonlinear second order differential equation for derived from (17) has predetermined boundary conditions at () and depends only on one parameter – the radial distance above the origin of AdS where the D5brane closes. Solving perturbatively for near and feeding this solution into a numerical shooting technique, we can generate all possible D5brane embeddings.
A plot of the condensate as a function of the bare mass for various values of the magnetic field has been presented in figure 4. As one can see, there is a multivalued loop centred at the origin of the plane. This loops is similar to the spiral structure, which generally appear in holographic gauge theories dual to the “flat” Dp/Dq intersection, when subjected to an external magnetic field Filev:2009xp (). However, unlike the scenario described in refs. Filev:2009xp (),Filev:2007qu () in the present case the step of the spiral is extremely small and it seems that the spiral approaches a limiting loop (rather than shrinking to zero radius). In practice, in the plots we cannot distinguish between the spiral and the limiting loop and the curve of the condensate winds along the loop. The origin of this behaviour in the holographic set up can be understood after one parametrises the D5brane embeddings by and verifies that is a fixed point of the equation of motion for . This causes embeddings, which start at infinitesimal radial distances () to bend and fold along the near the origin of AdS, with the function changing rapidly form to , before they enter a steady regime of slowly varying at some still small but larger . In the steady regime the asymptotic properties of the D5branes depend very weakly on changes of the initial (IR) conditions and as a result the “folding” solutions become extremely close to some embeddings which started at larger and give rise to practically the same values of the bare mass condensate at infinity. A particular case of such behaviour is plotted in figure 5. To summarise: the condensate is a singe valued function of the bare mass at large bare masses (large initial parameter ). At very small initial parameters the D5brane embeddings are affected by the fixed point and start folding along the part of the geometry near the origin of AdS, as a result the condensate is multivalued function of the bare mass and as the initial parameter approaches zero the condensate curve winds clockwise along a very tight spiral approaching a limiting loop.
Similarly to the spiral structures of refs. Filev:2009xp (),Filev:2007qu () the “folding” solutions have higher free energy (they are also unstable) and only the first (lower) branch of the spiral remains. Furthermore, since negative bare masses can be mapped to positive bare masses with flipped condensate, one can also show that the (lower) positive branch of the condensate curve has the lowest free energy this is why we focus on its properties.
Finally we point out that at zero bare mass the condensate has a finite negative value (for the stable branch) suggesting that a spontaneous breaking of the global symmetry corresponding to rotations along the transverse two sphere is taking place. Therefore our system realises the effect of magnetic catalysis of spontaneous symmetry breaking magcat ().
3.4 Phase Diagram
The most interesting property of the condensate curve is the existence of a first order phase transition pattern for sufficiently small magnetic fields, which ends at a critical point of a second order quantum phase transition for . In figure 6 we have presented plots of the condensate zoomed in near the region of the phase transition. One can see that for for there is a first order phase transition within the confined phase, which mimics the confinement/deconfinement phase transition at vanishing magnetic field. The critical parameter , at which the phase transition takes place, can be found using Maxwell’s equal are law, or alternatively by computing the density of the regularised wick rotated onshell action, which is identified with the free energy density of the system. The two approaches give the same result, because the condensate is a derivative of the free energy density with respect to the bare mass (at fixed radius of the two sphere and fixed magnetic field). One can also see that at the phase transition ends on a critical point of a second order phase transition and for the phase transition is replaced by a smooth cross over.
The observations from figure 6 can be summarised in the phase diagram presented in figure 7, where we have presented the bare mass parameter in physical notations . The blue curve in the figure represents a critical line of a first order phase transition. The critical line ends on a critical point of a second order phase transition (the solid dot). In the next subsection we focus on the properties of the theory at this critical point.
3.5 Magnetisation
In this subsection we study the magnetisation of the system. Identifying the free energy (spacial) density with the ( dimensional) density of the regularised euclidean “onshell” action, for the free energy density we obtain:
(19)  
where , we used that the volume of the field theory is and we subtracted the counter terms prescribed in refs. Karch:2009ph (), Karch:2005ms (). Using equation (19) and that , one easily arrives at the following expression for the magnetisation:
(20) 
where we have introduced the dimensionless quantity . In figure 8 we have presented a plot of the parameter at fixed magnetic field as a function of the bare mass parameter . The plots represent decreasing values of the magnetic field from to . Note that the magnetisation is negative (for positive ) and hence the system is diamagnetic. The innermost (blue) curve corresponds to and is in practise indistinguishable in the plot from the curve corresponding to the limit . The most interesting property of this curve is that for the magnetisation has a nonzero value^{8}^{8}8In the holographic set up the source of the finite magnetisation in the limit of vanishing magnetic field is the D3brane throat whose tension is proportional to (see (15)). in the limits and (the lower and upper branch of the blue curve). Therefore, in this regime the theory has a persistent diamagnetic response, which is independent of the magnetic field. This effect is very similar to the effect of persistent diamagnetic current in mesoscopic systems, such as nano tubes and quantum dots Zelyak (). Although such systems contain large namer of microscopic elements (electrons, molecules, etc.) so that macroscopic quantities are well defined, they are still relatively small and experience strong finite size quantum effects, which are the main cause of their peculiar properties. Recall that the bare mass parameter is proportional to the radius of the field theory two sphere and hence to the size of the system. Therefore, the persistent diamagnetic response is realised when the systems is sufficiently small () and finite size quantum effects (such as the existence of Casimir energy) are important, this makes the analogy to mesoscopic systems even stronger. As can be seen from the second plot in figure 8 the phase transition from ordinary to “persistent” diamagnetism is of a first order, where the critical parameter has the same value as for the confinement/deconfinement phase transition at strictly vanishing magnetic field reviewed in section 2.
Another interesting property of the magnetisation is its behaviour near the quantum critical point. Plots of the magnetisation and the corresponding magnetic susceptibility as functions of the magnetic field and for are presented in figure 9. From the second plot one can see that the magnetic susceptibility diverges at the phase transition, as expected for a second order phase transition. The question for the critical exponents that govern this behaviour will be addressed in the next subsection where we focus on the critical regime of the theory.
3.6 Critical behaviour
In this subsection we explore the properties of the theory at the critical point, in particular we calculate the critical exponents of various observables. For all observables that we study we find with a very high numerical accuracy that the critical exponent is one third, and for the corresponding diverging susceptibilities we find critical exponent . Generally, it is expected that at a second order phase transition the second derivatives of the free energy will diverge. In our considerations we will keep the radius of the two sphere fixed. The variation of the free energy density is then given by:
(21) 
where is the magnetisation of the theory and is the physical magnetic field.

We begin by studying the second derivative of the free energy with respect to the bare mass at fixed magnetic field . At a fixed radius of the two sphere we have: and , and it is sufficient to study the behaviour of the derivative at a fixed magnetic field (). We will show that as the bare mass parameter approaches the critical value the first derivative of the condensate diverges as , where is the critical exponent characterising the behaviour of the condensate near the phase transition, that is near the phase transition we have:
(22) Note that in general the crucial exponents on both sides of the transition need not agree, however this is often the case in critical phenomena and our system is not an exception. To determine numerically the value of the critical exponent we study the logarithmic derivative whose value in the limit coincides with . Refining our numerical techniques and focusing on the region near the phase transition, we obtain the data presented in figure 10.
Note that the points in a very close proximity to the phase transition have been removed to avoid the numerical evaluation of a term behaving like , which would inevitably generate a huge numerical error. Note also that we have evaluated . A nonlinear fit (the solid curve) to the data, provided a value for the critical exponent of with 95% confidence. We consider this result as a very strong indication that the critical exponent is and the second derivative of the condensate diverges at the phase transition as:
(23) 
Our next focus is the magnetic susceptibility related to the second derivative of the free energy with respect to the magnetic field: .
Figure 11: Plots of the quantities versus at , the dashed curves represent fits. For numerical studies it is convenient to introduce the dimensionless quantity , where is defined in equation (20). Plots of the parameters and have already been presented in figure 9. In this subsection we are interested in the behaviour of the parameters across the phase transition. This is why we generate a plot of the magnetisation at for a small range of the magnetic field around the critical value: . In figure LABEL:fig:fig11 we have presented plots of the quantities versus , the dashed curves represent fits. One can see the excellent agreement on both sides of the phase transition, suggesting that the magnetisation approaches its critical value as and hence for the critical behaviour of the susceptibility we obtain:
(24) We see that the magnetic susceptibility has the same critical exponent as the second derivative of the free energy with respect to the bare mass .

Finally, we are interested in the properties of the mixed derivative across the phase transition. To study its behaviour we generate a plot of the magnetisation parameter at and for a small range of the bare mass parameter around the critical one: . As can be seen from the plots in figure 12, there is an excellent agreement with the fits, suggesting that the second derivative diverges as:
(25) across the phase transition. We will return to the critical behaviour of the theory in the next section when we study the behaviour of the meson spectrum of the theory near the phase transition and show that at the phase transition the theory has a massless mode, which is also approached with a critical exponent .
4 Meson spectrum
In this section we analyse the spectrum of the quadratic fluctuations of the system corresponding to the spectrum of the mesonic bound states in the dual field theory. To obtain the equations of motion for the fluctuations we expand near the classical profile of the D5brane:
(26) 
where parameterise the transverse to the D5brane and the perturbations are allowed to depend on all word volume coordinates. The next step is to expand to second order in the total action of the D5brane:
(27) 
where . Using the equations of motions for the classical embedding of the D5brane, one can show that the first order term in the expansion of the action vanishes, as it should since we are expanding near a local extremum. Defining the matrix and expanding , for the second order expansion of the DBI action quite generally we obtain:
(28) 
On the other hand, since the supergravity background that we consider has a nonzero RR potential , for the second order term of the Wess–Zumino action we obtain:
(29) 
where we have used a gauge in which the magnetic component of the background RR four potential is given by:
(30) 
and that . We see that the Wess–Zumino part of the action couples the fluctuations along and , this effect has been used in the flat case analysed in ref. Filev:2009xp () to show that the Goldstone modes of the theory obey a nonrelativistic linear dispersion relation and to verify the modified counting rule valid for such modes. However, in the present work we will focus on the properties of the theory near the quantum critical point and in particular on the spectrum of fluctuations along , since the spectrum of this mode contains a massless particle at the critical point responsible for the divergent correlation length. To this end we have to write the explicit expressions for the second order action derived from equations (28), (29) and derive the corresponding equations of motion.
4.1 Quadratic Lagrangian
Our results for the second order Lagrangian derived from equations (28), (29) are:
(31)  
(32)  
(33)  
(34)  
(35)  
where are the components of the then dimensional metric and the matrices and are given by:
(36)  
(37) 
Note that in deriving the last expression in equation (35) we have used the equation of motion for and the Bianchi identity for the gauge field (). In the next subsection we use equations (31),(33) and (35) to show that for modes with vanishing angular momentum on the two sphere the equations of motion for and the gauge field decouple. Next we obtain the spectrum of the fluctuations along .
4.2 Spectrum along
In this subsection we analyse the spectrum of the fluctuations along . The reason to focus on this mode is that it contain the massless mode corresponding to the divergent correlation length at the critical point. As one can see from equation (35), the fluctuations along couple to the fluctuations of the gauge field. To obtain the corresponding equations of motion we vary the combined Lagrangian . We observe that if one restricts the fluctuations to the “Smode” by suppressing the angular momentum of the two sphere, the equations of motion for the fluctuations along decouple from those of the gauge field. To find the bound states we consider a consistent ansatz:
(38) 
and solve numerically the resulting second order equation for . One can show that the general solution for at infinity behaves as:
(39) 
Furthermore, substituting the asymptotic solution in the “onshell” quadratic action one can show that only modes with are renormalisable and correspond to bound states with finite energy – this is our quantisation condition (). One can also show that the regular solution for has predetermined boundary conditions for the first derivative at (the point where and the D5brane closes). Therefore we are left with just one parameter specifying the solution, which corresponds to an overall amplitude (the equation is linear) which we set to one. Finally, we define the dimensionless parameters:
(40) 
We are now ready to generate the numerical solution for . For each value of the starting parameter we select those modes for which the coefficient in equation (39) vanishes. This gives us a discrete set of values for corresponding to the meson spectrum associated to the fluctuations along . As one can expect, the spectrum depends crucially on the magnitude of the magnetic field . Plots of the spectrum of the first three states as a function of the bare mass parameter (at fixed ) are presented in figures 13, 14, 15.
One can see that for large bare mass parameter (and at vanishing magnetic field) the spectrum approaches the flat result first obtained in refs. Arean:2006pk (), Myers:2006qr ():
(41) 
which is represented by dashed lines in the figures. Another interesting feature of the spectrum is that at