Ah. Dhahri*ab,
E. Dhahria and
E. K. Hlilc
aLaboratoire de Physique Appliquée, Faculté des Sciences de Sfax, Université de Sfax, BP 1171, 3000, Tunisia. E-mail: dhahridhahri14@gmail.com; Tel: +216 20 20 45 55
bFaculté des Sciences de Monastir, Avenue de l'environnement 5019, Monastir, Tunisia
cInstitut Néel, CNRS et Université J. Fourier, BP 166, 38042 Grenoble, France
First published on 13th February 2019
La0.67−xBixBa0.33MnO3 (x = 0 and 0.05) ceramics were prepared via the sol–gel method. Structural, magnetic and magnetocaloric effects have been systematically studied. X-ray diffraction shows that all the compounds crystallize in the rhombohedral structure with the Rc space group. By analyzing the field and temperature dependence of magnetization, it is observed that both samples undergo a second order magnetic phase transition near TC. The value of TC decreases from 340 K to 306 K when increasing x from 0.00 to 0.05, respectively. The reported magnetic entropy change for both samples was considerably remarkable and equal to 5.8 J kg−1 K−1 for x = 0.00 and 7.3 J kg−1 K−1 for x = 0.05, respectively, for μ0H = 5 T, confirming that these materials are promising candidates for magnetic refrigeration applications. The mean-field theory was used to study the magnetocaloric effect within the thermodynamics of the model. Satisfactory agreement between experimental data and the mean-field theory has been found.
From this viewpoint, this paper reports the structural, magnetic and magnetic entropy change of Bi-substituted perovskite manganites La0.62Bi0.05Ba0.33MnO3. It is found that these materials show quite large magnetic entropy changes induced by low magnetic field changes.
Fig. 1 A flow chart illustrating the processing procedure for La0.67−xBixBa0.33MnO3 powders preparation. |
The phase purity and structure of sample were identified by X-ray powder diffraction at room temperature using a Siemens D5000 X-ray diffractometer with a graphite monochromatized CuKα radiation (λCuKα = 1.5406 Å) and 20° ≤ 2θ ≤ 90° with steps of 0.02° and a counting time of 18 s per step. According to our measurements, this system is able to detect up to a minimum of 3% of impurities. The structure analysis was carried out using the Rietveld method with FULLPROF software (version 0.2-Mars 1998-LLB-JRC).8 Scanning electron microscopy (SEM) using a Philips XL30 equipped with a field emission gun at 20 kV was used to characterize La0.67−xBixBa0.33MnO3 morphologies.
Magnetization (M) versus temperature (T) and magnetization versus magnetic field (μ0H) were performed by using BS1 and BS2 magnetometers developed in Louis Neel Laboratory at Grenoble. The isothermal curves were determined in the magnetic field range of 0–5 T. The temperature interval is fixed to 2 K in the vicinity of the Curie temperature (TC). The temperature steps were smaller near TC and larger further away.
Fig. 2 X-ray powder diffraction patterns and Rietveld refinement for the compounds La0.67−xBixBa0.33MnO3 (x = 0 and 0.05) at room temperature. |
Crystallographic data | Refinement conditions | Average crystallite size | Strain ε (%) | Average grain size D (nm) | ||||||
---|---|---|---|---|---|---|---|---|---|---|
Lattice parameters | Cell volume | Discrepancy factors | Debye–Scherrer technique Ds (nm) | Williamson–Hall technique Dw (nm) | ||||||
a = b (Å) | c (Å) | V (Å3) | Rwp % | Rp % | RF % | χ2 % | ||||
a a and c: hexagonal cell parameters; V: cell volume; Rwp, Rp; RF: the residuals for, respectively, the weighted pattern, the pattern and the Bragg structure factor; χ2: the goodness of fit. The numbers in parentheses are estimated standard deviations to the last significant digit. | ||||||||||
x = 0 | ||||||||||
5.5160 (3) | 13.5023(1) | 355.78(2) | 4.12 | 3.43 | 2.78 | 1.85 | 38 | 49 | 0.16 | 220 |
x = 0.05 | ||||||||||
5.50184 (2) | 13.5141(4) | 356.40(1) | 5.23 | 3.42 | 3.42 | 1.25 | 35 | 40 | 0.14 | 260 |
In order to quantitatively discuss the ionic match between A and B sites in perovskite compounds, a geometrical quantity, noted Goldschmidt tolerance factor (t), is usually introduced and is defined as:10
The value of average crystallite size was estimated from the full width at half maximum (FWHM) of X-ray diffraction peaks. The effects of synthesis, instruments and processing conditions were taken into consideration while making the calculation of crystallite size. The dependence of the size effect is given by Scherrer's formula: , where λ is the wavelength of CuKα radiation (λ = 1.5406 Å), K is grain shape factor (=0.89) and Ds is the thickness of the crystal. Using the Williamson–Hall (W–H) method,13 the average values of both Dw and lattice strain (ε) can be obtained from the intercept and the slope of the following relation, respectively,
Fig. 3 shows the SEM photograph of the compounds. The samples contained connected particles with hexagonal shape and clear grain boundaries. These particles are largely agglomerated with a broad size distribution. The average value of thickness of both compounds is listed in Table 1.
Fig. 3 Scanning electron micrograph for La0.67−xBixBa0.33MnO3 (x = 0 and 0.05) manganite. The inset: dispersion histogram. |
After measuring the diameters of all the particles in SEM image, the size distribution histogram is fitted with the log-normal function expressed as:
〈D〉 = D0exp(σ2/2) |
σD = 〈D〉[exp(σ2) − 1]1/2 |
The results analysis showed 〈D〉 ≥ 397.48 μm and σD = 291.53 μm.
The inset of Fig. 4 shows the temperature dependence of the inverse magnetic susceptibility of x = 0 and x = 0.05. It could be fitted to the Curie–Weiss law just above TC (the PM region): χ = C/T − θCW, where θCW is the Weiss temperature and C is the Curie constant defined as: , where μ0 = 107 H m−1 is the permeability, g is the Landé factor, μB = 9.27 × 1024 J T−1 is the Bohr magneton, kB = 1.38 × 1023 J K−1 is the Boltzmann constant, J = L + S is the total moment and μeff is the effective paramagnetic moment. We can determine the effect of paramagnetic moment (μexpeff) from the curie constant. The theoretical μcalceff is estimated using the following expression: . The parameters μexpeff and μcalceff are summarized in Table 2. It is found that the μexpeff is greater as compared to μcalceff. This discrepancy validates the formation of ferromagnetic spin clusters within the paramagnetic state.19 A linear fit yields positive Curie–Weiss temperature θCW = 312 K (x = 0.05). This result confirms a mean FM interaction between spins for all samples (Table 2). Moreover, this value is higher than TC, which may be due to the existence of short range FM ordering.19
x | TC (K) | W/W0 (10−2) | θCW (K) | θMn–O–Mn (°) | dMn–O (Å) | μexpeff (μB) | μcalceff (μB) |
---|---|---|---|---|---|---|---|
0 | 340 | 4.73 | 348 | 165.32 | 1.959 | 5.32 | 4.586 |
0.05 | 306 | 4.66 | 312 | 165.12 | 1.965 | 5.81 | 4.586 |
The structure analysis shows that the unit cell becomes slightly larger as the 6s2 lone pair character becomes dominant, it has been shown that the Bi–O bond is shorter than the La–O, despite of the similar ionic radius of La3+ and Bi3+ ions.20 This can be interpreted as arising from the rather covalent character of the Bi–O bonds. The electronegativity of Bi enhances hybridisation between 6s2 of Bi3+ orbitals and 2p of O2− orbitals and this hybridisation produces a local distortion. It is observed that transition temperature TC decreases with increase in Bi ratio. This is presumably due to tilts the MnO6 octahedra, resulting in a reduced overlap between the Mn-3d and O-2p orbitals.21 It should also be noted that the La0.67Ba0.33MnO3 sample is ferromagnetic while Bi0.67Ba0.33MnO3 is antiferromagnetic, indicating a competition between the double exchange and the antiferromagnetic super exchange in these compounds can decrease TC. This phenomenon has been observed in the compound Bi0.6−xLaxCa0.4MnO3.22
The magnetic entropy is related to the magnetization M, magnetic field strength μ0H and absolute temperature T through the Maxwell relation:
In the case of magnetization measurement in small discrete magnetic fields and temperature interval ΔT, ΔSM can be approximated to:
The −ΔSM(T) data calculated from the M(μ0H) curves (inset in Fig. 5) at different magnetic fields for the La0.67−xBixBa0.33MnO3 (x = 0 and 0.05) are plotted in Fig. 5. The compounds exhibit large changes in magnetic entropy around Curie temperature (TC), which is a characteristic property of simple ferromagnets due to the efficient ordering of magnetic spins at the temperature induced by magnetic field.25 Large magnetic entropy changes ΔSmaxM are reported for all the samples and are summarized in Table 3. The magnitude of (−ΔSmaxM(T)) for all samples increases with increasing the applied magnetic field (inset of Fig. 5). For example, the maximum magnetic-entropy value increases from 2.37 J kg−1 K−1 for x = 0.00 to 2.8 J kg−1 K−1 (2T) and 5.8 J kg−1 K−1 for x = 0.00 to 7.3 J kg−1 K−1 for x = 0.05 respectively (5T). Guo et al.26 indicated that the large magnetic entropy change in perovskite compounds could originate from the spin–lattice coupling in the magnetic ordering process. Strong coupling between spin and lattice is corroborated by the observed significant lattice change accompanying magnetic transition in perovskite manganites.27 The lattice structural change in the Mn–O bond distance as well as in the 〈Mn–O–Mn〉 bond angle would in turn favor the spin ordering. Thus a more abrupt variation of magnetization near Curie temperature (TC) occurs, resulting in a large magnetic entropy change as a large MCE.
Composition | TC (K) | (−ΔSmaxM) (J kg−1 K−1) | RCP (J kg−1) | μ0H (T) | Ref. |
---|---|---|---|---|---|
Gd5(Sr2Ge2) | 275 | 18.5 | 535 | 5 | 24 |
Gd | 294 | 10.2 | 410 | 5 | 25 |
La0.5Sm0.1Sr0.4Mn0.95In0.05O3 | 308 | 4.50 | 193.48 | 5 | 26 |
La0.67Sr0.33Mn0.9Cr0.1O3 | 328 | 5 | — | 5 | 27 |
La0.67Ba0.33MnO3 | 346 | 5.80 | 151 | 5 | This work |
La0.62Bi0.05Ba0.33MnO3 | 310 | 7.30 | 209 | 5 | This work |
Gd | 297 | 4 | 120 | 2 | 28 |
MnFeP0.45As0.5 | 300 | 14.5 | 188 | 2 | 28 |
La0.7Sr0.3Mn0.95Ti0.05O3 | 308 | 2.2 | 90 | 2 | 29 |
La0.7Sr0.3Mn0.9Fe0.1O3 | 260 | 1.7 | 83 | 2 | 30 |
La0.67Ba0.33MnO3 | 343 | 2.37 | 39 | 2 | This work |
La0.62Bi0.05Ba0.33MnO3 | 308 | 2.80 | 80 | 2 | This work |
The change of magnetic entropy can be also calculated from the field dependence of the specific heat by the following integration:
From this equation, it determine the change of specific heat induced by the external magnetic field zero to μ0H as:
Fig. 6 shows the temperature dependence of ΔCp(μ0H, T) under different field variations for the samples (for example x = 0.05) calculated from the ΔSM(μ0H, T). The ΔCp(μ0H, T) undergoes a sudden change from positive to negative around TC with a positive value above TC and a negative value below TC. The maximum/minimum value of ΔCp(μ0H, T) observed at 320/300 K, exhibits an increasing trend with applied field and is obtained to be 122.4/−115.43 J kg−1 K−1 for x = 0.05 at 5 T.
Fig. 6 Change of specific heat of the samples as a function of temperature at different magnetic field. |
It should be noted that (−ΔSmaxM) is not the only parameter deciding about an applicability of material. To estimate if a material can be a good candidate for magnetic refrigeration (MR), Gschneidner and Pecharsky28 defined the relative cooling power (RCPS), which is the important index which is used to evaluate the cooling efficiency of a magnetic refrigerant. It is defined as the product between the maximum values of the magnetic entropy change (−ΔSmaxM) and the full width at half maximum δTFWHM of the magnetic entropy change curve (RCP(S) = −ΔSmaxM × δTFWHM).29 This parameter corresponds to the amount of heat that can be transferred between the cold and hot parts of the refrigerator in one ideal thermodynamic cycle. The results are summarized in Table 3. Fig. 7 shows the absolute value of RCPS and δTFWHM for the sample (x = 0.05) versus applied field at 310 K. It can be seen from this figure that RCPS (δTFWHM) increases monotonically as the field increases. The value of RCP is about 51% of Gd at 294 K for μ0H = 5 T.30 To evaluate the applicability of ours samples as a magnetic refrigerant, the obtained values of ΔSM in our study are compared in Table 3 with those reported in the literature for several other magnetic compounds.31–37
Heff = H + Hexch = H + λM |
Amaral et al. proposed a model based on mean field theory and presented an approach of applying this method scenario to isotherm magnetization M(T, H) measurements.39 In our study, it consider the general mean field law:40 M(H, T) = BJ[(H + Hexch)/T], the Brillouin function BJ is written as: , where , J is the total angular momentum in the lattice, g is the gyromagnetic factor (landé factor), μB is the Bohr magnetron and kB is the Boltzmann's constant. The mean field exchange parameter λ is not predetermined. Then for corresponding values with the same (H + Hexch)/T, M is also the same, the value of the inverse BJ−1(M) function,41
The study of the exchange field induced by the magnetization change makes it possible to find the value of the average field exchange parameter λ. Fig. 8 shows H/T versus 1/T for some of the values of M (5 emu per g per step) from 266 K to 342 K for x = 0.05. According to the mean scaling method such H/T versus 1/T curves should show a series of straight lines at different temperature. The linear relationship between H/T and 1/T is kept. Linear fits are then easily made to each isomagnetic line. Typically, the interpolation step was of 1 emu g−1. The slope of this isomagnetic line, will then give the exchange field (Hexch).
Fig. 8 H/T versus 1/T curves with constant values of magnetization per curve for La0.62Bi0.05Ba0.33MnO3 compound. |
For all compounds in the paramagnetic domain or the materials of domain ordered such as anti-ferromagnetic, it can always expand increasing M in powers of H, or H in powers of M. In this latter approach it stop at the third order and considering that the magnetization is an odd function of field, it can write:42
Hexch = λ1M + λ3M3 |
Fig. 9 shows the evolution of the exchange field versus the magnetization for the La0.67−xBixBa0.33MnO3 (x = 0.05 for example). The experimental data should be included for the fit by eqn (Hexch = λ1M + λ3M3). The results show a very small dependence on M3 (λ3 = −1.3984 × 10−5 (T g emu−1)3), is found for this second order transition system, thus Hexch = λ1M with λ1 = 1.25 T g emu−1. After obtaining the mean field exchange parameter the next step of this method consists on building the scaling plot of M vs. (H + Hexch)/T (Fig. 10). It has successfully fitted the scaled magnetization data with the Brillouin function. From the scaling plot and the subsequent fit with the saturation magnetization equal to 72 emu g−1 (this value is close to the experimental one (Ms = 69 emu g−1 at 10 K) (inset b of Fig. 4), and the value of the total angular momentum of the manganite is J = 1.9.
Fig. 9 Exchange field versus magnetization for La0.62Bi0.05Ba0.33MnO3 sample, with the function λ1M + λ3M3 fit. |
Fig. 10 Mean-field scaling plot and Brillouin function fit. Inset 1/χ versus T of La0.62Bi0.05Ba0.33MnO3 compound. |
The magnetization measurements in the law temperature range show that the saturation magnetization is about Mtheors = JgμB = 3.67μB/f.u. This value is close to the experimental value Mexpers = 3.7μB at 10 K. The value of Mtheors per formula unit is given by: Mtheors = (nMn3+)MsMn3+ + (nMn4+)MsMn4+, where MsMn3+ = 4μB and MsMn4+ = 3μB are the magnetic moments, nMn3+ = 0.67 and nMn4+ = 0.33 are the contents of Mn3+ and Mn4+ ions respectively. It deduce that the total angular momentum of the compound is J = 1.835 (it was assumed that g = 2).
From a linear approximation of the susceptibility . Inset of the Fig. 10 shows the evolution of the inverse susceptibility (1/χ) versus temperature for the La0.62Bi0.05Ba0.33MnO3 compound. The straight line represents the Curie–Weiss law with λ = 1.25 T g emu−1 and TC = 306 K. The intersection is obtained at a temperature value equals to the critical one. At this point, the susceptibility becomes infinite, which corresponds to the ferromagnetic–paramagnetic transition. Fig. 11 shows how experimental data can be described using the mean-field method. A good agreement between the experimental M(T, μ0H) curves and the mean-field generated curves with the obtained parameters, except near the paramagnetic-ferromagnetic transition (TC), which are not adequately described. This result is probably due to the formation of magnetic domains and critical effects.
Fig. 11 Experimental magnetization versus μ0H (black symbols) of La0.62Bi0.05Ba0.33MnO3 sample and the interpolation using the mean field method (red lines). |
The linearity relation between H/T and 1/T in Fig. 8 improves that the mean field exchange parameter λ is independent of the temperature T. So the following equation will be simplified:
Fig. 12 shows the evolution of the magnetic entropy change (−ΔSM) data as a function of temperature at several magnetic applied fields for the La0.62Bi0.05Ba0.33MnO3 compound, by using the Maxwell relation and that basing on mean field theory. Both results are in good agreement, except close to TC, here an excepted small difference appears, due to the formation of magnetic domains and critical effects.
This journal is © The Royal Society of Chemistry 2019 |