arxiv: v1 [cond-mat.stat-mech] 19 May 2016

Transcription

1 The magnetic susceptibility on the transverse antiferromagnetic Ising model: Analysis of the reentrant behaviour arxiv: v1 [cond-mat.stat-mech] 19 May 2016 Minos A. eto, J. Ricardo de Sousa, Igor T. Padilha, Octavio R. Salmon, J. Roberto Viana Departamento de Física, Universidade Federal do Amazonas, 3000, Japiim, , Manaus-AM, Brazil F. Dinóla eto Cavendish Laboratory, University of Cambridge; J. J. Thomson Avenue, Cambridge, CB3 0HE, Engl. (Dated: May 20, 2016) Abstract ABSTRACT We study the three-dimensional antiferromagnetic Ising model in both uniform longitudinal (H) transverse (Ω) magnetic fields by using the effective-field theory with finite cluster = 1 spin (EFT-1). We analyzed the behavior of the magnetic susceptibility to investigate the reentrant phenomena we have seen the same phase diagram previously obtained in another papers. Our results shows the presence of two divergences in the susceptibility that indicates the existence of a reentrant behaviour. PACS numbers: Ak; Fr; Rh PACS numbers: Electronic address: minos@pq.cnpq.br 1

2 I. ITRODUCTIO The Ising antiferromagnet (AF) in transverse longitudinal magnetic fields divided into two equivalent interpenetrating sublattices A B, which is described by the following Hamiltonian H = J <i,j> σ z iσ z j H i σ z i Ω i σ x i, (1) where J is a positive coupling constant, which sts for the antiferromagnetic interaction between the magnetic moments (electrical origin), i, j represents the sum over all pairs of nearest-neighbour spins on a cubic lattice, σ ν i is the ν(= x,z) component, at site i, of the spin-1/2 Pauli operator, Ω H st for the transverse longitudinal magnetic fields, respectively. The operators σ x i σ z i do not commute, thus the field (Ω), which is transverse to the spin direction causes quantum spin fluctuations quantum tunnelling between the spin-up spin-down eigenstates of σi z. The fluctuations reduce the critical temperature T c below which the spins maintain a long-range order. The critical temperature vanishes at a critical field Ω c, a quantum phase transition occurs between the antiferromagnetic paramagnetic (P) quantum state. The magnetic susceptibility measures the capacity of a material to magnetise under the action of an external magnetic field. This property has been studied extensively in recent years. Some new results for the magnetic susceptibility for crystal fields Kondo effect were obtained by Desgranges [1]. Kashif, et al. [2] have studied the effects of copper addition on density magnetic susceptibility of lithium borate glasses: (100 x)(li 2 O 2B 2 O 3 ) xcuo, where x = 0, 5, 10, 15, mol%. At low temperatures, thedc magnetic susceptibility measurements for the heavy fermion CeCu 6 x Au x show an anomalous behaviour for the quantum critical concentration x = 0.1: The temperature dependence near to the quantum phase transition χ 1 θ+t α shows a coefficient α = 0.8 [3]. These results have been confirmed in some literatures [4, 5]. Recently, an intrinsic magnetic susceptibility of height purified single-wall carbon nanotubes[6]. The effect of fast neutron irradiation on the magnetic susceptibility magnetoresistance of Si whiskers with impurity concentration in metal-insulator transition they were observed by Druzhinin, et al. [7]. Inelastic quasi-elastic excitations in the neutron 2

3 scattering on the compounds CeCo 1 x Cu x Al 4 (x = 0, 0.05, 0.1) have been studied the antiferromagnetic ordering temperature, where the crystal field manifestation in inelastic neutron scattering, magnetic susceptibility specific heat can be observed on the compounds CeCoAl 4 [8]. Experimentally, another interesting analysis is the magnetism of linear chains antiferromagnetic at high temperatures by study of the dynamics susceptibility [9] of the hydrated compounds RbCoCl 3 2H 2 O. For these several years [10], magnetic phenomena with frustration or competing interations have been studied, because it is considered to be one of the most essential mechanisms for spin glass (SG) phenomena [11]. One of the characteristics that appear in this type of phenomena is the reentrant behaviour. This phenomenon can be seen in amorphous materials such as e.g. (Fe 1 x Mn x ) 75 P 16 B 6 Al 3, which displays the reentrance for several concentration x [12]. More recently, a reentrant spin glass behaviour in CE-type AFM Pr 0.5 Ca 0.5 MnO 3 manganite in the sequence of multiple magnetic transitions it was observed [13], showing a transition: P F AF SG. The prove of the existence of the SG phase, in this polycrystalline material, can be seen in behaviour of the susceptibility versus temperature. Another SG material that has the reentrant behaviour is the IrSr 2 Sm 1.15 Ce 0.85 Cu O 10 [14]. The reentrant behaviour, in this compounds, it can be seen via magnetic property measurements of AC DC susceptibility. The pure transverse Ising model (TIM) has been used to describe a variety of physical systems, for example, correlation functions for lattice gauge theories with action Boltzmann factor [15], thermodynamical properties of the mixed with four-spin interactions [16], critical behavior in two dimensional of the fidelity susceptibility [17], effect of surface dilution [18], description ferromagnetism in a transverse Ising antiferromagnet [19], reentrant phenomenona in nanosystems [20] effects of the romly in thin film [21]. Using EFT-1 with finite cluster =1 spin, in the present paper, we investigate the magnetic susceptibility analysis on the transverse of the Ising antiferromagnet in both external longitudinal transverse fields. This work is organized as follows: In Sec. II we outline the formalism its application to the transverse Ising antiferromagnetic longitudinal magnetic field; in Sec. III we discuss the results; finally, in Sec. IV we present our conclusions. 3

4 II. FORMALISM It was originally introduced by de Gennes[22] as a pseudospin model for hydrogen-bonded ferroelectricsuchaskh 2 PO 4 [23,24]. Theoretically, varioustechniqueshavebeenusedinthe transverse Ising model (TIM) as: renormalization group (RG) method [25], cluster variation method (CVM) [26], mean field theory (MFT) [27], pair approximation (PA) [28], Monte Carlo (MC) simulations [29] effective-field theory (EFT) [30 36] (we ll deal with here in this work). In what the ground-state of the model concerns (1), it is characterized by an antiparallel spin alignment in the horizontal vertical directions, thus it exhibits éel ordering within the initial sublattices A B. This is the AF state. This state must not be confused with the superantiferromagnetic (SAF) state, because the values of the longitudinal field H AF H SAF are different in the ground state. On the other h, the phase diagram exhibits interesting properties due to the competition between the antiferromagnetic exchange interaction the longitudinal field H. Particularly, the model (1) has an AF phase in the presence of a field, where the temperature of the transition decreases as the intensity of the fields increases, so at T = 0, a second-order phase transition occurs at critical values of H c Ω c. In order to treat the model (1) on a simple cubic lattice through the EFT-1 approach, we consider a simple example of cluster on a lattice consisting of a central spin surrounded by z spins. Then, these z spins are substituted by an effective field, which is proportional to the thermal average of the central spin. Accordingly, the Hamiltonian for this cluster is given by H 1A = H 1B = ( ( J J σ(1+δ)b )σ z H 1A z Ωσx 1A, (2) δ σ(1+δ)a )σ z H 1B z Ωσx 1B. (3) δ From the Hamiltonians (2) (3), we obtain the average magnetizations in sublattices A, m A = σ1a z, B, m B = σ1b z, using the approximate Callen-Suzuki [37], which are given by 4

5 H a 1A m A = (H a1a ) 2 +Ω tanhβ (H a 1A ) 2 +Ω 2, (4) 2 H a 1B m B = (H a1b ) 2 +Ω tanhβ (H a 1B ) 2 +Ω 2, (5) 2 where a 1A = J z σ(1+δ)b z a 1B = J z σ(1+δ)a z. δ δ ow, after using the identity exp(αd x )F(x) = F(x+a) (where D x = x is the differential operator) the van der Waerden identity for the two-state spin system (i. e., exp(aσ z i ) = cosh(a)+σi z sinh(a)) the Eqs. (4) (5) are rewritten as follows with F(x) = z m A = (α x +σ(1+δ)b z β x) F(x) x=0, (6) δ 0 z m B = (α x +σ(1+δ)aβ z x ) F(x) x=0, (7) δ 0 H x (H x)2 +Ω 2 tanhβ (H x) 2 +Ω 2, (8) where β x = sinh(jd x ) α x = cosh(jd x ). The Eqs. (6) (7) are rewritten in terms of multiple spin correlation functions. We remark that it is very difficult to address all the correlations between the spins in Eqs. (6) (7). So, in this work we resort to decouple the right-h sides of Eqs. (6) (7), in such a way that σ z ia σjb z...σz la ma m B...m A, (9) where i j... l m µ = σiµ z (µ = A,B). Although the approximation (9) (σ ) z 2 disregards correlations between different spins, it treats relations such as iµ = 1 exactly, whereas in the usual MFT all the self- multi-spin correlations are neglected. Then, by the use of the approximation (9), the Eqs. (6) (7) are rewritten as follows m A = m B = A p (T,h,δ)m p B, (10) p=0 A p (T,h,δ)m p A, (11) p=0 5

6 with A p (T,h,δ) = z! p!(z p)! αz p x β p x F(x) x=0, (12) where the coefficients A p (T,h,δ) are obtained by means of the mathematical property exp(αd x )F(x) = F(x+a). Close to the critical point, in terms of the uniform m = 1 2 (m A +m B ) staggered m s = 1 2 (m A m B ) magnetizations, we have m s 0 m m 0, thus, the expansion of the sublattice magnetization m A (up to linear order in m s ) is given by with m A = ω 0 (T,h,δ,m 0 )+ω 1 (T,h,δ,m 0 )m s, (13) ω 0 (T,h,δ,m 0 ) = ω 1 (T,h,δ,m 0 ) = A p (T,h,δ)m p 0, (14) p=0 p=0 pa p (T,h,δ)m p 1 0. (15) On the other h, due to the fact that we only observed second-order phase transitions, we analysed only the Eqs. (14) (15) in the limit of m s 0, so as to obtain the phase diagram. So we could locate the second-order line using the fact that m A = m 0 +m s in Eq. (13), thus: ω 0 (T,h,δ,m 0 ) = m 0 (16) ω 1 (T,h,δ,m 0 ) = 1, (17) at the critical point in which m s = 0, δ = Ω/J h = H/J. Therefore, we could find a solution for the second-order critical frontier, being m s the order parameter, which is suitable for describing the phase transition of the model (1). It is important to mention that the magnetizations of the two sublattices are not equal for m s 0, the system is in the antiferromagnetic phase. ervertheless, the magnetizations of the two sublattices are equal for m s = 0, so the system is in the saturated P phase. The initial magnetic susceptibility χ s is calculated by the derivative of the magnetization function m s (H) given by the subtraction of the equations (10) (11), i. e., 6

7 where Θ(m,m s ) (m,m s ). χ s = ( ) ms = Θ H H=0 1+, (18) III. RESULTS AD DISCUSSIOS In Figs. 1. (a-c), we show the behaviour of the magnetic susceptibility as a function of temperature by varying the values of the magnetic fields h δ. The numerical determination of the phase boundary (second-order phase transition) is obtained by solving simultaneously the set of Eqs. (16) (18). In Figs. (1.a) (1.b). We set the values of the magnetic fields δ (fixed h) h (fixed δ), respectively, where we can see clearly the existence of the reentrant phase transition (P AF P). It is a well-known example is the cuprate material, IrSr 2 Sm 1.15 Ce 0.85 Cu O 10, shows in measures of DC-susceptibility an evidence of complex set of magnetic transitions upon cooling that are characteristic of a reentrant SG ground-state. In Fig. 1.(a), for h = 6.20 δ = 0.7, the phase transition occurs at only one point (divergence of χ s ) in T (1) = However, inside the reentrant region when we increased the value of h, we find another transition peak. In the same figure, we find that for h = 6.30 δ = 0.7 we have T (1) T (1) = T(2) = T(2) = 3.082, for h = 6.40 δ = 0.7 we have = 3.082, respectively. In the limit of high temperature the magnetic susceptibility as a function of the temperature presents a behaviour kind 1/T. In Fig. 1. (b), the same behavior appears when we fix the value of h increase the value of δ. For case h = 6.20 δ = 0.2, the phase transition occurs at two one point in T (1) = T(2) occurs in the points T (1) = 3.303, for case h = 6.20 δ = 0.4, the phase transition = T(2) = We note that, fixing the value of h varying δ the values of the temperatures of transition decreases. A possible, degeneracy increase in the energy of the ground-state these regions may occur in magnetic models with competitive interactions are sensitive the these perturbation in low-temperature. The two transition points can also be seen in Fig. 1. (c), where we observe three behavior of the susceptibility as a function of temperature with some values of h δ. For case h = 6.05 δ = 0.01, the phase transition occurs at two one point in T (1) = T (2) = 3.473, the case h = 6.35 δ = 0.7, the phase transition occurs at two one 7

8 point in T (1) = 0.7 T(2) = the case h = 6.20 δ = 0.9, the phase transition occurs in the one point T (1) = For the occurrence of these reentrant phenomenon in a certain appropriate range of h δ, explanation of mechanisms has been given by eto de Sousa [38]. In this work the phase diagram in the h T plane is presented with selected values of δ that show the presence of reentrant behaviour in the region of low temperature, see Fig. (1.d). IV. COCLUSIOS These magnetic phenomena in which frustration competitive interactions are involved, have been studied due to the fact that they are considered one of the most essential mechanism for spins glass phenomena more recently in behaviours of multiple reentrant glass transitions in confined hard-sphere glasses [39]. This phenomenon appears in glass transition in a colloid polymer mixture with depletion attractions [40] showing that in the introduction in the short-range attraction to a colloid suspension of nearly hard spheres by addition of a free polymer produces this phenomenon. In both cases it appears the two transition points near the second transition point the derivative dχ s /dt > 0 is positive at high-temperature the entropy is the predominant factor the system is then in the disordered P phase but with an AF bias due the applied fields. This is a good model that could represent multiple magnetic transitions as it is seen experimentally in some materials [13, 14]. As expected, χ s increases rapidly with increasing temperature diverges at the critical point T c behavior is presented as 1/T in the limit of high-temperature, as expected. In this work we analyze the reentrant behavior on the magnetic susceptibility on the region of the reentrant phenomenon using the EFT-1. Our results are consistent with secondorder transitions from the AF ordered to the P disordered phase at zero transverse field. Furthermore, the investigations of this model with quenched disorder is expected to show many characteristic phenomena, as already analysed in models of metal-insulator transition [41]. This will be discussed in future work. ACKOWLEDGEMET 8

9 This work was partially supported by FAPEAM CPq (Brazilian Research Agency). [1] H. -U. Desgranges, Physica B, (2015). [2] I. Kashifet al., Physica B, (2008). [3] T. Tomanic et al., Physica B, (2008). [4] A. Schröder et al., Phys. Rev. Lett., (1998). [5] A. Schröder et al., ature, (2008). [6] Y. akai et al., Phys. Rev. B, (R) (2015). [7] A. A. Druzhinin et al., J. Magn. Magn. Mater., (2015). [8] T. Toliński et al., J. Magn. Magn. Mater., (2013). [9] E. J. ijhof et al., Phys. Rev. B, (1986). [10] G. Toulouse, Commun. Phys., (1977); J. Villian J, J. Phys. C, (1977); T. Oguchi T, J. Phys. Soc. Jpn., (1983); T. Morita, Physica, 98A 566 (1979). [11] H. Kitatani et al., J. Phys. Soc. Jpn., (1986) their references [12] I. Mirebeau et al., Phys. Rev. B, (1990). [13] G. Cao et al., J. Magn. Magn. Mater., (2006). [14] R. H. Colman A. C. Mclaughlin, Phys. Rev. B, (2012). [15] F. C. Sá Barreto A. L. Mota, Physica A, (2015). [16] Ghliyem M et al. 2014, Physica A, [17] Y. ishiyama, Physica A, (2013). [18] T. Kaneyoshi, Physica A, (2012). [19] T. Kaneyoshi, J. Magn. Magn. Mater., (2016). [20] T. Kaneyoshi, J. Magn. Magn. Mater., (2013). [21] Ü. Akinci, J. Magn. Magn. Mater., (2013). [22] P. G. de Gennes, Solid State Commun., (1963). [23] R. Blinc B. Zeks, Adv. Phys., (1972). [24] R. J. Elliot A. P. Young, Ferroeletrics, 7 23 (1974). [25] D. S. Fisher, Phys. Rev. Lett., (1992). [26] V. K. Saxena, Phys. Rev. B, (1983). [27] V. K. Saxena, Phys. Lett. A, (1982). 9

10 [28] O. Canko et al., J. Magn. Magn. Mater, (2005). [29] R. J. Creswick et al., Phys. Rev. B, (1988). [30] A. Saber et al., J. Phys. Condens. Matter, (1999). [31] E. F. Sarmento et al., J. Magn. Magn. Mater., (1992). [32] T. Bouziane M. Saber, J. Magn. Magn. Mater., (2009). [33] R. Honmura T. Kaneyoshi, J. Phys. C, 12, 3979 (1979). [34] D. A. do ascimento et al., Physica A, (2013). [35] D. A. do ascimento et al., J. Magn. Magn. Mater., (2012). [36] A. C. Paz et al., Sci. Amazon. 3, 78 (2014). [37] F. C. Sá Barreto et al., Ferroelectrics, (1981). [38] M. A. eto J. R. de Sousa, Phys. Scr. 87, (2013). [39] S. Mal et al., ature Commun. 5, 4435 (2014). [40] T. Eckert E. Bartsch, Phys. Rev. Lett. 89, (2002). [41] E. S. Bǒzin et al., Scientific Reports 4, 4081 (2014). 10

11 20 (a) 20 (b) 15 h=6.20; δ=0.7 h=6.30; δ=0.7 h=6.40; δ= h=6.20; δ=0.2 h=6.20; δ=0.4 h=6.20; δ=0.7 χs/j 10 χs/j k B T/J k B T/J χs/j (c) h=6.05; δ=0.01 h=6.35; δ=0.7 h=6.20; δ= k B T/J h δ=0.1 δ=0.3 δ=0.5 δ=0.7 δ=0.9 (d) k B T/J FIG. 1: Dependence of the susceptibility χ s /J y, as a function of the temperature. In the Fig. (1.a) we show this dependence by setting the value of the transverse magnetic field δ = 0.7 we vary the longitudinal magnetic field. In the Fig. (1.b) we show this dependence by setting the value of the longitudinal magnetic field δ = 6.20 we vary the transverse magnetic field. In the Fig. (1.c), we vary both the fields. In the Fig. (1.d), we reproduce the h T phase diagram obtained by eto de Sousa [38] we show for the reentrant region. 11