Mechanical behavior, enhanced dc resistivity, energy band gap and high temperature magnetic properties of Y-substituted Mg-Zn ferrites
M. A. Ali, M. N. I. Khan, M. M. Hossain, F.-U.-Z. Chowdhury, M. N. Hossain, R. Rashid, M. A. Hakim, S. M. Hoque, M. M. Uddin
11 Mechanical behavior, enhanced dc resistivity, energy band gap and high temperature magnetic properties of Y-substituted Mg-Zn ferrites
M. A. Ali , M. N. I. Khan , M. M. Hossain , F.-U.-Z. Chowdhury , M. N. Hossain , R. Rashid , M. A. Hakim ,S. M. Hoque ,and M. M. Uddin Department of Physics, Chittagong University of Engineering and Technology (CUET), Chattogram4349,Bangladesh. Materials Science Division, Atomic Energy Center, Dhaka 1000, Bangladesh. Department of Glass and Ceramic Engineering, Bangladesh University of Engineering and Technology (BUET), Dhaka 1000, Bangladesh.
ABSTRACT
We report the synthesis of Y-substituted Mg-Zn [Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05)] ferrites using conventional standard ceramic technique. The samples were characterized by X-ray diffraction (XRD) analysis, field emission scanning electron microscopy (FESEM), FTIR spectroscopy, UV–Vis spectroscopy and quantum design physical properties measurement system (PPMS).XRD patterns confirm the single phase cubic spinel structure up to x = 0.03 and appearance of a secondary phase of YFeO for higher Y contents. FESEM images depict the distribution of grains and EDS spectra confirmed the absence of any unwanted element. Completion of solid state reaction and formation of spinel structure has been revealed from FTIR spectra. The FTIR data along with lattice constant, bulk density and porosity were further used to calculate the stiffness constant ( C ij ), elastic constant and Debye temperatures. Mechanical stability of all studied compositions is confirmed from C ij using Born stability conditions. Brittleness and isotropic nature are also confirmed using Poisson’s ratio and anisotropy constants, respectively. The enhancement of dc electrical resistivity (10 Ω-cm to 10 Ω-cm) with Y content is observed. The energy band gap (increased with Y contents) is found in good agreement with dc electrical resistivity. Ferrimagnetic to paramagnetic phase change has been observed from the field dependent high temperature magnetization curves. The magnetic moments and saturation magnetization were found to be decreased with increasing temperature. The Curie temperature (T c ) has been measured from temperature dependent magnetic moment (M-T) and initial permeability ( µ ʹ i -T) and found to be in good agreement with each other. Decrease in T c with Y content is due to redistribution of cations and weakening of the exchange coupling constant. The magnetic phase transition has been analyzed by Arrott plot and found to have second order phase transition. The dc resistivity endorses the prepared ferrites are suitable for high frequency and high temperature magnetic device applications as well. Keywords:
Y-substituted Mg-Zn ferrites, dc resistivity, FTIR spectroscopy, mechanical properties, UV-Vis spectroscopy, high temperature hysteresis loop, Curie temperature. Introduction
Since last century, the spinel type ferrites have come out as one of the leading materials in the field of materials science due to their vast applications. The physics involved in spinel structure is still interesting as their properties can easily be tuned, consequently their applications as well. In recent years, new opportunities have been disclosed and numerous directions have been developed for research on these materials by manipulating their physical properties. The long-established applications of the ferrites (in thin film, bulk or nano form) are in transformer applications, solar hydrogen production, multi-layer chip inductors, high density magnetic recording, telecommunications, microwave devices, as contrast agents in magnetic resonance imaging (MRI), sensors, catalysts etc. [1-8]. Moreover, the new application of spinel ferrites has been listed as photocatalysts for photocatalytic degradation of dyes [9-12]. Mg-Zn ferrite is one of the most used soft ferrites because of its high electrical resistivity, low cost, low dielectric loss, high mechanical hardness and superior environmental stability [13]. Mg-Zn ferrites with higher resistivity (10 –10 Ω-cm) make its suitability in high frequency applications [13]. Most of ( ions are located on the B -sites and small fraction ( A -sites [14, 15]. However, the Zn ions have a preference to occupy A -sites in the spinel lattice [16, 17]. The Fe ions are distributed on both A - and B -sites [18]. Therefore, the physical properties of the Mg-Zn ferrites are determined by the cations distribution over the A - and B -sites. The physical properties can be altered by introducing the different metallic ions results the cations distribution modification on the A - and B -sites. The factors that determine cation distributions are ionic radius, charge, site preference and level of substituent ion, methods, conditions and sintering temperature for the preparation of ferrites [19-22]. Study of ferrites usually in most cases has been limited to their structural, electrical and magnetic properties [9-25]. However, study of mechanical properties of ferrites has scientific importance due to their potential applications in industry and in research field. Elastic constants are closely related to many important physical properties of solids, for instance acoustic-phonon vibration and internal stress, Debye temperature. This type of study provides the information about bonding nature within the materials. Also to understand the thermal properties and minimum changes in devices performance with mechanical shock, temperature, vibration, etc. during the manufacturing and testing of ferrite materials [26] Due to the mentioned interest the study of mechanical properties have already been considered by researchers [26-29]. Moreover, the study of magnetic properties of Mg-Zn/other ferrites have also been usually limited in range of room temperature to as low as 4 K [30]. But the most technological uses of these ferrites are either at or above room temperature. The hysteresis loop, magnetization as well as magnetic moment are strongly dependent on operating temperature. The room temperature measurement of magnetic properties cannot provide exact behavior of magnetic material at higher temperature since the temperature dependent magnetic behavior is highly non-linear. Therefore, to design a magnetic device using soft ferrites, thermal effects on magnetic properties need to be considered [31-32]. This prospect also motivates the researchers and some reports on the high temperature magnetic properties also available [30, 33-34]. Furthermore, the study of dc electrical resistivity provides important information regarding suitability of materials in high frequency application. In addition, determination of band gap of semiconducting materials is very important for their practical application. Therefore, the present paper reports the investigation of mechanical properties, high temperature magnetic properties along with FTIR spectra and UV-Vis measurement of Mg-Zn ferrites. Mg-Zn ferrites have been subjected by a large number of research group to modify the different properties using different synthesis routes and conditions [12, 14, 35-49], where some of them are dealt with tailoring electrical properties and/or magnetic properties of Mg-Zn ferrites. The motivation of our present research is to enhance the electrical resistivity of the Mg-Zn ferrites by Y substitution for Fe with favorable magnetic properties that is one of the significant characteristics of ferrites for applications in high frequency devices. Melagiriyappa et al. [35] have been studied the electrical properties of Sm substituted Mg-Zn ferrites where increased of electrical resistivity has been reported. The initial permeability and saturation magnetization (M s ) are focus and the variation is also being explained with the variation of Cu contents as well as density of the samples [36]. Mn substituted Mg-Zn ferrites have been investigated by Mohseni et al. [37] in which the structural, vibrational, micro-structural and magnetic properties are considered. The increase in M s and decrease in coercivity due to Mn substitution is reported. The effect of Cr substitution in the Mg–Zn ferrite has been studied by Haralkar et al. [42] and Xia et al. [43]; they have reported on the decrease of M s with Cr contents. Mukhtar et al. [48] have investigated the Pr substituted Mg-Zn ferrites where the samples are characterized by XRD, FTIR spectroscope, SEM, EDX and VSM measurements. The decrease in M s and increase in coercive field ( H c ) has been reported with increase in Pr contents. The improvement in Magnetic properties of Mg-Zn ferrites by Co substitution has also been reported [45]. Recently, the effect of Y substitution on the electrical resistivity and magnetic properties of different ferrites has been reported [34, 50-56]. But to the best of our knowledge, Y is not considered as a substituent in the Mg-Zn ferrites yet. We have studied the effect of Y substitution on structural properties, morphological study, ac conductivity, dielectric properties, permeability and room temperature magnetic properties of Mg-Zn ferrites and results are available elsewhere [57]. Therefore, in this article, FTIR spectroscopy, mechanical properties, dc resistivity, UV-Vis spectroscopy and temperature dependent magnetic properties of Y-substituted Mg-Zn ferrites prepared from nano-sized raw materials have been presented for the first time. Materials and methods
A series of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05) ferrites, in the steps of 0.01, have been synthesized using analytical grade nano-sized materials by conventional ceramic technique. The use of nano powders as raw materials is inspired by successful completion some other projects that can be found elsewhere [58-61]. The detail of the sample preparation is presented elsewhere [57]. The sample preparation method comprises the following operations as shown in Fig. 1. The X-ray diffraction (XRD) spectra of the prepared samples has been recorded by a PHILIPS X’Pert PRO X-ray diffractometer in the range 2θ = 15–70 using Cu K α radiation (λ = 1.54 Å). The microstructure images have been taken by the Field Emission Scanning Electron Microscope (FESEM) (model: JEOL JSM-7600F) equipped with an energy dispersive spectrometer (EDS). A PerkinElmer FTIR spectrometer has been used to compute the Fourier Transform Infrared (FTIR) spectra using the attenuated total reflection (ATR) technique in the range 400–1000 cm -1 . Elastic properties were calculated using FTIR data along with structural parameters. For dc electrical resistivity measurements, the as-prepared powders were pressed to pellet (diameter 8.4 mm, thickness 2.4 mm) and sintered. Conductive silver paint was applied and two probe method was used to measure the dc electrical resistivity. Temperature dependent M vs H curves and magnetic moment measurement were done by a Quantum Design PPMS VSM magnetometer. Fig. 1.
Flow chart of sample preparation method. Results and discussion
Fig. 2 depicts the XRD patterns of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05) ferrites. The peaks are indexed for different atomic planes and compared with reported peaks of different ferrites [34, 50-56]. The XRD patterns reveal synthesis of spinel ferrites in which the single phase noted up to x =0.03. A secondary phase is appeared at 2 =33.2° for higher Y contents. This extra phase (marked as red circle) results from YFeO because of highly reactive tendency of Fe ions with Y ions [62]. This type of secondary phase is also shown in case of Y substitution in other ferrites [55-56]. The structural parameters have been calculated using the XRD data and can be found elsewhere [57]. Oxide of raw materials (Fe O , MgO, ZnO and Y O ) Weighing by different mole Percentage Dry mixing by Milling, 6 hrs Pressing to desired shape Milling 2.5 hrs Pre-sintering (850 °C, 4 hrs) Sintering (1250 °C, 4 hrs) Finished Product Fig. 2.
The XRD patterns of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05). The physical properties, especially the magnetic properties, are sensitive to the microstructure of ferrites and by manipulation of microstructure the properties can also be tailored. Fig. 3 illustrates the FESEM images of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05) where considerable modification in microstructure is observed due to the Y substitution. The EDS spectra for each composition are presented in Fig. 4(a-f). The qualitative analysis of EDS spectra endorses the existence of Mg, Zn, Fe, Y and O ions and nonattendance of any undesired element. The quantitative analysis of EDS spectra provides the ratio of metal cations and anions present in the samples which are also calculated [Table 1] and found in good agreement with other reported results [14]. Table 1
Cations-anions % (At) and cations-anions ratio of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05). Y contents ( x ) Grain size (µm) cations (At %) anions (At %) Ratio of cation: anion 0.00 1.22 40.03 59.97 2.802:4.197 2.847:4.152 2.820:4.179 2.946:4.053 3.355:3.775 2.688:4.312 0.01 3.65 40.68 59.32 0.02 3.75 40.29 59.71 0.03 1.28 42.09 57.91 0.04 1.07 47.94 53.94 0.05 1.87 38.40 61.60
20 30 40 50 60 70 I n t e n s i t y ( a r b . u . ) x = 0.05 x = 0.04 x = 0.03 x = 0.02 x = 0.01 ( degree ) x = 0.00 (220)(311) (222) (400) (422)(551) (440) Fig. 3 . FESEM images of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05) [57]. (b) x =0.01 (a) x =0.00 (d) x =0.03 (c) x =0.02 (e) x =0.04 (f) x =0.05 FTIR spectroscopy
FTIR spectra of ferrites generally exhibit two major peaks for the M-O vibrations at A- and B - sub-lattice [63]. Fig. 5 illustrates the FTIR spectra of the considered ferrite compositions in which two major peaks are noted within the expected range of frequency: 365 to 1100 cm -1 . Table 2 shows the peaks points for the studied compositions, the completion of solid state reaction has been confirmed from the position of the band within the expected range of frequency for ferrites with spinel structure [64]. The bands at higher frequency (ν = 560 cm -1 ) Fig. 4 . EDS spectra of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05). x =0.00 x =0.01 x =0.02 x =0.03 x =0.04 x =0.05 correspond to the stretching vibrations of the Metal–Oxygen bunch in the A -sites and at lower frequency (ν = 396-402 cm -1 ) are allocated for the stretching vibration of the Metal–Oxygen bunch in the B -sites, respectively [65] . The band frequency ν remains constant for all composition. The band frequency ν slightly shifts to lower frequency side for x = 0.04 and higher frequency side at x = 0.05, this shift in peak position is might be related to the appearance of impurity phase. These measurements (XRD, FESEM and FTIR) confirmed the successful incorporation of Y into Mg-Zn ferrites. Table 2
FTIR absorption bands for Mg Zn Y x Fe x O . Y ( x ) ν (cm -1 ) ν (cm -1 ) FcT×10 (dynes/cm) FcO×10 (dynes/cm) K av ×10 (dynes/cm)
560 400 2.29 1.17 1.73
560 400 2.29 1.17 1.73
560 400 2.29
560 400 2.29
560 396 2.29
560 402 2.29
The stability and bonding characteristics among the adjacent atomic plane can be understood from the stiffness constants. For cubic crystal system, there are three stiffness constants: C , C and C . The FTIR spectroscopy data can be used to determine the stiffness constant, elastic constants and Debye temperature etc. [66, 67]. The equations used to calculate C is given by [68] 𝐶 = 𝐾 𝑎𝑣 𝑎 , where K a v is the average of tetrahedral and octahedral force constant calculated Wave number (cm -1 ) T r a n s m i tt a n ce ( % ) x =0.00 x =0.01 x =0.02 x =0.03 x =0.04 x =0.05 Fig. 5.
FTIR spectra of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05) . 0 [Table 2] from FTIR spectra using the following equation [69]: 𝐹 𝑐 = 4𝜋 𝑐 𝜈 𝑚, where, c, υ and m are the speed of light in free space, vibrational frequency and the reduced mass of Fe and O ions, respectively. The equivalent reduced mass is about 2.061×10 -23 g and a is the lattice constant. The relation [68]: 𝐶 = 𝜍𝐶 is used to estimate the C , where σ is the Poisson’s ratio [70]: 𝜍 = 0.324 × [1 − 1.043 × 𝑃 ] , P is the pore fraction. The C is calculated from the relation: 𝐶 = 𝐶 −𝐶 . The obtained values of C ij [Table 3] confirmed that the ferrites compositions under consideration satisfy the Born criteria for mechanical stability [71, 72]: 𝐶 − 𝐶 > 0; 𝐶 + 2𝐶 > 0; 𝐶 > 0. Therefore, the synthesized compositions are mechanically stable. The stiffness constant C represents the elasticity of length and provides the stiffness of materials along the (100) direction while the constants C and C born the elasticity in shape. The variation in stiffness constants with Y contents is shown in Fig. 6 (a). Fig. 6.
The variation of (a) stiffness constants ( C ij ) and elastic moduli with Y contents. The values of C and C increase up to x = 0.02 and then C gradually decreases while C decreases up to x = 0.04 and increases further. The improvement in C and C endorse the good solubility of Y ions into B -sites and strengthen of the inter-atomic bonding. While the decline trend in C and C curves indicates the decrease in crystallization process. Microstructure images [Fig. 3] also ratify the crystallization of prepared samples. In this case, the variation in C and C and bulk density [Table 3] may be attributed from the variation of grains size that can be observed from Fig. 3. [73]. The equations [70]: 𝐸 = 𝐶 −𝐶 (𝐶 +2𝐶 )(𝐶 +𝐶 ) , 𝐵 = +2𝐶 )3 , C a nd C ( G P a ) Y contents ( x ) C C C ( G P a ) C (a) EG Y contents ( x ) B a nd G ( G P a ) B (b) E ( G P a ) 𝐺 =
𝐸2(𝜍+1) are used to calculate the Young’s modulus ( E ), Bulk modulus ( B) and shear modulus ( G ) and presented in Table 3. Fig. 6 (b) represents the variation in different elastic moduli [ B , G and E ] with Y contents ( x ). The bulk modulus provides the value of required resistance for volume deformation, an increase is observed up to x = 0.02 due to improvement in microstructure and then decreases gradually. This variation might be related to the variation in bulk density as well as the grains size ( x ) [73]. The shear and Young’s modulus are noted to decrease for x = 0.01 and then increase gradually with the variation of x . The value of E provides the stiffness of solids where the high value of E indicates the stiffer solids. The values of E and G of substituted compositions are lower than that of the parent one [ x = 0.00]. The change in inter-atomic bonding is responsible for the variation in elastic moduli [26]. Although, the variation of elastic moduli [Fig. 6] is shown as a function of Y contents in magnified scale, but if we have a look on the values presented in Table 3, an insignificant variation is observed. Table 3 shows the Poisson’s ratio [ σ ] for all compositions. The Poisson’s ratio is extremely important in materials application point of view which determines the ductile behavior or brittleness of materials with a decisive factor of 0.26 [74, 75]. A value of σ less than 0.26 indicates the brittleness of materials otherwise ductile behavior of materials. Like other ceramics materials, brittleness of studied compositions is confirmed and the values of σ for all compositions lie in the expected range from -1 to 0.5, confirming excellent elastic behavior with the theory of elasticity [76]. The anisotropy constant A ( A = 2 C /( C - C )) [77-78] also calculated [Table 3] and found to have the value of 1 (one) for all composition, indicating isotropic nature of cubic ferrites compositions. The Debye temperature ( Θ D ), the most fundamental property of solids, used to resolve many of characteristics thermal properties; determine the lower temperature region and higher temperature region for solids. Moreover, the Θ D also provides the information regarding bonding strength. The higher value of Θ D is associated with stronger bonds in solids. The average sound velocity can be used to calculate Θ D by the following equation [79]: Θ D = h k B
3n 4π N A ρ/M 𝑣 m , where the symbols signify the usual meanings [75]. The average sound velocity is calculated by the equation [79]: 𝑣 m = 1 3 1 𝑣 l3 + 2 𝑣 s3 −1 3 , where v l [ 𝑣 l = 3B + 4G 3ρ ] and v s [ 𝑣 s = G ρ ] are the longitudinal and transverse sound velocities in an isotropic material. Fig. 7.
The variation of sound velocities and Debye temperature with Y contents in the composition of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05). Table 3
The lattice constant ( a ), bulk density (ρ b ), porosity ( P ), average force constant ( K av ), stiffness constant ( C ij ) elastic moduli ( E, B, G ), Poisson’s ratio ( σ ), anisotropic constant ( A ), sound velocities ( v l , v s , v m ) and Debye temperature ( Θ D ). Parameters x = x = x = x = x = x = a (Å) ρ b (kg/m ) 4006 4172 4197 4166 4119 4084 P (%) 18.08 14.76 14.25 15.06 16.38 17.17 K av. (N/m) 173.40 173.40 174.58 173.40 172.24 173.99 C (GPa) 205.55 205.59 206.86 205.38 204.01 206.14 C (GPa) 32.19 42.10 43.99 41.15 36.89 34.93 C (GPa) 86.68 81.70 81.43 82.11 83.56 85.60 E (GPa) 196.83 191.18 191.43 191.64 192.71 196.01 B (GPa) 89.98 96.57 98.28 95.89 92.60 92.00 G (GPa) 86.68 81.70 81.43 82.11 83.56 85.60 σ A v (m/s) 6918.9 6742.7 6740.1 6748.3 6777.8 6851.1 v s (m/s) 4365.4 4091.0 4063.0 4109.9 4194.3 4278.7 v m (m/s) 4803.9 4520.7 4492.6 4539.8 4625.3 4713.8 Θ D (K) 607.43 579.12 576.38 580.70 589.12 598.39 The calculated sound velocities mentioned are presented in Table 3 along with the calculated Θ D . It is observed [Table 3] that the transverse sound velocity is lower than the longitudinal sound velocity. Energy is transferred from particle to particle by vibration of particle during v s v l v l , v s , v m ( m / s ) Y contents ( x ) v m D D ( K ) propagation of wave through a medium. During the propagation of transverse wave, the vibration of particle is at right angles with the direction of wave transmission, caused a higher energy requirement to vibrate the adjacent particle. As a result the energy of the wave reduced and hence the velocity of transverse wave is lower than the longitudinal wave and is roughly half of longitudinal wave velocity [80]. Fig. 7 illustrates the variation of Θ D with Y contents in which it seen that the Θ D decreases with Y contents up to x = 0.02 and then increases with Y contents. The similarity in the variation of sound velocities and Θ D is observed similar to other ferrites [26, 73, 76]. The variation of Θ D roughly follows the variation of rigidity modulus [76] as can be shown in Fig. 6 (b). 3.5 dc resistivity High resistivity of ferrites is essential for most electronic application and therefore, the study of dc electrical resistivity is important for prediction of their application in electronic and high frequency devices. The room temperature dc electrical resistivity ( 𝜌 𝑅𝑇 ) of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05) is estimated by the following equation: 𝜌 𝑅𝑇 = 𝑅𝐴𝐿 , where,
R, A and L is resistance, area and thickness, respectively of the pellet. The value of 𝜌 𝑅𝑇 is noted to increased due to Y substitution unlike for x = 0.01 (Fig. 8 (a)), can be explained by Verwey mechanism based on electrons exchange involving the ions of the element with different valence states positioned octahedral ( B -site) [81-82]. In the present case, the Fe and Fe is responsible for the conduction because of the hopping tendency between them [83, 84] and the Fe ions are the product of sintering process [42]. The value of 𝜌 𝑅𝑇 is found to be 0.45 × 10 Ω-cm for x = 0 and 0.40 × 10 Ω-cm, 1.57 × 10 Ω-cm, 1.78 × 10 Ω-cm, 2.19 × 10 Ω-cm and 1.43 × 10 Ω-cm for x = 0.01, 0.02, 0.03, 0.04 and 0.05, respectively. Fig. 8 (a) clearly indicates the increase in resistivity for substituted composition (except for x = 0.01) where the resistivity decreases slightly. The ac resistivity was found to decrease for x = 0.01 [57]. The Y ions replace Fe ions at B -sites. [50], simultaneously small fraction of Zn ions forfeit balanced by forcing some Fe to move at the A -sites. The Fe ions have higher electronic valence than Zn ions; hence the charge carriers (here, metallic vacancy) have been increased by balancing the electrical charge, results decrease in resistivity. Two mechanisms might be involved in increasing resistivity, firstly, the enhancement in resistivity due to Y substitution is attributed by forming the stable bonds (electric) of Y and Fe , causes localization of Fe , limiting the Verwey mechanism [81], hence upturn the resistivity. Secondly, the preference of Y ions to reside in the B -sites, lead to a replacement of Fe ions at B -sites, results a reduction of Fe formation. Since, Y ions do not take place in the conduction due to its stable oxidation state and electrons hopping takes place between Fe and Fe , therefore, the degree of electrons change (Fe ↔ Fe ) is reduced and consequently the resistivity is increased. The resistivity of materials greatly depends on its microstructure: the porosity and the distribution of grains. The lowering of resistivity for x = 0.05 contents may be explained using the microstructure [Fig. 3(f)] of substituted compositions and the improved connectivity between the grains [85]. From Fig. 3 (d-e), it is seen that the grains distribution are more homogeneous for x = 0.05 than that for x = 0.03 and 0.04; and expected to have the improved connectivity between the grains. Therefore, a slight decrease in resistivity is also expected. Fig. 8. (a) variation dc resistivity with Y contents, and (b) plot of [
F(R ∞ )hν ] against hν of Mg Zn Y x Fe x O (0≤ x ≤0.05) ferrites. 3.6 Diffuse reflectance results
To know the optical properties of Y-substituted Mg-Zn ferrites, the UV–Vis diffuse reflectance study has been performed. The optical reflectance data can be used to calculate the band gap using the equation
𝐹 𝑅 ∞ = 𝐴 ℎ𝜈−𝐸 𝑛 ℎ𝜈 , where, F(R ∞ )=(1-R) /2R ∞ , is the Kubelka-Munk Function [86]. Therefore, the plot of [F(R ∞ )hν] against photon energy (hν) gives the way of calculating the energy band gap E g (eV), by extrapolating the linear part to the x-axis (energy axis). The Y contents ( x ) ( - c m ) (a) h (eV) x =0.00 E g = 2.56 eVE g = 2.59 eV x =0.01 E g = 2.71 eV x =0.02 E g = 2.73 eV x =0.03 [ F ( R ) h ] E g = 2.78 eV x =0.04 E g = 2.58 eV x =0.05 (b) band gap is found to be increased due to the Y substitution that is consistent with the Y-substituted Co ferrites [53]. The variation in band gap indicates the change in structural parameters attributed from the change in Y contents [53]. The band gap enhancement is associated with the variation of dc resistivity [Fig. 8 (b)]. Magnetic properties
High temperature magnetic properties mapping
The magnetization vs applied magnetic field (upper part of hysteresis loops) at different measuring temperatures, as shown in Fig. 9 (a-f), exhibiting typical magnetic nature at elevated temperature, revealing the existence of magnetic ordered structure near T c (in section 3.7.3). Beyond T c , destruction of magnetic order is also observed i.e., transition from ferrimagnetic (below T c ) to paramagnetic (above T c ) nature. It is [Fig. 9 (a-f) also expected that the coercivity (H c ) and remanent magnetization (M r ) decrease with increasing temperature and tend to be zero near T c . It is very important that the ferrimagnetic ordered structured remains up to T c indicating the suitability of these materials to use them at higher temperature. M ( e m u / g ) Applied field (Oe) C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C (b) (a) Applied field, H (Oe) M ( e m u / g ) C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C Fig. 9 (a-f)
The M vs H plot of Mg Zn Y x Fe x O (0 ≤ x ≤ 0.05) measured at different temperatures. The value of M s for each composition at different selected temperatures is obtained by the following: at first the M s vs 1/H curves are plotted (figures not shown here) and are extrapolated to zero i.e., 1/H = 0 [53]. The obtained M s values are plotted in Fig. 10 (a) for different Y contents at different temperature. The variation of M s with Y contents has been explained elsewhere [57]. A decrease in the M s with increasing temperature is observed due to more randomness of magnetic moments that attributed from the increase in thermal energy. A change in site occupancy of cation may also induce due to thermal fluctuation; results reduced M s at higher temperatures. C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C Applied field, H (Oe) M ( e m u / g ) (c) C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C Applied field, H (Oe) M ( e m u / g ) (d) M ( e m u / g ) Applied Field, H (Oe) C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C (e) (f) C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C Applied Field, H (Oe) M ( e m u / g ) Fig. 10.
Variation of (a) saturation magnetization at different measuring temperatures, (b) temperature dependent magnetic moment and (c) dM/dT vs temperature graph for different Y contents. It is interesting to note that the width of the magnetic transition, inferred from the magnitude of the FWHM of the dM/dT vs T curve [Fig. 10 (c)] and the magnetic ordering temperature (taken at the peak of the dM/dT versus T graph) show features strongly dependent on the level of Y substitutions. For the Y substituted compounds [Fig. 10 (c)], the transition width does not show any systematic variation with Y content. The Y-free compound is characterized by a very large transition width. The transition temperature lies within a narrow range of 341 ± 4 K. This implies that the average mean-field exchange energy among the spins remains insensitive to the Y content. It should be noted that, substitutional disorder always introduces local strains inside the x = 0.00 Temperature ( C) x = 0.01 x = 0.02 x = 0.03 M ( e m u ) x = 0.04 x = 0.05 (b) (a) C C C ) C C M s ( e m u / g ) Y contents nge interactions inside these samples.
Fig. 5.2.29 dM/dT vs temperature graph for Y substituted compositions.
360 420 480 540 600-0.0030.000
Y=0.02 d M / d T T (K)
360 420 480 540 600-0.0030.000
Y = 0.03 d M / d T T (K)
360 420 480 540 600-0.0030.000
Y = 0.04 d M / d T T (K)
360 420 480 540 600-0.006-0.0030.000
Y = 0.05 d M / d T T (K)
360 420 480 540 600-0.0030.000
Y = 0.00 d M / d T T (K) 360 420 480 540 600-0.0030.000
Y = 0.01 d M / d T T (K) compounds. These local distortions and inhomogeneities lead to a distribution in the exchange energies. Within the mean-field approximation, the Curie temperature is directly proportional to the exchange energy [87]. Structure specific corrections often lower the theoretically predicted transition temperature, but this link between transition temperature and exchange energy remains valid. Within this framework, width of the magnetic transition is directly related to the distribution and homogeneity of atoms inside the solid. Therefore, different transition widths in the Y substituted compounds reflect the distribution and inhomogeneity of the exchange interactions inside these samples. 3.7.2 Curie temperature
The magnetic moment and initial permeability and hence the magnetic properties are very sensitive to temperature. Both exhibit an abrupt dropping near the transition point from ferrimagnetic to paramagnetic state. Therefore, the Curie temperature (T c ) has been measured by applying the two basic characteristics: (i) the samples are subjected to an applied magnetic field of 1000 Oe and the magnetic moment is calculated at different temperature. The heating rate is kept slow to avoid the nucleation and particle growth at elevated temperatures. The variation of magnetic moment with temperature is shown in Fig.10 (b) where the magnetic moment is found to decrease with increasing temperature as expected and tends to zero at different measuring temperature for different compositions. The T c value is then calculated from the first derivative of magnetic moment against temperature graph, Fig. 10 (c). Moreover, the T c is also calculated from the temperature dependent initial permeability, µ i plot as shown in Fig.11 (a). A peak, also known as Hopkinson peak, is observed from the graph where the value of initial permeability dropped sharply. The characteristic temperature at which the peak observed is known as the T c where the spin magnetic moments are completely disordered and the conversion from ferrimagnetic to paramagnetic materials is taken place. The obtained T c is found to be decreased with the Y substation. The value of T c depends on the cation distribution on A- and B- site and the exchange coupling constant [53]. The coupling constant due to exchange interaction between ions in sub-lattice ( A and B ) is much stronger than the coupling constant due to exchange interaction among the ions in the same sub-lattice sites ( A/B ). The substitution of Y ions in the spinel structure leads to increase in the bond lengths and distance between ions presents at the A - and B –sites that results a decrease in exchange constants and hence the T c values [53]. The T c value for the Mg Zn Fe O ( x = 0.00) is obtained from the M-T graph and µʹ i -T graph and found to be 77°Cand 75°C, respectively which is higher than that of (T c ~47°C) reported by Khot et al. [88] and Hashim et al . [44] and lower than that of (T c ~97°C) obtained by Mazen et al . [89]. Fig. 11. (a) The temperature dependent initial permeability, and (b) the variation of Curie temperature with Y contents. 3.7.3 Phase transition
In section 3.7.1, the ferrimagnetic to paramagnetic phase transition has been observed and the transitions temperatures have also been obtained (section 3.7.2). The determination of the nature of phase transition (first or second order) is also significant. The Arrott plot (M (emu/g) against H/M (Oe.g/emu)) has been used to see the nature of phase transition. Banerjee [90] has proposed the criteria for first and second order transition. The positive slope of Arrott plot indicates the second order phase transition otherwise the transition is first order. The curves intercept at M (Y-axis) indicate the presence of spontaneous magnetization while the curve passing through the origin measure the T c [91]. Fig. 12 (a-f) illustrates that the second order phase transition has been occurred in the prepared samples. It is also observed that the spontaneous magnetization is present up to ≤ 75 °C for x = 0.00, ≤ 69 °C for x = 0.01, ≤ 67 °C for x = 0.02, ≤ 65 °C for x = 0.03 and ≤ 63 °C for x = 0.04 and 0.05 and thereafter the curves (after T c for each compositions) are straight upwards from H/M axis ( x -axis) which indicates the paramagnetic behavior.
40 60 80 100200300400500600700 x =0.00 x =0.01 x =0.02 x =0.03 x =0.04 x =0.05 i Temperature ( C) (a) T c ( C ) Y contents
From -T From M-T(b) Fig. 12 . The Arrott plot of Mg Zn x Fe x O (0≤ x ≤0.05) ferrites. C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C H/M (Oe.g/emu) M ( e m u / g ) (a) H/M (Oe.g/emu) (b) M ( e m u / g ) C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C H/M (Oe.g/emu) M ( e m u / g ) (c) C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C H/M (Oe.g/emu) M ( e m u / g ) (d) M ( e m u / g ) H/M (Oe.g/emu) C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C (e) H/M (Oe.g/emu) M ( e m u / g ) (f) C 47 C 52 C 57 C 62 C 67 C 69 C 71 C 73 C 75 C 77 C 79 C 81 C 83 C 85 C 87 C 92 C 97 C 102 C 127 C 147 C 167 C 187 C Conclusions
The Yttrium-substituted Mg-Zn [Mg Zn x Fe x O (0 ≤ x ≤ 0.05)] ferrites, with spinel structure as confirmed from the XRD and FTIR studies, have been synthesized using ceramic technique. The calculated stiffness constant revealed the mechanical stability of the compositions. The bulk modulus is increased for substituted compositions up to x = 0.2 and decreased thereafter. The shear modulus and Young’s modulus are found to be decreased at x = 0.01 and then increased with further increase in Y contents. The values of Poisson’s ratio (0.13 to 0.17) indicating the brittleness character of studied compositions. The value of anisotropic constants ( A = 1) be a sign of isotropic nature of studied spinel ferrites. The variation of Θ D with increasing Y contents is associated with the variation in rigidity of the compositions. The dc resistivity increases from 0.45×10 Ω-cm (for x = 0.00) to 1.57×10 Ω-cm, 1.78×10 Ω-cm, 2.19×10 Ω-cm and 1.43×10 Ω-cm for x =0.02, 0.03, 0.04 and 0.05, respectively. The obtained optical energy band gap are found to be 2.56 eV (for x = 0.00) to 2.59 eV, 2.71 eV, 2.73 eV, 2.78 eV and 2.58 eV for x = 0.01, 0.02, 0.03, 0.04 and 0.05, respectively. Both the ferrimagnetic and paramagnetic nature have been confirmed from temperatures dependent magnetization curves. Very close values of T c obtained from the M-T and µ ʹ i -T is noticed which is decreased due to the weakening of exchange coupling constant with Y substitution. The second order phase transition (ferrimagnetic to paramagnetic) has been occurred. The prepared ferrites are suitable for high frequency application as well as at temperature higher than room temperature (up to T c ). Acknowledgments
The authors are grateful to the Directorate of Research and Extension, Chittagong University of Engineering and Technology, Chattogram-4349, Bangladesh under grant number CUET DRE (CUET/DRE/2016-2017/PHY/005) for arranging the financial support for this work. We are thankful for the laboratory support of the Materials Science Division, Atomic Energy Center, Dhaka-1000, Bangladesh.
References A. Goldman, Modern Ferrites Technology, 2nd ed. (Springer, New York, USA, 2006). M. A. Ahmed, H. E. Hassan, M. M. Eltabey, K. Latka, T. R. Tatarchuk, Phys. B: Phys. Cond. Mater. 530 (2018) 195. B. R. Babu, T. Tatarchuk, Mater. Chem. Phys. 207 (2018) 534. T. Tatarchuk, M. Bououdina, W. Macyk, O. Shyichuk, N. Paliychuk, I. Yaremiy, B. AlNajar, M. Pacia, Nanoscale res. lett. 12(2017) 1. T. Tatarchuk, M. Bououdina, N. D. Paliychuk, I. P. Yaremiy, V. V. Moklyak, J. Alloys Compd. 694 (2017) 777. D. K. Pradhan, S. Kumari, V. S. Puli, P. T. Das, D. K. Pradhan, A. Kumar, J. F. Scott, R. S. Katiyar, Phys. Chem. Chem. Phys. 19 (2017) 210. A. Yadav, D. Varshney, Superlattice Microst. 113 (2017) 153. S.K. Rashmi, H.S. B. Naik,, H. Jayadevappa, R. Viswanath, S.B. Patil, M. M. Naik, Mater. Sci. Eng. B 225 (2017) 86. S.B. Patil, H.S.B. Naik, G. Nagaraju, R. Viswanath, S.K. Rashmi, Eur. Phys. J. Plus, 132 (2017) 328.
M. Li, Q. Gao, T. Wang, Y-S. Gong, B. Han, K-S. Xia, C-G. Zhou, Mater. Design, 97 (2016) 341.
M. Sundararajan, V. Sailaja, L. J. Kennedy, J. J. Vijaya, Ceram. Int. 43 (2017) 540.
S.F. Mansour, M.A. Abdoa, S.I. El-Dek, J. Magn. Magn. Mater. 422 (2017) 105.
B. Skołyszewska, W. Tokarz, K. Przybylski, Z. Ka̧ kol, Physica C 387 (2003) 290.
R. Sharma, P. Thakur, P. Sharma, V. Sharma, J. Alloys Compd. 704 (2017) 7.
M. U. Rana, M. U. slam, T. Abbas, Pakistan J. Appl. Sci. 2 (2002) 1110.
M. A. Ali, M.N.I. Khan, F.-U.-Z. Chowdhury, M.M. Hossain, M.Z. Rahaman, S.M. Hoque, M.A. Matin, M.M. Uddin, Results in Physics 14 (2019) 102517.
H. Kavas, A. Baykal, M. S. Toprak, Y.Koseoglu, M.Sertkol, B.Aktas, J. Alloys Comp. 479 (2009) 49.
P. Thakur, R. Sharma, M. Kumar, S.C. Katyal, N.S. Negi, N. Thakur, Vineet Sharma, P. Sharma, Mater. Res. Exp. 3 (2016) 075001.
K. El-Sayed, M. B.Mohameda, Sh. Hamdy, S.S. Ata-Allah, J. Magn. Magn. Mater. 423 (2017) 291.
K. A. M. Khalaf, A. D. Al-Rawas, H. M. Widatallah, K. S. Al-Rashdi, A. Sellai, A. M. Gismelseed, M. Hashim, S. K. Jameel, M. S. Al-Ruqeishi, K. O. Al-Riyami, M. Shongwe, A. H. Al-Rajhi, J. Alloys Compd. 657 (2016) 733.
M. Hashim, Alimuddin, S. Kumar, S. E. Shirsath, R.K. Kotnala, J. Shah,R. Kumar, Mater. Chem. Phys. 139 (2013) 364.
M. D. Hossain, M. N. I. Khan, A. Nahar, M. A. Ali, M. A. Matin, S. M. Hoque, M. A. Hakim, A. T. M. K. Jamil, J. Magn. Magn. Mater. 497 (2020) 165978.
C. C. Naik, A V Salker, Mater. Res. Express 6 (2019) 066112.
M. Atif, Mater. Res. Express 6 (2019) 076104.
K. Bindu, K. M. Ajith, H. S. Nagaraja, Mater. Res. Express 6 (2019) 045011.
A. Bhaskar, S.R. Murthy, J. Magn. Magn. Mater. 355 (2014) 100.
V.K. Lakhani, K.B. Modi, Solid State Sci. 12 (2010) 2134.
Mukhlis M. Ismail, Nasma A. Jaber, J. Braz. Soc. Mech. Sci. Engg. 40 (2018) 250.
S. E. Shirsath, S.M. Patange, R.H. Kadam, M.L. Mane, K.M. Jadhav, J. Mol. Struct. 1024 (2012) 77.
A. Franco, F. C. e Silva, Appl. Phys. Lett. 96 (2010) 172505.
L. Zegadi, J. J. Rousseau, B. Allard, P. Tenant, D. Renault, IEEE Trans. Magn. 36 (2000) 2022.
P. Andrei, A. Stancu, H. Hauser, P. Fulmek, J. Optoelectron. Adv. Mater. 9 (2007) 1137.
A. Franco, M. S. Silva, J. Appl. Phys. 109 (2011) 07B505.
T.E.P. Alves, H.V.S. Pessoni, A. Franco Jr. Phys. Chem. Chem. Phys., 19 (2017) 16395.
E. Melagiriyappa, H.S. Jayanna, B.K. Chougule, Mater. Chem. Phys. 112 (2008) 68.
M. M. Haque, M. Huq, M.A. Hakim, Mater. Chem. Phys. 112 (2008) 580.
H. Mohseni., H. Shokrollahi, I. Sharifi, Kh. Gheisari, J. Magn. Magn. Mater. 24 (2012) 3741.
U.R. Ghodake, Rahul. C. Kambale, S.S. Suryavanshi, Ceram. Int. 43 (2017) 1129.
A. Yadav, D. Varshney, J Supercond Nov Magn 30 (2017) 1297.
A. Ashrafizadeh, Ali Ghasemia, A. Paesano Jr., C. F. C. Machado, X. Liu, A. Morisako, J. Alloys Compd. 506 (2010) 279.
A. Yadav, R. Rajpoot, M.A. Dar, D. Varshney, AIP Conf. Proc. 1731 (2016) 050008.
S.J. Haralkar, R.H. Kadam, S.S. More, S.E. Shirsath, M.L. Mane S. Patil, D.R. Mane, Physica B 407 (2012) 4338.
A. Xia, S. Su, S. Liu, C. Jin, J. Mater. Sci.: Mater Electron 24 (2013) 4166. M. Hashim, S. S. Meena, R. K. Kotnala, S. E. Shirsath, P. Bhatt, S Kumar, E. Şentürk, R. Kumar, A. Alimuddin, N. Gupta, J. Magn. Magn. Mater. 360 (2014) 21.
R. Sharma, P. Thakur, M. Kumar, N. Thakur, N.S. Negi, P. Sharma, V. Sharma, J. Alloys Compd. 684 (2016) 569.
F. Hirosawa, T. Iwasaki, S. Watano, Appl. Nanosci. 7 (2017) 209.
R. A. Bugad, T. R. Mane, B. B. Navale, J. V. Thombare, A. R. Babar, B. R. Karche, J Mater Sci: Mater Electron. 28 (2017) 1590.
M. W. Mukhtar, M. Irfan, I. Ahmad, I. Ali, M. N. Akhtar, M. A. Khan, G. Abbas, M.U. Rana, A. Ali, M. Ahmad, J. Magn. Magn. Mater.381 (2015) 173.
R. Sharma, P. Thakur, M. Kumar, P. Sharma, V. Sharma, J. Alloys Compd. 746 (2018) 532.
V. Jančárik, E. Ušák, M. Šoka, M. Ušáková. Acta Physica Polonica Series A 126 (2014) 90.
S. Kumari, V. Kumar, P. Kumar, M. Kar, L. Kumar, Adv. Powd. Tech. 26 (2015) 213.
S. E. Jacobo, P. G. Bercoff, Ceram. Int. 42 (2016) 7664.
T.E.P. Alves, H.V.S. Pessoni, A. Franco Jr. Phys. Chem. Chem. Phys., 19 (2017) 16395.
S. Chakrabarty, A. Dutta, M. Pal, J. Magn. Magn. Mater. 461 (2018) 69.
M. Ishaque, M. A. Khan, I. Ali, M. Athair, H. M. Khan, M. A. Iqbal, M.U. Islam, M. F. Warsi, Mater. Sci. Semi.Proces. 41 (2016) 508.
M. Ishaque, M. A. Khan, I. Ali, H. M. Khan, M. A. Iqbal, M.U. Islam, M. F. Warsi, Ceram. Int. 41 (2015) 4028.
M.A. Ali, M.N.I. Khan, F.-U.-Z. Chowdhury, M.M. Hossain, A.K.M. Akhter Hossain, R. Rashid, A. Nahar, S.M. Hoque, M.A. Matin, M.M. Uddin, J. Mater. Sci.: Mater. Electr. 30 (2019) 13258.
M. A. Ali, M. M. Uddin, M.N.I. Khan, F.-U.-Z. Chowdhury, S.M. Hoque, S.I. Liba, Chin. Phys. B 26 (2017) 077501.
M. A. Ali, M. M. Uddin, M. N. I. Khan, F.-U.-Z. Chowdhury, S. M. Hoque, J. Magn. Magn. Mater. 424 (2017) 148.
M. A. Ali, M. N. I. Khan, F.-U.-Z. Chowdhury, S. Akhter, and M. M. Uddin, J. Sci. Res. 7 (2015) 65.
M A Ali, M. M. Uddin, M. N. I. Khan, F -U -Z Chowdhury, D. K. Saha, S. M. Hoque, S. I. Liba, S. Akhter; arxiv:1506.06280.
A.B. Gadkari, T.J. Shinde, P.N. Vasambekar, Mater. Chem. Phys. 114 (2009) 505.
R. D. Waldron, Phys. Rev.99 (1955) 1727.
O.M. Hemeda, M.Z. Said, M.M. Barakat, J. Magn. Magn. Mater. 224 (2001) 132.
W. B. White, W. A. De Angelis, Spectrochim. Acta, Part A, 23 (1967) 985.
K.B. Modi, J. Matter. Sci. 39 (2004) 2887.
K.B. Modi, M.C. Chhantbar, H.H. Joshi, Ceram. Int. 32 (2005) 111.
K.B. Modi, J.D. Gajera, M.P. Pandya, H.G. Vora, H.H. Joshi, Pramana 62 (2004)1173.
S.J. Keny, J. Manjanna, G. Venkateswaran, R. Kameswaran, Corr. Sci. 48 (2006) 2780.
S.S. Bhatu, V.K. Lakhani, A.R. Tanna, N.H. Varsoya, J.U. Buch, P.U. Sharma, U.N. Trivedi, H.h. Joshi, K.B. Modi, Ind. J. Pure Appl. Phys. 45 (2007) 596.
G. Grimvall, Thermophysical Properties of Materials (North-Holland, Amsterdam,1986).
M. A. Ali, M. Roknuzzaman, M. T. Nasir, S. H. Naqib, A. K. M. A. Islam, Int. J. Mod. Phys. B 30 (2016) 1650089.
K. Sadhana, K. Praveena, S.R. Murthy, J. Magn. Magn. Mater. 323 (2011) 2977.
I. N. Frantsevich, F. F. Voronov, and S. A. Bokuta, Elastic Constants and Elastic Moduli of Metals and Insulators Handbook (Naukova Dumka, Kiev, 1983) 60.
M. A. Ali, M. A. Hadi, M. M. Hossain, S. H. Naqib, A. K. M. A. Islam, Phys. Status Solidi B 254 (2017) 1700010.
S.M. Patange, Sagar E. Shirsath, K.S. Lohar, S.G. Algude, S.R. Kamble, Nilesh Kulkarni, D.R. Mane, K.M. Jadhav, J. Magn. Magn. Mater. 325 (2013) 107.
C. M. Zener S. Siegel, J. Phys. Chem. 53 (1949) 53.
M. A. Ali, A. K. M. A. Islam, Phys. B 407(2012) 1020. O. L. Anderson, J. Phys. Chem. Solids 24 (1963) 909.
V.K. Lakhani, K.B. Modi, Solid State Sci. 12 (2010) 2134.
E. Verwey, J. De Boer, Recl. Trav. Chim. Pays-Bas, 55 (1936) 531.
H. S. Aziz, S. Rasheed, R. A. Khan, A. Rahim, J. Niser, S. M. Shah, F. Iqbal, A. R. Khan, RSC Adv. 6 (2016) 6589.
M. A. Matin, M. N. Hossain, F. A. Mozahid, M. A. Ali, M. A. Hakim, M. F. Islam, SN Appl. Sci. 1 (2019) 1479.
N. Rezlescu, D. Condurache, C. Naum, E. Luca, Rev. Roum. Phys. 18 (1973) 727.
L. Chauhan, A. K. Shukla, K. Sreenivas, Ceram. Int. 41 (2015) 8341.
P. Kubelka, F. Munk, An article on optics of paint layers. Z. Tech. Phys. 12 (1931) 593.
J. S. Kouvel, Phys. Rev. 102 (1956) 11489 and references therein.
S. S. Khot, N. S. Shinde, B. P. Ladgaonkar, B. B. Kale, S. C. WataweAdv. Appl. Sci. Res. 2 (2011) 460.
S. A. Mazen, S. F. Mansour, H. M. Zaki, Cryst. Res. Technol. 38 (2003) 471.
B. K. Banerjee, Phys. Lett., 12 (1964) 16.