arxiv: v1 [hep-th] 19 Jan 2015

Transcription

1 A holographic model for antiferromagnetic quantum phase transition induced by magnetic field arxiv: v1 [hep-th] 19 Jan 215 Rong-Gen Cai, 1, Run-Qiu Yang, 1, and F.V. Kusmartsev 2, 1 State Key Laboratory of Theoretical Physics,Institute of Theoretical Physics, Chinese Academy of Sciences,Beijing 119, China. 2 Department of Physics, Loughborough University, Loughborough, Leicestershire, LE11 3TU, United Kingdom We propose a gravity dual of antiferromagnetic quantum phase transition (QPT) induced by magnetic field and study the criticality in the vicinity of quantum critical point (QCP). Results show the boundary critical theory is a strong coupling theory with dynamic exponent z = 2. The hyperscaling law is violated and logarithmic corrections appear near the QCP. We compare our theoretical results with experimental data on variety of materials including low-dimensional magnet, BiCoPO 5 and pyrochlores, Er 2 2xY 2xTi 2O 7. Our model describes well the existing experiments and predicts QCP and other high field magnetic properties of these compounds. Introduction. Quantum phase transition (QPT) and the behavior of quantum systems in the vicinity of the corresponding quantum critical point (QCP) have attracted a lot of attention both in theory and in experiment recently [1 3]. In contrast to their classical counterparts induced by thermal fluctuations arising at finite temperature T >, QPTs happen at zero temperature and are governed by quantum fluctuations associated with the Heisenberg uncertainty and driven by a certain control parameter rather at zero temperature, e.g., composition, magnetic field or pressure and commonly separate an ordered sate from a disordered phase. In condensed matter physics, such quantum criticality is considered to play an important role in some of the most interesting phenomena associated with itinerant electronic systems [4, 5]. One of intensively discussed QPTs is ordereddisordered QPT in antiferromagnetic materials induced by magnetic field (see for example, Refs. [6 8].) To study quantum critical behavior in the heavy-fermion systems are very suitable since they can be tuned continuously from an antiferromagnetic (AF) to a paramagnetic (PM) metallic state by varying a single parameter [4]. In these materials, QPT naturally belongs to the phenomenon involving strongly correlated many-body systems [9, 1]. However, the complete theoretical descriptions valid at all the energy (or temperature) region are still lacking. In order to study and characterize strongly coupled quantum critical systems, some new methods are called for. Therefore to develop a proper theoretical description of such a system a new method based on AdS/CFT duality is needed. This duality relates a weak coupling gravitational theory in a (d + 1)-dimensional asymptotically anti de-sitter (AdS) space-time to a d-dimensional strong coupling conformal field theory (CFT) in the AdS boundary [11 13]. In recent years, we have indeed witnessed that the duality has been extensively applied into condensed matter physics and a great deal of progress have been made [14 17]. In Ref. [18] we realized the ferromagnetic/paramagnetic phase transition in a holographic setup, and in Ref. [19] the holographic antiferromagnetic/paramagnetic phase transition was studied. We showed that the antiferroamgnetic transition temperature T N is indeed suppressed by an external magnetic field and tends to zero when the magnetic field reaches its critical value B c. This way the antiferromgnetic QPT induced by magnetic field is realized. In this Letter we will elaborate this QPT and study the critical behavior near the QCP. Holographic model. In order to describe spontaneous staggered magnetization which breaks the time reversal symmetry, we introduce two antisymmetric tenser fields coupling with U(1) Maxwell field strength. The bulk action takes the following form [19] S = 1 2κ 2 where d 4 x g[r+6/l 2 F µν F µν +λ 2 (L 1 +L 2 +L 12 )], L 12 = k 2 M(1)µν M (2) µν, L (a) = 1 4 µ M (a)ντ µ M ντ (a) 1 4 m2 M (a)µν M µν (a) 1 2 M(a)µν F µν 1 8 JV(M(a) µν ), V(M (a) µν ) = M(a)µ νm (a)ν τm (a)τ σm (a)σ µ, a = 1,2. L is the radius of AdS space, 2κ 2 = 16πG with G the Newtonian gravitational constant, k, m 2 and J are all model parameters with J <, λ 2 characterizes the back reaction of the two polarization fields M µν (a) to the background geometry, and L 12 describes the interaction between two polarization fields. The equations of motion for polarization fields read 2 M (a) µν m2 M (a) µν km(b) µν JM(a) δ (a) τ µ M δ M (a) (1) (2) τν F µν =. (3) Here (a,b) = (1,2) or (2,1). In the probe limit of λ, we can neglect the back reaction of the two polarization fields on the background geometry. The background we

2 2 will consider is a dyonic Reissner-Nordström-AdS black brane solution of the Einstein-Maxwell theory with a negative cosmological constant, and the metric reads [2] ds 2 = r 2 ( f(r)dt 2 +dx 2 +dy 2 )+ dr2 r 2 f(r), f(r) = 1 1+µ2 +B 2 r 3 + µ2 +B 2 r 4. (4) Here both the black brane horizon r h and AdS radius L have been set to be unitary. The temperature of the black brane is T = (3 µ 2 B 2 )/4π. (5) For the solution (4), the corresponding gauge potential is A µ = µ(1 1/r)dt + Bxdy. Here µ is the chemical potential and B can be viewed as the external magnetic field of dual boundary field theory. The zero temperature limit then corresponds to B 2 +µ 2 = 3. Critical magnetic field and Néel temperature. Let us first consider the influence of the external magnetic field B on the antiferromagnetic critical temperature T N. We consider a self-consistent ansatz with nonvanishing M (a) tr,m xy (a) (a=1,2) and define α = (M (1) xy +M (2) xy )/2, β = (M (1) xy M (2) xy )/2. (6) It is found that the equations for M (a) tr decouple with those of α and β and therefore we pay main attention on αandβ. ThebehaviorofthesolutionsofequationsinUV regiondepends on the value ofm 2 +k+4. When m 2 +k+ 4 =, the asymptotic solutions will have a logarithmic term, we will not consider this case in the present paper. When m 2 +k +4, we have the asymptotic solution at the AdS boundary: α UV = α + r (1+δ1)/2 +α r (1 δ1)/2 β UV = β + r (1+δ2)/2 +β r (1 δ2)/2, B m 2 +k +4, δ 1 = 17+4k +4m 2, δ 2 = 17 4k +4m 2, (7) where α ± and β ± are all finite constants. The antiferromagnetic condensation phase requires the parameters m 2 and k satisfy the following conditions [19] 3/2+m 2 < k < 4+m 2. (8) and α + = β + = in the spirit of the AdS/CFT correspondence. Near the critical temperature, the staggered magnetization is very small, i.e., β is a small quantity. In that case we can neglect the nonlinear terms of β and obtain α + f α f β + f β f ( + Jα3 2f r 6 f rf + 4 r 2 m2 +k r 2 f ( 3Jα 2 + r 6 f 2f rf 4 r 2 m2 k r 2 f ) α = B r 2 f, ) β =. (9) Near the IR region, the regular initial conditions are imposed. Without loss of generality, we can set β(r h ) = 1. As pointed out in Ref. [19], with increasing the magnetic field B from zero, the value of α and the effective mass of β increase, so that the critical temperature T N decreases. When T N is near zero, the regular initial conditions for (9) imply that B c = Jα 3 α (m 2 +k), (1) with α = ± (m 2 k)/3j at the horizon. Here α is the initial value of α at the IR fixed point and B c is the critical magnetic field where the QPT happens. Since we need that the value of α monotonically increases with the increasing of B in order that the Néel temperature monotonically decreases, we need an additional restriction, k m 2. (11) The result (1) can also be understood in the near horizon geometry [21]. In the zero temperature limit, an AdS 2 geometry emerges in the IR region. The solution for α in this region has the following form α(r) = c + r (1+δ )/2 +c r (1 δ )/2 + Jα 2 B (12) k m2, with δ = 1 2(Jα 2 k m2 )/3, r = 1/(r 1) and r 1 +. The existence of the IR fixed point requires δ > 1 and c + =. As a result, we have α = lim r 1 α(r) = B/(Jα 2 k m2 ). The asymptotic solution of β in the IR region takes the form β = d + r (1+δ )/2 +d r (1 δ )/2, (13) with δ = 9 6(3Jα 2 +k m2 )/3. The existence of non-trivial IR fixed point needs δ = 1 and d + =, thus the condition 3Jα 2 +k m2 = has to be satisfied. Spectrum and energy gap near the QCP In order to investigate the magnetic fluctuations in the vicinity of QCP, we need to consider the perturbations of two polarization fields. Because of the non-linear term in equation (3), the perturbations for all the components have to be included δm (a) µν = ǫc(a) µν e i(ωt+qx), (µ,ν) (r,y),(t,x) δm (a) µν = iǫc(a) µν e i(ωt+qx), (µ,ν) = (r,y),(t,x). (14) Put this perturbations into the equation (3) and compute to the 1st order for ǫ, we can get their equations of the perturbations. In general, all the components couple with each others. However, in the limit of q, we find the the perturbations for the xy components decouple with others [21]. Let β = (C xy (1) C xy (2) )/2 and α = (C xy (1) + C xy (2) )/2. Near the horizon, we impose the

3 3 ingoing conditions for α and β. Near the AdS boundary, we have the following asymptotic solutions for α and β α α + r (1+δ1)/2 + α r (1 δ1)/2, β β + r (1+δ2)/2 + β r (1 δ2)/2. (15) According to the dictionary of AdS/CFT, the terms α + and β + are the sources for the dual operators and the terms α and β are the vacuum expectation values in the boundary theory. By the definition, the retarded Green s functions for α and β read 1 ] [ ] [ α G α α,g α β ][ α+ =. (16) β β + G β α,g β β We can see that the magnetic quasi-particle excitations are the mixture of two polarization fields. In the paramagnetic phase or the antiferromagnetic phase without magnetic field, i.e. β = or α =, these two quasiparticle excitations decouple into paramagnetic and antiferromagnetic excitations [21]. In these two cases, we can call the excitations α and β as para-magnon and antiferro-magnon respectively. The Green s functions correspond to them then are, G αα = α / α +, G ββ = β / β +. (17) Using retarded Green s function, we can define spectrum function as P(ω, q ) = ImG(ω, q )/π. When we turn on a small momentum q, the peak of P(ω, q ) for longlife quasi-particle can be given by following dispersion relation, ω = +ǫ q, ǫ q = =. (18) Here is the energy gap of quasi-particle excitation. In the vicinity of QCP, for the case of ω =, the retarded Green s function usually has the form of G 1/(q 2 + 1/l 2 ), where l is called correlation length. At QCP, in general, the energy gap vanishes. Thus we have ω = ǫ q. In addition, for small frequency and wave vector, we can define the dynamic exponent z in the way as ω q z. As the equations involved here are nonlinear, we have to solve them numerically. The different parameters satisfying restrictions (8) and (11) give similar results, we here therefore just show typical examples by taking k = m 2 = 3/2, J = 6 in Fig. 1 and J = 1 in Figs. 2, 3 and 4. ItisinterestingtocompareourresultsontheNéeltemperature T N versusmagnetic field B (B < B c ) with some real materials. In Fig. 1 we present the experimental data for BiCoPO 5 [22] by red circles. The green dashed 1 Here we drop off some constants associated with the conformal dimension. TN AF Holographic Empirical fit Data PM B/Bc FIG. 1. The antiferromagnetic critical temperature T N versus the external magnetic field B calculated at the model parameters k = m 2 = 3/2, J = 6 (solid curve) and compared with experimental data from BiCoPO 5 (red circles) and mean field theory (green dashed line). The experimental data are from [22] and rescaled. TN AF Holographic x= x=.31 x=.85 PM B/Bc FIG. 2. The antiferromagnetic critical temperature T N versus the external magnetic field B calculated at the model parameters k = m 2 = 3/2, J = 1 (solid curve) and compared with experimental data for pyrochlore compounds: Er 2 2xY 2xTi 2O 7. Black circles corresponds to zero doping, x = ; magenta triangulars - to x =.31; green crosses - to x =.85. The experimental data are from [23]. curve is the fitting proposed by [22]. Here we rescale the original data. The authors in [22] argued that these data can be described in terms ofamean-field power law. However, compared with our model curve, they does not satisfy conventional mean field solution but can be fitted well by our holographic model. In the Fig.2 we compare our results (the solid curve) with experimental data taken from [23] for other class of materials - the pyrochlore compounds, Er 2 2x Y 2x Ti 2 O 7. Dataarepresentedbyblackcirclesforzerodoping,x = ; by magenta triangulars - for x =.31; by green crosses - for x =.85. This class of compounds is very different from the BiCoPO 5 and the low temperature antiferromagnetic state is residing there on the pyrochlore lattice [23]. Again as in previous case the model describes well the Nel temperature lowering with magnetic field (see, Fig.2) and predicts here the Quantum Critical Point at 2.5 Tesla. There magnetic spins become partially aligned in the direction of the magnetic field. Therefore the system requires less thermal energy to destroy the remaining magnetic spins order. It is also interesting here that the

4 ImG B=1.1Bc B=1.5Bc B=1.11Bc G l ln(ω/tn) B/Bc q 2 /Tc B/Bc 1 FIG. 3. Left: The antiferromagnetic spectrum function in the case with different magnetic field when B > B c. Right: The gap energy versus the external magnetic field when B/B c FIG. 4. Left: G 1 as a function of q 2 when ω = in the case of B = 1.1B c. The solid line is the fitting curve by G 1 q 2 +1/l 2. Right: The correlation length l versus the magnetic field when B/B c critical temperature versus critical field curve is independent of the diamagnetic doping at the magnetic site, i.e. the influence of dilution in Er 2 2x Y 2x Ti 2 O 7 solid solutions is similar to the increase of magnetic field for pure Er 2 Ti 2 O 7. This indicates the universality of holographic model in ability to describe complex interactions in these pyrochlores which local magnetic structure changes with doping. Thus, the model shows that with increasing magnetic field from zero to B c, the Néel temperature decreases from T N to zero. For small B, numerical results show that T N T N B 2. When magnetic field is close to B c, we find that Néel temperature is fitted well by following relation T N /ln T N (1 B/B c ), (19) where T N = T N /T N. We will present the relation (19) by considering the geometry AdS 2 at IR limit [21]. When magnetic field B is larger than the critical value B c, the antiferromagnetic phase disappears even at the zero temperature. In this case, the system comes into quantum disordered phase at zero temperature, in which there is a gapped magnetic excitation. In the left plot of Fig. 3, we show ImG with respect to the frequency of anferromagnetic excitation in the case with different magnetic field. Since we are mainly interested in the antiferromagnetic excitation in the paramagnetic phase, here we just show the results of G ββ. In the case of < B/B c 1 1, there is a distinct peak which gives the energy gap for the excitation. With increasing magnetic field, the peak moves towards higher energy and becomes more and more indistinct. This means that the gap increases but the lifetime decreases when magnetic field increases. At the critical magnetic field B = B c, we see ω =, which corresponds to a gapless longlifetime antiferromagnetic excitation. In the region of B/B c 1 +, we find the energy gap is fitted well by following equation (see the right plot of Fig. 3) (B/B c 1),with = /T N. (2) In the left plot of Fig. 4, we plot the inverse of Green s function G 1 ββ (q) in the case of ω = and B = 1.1B c. We can see that it obeys the behavior of G 1 q 2 +1/l 2 as we expected before. Thus the Green s function can give the correlation length by fitting the curve of G 1 as a function of q 2. The correlation length as a function of B is plotted in the right plot of Fig. 4. Results show that the correlation length l as a function of the tuning parameter B obeys the following relation l (B/B c 1) ν, with ν 1/2. (21) As to the dynamical exponent z, we can obtain by using the similar numeric method. In antiferromagnetic metal, z = 2 [5]. In our model, the results indeed show that z 2. This numerical result can be confirmed by the emergent AdS 2 geometry in the IR region [21]. The same result from the holographic frame but in a different model can be found in Ref. [24]. Furthermore from the energy gap (2), we see that this energy gap satisfies the universal scaling relation B B c zν. Discussion The relation we found in this paper of the Néel temperature with respect to the magnetic field in the vicinity of QCP is quite non-trivial and can fit well experimental results[22, 23]. The relation (19) is not the usual power-law behaviors or square-root form in the mean field theory. But it is the expected result in the 2-D QPTs in strong coupling case [5]. This non-trivial coincidence strongly implies a connection between these two different theories. In the antiferromagnetic metal where magnetic ordered is dominated by itinerant electrons, dynamical exponent z is 2 near the QCP. Since our dual boundary theory is a 2-D theory, the effective dimension is thus d eff = d+z = 4, which is just the upper critical dimension of the Hertz field theory [2, 25]. In this case, the hyperscaling is violated and logarithmic correction behavior appears. In fact, the d = z = 2 quantum critical theory is in general not a weak coupling theory at any T >. Instead, a strongly coupled effective classical model emerges that can be used to determine the critical dynamics [26]. Our results show that it can be described well by AdS/CFT correspondence and this provides a new example of the applicability of the gravity/gauge duality in condensed matter theory.

5 5 Acknowledgements This work was supported in part by the National Natural Science Foundation of China (No , No.11358, No , and No ). [1] S.SachdevandB.Keimer, Phys.Today64(2), 29(211). [2] J. A. Hertz, Phys. Rev. B 14, (1976). [3] S. Sachdev, Quantum Phase Transitions, Cambridge University Press, Cambridge (1999). [4] P. Gegenwart, Q. Si, and F. Steglich, Nat. Phys. 4, (28). [5] H. v. Lohneysen, et al., Rev. Mod. Phys. 79, 115 (27). [6] G. Chaboussant, et al., Phys. Rev. B 55, 346 (1997). [7] W. Shiramura, et al., J. Phys. Soc. Jpn. 66, 19 (1997); A. Oosawa, et al., Phys. Rev. B 66, 1445 (22). [8] Ch. Rüegg, et al., Nature 423, 62 (23); T. Nikuni, et al., Phys. Rev. Lett (2). [9] Q. Si, S. Rabello, et al., Nature 413, (21). [1] T. Senthil, et al.,science 33, (24); H. Tanaka, et al., J. Phys. Soc. Jpn. 7, 939 (21). [11] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999)] [12] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428, 15 (1998) [13] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998) [14] S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, Phys. Rev. Lett. 11, 3161 (28) [15] S. S. Lee, Phys. Rev. D 79, 866 (29) [16] H. Liu, J. McGreevy and D. Vegh, Phys. Rev. D 83, 6529 (211) [17] M. Cubrovic, J. Zaanen and K. Schalm, Science 325, 439 (29) [18] R. -G. Cai and R. -Q. Yang, Phys. Rev. D 9, 8191 (214) [19] R. -G. Cai and R. -Q. Yang, [arxiv: [hep-th]]. [2] R. -G. Cai and Y. -Z. Zhang, Phys. Rev. D 54, 4891 (1996) [21] R. -G. Cai, R. -Q. Yang, and F.V. Kusmartsev, in preparation. [22] E. Mathews et al. Solid State Communications 154(213) [23] J.F. Niven et al. Proc. R. Soc. A 47: (214). [24] N. Iqbal, H. Liu, M. Mezei and Q. Si, Phys. Rev. D 82, 452 (21) [25] Ar. Abanov, A. Chubukov, Phys. Rev. Lett (24); [26] S. Sachdev, and E. R. Dunkel, Phys. Rev. B 73, (26)