About this Journal Submit a Manuscript Table of Contents
Advances in High Energy Physics
Volume 2011 (2011), Article ID 917127, 17 pages
doi:10.1155/2011/917127
Research Article

A Note on Holography and Phase Transitions

1UPMC Paris 06, UMR 7589, LPTHE, 75005 Paris, France
2CNRS, UMR 7589, LPTHE, 75005 Paris, France
3Department of Physics, Northeastern University, Boston, MA 02115, USA
4Departamento de Física, Universidad Nacional de La Plata, C.C. 67, 1900 La Plata, Argentina

Received 23 February 2011; Revised 18 June 2011; Accepted 1 July 2011

Academic Editor: Anastasios Petkou

Copyright © 2011 Marc Bellon et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Focusing on the connection between the Landau theory of second-order phase transitions and the holographic approach to critical phenomena, we study diverse field theories in an anti de Sitter black hole background. Through simple analytical approximations, solutions to the equations of motion can be obtained in closed form which give rather good approximations of the results obtained using more involved numerical methods. The agreement we find stems from rather elementary considerations on perturbation of Schrödinger equations.

1. Introduction

Much activity has been centered in the last few years on the application of the gauge/gravity duality, originally emerged from string theory [13], to analyze strongly interacting field theories by mapping them to classical gravity. More recently such duality was successfully applied to describe, holographically, systems undergoing phase transitions like superconductors, superfluids, and other strongly interacting systems [417] (for a wider list of references see e.g., [18]).

In the holographic approach to the study of phase transitions one starts, on the gravity side, with a field theory in asymptotically anti de Sitter (AdS) space-time with temperature arising from a black hole metric, either introduced as a background or resulting from back reaction of matter on the geometry. Then, using the AdS/CFT correspondence, one can study the behavior of the dual field theory defined on the boundary, identify order parameters, and analyze the phase structure of the dual system.

Several models with scalar and gauge fields in a bulk which corresponds asymptotically to an AdS black hole metric have been studied [417]. The rather involved systems of nonlinear coupled differential equations that have to be solved require general application of numerical methods. In this way one often finds nontrivial hairy solutions that cease to exist for 𝑇 > 𝑇 𝑐 , where 𝑇 is the Hawking temperature associated with the black hole and 𝑇 𝑐 a critical temperature depending on the space-time dimensions and the parameters of the model on the bulk.

The main point in this AdS/CFT based calculation is that the asymptotics of the solution in the bulk encodes the behavior of the QFT at finite temperature defined on the border. One finds in general a typical scenario of second-order phase transition. The critical exponents can be computed with a rather good precision and coincide with those obtained within the mean field approximation in a great variety of models.

It is the purpose of this paper to get an analytical insight complementing the numerical results, with a focus on the connection between the holographic approach and mean field theory for the calculation of critical exponents. To this end, it will be important to first stress some connections, already signaled in [6], between the Landau phenomenological theory of second-order phase transitions and the gauge/gravity approach. We will expose the importance of analyticity to explain the similitude in the results for relevant physical quantities (critical exponents, dependence of the critical temperature on the charge density, and magnetic field, etc.) in very diverse models. This will be done first by stressing the relationship between the equations of motion and a Schrödinger problem, so that usual perturbative techniques allow to prove the critical behavior. We then show that simple matching conditions lead to results that broadly agree with elaborate numerical calculations.

It should be noted that our calculations assume that the backreaction of dynamical fields is negligible (probe approximation) which is valid when the gauge coupling constant is large. This approximation is useful to study the behavior near the phase transition which is precisely the domain we will analyze, comparing our analytical results with those obtained numerically. In fact, the holographic results we have previously obtained solving numerically the Einstein-Yang-Mills-Higgs equations of motion both considering the probe approximation [7] and the case in which matter backreacts on the geometry [5] show that when parameters are such that condensation takes place the solutions are very similar and the nature of the phase transition is identical.

With this in mind and using an analytic approach proposed in [8], we study the equations of motion for different models, showing how one can determine in a very simple way the critical behavior of the systems defined on the border with a good agreement with numerical results.

2. Holography and Mean Field Results

In the Landau approach to second-order phase transitions one considers an order parameter 𝒪 and assumes analyticity of the free energy 𝐹 = 𝐹 ( 𝑇 , 𝒪 ) , which can be expanded in even powers of 𝒪 𝐹 ( 𝑇 , 𝒪 ) = 𝐹 = 𝐹 0 + 𝐹 2 [ 𝑇 ] 𝒪 2 + 𝐹 4 [ 𝑇 ] 𝒪 4 + . ( 2 . 1 ) The dependence of the order parameter on 𝑇 is obtained by minimizing 𝐹 as a function of 𝒪 . The next step is to expand the coefficients 𝐹 2 and 𝐹 4 , in powers of 𝑇 𝑇 𝑐 . To ensure stability and a change of behavior at 𝑇 𝑐 , one takes 𝐹 2 = 𝑎 ( 𝑇 𝑇 𝑐 ) and 𝐹 4 = 𝑏 / 4 with 𝑎 , 𝑏 positive constants. In this way, one finds that the minimum of 𝐹 is at 𝒪 = 0 for 𝑇 > 𝑇 𝑐 (disordered phase) while for 𝑇 < 𝑇 𝑐 one has 𝑇 𝒪 𝑐 𝑇 1 / 2 . ( 2 . 2 )

Let us give a brief description of gauge/gravity approach to phase transitions to connect it to the Landau theory. On the gravity side one considers a classical field theory in a Schwarzschild-AdS black hole background (or one in which backreaction of fields on an asymptotically AdS space leads to a black hole solution). The choice of such a geometry is dictated by the fact that the warped AdS geometry prevents massive charged particles to be repelled to the boundary by a charged horizon and as a result a condensate floating over the horizon can be formed. To find such condensate one should look for nontrivial static solutions for the fields outside the black hole by imposing appropriate boundary conditions. The behavior of fields at infinity then allows to determine the dependence of the order parameter on temperature, identified with the Hawking temperature of the black hole.

Basically, the free energy 𝐹 is identified with the minimum of the action in the gravitational theory with prescribed boundary condition for the fields. Regularity of the fields at the horizon and smoothness of the background geometry, basic assumptions within the gauge/gravity duality, imply analyticity of this action. This is a first point of contact with Landau theory of second-order phase transitions leading to a mean field behavior. We will now see that if continuity and smoothness of fields are imposed in the region between the horizon and infinity, the mean field behavior found using a numerical approach is reproduced with a good precision, as can be seen following a very simple analytic approach proposed in [8] which we apply below for different models.

3. The Abelian Higgs Model in 𝑑 Space Dimensions

Dynamics of the system is governed by the action 𝑑 𝑆 = 𝑑 + 1 𝑥 | | 𝑔 | | 1 4 𝐹 𝜇 𝜈 𝐹 𝜇 𝑣 | | 𝜇 𝑖 𝐴 𝜇 Ψ | | 2 𝑚 2 | | Ψ | | 2 . ( 3 . 1 ) The background metric is the standard A d S 𝑑 + 1 -Schwarschild black hole 𝑑 𝑠 2 = 𝑓 ( 𝑟 ) 𝑑 𝑡 2 + 1 𝑓 ( 𝑟 ) 𝑑 𝑟 2 + 𝑟 2 𝑑 𝑥 𝑖 𝑑 𝑥 𝑖 , 𝑖 = 1 , 2 , . . . , 𝑑 1 , 𝑓 ( 𝑟 ) = 𝑟 2 𝑟 1 𝑑 𝑟 𝑑 , ( 3 . 2 ) with 𝑖 = 1 , 2 , , 𝑑 1 and we have chosen the AdS radius to be unity. Different values for the mass term can be considered [9, 10]: 𝑚 2 = 0 , 2 and 𝑚 2 = 𝑚 2 B F = 9 / 4 in 𝑑 = 3 and 𝑚 2 = 0 , 3 , 4 in 𝑑 = 4 , the last value corresponding to the Breitenlohner-Freedman bound 𝑚 2 B F < 0 marking the boundary of stability for a scalar field in AdS [19]. The black hole temperature is given by 𝑑 𝑇 = 𝑟 4 𝜋 . ( 3 . 3 )

In order to look for simple classical solutions for this model one can propose the following ansatz: | | Ψ | | Ψ = = 𝜓 ( 𝑟 ) , 𝐴 𝜇 = 𝜙 ( 𝑟 ) 𝛿 𝜇 0 . ( 3 . 4 ) Imposing regularity of the solution at the horizon ( 𝑟 = 𝑟 ) one gets 𝜓 𝑟 = 𝑑 𝑟 𝑚 2 𝜓 𝑟 , 𝜙 𝑟 = 0 . ( 3 . 5 ) Concerning the asymptotic behavior of the scalar potential 𝜙 and the scalar field 𝜓 one has 𝜓 𝜓 = 𝑟 𝜆 + 𝜓 + 𝑟 𝜆 + , 𝜌 𝜙 = 𝜇 𝑟 𝑑 2 + , ( 3 . 6 ) with 𝜆 ± = 1 2 𝑑 ± 𝑑 2 + 4 𝑚 2 . ( 3 . 7 ) According to the gauge/gravity correspondence, 𝜇 corresponds to the chemical potential in the dual theory defined on the boundary and 𝜌 to the charge density. Concerning the scalar field 𝜓 , both falloffs are acceptable provided the following condition holds [20]: 𝑑 2 4 < 𝑚 2 𝑑 < 2 4 + 1 , ( 3 . 8 ) One can then pick either 𝜓 + = 0 or 𝜓 = 0 leaving a one parameter family of solutions so that the condensate of the dual operator 𝒪 𝜎 will be given by 𝒪 𝜎 = 𝜓 𝜎 , ( 3 . 9 ) with the boundary condition written as 𝜀 𝜎 𝜌 𝜓 𝜌 = 0 , ( 3 . 1 0 ) where 𝜎 , 𝜌 = + , and 𝜀 + = 1 . For definiteness we will impose the condition 𝜓 = 0 .

It will be convenient for the analysis that follows to change variables according to 𝑟 𝑧 = 𝑟 , ( 3 . 1 1 ) so that the horizon is fixed at 𝑧 = 1 and the boundary at 𝑧 = 0 . In terms of this new variable 𝑧 , the equations for the functions in ansatz (3.4) become 𝜓 𝑑 1 + 𝑧 𝑑 𝑧 1 𝑧 𝑑 𝜓 + 𝜙 2 𝑟 2 1 𝑧 𝑑 2 𝑚 2 𝑧 2 1 𝑧 𝑑 𝜙 𝜓 = 0 , ( 3 . 1 2 ) 𝑑 3 𝑧 𝜙 2 𝜓 2 𝑧 2 1 𝑧 𝑑 𝜙 = 0 , ( 3 . 1 3 ) where the prime denotes 𝑑 / 𝑑 𝑧 .

Conditions (3.5) and (3.6) now read 𝜓 𝑚 ( 1 ) = 2 𝑑 𝜙 𝜓 𝜓 ( 1 ) = 0 , ( 1 ) = 0 , ( 3 . 1 4 ) 𝐵 ( 𝑧 ) 𝑧 0 𝐷 + 𝑧 𝜆 + + 𝐷 𝑧 𝜆 , 𝜙 𝐵 ( 𝑧 ) 𝑧 0 𝜇 𝑞 𝑧 𝑑 2 , ( 3 . 1 5 ) with 𝐷 ± = 𝑟 𝜆 ± 𝜓 ± and 𝑞 = 𝜌 / 𝑟 𝑑 2 .

In view of (3.3), the system (3.12)-(3.13) only depends on the black-hole temperature 𝑇 through the nonlinear 𝜙 2 𝜓 term in the scalar field equation (3.12): this describes the coupling of the scalar to the electric potential and gives an effective negative mass.

In the limit 𝑇 the nonlinear term in (3.12) vanishes so that the equation becomes linear with no nontrivial soliton solutions. Now, as already pointed out in [11] the coupling of the scalar to the Maxwell field is responsible for producing a negative effective mass for 𝜓 , and this effect becomes more important at low temperatures. This indicates that one should expect an instability taking place at some point towards forming scalar hair. At such point, the stabilizing effect of the 𝑚 -term is overcome. If a nontrivial solution exists at low temperatures, with an asymptotic behavior associated with a nonzero v.e.v. of an order parameter 𝑂 𝜓 in the dual theory, it is natural to expect that a critical temperature 𝑇 𝑐 should exist such that 𝑂 𝜓 𝑇 > 𝑇 𝑐 = 0 . This was indeed verified numerically for different values of the mass 𝑚 (including the Breitenlohner-Freedman bound value) and various space-time dimensions [9, 10].

Now, (3.12) for fixed potential 𝜙 can be viewed as a Schrödinger equation with zero energy: in order to have nontrivial solutions, the resulting potential has to adjust itself for the existence of a unique zero eigenvalue for every temperature below the critical one. In the vicinity of the critical point, the variation of the potential is linear in the temperature variation, but quadratic in the normalization of the field through the effect of the nonlinear term in (3.13). The eigenvalue of the Schrödinger equation is, to first order, simply proportional to some integral of the potential as can be seen from its variational evaluation: with a suitable weight function 𝑎 ( 𝑧 ) , the integral of 𝑎 ( 𝑧 ) 𝜓 ( 𝑧 ) times the left hand side of (3.12) gives the eigenvalue for a normalized 𝜓 and its variations are of second order with respect to the variations of 𝜓 . Fixing 𝜓 to the eigenfunction for 𝑇 = 𝑇 𝑐 , the two sources of the variation of this integral through the potential must compensate themselves, and one obtains a Δ 𝑇 𝜙 2 relation. It is precisely this type of behavior that has been found from the numerical solution to the differential equations on the gravity side.

To go further in the analysis without resorting to a numerical analysis, we will consider expansions of the fields in the bulk near 𝑧 = 1 and 𝑧 = 0 . Imposing the conditions of continuity and smoothness of the solutions at a point 𝑧 𝑚 intermediate between the boundary ( 𝑧 = 0 ) and the horizon ( 𝑧 = 1 ) will give algebraic equations between the parameters of the solution.

For the solution near the horizon ( 𝑧 = 1 ) we have, up to order ( 𝑧 1 ) 2 , 𝜓 𝐻 ( 𝑧 ) = 𝜓 0 + 𝜓 1 1 ( 𝑧 1 ) + 2 𝜓 2 ( 𝑧 1 ) 2 , 𝜙 𝐻 ( 𝑧 ) = 𝜙 0 + 𝜙 1 1 ( 𝑧 1 ) + 2 𝜙 2 ( 𝑧 1 ) 2 , ( 3 . 1 6 ) with 𝜓 0 , 𝜓 1 , 𝜓 2 , 𝜙 0 , 𝜙 1 , 𝜙 2 constants. The boundary conditions (3.14) at 𝑧 = 1 imply 𝜓 1 𝑚 = 2 𝑑 𝜓 0 , 𝜙 0 = 0 . ( 3 . 1 7 ) Substituting these values in (3.16) and using the differential (3.12)-(3.13) we can obtain 𝜙 2 and 𝜓 2 as a function of 𝜙 1 and 𝜓 0 . We get 𝜓 𝐻 ( 𝑧 ) = 𝜓 0 + 𝑚 2 𝑑 𝜓 0 ( 1 𝑧 ) + 2 𝑑 𝑚 2 𝑟 2 + 𝑚 4 𝑟 2 𝜙 2 1 4 𝑑 2 𝑟 2 𝜓 0 ( 1 𝑧 ) 2 , 𝜙 𝐻 ( 𝑧 ) = 𝜙 1 1 ( 1 𝑧 ) + 2 𝑑 3 2 𝜓 2 0 𝑑 𝜙 1 ( 1 𝑧 ) 2 . ( 3 . 1 8 )

As announced, imposing the conditions of continuity and smoothness at an intermediate point 𝑧 𝑚 allows to obtain a solution. Interestingly enough, as first observed in [8] the result of this crude approximation is quite stable with respect to the intermediate point, so we will consider the case 𝑧 𝑚 = 1 / 2 . This can be understood from the fact that 𝜓 is the ground state of a Schrödinger equation, so that it cannot have nodes: additional terms in the expansion of 𝜓 should rapidly fade away.

We will first analyze the boundary condition 𝐷 = 0 , so we have the set of equations 𝜓 𝐻 1 2 = 𝜓 𝐵 1 2 , 𝜓 𝐻 1 2 = 𝜓 𝐵 1 2 , 𝜙 𝐻 1 2 = 𝜙 𝐵 1 2 , 𝜙 𝐻 1 2 = 𝜙 𝐵 1 2 , ( 3 . 1 9 ) to be solved for 𝜓 0 , 𝜙 1 , 𝐷 + , and 𝜇 . We obtain 𝜓 2 0 = 𝑑 2 𝑑 + 1 1 6 ( 𝑑 2 ) 𝐵 𝑞 + 2 𝑑 ( 𝑑 5 ) 𝑟 𝐴 𝑟 𝐴 , 𝐷 + = 2 𝜆 + 1 4 𝑑 + 𝑚 2 𝑑 𝐵 2 𝜓 0 , 𝜙 1 𝑟 = 𝐴 𝐵 , 1 𝜇 = 2 𝑑 + 2 8 𝑑 𝑞 𝐵 + 2 𝑑 𝑟 𝐴 𝐵 ( 3 . 2 0 ) (for simplicity 𝐷 + is written in terms of 𝜓 0 ). Here 𝐴 = 1 6 𝑑 2 𝜆 + + 𝑚 4 𝜆 + + 2 + 2 𝑑 𝑚 2 6 + 5 𝜆 + , 𝐵 = 𝜆 + + 2 . ( 3 . 2 1 ) Using the AdS/CFT dictionary (3.9), we can identify the v.e.v. 𝒪 of the operator 𝒪 dual to the scalar field with the asymptotic coefficient 𝜓 + , 𝒪 𝜓 + = 𝑟 𝜆 + 𝐷 + . Now, one can write 𝐷 + (or 𝜓 0 since they are proportional) as a function of 𝑇 = ( 𝑑 / 4 𝜋 ) 𝑟 . Remembering that 𝑞 = 𝜌 / 𝑟 𝑑 2 one has 𝜓 2 0 = 𝑑 ( 5 𝑑 ) 2 𝑇 𝑐 𝑇 𝑑 1 𝑇 1 𝑇 𝑐 𝑑 1 , ( 3 . 2 2 ) where we have defined 𝑇 𝑐 𝑑 1 = 2 4 𝑑 ( 𝑑 2 ) ( 5 𝑑 ) ( 4 𝜋 / 𝑑 ) 𝑑 1 𝐵 𝐴 𝜌 . ( 3 . 2 3 ) We then have, for the order parameter 𝒪 = 𝐶 𝑑 𝑇 𝑇 𝑐 ( 𝑑 2 + 4 𝑚 2 + 1 ) / 2 𝑇 1 𝑇 𝑐 𝑑 1 , ( 3 . 2 4 ) with 𝐶 𝑑 = 4 𝑑 + 𝑚 2 𝜆 + 𝑑 + 2 2 𝜋 𝜆 + ( 5 𝑑 ) 8 𝑑 1 / 2 . ( 3 . 2 5 )

One can see from (27) that for 𝑇 close to 𝑇 𝑐 one has the typical second-order phase transition behavior 𝒪 ( 1 𝑇 ) / 𝑇 𝑐 . Note that 𝑇 𝑐 𝜌 1 / ( 𝑑 1 ) in agreement with the change of dimensions of the charge density for different 𝑑 ’s.

Our results for the critical temperature, as inferred from (3.23)–(3.25) in the cases 𝑑 = 3 , 4 , 𝑚 2 = 2 (with 𝑧 𝑚 = 1 / 2 ) are 𝑇 𝑐 = 0 . 1 5 𝜌 1 / 2 and 𝑇 c = 0 . 2 𝜌 1 / 3 , respectively. They can be compared with those obtained using a different analytical approximation based on perturbation theory near the critical temperature. For the case 𝑚 = 𝑚 B F , we obtain 𝑇 𝑐 = 0 . 1 2 𝜌 1 / 2 and 𝑇 𝑐 = 0 . 2 5 𝜌 1 / 3 [21, 22] in very good agreement with the exact numerical results, 𝑇 𝑐 = 0 . 1 5 𝜌 1 / 2 and 𝑇 𝑐 = 0 . 2 5 𝜌 1 / 3 . One can conclude that there is a good quantitative agreement between the three sets, which is not much affected by the choice of the point 𝑧 𝑚 in the method. This last fact was already observed in [8] for the particular case 𝑑 = 3 compared with the numerical results given in [9].

4. Non-Abelian Gauge Field in 𝑑 = 3 Space Dimensions

We now consider the case of an S U ( 2 ) Yang-Mills theory in 3 + 1 dimensional Antide Sitter-Schwarzschild background, as a prototype for the gauge/gravity duality in the case of pure gauge theories [12, 13]. The action is 1 𝑆 = 4 𝑑 4 𝑥 | | 𝑔 | | 𝐹 𝑎 𝜇 𝜈 𝐹 𝑎 𝜇 𝜈 , ( 4 . 1 ) with the field strength defined as 𝐹 𝑎 𝜇 𝜈 = 𝜕 𝜇 𝐴 𝑎 𝜈 𝜕 𝜈 𝐴 𝑎 𝜇 + 𝑔 𝜖 𝑎 𝑏 𝑐 𝐴 𝑏 𝜇 𝐴 𝑐 𝜈 , 𝑎 = 1 , 2 , 3 . ( 4 . 2 ) Writing the gauge field as an isospin vector 𝐴 𝐴 = 1 , 𝐴 2 , 𝐴 3 = 𝐴 𝑎 𝜇 𝑑 𝑥 𝜇 , ( 4 . 3 ) we consider the following ansatz for solving the equations of motion: 𝐴 = 𝐽 ( 𝑟 ) 𝑑 𝑡 𝑒 3 𝐾 ( 𝑟 ) 𝜌 𝑑 𝜑 𝑒 𝜌 + 𝐾 ( 𝑟 ) 𝑑 𝜌 𝑒 𝜑 , ( 4 . 4 ) where 𝑟 is the radial variable in spherical coordinates and 𝑥 = 𝜌 c o s 𝜑 , 𝑦 = 𝜌 s i n 𝜑 . ( 4 . 5 ) The background metric is 𝑑 𝑠 2 = 𝑓 ( 𝑟 ) 𝑑 2 1 𝑡 + 𝑓 𝑑 ( 𝑟 ) 2 𝑟 + 𝑟 2 𝑑 𝑥 2 + 𝑑 𝑦 2 , ( 4 . 6 ) with 𝑓 ( 𝑟 ) = 𝑟 2 𝑟 1 3 𝑟 3 . ( 4 . 7 ) In terms of the 𝑧 variable defined in (3.11) the equations of motion read 1 𝑧 3 𝐾 = 1 𝑟 2 𝐾 2 𝐽 2 𝑧 2 1 𝑧 3 𝐽 𝐾 , 2 ( 𝑧 ) = 𝑟 2 𝐽 ( 𝑧 ) 𝐾 2 ( 𝑧 ) 1 𝑧 3 . ( 4 . 8 )

At the horizon ( 𝑧 = 1 ), 𝐽 must vanish, so we have the following expansions for the fields, up to order ( 𝑧 1 ) 2 , 𝐾 𝐻 ( 𝑧 ) = 𝐾 0 + 𝐾 1 1 ( 𝑧 1 ) + 2 𝐾 2 ( 𝑧 1 ) 2 , 𝐽 𝐻 ( 𝑧 ) = 𝐽 1 1 ( 𝑧 1 ) + 2 𝐽 2 ( 𝑧 1 ) 2 . ( 4 . 9 ) Using the equations of motion (4.8), the coefficients 𝐾 1 , 𝐾 2 , and 𝐽 2 can be written in terms of 𝐾 0 and 𝐽 1 as 𝐾 𝐻 ( 𝑧 ) = 𝐾 0 𝐾 3 0 3 𝑟 2 ( 𝑧 1 ) + 𝐽 2 1 𝐾 0 𝑟 2 + 6 𝐾 3 0 𝑟 2 + 3 𝐾 5 0 3 6 𝑟 4 ( 𝑧 1 ) 2 , 𝐽 𝐻 ( 𝑧 ) = 𝐽 1 𝐽 ( 𝑧 1 ) 1 𝐾 2 0 3 𝑟 2 ( 𝑧 1 ) 2 . ( 4 . 1 0 )

At the 𝑧 = 0 boundary one has the asymptotic expansions 𝐾 𝐵 𝐶 ( 𝑧 ) = 1 𝑟 𝑧 , 𝐽 𝐵 ( 𝑧 ) = 𝐷 0 𝐷 1 𝑟 𝑧 , ( 4 . 1 1 ) where 𝐷 1 can be associated with the charge density. Coefficient 𝐶 1 in (4.11) should be identified with the order parameter 𝒪 𝐾 for the theory on the border. Since the order parameter is related to a vector field ( 𝐴 𝑖 ), the associated theory on the border is a p-wave superconductor [13].

As in the scalar field case, we will match both solutions at an intermediate point which we again choose as 𝑧 𝑚 = 1 / 2 , 𝐾 𝐻 1 2 = 𝐾 𝐵 1 2 , 𝐾 𝐻 1 2 = 𝐾 𝐵 1 2 , 𝐽 𝐻 1 2 = 𝐽 𝐵 1 2 , 𝐽 𝐻 1 2 = 𝐽 𝐵 1 2 , ( 4 . 1 2 ) and solve these equations for 𝐾 0 , 𝐽 1 , 𝐶 1 , and 𝐷 0 . From the first two identities (4.12), we obtain 𝐶 1 = 4 3 𝑟 𝐾 0 + 𝐾 3 0 9 𝑟 , 𝐽 1 = 4 8 𝑟 2 + 2 2 𝐾 2 0 𝐾 + 3 4 0 / 𝑟 2 ( 4 . 1 3 ) (we chose 𝐽 1 < 0 so 𝐽 ( 𝑧 ) > 0 ).

Substituting these values in (4.12), we can obtain 𝐾 0 as the root of the following polynomial in 𝐾 2 0 : 3 𝑟 2 𝐾 8 0 + 4 0 𝐾 6 0 + 2 0 7 𝑟 2 𝐾 4 0 + 4 8 6 𝑟 4 𝐾 2 0 9 𝑟 2 𝐷 2 1 4 8 𝑟 4 = 0 . ( 4 . 1 4 )

This equation implies that 𝐾 2 0 vanishes as 𝐾 2 0 = 𝐷 2 1 4 8 𝑟 4 5 4 𝑟 2 𝐷 + 𝑂 2 1 4 8 𝑟 4 2 . ( 4 . 1 5 )

Finally, introducing the temperature 𝑇 = 3 / ( 4 𝜋 ) 𝑟 , we can write 𝐾 2 0 = 1 2 8 𝜋 2 𝑇 8 1 4 𝑐 𝑇 2 𝑇 1 𝑐 𝑇 4 f o r 𝑇 n e a r 𝑇 𝑐 , ( 4 . 1 6 ) where 𝑇 2 𝑐 = 3 3 6 4 𝜋 2 𝐷 1 . ( 4 . 1 7 ) Then, for 𝑇 close to 𝑇 𝑐 , we have the typical second-order phase transition behavior for 𝐶 1 = 𝒪 𝐾 ( 1 𝑇 ) / 𝑇 𝑐 , in good agreement with [12, 13]. Our numerical value for 𝑇 𝑐 from (4.17), 𝑇 𝑐 / 𝐷 1 = 0 . 0 9 1 , can be compared with the numerical value obtained in [13], 𝑇 𝑐 / 𝐷 1 = 0 . 1 2 5 .

5. The 𝑑 = 3 Scalar Case in the Presence of an Applied Magnetic Field

We will consider here a system with dynamics governed by the action (3.1) in the 𝑑 = 3 , 𝑚 2 = 2 case, when an external magnetic field 𝐻 is applied. This corresponds to the case where the background metric is an AdS 3 + 1   magnetically charged black hole [14, 15]: 𝑑 𝑠 2 = 𝑓 ( 𝑟 ) 𝑑 𝑡 2 + 1 𝑓 ( 𝑟 ) 𝑑 𝑟 2 + 𝑟 2 𝑑 𝑥 𝑖 𝑑 𝑥 𝑖 𝑟 , 𝑖 = 1 , 2 , 𝑓 ( 𝑟 ) = 2 𝐿 2 𝑀 𝑟 + 𝐻 2 𝑟 2 . ( 5 . 1 ) The 𝑓 ( 𝑟 ) = 0 condition for having event horizons gives a quartic algebraic equation with 4 roots which can be explicitly written in terms of surds. There are two complex conjugate roots which we will call 𝑟 1 and 𝑟 2 and two real roots 𝑟 3 < 𝑟 4 . We will then call 𝑟 𝑟 4 the external black hole horizon. Taking 𝐿 = 1 from here on, one has 𝑟 = 1 2 6 1 / 3 𝐵 + 1 2 𝑀 𝐵 B 2 , ( 5 . 2 ) where 𝐵 = 8 3 1 / 3 𝐻 2 𝐷 + 9 𝑀 2 1 / 3 + 2 1 / 3 𝐷 + 9 𝑀 2 1 / 3 , 𝐷 = 8 1 𝑀 4 7 6 8 𝐻 6 . ( 5 . 3 ) The actual form of the other roots 𝑟 𝑎 ( 𝑎 = 1 , 2 , 3 ) is not necessary since one will only need the standard relationships between roots and coefficients 𝑟 1 + 𝑟 2 + 𝑟 3 = 𝑟 , 𝑟 1 𝑟 2 + 𝑟 1 𝑟 3 + 𝑟 2 𝑟 3 = 𝑟 2 , 𝑟 1 𝑟 2 𝑟 3 = 𝐻 2 𝑟 , ( 5 . 4 ) together with the relation 𝐻 2 𝑟 4 𝑀 𝑟 3 + 1 = 0 . ( 5 . 5 ) The Hawking temperature associated to the black hole is 𝑓 𝑇 = 𝑟 = 1 4 𝜋 4 𝜋 3 𝑟 2 𝑟 3 . ( 5 . 6 )

In terms of the variable 𝑧 = 𝑟 / 𝑟 , the equations of motion read 𝜓 𝑟 ( 𝑧 ) + 𝑑 𝑓 / 𝑧 / 𝑑 𝑧 𝑓 𝑟 / 𝑧 𝜓 𝑟 ( 𝑧 ) + 2 𝑧 4 𝜙 2 ( 𝑧 ) 𝑓 𝑟 / 𝑧 2 + 2 𝑓 𝑟 𝜙 / 𝑧 = 0 , ( 𝜓 𝑧 ) 2 2 ( 𝑧 ) 𝑓 𝑟 𝑟 / 𝑧 𝜙 ( 𝑧 ) 2 𝑧 4 = 0 , ( 5 . 7 ) where 𝜓 ( 𝑧 ) 𝜓 ( 𝑟 / 𝑧 ) and 𝜙 ( 𝑧 ) 𝜙 ( 𝑟 / 𝑧 ) (but from now on, the tilde will be omitted). In terms of the new variables, we have 𝑓 𝑟 𝑧 = 𝑟 𝑧 2 ( 1 𝑧 ) 1 ̃ 𝑟 1 𝑧 1 ̃ 𝑟 2 𝑧 1 ̃ 𝑟 3 𝑧 , 𝑟 𝑑 𝑓 / 𝑧 / 𝑑 𝑧 𝑓 𝑟 = 𝑧 / 𝑧 2 𝑀 3 / 𝑟 3 + 2 𝐻 2 𝑧 4 / 𝑟 4 𝑧 ( 1 𝑧 ) 1 ̃ 𝑟 1 𝑧 1 ̃ 𝑟 2 𝑧 1 ̃ 𝑟 3 𝑧 , ( 5 . 8 ) where ̃ 𝑟 𝑎 = 𝑟 𝑎 / 𝑟 , 𝑎 = 1 , 2 , 3 .

The boundary conditions for system (5.7) at the horizon are 𝜓 𝜙 ( 1 ) = 0 , ( 5 . 9 ) ( 1 ) = 2 𝑟 3 𝑟 𝐻 2 / 𝑟 3 𝜓 ( 1 ) , ( 5 . 1 0 ) while asymptotically one has 𝜓 ( 𝑧 ) = 𝐷 1 𝑧 + 𝐷 2 𝑧 2 , 𝜌 𝜙 ( 𝑧 ) = 𝜇 𝑟 𝑧 . ( 5 . 1 1 )

For the solution of system (5.7) near the horizon, we have the expansions, up to order ( 𝑧 1 ) 2 , 𝜓 𝐻 ( 𝑧 ) = 𝜓 0 + 𝜓 1 1 ( 𝑧 1 ) + 2 𝜓 2 ( 𝑧 1 ) 2 , 𝜙 𝐻 ( 𝑧 ) = 𝜙 0 + 𝜙 1 1 ( 𝑧 1 ) + 2 𝜙 2 ( 𝑧 1 ) 2 , ( 5 . 1 2 ) with 𝜓 0 , 𝜓 1 , 𝜓 2 , 𝜙 0 , 𝜙 1 , and 𝜙 2 constants. Using (5.10) we get 𝜙 0 𝜓 = 0 , 1 = 2 𝑟 3 𝑟 𝐻 2 / 𝑟 3 𝜓 0 . ( 5 . 1 3 ) Substituting these values in (5.12) and using the differential equations, we can obtain 𝜙 2 and 𝜓 2 as functions of 𝜙 1 and 𝜓 0 . We get 𝜓 ( 𝑧 ) = 𝜓 0 + 2 𝑟 4 𝐻 2 3 𝑟 4 𝜓 0 𝑟 ( 1 𝑧 ) 4 8 𝑟 4 + 𝑟 2 𝜙 2 1 1 2 𝐻 2 4 𝐻 2 3 𝑟 4 2 𝜓 0 ( 1 𝑧 ) 2 , 𝜙 ( 𝑧 ) = 𝜙 1 𝑟 ( 1 𝑧 ) + 4 𝜓 2 0 𝜙 1 𝐻 2 3 𝑟 4 ( 1 𝑧 ) 2 . ( 5 . 1 4 )

As before, we impose matching conditions at 𝑧 = 1 / 2 : 𝜓 𝐻 1 2 = 𝜓 𝐵 1 2 , 𝜓 𝐻 1 2 = 𝜓 𝐵 1 2 , 𝜙 𝐻 1 2 = 𝜙 𝐵 1 2 , 𝜙 𝐻 1 2 = 𝜙 𝐵 1 2 , ( 5 . 1 5 ) and solve for 𝜓 0 , 𝜙 1 , 𝐷 2 , and 𝜇 . We look for solutions with 𝐷 1 = 0 and 𝐷 2 corresponding to the order parameter of the 2 + 1 system defined on the boundary. We obtain 𝜙 1 𝑟 = 2 𝑅 , 𝐻 , 𝐻 𝜇 = 2 𝜌 / 𝑟 𝑟 + 2 𝑅 , 𝐻 𝑟 4 3 𝜌 / 𝑟 𝑟 + 𝑅 𝜓 , 𝐻 2 0 + 6 2 𝐻 2 3 𝑟 4 , 𝜓 2 0 = 3 𝑟 4 𝐻 2 𝜌 / 𝑟 𝑟 2 𝑅 , 𝐻 2 𝑟 4 𝑅 𝑟 , 𝐷 , 𝐻 2 = 8 8 𝑟 8 𝜙 2 1 𝑟 6 6 8 𝐻 2 𝑟 4 + 1 6 𝐻 4 𝜓 0 4 𝐻 2 3 𝑟 4 2 , ( 5 . 1 6 ) where 𝑅 𝑟 = , 𝐻 7 𝑟 8 6 𝐻 2 𝑟 4 + 2 𝐻 4 𝑟 6 . ( 5 . 1 7 )

The equation for 𝜓 2 0 in terms of the dimensionless variable 𝑢 = 𝑟 / 𝐻 takes the form 𝜓 2 0 = 3 𝑢 4 1 ( 𝜌 / 𝐻 ) 2 2 / 𝑢 4 6 + 7 𝑢 4 2 𝑢 4 2 / 𝑢 4 6 + 7 𝑢 4 , ( 5 . 1 8 ) with the temperature (5.6) given by 𝑇 𝐻 = 1 1 4 𝜋 3 𝑢 𝑢 3 . ( 5 . 1 9 )

The minimum value that 𝑢 can take is the one for which 𝑇 = 0 , 𝑢 = 3 1 / 4 , and corresponds to the condition 3 3 𝑀 4 = 2 7 𝐻 6 . From this, we see that in order to have a nontrivial solution the following inequality should hold: 𝜌 𝐻 2 2 / 𝑢 4 6 + 7 𝑢 4 . ( 5 . 2 0 ) The maximum of the r.h.s. is attained for 𝑢 𝑚 = ( 2 / 7 ) 1 / 8 > 𝑢 0 so that there is a critical value 𝐻 𝑐 of the magnetic field beyond which no nontrivial solution exists, 𝐻 𝑐 = 0 . 4 1 𝜌 . ( 5 . 2 1 )

Using (5.18)-(5.19), one can determine the critical temperature 𝑇 𝑐 as a function of 𝐻 . We give in Figure 1 the resulting 𝑇 𝑐 = 𝑇 𝑐 ( H ) curve. Interestingly, we find that in the range 𝐻 𝑐 𝜌 > 𝐻 > 2 2 / 𝑢 4 0 6 + 7 𝑢 4 0 = 0 . 3 2 7 𝜌 , ( 5 . 2 2 ) the curve 𝑇 𝑐 ( 𝐻 ) becomes double valued so that a nontrivial solution only exists in the range 𝑇 𝑐 2 > 𝑇 > 𝑇 𝑐 1 . According to the gauge/gravity duality 𝐷 2 should be identified with the order parameter, 𝒪 𝜓 = 𝐷 2 . We obtain the following expression: 𝐷 2 = 5 𝑢 4 2 6 𝑢 4 2 ( 𝜌 / 𝐻 ) 2 2 / 𝑢 4 6 + 7 𝑢 4 𝑢 2 2 / 𝑢 4 6 + 7 𝑢 4 1 / 4 . ( 5 . 2 3 )

917127.fig.001
Figure 1: The phase diagram of 𝑇 𝑐 against the magnetic field 𝐻 . The condensed phase ( 𝐷 2 0 ) corresponds to the lower left part below the line. The critical temperature decreases as the magnetic field grow up to the critical value 𝐻 𝑐 .

Note that 𝐷 2 becomes negative for 𝑢 0 < 𝑢 < ( 2 / 5 ) 1 / 4 or, equivalently, for 0 < 𝑇 < 0 . 0 3 1 6 𝐻 . However, when 𝐻 0 the numerator and denominator coefficients of 𝐷 2 that multiply 𝜓 0 cancel out. This is consistent with the result of the nomagnetic field model. (In fact, we have checked that the results of both models, 𝐻 = 0 , 𝑑 = 3 , and 𝑚 2 = 2 , are identical. We present in Figure 2 curves for 𝐷 2 as a function of temperature for different values of the external magnetic field.

917127.fig.002
Figure 2: A plot of the order parameter 𝐷 2 as a function of temperature 𝑇 . The charge density is 𝜌 = 1 . The dashed line corresponds to a magnetic field 𝐻 = 0 . 1 , the dash-dot one to 𝐻 = 0 . 2 , and the solid one to 𝐻 = 0 . 3 . (the critical value is 𝐻 𝑐 = 0 . 6 ).

The results are qualitatively in agreement with those described in [14, 15] where the reported curves obtained by numerically solving the equations of motion are similar to those in Figure 2.

6. Summary and Discussion

We have analyzed a number of models which have been proposed to study phase transition through the AdS/CFT correspondence. The common feature of all three models we discussed was that the space time bulk geometry was an Anti de Sitter black hole. Although the dynamical field content was very different—a charged scalar coupled to an electric potential, the same model in an external magnetic field, and a pure non-abelian gauge theory—the emerging scenarios are very similar and always include a second-order phase transition with mean-field critical exponents.

On general grounds, we were able to explain why the highly symmetric ansätze generally used produce the critical behaviors seen in mean field theory or the Landau approach. Founded on basic principles as the connection between the equations of motion and the Schrödinger equation, we clarify the similarity between several relevant quantities along a variety of models. In particular we showed that resorting to simple matching conditions, we obtain closed form solutions that significantly agree with the results obtained by numerically solving the exact set of equations of motion. This uncovers the important role played by analyticity to explain the universal behavior of certain physical constants.

The method seems to work very well near the critical temperature, though it deviates from the numerical results as we approach 𝑇 0 . In this regime our approach should be refined.

Alternative analytic calculations have been recently presented in [23, 24] where the phase transition vicinity is studied solving the equations of motion in terms of a series expansion near the horizon. Although the approach in these works is close to the one proposed in [8] and applied here, the possibility in the latter of varying the intermediate point 𝑧 𝑚 at which the matching is performed allows to obtain better solutions at fixed order 𝑁 in the expansion, as already pointed out in [24]. As we have seen that, in the matching approach, the problem reduces to find the solution of an algebraic equations system and this can be done, to the order we worked here, in a straightforward way. Increasing the order will of course complicate the algebraic system but in view of its main features, it can be handled by a simple computational software like Mathematica, at least for the next few orders. For a large-order expansion the method followed in references [23, 24] seems to be more appropriate.

Although the matching method works very well near the critical temperature, it deviates from the numerical results as 𝑇 approaches 0. In this regime the method should be refined. In particular it is to be expected that taking into account the quantum fluctuations of the gravity theory, one should be able to go beyond mean field approximation results. Also, one should consider generalized Lagrangians (like the Stückelberg one considered in [16]) leading to various types of phase transitions (first or second order with both mean and nonmean field behavior) as parameters are changed. There is also the possibility that including fermions in the bulk model could substantially change the critical behavior of the theory in the bulk (see [25] and references therein). The simplicity of the approach presented here, not requiring refined numerical calculations, should be an asset when trying to explore these more complex situations.

Acknowledgment

F. A. Schaposnik is associated with CICBA.

References

  1. J. Maldacena, “The large N limit of superconformal field theories and supergravity,” Advances in Theoretical and Mathematical Physics, vol. 2, no. 2, pp. 231–252, 1998. View at Zentralblatt MATH
  2. S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, “Gauge theory correlators from non-critical string theory,” Physics Letters B, vol. 428, no. 1-2, pp. 105–114, 1998. View at Publisher · View at Google Scholar · View at MathSciNet
  3. E. Witten, “Anti de Sitter space and holography,” Advances in Theoretical and Mathematical Physics, vol. 2, no. 2, pp. 253–291, 1998. View at Zentralblatt MATH
  4. S. S. Gubser, “Phase transitions near black hole horizons,” Classical and Quantum Gravity, vol. 22, no. 23, pp. 5121–5143, 2005. View at Publisher · View at Google Scholar · View at MathSciNet
  5. A. R. Lugo, E. F. Moreno, and F. A. Schaposnik, “Holography and AdS4 self-gravitating dyons,” Journal of High Energy Physics, vol. 2010, no. 11, pp. 1–19, 2010. View at Publisher · View at Google Scholar
  6. N. Iqbal, H. Liu, M. Mezei, and Q. Si, “Quantum phase transitions in holographic models of magnetism and superconductors,” Physical Review D, vol. 82, no. 4, Article ID 045002, 2010. View at Publisher · View at Google Scholar
  7. A. R. Lugo, E. F. Morenob, and F. A. Schaposnika, “Holographic phase transition from dyons in an AdS black hole background,” Journal of High Energy Physics, vol. 2010, no. 3, article 013, 2010. View at Publisher · View at Google Scholar
  8. R. Gregory, S. Kanno, and J. Soda, “Holographic superconductors with higher curvature corrections,” Journal of High Energy Physics, vol. 2009, no. 10, article 010, 2009. View at Publisher · View at Google Scholar
  9. S. A. Hartnoll, C. P. Herzog, and G. T. Horowitz, “Building a holographic superconductor,” Physical Review Letters, vol. 101, no. 3, 4 pages, 2008. View at Publisher · View at Google Scholar
  10. G. T. Horowitz and M. M. Roberts, “Holographic superconductors with various condensates,” Physical Review D, vol. 78, no. 12, Article ID 126008, 2008. View at Publisher · View at Google Scholar
  11. S. S. Gubser, “Breaking an Abelian gauge symmetry near a black hole horizon,” Physical Review D, vol. 78, no. 6, Article ID 065034, 2008. View at Publisher · View at Google Scholar
  12. S. S. Gubser, “Colorful horizons with charge in anti-de Sitter space,” Physical Review Letters, vol. 101, no. 19, p. 191601, 4, 2008. View at Publisher · View at Google Scholar
  13. S. S. Gubser and S. S. Pufu, “The gravity dual of a p-wave superconductor,” Journal of High Energy Physics, vol. 2008, no. 11, article 033, 2008. View at Publisher · View at Google Scholar
  14. E. Nakano and W.-Y. Wen, “Critical magnetic field in a holographic superconductor,” Physical Review D, vol. 78, no. 4, Article ID 046004, 2008. View at Publisher · View at Google Scholar
  15. T. Albash and C. V. Johnson, “A holographic superconductor in an external magnetic field,” Journal of High Energy Physics, vol. 2008, no. 9, article 121, 2008. View at Publisher · View at Google Scholar
  16. S. Franco, A. M. García-García, and D. Rodríguez-Gómez, “Holographic approach to phase transitions,” Physical Review D, vol. 81, no. 4, Article ID 041901, 5 pages, 2010.
  17. S. A. Hartnoll, C. P. Herzog, and G. T. Horowitz, “Holographic superconductors,” Journal of High Energy Physics, vol. 2008, no. 12, article 015, 2008. View at Publisher · View at Google Scholar
  18. G. T. Horowitz, “Introduction to holographic superconductors,” Lecture Notes in Physics, vol. 828, pp. 313–347, 2011. View at Publisher · View at Google Scholar
  19. P. Breitenlohner and D. Z. Freedman, “Stability in gauged extended supergravity,” Annals of Physics, vol. 144, no. 2, pp. 249–281, 1982. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  20. I. R. Klebanov and E. Witten, “AdS/CFT correspondence and symmetry breaking,” Nuclear Physics. B, vol. 556, no. 1-2, pp. 89–114, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  21. G. Siopsis and J. Therrien, “Analytic calculation of properties of holographic superconductors,” Journal of High Energy Physics, vol. 2010, no. 5, article 013, 2010. View at Publisher · View at Google Scholar
  22. G. Siopsis, J. Therrien, and S. Musiri, “Holographic conductors near the Breitenlohner-Freedman bound,” http://arxiv.org/abs/1011.2938.
  23. C. P. Herzog, “Analytic holographic superconductor,” Physical Review D, vol. 81, no. 12, Article ID 126009, 11 pages, 2010.
  24. F. Aprile, S. Franco, D. Rodríguez-Gómez, and J. G. Russo, “Phenomenological models of holographic superconductors and Hall currents,” Journal of High Energy Physics, vol. 2010, no. 5, article 102, 2010. View at Publisher · View at Google Scholar
  25. T. Faulkner, N. Iqbal, H. Liu, J. McGreevy, and D. Vegh, “Holographic non-Fermi-liquid fixed points,” Philosophical Transactions of the Royal Society A, vol. 369, no. 1941, pp. 1640–1669, 2011. View at Publisher · View at Google Scholar · View at PubMed