EPJ Web of Conferences 40, 03002 (2013) DOI: 10.1051/epjconf/20134003002 © Owned by the authors, published by EDP Sciences, 2013 Magnetic properties of Ce3+ in PbCeA (A= Te, Se, S) S. Isber1, X. Gratens2, S. Charar3 and Z. Golacki4 1 American University of Beirut, Department of Physics, Bliss Street, PO Box 11-0236 Beirut, Lebanon Instituto de Física, Universidade de São Paulo, 05315-970, São Paulo, Brazil 3 Laboratoire Charles Coulomb UMR 5221, Université Montpellier 2, Place Eugène Bataillon - CC069, F-34095 Montpellier Cedex 5 - France 4 Institute of Physics, Polish Academy of Sciences, Pl.02-668 Warsaw, Poland 2 Abstract. The magnetic susceptibility of Pb1-xCexA (A = S, Se and Te) crystals with 0.006 ≤ x ≤ 0.036 were studied in the temperature range from 20 mK up to room temperature. X-band (~9.5 GHz) Electron Paramagnetic Resonance (EPR) showed small shifts in the effective Landé factors that were attributed to crystal-field admixture. The EPR measurements were correlated with the magnetic susceptibility data and resulted in estimating the crystal-field splitting E(8) - E(7) of the lowest 2F5/2 manifold for Ce3+ ions in PbA (A = S, Se and Te) of about 340 K, 440 K and 540 K for Pb1-xCexTe, Pb1-xCexSe, and Pb1-xCexS, respectively. The values for the crystal-field splitting deduced from the magnetic data were found to be in agreement with the calculated ones based on the point charge model. Moreover, the deHaas van-Alphen magnetic oscillations in the susceptibility measurements of Pb1-xCexTe (x~ 0.05 and 0.07) were observed at ultra-low temperature (20 mK); The oscillations were investigated and the values of the oscillatory period for Pb1-xCexTe (x = 0.0048 and 0.007) are reported. 1 Introduction Lead salts PbS, PbSe, and PbTe establish an significant group of Diluted magnetic semiconductors (DMS) consisting of a magnetic semiconducting host, in which fraction of cations is substituted by magnetic rare-earth or transition metal ions. In past few decades, various compound semiconductors from III-V, II-VI, IV-VI, etc. have been identified as potential candidates in thermoelectric materials with a high figure-of-merit, optoelectronic and spintronic applications [1-4]. The amount of knowledge, which has been accumulated, is less in IV-VI DMS than compared to the situation of IIVI DMS [5-13]. Magnetization, Magnetic Susceptibility and Electron Paramagnetic Resonance (EPR) measurements are poweful tools to study magnetic proprieties and electronic structures of DMS. Most of the investigated DMS systems in the literature are doped with transition metal and rare-earth ions having S (L = 0) states, such as Mn2+, Gd3+ or Eu2+ [6-13]. However, the study of IV-VI DMS containing a fraction of non S-state rare-earth ions is much less explored [14-16]. For this class of DMS, the magnetic properties could be strongly affected by the presence of crystal-field splitting of the rare-earth states. For the rock salt structure of PbA (A = S, Se and Te), the crystal field at the cation site has cubic (Oh) symmetry and a crystal-field splitting of the 2 F 5/2 manifold of Ce3+ ion is expected into a 7 doublet and a 8 quadruplet, the doublet being the ground state. Electron Paramagnetic Resonance (EPR) of Ce3+ ions located in cubic symmetry site was observed at 4.2 K in Pb1-xCexA (A = S, Se and Te). The analysis of the EPR spectra lead to Landé g = 1.331, 1.361 and 1.400 for Pb1xCexS, Pb1-xCexSe and Pb1-xCexTe, respectively [17]. The theoretical value for 7 is g = 1.429. The present work is focused on the magnetic susceptibility of Pb1-xCexA (A = S, Se and Te). The present data were used to estimate the cubic crystal-field splitting of the 2F5/2 manifold of Ce3+ in the three compounds. 2 Experimental Experiments were carried out on an n-type single crystal of Ce doped PbA (A = Te, Se, S) grown by the Bridgman method [18]. The crystals were not intentionally doped with any other element. The solubility limit of Ce in the PbA matrix is about 10%. X-ray diffraction analysis indicated a fcc structure for all samples. The EPR measurements in the X-band (~9.56 GHz) were performed on a parallelepidic sample with the longest axis along the [001] [17]. These Pb1-xCexA crystals of few mm3 were pulled out from the as grown ingots. The extracted samples were cut along the (001) planes in form This is an Open Access article distributed under the terms of the Creative Commons Attribution License 2.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Article available at http://www.epj-conferences.org or http://dx.doi.org/10.1051/epjconf/20134003002 EPJ Web of Conferences of small cubes. The Cerium content (x) in the samples was deduced from the saturation of the magnetization curves at high magnetic fields and by electron microprobe measurements. The magnetization was measured at 2 K in magnetic fields H up to 7 Tesla using a SQUID magnetometer. The magnetic susceptibility was measured at 0.1 Tesla over the temperature range from 2 to 300 K. The errors on the data were estimated to 2%. The samples of Pb1-xCexA (A = Te, Se and S) were all ntype with carrier concentrations ranging from 6×10 18 to 1020 cm-3. The magnetic susceptibility at 20 mK was measured in a system using a force magnetometer, a plastic dilution refrigerator, and a 9 Tesla superconducting magnet. The fcc structure of the samples were checked by X-Ray diffraction and the cerium content were determined from the magnetization saturation at high external magnetic field; the calculation of the Ce concentration (x) in the sample were corrected for the lattice diamagnetism (D = -3.77×10-9 m3/Kg for PbS [19], D = -4.52×10-9 m3/Kg for PbSe [20], and D = -3.77×10-9 m3/Kg for PbTe [17]); The obtained average concentrations of cerium atoms were also checked by Microprobe measurements, as given in Table I; An excellent agreement between the Ce concentrations values by microprobe to those obtained from magnetic measurements indicated that cerium is present only in its trivalent (Ce3+) state in our Pb1−xCexA samples. The temperature variation of  for our Pb1xCexA (A = S, Se and Te) samples were found to strongly deviates from the free Ce ion paramagnetic behavior; it was also noticed that the experimental magnetic susceptibility fits the predicted values of the susceptibility for the doublet (7) only at temperatures below 30 K indicating that the Kramer doublet (7) lies below the quadruplet (8). T able I: Values of x, D and  for PbCeA (A= Te, Se, S). The concentration of Ce3+ ions (x) was obtained by three different techniques, (a)Modified Brillouin function, (b)Three-parameters fit of the magnetic susceptibility, (c)Microprobe measurements, (d) The agreement factor, defined as Se = [ /(N-L)]1/2, N is the number of data points, K is the number of parameters and  = (Robs – Rcalc)2 is the sum of the squares of the residuals. x(a)/ x(c)/x(b) (%)  (K) D (m3/Kg) Se(d) Se Se Se S S S Te 1.2/--/1.1 2.0/1.9/1.8 3.4/3.1/3.6 0.9/--/0.9 1.9/--/1.8 2.1/2.2/1.9 0.6/--/0.6 407 442 420 545 542 529 340 -3.8110-9 -3.6810-9 -4.2810-9 -3.6710-9 -3.0210-9 -3.310-9 -3.6410-9 7.210-3 3.210-3 7.310-3 2.810-3 3.510-3 3.210-3 6.710-3 2   N This result was also confirmed by the values of the Landé factors obtained by EPR measurements [17]. In order to simulate our experimental magnetic susceptibility data, we accounted for the presence of the quadruplet (8) 2  eff  B k BT x D (1) Where B and k B are the Bohr magneton and Boltzmann constant, respectively, and N is the number of cation per gram. The effective magnetic moment eff of Ce3+ is given by: 2  eff 3 Experimental results and discussions A laying at energy  E(8) - E(7) above the doublet (7). In such a case, the magnetic susceptibility of Pb1-xCexA can be described by the Van-Vleck expression given by:  2  eff 32 k B T   1  exp(   / k B T )    5  26 exp(   / k B T )   ( free )   21 1  2 exp(   / k B T )      (2) Where 2  eff ( free )  1 3 2 g J J ( J  1) is the effective 3+ magnetic moment of the Ce free ion which was obtained from the experimental g-value of the 7 doublet and gJ = 3/5 g. The values for  eff2 ( free ) were found to be 2.33, 2.38 and 2.45 for Pb1-xCexS, Pb1-xCexSe and Pb1respectively. In Eq. (2),  represents the cubic Crystal-Field Splitting (CFS), (i.e. the energy difference between 8 and 7). In our simulation of the magnetic susceptibility data, the effect of the exchange integral between Ce3+ ions was ignored as its magnitude is very small (|J/kB| ≤ 1 K) when compared to the CFS. In order to obtain the CFS value , the effective magnetic moment for each of the three compounds was fitted to Eq. (1) in the temperature range from 2 to 300 K with three fitting parameters x,  and D. The obtained values of x,  and D are given in Table I. Figure 1 shows the xCexTe, 2 temperature variation of R =  eff 2  eff ( free ) for the three compounds. The obtained values of the fitting parameters x and D are given in Table I. The values of CFS ( were found to be 536 K, 443 K and 340 K for Pb 1-xCexS, Pb1xCexSe and Pb1-xCexTe, respectively. We interpreted the results of our measurements on magnetic susceptibilities and electron paramagnetic resonance of Ce doped lead chalcogenides by fitting the crystal field and diamagnetic susceptibility parameters. The temperature dependence of the single ion susceptibilities was computed according to Eq. (2). Furthermore, the obtained values of the CFS were compared to the ones calculated from the point charge model [21, 22]. In our calculations, we accounted for the six doubly negative charges of the nearest neighbor anions located at a /2 from the Ce3+ ions, where a = 0.6454 nm, 0.6124 nm and 0.5936 nm for PbTe, PbSe and PbS, respectively. We used the Hartree-Fock method with a <r4> value of 3.5 a.u. for the fourth power of the 4f electron radius [23]. 03002-p.2 Joint European Magnetic Symposia 2012 1.0 Pb1-xCexSe 0.6 x = 0.012 x = 0.020 x = 0.034 Pb1-xCexS 2  eff / eff ( Free ) 2 0.8 0.4 x = 0.009 x = 0.021 Pb1-xCexTe x = 0.006 0.2 0 50 100 150 200 250 300 T (K) F ig. 1: Experimental data (symbols) and theoretical (continuous lines) reduction of the effective magnetic moment of Ce3+ in Pb1-xCexS, Pb1-xCexSe and Pb1-xCexTe samples. The effective magnetic moment was normalized to its value for the free Ce3+ ion. Fewer points of the data were displayed for a better view. The lines are the fit between 2 K and 300 K using the expressions in Eqs. (1) and (2). The values of the CFS obtained from the point charge model were  = 400 K, 340 K and 262 K for Pb 1-xCexS, Pb1-xCexSe, and Pb1-xCexTe, respectively which are about 100 K below the values obtained from the magnetic susceptibility results. In our calculation, the 0.2 12 13 H  <001> d M /d H T = 20 mK 0.0 -0.1 (a) -0.2 0.12 0.18 0.24 -1 0.30 0.36 0.42 -1 d M /d H H (Tesla ) x = 0.005 FFT Amplitude of H  <001> T = 20 mK 90 (b) 0 40 80 120 160 4 Conclusion On the basis of the simplified Van Vleck paramagnetic susceptibility model of the Ce3+ ion in a octahedral site symmetry, the temperature dependence of the magnetic susceptibilities of three Ce3+ chalcogenide compounds PbA (A = Te, Se, and S), within the temperature range 4.2-300 K have been investigated and interpreted. In our study, the CFS parameter  were found to be 536 ± 50 K, 440 ± 50 K and 340 ± 50 K for Pb1-xCexS, Pb1-xCexSe and Pb1-xCexTe, respectively. It was also shown that the doublet 7 lies below the quadruplet 8. A good explanation of our observed magnetic data shows that the experimental magnetic susceptibilities including the comparison of the Crystal-Field Splitting  to: (i) the spin-orbit splitting (~1600 K) between the 2F5/2 and 2F7/2 manifolds and (ii) to the point charge model. This simple model can be further used along with other physical measurements such as EPR, optical spectroscopy or low temperature magnetization, as basis for reconstruction the magnetic functions of the compound. The comparison of the obtained CFS values to those calculated from the point charge model was satisfactory; we expect that this agreement can be further improved by accounting for the magnetic dipole interactions, the exchange interactions and the presence of distorted Ce3+ sites. Further analysis on the dHvA results will follow in future communications. x = 0.005 14 0.1 spin-orbit constant was taken from the literature (1600 K) while the Landé factors were deduced from low temperatures EPR measurements. The existence of four types of sites was attributed to the presence of Pb vacancies acting as charge compensators near the Ce3+ ions. The results were explained by the formation of two lead vacancies per five Ce3+ ions. The presence of these vacancies could explain the existence of two types of antiferromagnetic nearest neighbor (NN) Ce3+ pairs previously reported by magnetization steps measurements [15]. We propose that cerium pairs with isotropic exchange interaction are due to two NN Ce3+ ions located into cubic sites with no lead vacancy in the near vicinity. The other group of NN pairs is ascribed to pairs with a nearby lead vacancy affecting both strength of the exchange interaction and its anisotropy. Figure 2 represents the Field dependence of the oscillatory part of the magnetic susceptibility at 20 mK for H || [001]; in this figure the de Haas-van Alphen (dHvA) oscillatory part of the magnetic susceptibility were clearly observed (upper figure); These oscillations are periodic in 1/H, and the Fourier transform of the oscillatory part revealed one frequency f ~ 90 Tesla as shown in the lower part of figure 2. The analysis of the oscillatory component in the direction H || <001> show periods of 1.12×10-2 and 1.23×10-2 Tesla-1 for Pb1-xCexTe samples with x = 0.0048 and 0.007, respectively. 200 H (Tesla) F ig. 2: Field dependence of the oscillatory part of the magnetic susceptibility at 20 mK for H || [001] (upper figure); Fourier transform of the oscillatory part of the susceptibility at 20 mK between 0 and 200 Tesla, revealing a single frequency f~ 90 Tesla (lower figure). References 03002-p.3 EPJ Web of Conferences 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. G. Springholz, T. Schwarzl, W. Heiss, T. Fomherz, G. Bauer, M. Aigle, H. Pasher, and I. Vavra, Physica E (Amsterdam) 13, 876 (2002). G. Grabecki, J. Wrobel, T. Dietl, E. Papis, E. Kaminska, A. Piotrowska, G. Springholz, and G. Bauer, Physica E (Amsterdam)13, 649 (2002). P. Banerjee, J. Pal, and B. Ghosh, J. Mater. Sci. 40, 1333 (2005). G. Bauer and H. Pascher, in Diluted Magnetic Semiconductors, edited by M. Jain (World Scientific, Singapore, 1991). A. Lewicki, J. Spalek and A. Mycielski, Solid State Phys. 20, 2005 (1987). A. Jwardowski, J. Appl. Phys. 67, 5108 (1990). S. Isber, P. Masri, S. Charar, X. Gratens and S. K. Misra, J. of Physics. C9, 10023 (1997). V. Bindilatti, E. ter Haar, N.F. Oliviera, Jr., M.T. Liu, Y. Shapira, X. Gratens, S. Charar, S. Isber, P. Masri, M. Averous, Z. Golacki, and E.J. McNiff, Jr. Phys. Rev B57, 7854 (1998). S. Isber, S. Misra, S. Charar, X. Gratens, M. Averous, Z. Golacki, Phys. Rev. B56, 13724 (1997). X. Gratens, E. ter Haar, V. Bindilatti, NF. Oliveira, Y. Shapira, MT. Liu, Z. Golacki Z, S. Charar, A. Errebbahi, J. PHYS- Cond. Matt. 12, 3711 (2000). V. Bindilatti, E. ter Haar, N.F. Oliveira, M.T. Liu, Y. Shapira, X. Gratens, S. Charar, S. Isber, P. Masri, M. Averous, Z. Golacki, E. J. McNiff, Phys. Rev. B57, 7854 (1998). X. Gratens, S. Isber, S. Charar, C. Fau, M. Averous, S. K. Misra, Z. Golacki, M. Ferhat and J. Tedenac, Phys. Rev. B55, 8075 (1997). M. Gorska and J. R. Anderson, Phys. Rev. B38, 9120 (1988). S. Isber, S. Charar, X. Gratens, C. Fau, M. Averous, S. K. Misra and Z. Golacki, Phys. Rev. B54, 7634 (1996). X. Gratens, E. ter Haar, V. Bindilatti, N.F. Oliveira, Y. Shapira, Z. Golacki, J. Magn Magn Mater 226, 2036 (2001); X. Gratens, E. ter Haar, V. Bindilatti, N.F. Oliveira, Y. Shapira, Z. Golacki, Physica B 284-288, 1519 (2000). X. Gratens, V. Bindilatti, N.F. Oliveira and Z. Golacki, Physica B 329–333, 1245 (2003). A.M. Gennaro, G. B. Martins, C. Rettori, G. E. Barberis and C.Y. An in Rare Earth Doped Semiconductors, Edited by G.S. Pomrenke, P.B. Klein and D.W. Langer, MRS proceedings V 301, 213 (Pittsburgh, Pennsylvania, 1993). X. Gratens, V. Bindilatti, V. A. Chitta, and N. F. Oliveira, Jr., S. Isber and Z. Golacki Phys. Rev. B79, 075207 (2009). N. Senguttuvan, N. Kidokoro, K. Ootsuka, M. Ishii, M. Kobayashi, T. Taira, Y. Sato, S. Kurimura J. Crystal Growth 229, 188 (2001). Deduced from magnetization measurements at 2 K. J. R. Anderson, G. Kido, Y. Nishina, M. Gorska, L. Kowalczyk and Z. Golacki, Phys. Rev. B 41, 1014 (1990). R. Visser, J. Andriessen, P. Dorenbos and C W E Van Eijk, J. Phys: Condens. Matter 5, 5887 (1993). 22. 23. 03002-p.4 K. R. Lea, M. J.M Leask, W. P. Wolf, J. Phys. Chem. Solids 23, 1381 (1962). A. J. Freeman and R.E. Watson, Phys. Rev. 127, 2058 (1962).