3.1. X-ray diffraction, surface morphology and Chemical Composition
Figure 1a depicts the X-ray diffraction patterns of the pure CoFe
2O
4 that was dried using infrared (IR) lamp, dried with oven at 200 °C (natural dried), and the Co
1-2xNi
xMn
xFe
2-yCe
yO
4 ferrites spinels resulting from the substitution of tetrahedral A-sites and octahedral B-sites of CoFe
2O
4 with Ni, Mn, and Ce atoms in a linearly increasing molar amount.
The substitution was done at equal molar amount of x and y, where 0 ≤ x = y ≤ 0.3. These substituted spinel nanoferrites were dried naturally at 200 °C for 4 hours and subjected to annealing as described in the experimental section. The final products of the appropriate substitution of Ni, Mn and Ce atoms are Co0.8Ni0.1Mn0.1Fe1.9Ce0.1O4, Co0.6Ni0.2Mn0.2Fe1.8Ce0.2O4, and Co0.8Ni0.1Mn0.1Fe1.9Ce0.1O4, including the unsubstituted CoFe2O4 spinel ferrites. The samples are labelled S1 (IR lamb dried dried CoFe2O4), S2 (natural-dried CoFe2O4), S3 (Co0.8Ni0.1Mn0.1Fe1.9Ce0.1O4), S4 (Co0.6Ni0.2Mn0.2Fe1.8Ce0.2O4), and S5 (Co0.4Ni0.3Mn0.3Fe1.7Ce0.3O4).
The diffraction patterns of the samples having peaks at 2θ positions can be indexed as follows: 30.18° (220), 35.58° (311), 43.39° (400), 53.82° (422), 57.23° (511) and 62.91° (440). The position of the XRD peaks for the samples matched well with the standard XRD pattern for CoFe
2O
4 spinel ferrites with JCP2 card No. 96-154-0974, and it is in good agreement with the XRD results previously reported elsewhere [
19,
24,
25,
26,
27]. It is noteworthy to mention that at x = y ≥ 0.2, additional peaks at 2θ positions, 28.63, 33.10 and 47.87° were observed. These peaks correspond respectively to (111), (200) and (220) planes of CeO
2 with PDF card no. 96-721-7888 [
19,
28,
29,
30]. This indicates the formation of CeO
2 secondary phase in the spinel structure. No other phase attributable to Ni and Mn oxides was observed; this is owing to the fact that ionic radii of Ni
2+, Mn
2+ and Co
2+ are very close to each other with 0.083, 0.081, and 0.079 nm, respectively. However, the ionic radius of Ce
3+ (0.101 nm) is greater than that of Fe
3+ (0.064 nm). As a results, cerium dioxides secondary phase was formed as Ce content increased in the sample. Mkwae et al. [
31] reported the formation of the CeO
2 as a secondary phase at higher Ce dopant content during the synthesis of MgCe
xFe
2-xO
4 by hydrothermal synthesis method. Many other researchers have previously observed the CeO
2 secondary phase when Ce
3+ was substituted in their spinel nanoferrites [
19,
32,
33,
34]. In addition, the crystallinity of the sample with x = y = 0.1 is less than that of the other samples as some peaks corresponding to the spinel appears to be absent or inconspicuous as observed in their low intensity.
Table 1 shows the crystallite sizes of the samples, calculated using Debye-Scherrer’s equation [
35].
Figure 1b–f shows the scanning electron microscope (SEM) images of the spinel nanoferrites samples. The two CoFe
2O
4 samples- infrared lamb dried and natural-dried in
Figure 1b,c, possess similar surface morphology with predominantly fine particles and few agglomerates. Upon addition of impurities to the spinel structure, difference in surface morphology becomes distinct especially on the x = y = 0.1 shown in
Figure 1d. This sample contains more agglomerates than fine nanoparticles as compared to the previous two samples with no impurities- Ni, Mn and Ce cations. The formation of these closely packed agglomerates could be due to the introduction of the Ni, Mn, and Ce dopants into the spinel, which might have led to the complication of the spinel structure. At the substitution of 0.1 molar concentration of the impurities, i.e., when x = y = 0.1, the crystallinity of the spinel decreased as observed earlier in the x-ray diffraction (XRD) analysis. At x = y = 0.2, the spinel looks like solidified magma sparsely covered by agglomerates (see
Figure 1e). This is more crystalline than the sample with x = y = 0.1 as can be observed in the high XRD peak intensities in
Figure 1a with secondary phases corresponding to the CeO
2. The sample depicted in
Figure 1f has the appearance of tightly bound particles.
Pore sizes and pore size distribution were estimated by conducting nitrogen adsorption-desorption measurement on the sample. The samples are generally mesoporous, having pore size p, where 2 nm < p < 50 nm.
Figure 2 depicts the nitrogen adsorption-desorption isotherm of the Co
0.8Ni
0.1Mn
0.1Fe
1.9Ce
0.1O
4 double-substitution spinel. It is a type IV isotherm [
36] which is characteristic of mesoporous materials. The inset is the pore size distribution.
Table 1 presents the samples’ crystallite sizes of the samples, calculated from the most prominent peak located at (311) plane of the spinel, the pore sizes, and the specific surface areas. The crystallite sizes are observed to vary with dopant amount, while the pore sizes and specific surface areas do not follow a particular order as regards the change in their values upon dopant addition.
Figure 3a–e show the High-Resolution Transmission Electron Microscope (HRTEM) images of the spinel nanoferrites. These spinels are generally crystalline nanoparticles. This is evident in the clear fringes observed on the HRTEM micrographs. The interplanar
d-spacing (
dhkl) fringes of the crystalline plane (
Figure 3a) can be attributed to the most prominent ring indexed (311) in the corresponding selected area electron diffraction (SAED) pattern and was calculated to be
d311 = 0.255 nm. This value varies slightly between samples due to changes in lattice parameters resulting from sites loading. The cubic spinel structure was highly expected to undergo both expansion and shrinkage between 0 ≤ x = y ≤ 0.3. The SAED pattern observed for these nanoferrites show bright crystalline spotted rings with the (311) exhibiting prominence. The SAED ring corresponding to the CeO
2 can also be observed very close to the most prominent ring of 311 peak. This agrees well with XRD result and the ring corresponds to (111) plane of CeO
2 fluorite structure [
37] whose inter-planar distance (
d-spacing) is 0.32 nm. It is noteworthy that the secondary phase of CeO
2 was only formed as Ce
3+ content increased in the spinel structure.
Figure 4a–h depicts the x-ray photoelectron spectroscopy (XPS) of the Co
1-2xNi
xMn
xFe
2-yCe
yO
4 nanoferrites with 0 ≤ x = y ≤ 0.3. The survey spectra shown in
Figure 4a,b before and after sputtering the surface with Ar
+ ions were identical with peaks confirming the presence of Co, Ni, Mn, Fe, Ce, O, and C elements in the Co
1-2xNi
xMn
xFe
2-yCe
yO
4. The sputtering took place for a minute at the sputtering rate of 15 nm per minute. The similarity of the survey spectra suggests the homogeneity of these nanoferrites. Observing from the survey, Ni 2p peak appears between the edge and trough of a peak. This may suggest the scarcity or total absence of Ni atoms in the sample composition. Furthermore, Ni 2p scan (
Figure 4d) confirms the scarcity of Ni atoms as the peak situated at 848.66 eV and 886.12 eV could not be ascribed to Ni.
Figure 4c shows Co 2p core level spectrum. The fitting peaks at 779.96 eV ascribed to B-site Co
2+ and 781.62 eV ascribed to A-site Co
2+ correspond to Co 2p
3/2 with its satellites at 785.11 eV, while the peak at 796.06 eV corresponds to Co 2p
1/2 with its satellite at 803.25 eV. It is obvious that there are no Co
3+ cations in the sample since they normally exhibit peaks centred at 798.5 eV and 783.3 eV while occupying tetrahedral and octahedral sites respectively [
38].
Figure 4e shows the Mn 2p core level spectrum consisting of spin-orbit doublets with their respective satellites.
The Mn 2p
3/2 peak at 640.68 eV ascribed to Mn
2+ cations and 642.06 eV ascribed to Mn
3+ cations has apparently combined satellites located at 645.75 eV. In addition, the peak at 652.96 eV is attributed to Mn 2p
1/2 with its satellite at 657.54 eV. The Ce 3d core level spectrum depicted in
Figure 4f exhibits peaks at 882.72 eV and 888.47 eV which are respectively attributable to Ce
3+ and Ce
4+ cations of Ce 3d
5/2. The peak at 901.88 eV is ascribed to Ce 3d
3/2 with its satellite peak at 916.21 eV. The Fe 2p spectrum shown in
Figure 4g comprises two spin-orbit doublets situated as follows: Fe 2p3/2 (710.43 eV, Fe3+ in B-sites and 712.45 eV, Fe3+ in A-sites) with its combined shake-up satellites at 718.20 eV and Fe 2p1/2 (724.10 eV, Fe3+ in B-sites and 725.61 eV, Fe3+ in A-sites) with its satellites at 732.32 eV. The positions of the binding energy of Co 2p and Fe 2p peaks and their satellites confirms the oxidation state of Co and Fe (Co2+ and Fe3+) in the sample. The spectra also indicate the presence of Co and Fe at the octahedral and tetrahedral sites [
19]. The O 1s spectrum (
Figure 4h) consists of three peaks ascribed to OL (centred at 529.00 eV), OV (centred at 529.9 eV) and OC (centred at 531.23 eV). OL is the oxygen inherent in the sample, OV signifies the oxygen deficient sites in the sample. OC is connected to surface oxygen on the sample owing to its exposure to air.
Various deconvolutions of O 1s peaks are presented in
Table 2. Chemisorbed oxygen is only present in two of the samples while the existence of oxygen vacancy and lattice oxygen is generally observed.
3.2. Gas sensing properties
The sensors were exposed to ethanol, propanol, ammonia, and LPG at operating temperature of 225 ºC. The concentration range for each of the target gases is based on the recommended exposure limits. Dry-air, helium, nitrogen, and argon were used as carrier-gases and for analyte gas dilution to produce different concentrations.
Figure 5 depicts the selectivity plots of Co
1-2xNi
xMn
xFe
2-yCe
yO
4 ,0 ≤ x=y ≤ 0.1, double-substitution spinel (i.e., S3) towards ethanol, propanol, ammonia, and LPG using dry-air and helium gas as carrier gases. It should be mentioned that the nitrogen and argon carrier-gas driven results are not included here as these did not yield good gas sensing performance as expected. These results are shown in
Figure S5 of the ESI.
Figure 5a shows that when dry air was used as a carrier and to dilute the analytes, sensors S1 and S3 were favourable to ethanol vapour. However, propanol vapour could be highly detected by S1 sensor. Moreover, sensor S3 is highly sensitive to ammonia. Furthermore, LPG could be detected by S1. Overall, sensor S3 gave the highest response of 7.95 to ammonia. The whole situation changed unexpectedly when helium is used as a carrier gas and to dilute the analytes. The expectation was that helium is an inert gas and will create an oxygen deficient atmosphere inside the sensing chamber resulting in the decrease in sensors’ sensitivity towards the analyte gas. The opposite was observed in
Figure 5b. The sensor S3 was extremely sensitive and selective to LPG at 6000 ppm concentration with a response of 116. Interestingly, these experiments that was designed to demonstrate the performance of semiconductor gas sensors in different atmospheres, operated with some accompanied anomalies.
Figure 5c depicts the transient current plot with time of the LPG concentrations on S3 sensor. The results show an increasing response with LPG concentration until 6000 ppm is reached thereafter the response begins to drop with further increase in concentration. However, on the 10000 ppm concentration the response is slightly higher than that towards the 9000 ppm. This can be attributed to the fact that on the 10000 ppm, the carrier gas is not present, and this is a maximum concentration. All the maximum gas analyte concentrations should look the same regardless of the carrier gas. The reason for dropping of the response after 6000 ppm is not clear as well as the high response when helium was used compared to dry air.
In
Figure 5d the transient current plot of the sensor toward LPG in dry air atmosphere, i.e., using dry air as carrier gas is shown. The sensor exhibits n-type characteristics towards the first three concentration cycles of LPG: 1000, 2000, and 3000 ppm. On the peak corresponding to the 3000 ppm concentration, a large split is observed, which is attributed to fluctuation in operating temperature as once observed by Mkwae et al. [
31]. However, from 4000 ppm onward, charge carrier transition began to occur in the sensor upon LPG sensing. The response of this sensor towards 3000 ppm of LPG is 3.35.
The LPG sensing of the sensor S3 in dry air atmosphere was characterized by two distinct phenomena, namely: current oscillation resulting from fluctuation in operating temperature, and charge carrier transition, from n- to p- type. The response, response time, and the recovery time of the sensor are 5.87 (10000ppm), 5.83 minutes, and 4.00 minutes respectively. The arrow at the bottom of the curve indicates the moment LPG was allowed into the chamber, while the arrow at the top of the curve indicates the moment LPG was stopped from entering the chamber. In helium gas atmosphere, there was no charge carrier transition observed in the sensors upon LPG detection, rather, there was a decrease in sensors’ response from 7000 ppm concentration of LPG in addition to the observed current oscillation which occurred at the peak as well as at the sensors’ relaxation. Generally, the sensors behaved abnormally upon LPG detection when helium was used as a carrier gas (i.e., in helium dominant chamber or atmosphere). The response increased from 1000 ppm to 6000 ppm then decreased the same way from 7000 ppm to 10000 ppm. This made the response of S3 at 10000 ppm (81.73) less than at 6000 ppm (116.43). The response could have been much greater than 116.43 had the abnormal decrease in response not occurred. Moreover, while the sensors’ response pattern towards other gases was normal, it was generally abnormal towards LPG (and occasionally, towards NH3) irrespective of the carrier gas.
Figure 6 shows the response patterns of the sensors towards LPG in dry air atmosphere and in helium gas atmosphere. The peaks labelled 1-10 corresponds to 1000-10000 ppm LPG concentration with 10 equal steps.
Figure 6a is the response pattern of the naturally dried pristine CoFe
2O
4-based sensor towards LPG in dry air atmosphere. This pattern is like that of the IR-dried sample. The pristine samples did not undergo charge carrier transition, rather, their LPG sensing performance was accompanied by current oscillation which becomes prominent with gas concentration. At concentrations of 5000 to 7000 ppm, the frequency of the oscillation was high, while the amplitude was higher from 8000 to 10000 ppm concentrations as seen in
Figure 6a. This current oscillation is caused by the fluctuating operating temperature during the LPG sensing. This anomaly was discussed elsewhere [
39]. The response pattern of the sensor (S5), fabricated from the sample with the chemical formula Co
0.4Ni
0.3Mn
0.3Fe
1.7Ce
0.3O
4 is depicted in
Figure 6b. This pattern is similar for all the doped samples. They exhibited charge carrier transition (n- to p- type) upon LPG detection.
Figure 6c,d show the response patterns of S2 and S5 respectively towards LPG in helium gas atmosphere. S1 and S2 were fabricated from the pure samples, and they exhibit similar sensitivity pattern, while S3, S4, and S5 are doped samples, and they also exhibit similar sensitivity pattern. In nitrogen gas atmosphere, the sensitivity pattern of the sensors towards LPG is similar to that of the dry air atmosphere - there was charge carrier transition of the sensors from n- to p-type. Also, there were series of current oscillation at the peaks of the response. In argon gas atmosphere, the performance of the sensors was generally poor and unquantifiable. The gas sensing performance of the sensor (S3) towards LPG in dry air atmosphere, nitrogen gas atmosphere, helium gas atmosphere, and argon gas atmosphere are compared in
Table 3. The resistance values presented in
Table 3 are the sensor resistance in the carrier gas before the target gas, i.e., LPG was introduced into the chamber. The idea of gas sensing using different carrier gases (as carried out in this work) is unprecedented. The outstanding response of the Co0.8Ni0.1Mn0.1Fe1.9Ce0.1O4 -based sensor (S3) towards LPG in helium gas atmosphere is noteworthy and could be suggested that helium gas created a more conduction-enhancing condition for the sensor than the dry air, nitrogen gas, and argon gas did. It probably was acting as a catalyst in the reaction between the sensing layer and the gas. Moreover, the optimal resistance for best performance appears to be attained in helium gas as can be observed in
Table 3. This could be a reason for higher conductivity in the sensor when the LPG was allowed into the chamber with helium as carrier gas. The responses presented in
Table 3 are towards 3000 ppm of LPG in dry air atmosphere, 6000 ppm of LPG in helium gas atmosphere, and 3000 ppm of LPG in nitrogen gas atmosphere because of the charge carrier transition that was observed on the sensor during LPG sensing in dry air atmosphere, and the decrease in response observed starting from 7000 ppm of LPG in helium gas atmosphere. The highest response recorded when dry air was used as a carrier gas was towards ammonia on the sensor S2. The response is 9 towards 1000 ppm; the transient current curve of the sensor is shown in
Figure 6e.
3.2.1. Gas sensing Mechanism
Figure 7 illustrates the gas sensing mechanism of the Co
0.8Ni
0.1Mn
0.1Fe
1.9Ce
0.1O
4-based sensor in the different ambient conditions (oxygen abundant and oxygen deficient) it was operated while determining its gas sensing performance. At the operating temperature of 225 °C, the atmospheric oxygen is chemisorbed on the surface of the sensor as
O - oxygen species (
Figure 7a) as illustrated in the equations 2-5. These oxygen species captured the free electrons at the conduction band, creating a depletion layer on the sensor [
40,
41], and thus increase the resistivity of the n-type sensor. When the target gas (LPG- a reducing gas) was introduced into the chamber, it reacted with the chemisorbed oxygen on the sensor surface to produce water and carbon dioxide (
Figure 7b). The surface reduction of the sensor’s sensing layer led to the release of the captured electrons, causing the thinning of the depletion layer, and resulting in increased conductivity of the sensor [
40,
41,
42]. In helium gas atmosphere (
Figure 7c), the sensor resistance (27.95 kΩ) was observed to be far less than observed in the air (see
Table 3). This implies the availability of more electrons at the conduction band of the n-type semiconductor-based sensor as a result of a very thin depletion region. This could be owing to the displacement of oxygen from the chamber where, at the operating temperature, it could have been adsorbed again to the sensor surface and reduced the quantity of electrons or impede the mobility of electrons in the conduction band. Nevertheless, it could not be explicitly stated that some chemisorbed species of oxygen are no longer available at all on the sensor surface, but rather that there is limited amount of them, otherwise, the current would have remained the same when LPG was allowed into the chamber, as adsorption/desorption of oxygen is currently understood to be the principle upon which the operation of electrochemical gas sensors is based [
20,
43,
44]. However, the sensor’s response increased with increasing concentration of the target gas (LPG). The effect of helium was more noticeable when used as a carrier gas for LPG (depicted in
Figure 7d). The relatively high sensor response towards helium-borne LPG indicates the catalytic effect of helium in the process. While helium did not actually participate in the chemical reaction (being inert as shown in
Figure 7c), it is obvious that it facilitated the reaction of LPG with the residual oxygen species on the sensing layer by easily bringing LPG into close contact with the sensing layer, and thus aiding further thinning of the depletion layer as shown in
Figure 7d.
Figure 8a shows the transient current curve resulting from the detection of 3 cycles each of 5000 ppm, 6000 ppm, and 7000 ppm LPG concentrations by the sensor (S3). The sensor could not maintain the same response towards three cycles of 5000 ppm concentration of LPG. Its response towards three cycles of 6000 ppm concentration is consistent, while its response to three cycles of 7000 ppm concentration is as observed for 5000 ppm concentration but with lower response. The decrease in the response that occur towards 7000 ppm, in comparison with the 6000-ppm concentration, is consistent with the result obtained earlier. The test was carried out over 10 months after the initial measurement. The response of the sensor has decreased to 12.39 over this period.
Figure 8b shows the sensitivity plot of the sensor (S3) towards 5000 ppm and 6000 ppm of LPG using helium as carrier gas in the presence of 50% relative humidity. The response of the sensor under this condition is 10.98 as compared to 12.39 under dry ambient condition. The sensor is not much affected by the high humidity; this indicates its potential for LPG detection and monitoring in an atmosphere with high relative humidity.
The long-term stability of the sensor (S3) in helium gas atmosphere is presented in
Figure 8c. The sensor retained about 50 % of its initial sensitivity towards LPG in 6 months. Surprisingly, in the following 5 months the response has drastically diminished relative to its initial response. This drastic diminishment in response could be the effect of series of current instability the sensor underwent while detecting LPG over this period. This instability, otherwise called current oscillation, occurred at the peaks (marked with red arrow in
Figure 6d) as well as at the relaxation (marked with green arrow in
Figure 6d). This could have caused an enormous stress on the sensor, and thereby rapidly reduced its efficiency.
The result of this work in comparison to other works on LPG detection is summarised in
Table 4. The performance of the sensor with helium as carrier is quite remarkable.