Ahmed Dhahri*ab,
J. Laific,
Soumaya Gouadriad,
M. Elhadib,
E. Dhahria and
E. K. Hlile
aLaboratoire de Physique Appliquée, Faculté des Sciences de Sfax, Université de Sfax, BP 1171, 3000, Tunisia. E-mail: dhahridhahri14@gmail.com
bDepartment of Physics, College of Science and Humanities – Dawadmi, Shaqra University, Riyadh, Saudi Arabia
cPhysics Department, College of Science, Jouf University, P.O. Box: 2014, Sakaka, Saudi Arabia
dDepartment of Physics, College of Science, Princess Nourah Bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
eUniv. Grenoble Alpes, CNRS, Grenoble INP, Institut Néel, 38000 Grenoble, France
First published on 1st February 2022
We present a detailed study on the physical properties of La0.6Ba0.2Sr0.2Mn1−xNixO3 samples (x = 0.00, 0.05 and 0.1). The ceramics were fabricated using the sol–gel route. Structural refinement, employing the Rietveld method, disclosed a rhombohedral Rc phase. The magnetization vs. temperature plots show a paramagnetic–ferromagnetic (PM–FM) transition phase at the TC (Curie temperature), which decreases from 354 K to 301 K. From the Arrott diagrams M2 vs. μ0H/M, we can conclude the phase transition is of the second order. Based on measurements of the isothermal magnetization around TC, the magnetocaloric effects (MCEs) have been calculated. The entropy maximum change (−ΔSM) values are 7.40 J kg−1 K−1, 5.6 J kg−1 K−1 and 4.48 J kg−1 K−1, whereas the relative cooling power (RCP) values are 232 J kg−1, 230 J kg−1 and 156 J kg−1 for x = 0.00, 0.05 and 0.10, respectively, under an external field (μ0H) of 5 T. Through these results, the La0.6Ba0.2Sr0.2Mn1−xNixO3 (0 ≤ x ≤ 0.1) samples can be suggested for use in magnetic refrigeration technology above room temperature. The electrical resistivity (ρ) vs. temperature plots exhibit a transition from metallic behavior to semiconductor behavior in the vicinity of TM–SC. The adiabatic small polaron hopping (ASPH) model is applied in the PM-semiconducting part (T > TMS). Throughout the temperature range, ρ is adjusted by the percolation model. This model is based on the phase segregation of FM-metal clusters and PM-insulating regions.
Gadolinium (Gd) is a pure lanthanide element and is the first material that has a high MCE with a Curie temperature (TC) near room temperature (RT).5 Interestingly, further MCE investigations were performed for binary Gd–M compounds, such as Gd3SiGe2,6 which shows a MCE twice that of Gd.7 The researchers focused on finding new cheaper materials with larger MCEs. In this context, manganites with perovskite structure have certain advantages over Gd: their elements are not expensive, they are chemically stable, they have high resistivity and they exhibit a good MCE under low magnetic fields.8,9 Among these, ABO3 compounds are known as perovskite manganites. These materials have a general formula Re1−x3+Ax2+(Mn1−x3+Mnx4+)O32− where Re3+ is a rare earth element (Nd3+, La3+, Pr3+, Sm3+, …) and A2+ is an alkaline earth ion (Sr2+, Ba2+, Ca2+).10,11 They present some interesting properties, which make them very attractive materials for industrial applications.
The pure stoichiometric lanthanum manganite LaMnO3 is antiferromagnetic, insulating at 150 K, and the substitution of the rare earth element by a lower valence ion causes the oxidation of Mn3+ into Mn4+ to ensure electroneutrality in the material. It is followed by the appearance of macroscopic magnetization, i.e. a ferromagnetic coupling between the Mn3+ ions (t2g3e1g) and the Mn4+ ion (t2g3e0g). Substitution of La3+ by a divalent or monovalent ion can result in a wide Curie temperature range which can vary from 150 K to 375 K. Experimentally, manganites, in particular manganese oxides La1−xSrxMnO3 with x = 0.3, are well studied systems. They present an FM–PM transition accompanied by a metal–semiconductor transition close to TC. Several studies have been performed on the magnetocaloric properties of this compound, which exhibits a large change in magnetic entropy, with a narrow range of working temperatures in the vicinity of TC. In addition, several investigations have been carried out to estimate the substitution effects of the Mn site and these have shown that, for the La1−xSrxMnO3 family, even a low rate of substitution at the Mn site would induce a significant change in the properties of magnetotransport.12–14 The values of TC and ΔSM are generally affected by the partial substitution of manganese ions by certain transition metals, for example the In ion,15,16 the Al ion,17 etc. On the other hand, the substitution of the rare earth La3+ by certain metals lead to a significant change in the magnetic and electrical properties.18,19
Indeed, the substitution of the Re site with divalent ions20 proves the oxidation of Mn3+ to Mn4+, which is the origin of the ferromagnetic character.21,22 The magnetic coupling between Mn4+ and Mn3+ is usually governed by the movement of the electron, for example between the two partly filled d layers with a strong Hund’s coupling on site.
The issue of replacing magnetic and non-magnetic ions in the Mn site is very important. For example, the partial substitution of Mn3+ ions with Ni2+ ions modifies the ratio of the Mn4+–O2–Mn3+ network and leads to a decrease of the double exchange (DE) interactions.23 Several studies have been carried out24–29 to explain the relationship between the magnetotransport and magnetic properties of Re1−xAxMnO3 substituted by different elements on the Mn site. Based on the information given in this article, we have carefully discussed the physical properties in La0.6Ba0.2Sr0.2Mn1−xNixO3 (0 ≤ x ≤ 0.1) compounds.
The patterns of our samples were indexed in the rhombohedral (Rc) symmetry (no. 167), with (La, Ba, Sr): 6a (0, 0, 0.25), (Mn, Ni): 6b (0, 0, 0) and O: 18 (x, 0, 0.25). The different structural parameters are tabulated in Table 1. When the Ni substitution increases, the volume and the lattice parameters decrease. A similar behavior has also been previously observed.31
x | |||
---|---|---|---|
0.00 | 0.05 | 0.10 | |
Rc phase | |||
a (Å) | 5.4971(4) | 5.4824(3) | 5.4785(1) |
c (Å) | 13.4713(5) | 13.4572(1) | 13.4491(3) |
V (Å3) | 352.54(1) | 350.30(2) | 349.58(4) |
(O)Biso (Å2) | 0.85(2) | 0.54(3) | 2.61(6) |
(O)x | 0.4472(1) | 0.4526(3) | 0.4435(5) |
(La, Sr, Ba)Biso (Å2) | 0.967(1) | 0.352(3) | 0.853(5) |
(Mn, Ni)Biso (Å2) | 0.368(2) | 0.267(7) | 1.812(6) |
Discrepancy factors | |||
Rwp (%) | 2.62 | 3.51 | 4.11 |
Rp (%) | 3.22 | 1.09 | 4.20 |
RF (%) | 5.172 | 6.532 | 2.235 |
χ2 (%) | 2.42 | 1.61 | 1.37 |
This decrease is explained by the fact that the average ionic radius of the manganese site (0.599 ≤ rMn+Ni ≤ 0.592) decreases, which can be assigned to the formulation of a higher level of Mn4+, compared to Mn3+.
Taking into account the neutrality of La0.63+Ba0.22+Sr0.22+(Mn0.6−y3+Mn0.4−x+y4+)Nix2+O32− (0 ≤ x ≤ 0.1), these results can be justified as follows: whenever the x ratio of Ni is increased, the proportion of Mn3+ weakens by y = 2x, while the ratio of Mn4+ increases by x. In Table 2, we have listed the (Mn/Ni)–O–(Mn/Ni) bond angles (θ(Mn/Ni)–O–(Mn/Ni)) and the (Mn/Ni)–O bond lengths (d(Mn/Ni)–O). It is noted that θ(Mn/Ni)–O–(Mn/Ni) decreases linearly with the increase of x, whereas d(Mn/Ni)–O increases, leading to a tilting of the BO6 octahedrons.
x | 0.00 | 0.05 | 0.10 |
---|---|---|---|
θ(Mn/Ni)–O–(Mn/Ni) (°) | 166.2(1) | 164.98(9) | 162.82(6) |
d(Mn/Ni)–O (Å) | 1.9571(1) | 1.9587(2) | 1.9633(1) |
〈rB〉 (Å) | 0.645 | 0.647 | 0.649 |
tG | 0.974 | 0.973 | 0.972 |
D | 0.286 | 0.285 | 0.283 |
W (10−2) (arb. units) | 4.73 | 4.71 | 4.66 |
DSC (nm) | 55 | 50 | 49 |
〈DSEM〉 ± σD (μm) | 1.33 ± 0.44 | 0.95 ± 0.38 | 0.79 ± 0.31 |
DWH (nm) | 130 | 121 | 118 |
p (%) | 7.26 | 7.32 | 7.35 |
The perovskite structure can be distorted from the ideal cubic structure, which greatly affects the properties. These distortions are principally given by the relationship between the ionic radius of the cations, defined by the tolerance factor tG:32
(1) |
In our work, the tolerance factor tG decreases with the increase of Ni (Table 2).
The rate of rhombohedric deformation D% can be calculated employing the expression: where , a1 = a2 = a and 35,36
From Table 2, the value of D decreases with the decrease in mean rMn+Ni.
On the other hand, the average crystallite size DSC can be calculated using Scherrer’s relationship:37
(2) |
As in Scherrer’s method, the crystallite size values were determined from the Williamson–Hall equation:38
(3) |
Fig. 3 Typical SEM images for x = 0.00 and 0.10, (a) and (b) present the histograms of particle size and (c) shows the EDX analysis for x = 0.10. |
ImageJ software was employed to determine a statistical count of the grain size, which was performed on the SEM images. This technique consists of measuring the diameters of all the particles in the SEM image. Then, we adjusted these data using the log-normal function.
(4) |
Where, σ is the median diameter obtained from the data dispersions and D0 is the median diameter obtained from the SEM images. Fig. 3(a) and (b) present the grain number (counts) versus the particle size. Using the fit results, the mean diameter and standard deviation σD = 〈D〉[exp(σ2) − 1]1/2 were determined (Table 2).
It is remarkable that the average size of the particles obtained is greater than the average size of the crystallites determined by XRD. This can be explained by the fact that each particle observed by SEM is made up of several crystallites. The energy dispersive X-ray microanalysis (EDX) spectrum of La0.6Ba0.2Sr0.2Mn0.9Ni0.1O3 is presented in Fig. 3(c) as an example. This technique confirms the composition and purity of the samples. The spectra reveal the homogeneous distribution of La, Ba, Sr, Mn, Ni and O atoms over a wide surface area.
To evaluate the porosity of the compounds, we calculated the X-ray density, , where d is the bulk density, a is the lattice constant, M is the molecular weight and Na is Avogadro’s number39 (see Table 2).
Fig. 4 Evolution of magnetization vs. temperature under a magnetic field μ0H = 0.05 T for the compounds LBSMNO. The inset: dM/dT curve vs. T. |
The most important cause of the decrease in TC is the reduction of the one-electron bandwidth Wd. It is given as:43
(5) |
The increase in the Ni content results in a greater length of the d(Mn/Ni)–O bond and a reduction in the bond angle θ(Mn/Ni)–O–(Mn/Ni), and therefore a weaker bandwidth W (Table 2) which causes the decrease of TC.44,45 In fact, this reduction in W leads to a reduction in the FM coupling between neighboring manganese atoms.
In addition, it is interesting to understand the behavior of M vs. temperature in the FM region. In this context, and according to Lonzarich and Taillefer,46 magnetization obeys the theory of spin waves. At low temperatures, this theory is that the magnetization has multiplied in T3/2 (Bloch’s law) and over a wide range of temperatures in T2, yet in the vicinity of TC it varies as follows: (1 − T4/3/TC4/3)1/2.
In the FM region, the M data has been adjusted by this relation:
M(T) = M0 + M3/2T3/2 + M2T2 | (6) |
For T > TC, the PM region, the Curie–Weiss law was used to analyse the inverse of magnetic susceptibility (χ−1 = H/M):
(7) |
Fig. 5 The inverse magnetic susceptibility χ−1 versus temperature. The inset shows the field dependence of the magnetization curves at 10 K. |
x | |||
---|---|---|---|
0.00 | 0.05 | 0.10 | |
TC (K) | 354 | 320 | 301 |
θCW (K) | 358 | 326 | 308 |
μexpeff (μB) | 4.92 | 4.77 | 4.64 |
μtheeff (μB) | 4.48 | 4.34 | 4.19 |
μ0HC (10−3 T) | 9.7 | 8.9 | 7.7 |
Ms (emu g−1) | 84 | 79 | 75 |
Mr (emu g−1) | 13 | 8 | 5 |
R | 0.066 | 0.10 | 0.15 |
The experimental effective moment μexpeff can be calculated using the parameter C and the values are given in Table 3.
The theoretical effective paramagnetic moment for La0.6Ba0.2Sr0.2Mn1−xNixO3 compositions could be calculated using the following expression:
The percentage of Mn3+ and Mn4+ ions was calculated and checked by the conventional chemical method. The orbital moment is frozen (L = 0) for Mn3+ and Mn4+, so, the theoretical effective moment can be given as: with S = 1.5 for (Mn4+, 3d3), 2 for (Mn3+, 3d4) and 4 for (Ni2+, 3d8) and g = 2. The calculated values of μtheeff(Mn3+), μtheeff(Mn4+) and Ni2+ are 4.90 μB, 3.78 μB and 2.828 μB, respectively. From Table 3, a difference between the experimental and theoretical values of μeff can be observed. This can be clarified by the presence of FM clusters within the PM phase.48
The evolution of magnetization vs. μ0H for La0.6Ba0.2Sr0.2Mn1−xNixO3 (0 ≤ x ≤ 0.1) at 10 K is depicted in the inset of Fig. 5. At μ0H = 1.5 T, the compounds exhibit a constant value of M. The magnetic moments determined by the magnetization data are obtained to be 3.56 μB per formula unit for x = 0.00, 3.35 μB per formula unit for x = 0.05 and 3.1 μB per formula unit for x = 0.1.
The calculated magnetic moment can be determined by: Msp = ((0.67 − 2x) × 4 + (0.4 + x) × 3 + 2x) μB = (3.88 − 3x) μB/f.u.
The magnetic moments of Ni2+, Mn3+ and Mn4+ have 2, 4 and 3 μB, respectively. The Msp values are 3.88 μB for x = 0.00, 3.73 μB for x = 0.05 and 3.58 μB for x = 0.10. This reduction can be attributed to competition between the ferromagnetic and antiferromagnetic interactions. In addition, the Ni2+ ion at the M site influences the valence states of the manganese ions, i.e. it decreases the Mn3+/Mn4+ ratio. This proves a reduction in the double Zener exchange (DE), which leads to a decrease in magnetization.
We analyzed the hysteresis loops at 10 K (μ0H = ±5 T) to better understand the magnetic properties at low temperatures (Fig. 6). The curves are similar, with hysteresis loops which are weak, with reasonably good coercive fields (μ0HC), which decrease with the increase of Ni content . This reduction can be assigned to the decrease in spin dependent electron hopping. In the weak μ0H region, M grew significantly and reached saturation as the field increased. The value of Ms (saturation magnetization) can be estimated at high μ0H at about 5 T. From Table 3, it was found that Ms decreased, which may be due to antiferromagnetic alignments. Likewise, the insertion of Ni in the Mn site modifies the valence states of the manganese ions and decreases the level of Mn3+/Mn4+, which in turn weakens the DE interaction. The inset of Fig. 6 shows a zoomed in view of the central portion of M versus μ0H, at small μ0H. The determined values of μ0HC with Ni substitution are summarized in Table 3. The weak hysteresis loop with large saturation values confirms the characteristic soft FM behavior of the compounds. In this context, it can be concluded that our compounds may be applicable to read and write processes in high density recording media or for information storage.49
Fig. 6 Magnetic hysteresis curves at T = 10 K for the LBSMNO samples. The inset shows a zoomed in view of the central portion of the hysteresis at a low field. |
The remanence ratio (R) is utilized to comprise the isotropic nature of our investigated compounds. The values of remanence varied in the range of 5–13 emu g−1. The ratio (R) is given by: R = Mr/Ms. We have summarized the values of R in Table 3. The small obtained values confirm the isotropic natures.50 In this context, for magnetic recording and memory devices,51 it is advantageous to have higher remanence ratios. The values of R reveal an increasing trend with Ni2+ substitution.
We have presented in Fig. 7(a–c), the external μ0H variation of the isothermal magnetization at different temperatures around TC. M increases rapidly at low μ0H and then achieves saturation, which shows ferromagnetic behavior. Above TC, thermally unsettled magnetic moments yield to increase magnetizations linearly at high temperatures, which means PM behavior. This phenomena demonstrates the magnetic phase transition.
Fig. 7 Magnetization versus field (M vs. μ0H) curves for LBSMNO compounds at 5 K. The inset: M2 vs. μ0H/M plots around TC. |
In the inset of Fig. 7(a–c), we have presented the Arrott plots (M2 vs. μ0H/M) for the La0.6Ba0.2Sr0.2Mn1−xNixO3 (0 ≤ x ≤ 0.1) compounds, from these curves we can conclude the nature of the magnetic phase transition. The slope of the curves is positive, so the transition is second order, according to Banerjee’s criteria.52 The MCE is an intrinsic characteristic of magnetic materials.53–55 Its principle is based on cooling or heating compounds when subjected to a magnetic field under adiabatic conditions, which is maximized when materials are close to their magnetic control temperature.
|ΔSM| can be determined, using Maxwell’s equations, by the formula:
(8) |
The peaks of the elaborated samples are the same as those of the experiment in the vicinity of TC where is the experimental value from the M(T) curve at μ0H.
We can use the following expression:
(9) |
The predicted (−ΔSM) vs. temperature plots are shown in Fig. 8 at different μ0H values. It can be seen that (−ΔSM) depends on μ0H and the temperature until a maximum value (−ΔSmaxM) is reached around TC. (−ΔSM) increases with the increase of μ0H for each compound due to the spin effect, which becomes important with the increase of μ0H. The (−ΔSmaxM) for x = 0, 0.05 and 0.1 are 7.65, 5.44 and 4.45 J kg−1 K−1, respectively, at μ0H = 5 T. Although, when increasing the Ni ratio, (−ΔSmaxM) and TC decrease. This can be explicated by the lowering of the Mn3+/Mn4+ ratio, which goes from 1.5 (x = 0) to 0.8 (x = 0.1) and subsequently favors the DE interaction of Mn3+–O–Mn4+ over the superexchange (SE) interaction of Mn4+–O–Mn4+, Mn3+–O–Mn3+ and Ni2+–O–Ni2+.56 This result is strongly affected by structural parameters, such as decreasing the bond angles (Mn/Ni)–O–(Mn/Ni) and increasing the bond distances (Mn/Ni)–O.
Fig. 8 Temperature dependence of the magnetic entropy change under different external fields for LBSMNO manganites. |
(−ΔSmaxM) is not the only factor that determines the applicability of such a material, but also the temperature range over which it remains considerable is significant.
The relative cooling power (RCP) is another very important parameter along with (−delta), which defines the amount of heat that can be released between cold and hot sinks in an ideal refrigeration cycle and can be given by the following formula:57
RCP = −ΔSmaxM × δTFWHM | (10) |
Fig. 9 present the relative cooling power as a function of μ0H for our compounds. The values of RCP increase with increasing μ0H and reache about 214 J kg−1 for x = 0, 230 J kg−1 for x = 0.05 and 285 J kg−1 for x = 0.1 at μ0H = 5 T. To better understand the performance of the MCE of our compounds, the values of (−ΔSmaxM) and RCP are compared to other manganites, as given in Table 4.58–62 It can be noted that our samples, especially for x = 0.1, have a suitable TC value, close to RT, and a relatively large magnetic entropy change to include other materials. This proves that our materials can be used in the field of magnetic refrigeration.
Composition | TC | |ΔSmaxM| (J kg−1 K−1) | RCP (J kg−1) | μ0H (T) | Ref. |
---|---|---|---|---|---|
Gd | 293 | 5 | 153 | 2 | 46 |
Gd | 293 | 9.5 | 410 | 5 | 47 |
Gd5(Sr2Ge2) | 275 | 18.5 | 535 | 5 | 47 |
La0.7Sr0.3Mn0.95Ti0.05O3 | 308 | 2.2 | 90 | 2 | 48 |
La0.6Ba0.2Sr02MnO3 | 354 | 7.65 | 214 | 5 | This work |
La0.6Ba0.2Sr0.2Mn0.95Ni0.05O3 | 320 | 5.44 | 230 | 5 | This work |
La0.6Ba0.2Sr0.2Mn0.9Ni0.1O3 | 301 | 4.45 | 285 | 5 | This work |
La0.6Ba0.2Sr02MnO3 | 354 | 3.88 | 82 | 2 | This work |
La0.6Ba0.2Sr0.2Mn0.95Ni0.05O3 | 320 | 2.69 | 98 | 2 | This work |
La0.6Ba0.2Sr0.2Mn0.9Ni0.1O3 | 301 | 2.37 | 121 | 2 | This work |
La0.7Sr0.3Mn0.9Fe0.1O3 | 260 | 1.7 | 83 | 2 | 49 |
La0.67Sr0.33Mn0.9Cr0.1O3 | 328 | 5 | — | 5 | 50 |
On the other hand, to affirm the nature of phase transition, Franco et al.63 proposed a phenomenological universal curve for the field dependence of ΔSM. We can construct the universal curve by normalizing all the ΔSmaxM: ΔSM(T, μ0H)/ΔSmaxM below and above TC, by imposing that the position of two additional reference points in the curve correspond to θ = ±1.
(11) |
By referring to Banerjee’s criteria for a 2nd order PM–FM transition, all ΔSM curves at different μ0H values should merge into one curve with temperature scaling. If not, the samples follow a 1st order phase transition. Fig. 10 represents the evolution of ΔSM(T, μ0H)/ΔSmaxM vs. temperature θ at different μ0H values for x = 0.05, for example. In this figure, all the data collapses into a single master curve around TC, indicating the 2nd order nature of this phase transition. These results are in good accordance with those obtained by the Banerjee criterion discussed before. In addition, this universal curve can be used for practical purposes, such as extrapolating results to fields or temperatures not available in the laboratory, improving data resolution and deconvoluting the response of overlapping magnetic transitions.64
Fig. 10 Universal behavior of the scaled ΔSM curves of the LBSM0.95N0.05O sample under various fields. |
We can adjust this curve by the Lorentzian function:
(12) |
• α = 1.85 ± 0.03, β = 1.55 ± 0.02 and γ = 0.53 ± 0.04, for T ≤ TC.
• α = 1.16 ± 0.01, β = 1.12 ± 0.02 and γ = − 0.07 ± 0.03, for T > TC.
From eqn (7), the stance and magnitude of the peak, namely, (ΔSmaxM, TC) and Tr1 and Tr2 are the only ones that are to describe ΔS, where Tr2 > TC and Tr1 < TC. In the end, to transpose ΔS(θ) into the real ΔSM(T), we only use these values, which are fixed by the properties of the compounds.
It is essential to study how the MCE evolves over the ranges of applied magnetic fields and the desired temperatures, taking into account that the 2nd order transition has been proven for all samples.
The evolution of the ΔSM vs. the field is given by the expression, according to Parker and Oesterreicher:65 ΔSmaxM = b(μ0H)n, b is a constant and n is an exponent, which depends on the magnetic state of the sample. The exponent n can be expressed by the following expression:
(13) |
According to the mean field approach for conventional ferromagnetic compounds, a minimum value of “n” is 2/3 at TC. Below TC, n has been predicted to be 1 and the materials are in the FM state. However, above TC, it is equal to 2 in the PM zone, according to the Curie–Weiss law. Yet, recent experimental data indicates a deviation from n = 0.66, in the case of a few soft magnetic amorphous compounds. The temperature dependence of n is shown in Fig. 11. The values of n are found to be 0.67, 0.45 and 0.32 for x = 0.00, 0.05 and 0.10, respectively. For x = 0, the value is close to the values of the mean field model. However, for x = 0.05 and 0.1, these values do not coincide with the predicted value of the mean field of 0.66. This difference is probably due to local inhomogeneities around TC.66
x = 0.00 | x = 0.05 | x = 0.10 | |
---|---|---|---|
ρ0 (Ω cm) | 0.0011 | 0.0015 | 0.0022 |
ρ2 (× 10−7 Ω cm K−2) | 74.39 | 94.60 | 1.13 |
ρ5 (× 10−15 Ω cm K−5) | 72.73 | 1.28 | 2.32 |
A (× 10−7 Ω cm) | 6.47 | 4.65 | 3.45 |
Ea/kB (K) | 2100 | 1390 | 780 |
U0/kB (K) | 8224 | 7822 | 6993 |
TmodC (K) | 305 | 297 | 275 |
R2 | 0.9998 | 0.9995 | 0.9999 |
To better understand the contribution of the different factors causing the conduction mechanism below the transition temperature (T < TM–SC), the ρ(T) curve was fitted using different theoretical models.
Conduction electrons meet different competitors, including scattering of the grain/domain boundary, electron–magnon scattering and electron–electron scattering. Using the following empirical relation,67 the electrical resistivity data is analyzed:
ρFM = ρ0 + ρ2T2 + ρ5T5 | (14) |
The values of the electrical resistivity were adjusted using eqn (14) and in the inset of Fig. 13, we have given the best fit. In Table 5 we have grouped together the estimated values of the adjusted parameters.
Above (T > TM–SC), the resistivity is simulated by the SPH mechanism.69
ρ is expressed as in the adiabatic SPH model:
(15) |
To understand the transport mechanism of the total resistivity over the whole temperature range, we used a phenomenological percolation model.70 For this model, resistivity is defined based on the contributions of the FM clusters in the PM region. Thus, the resistivity is expressed by the following expression:
ρ = ρFMf + ρPM(1 − f) | (16) |
The volume fraction follows the Boltzmann distribution and this is expressed via the following equation:
With T = TmodC, f = fc = 0.5 where fc is called a percolation threshold.71 Where f < fc, the sample remains semiconducting and for f > fc it acquires a metallic phase.72
Hence, in the entire temperature range eqn (16), can be given as:
(17) |
The data evaluated from eqn (17) are in agreement with the experimental results. It can be seen that the percolation model adequately describes the resistivity behavior over a wide range of temperatures, including the phase transition region. We have grouped the most suitable parameters in Table 5 and in Fig. 13, we have presented the fit of the data.
Fig. 13 Temperature dependence of resistivity fitted according to eqn (16). The inset shows f(T) as function of temperature. |
The inset in Fig. 13 shows f(T) vs. temperature for all samples. Below TM–SC, the volume concentration of the ferromagnetic phase remains equal to 1. This proves the strong dominance of the FM part in this zone. Afterwards, f(T) starts to lower to zero, since the metallic state (FM) moves to a semiconductor state (PM). This result proves the validity of the percolation approach.
This journal is © The Royal Society of Chemistry 2022 |