Ali Mleiki1,2, R Hanen1, H Rahmouni1, N Guermazi3, K Khirouni4, E K Hlil5, A Cheikhrouhou2. 1. Unité de Recherche Matériaux Avancés et Nanotechnologies (URMAN), Institut Supérieur des Sciences Appliquées et de Technologie de Kasserine, Kairouan University BP 471 Kasserine 1200 Tunisia alimlk1986@hotmail.com. 2. LT2S, Digital Research Center of Sfax Sfax Technoparc 3021 Sfax Tunisia. 3. Laboratoire Génie des Matériaux et Environnement (LGME), Ecole Nationale d'Ingénieurs de Sfax (ENIS), Université de Sfax Sfax Tunisia. 4. Laboratoire de Physique des Matériaux et des Nanomatériaux Appliquée à l'Environnement, Faculté des Sciences de Gabès cité Erriadh, Université de Gabès 6079 Gabès Tunisia. 5. Institut Néel, CNRS et Université Joseph Fourier BP 166 F-38042 Grenoble Cedex 9 France.
Manganese perovskites with the general formula AMnO3, where A = lanthanide or divalent alkaline earth, are of considerable technological importance, particularly with regard to electrical and magnetic properties. The physical properties of these materials can be changed using magnetic or electric fields, which opens up the field to important technological applications. The chemical and physical properties of perovskite oxides vary with composition from dielectric to superconducting, and they find application as magnetic sensors, spintronic devices, catalysts, magnetic recording devices, and magnetoresistive transducers and in magnetic refrigeration systems.[1-4]The structural, electrical, and magnetic properties of these compounds are sensitive to both the manganese valence and internal chemical factors, such as the average cationic radius 〈rA〉 in the A-site and the A-site cationic disorder. Introducing quenched disorder in manganites as different ions at A-site can leads to a change in structure, which modifies the original system completely. A-site doping can be strongly influenced by the average ionic radius, which exhibits a close relation between the distortion of MnO6 octahedra and the narrowing of the electronic bandwidth.The low bandwidth manganites Pr1−CaMnO3 have been the subject of several studies with calcium concentration due to their diverse physical properties such as charge-ordering (CO),[5] metamagnetic transition,[6] colossal magnetoresistance[7,8] and large magnetocaloric effects.[9] In the case of a 35% Ca-doped sample,[10] it undergoes a charge-ordering (CO) transition at about 230 K, an antiferromagnetic (AFM) transition in the vicinity of 160 K. Without applied magnetic field, its electrical conduction shows insulator characteristics. Below the AFM transition temperature, application of a high magnetic field (above μ0H = 4 T) leads this one to metallic and ferromagnetic (FM) state. In addition, it is worth mentioning for the Pr0.65(Ca1−Sr)0.35MnO3 compounds the partial substitution of Sr2+ ions for Ca2+ ions, having higher ionic radius, shows up many physical interesting phenomena,[11-13] particularly with x ≈ 0.3 the charge and orbital ordering appear below 200 K, but above x ≈ 0.3 collapse into the metallic state has been observed below 100 K.[14]Study of A-site doped Pr0.65Ca0.35MnO3 sample would offer a complementary understanding of the structural, magnetic and electrical properties. To the best of our knowledge, there are no reports on the physical properties of Ba doped Pr0.65Ca0.35MnO3 compounds. Substitution of calcium by barium is interesting as the difference between ionic sizes of Ca and Ba is significantly large. The selection of Ba as a dopant was inspired by an intention to induce chemical pressure inside the lattice to favor the electron/hole transfer between the neighboring Mn sites. In this context, we report how a direct substitution of barium ion in the A-site affects the physical properties as compared to pure Pr0.65Ca0.35MnO3. Our results show that at x = 0.10 doping level the charge ordering (CO) is destroyed and drastic changes of properties have been observed.
Experimental details
Polycrystalline perovskite Pr0.65Ca0.25Ba0.1MnO3 sample was fabricated by standard solid-state reaction method. After thoroughly mixing the stoichiometric reactants of high purity (99.99%) [Aldrich] Pr6O11, CaCO3, BaCO3 and MnO2, the sample is fired at 900 °C for 12 h. Then, the mixture was pressed into pellets and sintered at 1000–1200 K for 48 h in air with several intermediate grindings. The powder X-ray diffractometer (“Panalytical X pert Pro” diffractometer with Cu-Kα radiation (λ = 1.5406 Å)) was employed to identify the phase purity and crystalline structure at room temperature. The unit cell parameters were refined using the Rietveld analysis of the X-ray powder diffraction data with the help of FULLPROF software.[15] The static magnetic measurements were carried out in magnetic fields up to 5 T at temperatures between 5 and 285 K on a BS2 magnetometer equipped with a super conducting coil and developed at Neel Institute. On the both side of the pellet, we deposited a thin aluminum film through a circular mask of 6 mm of diameter. Then, we obtained a configuration of a plate capacitor used to measure both the capacitance and the conductance. The sample is mounted in a cryostat which allows the variation of temperature from 80 to 300 K. Conductance and capacitance measurement are performed using an Agilent 4294A analyzer.
Results and discussion
X-ray diffraction
In order to determine the crystalline structure of Pr0.65Ca0.25Ba0.1MnO3, room temperature powder X-ray diffraction (XRD) pattern was collected. The XRD pattern of Pr0.65Ca0.25Ba0.1MnO3 sample is shown in Fig. 1. It is seen from the results that all peaks can be indexed according to a perovskite structure with orthorhombic Pnma symmetry. Rietveld refinement confirms that there are no detectable impurities in this compound. From the XRD pattern lattice parameters obtained from refinements are: 5.4565(3) Å, b = 7.7240(4) Å and c = 5.4519(5) Å.
Fig. 1
Rietveld refined powder XRD patterns of the Pr0.65Ca0.25Ba0.1MnO3 compound.
Magnetic and magnetocaloric effect
The temperature dependence of the magnetization, M(T), measured via field-cooled (FC) protocol at 0.05 T applied magnetic field is shown in Fig. 2(a). With decreasing temperature the Pr0.65Ca0.25Ba0.1MnO3 compound displays a typical paramagnetic–ferromagnetic (PM–FM) transition. The Curie temperature (TC), corresponding to the minimum in dM/dT − T curve (see inset of Fig. 2(a)), is determined to be 85 K. Apparently, as compared to Pr0.65Ca0.35MnO3,[16] a 10% barium substitution at the A-site destroys the charge ordering (CO) in the pristine compound. Because the larger ionic size of Ba2+ compared to that of Ca2+, the destruction of CO state can be attributed to the increase of the average ionic radius of the A-site (〈rA〉) which leads to the increase of eg bandwidth, favors the electrons movement through the lattice and promotes the FM metallic state. The isoelectronic A-site substitution controls the structure by the distortion of MnO6 octahedra, i.e., tilting and rotation of the MnO6 octahedra, which modifies local structural parameters such as the Mn–O bond distance and Mn–O–Mn bond angle, and then directly influence the case of electron hopping between Mn ions.[17] Several works on manganites show that the FM state as well as the CO are affected due to the chemical pressure at A-site.[18-20]
Fig. 2
(a) Temperature dependence of FC magnetization measured at a field of 0.05 T, the inset shows the derivative of magnetization with respect to temperature (dM/dT vs. T). (b) Inverse susceptibility versus temperature plot for Pr0.65Ca0.25Ba0.1MnO3 bulk.
In Fig. 2(b) we have reported the inverse susceptibility 1/χ as a function of temperature for Pr0.65Ca0.25Ba0.1MnO3. A deviation has been observed which starts around 225 K far above TC, which may be due to growth of short range magnetic interactions in the compound.[17,19] However, the curve in PM regime from 230 to 285 K obeys the Curie–Weiss (CW) formula χ = C/(T − θP), where C is the Curie constant and θP is the PM Curie temperature. The fitting of the linear region to the CW law yields positive θP, θP = 66.5 K, indicating dominant FM interactions. On another hand, the detailed calculation of the theoretical effective moment (μtheoeff) can be found elsewhere.[21] Moreover, the experimental effective moment is found to be μexpeff = 14.5μB, which is much larger than its theoretical counterpart μtheoeff = 5.53μB, confirming clearly the magnetic inhomogeneity of our sample. This behavior can be explained in terms of the presence of FM clusters of Mn3+, Mn4+ double exchange pairs together in the PM region[22]Fig. 3(a) shows the evolution of isothermal magnetization as a function of the magnetic applied field, M(H), for Pr0.65Ca0.25Ba0.1MnO3 sample. At low temperature, we can see that magnetization increases sharply with the applied field for μ0H < 0.5 T indicating the presence of FM behavior. The magnetic nature of the sample above 156 K is clearly PM as the magnetization increases linearly with increasing magnetic field.
Fig. 3
(a) Isothermal magnetization around TC for Pr0.65Ca0.25Ba0.1MnO3 compound at different temperatures. (b) Standard Arrott plots (isotherms M2vs. H/M).
Based on the Arrott plots, the isotherms of M2 as a function of H/M are shown in Fig. 3(b). The order of the magnetic transition was analyzed based in the criterion suggested by Banerjee.[23] As can be observed, the positive slopes of the Arrott plots for Pr0.65Ca0.25Ba0.1MnO3 sample indicate a second-order magnetic transition. According to the mean field theory, the Arrott plots should give a series of parallel straight lines around TC.[24] However, all M2versus H/M curves show non-linear behavior in high-magnetic field having downward curvature which indicate the non-mean-field magnetic interactions in this sample.The magnetocaloric effect (MCE) can be related to the magnetic properties of Pr0.65Ca0.25Ba0.1MnO3. The MCE have been estimated from the temperature and field dependence of the magnetization M(H, T), (see Fig. 3(a)), through the thermodynamics Maxwell's relation:[25]It can be found that the maximum value of magnetic entropy (−ΔSM) increases monotonically with increasing magnetic field change. In addition, we can see in Fig. 4(a) the presence of a transition temperature shift induced by the applied magnetic field. Therefore, the −ΔSM is maximized around TC and its maximum values (−ΔSmaxM) are found to be approximately 2.5 and 4.4 J kg−1 K−1 under a magnetic field changes of 2 and 5 T. Other important parameter of refrigerant materials is the relative cooling power (RCP) which is defined by the product of the peak entropy change and the full width at half maximum of the peak, i.e., RCP = −ΔSmaxMδTFWHM (where δTFWHM is the full width at half-maximum of ΔSM(T)). Thus, the resulting RCP values are found to be 104.2 and 263.5 J kg−1 under a magnetic field changes of 2 and 5 T, respectively. The obtained values are comparable with those of several manganite systems with a large MCE,[26-32] suggesting that our compound may be suitable candidate for magnetic refrigeration. The large MCE in Pr0.65Ca0.25Ba0.1MnO3 can be understood by its sharp magnetization change at the FM–PM transition and its high magnetization saturation.
Fig. 4
(a) Temperature dependence of the magnetic entropy changes ΔSM under various external magnetic fields. (b) Universal curves of Pr0.65Ca0.25Ba0.1MnO3 sample.
Additionally, a phenomenological universal curve for ΔSM(T) in the vicinity of a second-order phase transition has been theoretically established[33,34] and successfully applied to several families of magnetic material,[35-38] which can be used as a further criterion to reveal the order of the magnetic phase transition. In addition, this scaling can be particularly helpful for investigating the refrigerant capacity of a similar material. However, a breakdown of the universal curve can be attributed to inhomogeneity within the material.[39] The temperature axis was transformed according to the relation:where T1 and T2 are the temperature of the two reference points. For present study, T1 and T2 are chosen such that ΔSM(T1)/ΔSmaxM = ΔSM(T2)/ΔSmaxM = 0.5. The rescaled ΔSM(T) curves for different maximum applied magnetic fields ranging from 1 up to 5 T with increments of 1 T are shown in Fig. 4(b). It has been reported that the essential region for collapse to determine the order of the transition is θ < −1 (T < TC). However, the deviation of collapse for first-order phase transition may not be obvious around −1 < θ < 0. The divergence of the curves is clear in the Pr0.65Ca0.25Ba0.1MnO3 compound, particularly far above and below the TC, indicating clearly that this compound does not follows an universal curve and the magnetic phase transition is first order. Thus, we conclude that there is not perfect agreement of the rescaled curves below TC, however a check of the Arrott plots (M2(H/M)) constructions suggested the second order quality of the transition in our sample.
DC-resistivity study
The variation of the dc-resistivity (ρ) as a function of temperature of Pr0.65Ca0.25Ba0.1MnO3 is represented schematically in Fig. 5. Across all the studied temperature ranges, the ρ decreases by increasing temperature which indicates that this compound exhibits a semiconductor behavior. The decrease in ρ when the temperature increases could be due to the increase the density of free carriers that acquire a sufficient thermal energy thermal energy to easily cross the barriers encountered.
Fig. 5
Variation of the electrical resistivity (ρ) versus temperature deduced from the dc conductivity σdc(T), the inset shows the variation of conductivity (σ) as a function of frequency at various temperatures for Pr0.65Ca0.25Ba0.1MnO3.
AC-conductivity study
The inset of Fig. 5 shows the frequency dependence of conductivity (σ) at different temperatures for Pr0.65Ca0.25Ba0.1MnO3. As we can observe, ac-conductivity (σac) spectrum exhibits the presence of several behaviors. At low frequency (f < 104 Hz), a large plateau, corresponds to dc-conductivity (σdc), appears and the σ increases by increasing temperature. Such behavior indicates that conduction process is thermally activated by increasing the density of free carriers in the structure. Above 105 Hz and for T > 180 K, a metallic behavior appears accompanied by a decrease in σ. However, for T < 180 K the σ monotonically increases with increasing frequency, indicating a superposition of different relaxation processes.[40] The range of temperature, where σ monotonically increases, could be described by the law of Jonscher:where A is a constant determining the strength of polarizability and s is the frequency exponent representing the degree of interaction between mobile charge and the lattices around them. To get an information about the conduction mechanism in our sample, we will study the variation of the exponent s depending on the temperature. From the region where the σ is proportional to ω the value of the exponent s is determined from the slope of the linear parts of σ versus frequency. As shown in Fig. 6(a), the exponent s decreases by increasing temperature. The result are in good agreement with a hopping mechanism in disordered material[41] and confirm that conduction process is thermally activated. Such behavior is observed in different oxides and amorphous compounds and is discussed using Pike model.[42] Also such variation is in good agreement with the model of Mott and Davis:[43]where kB is a Boltzmann's constant, T is the absolute temperature and Wm is defined as the energy required for an electron to pass from one site to another. The observed decrease in the frequency exponent s with increasing temperature suggests that the correlated barrier hopping (CBH) model best describes the ac conduction mechanism.[44] In the Pr0.65Ca0.25Ba0.1MnO3 compound, large number of crystal boundaries acts as barriers which are overcome by the charge carriers, resulting in conduction across the samples. According to Funke[45] in Joncher's equation, if s ≤ 1, it indicates that the hopping motion due to translational motion.
Fig. 6
(a) Temperature dependence of the frequency exponent s, the inset shows the plot of ‘τ’ versus 1000/T. (b) Arrhenius plot of (σacT) at different frequencies for Pr0.65Ca0.25Ba0.1MnO3.
Fig. 6(b) shows the product ‘conductivity temperature’ (σT) versus (1000/T) at different frequencies (1 kHz, 68 kHz, 200 kHz and 688 kHz). The high temperature conductance data were well fitted to the following equation:[43]where σ0 is a pre-exponential factor, T is the absolute temperature, Eac represents the activation energy of conduction mechanism and kB is the Boltzmann constant. It is clear, at high temperatures, that σ is follows this equation, suggesting that the electronic conduction is dominated by thermally activated hopping of small polaron (SPH).[46] From this figure, it is found that activation energy decreases when increasing frequency (see the inset of Fig. 6(b)). This can be explained by the fact that the charge carriers choose the easiest path between sites.[47] Therefore, the hops become only between closer sites and the activation energy decreases with frequency.
Complex impedance analysis
The complex impedance Z* of the network is defined as:where Z′ is the real part and Z′′ the imaginary part of the complex impedance. Fig. 7(a) displays the evolution of Z′ with frequency at different temperatures for Pr0.65Ca0.25Ba0.1MnO3. At low frequencies, the spectrum of Z′ has higher values. Such behavior is widely observed in perovskite systems.[48,49] At higher frequencies, however, the value of Z′ merges for all temperatures. Such observation indicates the presence of space charge polarization,[50] which is released at high temperatures. At low frequencies and when temperature increases, Z′ decreases. This behavior can be related mainly to the lowering in the density of trapped charges and then a rise in the mobility of charges carriers.[51] The behavior of Z′ observed for this compound at lower and higher frequencies is in good agreement with the reported results in the literature.[52,53] The variation of Z′′ as a function of frequency and temperature is shown in Fig. 7(b). This spectrum is characterized by the presence of a peak at a specific frequency, which is conventionally known as the “relaxation frequency” (fr). When the temperature increases, the relaxation frequency shifts to the higher frequency region (see inset of Fig. 7(b)). This evolution indicated the presence of a relaxation phenomenon in the studied compound.[54] In addition, it is noticed that the center of each peak, i.e., at the maximum of Z′′, does not lie at the identical frequency indicating a deviation from the typical Debye's model in the material.[55] Also, it indicates that the peak intensity of Z′′ decreases with temperature, confirming the semiconductor in this compound.[48] By determining the value of the relaxation frequency, which determined for the low frequency peak centered at fr, and using the following relation: 2πfrτ = 1, the relaxation time (τ) can be calculated. Using Arrhenius model, the inset of Fig. 6(a) exhibits the variation of log τ versus 1000/T. It is observed that the value of τ decreases by increasing the temperature, which indicate a thermally activated process. This result shows the presence of one conduction mechanism where the slope does not changes across all the studied temperature ranges.
Fig. 7
(a) Frequency dependence of the real part (Z′) of the impedance at different temperatures. (b) Frequency dependence of the imaginary part (Z′′) of the impedance at different temperatures for Pr0.65Ca0.25Ba0.1MnO3, The inset shows zoom-in of the high-temperature region.
Fig. 8 shows the Nyquist diagram of the system Pr0.65Ca0.25Ba0.1MnO3 obtained at T = 80 K. The spectrum is characterized by the appearance of a semicircle arc. To connect the microstructure of the compound with its electrical response, an equivalent electrical circuit can be established by modelling the impedance data using Zview software. As can be seen in Fig. 8, the agreement between the experimental and theoretical data is satisfactory. Then, the equivalent configuration is of the type ((R1//ZCPE3) + (Rs//ZCPE2) + (Rp//ZCPE1)) as shown in the inset of Fig. 8.
Fig. 8
Complex impedance spectrum (Nyquist plot) at T = 80 K with the electrical equivalent circuit.
Dielectric properties
The study of the dielectric properties is another important source of information regarding the conduction process. The complex dielectric permittivity according to the Mott's theory[43] is defined as:where ε′ and ε′′ are the real and imaginary parts of ε*, respectively. ε′ describes the stored energy and ε′′ describes the dissipated energy in the dielectric material, due to the presence of an applied electric field. The dissipation factor (D) and the capacitance (C) of the compound are measured over a wide frequency range (40 Hz to 10 MHz). Thus, ε′ and ε′′ can be determined using the following expressions: ε′ = Ct/(ε0A) and ε′′ = ε′D, where ε0 is the free space permittivity, t is the thickness of the sample and A is the area of the electrodes. The evolution of ε′ and ε′′ with frequency in the temperature range 80–160 K for Pr0.65Ca0.25Ba0.1MnO3 are shown in Fig. 9(a and b). As can be seen, at low frequencies, ε′ has higher values which can be related to the presence of space charge produced via the electrode-sample contact and interface between the grain and their boundaries.[56] According to Singh et al.[57] and Kar et al.,[58] such behavior can be attributed to the presence of several parameters of polarizations such as dipolar, atomic, ionic, interfacial and electronic polarization in the compound. Therefore, under the applied alternative electric field all electric dipole moments will be oriented in his direction. Then, all the polarization mechanisms can participate to polarization of the system. At higher frequencies, the decrease in the permittivity ε′ could be explained that the heavy particles cannot follow the applied electric field and then dipole response is weak, inducing a reduction of ionic and dipolar polarization contributions. Such behavior is observed when the charge carrier hopping cannot follow the alternating current beyond a certain frequency, which decreases the probability of electrons reaching the grain boundary and as a result, polarization decreases. However, at high frequencies, ε′ remains important which may be due to a localized accumulation of charges involving the presence of interfacial polarization.[59]
Fig. 9
Variation of real part (a) and imaginary part (b) of permittivity as a function of frequency, (c) variation of loss tangent with frequency at different temperature for Pr0.65Ca0.25Ba0.1MnO3.
Fig. 9(b) clearly shown that ε′′ decreases with increasing the frequency and no peak is observed, confirming that the polarization phenomenon in this sample is governed by a hopping process.[60] This result may be due to space charge polarization and interface effects.[61] At low frequencies, ε′′ has higher values, suggesting the contribution of all polarization effects particularly the interfacial/space polarization due to an inhomogeneous dielectric structure. At high frequencies, where the electric dipoles cannot follow the applied electric field, the low value attained may be due to the gradual decrease of the dielectric loss. On another hand, the variation of tang(δ) (=ε′′/ε′) as a function of frequency is shown in Fig. 9(c). At constant frequency the dielectric loss increases by increasing temperature which may be due to the increase in conductivity. However, by increasing frequency, the loss tangent first decreases and then increases. Since the charge carriers cannot follow the electric field, at higher frequencies, the value of tang(δ) remains low.
Conclusion
Structural, magnetic and electrical properties of Pr0.65Ca0.25Ba0.1MnO3 manganite were investigated. A Rietveld analysis shows that the sample crystallizes in the orthorhombic structure with Pnma space group. Magnetic measurements show a sharp paramagnetic to ferromagnetic transition around Curie temperature TC. A large maximum of magnetic entropy change (−ΔSmaxM) value has been found with a large relative cooling power (RCP) value. The variation of resistivity as a function of the temperature shows the presence of a semiconductor behavior. The ac-conductivity study is obeyed to the Jonscher law. The dominance of hopping mechanism was confirmed by the variation of the exponent s. Complex impedance analysis permits to model the sample by an electrical equivalent circuit and displays the presence of a relaxation phenomenon. The dielectric study shows the contribution of different types of polarization.