M. Jeddi*a,
J. Massoudia,
H. Gharsallahab,
Sameh I. Ahmedc,
E. Dhahria and
E. K. Hlild
aLaboratoire de Physique Appliquée, Faculté des Sciences, Université de Sfax, B. P. 1171, 3000 Sfax, Tunisia. E-mail: marwajeddi@gmail.com
bInstitut Préparatoire aux Études d'Ingénieur de Sfax, Université de Sfax, B. P. 1172, 3018 Sfax, Tunisia
cDepartment of Physics, College of Science, Taif University, P. O. Box 11099, Taif 21944, Saudi Arabia
dInstitut Néel, CNRS, Université J. Fourier, B. P. 166, 38042 Grenoble, France
First published on 11th February 2021
Nd0.6Sr0.4MnO3 polycrystalline manganite was synthesized by two different methods: the auto-combustion reaction (NSMO-AC) and the sol–gel method (NSMO-SG). The structural, magnetic, magnetocaloric and critical behavior of the samples were examined. Rietveld refinements of the XRD patterns revealed that both compounds are pure single phase indexed to the orthorhombic system adopting the Pnma space group. The nanometric size estimated using the Williamson–Hall method was confirmed by TEM micrographs. Magnetic measurements as a function of temperature indicated that both samples underwent a second order ferromagnetic (FM)–paramagnetic (PM) phase transition at Curie temperature (TC). The relative cooling power was observed to be around 95.271 J kg−1 for NSMO-AC and 202.054 J kg−1 for NSMO-SG at μ0H = 5 T, indicating that these materials are potential candidates for magnetic refrigeration application close to room temperature. The critical behavior was estimated using diverse techniques based on the isothermal magnetization data recorded around the critical temperature TC. The calculated values are fully satisfactory to the requirements of the scaling theory, implying their reliability. The estimated critical exponents matched well with the values anticipated for the mean-field model and the 3D Ising model for NSMO-AC and NSMO-SG, respectively, showing that the magnetic interactions depended on the process of elaboration.
Nowadays, intensive research activities attempt to explore new magnetic materials presenting optimal magnetocaloric properties near room temperature as the most efficient refrigerant materials in the field of magnetic cooling technologies.7–9 Perovskite-based manganites of compositional formula (RE1−xAx)MnO3, where RE is a trivalent rare earth atom and A is a divalent alkaline earth one,10–12 can be considered as good competitors for other magnetocaloric materials. This interest is due to their powerful ability to be exploited in different fields.13–20 There are several robust arguments that emphasize the fact that this class of compounds will fulfill an impressive contribution in the emerging technologies21 given that their critical temperature can be adjusted by substitution; their preparation is relatively easy and their cost is lower than that of gadolinium, which is considered as a reference.22–24 Among these materials, only lanthanum manganites have been fully characterized with a deeper investigation into the structural and magnetic properties.25–27 Neodymium-based manganites seem to be very promising magnetocaloric candidates28,29 owing to the reason that the magnitude of their magnetic entropy change and the value of their transition temperature can be strongly tuned by controlling their microstructure and their chemical composition following the preparation method.30
As the magnetocaloric properties of manganites were examined only to a limited extent, the exploration of other systems and new elaboration process inducing an improvement in the magnetocaloric effect is still under the initial stage.31,32 Different methods, such as solid state reaction,33–36 chemical co-precipitation,37–39 and sol–gel method40–43 are available for the synthesis of RE manganites with a perovskite-like structure. The methods of synthesis play a significant role in the formation of the required crystal structure, the control of the crystallite size, and the contribution of magnetic interactions.
The sol–gel route is a classic soft chemistry process, mainly explored in the context of nanometric manganites, which lead to convincing results in terms of particle size. However, calcination is an important step to induce the good crystallization of materials and to remove foreign impurities. This step determines the final dimension of objects and generates sintering phenomena.44,45 The auto-combustion reaction cannot be described as a soft chemistry process as the temperatures reached during autoignition are very high. However, the fact that the transition from the liquid state to the solid state is very fast allows access to crystallized and single-phase materials with very small sizes, which eliminate the need for a subsequent annealing step.46,47 Therefore, this phenomenon would constitute a major advantage as far as the sol–gel process is concerned. The synthesis by auto-combustion is a promising technique because it is fast and requires simple equipment and low processing time. It favors high product crystallinity, homogeneous chemical composition, small particle size and narrow particle size dispersion.48 The reagents are inexpensive, the dopants can be easily integrated into the final product, and multiple steps are not involved.
Motivated by these considerations, a detailed investigation of the structural, magnetic, magnetocaloric, and critical behavior of strontium-doped neodymium manganite Nd0.6Sr0.4MnO3 synthesized by the auto-combustion process and the sol–gel method is reported.
Fig. 2 (a) Rietveld refinement of the X-ray diffraction patterns, (b) crystal structure and (c) Williamson–Hall plot of NSMO-AC and NSMO-SG compounds. |
Compound | NSMO-AC | NSMO-SG |
Space group | Pnma | Pnma |
a (Å) | 5.44243 | 5.44086 |
b (Å) | 7.66322 | 7.68139 |
c (Å) | 546654 | 5.46634 |
V (Å3) | 227.990 | 228.457 |
〈dMn–O〉 (Å) | 1.96 | 1.95 |
〈θMn–O–Mn〉 (°) | 160.30 | 161.08 |
DWS (nm) | 18 | 25 |
DTEM (nm) | 32 | 61 |
χ2 (%) | 1.68 | 1.34 |
The average crystallite size was determined according to the Williamson–Hall approach using the following relation:53–55
(1) |
Fig. 2c represents the (βcosθ) (y-axis) vs. (4sinθ) (x-axis) data related to the most intense peaks of NSMO-AC and NSMO-SG samples. The crystallite size DWH was deduced from the y-intercept of the linearly fitted data, and the microstrain ε is deduced from the slope. The DWH and ε values were, respectively, equal to 18 nm and 0.0008 for NSMO-AC and 25 nm and 0.0003 for NSMO-SG.
The surface TEM micrographs of the aforementioned compounds are exemplified in the inset of Fig. 3a. It is worth noticing that the shape of the particles differed based on the synthesis route. The size distribution of the grains analyzed quantitatively by fitting the histogram according to a Lorentzian function is depicted in Fig. 3a. The mean diameters of NSMO-AC and NSMO-SG compounds were around 32 and 61 nm, respectively. The particle size determined from the TEM image was higher than that estimated by the XRD experiment, which can be elucidated by the fact that each particle analyzed by TEM was composed of numerous different crystallized grains.
Fig. 3 (a) Size distribution histogram of NSMO-AC and NSMO-SG compounds. The inset presents the TEM image, (b) EDAX analysis spectrum. |
Fig. 3b exhibits the EDAX spectra of NSMO-AC and NSMO-SG compounds, respectively. The analysis was performed on several zones. All the ingredients incorporated within the synthesis (Nd, Sr, Mn, and O) were present. One can see that there are no impurities. The composition of the compound was identical to the desired one.
Fig. 4 Magnetization (M) vs. temperature and the plot of dM/dT at μ0H = 0.05 T of NSMO-AC and NSMO-SG compounds. |
Fig. 5 presents the temperature dependence of the inverse magnetic susceptibility χ−1(T) data. The solid blue line was the best fit of Curie Weiss law in the PM region defined as follow:56,57
(2) |
For the NSMO-AC sample (Fig. 5a), the χ−1(T) curve showed a downturn deviation from the Curie Weiss law indicating the existence of a Griffiths phase (GP),58–60 which corresponded to the appearance of the FM clusters within the PM region in the temperature range TRand ≤ T ≤ TGP, thereby defining the Griffiths regime. TGP was the temperature at which the χ−1(T) curve deviated from the Curie Weiss law called the Griffiths temperature.61 This was also illustrated by the dM/dT plot (Fig. 4). TRand was the critical temperature of a random FM state called the random transition temperature.62 TRand was equivalent to the Curie Weiss temperature (θCW).
The GP was characterized by an exponent λ, defined by the following equation:
χ−1(T) ∝ (T − TRand)1−λ; (0 ≤ λ < 1) | (3) |
The double logarithmic plot of (χ−1) vs. tm, tm = T/TRand − 1, reproduced in Fig. 5a2, shows a linear behavior at low tm values. This result confirms the proposed GP.64
The observed slope change, as marked in Fig. 5a, is ascribed to the presence of a charge ordering state (CO).65,66 As reported in a previous study,67 the CO transition refers basically to the appearance of an antiferromagnetic (AFM) ordering phase which makes the inverse susceptibility curves plotted upward. TCO is the CO transition temperature (Table 2). It could also be determined from the conspicuous kink in the dM/dT plot (Fig. 4).
However, for the NSMO-SG sample (Fig. 5b), there was no obvious observable downturn or upward deviation. The χ−1(T) followed the Curie Weiss law. The obtained θCW value was positive (Table 2), proving the FM behavior of the studied sample.68,69
For the NSMO-SG compound, the disappearance of the GP and the CO state was attributed to the enhancement of FM interactions. Each particle was supposed to be composed of two parts.70 The first one was the core where the double exchange interactions dominated and promoted an FM character. As for the other part, it was a disordered shell where the magnetic interactions were modified by the vacancies, defects as well as broken bonds and disfavored the FM behavior.71 Therefore, the appearance of the GP and the CO state until the reduction of the crystallite size arose from the loss of the FM order, since the surface contribution was larger.
Isothermal magnetizations versus the applied magnetic field M(μ0H,T) for NSMO-AC and NSMO-SG samples performed at several temperatures are plotted in Fig. 6. At low temperature values, the M(μ0H,T) data of NSMO-AC increased rapidly in the low external magnetic field region and were not as yet saturated for an applied magnetic field as great as 5 T. Contrary to the NSMO-SG compound, the M(μ0H,T) data increased sharply in the low magnetic field region and subsequently saturated as the field value increased, which corroborated to the high FM character of the NSMO-SG sample compared to that of the NSMO-AC sample. At high temperature values, the M(μ0H,T) data changed linearly with the applied magnetic field. Both compounds showed a PM behavior.
Fig. 6 Isothermal magnetization curves recorded at various temperature values around TC for NSMO-AC and NSMO-SG compounds. |
(4) |
The temperature dependence of the magnetic entropy change (−ΔSM(T)) of NSMO-AC and NSMO-SG compounds, established at various external magnetic fields, is shown in Fig. 7. It was worth noting that the magnitude of ΔSM increased gradually when the external magnetic field was increased and achieved its maximum around TC, which was indicative of a much considerable magnetic entropy change to be predicted at higher magnetic field values. For μ0H = 2 T, the NSMO-AC sample presented a maximum ΔSmaxM value of 0.247 J kg−1 K−1, and the NSMO-SG sample reached a maximum ΔSmaxM value of 1.689 J kg−1 K−1. For μ0H = 5 T, the maximum ΔSmaxM value was 0.713 J kg−1 K−1 for the NSMO-AC compound and 3.594 J kg−1 K−1 for the NSMO-SG compound. Thus, the NSMO-SG compound achieved an important MCE near room temperature with a considerable ΔSmaxM value making the system potentially beneficial for magnetic refrigeration applications.
Fig. 7 Magnetic entropy change (ΔSM) vs. temperature at several applied magnetic fields for NSMO-AC and NSMO-SG compounds. |
As can be seen, the ΔSM peak of the NSMO-AC sample was larger than that of the NSMO-SG sample. Such an enlargement of the ΔSM peak was assigned to the increase in the surface/volume ratio. When the crystallite size reduced, the surface proportion increased, which resulted in a weakness of the FM coupling compared to that of the core. This fact led to a scattering of the Curie temperature values, and thus a broad magnetic transition.73,74
The magnetocaloric effectiveness could be assessed through the relative cooling power (RCP),75 which was defined as the quantity of heat transferred between the hot and the cold sinks in an ideal refrigeration cycle. The latter takes into account the magnitude of (ΔSM) and its FWHM (δTFWHM). It was calculated as follow:
RCP = (−ΔSmaxM) × δTFWHM | (5) |
The obtained RCP was 95.271 J kg−1 for NSMO-AC and 202.054 J kg−1 for NSMO-SG at μ0H = 5 T, which accounted for about 23% and 50% of that estimated in pure gadolinium (Gd), respectively.
To further estimate the efficacy of our compounds in the magnetic cooling systems, the calculated ΔSmaxM and the RCP values found in the present study are compared with that of other magnetic materials26,29,30,76–79 and tabulated in Table 3.
Compound | μ0H (T) | TC (K) | −ΔSmaxM (J kg−1 K−1) | RCP (J kg−1) | Ref. |
---|---|---|---|---|---|
Nd0.6Sr0.4MnO3 (auto-combustion method) | 1 | 100 | 0.083 | 7.940 | This work |
2 | 0.247 | 27.481 | |||
3 | 0.408 | 49.929 | |||
4 | 0.56 | 70.632 | |||
5 | 0.713 | 95.271 | |||
Nd0.6Sr0.4MnO3 (sol–gel method) | 1 | 275 | 0.853 | 28.322 | This work |
2 | 1.689 | 72.105 | |||
3 | 2.417 | 116.637 | |||
4 | 3.033 | 163.900 | |||
5 | 3.594 | 202.054 | |||
Gd | 2 | — | 5.5 | 164 | 76 |
5 | 10.2 | 410 | 77 | ||
Nd0.6Sr0.4MnO3 (solid-state method) | 2 | 245 | 1.71 | 74.88 | 30 |
5 | 3.68 | 216.03 | |||
La0.6Sr0.4MnO3 | 2 | 365 | 0.97 | 98 | 26 |
5 | 2.14 | 264 | |||
Nd0.50-0.17Sr0.33MnO3 | 1 | 220 | 0.19 | — | 29 |
Nd0.7Sr0.3MnO3 | 2 | 240 | 5.3 | 80 | 79 |
Nd0.67Ba0.33MnO3 | 5 | 145 | 3.91 | 265 | 78 |
Nd0.67Ba0.33Mn0.98Fe0.1O3 | 5 | 134 | 2.97 | 242 | 78 |
Fig. 8 Modified Arrott plots (MAP): M1/β vs. (μ0H/M)1/γ with (a) mean-field model, (b) tri-critical mean-field model, (c) 3D-Heisenberg model and (d) 3D-Ising model for NSMO-AC and NSMO-SG samples. |
Based on the universal scaling relations, a second-order magnetic transition near the Curie temperature TC was specified by means of interconnected critical exponents, namely β is the exponent related to the spontaneous magnetization MS under TC, γ is the exponent associated with the magnetic susceptibility χ0−1 beyond TC and δ is the exponent linked to the critical isothermal magnetization at TC. The critical exponent values were mathematically approximated from the magnetization data via the asymptotic formulas as follow:82
MS(T < TC,μ0H → 0) = M0|ε|β | (6) |
(7) |
M(T = TC,μ0H) = D(μ0H)1/δ | (8) |
According to the Arrott–Noakes equation of state,83 the experimental data were examined through the modified Arrott-plot (MAP) formula:
(9) |
The appropriate model brought about a set of reasonably good parallel straight lines.
Fig. 8 displays the plot of M1/β vs. at various temperatures using different theoretical models: (a) the mean-field model (β = 0.5 and γ = 1), (b) the tri-critical mean-field model (β = 0.25 and γ = 1), (c) the 3D-Heisenberg model (β = 0.365 and γ = 1.336) and (d) the 3D-Ising model (β = 0.325 and γ = 1.240). All models presented nearly parallel lines and quasi straight in the region of high field values. It was complicated to identify the most adequate model to characterize our samples.
Herein, it was essential to calculate the relative slope (RS), which was considered as a new specific indicator to select the suitable model. The RS was expressed at the critical temperature as follow:
(10) |
Fig. 9 presents the RS vs. T data of both the systems for the following four models: the mean-field model, the tri-critical mean-field model, the 3D-Heisenberg model and the 3D-Ising model. The most ideal model was the one that presented an RS value equal to 1.84
It was clear that the mean-field model was the most satisfactory one to define the critical properties of the NSMO-AC compound, whereas the 3D Ising model was the ideal one to specify the critical behavior of the NSMO-SG compound. The change in the universality class could be associated with the change in the elaboration process, which influenced the crystallite size. Similar results were reported in ref. 85 and 86, where the authors suggested that the reduction in the crystallite size presented an intensive impact on the magnetic interactions and strongly influenced the universality class.
Departing from the MAP, the linear extrapolation of the isotherm in the high magnetic field region provided the spontaneous magnetization MS values as well as the inverse susceptibility χ0−1 values as intercepts on the coordinate axes M1/β and , respectively. The MS(T) and χ0−1(T) data are demonstrated in Fig. 10. Using eqn (6) and (7), the best fit results of MS(T) and χ0−1(T) plots, respectively, provide new values of β, γ, and TC (Table 4).
Fig. 10 Spontaneous magnetization MS(T) and the inverse initial susceptibility χ0−1(T) data for NSMO-AC and NSMO-SG compounds. Solid lines correspond to the fitting curves. |
Model/compound | Technique | TC (K) | β | γ | δ |
---|---|---|---|---|---|
Mean-field model | 0.5 | 1 | 3 | ||
Tri-critical mean-field model | 0.25 | 1 | 5 | ||
3D-Heisenberg model | 0.365 | 1.336 | 4.80 | ||
3D-Ising model | 0.325 | 1.240 | 4.82 | ||
NSMO-AC | MAP | 105.195 ± 0.006 | 0.461 ± 0.001 | 1.209 ± 0.003 | |
KF | 104.167 ± 0.008 | 0.621 ± 0.002 | 1.123 ± 0.005 | ||
CIA (exp.) | 2.516 ± 0.038 | ||||
CIA (cal.) | 2.80 | ||||
NSMO-SG | MAP | 274.150 ± 0.006 | 0.320 ± 0.001 | 1.239 ± 0.002 | |
KF | 273.980 ± 0.001 | 0.323 ± 0.002 | 1.235 ± 0.004 | ||
CIA (exp.) | 4.831 ± 0.016 | ||||
CIA (cal.) | 4.82 |
From Table 4, it is worth highlighting that the obtained critical exponents of NSMO-SG are self-consistent and in accurate estimates when compared with the theoretical ones. However, the estimated critical exponents of NSMO-AC exhibited a slight deviation from the exponent values predicted for the theoretical interaction models. Generally, for the inhomogeneous FM materials, the values of the critical exponents were incompatible when compared with those of the traditional universality classes. Similar results were reported by Sudakshina et al.87
A further processing of MS(T) and χ0−1(T) was realized using the Kouvel–Fisher (KF) method88 so as to determine accurately the values of β, γ, and TC. The function construction was defined by the following expressions:
(11) |
(12) |
Under this method, the plots of and , shown in Fig. 11, appear as straight lines with 1/β and 1/γ slopes, respectively. The intercepts on the T axis of their extrapolation to the ordinate equal to zero corresponding to the TC values. It was noted that the critical exponent as well as the Curie temperature values determined from the KF method (Table 4) go in tandem with those deduced from the MAP. Consequently, we could conclude that the aforementioned methods were feasible and effective to investigate the critical properties.
Fig. 11 Kouvel–Fisher (KF) curves for the spontaneous magnetization MS(T) and the inverse initial susceptibility χ0−1(T) of NSMO-AC and NSMO-SG compounds. Solid lines correspond to the linear fits. |
Fig. 12 presents the critical isotherm (M vs. μ0H) curves at TC = 100 and 275 K for NSMO-AC and NSMO-SG, respectively, plotted on a log–log scale. By fitting the M(μ0H) with eqn (8), the third exponent δ value is achieved. The obtained values, recapitulated in Table 4, are similar to those calculated theoretically through the Widom scaling law:89
(13) |
Fig. 12 Critical isotherm (M vs. μ0H), plotted on a log–log scale, for NSMO-AC and NSMO-SG compounds. |
These findings validated the exactness of the calculated values of β and γ.
For more certainity, the reliability of the acquired critical exponents could be determined through the universal scaling theory around the critical region, which could be described by the equation below:
(14) |
Fig. 13 displays the plot of M|ε|−β vs. μ0H|ε|−β−γ by using β, γ, and TC values determined by the KF method. It can be obviously seen that all of the experimental data fell into two independent branches around TC. This observation ensured that eqn (14) is obeyed, which denotes the accuracy of the obtained critical exponents and that of the Curie temperature.
Fig. 13 Scaling plots M|ε|−β vs. μ0H|ε|−β−γ around TC for NSMO-AC and NSMO-SG samples. The inset depicts the same plots on the log–log scale. |
This journal is © The Royal Society of Chemistry 2021 |