Author s version 1997 arxiv:cond-mat/ v1 [cond-mat.stat-mech] 3 Apr 1998

Size: px

Start display at page:

Download "Author s version 1997 arxiv:cond-mat/ v1 [cond-mat.stat-mech] 3 Apr 1998"

Stewart Robertson
6 years ago
Views:

1 Z. Phys. B, 8 87 (997) ZEITSCHRIFT FÜR PHYSIK B Author s version 997 arxiv:cond-mat/98444v [cond-mat.stat-mech] Apr 998 Susceptibilities and correlation functions of the anisotropic spherical model D. A. Garanin I. Institut für Theoretische Physik, Universität Hamburg, Jungiusstr. 9, D-55 Hamburg, Germany ( garanin@physnet.uni-hamburg.de) Received: 8 August 996 Abstract. The static transverse and longitudinal correlation functions (CF) of a -dimensional ferromagnet are calculated for the exactly solvable anisotropic spherical model (ASM) determined as the limit D of the classical D-component vector model. The results are nonequivalent to those for the standard spherical model of Berlin and Kac even in the isotropic case. Whereas the transverse CF has the usual Ornstein- Zernike form for small wave vectors, the longitudinal CF shows a nontrivial behavior in the ordered region caused by spin-wave fluctuations. In particular, in the isotropic case below T c one has S zz (k) /k (the result of the spin-wave theory) for k< κ m T c T. PACS: 75..Hk; 75.4.Cx The spherical model (SM) of Berlin and Kac [] was believed many years to be the only model in the statistical theory of magnetism, which is exactly solvable in dimensions. The other model, which is the limit D of the classical D-component vector model introduced by Stanley [], has the same partition function in the homogeneous case [], and for this reason the latter was not considered as an independent model by many researchers and received relatively little attention. However, the D model is nonequivalent to the standard spherical model of Berlin and Kac and goes beyond it in many respects. First, it can be easily generalized for anisotropic systems [4]. Second, it does not use the global spin constraint, which leads to unphysical results in spatially inhomogeneous situations. In particular, in the SM the Curie temperature T c of a 4- dimensional ferromagnetic film with free boundary conditions, which is infinite in dimensions and finite in the 4th dimension, has been found by Barber and Fisher [5] tobe anon-monotonousfunctionofthenumberofitslayers N. To the contrast, an improved versionof the spherical model using separate spin constraints in each layer[6] leads in the case mentioned above to the monotonically increasing T c (N), as it should be. However, even this improved version of the SM fails on physically meaningful -dimensional ferromagnetic films, since the latter are -dimensional objects and without anisotropy their T c is zero for any finite thickness. An adequate description of spatially inhomogeneous and low-dimensional ferromagnetic systems on the spherical level can be archieved only with the use of the anisotropicd model, which can be called also anisotropic spherical model (ASM). This model was recently applied to study the dimensional crossover of T c (N) of anisotropic ferromagnetic films [7] and the effects of thermal fluctuations in the Bloch-wall phase transition [8]. Here it will be shown that the ASM deviates from the standard SM even in the spatially homogeneous case and even in the isotropic limit, if the spin-spin correlation functions are concerned. Indeed, below T c in the D model there are two different longitudinal and transverse correlation functions, whereas there is only one CF in the SM. A less trivial reason for the difference between two models relies on the fact that a wave-vectordependent CF is proportional to the appropriate susceptibility, which is the linear response to a spatially inho-

2 mogeneous sinusoidal field. In this case, as was argued above, the global spin constraint used by the SM modifies the results. As we shall see, the k-dependent longitudinal CF of the ASM has a nontrivial non-ornstein-zernike form below T c, which is similar to that following from the lowest-order spin-wave theory for T T c [9,,]. For the uniaxially anisotropic ferromagnetic model the ASM yields finite longitudinal and transverse susceptibilities and correlation lengths in the ordered region. The anisotropic generalization of the classical D- component vector model of Stanley [] can be described by the Hamiltonian H H i m i D J ij η α m αi m αj, () ij α where m i and η α are anisotropy coefficients. Here we consider the uniaxial model with η η z and η α η for α. The model () can be conveniently treated by the classical spin diagram technique [,,4], which allows classification of diagrams in powersof/d ford []. The equationofstate m(h,t) of the anisotropic spherical model, where m m z, is contained in the diagram series corresponding to the self-consistent Gaussian approximation (SCGA) [,4], which becomes exact in the limit D. In terms of the dimensionless variables θ T/Tc MFA and h H/J, where Tc MFA J /D is the mean-field transition temperature and J is the zero Fourier component of the exchange interaction, one comes to the system of equations for the magnetization m and the gap parameter G [,4] Here G m mh, θgp(ηg) m. () P(X) v (π) Xλ q, () where v is the unit cell volume and λ q J q /J. In the long-wavelength limit λ q αq, where α a and a is the lattice spacing. For -dimensional lattices with the nearest neighbour interactions the integral P(X) has the following properties: { P(X) X /z, X W c ( X) / (4), X, where z is the number of nearest neighbors, W is the Watson integral, and c v /(4πα / ). For the simple cubic (sc) lattice v a and α a /6, hence c (/π)(/) /. The solution of the system of equations () simplifies for zerofield, where below T c one has G (the zero spin-wave gap) and m ( θ/θ c ) /, θ θ c /P(η). (5) It can be seen that in the isotropic case, η, the value of the phase transition temperature θ c reduces to the well-known result θ c /W []. Using () one can calculate the longitudinal susceptibility χ z m/ H. The zero-field reduced susceptibility, χ z J χ z, has the form G χ z G, θ > θ c (6) θ m [ηp (η)p(η)], θ < θ c, where P (X) dp(x)/dx and G satisfies the equation θgp(ηg) above θ c. Solving this equation near θ c in the linear approximation in G, using (5), and introducing the reduced temperature variable ǫ θ c θ one can rewrite (6) in the form where I(η) ηp (η) P(η) χ z I(η) { ( ǫ), θ > θ c (ǫ), θ < θ c,, (7) η /z, η c, η. P(η) η (8) (9) The first line of (8) is valid for the weakly-anisotropic model, η, in the narrow temperature interval ǫ ǫ c η, () 4W whereas the second one is valid in the whole region below θ c. The latter is the reason to define ǫ in the nonstandard form (7). It can be seen from (8) that the critical behaviour of the longitudinal susceptibility in the ASM is the same as that in the mean field approximation (MFA), including the famous ratio of critical amplitudes. This equivalence holds for all the critical indices [], since for η < the square-root singularity of P(X) at X is suppressed. In the extreme case of the spherical Ising model, η, the mean field approximation for the ASM becomes exact, since the fluctuations of the transverse spin components die out with the transverse coupling η and the influence of longitudinal fluctuations vanishes in the spherical limit, D. For the anisotropic model, η <, the quantity χ z is finite in the ordered region, vanishing at θ and diverging at θ θ c. On the contrary, the longitudinal susceptibility of the isotropic model, η, diverges as χ z ( ǫ) aboveθ c andisinfiniteinthewholeregionbelowθ c dueto the Goldstone-mode fluctuations. In this case for h one has

3 5 Susceptibilities of the uniaxial spherical model: η.99 z z z z z z χ... 5 χ z χ z Tc T/T MFA c FIG.. Temperature dependences of the longitudinal and transverse susceptibilities of the uniaxial spherical model. m m(h) m() θc m /h/, () which coincides with the result of the standard spherical model and is in accord with the prediction of the lowest-order spin-wave theory [5] for θ θ c. A similar expression for a general D-component vector model, which contains the additional factor /D, was derived by Fisher and Privman [6] from the scaling arguments. Measuring the square-root singularity of magnetization () in systems with spontaneous breaking of a continuous symmetry is a difficult task, and it was done only resently on EuS [7] and on EuO [8]. In the following we shall study the correlation functions S αα (r) m α () m α (r), α,,...,d, () which are translationally invariant in the spatially homogeneous case considered here. The fluctuation-dissipation theorem relates wave-vector-dependent CFs to the appropriate susceptibilities as χ α (k) βs αα (k), where β /T. Both longitudinal and transverse CFs enter the self-consistent Gaussian approximation and were(approximately) calculated in [,4]. Since the SCGA becomes exact in the limit D, the transverse CFs describing in the SCGA fluctuations of D transverse spin components become exact, too, and with the use of the results of [4] one obtains χ (k) D θ S (k) G ηgλ k () in the whole temperature range [cf. (6)]. One can see that for h below θ c, where G, the homogeneous transverse susceptibility, χ χ (), diverges in the isotropic limit, η, and remains a finite constant for the uniaxially-anisotropic model. Comparing (8) and (), one finds that both susceptibilities become equal to each other at the characteristic temperature slightly below θ c, which is determined by ǫ ǫ [see (7) and FIG.. (a) The compact part ˆΛ zz(k) of the longitudinal CF S zz(k); (b) block summation for one-site spin cumulants Λ... renormalized by Gaussian fluctuations of the molecular field; (c) ladder equation for the transverse four-spin correlation line Ṽk; (d) Dyson equation for the renormalized transverse interaction ηβ J q. Unlabeled small circles denote transverse spin components. ()]. The temperature dependences of the longitudinal and transverse susceptibilities for the weakly-anisotropic spherical model are represented in Fig.. One can see that even a small anisotropy, as η, has a profound influence on χ z below θ c. Unlike the transverse CF given by (), the exact longitudinal CF, S zz (k), cannot be determined from the SCGA in the spherical limit, since for D the fluctuations of the single longitudinal spin component become nonessentialands zz (k) disappearsfromthescgaequations. Therefore, the wave-vector-dependent longitudinal CF in the anisotropic spherical model should be considered separately, which is the main purpose of this paper. With the help of the classical spin diagram technique [,,4] S zz (k) determined by () can be represented as S zz (k) ˆΛ zz (k) ˆΛ zz (k)βj k, (4) where ˆΛ zz (k) is the compact (irreducible) part of S zz (k) given by the diagrams, which cannot be cut by the one longitudinal interaction line βj k. The quantity ˆΛ zz (k) is in turn given in the limit D by the set of diagrams represented in Fig.. Such a choice of diagrams is based on the arguments of []. More technical details can be found in [9], where the transverse CF, S αα (k), α, was calculated up to the first order in /D for low-dimensional ferro- and antiferromagnets in magnetic field. The analytical form of ˆΛ zz (k) of Fig. a reads ˆΛ zz (k) Λ zz Λ ααzṽk. (5) Here in the limit D one has [9]

4 Λ zz Λ αα Λ ααββ ξ, Λααz Λ ααββ ξ, (6) where α β z, Λ αα θg D, Λααββ G satisfies the system of equations (), and ξ β(h mj ) D θ ( ) θg D θg/, (7) (hm) Dm θg (8) is the temperature-normalized molecular field. The quantity Ṽk in (5) is the solution of the ladder equation Fig. c and has the form Ṽ k where Λ ααββ is given by (7), V k D V k Λ ααββ V k, (9) v (π) ηβ J q ηβ J k q, () and the renormalized transverse interaction ηβ J q determined by Fig. d reads ηβ J q ηβj q Λ αα ηβj q. () The factor D in () results from the summation over transverse spin components in Fig. a or in Eq. (5), as well as in the ladder equation Fig. c. Such a factor does not appear, if one tries to take into account the similar diagrams for the transverse CF S (k) S αα (k), α, and thus these diagrams vanish in the spherical limit and S (k) has the trivial Ornstein-Zernike form (). Combining now Eqs. (4)-(), one comes after simplifications to the final result where χ z (k) G λ k m θg r k, () r k v (π) Gηλ q Gηλ k q. () The result similar to (), in which, however, only the last term is present, was obtained earlier [9,,] within the lowest-order spin-wave theory well below T c. On the other hand, for the standard spherical model the correlation function has the trivial Ornstein-Zernike for all temperatures [,]. The integral in () can be easily calculated for k and in the corner of the Brillouin zone, k b, for the sc lattice, where λ b q λ q. The result has the form { I(ηG), k r k P(ηG) (4), k b, where I(X) is given by (9). One can see that χ z () of () is in accord with the previously obtained expression (6) for h below θ c, where G. Under the same conditions χ z (b) θ/(θ c ), which can be used to control numerical calculations. For the weakly anisotropic model in the case ηg the integral () can be calculated analytically for small wave vectors, a k, which results in ( r k c arctan αk αk ηg ). (5) The last term in () modifies the k-dependence of χ z (k) below θ c in the gapless case ηg in the range of small wave vectors, which shrinks if θ c is approached from below. Instead of the Ornstein-Zernike form χ z (k) /k one has (cf. [9,,]) χ z (k) k, k κ m m θ c α/. (6) The latter defines a length scale ξ m, which exists only in the ordered region (m > ) and diverges at θ c due to the vanishing of magnetization (5): ξ m κ m θ m c α /. (7) For the sc lattice c α / π a and hence ξ m θ 4πm a. The length ξ m is analogous to the bare, i.e., the meanfield correlationlength below T c, which follows, in particular, from the Landau-Ginzburg phenomenological free energy. This analogy is, however, not complete, since ξ m diverges at the actual transition temperature T c and not at Tc MFA (θ ). The crossover of χ z (k) at k κ m for temperatures slightly below T c wasdescribed earlierwith the help of the renormalization group approach []. The k-dependences of χ z (k) of () in the isotropic case obtained with the help of numerical integration in () are represented in Fig. for different temperatures below T c. Curves of such a type for Heisenberg model systems such as EuO and EuS could be observed, in principle, in neutron scattering experiments, but such experiments were not carried out up to now. Our next task is to calculate explicitly the real-space correlation functions () in the small-anisotropy case η. Since the CFs themselves are proportional to /D and thus vanish in the spherical limit [see, e.g., ()], it is more convenient to deal with the appropriate susceptibilities, which are given by χ αα (r) v (π) eikr χ αα (k). (8) In the transverse case using () one comes to the wellknown result χ (r) c α / e r/ξ c, r ξ c α ηg (9) 4

5 4 χ z (k). T c.4 T c.6 T c 4 5 Correlation lengths of the uniaxial spherical model: η.99.8 T c T c 5 ξ c /α / ξ cz /α / ξ m /α / k/k max FIG.. The wave vector dependences of the inverse longitudinal susceptibility χ z (k) of the isotropic spherical ferromagnet along the [] direction of k for different temperatures in the ordered region. for r a. One can see that in the isotropic case, η, the transverse correlation length, ξ c, if infinite in the wholeregionbelowθ c,whereg. Inmorecomplicated situations, as for the longitudinal susceptibility below θ c, (), the correlation length can be determined as ξ cz /κ c, where κ c is the singularity point of χ z (k) on the imaginaryaxisk iκ, whichis closesttothe origin. This leads with the use of (5) to the transcedental equation for u c ξ c /(ξ cz ) having the form u c ln u c u c a, a ξ c ξ m m θc η, () where ξ m is given by (7). This equation coincides with Eq. (4.9) of [8], which determines the temperaturedependent width δ L ξ cz of the (linear) domain wall in the uniaxial spherical model. The asymptotic solutions of () read { u c a/, a (slightly below θc ) e a (), a (far below θ c ). Note that the transition between two regimes in () occurs at ǫ ǫ η [see ()]. The explicit results for the correlation length ξ cz itself, including that above θ c, read α G αi(η) ǫ, < ǫ ǫ ξ cz ξm ξ c αi(η), < ǫ ǫ () ǫ [ ( ξ c exp ξ )] c, ǫ ǫ, ξ m where I(η) is given by the second line of (9). The numerically calculated temperature dependences of the longitudinal and transverse correlation lengths are illustrated in Fig. 4. T c T/T MFA c FIG. 4. Temperature dependences of the transverse, ξ c, and longitudinal, ξ cz, correlation lengths, as well as the characteristic length ξ m, in the uniaxial spherical model. To calculate the Fourier-transform χ z (r) below θ c with the help of (8), (), and (5), it is convenient to perform at first the integration over the angle variables and then to deform the k-integration contour in the complex plane. In this way one comes to the well-behaved expression [ χ z (r) ac α / u e r/ξcz c r a( u c )u c ] du e ru/ξ c ( uln u u a ) π u, () where the first term is the contribution of the pole determining the longitudinal correlation length ξ cz [see ()] and the second one is the integral along the cut of the arctan-function of (5). The expression above simplifies in various limiting cases. For a (ǫ ǫ ), the cut term of () is small as a in comparison to the pole term and, additionally, it decays as exp( r/ξ c ), i.e., much faster than the first term of (). Thus, in this case only the pole term in () is relevant, and with the use of () one obtains χ z (r) c α / e r/ξcz, (4) r where ξ cz is given by the the middle line of (). In the opposite case, a (ǫ ǫ ), the pole sticks to the beginning of the cut, u, and the amplitude of the pole term becomes exponentially small because of the factor u c [see ()]. Thus, in this case only the cut term in () is essential, and the result can be simplified to χ z (r) c α / e r/ξcz πr dx x exp ( ) r x, (5) πξ m 5

6 where ξ cz is given by the the lower line of (). The asymptotic forms of (5) are (4) for r ξ m and ( ) χ z (r) m ξm e r/ξcz (6) θ r for r ξ m. One can say that for the isotropic model, where the true correlation length ξ cz is infinite below θ c, the bare one, ξ m of (7), still plays its role in some sence: For r> ξ m the slow decay χ z (r) /r changes to a faster one, χ z (r) /r. As a conclusion, the anisotropic spherical model (ASM) considered here is a good exactly solvable toy model for classical spin systems, which can be successfully applied in many situations where the standard spherical model of Berlin and Kac fails. More realistic models of phase transitions posess, naturally, the other, nonspherical, values of the critical indices, but this is, however, only a quantitative effect, which is less important in comparison to the profound role played by the gapless spin waves in the isotropic model below T c. The latter are properly taken into account in the ASM, which is thereby a very important step beyond the mean field approximation. The author thanks D. Görlitz, J. Kötzler, and Hartwig Schmidt for valuable discussions. The financial support of Deutsche Forschungsgemeinschaft under contract Schm 98/5- is greatfully acknowledged. [9] D. A. Garanin, J. Stat. Phys. 8, 97 (996). [] G. S. Joyce, in Phase Transitions and Critical Phenomena, edited by C. Domb and M. S. Green (Academic Press, New York, 97), Vol.. [] S. Ma, Modern theory of critical phenomena (Benjamin, Reading, 976). [] T. N. Berlin and M. Kac, Phys. Rev. 86, 8 (95). [] H. E. Stanley, Phys. Rev. Lett., 589 (968). [] H. E. Stanley, Phys. Rev. 76, 78 (968). [4] H. Okamoto, Phys. Lett. A, 5 (97). [5] M. N. Barber and M. E. Fisher, Ann. Phys. 77, (97). [6] G. Costache, D. Mazilu, and D. Mihalache, J. Phys. C 9, L5 (976). [7] D. A. Garanin, J. Phys. A 9, L57 (996). [8] D. A. Garanin, J. Phys. A 9, 49 (996). [9] H. Mori and K. Kawasaki, Prog. Theor. Phys. 7, 59 (96). [] A. Z. Patashinskii and V. L. Pokrovskii, Zh. Eksp. Teor. Fiz (97) [Sov. Phys. JETP 7, 7 (97)]. [] S. Lovesey and K. Trohidou, J. Phys.: Condensed Matter, 87 (99). [] D. A. Garanin and V. S. Lutovinov, Solid State Commun. 5, 9 (984). [] D. A. Garanin, J. Stat. Phys. 74, 75 (994). [4] D. A. Garanin, Phys. Rev. B 5, 59 (996). [5] T. Holstein and H. Primakoff, Phys. Rev. 58, 98 (94). [6] M. E. Fisher and V. Privman, Phys. Rev. B, 447 (985). [7] J. Kötzler, D. Görlitz, R. Dombrowski, and M. Pieper, Z. Phys. B 94, 9 (994). [8] A. Flosdorff, D. Görlitz, and J. Kötzler, J. Appl. Phys. 79, 464 (996). 6

Phase transitions and critical phenomena

Phase transitions and critical phenomena Classification of phase transitions. Discontinous (st order) transitions Summary week -5 st derivatives of thermodynamic potentials jump discontinously, e.g. (